PGS Documentation

Page tree

Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Analysis of variance (ANOVA) is a very powerful technique for identifying differentially expressed genes in a multi-factor experiment such as this one. In this data set, ANOVA will be used to generate a list of genes that are significantly different between Down syndrome and normal samples with an absolute difference bigger than 1.3 fold.

The ANOVA model should include Type because it is the primary factor of interest. From the exploratory analysis using the PCA plot, we observed that tissue is a large source of variation; therefore, Tissue should be included in the model. In the experiment, multiple samples were taken from the same subject, so Subject must be included in the model. If Subject were excluded from the model, the ANOVA assumption that samples within groups are independent will be violated. Additionally, the PCA scatter plot showed that the Downs syndrome and normal samples separated within tissue type, so the Type*Tissue interaction should be included in the model.

  • To invoke the ANOVA dialog, select Detect Differentially Expressed Genes in the Analysis section of the Gene Expression workflow
  • In the Experimental Factor(s) panel, select Type, Tissue and Subject by pressing <Ctrl> and left clicking each factor
  • Use the Add Factor > button to move the selections to the ANOVA Factor(s) panelTo specify
  • the interaction, select Select both Type and Tissue by pressing holding <Ctrl> on your keyboard and left clicking each factor.
  • Select the Add Interaction > button to add the Type * Tissue interaction to the ANOVA Factor(s) panel (Figure 1).  

Do NOT select OK or Apply. We will be adding contrasts to this ANOVA model in an upcoming section of the tutorial.

...

Numbered figure captions
SubtitleTextConfiguring ANOVA factors and interactions
AnchorNameANOVA Configuration

...

Image Added

Random vs.

...

fixed effects – mixed model ANOVA

Most factors in ANOVA are fixed effects, whose levels in a data set represent all the levels of interest. In this study, Type and Tissue are fixed effects. If the levels of a factor in a data set only represent a random sample of all the levels of interest (for example, Subject), the factor is a random effect. The ten subjects in this study represent only a random sample of the global population about which inferences are being made. Random effects are colored red on the spreadsheet and in the ANOVA dialog. When the ANOVA model includes both random and fixed factors, it is a mixed-model ANOVA.

...

You can specify which factors are random and which are fixed when you import your data or after importing by right-clicking on the column corresponding to a categorical variable, selecting Properties, and checking Random effectEffect. By doing that, the ANOVA will automatically know which factors to treat as random and which factors to treat as fixed.

Nested/Nesting Relationships

The subject factor in the ANOVA model is listed as “5. Subject (3. Type)”, which means that Subject is nested in Type. Partek Genomics Suite can automatically detect this sort of hierarchical design and will adjust the ANOVA calculation accordingly.

Linear Contrasts 

By default, an ANOVA only outputs a p-value for each factor/interaction. To get the fold change and ratio between Down syndrome and normal samples, a contrast must be set up.

  • Select Contrasts… to invoke the Configure dialog
  • Choose 6Choose 3.Type from the Select Factor/Interaction drop-down list. The levels in this factor are listed on the Candidate Level(s) panel on the left side of the dialog
  • Left click to select Down Syndrome from the Candidate Level(s) panel and move it to the Group 1 panel (renamed Down Syndrome) by selecting Add Contrast Level > in the top half of the dialog. 

Label 1 will be changed to the subgroup name automatically, but you can also manually specify the label

...

name.

  • Select Normal from the Candidate Level(s) panel and move it to the Group 2 panel (renamed Normal) 
  • The Add Contrast button can now be selected (Figure 2)

Numbered figure captions
SubtitleTextAdding a contrast of Down Syndrome and Normal samples
AnchorNameConfiguring ANOVA Contrasts

...

Image Added

Because the data is log2 transformed, Partek Genomics Suite will automatically detect this and will automatically select Yes for Data is already log transformed? in the top right-hand corner of the dialog. Partek Genomics Suite will use the geometric mean of the samples in each group to calculate the fold change and mean ratio for the contrast between the Down syndrome and normal samples.

  • Select Add Contrast to add the Down Syndrome vs. Normal contrast 
  • Select OK to apply the configuration
  • If successfully added, the Contrasts… button will now read Contrasts Included (Figure 3)

Numbered figure captions
SubtitleTextANOVA configuration with contrasts included
AnchorNameANOVA Configuration with Contrasts

...

Image Added

  • By default, Specify Output File is checked and gives a name to the output file. If you are trying to determine which factors should be included in the model and you do not wish to save the output file, simply uncheck this box
  • Select OK in the ANOVA dialog to compute the 3-way mixed-model ANOVA

Several progress messages will display in the lower left-hand side of the ANOVA dialog while the results are being calculated.

The result will be displayed in a child spreadsheet, ANOVA-3way (ANOVAResults). In this spreadsheet, each row represents a gene probe set and the columns represent the computation results for that gene probe set (Figure 4).  By Although not synonymous, probe set and gene will be treated as synonyms in this tutorial for convenience. By default, the genes are sorted in ascending order by the p-value of the first categorical factor. In this tutorial,Type is the first categorical factor, which means the most highly significant differently expressed gene between Down syndrome and normal samples is at the top of the spreadsheet in row 1.

 

Numbered figure captions
SubtitleTextANOVA spreadsheet
AnchorNameANOVA Spreadsheet

...

Image Added

For additional information about ANOVA in Partek Genomics Suite, see Chapter 11 Inferential Statistics in the User’s Manual (Help > User’s Manual).

Visualizing ANOVA

...

results

Deciding which factors to include in the ANOVA may be an iterative process while you decide which factors and interactions are relevant as not all factors have to be included in the model. For example, in this tutorial, Gender and Scan date were not included.  The Sources of Variation plot is a way to quantify the relative contribution of each factor in the model towards explaining the variability of the data.

  • View the sources of variation for each of the factors across the whole genome by clicking Plot Select View Sources of Variation from the Analysis section of the Gene Expression workflow with the ANOVA result spreadsheet activeactive 

