This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Seurat scRNAseq analysis: what is the mito gene list

Hi there, I am trying to analyze 10X genomics output with Seurat. After reading in data using "Read10X" function and create an Seurat object, I am stuck at this step:

> mito.genes <- grep(pattern = "^MT-", x = rownames(x = my_data@data), value = TRUE)

For the data I have, row names of the expression matrix are gene names ("Sox17" "Mrpl15" "Lypla1" "Tcea1" "Rgs20" etc) or ensembl ID, where no genes begin with "MT-".

Does anyone know the mitochondria gene list Seurat uses for quality control? Or I should generate a mitochondria gene list on my own? I am working on mouse samples.

Many thanks for your help!

rna-seq scrnaseq
• 17,765 views • link
updated by Lauren • Lauren • written by wenbinm • wenbinm

Mouse mt gene ids will start with "mt-"

• 0 views • link
written by Arup Ghosh Arup Ghosh

is rat mt gene ids also "mt-"?

• 0 views • link
written by Lauren • Lauren

I was told it is ^Mt for rat

• 0 views • link
written by Lauren • Lauren

2 answers

Looks like you have a mouse genome. Instead of ^MT-, you can look for ^Mt- or just add ignore.case = TRUE to the grep.

• 1 views • link
written by igor igor

It works. Thank you! It seems that Seurat is always using these 13 genes: 'Nd1','Nd2','Co1','Co2','Atp8','Atp6','Co3','Nd3','Nd4l','Nd4','Nd5','Nd6','Cytb'

• 1 views • link
written by wenbinm • wenbinm

At least in my case, mouse mitochondrial genes are flagged with ^mt-, rather than ^Mt-. But either will still work if you set ignore.case = TRUE as per igor's suggestion.

• 0 views • link
written by rogangrant • rogangrant

Log in to answer this question.

More posts like this