Efficient multivariate linear mixed model algorithms for genome-wide association studies
Open accessAbstract
Multivariate linear mixed models (mvLMMs) are powerful tools for testing associations between single-nucleotide polymorphisms and multiple correlated phenotypes while controlling for population stratification in genome-wide association studies. We present efficient algorithms in the genome-wide efficient mixed model association (GEMMA) software for fitting mvLMMs and computing likelihood ratio tests. These algorithms offer improved computation speed, power and P-value calibration over existing methods, and can deal with more than two phenotypes.
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?
In TwinsUK (females only), GWASs used a linear mixed model with an empirical genetic relatedness matrix to account for cryptic pedigree and population structure.
“GWASs in TwinsUK (females only) were performed using GEMMA [55], which implements an LME model with an empirical genetic relatedness matrix to account for cryptic pedigree and population structure.”
From Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci (2023) · cites this paper as [55]
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
PLINK: a tool set for whole-genome association and population-based linkage analyses.
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R., Bender, D., Maller, J., Sklar, P., de Bakker, P. I. W., Daly, M. J., et al. (2007) · The American Journal of Human Genetics · doi:10.1086/519795 · Cited by 34,089 · see citing works
GCTA: a tool for genome-wide complex trait analysis
Yang, J., Lee, S. H., Goddard, M. E., & Visscher, P. M. (2011) · The American Journal of Human Genetics · doi:10.1016/j.ajhg.2010.11.011 · Cited by 8,279 · see citing works
Yu, J., Pressoir, G., Briggs, W. H., Bi, I. V., Yamasaki, M., Doebley, J., McMullen, M. D., Gaut, B. S., Nielsen, D. M., Holland, J. B., et al. (2005) · Nature Genetics · doi:10.1038/ng1702 · Cited by 4,155 · 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 · cited 18× · see citing works
“Consequently, existing methods cannot be used to perform the LRT for mvLMMs in GWAS settings.”
BackgroundMethodInfluentialSelf-citationVariance 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
Mixed linear model approach adapted for genome-wide association studies
Zhang, Z., Ersoz, E. S., Lai, C., Todhunter, R. J., Tiwari, H. K., Gore, M. A., Bradbury, P. J., Yu, J., Arnett, D. K., Ordovás, J. M., et al. (2010) · Nature Genetics · doi:10.1038/ng.546 · Cited by 2,132 · see citing works
Efficient Control of Population Structure in Model Organism Association Mapping
Kang, H. M., Zaitlen, N., Wade, C. M., Kirby, A. W., Heckerman, D. E., Daly, M. J., & Eskin, E. (2008) · Genetics · doi:10.1534/genetics.107.080101 · Cited by 1,848 · see citing works
Applications of linear models in animal breeding
Henderson, C. R. (1984) · University of Guelph eBooks · Cited by 1,700 · see citing works
FaST linear mixed models for genome-wide association studies
Lippert, C., Listgarten, J., Liu, Y., Kadie, C. M., Davidson, R. I., & Heckerman, D. (2011) · Nature Methods · doi:10.1038/nmeth.1681 · Cited by 1,209 · see citing works
Estimating genetic parameters in natural populations using the ‘animal model’
Kruuk, L. E. B. (2004) · Philosophical Transactions of the Royal Society B Biological Sciences · doi:10.1098/rstb.2003.1437 · Cited by 1,196 · see citing works
Meyer, K. (2007) · Journal of Zhejiang University SCIENCE B · doi:10.1631/jzus.2007.b0815 · Cited by 924 · see citing works
Gilmour, A. R., Thompson, R., & Cullis, B. R. (1995) · Biometrics · doi:10.2307/2533274 · Cited by 915 · see citing works
Polygenic Modeling with Bayesian Sparse Linear Mixed Models
Zhou, X., Carbonetto, P. S., & Stephens, M. (2013) · PLoS Genetics · doi:10.1371/journal.pgen.1003264 · Cited by 801 · see citing works
Self-citationGenome-wide association analysis of metabolic traits in a birth cohort from a founder population
Sabatti, C., Service, S. K., Hartikainen, A., Pouta, A. M., Ripatti, S., Brodsky, J. F., Jones, C. G., Zaitlen, N., Varilo, T., Kaakinen, M., et al. (2008) · Nature Genetics · doi:10.1038/ng.271 · Cited by 721 · see citing works
Increased accuracy of artificial selection by using the realized relationship matrix
Hayes, B. J., Visscher, P. M., & Goddard, M. E. (2009) · Genetics Research · doi:10.1017/s0016672308009981 · Cited by 665 · see citing works
Lee, S. H., Yang, J., Goddard, M. E., Visscher, P. M., & Wray, N. R. (2012) · Bioinformatics · doi:10.1093/bioinformatics/bts474 · Cited by 635 · see citing works
Robust variance-components approach for assessing genetic linkage in pedigrees.
Amos, C. I. (1994) · PubMed · Cited by 520 · see citing works
Population Structure and Cryptic Relatedness in Genetic Association Studies
Astle, W. J., & Balding, D. J. (2009) · Statistical Science · doi:10.1214/09-sts307 · Cited by 455 · see citing works
Korte, A., Vilhjálmsson, B. J., Segura, V., Platt, A., Long, Q., & Nordborg, M. (2012) · Nature Genetics · doi:10.1038/ng.2376 · Cited by 440 · see citing works
Improved linear mixed models for genome-wide association studies
Listgarten, J., Lippert, C., Kadie, C. M., Davidson, R. I., Eskin, E., & Heckerman, D. E. (2012) · Nature Methods · doi:10.1038/nmeth.2037 · Cited by 403 · see citing works
MultiPhen: Joint Model of Multiple Phenotypes Can Increase Discovery in GWAS
O’Reilly, P. F., Hoggart, C. J., Pomyen, Y., Calboli, F. C. F., Elliott, P., Järvelin, M., & Coin, L. (2012) · PLoS ONE · doi:10.1371/journal.pone.0034861 · Cited by 358 · see citing works
Meyer, & Karin (2007) · Acta Scientiarum Naturalium Universitatis Sunyatseni · Cited by 327 · see citing works
A high-resolution association mapping panel for the dissection of complex traits in micePreprint
Bennett, B. J., Farber, C. R., Orozco, L. D., Kang, H. M., Ghazalpour, A., Siemers, N. O., Neubauer, M., Neuhaus, I., Yordanova, R., Guan, B., et al. (2010) · Genome Research · doi:10.1101/gr.099234.109 · Cited by 314 · see citing works
A multivariate test of association
Ferreira, M. A. R., & Purcell, S. (2008) · Bioinformatics · doi:10.1093/bioinformatics/btn563 · Cited by 247 · see citing works
A Unified Framework for Association Analysis with Multiple Related Phenotypes
Stephens, M. (2013) · PLoS ONE · doi:10.1371/journal.pone.0065245 · Cited by 241 · see citing works
Self-citationEstimating variances and covariances for multivariate animal models by restricted maximum likelihood
Meyer, K. (1991) · Genetics Selection Evolution · doi:10.1186/1297-9686-23-1-67 · Cited by 239 · see citing works
Statistical Estimation of Correlated Genome Associations to a Quantitative Trait Network
Kim, S. Y., & Xing, E. P. (2009) · PLoS Genetics · doi:10.1371/journal.pgen.1000587 · Cited by 218 · see citing works
Heritability and Genetic Correlations Explained by Common SNPs for Metabolic Syndrome Traits
Vattikuti, S., Guo, J., & Chow, C. C. (2012) · PLoS Genetics · doi:10.1371/journal.pgen.1002637 · Cited by 217 · see citing works
Price, A. L., Helgason, A., Þorleifsson, G., McCarroll, S. A., Kong, A., & Stefánsson, K. (2011) · PLoS Genetics · doi:10.1371/journal.pgen.1001317 · Cited by 200 · see citing works
Bayesian Quantitative Trait Loci Mapping for Multiple Traits
Banerjee, S., Yandell, B. S., & Yi, N. (2008) · Genetics · doi:10.1534/genetics.108.088427 · Cited by 121 · see citing works
Trzaskowski, M., Yang, J., Visscher, P., & Plomin, R. (2013) · Molecular Psychiatry · doi:10.1038/mp.2012.191 · Cited by 113 · see citing works
Pirinen, M., Donnelly, P., & Spencer, C. C. A. (2013) · The Annals of Applied Statistics · doi:10.1214/12-aoas586 · Cited by 99 · see citing works
Runcie, D. E., & Mukherjee, S. (2013) · Genetics · doi:10.1534/genetics.113.151217 · Cited by 72 · see citing works
Meyer, K., Johnston, D. J., & Graser, H. U. (2004) · Australian Journal of Agricultural Research · doi:10.1071/ar03164 · Cited by 55 · see citing works
Kostem, E., & Eskin, E. (2013) · The American Journal of Human Genetics · doi:10.1016/j.ajhg.2013.03.010 · Cited by 26 · see citing works
PX * AI: algorithmics for better convergence in restricted maximum likelihood estimation.
Meyer, K. (2006) · Cited by 10 · see citing works
Increased accuracy of artificial selection by using the realized relationship matrix: Erratum
(2009) · Genetics Research · doi:10.1017/s0016672309000111 · Cited by 4 · see citing works
Flagged first, then by citations
Cited by
Sorted by citations
Methods of integrating data to uncover genotype-phenotype interactions
Ritchie, M. D., Holzinger, E. R., Li, R., Pendergrass, S. A., & Kim, D. (2015) · Nature Reviews Genetics · doi:10.1038/nrg3868 · Cited by 1,085 · see citing works
Genetics and intelligence differences: five special findings
Plomin, R., & Deary, I. J. (2015) · Molecular Psychiatry · doi:10.1038/mp.2014.105 · Cited by 702 · see citing works
Adjusting for Heritable Covariates Can Bias Effect Estimates in Genome-Wide Association Studies
Aschard, H., Vilhjálmsson, B. J., Joshi, A. D., Price, A. L., & Kraft, P. (2015) · The American Journal of Human Genetics · doi:10.1016/j.ajhg.2014.12.021 · Cited by 287 · see citing works
Zhang, Y., Qi, G., Park, J.-H., & Chatterjee, N. (2018) · Nature Genetics · doi:10.1038/s41588-018-0193-x · Cited by 288 · see citing works
Genetic architecture of plant stress resistance: Multi-trait genome-wide association mapping
Thoen, M. P. M., Davila Olivas, N. H., Kloth, K. J., Coolen, S., Huang, P., Aarts, M. G. M., Bac‐Molenaar, J. A., Bakker, J., Bouwmeester, H. J., Broekgaarden, C., et al. (2017) · New Phytologist · doi:10.1111/nph.14220 · Cited by 188 · see citing works
Insights into the genetic architecture of the human face
White, J. D., Indencleef, K., Naqvi, S., Eller, R. J., Hoskens, H., Roosenboom, J., Lee, M. K., Li, J., Mohammed, J., Richmond, S., et al. (2021) · Nature Genetics · doi:10.1038/s41588-020-00741-7 · Cited by 166 · see citing works
Statistical methods to detect pleiotropy in human complex traits
Hackinger, S., & Zeggini, E. (2017) · Open Biology · doi:10.1098/rsob.170125 · Cited by 143 · see citing works
Hata, T., Miyata, N., Takakura, S., Yoshihara, K., Asano, Y., Kimura-Todani, T., Yamashita, M., Zhang, X.-T., Watanabe, N., Mikami, K., et al. (2019) · Endocrinology · doi:10.1210/en.2019-00408 · Cited by 116 · see citing works
Novel genetic loci affecting facial shape variation in humans
Xiong, Z., Dankova, G., Howe, L. J., Lee, M. K., Hysi, P. G., de Jong, M. A., Zhu, G., Adhikari, K., Li, D., Li, Y., et al. (2019) · eLife · doi:10.7554/elife.49898 · Cited by 92 · see citing works
Community-wide genome sequencing reveals 30 years of Darwin’s finch evolution
Enbody, E. D., Sendell-Price, A. T., Sprehn, C. G., Rubin, C.-J., Visscher, P. M., Grant, B. R., Grant, P. R., & Andersson, L. (2023) · Science · doi:10.1126/science.adf6218 · Cited by 69 · see citing works
Imputing Phenotypes for Genome-wide Association Studies
Hormozdiari, F., Kang, E. Y., Bilow, M., Ben-David, E., Vulpe, C., McLachlan, S., Lusis, A. J., Han, B., & Eskin, E. (2016) · The American Journal of Human Genetics · doi:10.1016/j.ajhg.2016.04.013 · Cited by 53 · see citing works
Liu, Z., & Lin, X. (2019) · Journal of the American Statistical Association · doi:10.1080/01621459.2018.1513363 · Cited by 36 · see citing works
Heritability of interpack aggression in a wild pedigreed population of North American grey wolves
vonHoldt, B. M., DeCandia, A. L., Heppenheimer, E., Janowitz‐Koch, I., Shi, R., Zhou, H., German, C. A., Brzeski, K. E., Cassidy, K. A., Stahler, D. R., et al. (2020) · Molecular Ecology · doi:10.1111/mec.15349 · Cited by 22 · see citing works
“We employed a linear mixed model (LMM) in gemma (Zhou & Stephens, 2014) and included a kinship relatedness matrix estimated for 391 individuals in related (see above section for more details).”
MethodSun, J., Wang, W., Zhang, R., Duan, H., Tian, X., Xu, C., Li, X., & Zhang, D. (2022) · Translational Psychiatry · doi:10.1038/s41398-022-02074-x · Cited by 18 · see citing works
Li, Y., Xiong, Z., Zhang, M., Hysi, P. G., Qian, Y., Adhikari, K., Weng, J., Wu, S., Du, S., Gonzalez-Jose, R., et al. (2023) · PLOS Genetics · doi:10.1371/journal.pgen.1010786 · Cited by 9 · see citing works
“PLOS Genetics | https://doi.org/10.1371/journal.pgen.1010786 July 17, 2023 16 / 25 Genetic correlation, heritability, GWASs, meta-analyses Genetic correlations were estimated from all SNPs using GEMMA [55] based on a multivariate linear mixed model.”
MethodInfluentialLiu, Z., Barnett, I., & Lin, X. (2020) · The Annals of Applied Statistics · doi:10.1214/19-aoas1312 · Cited by 7 · see citing works
- Authors
- Zhou, X., & Stephens, M.
- Affiliations
- University of Chicago (US)
- Corresponding author
- Zhou, Xiang
- Stephens, Matthew
- Venue
- Nature Methods
- Published
- 2014-02-16
- Volume
- 11
- Issue
- 4
- Pages
- 407-409
- ISSN
- 1548-7091 (ISSN-L)
- 1548-7105
- PMC
- PMC4211878
- OpenAlex
- W1952660301
- Cited by
- 841
- References
- 37
- Citations by year
- 2026: 63 · 2025: 63 · 2024: 77 · 2023: 76 · 2022: 83 · 2021: 70 · 2020: 74 · 2019: 80 · 2018: 85 · 2017: 58 · 2016: 57 · 2015: 40 · 2014: 14
- FWCI
- 37.12
- Type
- journal article
- Language
- en
- Open access
- green
- www.ncbi.nlm.nih.gov
- Licence
- www.springer.com/tdm (tdm)
- Locations
- Nature Methods (publishedVersion)
- PubMed (publishedVersion)
- Europe PMC (PubMed Central) (acceptedVersion)
- PubMed Central (acceptedVersion)
- Funders
- Broad Institute
- Terveyden ja hyvinvoinnin laitos
- National Institutes of Health
- National Heart, Lung, and Blood Institute (R01HL092206, HL092206)
- National Human Genome Research Institute (R56HG002585, HG02585, R01HG002585)
- Topics
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genetic Associations and Epidemiology
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genetic and phenotypic traits in livestock
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genetic Mapping and Diversity in Plants and Animals
- Keywords
- multivariate linear mixed model, genome-wide association study, Gemma, likelihood ratio test, single nucleotide polymorphisms, correlated phenotypes, population stratification, computational efficiency, p-value calibration
- MeSH terms
- Algorithms, Software, Multivariate Analysis, Likelihood Functions, Linear Models, Genome, Polymorphism, Single Nucleotide
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