<p>Indeed there's no satisfactory doublet/multiplet identification methods as of 2020.</p>
<p>The most reliable way to remove doublets is still manual removal of doublet subclusters guided by marker genes.</p>
<p>Not only trajectory analysis can be impacted by doublets, but clustering and embedding can also be affected for obvious reasons.</p>

CC BY-SA 4.0
2020-12-30T13:50:40.023Z
2020-12-30T13:50:40.023Z
Code42
<p>I have attached a histogram plot of the number of genes per cell in a single cell RNA seq data set of lung endothelial cells. I do not find a bimodal or multimodal distribution of the number of genes per cell within any particular condition. Doublets/multiplets should have an aberrantly high number of gene expression and if many, can lead to a bi or multimodal histogram plot of the number of genes detected per cell.  Can I conclude that I have an insignificant presence of doublets or multiplets in my data set from this plot and hence continue with downstream analyses?</p>
<p>I am reluctant to use automatic doublet-finding algorithms because they have a high possibility to remove cells in a transitional differentiating phase. I also think that the default threshold of 2500 genes per cell given by Seurat might sometimes lead to the loss of good cells depending on the cell type.</p>
<p><a href="https://i.sstatic.net/olCSB.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/olCSB.png" alt="nGene per cell" /></a></p>
<p>Apart from trajectory analyses, are there any other analyses on single cell RNA seq data that are sensitve to doublets/multiplets?</p>
<p>Thanks in advance for your help.</p>

CC BY-SA 4.0
2021-05-29T21:04:18.867Z
You should not expect a multi-modal distribution. See this previous discussion: https://bioinformatics.stackexchange.com/questions/3211/determining-doublets-in-single-cell-rna-seq
burger
<rna-seq><scrnaseq><seurat><single-cell>
What is the best way to address the question of doublets and multiplets in a single cell RNA seq data set?
# What is the best way to address the question of doublets and multiplets in a single cell RNA seq data set?

(Asked by: username_0 on 2020-08-04)

I have attached a histogram plot of the number of genes per cell in a single cell RNA seq data set of lung endothelial cells. I do not find a bimodal or multimodal distribution of the number of genes per cell within any particular condition. Doublets/multiplets should have an aberrantly high number of gene expression and if many, can lead to a bi or multimodal histogram plot of the number of genes detected per cell. Can I conclude that I have an insignificant presence of doublets or multiplets in my data set from this plot and hence continue with downstream analyses?

I am reluctant to use automatic doublet-finding algorithms because they have a high possibility to remove cells in a transitional differentiating phase. I also think that the default threshold of 2500 genes per cell given by Seurat might sometimes lead to the loss of good cells depending on the cell type.

[![nGene per cell](https://i.sstatic.net/olCSB.png)](https://i.sstatic.net/olCSB.png)

Apart from trajectory analyses, are there any other analyses on single cell RNA seq data that are sensitve to doublets/multiplets?

Thanks in advance for your help.

**Tags:** rna-seq, scrnaseq, seurat, single-cell
> **Comments:**
> **username_1**: You should not expect a multi-modal distribution. See this previous discussion: https://bioinformatics.stackexchange.com/questions/3211/determining-doublets-in-single-cell-rna-seq  

---

## 1 answers

### Answer by username_2 on 2020-12-30. Score: 2

Indeed there's no satisfactory doublet/multiplet identification methods as of 2020.

The most reliable way to remove doublets is still manual removal of doublet subclusters guided by marker genes.

Not only trajectory analysis can be impacted by doublets, but clustering and embedding can also be affected for obvious reasons.


---

