Considerations when using the significance analysis of microarrays (SAM) algorithm

• Published in: BMC bioinformatics Vol. 6; no. 1; p. 129
• Main Authors: Larsson, Ola, Wahlestedt, Claes, Timmons, James A
• Format: Journal Article
• Language: English
• Published: England BioMed Central 29.05.2005; BMC
• Subjects: Adult; Aged; Algorithms; Computational Biology - methods; Correspondence; Data Interpretation, Statistical; Gene Expression Profiling - methods; Humans; Microarray Analysis; Middle Aged; Models, Genetic; Models, Statistical; Muscle, Skeletal - metabolism; Oligonucleotide Array Sequence Analysis - methods; Pattern Recognition, Automated; Phenotype; Research Design; Software; Software Design
• ISSN: 1471-2105; 1471-2105
• DOI: 10.1186/1471-2105-6-129

Users of microarray technology typically strive to use universally acceptable data analysis strategies to determine significant expression changes in their experiments. One of the most frequently utilised methods for gene expression data analysis is SAM (significance analysis of microarrays). The impact of selection thresholds, on the output from SAM, may critically alter the conclusion of a study, yet this consideration has not been systematically evaluated in any publication. We have examined the effect of discrete data selection criteria (qualification criteria for inclusion) and response thresholds (out-put filtering) on the number of significant genes reported by SAM. The use of a reduced data set by applying arbitrary restrictions vis-à-vis abundance calls (e.g. from D-chip) or application of the fold change (FC) option within SAM (named the FC hurdle hereafter), can substantially alter the significant gene list when running SAM in Microsoft Excel. We determined that for a given final FC criteria (e.g. 1.5 fold change) the FC hurdle applied within Microsoft Excel SAM alters the number of reported genes above the final FC criteria. The reason is that the FC hurdle changes the composition of the control data set, such that a different significance level (q-value) is obtained for any given gene. This effect can be so large that it changes subsequent post hoc analysis interpretation, such as ontology overrepresentation analysis. Our results argue for caution when using SAM. All data sets analysed with SAM could be reanalysed taking into account the potential impact of the use of arbitrary thresholds to trim data sets before significance testing.