mitochondrial RNA percentage in bulk RNAseq data
1
1
Entering edit mode
sohail.mbio ▴ 10
@335cd804
Last seen 13 months ago
United States

Hi,

I am analysing bulk RNAseq data of 86 samples by DESeq2. I really liked one feature of Seurat in scRNAseq analysis "mitochondrial content percentage" . Is there a simple way to extract "mt percentage" information from deseq2 object. Like

plotCounts(dds, gene="MT-", intgroup="Group")

FYI: My deseq2 object have gene symbols, so all mitochondrial genes have MT- as suffix.

thanks

rnaseqGene RNASeqR DESeq2 RNAseq • 2.5k views
ADD COMMENT
3
Entering edit mode
@mikelove
Last seen 21 hours ago
United States

We don't have this functionality, but you can just do:

cts <- counts(dds, normalized=TRUE)[idx,]

Then plot as you like, e.g. you can do colSums() of this, etc.

idx is just the lines matching whatever pattern you want, see grep or grepl which you can run on rownames(dds). The dds object can be subset with numeric or logical indices.

ADD COMMENT

Login before adding your answer.

Traffic: 786 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