Sources of Variation tab will appear (Figure 5) with a bar chart showing the signal to noise ratio for each factor accross the whole genome. Sources of variation can also be viewed as a pie chart showing sum or squares by selecting the Pie Chart (Sum of Squares) tab in the upper left-hand side of the Sources of Variation tab.

 

Numbered figure captions
SubtitleTextSources of Variation tab showing a bar chart
AnchorNameANOVA Sources of Variation

...

Image Added

This plot presents the mean signal-to-noise ratio of all the genes on the microarray. All the non-random factors in the ANOVA model are listed on the X-axis (including error). The Y-axis represents the mean of the ratios of mean square of all the genes to the mean square error of all the genes. Mean square is ANOVA’s measure of variance. Compare the bar for each signal to the bar for error; if a factor's bar is higher than error's bar, that factor contributed significant variation to the data across all the variables. Notice that this plot is very consistent with the results in the PCA scatter plot. In this data, on average, Tissue is the largest source of variation.

...

  • Right-click on a row header in the ANOVA spreadsheet (Figure 6)
  • Select ANOVA Interaction Plot to generate an Interaction Plot tab for that individual gene

Numbered figure captions
SubtitleTextCalling an ANOVA Interaction Plot for a gene
AnchorNameCalling ANOVA Interaction Plot

...

Image Added

Generate these plots for rows 3 (DSCR3) and 8 (CSTB). If the lines in the interaction plot are not parallel, then there is a chance that there is an interaction between Tissue and Type. Error bars show standard error of the least squared mean. DSCR3 is a good example of this (Figure 7). We can look at the p-values in column 9, p-value(Type * Tissue) to check if this apparent interaction is statistically significant. 

 

Numbered figure captions
SubtitleTextInteraction Plot for DSCR3
AnchorNameInteraction Plot

...

Image Added

We can view the expression levels of a gene for each sample using a dot plot. 

  • Right click on the gene row header and select Dot Plot (Orig. Data) from the pop-up menu. This generates a Dot Plot tab for the selected gene (Figure 8)

Numbered figure captions
SubtitleTextDot Plot showing DSCR3 expression levels for each sample
AnchorNameDot Plot

...

Image Added

In the plot, each dot is a sample of the original data. The Y-axis represents the log2 normalized intensity of the gene and the X-axis represents the different types of samples. The median expression of each group is different from each other in this example. The median of the Down syndrome samples is ~6.3, but the median of the normal samples is ~6.0. The line inside the Box & Whiskers represents the median of the samples in a group. Placing the mouse cursor over a Box & Whiskers plot will show its median and range. 

Create Gene List

...

 

...

Numbered figure captions
SubtitleTextCreating a gene list from ANOVA results
AnchorNameList Manager

Image Removed

The spreadsheet Down_Syndrome_vs_Normal (A) will be created as a child spreadsheet under the Down_Syndrome-GE spreadsheet.

This gene list spreadsheet can now be used for further analysis such as hierarchical clustering, gene ontology, integration of copy number data, or be exported into other data analysis tools such as pathway analysis.

You can practice creating new gene list criteria of your own to become familiar with the List Manager tool. For more information, you can always click on the (Image Removed) buttons.

Generating Gene Lists from a Volcano Plot

Next, we will generate a list of genes that passed a p-value threshold of 0.05 and fold-changes greater than 1.3 using a volcano plot.

...

Numbered figure captions
SubtitleTextGenerating a Volcano Plot from ANOVA results
AnchorNameGenerating Volcano Plot

Image Removed

...

Numbered figure captions
SubtitleTextVolcano plot generated from ANOVA spreadsheet
AnchorNameVolcano Plot

Image Removed

In the plot, each dot represents a gene. The X-axis represents the fold change of the contrast (Down syndrome vs. Normal), and the Y-axis represents the range of p-values. The genes with increased expression in Down syndrome samples are on the right side of the N/C (no change) line; genes with reduced expression in Down syndrome samples are on the left. The genes become more statistically significant with increasing Y-axis position. The genes that have larger and more significant changes between the Down syndrome and normal groups are on the upper right and upper left corner. 

In order to select the genes by fold-change and p-value, we will draw a horizontal line to represent the p-value 0.05 and two vertical lines indicating the –1.3 and 1.3-fold changes (cutoff lines).

...

Numbered figure captions
SubtitleTextSetting cutoff lines for -1.3 to 1.3 fold changes and a p-value of 0.05
AnchorNameSet Cutoff Lines

Image Removed

  • Select OK to draw the cutoff lines
  • Select OK in the Plot Rendering Properties dialog to close the dialog and view the plot

The plot will be divided into six sections. By clicking on the upper-right section, all genes in that section will be selected.

...

Numbered figure captions
SubtitleTextCreating a gene list from a volcano plot
AnchorNameCreate List from Volcano Plot

Image Removed

Note: If no column is selected in the parent (ANOVA) spreadsheet, all of the columns will be included in the gene list; if some columns are selected, only the selected columns will be included in the list.

  • Specify a name for the gene list (example: volcano plot list) and write a brief description about the list. The description is shown when you right-click on the spreadsheet > Info > Comments. Here, I have named the list "volcano plot list" and described it as "Genes with >1.3 fold change and p-value <0.05" (Figure 14). 

 The list can be saved as a text file (File > Save As Text File) for use in reports or by downstream analysis software.

 

Numbered figure captions
SubtitleTextSaving a list created from a volcano plot
AnchorNameSave List Dialog

Image Removed

 

Page Turner
button-linkstrue

...