<p>The table pictured here <em>is</em> the result of a differential gene expression (DGE) test. In the table, each row is a gene (labeled with the gene ID and name) and each column (besides the first two) are the results of DGE tests. Each cluster has two sub-columns: L2FC (log2 fold-change) which tells you the average relative expression changes between cells of the cluster vs. others; and p-value, which is how significant the change is.</p>
<p>You can make this table yourself by performing DGE tests for all clusters. You've already done this yourself in your Rpubs script! To make it easier though, so you don't have to manually test each cluster, use the <a href="https://www.rdocumentation.org/packages/Seurat/versions/1.4.0/topics/FindAllMarkers" rel="nofollow noreferrer">FindAllMarkers</a> function. It will compare each cluster to all other clusters and report the L2FCs and p-values.</p>
<p>However, I recommend using the p_val_adj (adjusted p-value) instead of the p-value. This is due to a problem called multiple testing: you are performing many, many tests (each gene in each cluster is a separate statistical test), and we expect 5% of differentially expressed genes (DEGs) to be DE by random chance. Over something like 100,000 tests, that means you'd expect 5,000 of those to be false positives (which is huge). So p-value is adjusted to account for multiple tests.</p>
<p>Once you filter the DEG table (p-adj&lt;0.05), you can do other downstream analyses to describe what has changed between the conditions. You already have a good start in your Rpub script with the violin plots. Cluster 6 is clearly distinct in terms of those genes' expression levels. You could try to figure out what's special with that cluster. Good luck!</p>

CC BY-SA 4.0
2022-11-29T15:26:13.863Z
With my violin plots, does that mean that those two genes , c("MS4A1", "CD79A"), are found within multiple clusters?
Antonio
CC BY-SA 4.0
2022-11-30T16:21:52.647Z
@Antonio Each violin plot shows the expression of that gene in all cells (the black dots). As you see, most values are 0 (which is very common for scRNA-seq), along with some random low-level noise. The important bits to look at is anything above the noise level (in this case, ~0.5). Any violin with color is usually somewhat significant. For the MS4A1 gene, you can see the only clear violin is in cluster 6, while CD79A is in 6 and 20. If you look at your UMAP, cluster 6 is off to the side by itself. This may be a gene that is cell-type specific or only expressed in certain conditions.
Sean Robertson
CC BY-SA 4.0
2022-11-29T02:28:41.133Z
2022-11-29T02:28:41.133Z
Sean Robertson
<p>I want to perform Differential Gene Expression Analysis. I recently completed RNA Sequence Analysis using the Seraut Tutorial up to <em>this</em> point: <a href="https://satijalab.org/seurat/articles/pbmc3k_tutorial.html#finding-differentially-expressed-features-cluster-biomarkers-" rel="nofollow noreferrer">Finding Differentially Expressed Features (Cluster Biomarkers)</a>.</p>
<p>My project can be found at <em>this</em> <a href="https://rpubs.com/ap3340/SingleCellHumanPBMC" rel="nofollow noreferrer">Rpubs Link</a> which contains the link to the dataset on the 10x genomics website.</p>
<p><strong>My question is... is the following table derived from performing the steps found in the first link I provided?</strong>
<a href="https://i.sstatic.net/sr9p0.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/sr9p0.png" alt="enter image description here" /></a></p>
<p><strong>If so, once I derive that table ( ID, Name, L2FC, p-value, etc.), does that mean I can perform Differential Gene Expression Analysis?</strong></p>
<p>I am a complete newbie and am teaching myself all of this. Thanks, y'all.</p>
<p><strong>UPDATE:</strong>
Yes, it did provide me that table. With this table, am I able to perform gene expression?</p>

<r><rna-seq><gene-expression>
What are the conditions to perform Gene Expression Analysis on data you've performed RNA Sequence Analysis on?
# What are the conditions to perform Gene Expression Analysis on data you've performed RNA Sequence Analysis on?

(Asked by: username_0 on 2022-11-28)

I want to perform Differential Gene Expression Analysis. I recently completed RNA Sequence Analysis using the Seraut Tutorial up to *this* point: [Finding Differentially Expressed Features (Cluster Biomarkers)](https://satijalab.org/seurat/articles/pbmc3k_tutorial.html#finding-differentially-expressed-features-cluster-biomarkers-).

My project can be found at *this* [Rpubs Link](https://rpubs.com/ap3340/SingleCellHumanPBMC) which contains the link to the dataset on the 10x genomics website.

**My question is... is the following table derived from performing the steps found in the first link I provided?**
[![enter image description here](https://i.sstatic.net/sr9p0.png)](https://i.sstatic.net/sr9p0.png)

**If so, once I derive that table ( ID, Name, L2FC, p-value, etc.), does that mean I can perform Differential Gene Expression Analysis?**

I am a complete newbie and am teaching myself all of this. Thanks, y'all.

**UPDATE:**
Yes, it did provide me that table. With this table, am I able to perform gene expression?

**Tags:** r, rna-seq, gene-expression

---

## 1 answers

### Answer by username_1 on 2022-11-29. Score: 4

The table pictured here *is* the result of a differential gene expression (DGE) test. In the table, each row is a gene (labeled with the gene ID and name) and each column (besides the first two) are the results of DGE tests. Each cluster has two sub-columns: L2FC (log2 fold-change) which tells you the average relative expression changes between cells of the cluster vs. others; and p-value, which is how significant the change is.

You can make this table yourself by performing DGE tests for all clusters. You've already done this yourself in your Rpubs script! To make it easier though, so you don't have to manually test each cluster, use the [FindAllMarkers](https://www.rdocumentation.org/packages/Seurat/versions/1.4.0/topics/FindAllMarkers) function. It will compare each cluster to all other clusters and report the L2FCs and p-values.

However, I recommend using the p\_val\_adj (adjusted p-value) instead of the p-value. This is due to a problem called multiple testing: you are performing many, many tests (each gene in each cluster is a separate statistical test), and we expect 5% of differentially expressed genes (DEGs) to be DE by random chance. Over something like 100,000 tests, that means you'd expect 5,000 of those to be false positives (which is huge). So p-value is adjusted to account for multiple tests.

Once you filter the DEG table (p-adj<0.05), you can do other downstream analyses to describe what has changed between the conditions. You already have a good start in your Rpub script with the violin plots. Cluster 6 is clearly distinct in terms of those genes' expression levels. You could try to figure out what's special with that cluster. Good luck!

> **Comments:**
> **username_0**: With my violin plots, does that mean that those two genes , c("MS4A1", "CD79A"), are found within multiple clusters?  
> **username_1**: @Antonio Each violin plot shows the expression of that gene in all cells (the black dots). As you see, most values are 0 (which is very common for scRNA-seq), along with some random low-level noise. The important bits to look at is anything above the noise level (in this case, ~0.5). Any violin with color is usually somewhat significant. For the MS4A1 gene, you can see the only clear violin is in cluster 6, while CD79A is in 6 and 20. If you look at your UMAP, cluster 6 is off to the side by itself. This may be a gene that is cell-type specific or only expressed in certain conditions.  

---

