Souporcell: robust clustering of single-cell RNA-seq data by genotype without reference genotypes
Open accessAbstract
Methods to deconvolve single-cell RNA-sequencing (scRNA-seq) data are necessary for samples containing a mixture of genotypes, whether they are natural or experimentally combined. Multiplexing across donors is a popular experimental design that can avoid batch effects, reduce costs and improve doublet detection. By using variants detected in scRNA-seq reads, it is possible to assign cells to their donor of origin and identify cross-genotype doublets that may have highly similar transcriptional profiles, precluding detection by transcriptional profile. More subtle cross-genotype variant contamination can be used to estimate the amount of ambient RNA. Ambient RNA is caused by cell lysis before droplet partitioning and is an important confounder of scRNA-seq analysis. Here we develop souporcell, a method to cluster cells using the genetic variants detected within the scRNA-seq reads. We show that it achieves high accuracy on genotype clustering, doublet detection and ambient RNA estimation, as demonstrated across a range of challenging scenarios.
AI snapshot
Written from the paper’s full text.
TL;DR
souporcell clusters single-cell RNA-seq cells by genotype without reference genotypes, using variants detected directly in the reads. It jointly calls genotype clusters, cross-genotype doublets, and ambient RNA, and the paper reports that it outperforms demuxlet, vireo, and scSplit across synthetic, maternal/fetal, and *Plasmodium falciparum* benchmarks.
Key findings
- On a synthetic mixture of five human iPSC lines, souporcell assigned 6612/6622 singletons and 415/451 doublets correctly; four singletons were falsely labeled as doublets, 35 doublets were misidentified as singletons, and one doublet and four singletons were unassigned.
- For doublet detection on a synthetic mixture with 6% doublets and 10% ambient RNA, ROC AUC was 0.98 for souporcell and 0.91 for vireo; demuxlet’s posterior doublet probability was reported as poorly calibrated, with nearly 40% false positive doublets at its default threshold.
- Cell assignment versus demuxlet on three experimental iPSC mixture replicates gave Adjusted Rand Index 1.0 for souporcell and vireo and an average of 0.97 for scSplit.
- Across ambient RNA from 2.5%–50%, souporcell and vireo retained high accuracy, with souporcell more robust at calling doublets in high-ambient-RNA cases; all tools performed well down to a minority cluster of 1.2% (50 cells), but only souporcell and vireo correctly identified all minority sample singletons down to 0.6% of all cells.
- In maternal/fetal samples, souporcell clusters agreed with demuxlet and segregated with expected cell types, while vireo and scSplit had major discordances; in *P. falciparum* datasets, souporcell had superior ARI across the board, with vireo performing poorly on one dataset (ARI versus demuxlet 0.24).
Methods
- Reads are remapped with minimap2 rather than using STAR-aligned BAMs directly, putative SNPs are called with freebayes, and per-cell allele counts are obtained with vartrix.
- Clustering uses a sparse mixture model fit by a deterministic annealing variant of expectation maximization, with cluster centers represented as alternate allele fractions; doublets are detected with a beta-binomial model, and genotype plus ambient RNA are co-inferred in STAN by maximum likelihood/gradient descent.
- Validation used five HipSci human iPSC lines sequenced individually and in mixtures with 10x Chromium, synthetic mixtures with simulated doublets and ambient RNA, maternal/fetal samples from Vento et al., and three *P. falciparum* datasets containing six genetically distinct strains with 1893–2608 cells and median UMIs of ~1000.
Limitations
- The paper states that limitations include low signal-to-noise due to decreased UMI per cell and high numbers of donors causing increased local maxima.
- The text also makes evident that genotype accuracy from scRNA-seq is significantly lower than from whole-genome sequencing because most variants have very low coverage per cluster.
- The paper notes that no good generative models were available for batch effects, allele-specific expression, ambient RNA, and doublets in scRNA-seq, so realistic benchmarking required experimentally generated mixtures rather than fully in silico data.
Model-written summary. Check it against the paper before relying on it.
Integrity signals
Registry
Cites
Statistics
Claims the paper makes
4 of 5 cited sources checked · 2 supported · 0 flagged · accuracy not yet measured
- Findingown data
Demuxlet significantly overestimates doublets versus expectations based on the number of cells loaded, especially as ambient RNA increases.
“Demuxlet significantly overestimates doublets versus expectations based on the number of cells loaded11 (Table S2) especially as ambient RNA increases (Fig. 2j).”
Quote in the paperSubstantial claim1 of 1 sources checked- cites [11] Zheng, G. X. Y. et al. Massively parallel digital transcriptional profiling of single cells. Nat. Commun. 8, 14049 (2017). source
- Findingown data
The default doublet probability threshold for demuxlet gives nearly 40% false positive doublets.
“The default doublet probability threshold for demuxlet gives nearly 40% false positive doublets.”
Quote in the paperNot a substantial claim - Findingown data
Souporcell is able to accurately estimate the amount of ambient RNA in the experiment.
“With these data we also show that souporcell is able to accurately estimate the amount of ambient RNA in the experiment (Fig. 2k).”
Quote in the paperNot a substantial claim- The authors state this as an assumption: Allele counts for locus i of each cluster j are drawn from a binomial distribution with an alternative allele fraction of (1-ρ)fij+ρ*ai, where fij is 0, 0.5, or 1 (with a haploid mode limited to 0 and 1), ρ is the amount of ambient RNA, and ai is the average allele fraction in the experiment.
- Findingown data
The pattern of souporcell agreement with demuxlet and vireo/scSplit discordance is similar for the other samples tested.
“This is similar for the other samples tested (Fig. S2, Table S3).”
Quote in the paperNot a substantial claim - Findingown data
The ARI of scSplit and demuxlet suffered due to poor doublet detection.
“The ARI of scSplit and demuxlet suffered due to poor doublet detection.”
Quote in the paperSubstantial claim
- Findingown data
Souporcell surpasses both vireo and scSplit in genotype accuracy on a synthetically mixed sample with 6% doublets and 10% ambient RNA.
“Nevertheless, souporcell surpasses both vireo and scSplit in genotype accuracy on a synthetically mixed sample with 6% doublets and 10% ambient RNA (Fig. S1i).”
Quote in the paperSubstantial claim - Findingreasoning
Alignment parameter differences between STAR and aligners intended for variant calling constitute a set of artifacts.
“The next set of artifacts is alignment parameter differences between STAR and aligners intended for variant calling.”
Quote in the paperSubstantial claim - Findingown data
PCA of the normalized log likelihood matrix reveals clear separation of the five genotype clusters, with interspersed doublets.
“To visualize the five clusters identified by genotype we carried out a Principal Component Analysis (PCA) of the normalized log likelihood matrix, which reveals a clear separation of the clusters, with interspersed doublets (Fig. 2c and d).”
Quote in the paperNot a substantial claim - Findingown data
Comparing souporcell to demuxlet, there are 21 cells that demuxlet labels as maternal or fetal but which appear in the other individual's cell type clusters.
“Comparing souporcell to demuxlet, there are 21 cells that demuxlet labels as maternal or fetal but which appear in the other individual’s cell type clusters.”
Quote in the paperNot a substantial claim - Findingreasoning
The soft clip penalty being higher in STAR and not exposed as a parameter leads to false positive variants due to lack of soft clipping where other mappers would soft clip poorly matching read ends.
“The second type of artifact we found was due to the soft clip penalty being higher in STAR and not being exposed as a parameter to the user. This leads to false positive variants due to the lack of soft clipping where other mappers would soft clip poorly matching read ends.”
Quote in the paperSubstantial claim - Findingown data
STAR alignments introduce artifacts that are a significant source of false positive variants and reference bias.
“we have found several artifacts of the STAR alignments (methods) that are a significant source of false positive variants and reference bias.”
- Also stated in: § Supplementary Methods · plain
“We identified three different artifacts introduced by the STAR alignments resulting in false positive variants as well as reference bias that causes reads that do not support the reference allele to appear as though they do.”
Quote in the paperSubstantial claim - Also stated in: § Supplementary Methods · plain
- Findingreasoning
In the synthetic mixture, a low-dimensional representation of the expression matrix reveals relatively little variation because only one cell type is present.
“A low dimensional representation of the expression matrix, E, reveals relatively little variation as expected since there is only one cell type present (Fig. 2a).”
Quote in the paperNot a substantial claim - Interpretationown data
Demuxlet's poor ARI was due primarily to extremely high levels of false positive doublets.
“Again, demuxlet’s poor ARI was due primarily to extremely high levels of false positive doublets (Fig. 2l).”
Quote in the paperNot a substantial claim - Interpretationreasoning
An experimental design involving m overlapping mixtures for 2m-1 multiplexed individuals gives each individual a unique inclusion/exclusion signature across mixtures, allowing clusters to be assigned to individuals.
“To address this, we propose an experimental design involving m overlapping mixtures for 2m-1 multiplexed individuals (Table 1). Each individual is assigned a binary number from 1 to 2m, where each bit corresponds to the inclusion (1) or exclusion (0) from each of the mixtures. This gives each individual a unique signature of inclusion/exclusion across the mixtures.”
Quote in the paperSubstantial claim - Interpretationreasoning
Remapping reads with minimap2 produces alignments more conducive to accurate variant calling than STAR alignments.
“Instead, we remap the reads with minimap212 (Fig. 1a) which produces alignments more conducive to accurate variant calling.”
Quote in the paperSubstantial claim - Premisecitation
Ambient RNA shifts the observed allele fraction away from the underlying genotype allele fractions.
“The ambient RNA shifts the observed allele fraction away from the underlying genotype allele fractions9.”
Quote in the paperSubstantial claim0 of 1 sources checked- cites [9] Young, M. D. & Behjati, S. SoupX removes ambient RNA contamination from droplet based single cell RNA sequencing data. Preprint at bioRxiv https://doi.org/10.1101/303727(2018). source
- Interpretationreasoning
The Plasmodium sample where vireo performed poorly was more difficult due to sample skew caused by a clonal expansion of one of the six strains.
“This sample was more difficult due to sample skew caused by a clonal expansion of one of the six strains.”
Quote in the paperNot a substantial claim - Interpretationreasoning
Mixture model clustering has the advantage over hard clustering that cells can be partially assigned to multiple clusters, which naturally allows for both doublet cells and varying levels of ambient RNA.
“The advantage of mixture model clustering over hard clustering is that cells can be partially assigned to multiple clusters, which naturally allows for both doublet cells and varying levels of ambient RNA (Fig. 1d).”
Quote in the paperSubstantial claim- The authors state this as an assumption: Allele counts for locus i of each cluster j are drawn from a binomial distribution with an alternative allele fraction of (1-ρ)fij+ρ*ai, where fij is 0, 0.5, or 1 (with a haploid mode limited to 0 and 1), ρ is the amount of ambient RNA, and ai is the average allele fraction in the experiment.
- Interpretationreasoning
The default parameters of cellranger make its output poorly suited for variant calling.
“This is a parameter which is exposed to the user but, the default makes the output of cellranger poorly suited for variant calling.”
Quote in the paperSubstantial claim - Interpretationnone given
Souporcell explicitly models and estimates the amount of ambient RNA in the experiment.
“Furthermore, souporcell explicitly models and estimates the amount of ambient RNA in the experiment, which is a major confounder of scRNAseq analysis with regard to both expression and genotype.”
Quote in the paperSubstantial claim - Interpretationown data
Based on the position of the 21 discordant cells in the expression t-SNE plot, it is most likely that these are errors in the demuxlet assignments that are not made by souporcell.
“Based on the position of these cells in the expression t-SNE plot, it is most likely that these are errors in the demuxlet assignments that are not made by souporcell.”
Quote in the paperSubstantial claim - Interpretationreasoning
Souporcell enables donor multiplexing designs to be used more easily than was previously possible, including in situations when no WGS or genotyping data are available.
“souporcell enables donor multiplexing designs to be used more easily than was previously possible, including in situations when no WGS or genotyping data are available.”
Quote in the paperSubstantial claim
- May combine two claims. 21
- Classified as probably not a substantial claim. 3
May combine two claims.
PremisecitationVartrix avoids reference bias due to ambiguous support such as alignment end effects and corrects bases with duplicate reads via the UMI.
“Next, we count alleles per cell with vartrix14(Fig. 1c) which avoids reference bias due to ambiguous support such as alignment end effects and corrects bases with duplicate reads via the UMI.”
Quote in the paperSubstantial claim1 of 1 sources checked- cites [14] Petti, A. A. et al. A general approach for detecting expressed mutations in AML cells using single cell RNA-sequencing. Nat. Commun. 10, 3660 (2019). source
The Methods section (‘Single-cell mutation identification and analysis’; page unavailable) identifies VarTrix as a tool for automated cell-variant assignment. It also describes resolving discordant duplicate reads by UMI, but attributes that procedure to the authors’ Pysam-based workflow, not to VarTrix. Our copy does not say that VarTrix avoids reference bias or alignment-end effects. (LLM, conf=0.85)
Experimental — accuracy not yet measured
Probably not a substantial claim: may not pass S3 (load-bearing).
Interpretationnone givenIf using known variable sites such as from the 1k genomes project or if the genotypes of the donors are known, the remapping process is significantly less important.
“If using known variable sites such as from the 1k genomes project24 or if the genotypes of the donors are known, the remapping process is significantly less important.”
Quote in the paperNot a substantial claim1 of 1 sources checked- cites [24] 1000 Genomes Project Consortium. A global reference for human genetic variation. Nature 526, 68–74 (2015). source
Probably not a substantial claim: may not pass S3 (load-bearing).
PremisecitationGATK recommendations for variant calling on bulk RNAseq involve removing these regions of the alignments prior to variant calling.
“This has been noted by others in the past, and GATK recommendations for variant calling on bulk RNAseq involve removing these regions of the alignments prior to variant calling23.”
Quote in the paperNot a substantial claim1 of 1 sources checked- cites [23] Sahraeian, S. M. E. et al. Gaining comprehensive biological insight into the transcriptome by performing a broad-spectrum RNA-seq analysis. Nat. Commun. 8, 59 (2017). source
The Results section (‘Variant calling’; page unavailable) says GATK applies RNA-specific filtering steps to handle splicing in RNA-seq data. It does not say that GATK recommends removing the regions in question from alignments before variant calling, so that specific claim is unsupported. (LLM, conf=0.82)
Experimental — accuracy not yet measured
May combine two claims.
Findingown dataAnalysis of the expression profile of one Plasmodium dataset reveals that the genotypes are distributed across the Plasmodium intra-erythrocytic cycle while being well separated in normalized loss cluster space.
“Analysis of the expression profile of one of these (see Fig. S3 for the others) reveals that the genotypes are distributed across the Plasmodium intra-erythrocytic cycle (Fig. 3b) while being well separated in normalized loss cluster space (Fig. 3c,d).”
Quote in the paperSubstantial claimMay combine two claims.
Findingown dataFor the synthetic mixture, souporcell assigned 6612/6622 singletons and 415/451 doublets correctly; four singletons were falsely labeled as a doublet, 35 doublets were misidentified as singletons, and one doublet and four singletons were unassigned.
“For these data souporcell assigned 6612/6622 singletons and 415/451 doublets correctly; four singletons were falsely labeled as a doublet, 35 doublets were misidentified as singletons, and one doublet and four singletons were unassigned.”
Quote in the paperSubstantial claimMay combine two claims.
FindingreasoningSTAR handles spliced reads that do not match the reference well by introducing multiple splice events that force the read to fit the reference in a statistically spurious fashion, causing variant callers to count these reads as supporting the reference allele and introducing reference bias.
“The first artifact is due to the way STAR handles spliced reads when the read does not match the reference well. STAR will take such a read and introduce multiple splice events that force it to fit the reference in a statistically spurious fashion.”
- Also stated in: § Supplementary Methods · plain
“Consequently, for variants in these loci, the variant callers will count these reads as supporting the reference allele, thereby introducing reference bias and noise to the downstream clustering.”
Quote in the paperSubstantial claim- Also stated in: § Supplementary Methods · plain
May combine two claims.
Findingown dataOn a synthetic mixture with 6% doublets and 10% ambient RNA, the ROC AUC for doublet calls was 0.98 for souporcell and 0.91 for vireo.
“we calculated a receiver-operator characteristic (ROC) curve of the doublet calls (Fig. 2i) on a synthetic mixture with 6% doublets and 10% ambient RNA that showed the area under the curve values of 0.98 and 0.91 for souporcell and vireo, respectively.”
Quote in the paperSubstantial claimMay combine two claims.
Findingown dataThe Adjusted Rand Index of remaining cell assignments versus demuxlet was 1.0 for souporcell and vireo across the three replicates and an average of 0.97 for scSplit.
“The Adjusted Rand Index (ARI) of the remaining cell assignments versus demuxlet (Table S2) were 1.0 (fully concordant) for souporcell and vireo across the three replicates and an average of 0.97 for scSplit.”
Quote in the paperSubstantial claimMay combine two claims.
Findingown dataSouporcell and vireo retain high accuracy across ambient RNA levels, with souporcell being more robust at accurately calling doublets in high ambient RNA cases.
“We found that souporcell and vireo retain high accuracy with souporcell being more robust at accurately calling doublets in high ambient RNA cases (Fig. 2k).”
Quote in the paperSubstantial claimMay combine two claims.
Findingown dataSouporcell outperforms vireo and scSplit and surpasses demuxlet on both cell assignment and doublet accuracy.
“We show that souporcell not only outperforms these new methods, but also surpasses demuxlet on both cell assignment and doublet accuracy.”
- Also stated in: § Discussion · plain
“Our benchmarks show that souporcell can outperform all other currently available methods, including those that require genotypes a priori.”
Quote in the paperSubstantial claim- Also stated in: § Discussion · plain
May combine two claims.
Findingown dataDonor of origin is the most significant driver of expression, but donor cells overlap in expression patterns and it is not possible to assign a donor to each cell based solely on expression patterns.
“Indeed, the most significant driver of expression appears to be the donor of origin, but the donor cells overlap in expression patterns and it is not possible to assign a donor to each cell based solely on expression patterns.”
Quote in the paperSubstantial claimMay combine two claims.
Findingown dataSouporcell is highly accurate for challenging datasets involving closely related maternal/fetal samples, and varying mixtures of Plasmodium falciparum strains.
“Moreover, souporcell is highly accurate for challenging datasets involving closely related maternal/fetal samples, and varying mixtures of Plasmodium falciparum strains.”
Quote in the paperSubstantial claimMay combine two claims.
Findingown dataSouporcell clusters agree with demuxlet and segregate with the expected cell type clusters, while vireo and scSplit have major discordances with demuxlet.
“While souporcell clusters agree with demuxlet and segregate with the expected cell type clusters, vireo and scSplit have major discordances with demuxlet.”
Quote in the paperSubstantial claimMay combine two claims.
Findingown dataAll tools performed well down to a minority cell cluster comprising only 1.2% (50 cells) of total cells, but only souporcell and vireo were able to correctly identify all minority sample singletons as their own cluster down to 0.6% of all cells.
“We found that all tools performed well down to the minority cell cluster comprising only 1.2% (50 cells) of total cells (Fig. 2m), but only souporcell and vireo were able to correctly identify all minority sample singletons as their own cluster down to 0.6% of all cells.”
Quote in the paperSubstantial claimMay combine two claims.
FindingreasoningThe most common error mode for vireo and scSplit is calling homozygous reference loci as heterozygous variants, which is expected when ambient RNA is not accounted for.
“The most common error mode for vireo and scSplit is calling homozygous reference loci as heterozygous variants (Fig. S1j) which is expected when ambient RNA is not accounted for, as it is not in these two tools.”
Quote in the paperSubstantial claimMay combine two claims.
Findingown dataThe ARI for each method on the three Plasmodium datasets shows superior performance for souporcell across the board, with scSplit suffering on all datasets and vireo performing poorly on one dataset with an ARI versus demuxlet of 0.24.
“The ARI for each method (Table S4) on the three Plasmodium data sets show superior performance for souporcell across the board, with scSplit suffering on all datasets and vireo performing poorly on one, which had an ARI versus demuxlet of 0.24.”
Quote in the paperSubstantial claimMay combine two claims.
Findingown dataDemuxlet's posterior doublet probability output did not have enough significant digits and is 1.0 until it starts varying with 27% false positives.
“Demuxlet’s posterior doublet probability output did not have enough significant digits and is 1.0 until it starts varying with 27% false positives.”
Quote in the paperNot a substantial claimMay combine two claims.
FindingreasoningThe indel penalty relative to the mismatch penalty is much higher in STAR than other aligners, causing alignments to choose many mismatches over a single or few indels and thus create false positive variants.
“The final issue is that the indel penalty relative to the mismatch penalty is much higher in STAR than other aligners. This causes the alignments to choose many mismatches over a single or few indels when possible and thus create false positive variants.”
Quote in the paperSubstantial claimMay combine two claims.
Findingnone givenLimitations of souporcell include low signal to noise due to decreased UMI per cell and high numbers of donors causing increased local maxima.
“Limitations of souporcell include low signal to noise due to decreased UMI per cell and high numbers of donors causing increased local maxima.”
Quote in the paperNot a substantial claimMay combine two claims.
Interpretationown dataSouporcell is robust across a large range of parameters, and more so than any other currently available method.
“Using more realistic and challenging test cases than previous studies, we show that souporcell is robust across a large range of parameters, and more so than any other currently available method.”
Quote in the paperSubstantial claimMay combine two claims.
Interpretationown dataThe expression PCA and normalized cell-cluster loss PCA of the experimental mixture were similar to the synthetic mixture, indicating that the synthetic mixtures were an accurate approximation of real mixtures.
“The expression PCA (Fig. 2e) and normalized cell-cluster loss PCA (Fig. 2g,h) of the experimental mixture were similar to the synthetic mixture indicating that the synthetic mixtures were an accurate approximation of real mixtures.”
Quote in the paperSubstantial claimMay combine two claims.
Interpretationnone givenSouporcell clusters cells by genotype, calls doublet-cell barcodes, and infers the amount of ambient RNA in the experiment without the use of a genotype reference.
“We present souporcell, a method to cluster cells by genotype, call doublet-cell barcodes, and infer the amount of ambient RNA in the experiment without the use of a genotype reference.”
Quote in the paperSubstantial claimProbably not a substantial claim: may not pass S4 (based).
Interpretationnone givenFiltering loci to those with at least n cells supporting each allele provides fairly robust SNPs that have a good chance of aiding the clustering process.
“This provides us with fairly robust SNPs that have a good chance of aiding the clustering process.”
Quote in the paperNot a substantial claimMay combine two claims.
InterpretationreasoningSouporcell enables valuable genotype information to be extracted and ambient RNA estimation at no additional cost.
“souporcell also enables valuable genotype information to be extracted and ambient RNA estimation at no additional cost.”
Quote in the paperSubstantial claim
The authors state this as an assumption.
Allele counts for locus i of each cluster j are drawn from a binomial distribution with an alternative allele fraction of (1-ρ)fij+ρ*ai, where fij is 0, 0.5, or 1 (with a haploid mode limited to 0 and 1), ρ is the amount of ambient RNA, and ai is the average allele fraction in the experiment.
“we assume that the allele counts for locus i of each cluster j are drawn from a binomial distribution with an alternative allele fraction of (1-ρ)fij+ρ*ai, where fij is 0, 0.5, or 1 (with a haploid mode limited to 0 and 1), ρ is a parameter representing the amount of ambient RNA and ai is the average allele fraction in the experiment.”
The paper relies on this without support.
There are no good generative models available for batch effects, allele-specific expression, ambient RNA, and doublets in scRNAseq that can be used to generate in silico data for testing methods that cluster by genotype.
“Currently, there are no good generative models available for batch effects, allele-specific expression, ambient RNA, and doublets in scRNAseq that can be used to generate in silico data for testing methods that cluster by genotype.”
The paper relies on this without support.
Malaria infections often contain parasites from multiple different genetic backgrounds, and it is not possible to separate the strains prior to sequencing.
“Malaria infections often contain parasites from multiple different genetic backgrounds, and it is not possible to separate the strains prior to sequencing.”
The paper relies on this without support.
The decidua is primarily composed of maternal cells with some invading fetal trophoblasts, while the placenta is largely composed of cells of fetal origin with the exception of maternal macrophages.
“The decidua is primarily composed of maternal cells with some invading fetal trophoblasts, while the placenta is largely composed of cells of fetal origin with the exception of maternal macrophages.”
The paper relies on this without support.
In scRNAseq data most variants have very low coverage per cluster compared to what would be generated from WGS data, so genotype accuracy is significantly lower than one would attain with genome sequencing.
“In scRNAseq data most variants have very low coverage per cluster compared to what would be generated from WGS data, thus the genotype accuracy is significantly lower than one would attain with genome sequencing.”
The authors state this as an assumption.
The majority of barcodes contain a single cell.
“It is important for some assumptions in this model that the majority of barcodes contain a single cell.”
The authors state this as an assumption.
Only biallelic variants are considered.
“We will assume only biallelic variants.”
The authors state this as an assumption.
Ploidy is limited to 1 or 2.
“We assume ploidy is limited to 1 or 2.”
The authors state this as an assumption.
Each cluster is independent and each locus is independent.
“We treat each cluster as independent and each locus as independent, before marginalizing across the possible genotypes.”
The authors state this as an assumption.
The observed allele fractions are drawn from a binomial distribution with a probability skewed away from p = g/P by the level of ambient RNA ρ.
“The observed allele fractions are assumed to have been drawn from a binomial distribution with a probability that was skewed away from p = g/P by the level of ambient RNA ρ.”
The authors state this as an assumption.
Allele counts of each cell at each locus are drawn from a beta-binomial distribution from either a single cluster or a pair of clusters.
“We treat the allele counts of each cell at each locus as random variables drawn from a beta-binomial distribution from either a single cluster or a pair of clusters.”
The authors state this as an assumption.
The expected allele fractions of a doublet coming from cluster i and cluster j is the average of the allele fractions of the two clusters.
“The expected allele fractions of a doublet coming from cluster i, and cluster j is the average of the allele fractions of the two clusters.”
The authors state this as an assumption.
The beta-binomial is used to model uncertainty in the binomial parameter p.
“The beta-binomial is used to model our uncertainty in the binomial parameter p.”
The authors state this as an assumption.
Each locus in each cluster comes from one of three genotypes for diploid (0/0, 0/1, 1/1) and two genotypes for haploid (0, 1).
“The model treats each locus in each cluster as coming from one of three genotypes for diploid (0/0, 0/1, 1/1, here denoted by g = 0, 1, or 2) and two genotypes from haploid (0, 1).”
Claims other papers make about this one
Each is checked against this paper’s own text: does it say what the citing paper says it does?
Outside of human samples, SNP-based demultiplexing has been applied to demultiplex plasmodium samples and across mouse strains.
“Outside of human samples, SNP-based demultiplexing has been applied to demultiplex other species: plasmodium samples, and across mouse strains (Heaton et al, 2020; Mylka et al, 2022).”
From Evaluation of genetic demultiplexing of single-cell sequencing data from model species (2023) · cites this paper as Heaton et al, 2020
The Validation and Benchmarking section reports using souporcell to separate six genetically distinct Plasmodium falciparum strains, supporting the Plasmodium example. The cited paper does not report demultiplexing across mouse strains, so it supports only one of the two examples in the sentence. (LLM, conf=0.87)
Experimental — accuracy not yet measuredSoupOrCell uses the identification of conflicting genotypes to identify ambient RNA contamination, limiting its application to mixed-genotype experiments.
“SoupOrCell [21] uses the identification of conflicting genotypes to identify ambient RNA contamination, limiting its application to mixed-genotype experiments.”
From SoupX removes ambient RNA contamination from droplet-based single-cell RNA sequencing data (2020) · cites this paper as [21]
Self-check
Percentages, sums and reported differences re-computed from the paper’s own numbers.
Each contradiction between two passages listed here was confirmed by an independent check.
Reference-allele rate is defined inconsistently
MinorThe methods define g as an integer number of reference-allele copies, but also call g the expected reference-allele rate. A later passage gives that rate as g/P, which differs from g when the ploidy is two and g is one or two.
“The expected reference allele rate without ambient RNA is g and g is an integer value ∈ [0..P].”
“The observed allele fractions are assumed to have been drawn from a binomial distribution with a probability that was skewed away from p = g/P by the level of ambient RNA ρ.”
Cluster parameter switches between reference and alternate allele fractions
MinorThe definition says φ is the fraction of reference alleles, whereas Equation 1 applies φ to the alternate-allele count and (1−φ) to the reference-allele count. Those are opposite definitions of the same parameter.
“This is a real number representing the fraction of ref alleles in this cluster at this locus.”
“Equation 1: Cluster model Likelihood function ℒ(A)=∏c∈C∑k∈K1K∏l∈Lc(Al,c,0+Al,c,1Al,c,1)ϕk,lAl,c,1(1−ϕk,l)Al,c,0”
Cites
STAR: ultrafast universal RNA-seq aligner
Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., Batut, P., Chaisson, M., & Gingeras, T. R. (2013) · Bioinformatics · doi:10.1093/bioinformatics/bts635 · Cited by 51,709 · see citing works
A global reference for human genetic variation
The 1000 Genomes Project Consortium, authors, C., Auton, A., Abecasis, G. R., Steering committee, Altshuler, D. M., Durbin, R. M., Abecasis, G. R., Bentley, D. R., Chakravarti, A., et al. (2015) · Nature · doi:10.1038/nature15393 · Cited by 18,247 · see citing works
Self-citationMinimap2: Pairwise alignment for nucleotide sequences
Li, H. (2018) · Bioinformatics · doi:10.1093/bioinformatics/bty191 · Cited by 14,807 · see citing works
<i>Stan</i> : A Probabilistic Programming Language
Carpenter, B., Gelman, A., Hoffman, M. D., Lee, D., Goodrich, B., Betancourt, M., Brubaker, M., Guo, J., Li, P., & Riddell, A. (2017) · Journal of Statistical Software · doi:10.18637/jss.v076.i01 · Cited by 7,754 · cited 2× · see citing works
“This model is implemented in the domain-specific language for probabilistic models, STAN 16 , and it solves for the maximum-likelihood soup fraction with gradient descent.”
MethodMassively parallel digital transcriptional profiling of single cells
Zheng, G. X. Y., Terry, J. M., Belgrader, P., Ryvkin, P., Bent, Z. W., Wilson, R., Ziraldo, S. B., Wheeler, T. D., McDermott, G. P., Zhu, J., et al. (2017) · Nature Communications · doi:10.1038/ncomms14049 · Cited by 7,353 · see citing works
Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets.
Macosko, E. Z., Basu, A., Satija, R., Nemesh, J., Shekhar, K., Goldman, M., Tirosh, I., Bialas, A. R., Kamitaki, N., Martersteck, E. M., et al. (2015) · Cell · doi:10.1016/j.cell.2015.05.002 · Cited by 7,110 · see citing works
Haplotype-based variant detection from short-read sequencing
Garrison, E., & Marth, G. (2012) · arXiv · doi:10.48550/arxiv.1207.3907 · Cited by 3,700 · cited 1× · see citing works
“Until recently, a genotype reference obtained via whole-genome or exome sequencing has been required for each multiplexed individual before cell-sample categorization 6 .”
Single-cell reconstruction of the early maternal–fetal interface in humans
Vento-Tormo, R., Efremova, M., Botting, R. A., Turco, M. Y., Vento-Tormo, M., Meyer, K. B., Park, J.-E., Stephenson, E., Polański, K., Goncalves, A., et al. (2018) · Nature · doi:10.1038/s41586-018-0698-6 · Cited by 2,358 · see citing works
SoupX removes ambient RNA contamination from droplet based single-cell RNA sequencing data
Young, M. D., & Behjati, S. (2018) · doi:10.1101/303727 · Cited by 1,648 · see citing works
SoupX removes ambient RNA contamination from droplet-based single-cell RNA sequencing data
Young, M. D., & Behjati, S. (2020) · GigaScience · doi:10.1093/gigascience/giaa151 · Cited by 1,619 · see citing works
Multiplexed droplet single-cell RNA-sequencing using natural genetic variation
Kang, H. M., Subramaniam, M., Targ, S., Nguyen, M., Maliskova, L., McCarthy, E., Wan, E., Wong, S., Byrnes, L., Lanata, C. M., et al. (2018) · Nature Biotechnology · doi:10.1038/nbt.4042 · Cited by 1,339 · see citing works
Stoeckius, M., Zheng, S., Houck-Loomis, B., Hao, S., Yeung, B. Z., Mauck, W. M., Smibert, P., & Satija, R. (2018) · Genome Biology · doi:10.1186/s13059-018-1603-1 · Cited by 1,046 · see citing works
Haplotyping germline and cancer genomes with high-throughput linked-read sequencing
Zheng, G. X. Y., Lau, B. T., Schnall-Levin, M., Jarosz, M. V., Bell, J., Hindson, C. M., Kyriazopoulou-Panagiotopoulou, S., Masquelier, D. A., Merrill, L., Terry, J. M., et al. (2016) · Nature Biotechnology · doi:10.1038/nbt.3432 · Cited by 678 · see citing works
Self-citationMULTI-seq: sample multiplexing for single-cell RNA sequencing using lipid-tagged indices
McGinnis, C. S., Patterson, D. M., Winkler, J., Conrad, D. N., Hein, M. Y., Srivastava, V., Hu, J. L., Murrow, L. M., Weissman, J. S., Werb, Z., et al. (2019) · Nature Methods · doi:10.1038/s41592-019-0433-8 · Cited by 643 · see citing works
Common genetic variation drives molecular heterogeneity in human iPSCs
Kilpinen, H., Goncalves, A., Leha, A., Afzal, V., Alasoo, K., Ashford, S., Bala, S., Bensaddek, D., Casale, F. P., Culley, O. J., et al. (2017) · Nature · doi:10.1038/nature22403 · Cited by 586 · see citing works
Self-citationSahraeian, S. M. E., Mohiyuddin, M., Sebra, R., Tilgner, H., Afshar, P. T., Au, K. F., Bani Asadi, N., Gerstein, M. B., Wong, W. H., Snyder, M. P., et al. (2017) · Nature Communications · doi:10.1038/s41467-017-00050-4 · Cited by 376 · see citing works
Batch effects and the effective design of single-cell gene expression studies
Tung, P.-Y., Blischak, J. D., Hsiao, C. J., Knowles, D. A., Burnett, J. E., Pritchard, J. K., & Gilad, Y. (2017) · Scientific Reports · doi:10.1038/srep39921 · Cited by 356 · see citing works
Vireo: Bayesian demultiplexing of pooled single-cell RNA-seq data without genotype reference
Huang, Y., McCarthy, D. J., & Stegle, O. (2019) · Genome Biology · doi:10.1186/s13059-019-1865-2 · Cited by 348 · see citing works
Indexing Graphs for Path Queries with Applications in Genome Research
Sirén, J., Välimäki, N., & Mäkinen, V. (2014) · IEEE Transactions on Computational Biology and Bioinformatics · doi:10.1109/tcbb.2013.2297101 · Cited by 316 · see citing works
The Malaria Cell Atlas: Single parasite transcriptomes across the complete Plasmodium life cycle
Howick, V. M., Russell, A. J. C., Andrews, T. S., Heaton, H., Reid, A. J., Natarajan, K. N., Butungi, H., Metcalf, T., Verzier, L. H., Rayner, J. C., et al. (2019) · Science · doi:10.1126/science.aaw2619 · Cited by 302 · see citing works
Self-citationA general approach for detecting expressed mutations in AML cells using single cell RNA-sequencing
Petti, A. A., Williams, S. R., Miller, C. A., Fiddes, I. T., Srivatsan, S. N., Chen, D. Y., Fronick, C. C., Fulton, R. S., Church, D. M., & Ley, T. J. (2019) · Nature Communications · doi:10.1038/s41467-019-11591-1 · Cited by 215 · see citing works
Co‐evolution of NK receptors and HLA ligands in humans is driven by reproduction
Moffett, A., & Colucci, F. (2015) · Immunological Reviews · doi:10.1111/imr.12323 · Cited by 174 · see citing works
Genotype-free demultiplexing of pooled single-cell RNA-seq
Xu, J., Falconer, C., Nguyen, Q., Crawford, J., McKinnon, B. D., Mortlock, S., Senabouth, A., Andersen, S., Chiu, H. S., Jiang, L., et al. (2019) · Genome Biology · doi:10.1186/s13059-019-1852-7 · Cited by 94 · see citing works
The human-induced pluripotent stem cell initiative-data resources for cellular genetics
Streeter, I., Harrison, P. W., Faulconbridge, A., The HipSci Consortium, Flicek, P., Parkinson, H., & Clarke, L. (2017) · Nucleic Acids Research · doi:10.1093/nar/gkw928 · Cited by 92 · see citing works
Deterministic Annealing Variant of the EM Algorithm
Ueda, N., & Nakano, R. (1994) · neural information processing systems · Cited by 72 · see citing works
Flagged first, then by citations
Cited by
Sorted by citations
SoupX removes ambient RNA contamination from droplet-based single-cell RNA sequencing data
Young, M. D., & Behjati, S. (2020) · GigaScience · doi:10.1093/gigascience/giaa151 · Cited by 1,619 · see citing works
“SoupOrCell [21] uses the identification of conflicting genotypes to identify ambient RNA contamination, limiting its application”
Single-cell multi-omics analysis of the immune response in COVID-19
Cambridge Institute of Therapeutic Immunology and Infectious Disease-National Institute of Health Research (CITIID-NIHR) COVID-19 BioResource Collaboration, Stephenson, E., Reynolds, G., Botting, R. A., Calero-Nieto, F. J., Morgan, M. D., Tuong, Z. K., Bach, K., Sungnak, W., Worlock, K. B., et al. (2021) · Nature Medicine · doi:10.1038/s41591-021-01329-2 · Cited by 742 · see citing works
An integrated cell atlas of the lung in health and disease
Sikkema, L., Ramírez-Suástegui, C., Strobl, D. C., Gillett, T. E., Zappia, L., Madissoon, E., Markov, N. S., Zaragosi, L., Ji, Y., Ansari, M., et al. (2023) · Nature Medicine · doi:10.1038/s41591-023-02327-2 · Cited by 856 · see citing works
“Both annotation sources (for example, Ensembl or RefSeq) and annotation versions (for example, Ensembl release 84 or Ensembl release 91) contribute to the variation between different gene naming schemes.”
BackgroundCircuits between infected macrophages and T cells in SARS-CoV-2 pneumonia
Grant, R. A., Morales‐Nebreda, L., Markov, N. S., Swaminathan, S., Querrey, M., Rios-Guzman, E., Abbott, D. A., Donnelly, H. K., Donayre, A., Goldberg, I. A., et al. (2021) · Nature · doi:10.1038/s41586-020-03148-w · Cited by 703 · see citing works
SARS-CoV-2 infection triggers profibrotic macrophage responses and lung fibrosis
Wendisch, D., Dietrich, O., Mari, T., von Stillfried, S., Ibarra, I. L., Mittermaier, M., Mache, C., Chua, R. L., Knoll, R., Timm, S., et al. (2021) · Cell · doi:10.1016/j.cell.2021.11.033 · Cited by 528 · see citing works
“Demultiplexing of donors Donors of PBMCs were demultiplexed using the souporcell algorithm (Heaton et al., 2020) based on the BAM files & barcodes exported by CellRanger count.”
Spatially resolved multiomics of human cardiac niches
Kanemaru, K., Cranley, J., Muraro, D., Miranda, A. M. A., Ho, S. Y., Wilbrey-Clark, A., Patrick Pett, J., Polanski, K., Richardson, L., Litvinukova, M., et al. (2023) · Nature · doi:10.1038/s41586-023-06311-1 · Cited by 425 · see citing works
“Part of the SAN samples were mixed with different donors after the nuclei isolation for cost-efficient experimental design (Supplementary Table 8) and computationally demultiplexed (Soupercell, v.2.0) 55 on the basis of genetic variation between the donors.”
InfluentialSingle-nucleus cross-tissue molecular reference maps toward understanding disease gene function
Eraslan, G., Drokhlyansky, E., Anand, S. K., Fiškin, E., Subramanian, A., Slyper, M., Wang, J., Wittenberghe, N. V., Rouhana, J. M., Waldman, J., et al. (2022) · Science · doi:10.1126/science.abl4290 · Cited by 396 · see citing works
A human cell atlas of fetal chromatin accessibility
Domcke, S., Hill, A. J., Daza, R. M., Cao, J., O’Day, D. R., Pliner, H. A., Aldinger, K. A., Pokholok, D., Zhang, F., Milbank, J. H., et al. (2020) · Science · doi:10.1126/science.aba7612 · Cited by 375 · see citing works
Local and systemic responses to SARS-CoV-2 infection in children and adults
Yoshida, M., Worlock, K. B., Huang, N., Lindeboom, R. G. H., Butler, C. R., Kumasaka, N., Dominguez Conde, C., Mamanova, L., Bolt, L., Richardson, L., et al. (2022) · Nature · doi:10.1038/s41586-021-04345-x · Cited by 365 · see citing works
Mapping the developing human immune system across organs
Suo, C., Dann, E., Goh, I., Jardine, L., Kleshchevnikov, V., Park, J.-E., Botting, R. A., Stephenson, E., Engelbert, J., Tuong, Z. K., et al. (2022) · Science · doi:10.1126/science.abo0510 · Cited by 315 · see citing works
“0) (78) (for more details, see the supplementary materials).”
BackgroundCRISPR activation and interference screens decode stimulation responses in primary human T cells
Schmidt, R., Steinhart, Z., Layeghi, M., Freimer, J. W., Bueno, R., Nguyen, V., Blaeschke, F., Ye, C. J., & Marson, A. (2022) · Science · doi:10.1126/science.abj4008 · Cited by 285 · see citing works
Spatial multiomics map of trophoblast development in early pregnancy
Arutyunyan, A., Roberts, K., Troulé, K., Wong, F. C. K., Sheridan, M. A., Kats, I., García‐Alonso, L., Velten, B., Hoo, R., Ruiz-Morales, E. R., et al. (2023) · Nature · doi:10.1038/s41586-023-05869-0 · Cited by 309 · see citing works
Datlinger, P., Rendeiro, A. F., Boenke, T., Senekowitsch, M., Krausgruber, T., Barreca, D., & Bock, C. (2021) · Nature Methods · doi:10.1038/s41592-021-01153-z · Cited by 261 · see citing works
Cheroni, C., Caporale, N., & Testa, G. (2020) · Molecular Autism · doi:10.1186/s13229-020-00370-1 · Cited by 268 · see citing works
Maternal-fetal immune responses in pregnant women infected with SARS-CoV-2
Garcia‐Flores, V., Romero, R. J., Xu, Y., Theis, K. R., Arenas‐Hernandez, M., Miller, D., Peyvandipour, A., Bhatti, G., Galaz, J., Gershater, M., et al. (2022) · Nature Communications · doi:10.1038/s41467-021-27745-z · Cited by 233 · see citing works
The single-cell eQTLGen consortium
van der Wijst, M. G. P., de Vries, D., Groot, H. E., Trynka, G., Hon, C.-C., Bonder, M. J., Stegle, O., Nawijn, M. C., Idaghdour, Y., van der Harst, P., et al. (2020) · eLife · doi:10.7554/elife.52155 · Cited by 215 · see citing works
He, P., Lim, K., Sun, D., Pett, J. P., Jeng, Q., Polański, K., Dong, Z., Bolt, L., Richardson, L., Mamanova, L., et al. (2022) · Cell · doi:10.1016/j.cell.2022.11.005 · Cited by 221 · see citing works
Terekhova, M., Swain, A., Bohacova, P., Aladyeva, E., Arthur, L., Laha, A., Mogilenko, D. A., Burdess, S., Sukhov, V., Kleverov, D., et al. (2023) · Immunity · doi:10.1016/j.immuni.2023.10.013 · Cited by 199 · see citing works
Single-cell CUT&Tag analysis of chromatin modifications in differentiation and tumor progression
Wu, S. J., Furlan, S. N., Mihalas, A. B., Kaya-Okur, H. S., Feroze, A. H., Emerson, S. N., Zheng, Y., Carson, K., Cimino, P. J., Keene, C. D., et al. (2021) · Nature Biotechnology · doi:10.1038/s41587-021-00865-z · Cited by 201 · see citing works
Induction of trained immunity by influenza vaccination - impact on COVID-19
Debisarun, P. A., Gössling, K. L., Bulut, O. P., Kılıç, G., Zoodsma, M., Liu, Z., Oldenburg, M., Rüchel, N., Zhang, B., Xu, C., et al. (2021) · PLoS Pathogens · doi:10.1371/journal.ppat.1009928 · Cited by 185 · see citing works
“Next, we demultiplexed and clustered cells within each pool into four samples using souporcell, a genotype-free method [42].”
MethodFibroblast inflammatory priming determines regenerative versus fibrotic skin repair in reindeer
Sinha, S., Sparks, H. D., Labit, E., Robbins, H. N., Gowing, K., Jaffer, A., Kutluberk, E., Arora, R., Raredon, M. S. B., Cao, L., et al. (2022) · Cell · doi:10.1016/j.cell.2022.11.004 · Cited by 180 · see citing works
“Sample identify of each cell was computational demultiplexed based on natural genetic variation using souporcell (v2).(23) Briefly, souporcell analyses were performed with the number of clusters set to k = 3 and without the use of ‘common_variants’ or ‘skip_remap’ options as species-specific variants are not available…”
MethodDiscrete limbal epithelial stem cell populations mediate corneal homeostasis and wound healing
Altshuler, A., Amitai-Lange, A., Tarazi, N., Dey, S., Strinkovsky, L., Hadad-Porat, S., Bhattacharya, S., Nasser, W., Imeri, J., Ben-David, G., et al. (2021) · Cell stem cell · doi:10.1016/j.stem.2021.04.003 · Cited by 179 · see citing works
Cheng, Q., Butler, W., Zhou, Y., Zhang, H., Tang, L., Perkinson, K., Chen, X., Jiang, X. “., McCall, S. J., Inman, B. A., et al. (2022) · European Urology · doi:10.1016/j.eururo.2021.12.039 · Cited by 163 · see citing works
“[17] Heaton H, Talman AM, Knights A, et al. Souporcell: robust clustering of single-cell RNA-seq data by genotype without reference genotypes.”
MethodInfluentialLasry, A., Nadorp, B., Fornerod, M. W. J., Nicolet, D., Wu, H., Walker, C. J., Sun, Z., Witkowski, M. T., Tikhonova, A. N., Guillamot, M., et al. (2022) · Nature Cancer · doi:10.1038/s43018-022-00480-0 · Cited by 154 · see citing works
Signature of long-lived memory CD8+ T cells in acute SARS-CoV-2 infection
Adamo, S., Michler, J., Zurbuchen, Y., Cervia, C., Taeschler, P., Raeber, M. E., Sain, S. B., Nilsson, J., Moor, A., & Boyman, O. (2021) · Nature · doi:10.1038/s41586-021-04280-x · Cited by 141 · see citing works
Human SARS-CoV-2 challenge uncovers local and systemic response dynamics
Lindeboom, R. G. H., Worlock, K. B., Dratva, L. M., Yoshida, M., Scobie, D., Wagstaffe, H. R., Richardson, L., Wilbrey-Clark, A., Barnes, J. L., Kretschmer, L., et al. (2024) · Nature · doi:10.1038/s41586-024-07575-x · Cited by 143 · see citing works
“Souporcell (v.2.0) 39 was used to demultiplex each pool based on the genotype differences between the mixed samples.”
InfluentialA cell atlas foundation model for scalable search of similar human cells
Heimberg, G., Kuo, T. C.-Y., DePianto, D. J., Salem, O., Heigl, T., Diamant, N., Scalia, G., Biancalani, T., Turley, S. J., Rock, J. R., et al. (2024) · Nature · doi:10.1038/s41586-024-08411-y · Cited by 179 · see citing works
“Individual samples were genetically demultiplexing using the singularity container provided with Souporcell (v.2.0) 82 .”
InfluentialGenetic ancestry effects on the response to viral infection are pervasive but cell type specific
Randolph, H. E., Mu, Z., Fiege, J. K., Thielen, B. K., Grenier, J., Cobb, M., Hussin, J., Barreiro, L. B., Langlois, R. A., & Barreiro, L. B. (2021) · Science · doi:10.1126/science.abg0928 · Cited by 131 · see citing works
Molè, M. A., Coorens, T. H. H., Shahbazi, M. N., Weberling, A., Weatherbee, B. A. T., Gantner, C. W., Sancho‐Serra, C., Richardson, L., Drinkwater, A., Syed, N., et al. (2021) · Nature Communications · doi:10.1038/s41467-021-23758-w · Cited by 131 · see citing works
Senkowski, W., Gall-Mas, L., Falco, M. M., Li, Y., Lavikka, K., Kriegbaum, M. C., Oikkonen, J., Bulanova, D. R., Pietras, E. J., Voßgröne, K., et al. (2023) · Developmental Cell · doi:10.1016/j.devcel.2023.04.012 · Cited by 119 · see citing works
Oelen, R., de Vries, D. H., Brugge, H., Gordon, M. G., Vochteloo, M., Ye, C., Westra, H.-J., Franke, L., & van der Wijst, M. G. P. (2022) · Nature Communications · doi:10.1038/s41467-022-30893-5 · Cited by 117 · see citing works
Yu, L., Logsdon, D. M., Pinzón-Arteaga, C. A., Duan, J., Ezashi, T., Wei, Y., Orsi, A. E. R., Oura, S., Liu, L., Wang, L., et al. (2023) · Cell stem cell · doi:10.1016/j.stem.2023.08.002 · Cited by 117 · see citing works
Franken, A., Bila, M., Mechels, A., Kint, S., Dessel, J. V., Pomella, V., Vanuytven, S., Philips, G., Bricard, O., Xiong, J., et al. (2024) · Immunity · doi:10.1016/j.immuni.2024.02.007 · Cited by 117 · see citing works
A multi-omic atlas of human embryonic skeletal development
To, K., Fei, L., Pett, J. P., Roberts, K., Blain, R., Polański, K., Li, T., Yayon, N., He, P., Xu, C., et al. (2024) · Nature · doi:10.1038/s41586-024-08189-z · Cited by 98 · see citing works
Li, Q., Sharkey, A., Sheridan, M. A., Magistrati, E., Arutyunyan, A., Huhn, O., Sancho‐Serra, C., Anderson, H. E., McGovern, N. N., Esposito, L., et al. (2024) · Cell stem cell · doi:10.1016/j.stem.2023.12.013 · Cited by 120 · see citing works
A prenatal skin atlas reveals immune regulation of human skin morphogenesis
Gopee, N. H., Winheim, E., Olabi, B., Admane, C., Foster, A. R., Huang, N., Botting, R. A., Torabi, F., Sumanaweera, D., Le, A. P., et al. (2024) · Nature · doi:10.1038/s41586-024-08002-x · Cited by 108 · see citing works
“Possible maternal contamination (total of 118 cells) was identified using the souporcell pipeline (v.2.4.0) 88 as previously described 5,58 .”
InfluentialWang, M., Song, W., Chen, M., Wang, Q., Zhou, X., Xu, P., Krek, A., Yoon, Y., Ho, L., Orr, M. E., et al. (2022) · Molecular Neurodegeneration · doi:10.1186/s13024-022-00517-z · Cited by 100 · see citing works
“Genotype‐based approaches RNA reads from scRNA-seq provide a unique opportunity to infer SNVs per cell to demultiplex for individual samples [108–110], or cluster cells to trace clonal evolution [107] or genotype distributions [72].”
BackgroundMethodInfluentialYoon, Y., Kim, D., Tak, K. Y., Hwang, S., Kim, J., Sim, N. S., Cho, J.-M., Choi, D., Ji, Y., Hur, J. K., et al. (2022) · Nature Communications · doi:10.1038/s41467-022-30934-z · Cited by 101 · see citing works
“The processed fi les were demultiplexed with Souporcell 77 .”
A single-cell atlas of Plasmodium falciparum transmission through the mosquito
Real, E., Howick, V. M., Dahalan, F. A., Witmer, K., Cudini, J., Andradi-Brown, C., Blight, J., Davidson, M. S., Dogga, S. K., Reid, A. J., et al. (2021) · Nature Communications · doi:10.1038/s41467-021-23434-z · Cited by 93 · see citing works
Self-citationKhabirova, E., Jardine, L., Coorens, T. H. H., Webb, S., Treger, T. D., Engelbert, J., Porter, T., Prigmore, E., Collord, G., Piapi, A., et al. (2022) · Nature Medicine · doi:10.1038/s41591-022-01720-7 · Cited by 84 · see citing works
Cappuyns, S., Philips, G., Vandecaveye, V., Boeckx, B., Schepers, R., Brussel, T. V., Arijs, I., Mechels, A., Bassez, A., Lodi, F., et al. (2023) · Nature Communications · doi:10.1038/s41467-023-43381-1 · Cited by 92 · see citing works
Neavin, D., Senabouth, A., Arora, H., Lee, J. T. H., Ripoll-Cladellas, A., sc-eQTLGen Consortium, Franke, L., Prabhakar, S., Ye, C. J., McCarthy, D. J., et al. (2024) · Genome Biology · doi:10.1186/s13059-024-03224-8 · Cited by 90 · see citing works
“Seven of those are demultiplexing methods ( Demuxalot [6], Demuxlet [3], Dropulation [5], Freemuxlet [16], ScSplit [7], Souporcell [4], and Vireo [2]) which leverage the common genetic variation between individuals to identify cells that came from each individual and to identify heterogenic…”
BackgroundMethodSelf-citationLifelong restructuring of 3D genome architecture in cerebellar granule cells
Tan, L., Shi, J., Moghadami, S., Wright, C. P., Parasar, B., Seo, Y., Vallejo, K., Cobos, I., Duncan, L. E., Chen, R., et al. (2023) · Science · doi:10.1126/science.adh3253 · Cited by 84 · see citing works
“We found increasing interactions among certain human chromosomes—most prominently within a hub of Chr 1, 9, 11, 14, 15, 16, 17, 21, and 22, and between chromosome pairs such as Chr 2/9, 4/14, 8/11, 13/20; meanwhile, Chr 12 weakened its interactions with the hub (Fig.”
BackgroundLee, A., Kim, C., Park, S.-W., Joo, J., Choi, B., Yang, D., Jun, K., Eom, J., Lee, S., Chung, S. J., et al. (2023) · Science Advances · doi:10.1126/sciadv.abo2467 · Cited by 88 · see citing works
Evolutionary cell type mapping with single-cell genomics
Tanay, A., & Sebé-Pedrós, A. (2021) · Trends in Genetics · doi:10.1016/j.tig.2021.04.008 · Cited by 75 · see citing works
DropletQC: improved identification of empty droplets and damaged cells in single-cell RNA-seq data
Muskovic, W., & Powell, J. E. (2021) · Genome biology · doi:10.1186/s13059-021-02547-0 · Cited by 82 · see citing works
“problem, several methods have been developed to estimate and remove its contribution to gene expression [2 – 4].”
MethodGuo, C., Liu, Q., Zong, D., Zhang, W., Zuo, Z.-Q., Yu, Q., Sha, Q., Zhu, L., Gao, X., Fang, J., et al. (2022) · Cell Reports · doi:10.1016/j.celrep.2022.111606 · Cited by 86 · see citing works
Immune responses in checkpoint myocarditis across heart, blood and tumour
Blum, S. M., Zlotoff, D. A., Smith, N. P., Kernin, I., Ramesh, S., Zubiri, L., Caplin, J., Samanta, N., Martin, S., Wang, M., et al. (2024) · Nature · doi:10.1038/s41586-024-08105-5 · Cited by 96 · see citing works
“Next, the pooled scRNA-seq data were demultiplexed using Souporcell (v.2.5) 63 , and the resulting VCF files for each predicted donor were matched to the VCF files from our bulk RNA-seq data using Vireo (v.0.5.8, Python) 64 .”
InfluentialShannon, M. J., McNeill, G. L., Koksal, B., Baltayeva, J., Wächter, J., Castellana, B., Peñaherrera, M. S., Robinson, W. P., Leung, P. C. K., & Beristain, A. G. (2024) · Developmental Cell · doi:10.1016/j.devcel.2024.01.023 · Cited by 87 · see citing works
Vértesy, Á., Eichmüller, O. L., Naas, J., Novatchkova, M., Esk, C., Balmaña, M., Ladstätter, S., Bock, C., von Haeseler, A., & Knoblich, J. A. (2022) · The EMBO Journal · doi:10.15252/embj.2022111118 · Cited by 70 · see citing works
“For pooled 10 × GEX libraries, the donor cell-line of the assayed cells was determined by genotype-based demultiplexing using souporcell (Heaton et al , 2020).”
MethodEarly human lung immune cell development and its role in epithelial cell fate
Barnes, J. L., Yoshida, M., He, P., Worlock, K. B., Lindeboom, R. G. H., Suo, C., Pett, J. P., Wilbrey-Clark, A., Dann, E., Mamanova, L., et al. (2023) · Science Immunology · doi:10.1126/sciimmunol.adf9988 · Cited by 78 · see citing works
Garcia‐Flores, V., Romero, R. J., Tarca, A. L., Peyvandipour, A., Xu, Y., Galaz, J., Miller, D., Chaiworapongsa, T., Chaemsaithong, P., Berry, S. M., et al. (2024) · Science Translational Medicine · doi:10.1126/scitranslmed.adh8335 · Cited by 75 · see citing works
A village in a dish model system for population-scale hiPSC studies
Neavin, D. R., Steinmann, A. M., Farbehi, N., Chiu, H. S., Daniszewski, M., Arora, H., Bermudez, Y., Moutinho, C., Chan, C., Bax, M., et al. (2023) · Nature Communications · doi:10.1038/s41467-023-38704-1 · Cited by 72 · see citing works
“ThehiPSCvillageswerecapturedonDays 0, 1, 2, 3, 4, 5, 7, and 15 to capture the transcriptional changes that occur early during differentiation.”
Kim, S., Leem, G., Choi, J., Koh, Y., Lee, S., Nam, S., Kim, J. S., Park, C. H., Hwang, H. K., Min, K. I., et al. (2024) · Genome Medicine · doi:10.1186/s13073-024-01287-7 · Cited by 74 · see citing works
“To demultiplex the samples in a single 10X scRNA-seq library, we utilized Souporcell [22] (v2.”
MethodThe genetic architecture of cell type–specific cis regulation in maize
Marand, A. P., Jiang, L.-G., Gomezcano, F., Minow, M. A. A., Zhang, X., Mendieta, J. P., Luo, Z.-L., Bang, S., Yan, H., Meyer, C., et al. (2025) · Science · doi:10.1126/science.ads6601 · Cited by 51 · see citing works
Generating human bone marrow organoids for disease modeling and drug discovery
Olijnik, A.-A., Rodriguez-Romera, A., Wong, Z. C., Shen, Y., Reyat, J. S., Jooss, N. J., Rayes, J., Psaila, B., & Khan, A. O. (2024) · Nature Protocols · doi:10.1038/s41596-024-00971-7 · Cited by 71 · see citing works
Single-cell profiling of proteins and chromatin accessibility using PHAGE-ATAC
Fiskin, E., Lareau, C. A., Ludwig, L. S., Eraslan, G., Liu, F., Ring, A. M., Xavier, R. J., & Regev, A. (2022) · Nature Biotechnology · doi:10.1038/s41587-021-01065-5 · Cited by 68 · see citing works
Senabouth, A., Daniszewski, M., Lidgerwood, G. E., Liang, H. H., Hernández, D., Mirzaei, M., Keenan, S. N., Zhang, R., Han, X., Neavin, D. R., et al. (2022) · Nature Communications · doi:10.1038/s41467-022-31707-4 · Cited by 72 · see citing works
Bailey, J. I., Puritz, C. H., Senkow, K. J., Markov, N. S., Diaz, E., Jonasson, E., Yu, Z., Swaminathan, S., Lu, Z., Fenske, S. W., et al. (2024) · Nature Immunology · doi:10.1038/s41590-024-01975-x · Cited by 70 · see citing works
“41) for that library and two libraries that were prepared from these samples separately.”
BackgroundMMR vaccination induces trained immunity via functional and metabolic reprogramming of γδ T cells
Röring, R. J., Debisarun, P. A., Botey-Bataller, J., Suen, T. K., Bulut, O. P., Kılıç, G., Koeken, V. A. C. M., Sarlea, A. S., Bahrar, H., Dijkstra, H., et al. (2024) · Journal of Clinical Investigation · doi:10.1172/jci170848 · Cited by 66 · see citing works
SCITO-seq: single-cell combinatorial indexed cytometry sequencing
Hwang, B., Lee, D. S., Tamaki, W., Sun, Y., Ogorodnikov, A., Hartoularos, G. C., Winters, A. F., Yeung, B. Z., Nazor, K. L., Song, Y. S., et al. (2021) · Nature Methods · doi:10.1038/s41592-021-01222-3 · Cited by 63 · see citing works
The effect of background noise and its removal on the analysis of single-cell expression data
Janssen, P., Kliesmete, Z., Vieth, B., Adiconis, X., Simmons, S., Marshall, J. L., Aubin, C., Heyn, H., Levin, J. Z., Enard, W., et al. (2023) · Genome biology · doi:10.1186/s13059-023-02978-x · Cited by 64 · see citing works
Moerkens, R., Mooiweer, J., Ramírez-Sánchez, A. D., Oelen, R., Franke, L., Wijmenga, C., Barrett, R. J., Jonkers, I. H., & Withoff, S. (2024) · Cell Reports · doi:10.1016/j.celrep.2024.114247 · Cited by 64 · see citing works
MacDonald, L., Elmesmari, A., Somma, D., Frew, J., Mario, C. D., Madhu, R., Paoletti, A., Simakou, T., Hardy, O. M., Tolusso, B., et al. (2024) · Immunity · doi:10.1016/j.immuni.2024.11.004 · Cited by 65 · see citing works
Single-cell immune profiling reveals distinct immune response in asymptomatic COVID-19 patients
Zhao, X.-N., You, Y., Cui, X., Gao, H., Wang, G., Zhang, S., Lin, Y., Duan, L., Zhu, K.-L., Wang, Y., et al. (2021) · Signal Transduction and Targeted Therapy · doi:10.1038/s41392-021-00753-7 · Cited by 62 · see citing works
Johansen, M. D., Mahbub, R. M., Idrees, S., Nguyen, D. H., Miemczyk, S., Pathinayake, P. S., Nichol, K., Hansbro, N. G., Gearing, L. J., Hertzog, P. J., et al. (2022) · American Journal of Respiratory and Critical Care Medicine · doi:10.1164/rccm.202108-1901oc · Cited by 61 · see citing works
A single-cell view on host immune transcriptional response to in vivo BCG-induced trained immunity
Li, W., Moorlag, S. J. C. F. M., Koeken, V. A. C. M., Röring, R. J., de Bree, L. C. J., Mourits, V. P., Gupta, M. K., Zhang, B., Fu, J., Zhang, Z., et al. (2023) · Cell Reports · doi:10.1016/j.celrep.2023.112487 · Cited by 56 · see citing works
Age-specific nasal epithelial responses to SARS-CoV-2 infection
Woodall, M. N. J., Cujba, A.-M., Worlock, K. B., Case, K., Masonou, T., Yoshida, M., Polański, K., Huang, N., Lindeboom, R. G. H., Mamanova, L., et al. (2024) · Nature Microbiology · doi:10.1038/s41564-024-01658-1 · Cited by 58 · see citing works
“Cells that were identified as heterotypic doublets by Souporcell were discarded.”
Boháčová, P., Terekhova, M., Tsurinov, P., Mullins, R. D. Z., Husarcikova, K., Shchukina, I., Antonova, A. U., Echalar, B., Kössl, J., Saidu, A., et al. (2024) · Immunity · doi:10.1016/j.immuni.2024.08.019 · Cited by 62 · see citing works
Cheng, J., Liao, J., Shao, X., Lu, X., & Fan, X. (2021) · Advanced Science · doi:10.1002/advs.202101229 · Cited by 56 · see citing works
Guo, C., Wu, M., Huang, B., Zhao, R., Jin, L., Fu, B., Wang, P., Wang, D., Zheng, M., Fang, J., et al. (2022) · Genome Medicine · doi:10.1186/s13073-022-01049-3 · Cited by 54 · see citing works
“Heaton H, Talman AM, Knights A, Imaz M, Gaffney DJ, Durbin R, et al. Souporcell: robust clustering of single-cell RNA-seq data by genotype without reference genotypes.”
MethodPenter, L., Liu, Y., Wolff, J. O., Yang, L., Taing, L., Jhaveri, A., Southard, J., Patel, M., Cullen, N. M., Pfaff, K. L., et al. (2023) · Blood · doi:10.1182/blood.2022018246 · Cited by 53 · see citing works
“33 By applying souporcell 34 (Figure 2D, left) on 62 673 myeloid cells from the post-HSCT baseline time point from 8 patients (median 296 days post-HSCT, range 151-1290), we identi fi ed 1422 cells (2.3”
InfluentialMethods and Insights from Single-Cell Expression Quantitative Trait Loci
Kang, J. B., Raveane, A., Nathan, A., Soranzo, N., & Raychaudhuri, S. (2023) · Annual Review of Genomics and Human Genetics · doi:10.1146/annurev-genom-101422-100437 · Cited by 53 · see citing works
“Since then, advancements in the throughput of single-cell technologies and sample multiplexing (60, 70, 136) have enabled the sequencing of larger cohorts amenable to eQTL analyses.”
BackgroundLatent human herpesvirus 6 is reactivated in CAR T cells
Lareau, C. A., Yin, Y., Maurer, K., Sandor, K. D., Dániel, B., Yagnik, G., Peña, J., Crawford, J. C., Spanjaart, A. M., Gutierrez, J. C., et al. (2023) · Nature · doi:10.1038/s41586-023-06704-2 · Cited by 54 · see citing works
Cong, J., Wu, J., Fang, Y., Wang, J., Kong, X., Wang, L., & Duan, Z. (2024) · Environment International · doi:10.1016/j.envint.2024.108744 · Cited by 55 · see citing works
Rodríguez‐Ubreva, J., Arutyunyan, A., Bonder, M. J., Molina, L. D. P., Clark, S. J., de la Calle‐Fabregat, C., García‐Alonso, L., Handfield, L., Ciudad, L., Andrés‐León, E., et al. (2022) · Nature Communications · doi:10.1038/s41467-022-29450-x · Cited by 52 · see citing works
Recruitment of highly cytotoxic CD8+ T cell receptors in mild SARS-CoV-2 infection
Wagner, K. I., Mateyka, L. M., Jarosch, S., Grass, V., Weber, S., Schober, K., Hammel, M., Burrell, T., Kalali, B. N., Poppert, H., et al. (2021) · Cell Reports · doi:10.1016/j.celrep.2021.110214 · Cited by 45 · see citing works
“…et al., 2020 https://doi.org/10.1093/bioinformatics/ btaa611 scSplit 1.0 Xu et al., 2019 https://doi.org/10.1186/ s13059-019-1852-7 Souporcell 2.0 Heaton et al., 2020 https://doi.org/10.1038/ s41592-020-0820-1 uType software Invitrogen/ThermoFisher N/A Other Äkta pureSuperdeso 200 10/300GL GE…”
InfluentialRobust temporal map of human in vitro myelopoiesis using single-cell genomics
Alsinet, C., Primo, M. N., Lorenzi, V., Bello, E., Kelava, I., Jones, C. P., Vilarrasa‐Blasi, R., Sancho‐Serra, C., Knights, A. J., Park, J.-E., et al. (2022) · Nature Communications · doi:10.1038/s41467-022-30557-4 · Cited by 44 · see citing works
Self-citationHuman conjunctiva organoids to study ocular surface homeostasis and disease
Bannier-Hélaouët, M., Korving, J., Ma, Z., Begthel, H., Giladi, A., Lamers, M. M., van de Wetering, W. J., Yawata, N., Yawata, M., LaPointe, V. L. S., et al. (2024) · Cell stem cell · doi:10.1016/j.stem.2023.12.008 · Cited by 46 · see citing works
A single cell atlas of sexual development in Plasmodium falciparum
Dogga, S. K., Rop, J., Cudini, J., Farr, E., Dara, A., Ouologuem, D., DJIMDÉ, A. A., Talman, A. M., & Lawniczak, M. (2024) · Science · doi:10.1126/science.adj4088 · Cited by 43 · see citing works
Self-citationGeyer, M. B., DeWolf, S., Mi, X., et al. (2025) · Blood · doi:10.1182/blood.2025029532 · Cited by 43 · see citing works
Landscape and age dynamics of immune cells in the Egyptian rousette bat
Friedrichs, V., Toussaint, C., Schäfer, A., Rissmann, M., Dietrich, O., Mettenleiter, T. C., Pei, G., Balkema‐Buschmann, A., Saliba, A., & Dorhoi, A. (2022) · Cell Reports · doi:10.1016/j.celrep.2022.111305 · Cited by 41 · see citing works
Sample-multiplexing approaches for single-cell sequencing
Zhang, Y., Xu, S., Wen, Z., Gao, J., Li, S., Weissman, S. M., & Pan, X. (2022) · Cellular and Molecular Life Sciences · doi:10.1007/s00018-022-04482-0 · Cited by 40 · see citing works
“More recently, the computational tool, Souporcell, was designed for clustering scRNA-seq data by genotype without reference genotypes [31].”
MethodMapping the developing human immune system across organsPreprint
Suo, C., Dann, E., Goh, I., et al. (2022) · bioRxiv (Cold Spring Harbor Laboratory) · doi:10.1101/2022.01.17.476665 · Cited by 40 · see citing works
Heimberg, G., Kuo, T., DePianto, D., Heigl, T., Diamant, N., Salem, O., Scalia, G., Biancalani, T., Turley, S., Rock, J., et al. (2023) · doi:10.1101/2023.07.18.549537 · Cited by 40 · see citing works
BCG vaccination induces innate immune memory in γδ T cells in humans
Suen, T. K., Moorlag, S. J. C. F. M., Li, W., de Bree, L. C. J., Koeken, V. A. C. M., Mourits, V. P., Dijkstra, H., Lemmers, H., Bhat, J., Xu, C., et al. (2023) · Journal of Leukocyte Biology · doi:10.1093/jleuko/qiad103 · Cited by 41 · see citing works
Integrative genotyping of cancer and immune phenotypes by long-read sequencing
Penter, L., Borji, M., Nagler, A., Lyu, H., Lu, W. S., Cieri, N., Maurer, K., Oliveira, G., Al’Khafaji, A. M., Garimella, K. V., et al. (2024) · Nature Communications · doi:10.1038/s41467-023-44137-7 · Cited by 40 · see citing works
“Distinguishing normal and malignant hematopoiesis using expressed donor-and recipient-speci fi c single nucleotide polymorphisms (SNPs) ( souporcell ) 25 revealed two additional AML-derived ECs, mega-karyopoiesis and erythropoiesis, to be almost entirely recipient-derived in 6 of 8 analyzed cases…”
BackgroundGibson, A., Ram, R., Gangula, R., Li, Y., Mukherjee, E. M., Palubinsky, A. M., Campbell, C. N., Thorne, M., Konvinse, K. C., Choshi, P., et al. (2024) · Nature Communications · doi:10.1038/s41467-024-52990-3 · Cited by 41 · see citing works
Single cell profiling of COVID-19 patients: an international data resource from multiple tissuesPreprint
Consortia, C. Z. I. S.-C. C.-1., Ballestar, E., Färber, D. L., Glover, S. C., Horwitz, B. H., Meyer, K. B., Nikolić, M., Ordovás-Montañés, J., Sims, P. A., Shalek, A. K., et al. (2020) · medRxiv · doi:10.1101/2020.11.20.20227355 · Cited by 39 · see citing works
MALAT1 expression indicates cell quality in single-cell RNA sequencing dataPreprint
Clarke, Z. A., & Bader, G. D. (2024) · bioRxiv (Cold Spring Harbor Laboratory) · doi:10.1101/2024.07.14.603469 · Cited by 41 · see citing works
“For droplet-based methods, cells may release RNA transcripts into the solution of the experiment during cell processing creating ambient RNA which makes it difficult to tell which transcripts come from which cells (Young and Behjati 2018; Heaton et al. 2020).”
Krakow, E. F., Brault, M., Summers, C. N., et al. (2024) · Blood · doi:10.1182/blood.2024024105 · Cited by 38 · see citing works
Hendriks, D. F. G., Artegiani, B., Margaritis, T., Zoutendijk, I., Lopes, S. M. C. D. S., & Clevers, H. C. (2024) · Nature Communications · doi:10.1038/s41467-024-48550-4 · Cited by 39 · see citing works
Modeling human embryo implantation in vitro
Molè, M. A., Elderkin, S., Zorzan, I., et al. (2025) · Cell · doi:10.1016/j.cell.2025.10.027 · Cited by 38 · see citing works
Weber, L. M., Hippen, A. A., Hickey, P. F., Berrett, K. C., Gertz, J., Doherty, J. A., Greene, C. S., & Hicks, S. C. (2021) · GigaScience · doi:10.1093/gigascience/giab062 · Cited by 37 · see citing works
Murrow, L. M., Weber, R. J., Caruso, J. A., McGinnis, C. S., Phong, K. T., Gascard, P. D., Rabadam, G., Borowsky, A. D., Desai, T. A., Thomson, M., et al. (2022) · Cell Systems · doi:10.1016/j.cels.2022.06.005 · Cited by 35 · see citing works
“We used MULTI-seq lipid-based barcoding and in silico genotyping for sample multiplexing to minimize technical variability between samples (Figure S1C; Table S2; STAR Methods; Heaton et al., 2020; McGinnis et al., 2019b).”
SMARCB1 loss activates patient-specific distal oncogenic enhancers in malignant rhabdoid tumors
Liu, N.-Q., Paassen, I., Custers, L., Zeller, P., Teunissen, H., Ayyildiz, D., He, J., Buhl, J. L., Hoving, E., van Oudenaarden, A., et al. (2023) · Nature Communications · doi:10.1038/s41467-023-43498-3 · Cited by 36 · see citing works
“Reads were mapped to GRCh38 and cell genotypes were annotated by souporcell 75 .”
MethodCaporale, N., Castaldi, D., Rigoli, M. T., Cheroni, C., Valenti, A., Stucchi, S., Lessi, M., Bulgheresi, D., Trattaro, S., Pezzali, M., et al. (2024) · Nature Methods · doi:10.1038/s41592-024-02555-5 · Cited by 35 · see citing works
“1b) were demultiplexed using three different state-of-the-art algorithms available for genetic demultiplexing (demuxlet 27 , souporcell 28 , Vireo 29 ).”
Simple oligonucleotide-based multiplexing of single-cell chromatin accessibility
Wang, K., Xiao, Z., Yan, Y., Ye, R., Hu, M. S., Bai, S., Sei, E., Qiao, Y., Chen, H., Lim, B., et al. (2021) · Molecular Cell · doi:10.1016/j.molcel.2021.09.026 · Cited by 35 · see citing works
deMULTIplex2: robust sample demultiplexing for scRNA-seq
Zhu, Q., Conrad, D. N., & Gartner, Z. J. (2024) · Genome biology · doi:10.1186/s13059-024-03177-y · Cited by 32 · see citing works
“For multi-donor datasets, we ran SNP-based sample classification using vireo [11] or souporcell [12] when such genotype-based classification were not provided.”
- Authors
- Heaton, H., Talman, A. M., Knights, A., Imaz, M., Gaffney, D. J., Durbin, R., Hemberg, M., & Lawniczak, M. K. N.
- Affiliations
- Wellcome Sanger Institute (GB)
- Centre National de la Recherche Scientifique (FR)
- University of Cambridge (GB)
- Corresponding author
- Heaton, Haynes
- Hemberg, Martin
- Lawniczak, Mara K. N.
- ORCID
- Venue
- Nature Methods
- Published
- 2020-05-04
- Volume
- 17
- Issue
- 6
- Pages
- 615-620
- ISSN
- 1548-7091 (ISSN-L)
- 1548-7105
- PMC
- PMC7617080
- OpenAlex
- W3023290106
- Cited by
- 506
- References
- 25
- Citations by year
- 2026: 95 · 2025: 88 · 2024: 126 · 2023: 77 · 2022: 69 · 2021: 38 · 2020: 11 · 2018: 1
- FWCI
- 13.53
- Type
- journal article
- Language
- en
- Open access
- green
- www.ncbi.nlm.nih.gov
- Locations
- Nature Methods (publishedVersion)
- PubMed (publishedVersion)
- PubMed Central (submittedVersion)
- HAL (Le Centre pour la Communication Scientifique Directe) (submittedVersion)
- Apollo (University of Cambridge) (submittedVersion)
- Apollo
- Funders
- Wellcome Trust (206194/Z/17/Z)
- British Heart Foundation (RG/13/13/30194; RG/18/13/33946)
- Topics
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gene expression and cancer classification
- Keywords
- single-cell RNA sequencing, genotype clustering, cell lysis, batch effects
- MeSH terms
- RNA-Seq, Algorithms, Base Sequence, Cell Line, Genotype, Humans, RNA, Sensitivity and Specificity, Software, Cluster Analysis, Polymorphism, Single Nucleotide, Single-Cell Analysis
About these checks
Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.
Claims read by deepseek-v4.1-flash under definition version 3 (extract-v3) on 2026-10-13.
Classified by jev-1.13-20260917 (gold-guideline-v0.2) on 2026-10-13.
Cited sources are checked by gpt-6-sol.
Statistics and arithmetic are recomputed, not judged.
Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.
What counts as a claim ›How we measure confidence ›How we check this ›
Are you an author, or do you see a problem with this page? Request a review