Summary of Affymetrix Microarray Data Analysis

1. Affymetrix Microarray Data

• CEL files: contain intensity values, higher intensity (transcript abundance) more active genes
• CDF (chip description file) files: specify the probe and probe set to which each cell belongs
• Terms:
  • Probe: oligonucleotides of 25 base (pair) length used to probe RNA targets (25 base sequence)
  • Probe pair: a unit composed of a perfect match (PM) and its mismatch (MM)
  • Probe pair set: PMs and MMs related to a common affyID (a group of probe pairs corresponds to a particular gene or a fraction of a gene. Some genes are represented by more than one probe set.)
  • affyID: an identification for a probe set (which can be a gene or a fraction of a gene) represented on the array

2. Microarray Experimental Designs

• Biological and technical replicates
• Pooling (biological averaging), blocking, randomized
• Sample size determination

3. Data Exploration

• MA plots
  • M values are log fold changes, M=log2(T/C)=log2(T)-log2(C)
  • A values are average log intensities between two arrays, A=(log2(T)+log2(C))/2
• Images, residual images
• Histograms, boxplots
• RNA degradation plots

4. Data Preprocessing

• Approaches: background correction, normalization, PM correction, and summarization
  • Background correction methods:
    • rma: robust multiarray average method (Irizarry et al. 2003)
    • mas: Affymetrix Microarray Suite background correction method (2002)
    • GCRMA: modified RMA to estimate nonspecific binding (Wu et al. 2004)
  • Normalization methods:
    • quantile, contrast and loess: discussed and compared by Bolstad et al. (2003)
    • constant (scaling): taken by Affymetrix, usually done after summarization
    • invariantset: used in the dChip software (Li and Wong 2001)
    • qspline: normalized by fitting splines to the quantiles (Workman et al. 2002).
  • PM correction methods:
    • mas: an ideal mismatch subtracted from PM (Affymetrix 2002)
    • pmonly: no adjustment to the PM values.
    • subtractmm: subtract MM from PM (Affymetrix MAS 4.0 1999)
  • Summarization methods:
    • avgdiff: the average (Affymetrix MAS 4.0 1999)
    • mas: Tukey biweight on log2(PM-CM) (Affymetrix MAS 5.0 2002)
    • liwong: model-based expression index (MBEI) (Li and Wong 2001), fitting the following multi-chip model to each probeset:
      • y_ij = theta_i * phi_j + epsilon_ij, where y_ij = PM_ij - MM_ij
      • y_ij = mu_i + theta_i * phi_j + epsilon_ij, where y_ij = PM_ij
    • medianpolish: used in the RMA expression summary (Irizarry et al. 2003). A multichip linear model is fit to data from each probeset
      • y_ij = alpha_i + beta_j + epsilon_ij, where y_ij are the background-adjusted, normalized, and log-transformed PM intensities
    • playerout: Lazaridis et al. (2002)
• Popular methods

  Methods	Background correction	Normalization	PM correction	Summarization
  RMA	rma	quantile	pmonly	medianpolish (log2 scale)
  MAS5	mas	constant	mas	mas (log2 scale)
  MBEI	PM only	invariantset	pmonly or subtractmm	liwong

• Comparison of methods: compare the performance (power and FDR) of the methods using the data where the truth is known.
• R function expresso: combining the preprocessing methods together, but not every method can be combined.
  • rma background correction adjusting only PM probe intensities should only be used in conjunction with the pmonly PM correction.
  • subtractmm PM correction should not be used in conjunction with mas and medianpolish summarization methods because of likely negative corrections.

5. Analysis of Differentially Expressed Genes

• Approaches
  • Parametric test: t-test
  • Non-parametric tests: Wilcoxon sign-rank/rank-sum tests
• Linear models of microarrays (limma package):
  • linear models and design matrix
  • Contrasts and contrasts matrix
• ANOVA and MANOVA
• Multiple testings (p-value adjustments):
  • FWER: Bonferroni
  • FDR: Benjamini Hochberg

6. Clustering of Differentially Expressed Genes

• Annotation; Gene ontology
• Venn diagrams; clustering; classification
• Diagnostics

References

Bolstad (2014), bioc/vignettes/affy/inst/doc/builtinMethods.pdf

Affymetrix Microarray Suite (MAS) 5.0

Smyth (2004) Statistical Applications in Genetics and Molecular Biology

Wu et al. (2004) JASA

Bolstad et al. (2003) Bioinformatics

Irizarry et al. (2003) Biostatistics

Li and Wong (2001) PNAS

Lazaridis et al. (2002) Math Biosci.

Workman et al. (2002) Genome Biol.