<p>Mitochondrial gene filtering is done on a case-by-case basis, which is why there's a bit of discussion about choosing appropriate thresholds in the <a href="https://satijalab.org/seurat/articles/pbmc3k_tutorial.html" rel="nofollow noreferrer">Seurat PBMC3k vignette</a>.</p>
<p>In your specific case, there are a lot of cells with very high mitochondrial gene content. It would be a good idea to work out why that is before doing any more filtering.</p>
<p>However, you're asking specifically about <code>sctransform</code>, and <a href="https://satijalab.org/seurat/articles/sctransform_vignette.html" rel="nofollow noreferrer">the SCtransform vignette</a> does account for mitochondrial gene content in the transformation:</p>
<blockquote>
<p>During normalization, we can also remove confounding sources of variation, for example, mitochondrial mapping percentage</p>
</blockquote>
<pre><code># store mitochondrial percentage in object meta data
pbmc &lt;- PercentageFeatureSet(pbmc, pattern = &quot;^MT-&quot;, col.name = &quot;percent.mt&quot;)

# run sctransform
pbmc &lt;- SCTransform(pbmc, vars.to.regress = &quot;percent.mt&quot;, verbose = FALSE)
</code></pre>
<p>There's no filtering in that vignette, but there's also no mention of how clean the data was initially. It's still a good idea to look at your data for outliers, try to understand why they are outliers, and have a think about whether or not to exclude them.</p>

CC BY-SA 4.0
2021-08-06T23:55:55.733Z
2021-08-06T23:55:55.733Z
gringer
<p>I have a Seurat object which has a high expression of <code>mitochondrial genes</code></p>
<pre><code>&gt; pbmc1
An object of class Seurat 
36601 features across 18338 samples within 1 assay 
Active assay: RNA (36601 features, 0 variable features)
&gt; 
</code></pre>
<p><a href="https://i.sstatic.net/cgUpB.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/cgUpB.png" alt="enter image description here" /></a></p>
<p>I want to use <code>SCTransform in Seurat</code>, I don't know if still I should define the mitochondrial percentage or not, if not, how SCTransform function knows which percentage suits for my data?</p>
<p>I used SCTransform without clipping <code>MT genes</code> and I noticed I have the same number of cells</p>
<pre><code>&gt; pbmc_SCTransform &lt;- SCTransform(pbmc, method = &quot;glmGamPoi&quot;, vars.to.regress = &quot;percent.mt&quot;, verbose = FALSE)
&gt; pbmc_SCTransform
An object of class Seurat 
58629 features across 18338 samples within 2 assays 
Active assay: SCT (22028 features, 3000 variable features)
 1 other assay present: RNA
</code></pre>
<p>Is it fine or I should clip <code>MT genes</code> before running SCTransform function?</p>
<p>I am not too sure about the input composition but from 10X websummary two of my samples likely suffer from poor sample prep or poor beginning cell health. The gene expression median genes per cell are 41 and 50 and differential expression show only mitochondrial (MT) and MALAT1 genes. This is indicative of subpar lysis conditions.</p>

CC BY-SA 4.0
2021-08-06T14:42:22.180Z
Frustrating, why down voting? Do you know the answer? You believe in Seurat vignette  this has been already mentioned? No
Mahta Mira
<scrnaseq><seurat><10x-genomics>
sctransform - mitochondrial expression filtering / transformation
# sctransform - mitochondrial expression filtering / transformation

(Asked by: username_0 on 2021-08-06)

I have a Seurat object which has a high expression of mitochondrial genes
> pbmc1
An object of class Seurat 
36601 features across 18338 samples within 1 assay 
Active assay: RNA (36601 features, 0 variable features)
> 


I want to use SCTransform in Seurat, I don't know if still I should define the mitochondrial percentage or not, if not, how SCTransform function knows which percentage suits for my data?
I used SCTransform without clipping MT genes and I noticed I have the same number of cells
> pbmc_SCTransform <- SCTransform(pbmc, method = "glmGamPoi", vars.to.regress = "percent.mt", verbose = FALSE)
> pbmc_SCTransform
An object of class Seurat 
58629 features across 18338 samples within 2 assays 
Active assay: SCT (22028 features, 3000 variable features)
 1 other assay present: RNA

Is it fine or I should clip MT genes before running SCTransform function?
I am not too sure about the input composition but from 10X websummary two of my samples likely suffer from poor sample prep or poor beginning cell health. The gene expression median genes per cell are 41 and 50 and differential expression show only mitochondrial (MT) and MALAT1 genes. This is indicative of subpar lysis conditions.

**Tags:** scrnaseq, seurat, 10x-genomics
> **Comments:**
> **username_0**: Frustrating, why down voting? Do you know the answer? You believe in Seurat vignette  this has been already mentioned? No  

---

## 1 answers

### Answer by username_1 on 2021-08-06. Score: 1

Mitochondrial gene filtering is done on a case-by-case basis, which is why there's a bit of discussion about choosing appropriate thresholds in the Seurat PBMC3k vignette.
In your specific case, there are a lot of cells with very high mitochondrial gene content. It would be a good idea to work out why that is before doing any more filtering.
However, you're asking specifically about sctransform, and the SCtransform vignette does account for mitochondrial gene content in the transformation:

During normalization, we can also remove confounding sources of variation, for example, mitochondrial mapping percentage

# store mitochondrial percentage in object meta data
pbmc <- PercentageFeatureSet(pbmc, pattern = "^MT-", col.name = "percent.mt")

# run sctransform
pbmc <- SCTransform(pbmc, vars.to.regress = "percent.mt", verbose = FALSE)

There's no filtering in that vignette, but there's also no mention of how clean the data was initially. It's still a good idea to look at your data for outliers, try to understand why they are outliers, and have a think about whether or not to exclude them.


---

