Abstract
It has unambiguously been shown that genetic, environmental, demographic, and technical factors may have substantial effects on gene expression levels. In addition to the measured variable(s) of interest, there will tend to be sources of signal due to factors that are unknown, unmeasured, or too complicated to capture through simple models. We show that failing to incorporate these sources of heterogeneity into an analysis can have widespread and detrimental effects on the study. Not only can this reduce power or induce unwanted dependence across genes, but it can also introduce sources of spurious signal to many genes. This phenomenon is true even for well-designed, randomized studies. We introduce "surrogate variable analysis" (SVA) to overcome the problems caused by heterogeneity in expression studies. SVA can be applied in conjunction with standard analysis techniques to accurately capture the relationship between expression and any modeled variables of interest. We apply SVA to disease class, time course, and genetics of gene expression studies. We show that SVA increases the biological accuracy and reproducibility of analyses in genome-wide expression studies.
MeSH Terms
Algorithms
Breast Neoplasms/genetics
Computer Simulation
Data Interpretation, Statistical
Epigenesis, Genetic
Female
Gene Expression
Genes, BRCA1
Genes, BRCA2
Genetic Heterogeneity
Genetic Linkage
Genome, Fungal
Genome, Human
Humans
Kidney/metabolism
Linear Models
Mutation
Oligonucleotide Array Sequence Analysis
Quantitative Trait, Heritable
Reproducibility of Results
Saccharomyces cerevisiae/genetics,metabolism
Time Factors
Authors & Affiliations
2 authors, click to expand affiliations / ORCID
Leek Jeffrey T
Department of Biostatistics, University of Washington, Seattle, Washington, USA.
Storey John D
Conflict of Interest
Competing interests. The authors have declared that no competing interests exist.
References (33)
33 references, click to expand
-
A new approach to intensity-dependent normalization of two-channel microarrays.
Biostatistics. 2007 Jan;8(1):128-39
PMID: 16636140
-
Genetic interactions between polymorphisms that affect gene expression in yeast.
Nature. 2005 Aug 4;436(7051):701-3
PMID: 16079846
-
Molecular classification of familial non-BRCA1/BRCA2 breast cancer.
Proc Natl Acad Sci U S A. 2003 Mar 4;100(5):2532-7
PMID: 12610208
-
Genetic dissection of transcriptional regulation in budding yeast.
Science. 2002 Apr 26;296(5568):752-5
PMID: 11923494
-
Genetic analysis of genome-wide variation in human gene expression.
Nature. 2004 Aug 12;430(7001):743-7
PMID: 15269782
-
Significance analysis of time course microarray experiments.
Proc Natl Acad Sci U S A. 2005 Sep 6;102(36):12837-42
PMID: 16141318
-
Fundamental patterns underlying gene expression profiles: simplicity from complexity.
Proc Natl Acad Sci U S A. 2000 Jul 18;97(15):8409-14
PMID: 10890920
-
Correlation between gene expression levels and limitations of the empirical bayes methodology for finding differentially expressed genes.
Stat Appl Genet Mol Biol. 2005;4:Article34
PMID: 16646853
-
Principal components analysis corrects for stratification in genome-wide association studies.
Nat Genet. 2006 Aug;38(8):904-9
PMID: 16862161
-
Statistical significance for genomewide studies.
Proc Natl Acad Sci U S A. 2003 Aug 5;100(16):9440-5
PMID: 12883005
-
Genetics of gene expression surveyed in maize, mouse and man.
Nature. 2003 Mar 20;422(6929):297-302
PMID: 12646919
-
A transcriptional profile of aging in the human kidney.
PLoS Biol. 2004 Dec;2(12):e427
PMID: 15562319
-
Exploring the metabolic and genetic control of gene expression on a genomic scale.
Science. 1997 Oct 24;278(5338):680-6
PMID: 9381177
-
Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factors.
Nat Genet. 2003 Sep;35(1):57-64
PMID: 12897782
-
Estimation of false discovery proportion under general dependence.
Bioinformatics. 2006 Dec 15;22(24):3025-31
PMID: 17046978
-
Remarks on Parallel Analysis.
Multivariate Behav Res. 1992 Oct 1;27(4):509-40
PMID: 26811132
-
Multiple locus linkage analysis of genomewide expression in yeast.
PLoS Biol. 2005 Aug;3(8):e267
PMID: 16035920
-
Issues in cDNA microarray analysis: quality filtering, channel normalization, models of variations and assessment of gene effects.
Nucleic Acids Res. 2001 Jun 15;29(12):2549-57
PMID: 11410663
-
Integrative analysis of the cancer transcriptome.
Nat Genet. 2005 Jun;37 Suppl:S31-7
PMID: 15920528
-
Molecular heterogeneity of inflammatory breast cancer: a hyperproliferative phenotype.
Clin Cancer Res. 2006 Sep 1;12(17):5047-54
PMID: 16951220
-
Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation.
Nucleic Acids Res. 2002 Feb 15;30(4):e15
PMID: 11842121
-
Assessing stability of gene selection in microarray data analysis.
BMC Bioinformatics. 2006 Feb 01;7:50
PMID: 16451725
-
Some comments on instability of false discovery rate estimation.
J Bioinform Comput Biol. 2006 Oct;4(5):1057-68
PMID: 17099941
-
A reanalysis of a published Affymetrix GeneChip control dataset.
Genome Biol. 2006;7(3):401
PMID: 16563185
-
Singular value decomposition for genome-wide expression data processing and modeling.
Proc Natl Acad Sci U S A. 2000 Aug 29;97(18):10101-6
PMID: 10963673
-
Fluorescent cDNA microarray hybridization reveals complexity and heterogeneity of cellular genotoxic stress responses.
Oncogene. 1999 Jun 17;18(24):3666-72
PMID: 10380890
-
Experimental design for gene expression microarrays.
Biostatistics. 2001 Jun;2(2):183-201
PMID: 12933549
-
Cluster analysis and display of genome-wide expression patterns.
Proc Natl Acad Sci U S A. 1998 Dec 8;95(25):14863-8
PMID: 9843981
-
Treating expression levels of different genes as a sample in microarray data analysis: is it worth a risk?
Stat Appl Genet Mol Biol. 2006;5:Article9
PMID: 16646873
-
Genomic expression programs in the response of yeast cells to environmental changes.
Mol Biol Cell. 2000 Dec;11(12):4241-57
PMID: 11102521
-
Analysis of variance for gene expression microarray data.
J Comput Biol. 2000;7(6):819-37
PMID: 11382364
-
The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease.
Science. 2006 Sep 29;313(5795):1929-35
PMID: 17008526
-
Gene-expression profiles in hereditary breast cancer.
N Engl J Med. 2001 Feb 22;344(8):539-48
PMID: 11207349