Question: Seurat scRNAseq analysis: what is the mito gene list
0
gravatar for wenbinm
10 months ago by
wenbinm10
USA
wenbinm10 wrote:

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!

scrnaseq rna-seq • 1.5k views
ADD COMMENTlink modified 9 months ago by rogangrant0 • written 10 months ago by wenbinm10

Mouse mt gene ids will start with "mt-"

ADD REPLYlink written 5 months ago by arup1.7k
2
gravatar for igor
10 months ago by
igor8.3k
United States
igor8.3k wrote:

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

ADD COMMENTlink written 10 months ago by igor8.3k

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'

ADD REPLYlink written 10 months ago by wenbinm10
0
gravatar for rogangrant
9 months ago by
rogangrant0
rogangrant0 wrote:

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.

ADD COMMENTlink modified 9 months ago • written 9 months ago by rogangrant0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1166 users visited in the last hour