There are many ways to assess the differential expression of genes (DEG) between populations of single cells – Here we detail some of the methods available for researchers using SeqGeq.
Currently Volcano plots are the only truly statistically rigorous method of performing DEG analysis in SeqGeq. A volcano plot is named for the shape of the graph generated there, where more significant p-values tend to correspond with increased fold changes for a given comparison gene set.
After defining a set of populations you’re interested in determining DEGs for, you can use a volcano plot to ascertain genes expressed above a certain fold change and statistically significant p-value. To do so, open a gene view or pivoted plot, placing your primary population of interest (which you want the volcano plot to refer) on the y-axis there, against the population you’d like to compare against (on the x-axis), then click the DEG button at the top of the graph window there:
Note the y-axis here is an adjusted p-value, known as a “q-value”. The current p-value correction uses a popular “Bonferroni” method. The scale typically used to view this type of derived observation is inverse log transformed.
Inverse Log Transformation
A value of one will appear at the bottom of the scale indicating least significance between comparator population, while the more significant values (determined dynamically according to your data) will appear the top.
This statistic is the difference in expression for each gene between the populations under comparison. A positive fold change means the population of interest (on the y-axis) is expressing more than its comparator (on the x-axis), while negative values indicate relatively more expression in the comparator population.
Gating on the Volcano
You’ll see a dotted line indicates the commonly used 0.05 p-value by default (95% confidence level). This default value can be changed in the Graphs section of SeqGeq’s preferences. Though it is possible to create a gate within the volcano plot by hand, manually entering filter values is often highly desirable. In the Graph tab of the Volcano plot’s graph window you can choose Manual Gate:
The manual gating dialog will illustrate the minimum and maximum values for these parameters by default, so that researchers can enter their preferred significance and fold change gating values for either up or down regulated genes between the populations set in the original gene view graph.
If you have questions about filtering differentially expressed genes please contact tech support and we’ll be happy to help: email@example.com