Question: Removal of mRNA contamination from scRNA-seq data using SoupX: How to proceed when "Extremely high contamination estimated"?
Entering edit mode
AA ▴ 10
Last seen 11 months ago

I am using SoupX to quantify, profile, and remove ambient mRNA contamination from my scRNA-seq data.

I've tried running the automated workflow (from the vignette) on my dataset and am getting the below message.

sc = load10X("/path/to/output") 
sc = autoEstCont(sc)
# 127 genes passed tf-idf cut-off and 35 soup quantile filter.  Taking the top 35.
# Using 252 independent estimates of rho.
# Estimated global rho of 0.75
# Error in setContaminationFraction(sc, contEst, forceAccept = forceAccept) : 
# Extremely high contamination estimated (0.75).  This likely represents a failure in estimating the contamination fraction.  Set forceAccept=TRUE to proceed with this value.

I'm very new to bioinformatics and scRNA-seq analysis and am wondering how to proceed. What should I do to check if this is "real" before moving on and correcting expression profile.

I've been trying to do some of the visual sanity checks such as mentioned in the vignette but it seems I first need to do the "manual method" to estimate the contamination fraction. However after reading through the vignette several times I'm still confused on the exact code I need to run. I keep running into error "'x' must be an array of at least two dimensions".

scRNAseq SoupX • 282 views
Entering edit mode
Last seen 4 hours ago
United States

This support site is intended for support of Bioconductor packages. If you want support for non-Bioconductor packages, you should contact the author directly.


Login before adding your answer.

Traffic: 375 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6