text: user16549: 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).
  My project can be found at this Rpubs Link 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?
  
  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?
  
  user16565: 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 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!
  
original_text: 16549: <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.stack.imgur.com/sr9p0.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/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>
   16565: <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>
  