Mitochondrial genes in single cell nuclear RNAseq data
2
5
Entering edit mode
5.8 years ago
dalhoomist ▴ 50

I have done several single cell nuclear RNAseq experiments on human brain tissue. I followed this protocol outlined in Krishnaswami et al https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4941947/. The percentage of viable (intact) cells out of the total nuclei suspension is less than 5% - and these cells can be excluded in the analysis. The nuclear specific transcripts are enriched (for example MALAT-1). However, some mitochondrial specific genes like MT-ND1, MT-ND2, MT-ND3, MT-CO1, MT-CO2, and MT-CO3 are still detected at high levels (And often differentially expressed between clusters). Given that the RNA sequenced is nuclear RNA, how can I explain the presence of mitochondria-encoded RNA in my data?

RNA-Seq gene mitochondria • 7.7k views
ADD COMMENT
0
Entering edit mode

Since the explanation in this case is likely going to be related to experimental protocol/technique, you may want to post this at SeqAnswers and biology SE.

Informatically you could ignore the counts on mito genes if they are not of interest. But their presence seems to indicate that your nuclear preps were likely not clean or perhaps mitochondrial RNA's get transferred to nuclei? It also looks like mitochondrial genes are found in nuclei so that may also be another explanation?

ADD REPLY
0
Entering edit mode

Hi, I am facing the same question in a single nuclei dataset. The same genes that you are mentioning and they are also differential between my treatment groups. Did you ever get an idea on what is going on here?

ADD REPLY
0
Entering edit mode

Variable mitochondrial genes are pretty typical, and as genomax said, this could be indicative of contamination (either through intact cells or ambient RNA). What % of reads are in mitochondrial genes? If it's only a few percent, you can likely ignore them.

ADD REPLY
2
Entering edit mode
4.0 years ago
ldaveyl ▴ 20

I've read in some papers that use single nucleus datasets that mitochondrial genes are ignored for downstream analysis.

https://www.nature.com/articles/s41586-019-1404-z#Sec8

"Mitochondrial RNA genes were filtered out as well to exclude transcripts originating from outside the nucleus and avoid biases introduced by the isolation of nuclei and ultracentrifugation."

ADD COMMENT
0
Entering edit mode
11 months ago
Prakki Rama ★ 2.7k

Quoting from 10X Genomics documentation page:

Even in 10x assays with single nuclei as input (for example, Multiome assay), an elevated level of mitochondrial genes would mean that mitochondrial RNA could remain “stuck” to nuclear membranes or may get partitioned into GEMs and be detected by the gene expression biochemistry.

ADD COMMENT

Login before adding your answer.

Traffic: 1160 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6