variancePartition: interpreting drivers of variation in complex gene expression studies
Published version of a preprint
Abstract
Background: As large-scale studies of gene expression with multiple sources of biological and technical variation become widely adopted, characterizing these drivers of variation becomes essential to understanding disease biology and regulatory genetics. Results: We describe a statistical and visualization framework, variancePartition, to prioritize drivers of variation based on a genome-wide summary, and identify genes that deviate from the genome-wide trend. Using a linear mixed model, variancePartition quantifies variation in each expression trait attributable to differences in disease status, sex, cell or tissue type, ancestry, genetic background, experimental stimulus, or technical variables. Analysis of four large-scale transcriptome profiling datasets illustrates that variancePartition recovers striking patterns of biological and technical variation that are reproducible across multiple datasets. Conclusions: Our open source software, variancePartition, enables rapid interpretation of complex gene expression studies as well as other high-throughput genomics assays. variancePartition is available from Bioconductor: http://bioconductor.org/packages/variancePartition .
The claims and statistics of this paper haven’t been checked yet.
Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.
Integrity signals
Registry
Cites
Statistics
Claims the paper makes
The claims and statistics of this paper haven’t been checked yet.
Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.
Claims other papers make about this one
Each is checked against this paper’s own text: does it say what the citing paper says it does?
Donor identity explained the most variance per gene (mean 13.5%).
“Using a linear mixed model to estimate the variance explained for a set of factors per gene30, we found that donor identity explained the most variance per gene (mean 13.5%) (Fig. 2a).”
From Genetic analysis of the human microglial transcriptome across brain regions, aging and disease pathologies (2022) · cites this paper as [30]
Self-check
Percentages, sums and reported differences re-computed from the paper’s own numbers.
The self-check needs the paper’s full text.
Cites
Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2
Love, M. I., Huber, W., & Anders, S. (2014) · Genome Biology · doi:10.1186/s13059-014-0550-8 · Cited by 95,106 · see citing works
Fitting Linear Mixed-Effects Models Using <b>lme4</b>
Bates, D., Mächler, M., Bolker, B., & Walker, S. (2015) · Journal of Statistical Software · doi:10.18637/jss.v067.i01 · Cited by 86,835 · see citing works
Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., van Baren, M. J., Salzberg, S. L., Wold, B. J., & Pachter, L. (2010) · Nature Biotechnology · doi:10.1038/nbt.1621 · Cited by 15,556 · see citing works
ggplot2: Elegant Graphics for Data Analysis
Wickham, H. (2009) · doi:10.1007/978-0-387-98141-3 · Cited by 14,447 · see citing works
Smyth, G. K. (2004) · Statistical Applications in Genetics and Molecular Biology · doi:10.2202/1544-6115.1027 · Cited by 11,829 · see citing works
Near-optimal probabilistic RNA-seq quantification
Bray, N. L., Pimentel, H., Melsted, P., & Pachter, L. (2016) · Nature Biotechnology · doi:10.1038/nbt.3519 · Cited by 9,723 · see citing works
Genome-wide efficient mixed-model analysis for association studies
Zhou, X., & Stephens, M. (2012) · Nature Genetics · doi:10.1038/ng.2310 · Cited by 3,627 · see citing works
Variance component model to account for sample structure in genome-wide association studies
Kang, H. M., Sul, J. H., Service, S. K., Zaitlen, N. A., Kong, S.-Y., Freimer, N. B., Sabatti, C., & Eskin, E. (2010) · Nature Genetics · doi:10.1038/ng.548 · Cited by 2,817 · see citing works
Maximum Likelihood Approaches to Variance Component Estimation and to Related Problems
Harville, D. A. (1977) · Journal of the American Statistical Association · doi:10.1080/01621459.1977.10480998 · Cited by 2,667 · see citing works
Repeatability for Gaussian and non-Gaussian data: a practical guide for biologists.
Nakagawa, S., & Schielzeth, H. (2010) · Biological Reviews · doi:10.1111/j.1469-185x.2010.00141.x · Cited by 2,584 · see citing works
SEQC/MAQC-III Consortium (2014) · Nature Biotechnology · doi:10.1038/nbt.2957 · Cited by 983 · cited 3× · see citing works
“While technical variation in RNA-seq is known to depend on GC content [8, 9], variancePartition gives a clear illustration of how the effect of technical artifacts varies substantially across genes.”
Advantages and pitfalls in the application of mixed-model association methods
Yang, J., Zaitlen, N. A., Goddard, M. E., Visscher, P. M., & Price, A. L. (2014) · Nature Genetics · doi:10.1038/ng.2876 · Cited by 981 · see citing works
GC-Content Normalization for RNA-Seq Data
Risso, D., Schwartz, K., Sherlock, G., & Dudoit, S. (2011) · BMC Bioinformatics · doi:10.1186/1471-2105-12-480 · Cited by 867 · see citing works
Genetic variability in the regulation of gene expression in ten regions of the human brain
UK Brain Expression Consortium, Ramasamy, A., North American Brain Expression Consortium, Trabzuni, D., Guelfi, S., Varghese, V., Smith, C., Walker, R., De, T., Coin, L., et al. (2014) · Nature Neuroscience · doi:10.1038/nn.3801 · Cited by 644 · see citing works
Polarization of the Effects of Autoimmune and Neurodegenerative Risk Alleles in Leukocytes
Raj, T., Rothamel, K., Mostafavi, S., Ye, C., Lee, M. N., Replogle, J. M., Feng, T., Lee, M., Asinovski, N., Frohlich, I., et al. (2014) · Science · doi:10.1126/science.1249547 · Cited by 506 · see citing works
Reproducibility of high-throughput mRNA and small RNA sequencing across laboratories
The GEUVADIS Consortium, 't Hoen, P. A. C., Friedländer, M. R., Almlöf, J., Sammeth, M., Pulyakhina, I., Anvar, S. Y., Laros, J. F. J., Buermans, H. P. J., Karlberg, O., et al. (2013) · Nature Biotechnology · doi:10.1038/nbt.2702 · Cited by 250 · see citing works
Detecting and correcting systematic variation in large-scale RNA sequencing data
Li, S., Łabaj, P. P., Zumbo, P., Sykacek, P., Shi, W., Shi, L., Phan, J., Wu, P.-Y., Wang, M., Wang, C., et al. (2014) · Nature Biotechnology · doi:10.1038/nbt.3000 · Cited by 175 · see citing works
Flagged first, then by citations
Research works printed but not linked: 33
- [2]GTEx Consortium The Genotype-Tissue Expression (GTEx) pilot analysis: Multitissue gene regulation in humans Science 2015 348 6235 648 60 10.1126/science.1262110 25954001 PMC4547484 · DOI · Paper page
- [4]Lee MN Ye C Villani AC Raj T Li W Eisenhaure TM Imboywa SH Chipendo PI Ran FA Slowikowski K Ward LD Raddassi K McCabe C Lee MH Frohlich IY Hafler D. a Kellis M Raychaudhuri S Zhang F Stranger BE Benoist CO De Jager PL Regev A Hacohen N Common genetic variants modulate pathogen-sensing responses in human dendritic cells Science 2014 343 6175 1246980 10.1126/science.1246980 24604203 PMC4124741 · DOI · Paper page
- [5]Fairfax BP Humburg P Makino S Naranbhai V Wong D Lau E Jostins L Plant K Andrews R McGee C Knight JC Innate immune activity conditions the effect of regulatory variants upon monocyte gene expression Science 2014 343 6175 1246949 10.1126/science.1246949 24604202 PMC4064786 · DOI · Paper page
- [6]Lappalainen T Sammeth M Friedländer MR ‘t Hoen PAC Monlong J Rivas MA Gonzàlez-Porta M Kurbatova N Griebel T Ferreira PG Barann M Wieland T Greger L van Iterson M Almlöf J Ribeca P Pulyakhina I Esser D Giger T Tikhonov A Sultan M Bertier G MacArthur DG Lek M Lizano E Buermans HPJ Padioleau I Schwarzmayr T Karlberg O Ongen H Kilpinen H Beltran S Gut M Kahlem K Amstislavskiy V Stegle O Pirinen M Montgomery SB Donnelly P McCarthy MI Flicek P Strom TM The Geuvadis Consortium Lehrach H Schreiber S Sudbrak R Carracedo Á Antonarakis SE Häsler R Syvänen AC van Ommen G-J Brazma A Meitinger T Rosenstiel P Guigó R Gut IG Estivill X Dermitzakis ET Transcriptome and genome sequencing uncovers functional variation in humans Nature 2013 501 7468 506 11 10.1038/nature12531 24037378 PMC3918453 · DOI · Paper page
- [7]Zhang B Gaiteri C Bodea LG Wang Z McElwee J Podtelezhnikov AA Zhang C Xie T Tran L Dobrin R Fluder E Clurman B Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer’s disease Cell 2013 153 3 707 20 10.1016/j.cell.2013.03.030 23622250 PMC3677161 · DOI · Paper page
- [10]Alter O Brown PO Botstein D Singular value decomposition for genome-wide expression data processing and modeling Proc Natl Acad Sci U S A 2000 97 18 10101 6 10.1073/pnas.97.18.10101 10963673 PMC27718 · DOI
- [11]Eisen MB Spellman PT Brown PO Botstein D Cluster analysis and display of genome-wide expression patterns Proc Natl Acad Sci U S A 1998 95 25 14863 8 10.1073/pnas.95.25.14863 9843981 PMC24541 · DOI · Paper page
- [13]Robinson MD McCarthy DJ Smyth GK edgeR: a Bioconductor package for differential expression analysis of digital gene expression data Bioinformatics (Oxford) 2010 26 1 139 40 10.1093/bioinformatics/btp616 PMC2796818 19910308 · DOI · Paper page
- [14]Trapnell C Hendrickson DG Sauvageau M Goff L Rinn JL Pachter L Differential analysis of gene regulation at transcript resolution with RNA-seq Nat Biotechnol 2013 31 1 46 53 10.1038/nbt.2450 23222703 PMC3869392 · DOI · Paper page
- [15]Law CW Chen Y Shi W Smyth GK Voom: precision weights unlock linear model analysis tools for RNA-seq read counts Genome Biol 2014 15 2 29 10.1186/gb-2014-15-2-r29 PMC4053721 24485249 · DOI · Paper page
- [17]Mele M Ferreira PG Reverter F DeLuca DS Monlong J Sammeth M Young TR Goldmann JM Pervouchine DD Sullivan TJ Johnson R Segre AV Djebali S Niarchou A Consortium TG Wright F. a Lappalainen T Calvo M Getz G Dermitzakis ET Ardlie KG Guigo R The human transcriptome across tissues and individuals Science 2015 348 6235 660 5 10.1126/science.aaa0355 25954002 PMC4547472 · DOI · Paper page
- [18]Rouhani F Kumasaka N de Brito MC Bradley A Vallier L Gaffney D Genetic background drives transcriptional variation in human induced pluripotent stem cells PLoS Genet 2014 10 6 1004432 10.1371/journal.pgen.1004432 PMC4046971 24901476 · DOI · Paper page
- [19]Burrows CK Banovich NE Pavlovic BJ Patterson K Gallego Romero I Pritchard JK Gilad Y Genetic variation, not cell type of origin, underlies the majority of identifiable regulatory differences in iPSCs PLoS Genet 2016 12 1 1005793 10.1371/journal.pgen.1005793 PMC4727884 26812582 · DOI
- [20]Trabzuni D Thomson PC Analysis of gene expression data using a linear mixed model/finite mixture model approach: application to regional differences in the human brain Bioinformatics 2014 30 11 1555 61 10.1093/bioinformatics/btu088 24519379 · DOI
- [21]Listgarten J Kadie C Schadt EE Heckerman D Correction for hidden confounders in the genetic analysis of gene expression Proc Natl Acad Sci 2010 107 38 16465 10.1073/pnas.1002425107 20810919 PMC2944732 · DOI
- [22]Storey JD Madeoy J Strout JL Wurfel M Ronald J Akey JM Gene-Expression Variation Within and Among Human Populations Am J Hum Genet 2007 80 3 502 9 10.1086/512017 17273971 PMC1821107 · DOI
- [24]Yang J Benyamin B McEvoy BP Gordon S Henders AK Nyholt DR Madden P. a Heath AC Martin NG Montgomery GW Goddard ME Visscher PM Common SNPs explain a large proportion of the heritability for human height Nat Genet 2010 42 7 565 9 10.1038/ng.608 20562875 PMC3232052 · DOI · Paper page
- [26]Lippert C Listgarten J Liu Y Kadie CM Davidson RI Heckerman D FaST linear mixed models for genome-wide association studies Nat Methods 2011 8 10 833 5 10.1038/nmeth.1681 21892150 · DOI · Paper page
- [27]Pirinen M Donnelly P Spencer CCA Efficient computation with a linear mixed model on large-scale data sets with applications to genetic studies Ann Appl Stat 2013 7 1 369 90 10.1214/12-AOAS586 · DOI
- [30]Laird NM Ware JH Random-effects models for longitudinal data Biometrics 1982 38 4 963 74 10.2307/2529876 7168798 · DOI · Paper page
- [31]Pinheiro J Bates D Mixed-effects models in S and S-Plus 2000 New York Springer
- [38]Gelman A Analysis of variance – why it is more important than ever Ann Stat 2005 33 1 1 53 10.1214/009053604000001048 · DOI · Paper page
- [40]Munro SA Lund SP Pine PS Binder H Clevert DA Conesa A Dopazo J Fasold M Hochreiter S Hong H Jafari N Kreil DP Labaj PP Li S Liao Y Lin SM Meehan J Mason CE Santoyo-Lopez J Setterquist RA Shi L Shi W Smyth GK Stralis-Pavese N Su Z Tong W Wang C Wang J Xu J Ye Z Yang Y Yu Y Salit M Assessing technical performance in differential gene expression experiments with external spike-in RNA control ratio mixtures Nat Commun 2014 5 5125 10.1038/ncomms6125 25254650 · DOI
- [41]Li S, Tighe SW, Nicolet CM, Grove D, Levy S, Farmerie W, Viale A, Wright C, Schweitzer PA, Gao Y, Kim D, Boland J, Hicks B, Kim R, Chhangawala S, Jafari N, Raghavachari N, Gandara J, Garcia-Reyero N, Hendrickson C, Roberson D, Rosenfeld J, Smith T, Underwood JG, Wang M, Zumbo P, Baldwin DA, Grills GS, Mason CE. Multi-platform assessment of transcriptome profiling using RNA-seq in the ABRF next-generation sequencing study. Nat Biotechnol. 2014; 32(9). doi: 10.1038/nbt.2972. PMC4167418 25150835 · DOI · Paper page
- [44]Feng H Zhang X Zhang C mRIN for direct assessment of genome-wide and gene-specific mRNA integrity from large-scale RNA-sequencing data Nat Commun 2015 6 7816 10.1038/ncomms8816 26234653 PMC4523900 · DOI
- [45]Gallego Romero I Pai AA Tung J Gilad Y RNA-seq: impact of RNA degradation on transcript quantification BMC Biol 2014 12 1 42 10.1186/1741-7007-12-42 24885439 PMC4071332 · DOI
- [46]Liao Y Smyth GK Shi W FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features Bioinformatics 2014 30 7 923 30 10.1093/bioinformatics/btt656 24227677 · DOI · Paper page
- [47]Anders S Pyl PT Huber W HTSeq-A Python framework to work with high-throughput sequencing data Bioinformatics 2015 31 2 166 9 10.1093/bioinformatics/btu638 25260700 PMC4287950 · DOI · Paper page
- [49]Patro R Mount SM Kingsford C Sailfish enables alignment-free isoform quantification from RNA-seq reads using lightweight algorithms Nat Biotechnol 2014 32 5 462 4 10.1038/nbt.2862 24752080 PMC4077321 · DOI · Paper page
- [50]Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides accurate, fast, and bias-aware transcript expression estimates using dual-phase inference. bioRxiv. 2016. doi: 10.1101/021592. · DOI
- [51]Li B Dewey CN RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome BMC Bioinformatics 2011 12 1 323 10.1186/1471-2105-12-323 21816040 PMC3163565 · DOI · Paper page
- [53]Soneson C Love MI Robinson MD Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences F1000Research 2016 4 0 1521 10.12688/f1000research.7563.2 PMC4712774 26925227 · DOI
- [54]Frazee AC Pertea G Jaffe AE Langmead B Salzberg SL Leek JT Ballgown bridges the gap between transcriptome assembly and expression analysis Nat Biotechnol 2015 33 3 243 6 10.1038/nbt.3172 25748911 PMC4792117 · DOI
Cited by
Sorted by citations
Establishing Cerebral Organoids as Models of Human-Specific Brain Evolution
Pollen, A. A., Bhaduri, A., Andrews, M. G., Nowakowski, T. J., Meyerson, O. S., Mostajo-Radji, M. A., Di Lullo, E., Alvarado, B., Bedolli, M., Dougherty, M. L., et al. (2019) · Cell · doi:10.1016/j.cell.2019.01.017 · Cited by 600 · see citing works
Crowell, H. L., Soneson, C., Germain, P.-L., Calini, D., Collin, L., Raposo, C., Malhotra, D., & Robinson, M. D. (2020) · Nature Communications · doi:10.1038/s41467-020-19894-4 · Cited by 515 · see citing works
Dream: powerful differential expression analysis for repeated measures designs
Hoffman, G. E., & Roussos, P. (2021) · Bioinformatics · doi:10.1093/bioinformatics/btaa687 · Cited by 319 · see citing works
Self-citationLopes, K. D. P., Snijders, G. J. L., Humphrey, J., Allan, A., Sneeboer, M. A. M., Navarro, E., Schilder, B. M., Vialle, R. A., Parks, M., Missall, R., et al. (2022) · Nature Genetics · doi:10.1038/s41588-021-00976-y · Cited by 273 · see citing works
Li, P., Marshall, L., Oh, G., Jakubowski, J. L., Groot, D., He, Y., Wang, T., Petronis, A., & Labrie, V. (2019) · Nature Communications · doi:10.1038/s41467-019-10101-7 · Cited by 234 · see citing works
Wingo, A. P., Dammer, E. B., Breen, M. S., Logsdon, B. A., Duong, D. M., Troncosco, J. C., Thambisetty, M., Beach, T. G., Serrano, G. E., Reiman, E. M., et al. (2019) · Nature Communications · doi:10.1038/s41467-019-09613-z · Cited by 211 · see citing works
Girdhar, K., Hoffman, G. E., Jiang, Y., Brown, L., Kundakovic, M., Hauberg, M. E., Francoeur, N. J., Wang, Y.-C., Shah, H., Kavanagh, D. H., et al. (2018) · Nature Neuroscience · doi:10.1038/s41593-018-0187-0 · Cited by 139 · see citing works
Self-citationHumphrey, J., Venkatesh, S., Hasan, R., Herb, J. T., de Paiva Lopes, K., Küçükali, F., Byrska-Bishop, M., Evani, U. S., Narzisi, G., Fagegaltier, D., et al. (2023) · Nature Neuroscience · doi:10.1038/s41593-022-01205-3 · Cited by 138 · see citing works
The proteomic landscape of synaptic diversity across brain regions and cell types
van Oostrum, M., Blok, T. M., Giandomenico, S. L., tom Dieck, S., Tushev, G., Fürst, N., Langer, J. D., & Schuman, E. M. (2023) · Cell · doi:10.1016/j.cell.2023.09.028 · Cited by 135 · see citing works
Identity domains capture individual differences from across the behavioral repertoire
Forkosh, O., Karamihalev, S., Roeh, S., Alon, U., Anpilov, S., Touma, C., Nussbaumer, M., Flachskamm, C., Kaplick, P. M., Shemesh, Y., et al. (2019) · Nature Neuroscience · doi:10.1038/s41593-019-0516-y · Cited by 110 · see citing works
Kapoor, M., Wang, J.-C., Farris, S. P., Liu, Y., McClintick, J., Gupta, I., Meyers, J. L., Bertelsen, S., Chao, M., Nurnberger, J., et al. (2019) · Translational Psychiatry · doi:10.1038/s41398-019-0384-y · Cited by 92 · see citing works
A cross-species proteomic map reveals neoteny of human synapse development
Wang, L., Pang, K., Zhou, L., Cebrián-Silla, A., González-Granero, S., Wang, S., Bi, Q., White, M. L., Ho, B., Li, J., et al. (2023) · Nature · doi:10.1038/s41586-023-06542-2 · Cited by 82 · see citing works
A cross-species proteomic map reveals neoteny of human synapse development
Wang, L., Pang, K., Zhou, L., Cebrián-Silla, A., González-Granero, S., Wang, S., Bi, Q., White, M. L., Ho, B., Li, J., et al. (2022) · doi:10.1101/2022.10.24.513541 · Cited by 82 · see citing works
The three-dimensional landscape of cortical chromatin accessibility in Alzheimer’s disease
Bendl, J., Hauberg, M. E., Girdhar, K., Im, E., Vicari, J. M., Rahman, S., Fernando, M. B., Townsley, K. G., Dong, P., Misir, R., et al. (2022) · Nature Neuroscience · doi:10.1038/s41593-022-01166-7 · Cited by 65 · see citing works
Fancy, N. N., Smith, A. M., Caramello, A., Tsartsalis, S., Davey, K., Muirhead, R. C. J., McGarry, A., Jenkyns, M. H., Schneegans, E., Chau, V., et al. (2024) · Acta Neuropathologica · doi:10.1007/s00401-024-02727-9 · Cited by 48 · see citing works
“We then used the dream function from variancePartition package in R [31] to compare the expression changes across the categorical variables or regressed against β-amyloid densities.”
MethodSpatially resolved gene signatures of white matter lesion progression in multiple sclerosis
Alsema, A. M., Wijering, M. H. C., Miedema, A., Kotah, J. M., Koster, M., Rijnsburger, M., van Weering, H. R. J., de Vries, H. E., Baron, W., Kooistra, S. M., et al. (2024) · Nature Neuroscience · doi:10.1038/s41593-024-01765-6 · Cited by 38 · see citing works
Breen, M. S., Wingo, A. P., Koen, N., Donald, K. A., Nicol, M., Zar, H. J., Ressler, K. J., Buxbaum, J. D., & Stein, D. J. (2018) · Brain, Behavior, and Immunity · doi:10.1016/j.bbi.2018.05.016 · Cited by 35 · see citing works
Establishing the relationship between brain cellular senescence and brain structure
Lund, A. N., Kopell, B. H., Golestani, N., Vornholt, E., Thompson, R. C., Liharska, L. E., Luo, Y., Cheng, E., Park, Y. J., Fennessy, B., et al. (2026) · Cell · doi:10.1016/j.cell.2025.10.014 · Cited by 6 · see citing works
Self-citation
- Authors
- Hoffman, G. E., & Schadt, E. E.
- ORCID
- Hoffman, Gabriel E.
- Venue
- BMC Bioinformatics
- Published
- 2016-11-25
- Volume
- 17
- Issue
- 1
- ISSN
- 1471-2105
- PMC
- PMC5123296
- Cited by
- 831
- References
- 17
- Language
- en
- Funders
- National Heart, Lung, and Blood Institute (U01 HL107388-04)
- National Institute on Aging (U01 AG046170-02)
About these checks
Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.
Statistics and arithmetic are recomputed, not judged.
Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.
What counts as a claim ›How we measure confidence ›How we check this ›
Are you an author, or do you see a problem with this page? Request a review