{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,26]],"date-time":"2026-02-26T03:08:16Z","timestamp":1772075296108,"version":"3.50.1"},"reference-count":48,"publisher":"Oxford University Press (OUP)","issue":"4","license":[{"start":{"date-parts":[[2023,4,17]],"date-time":"2023-04-17T00:00:00Z","timestamp":1681689600000},"content-version":"vor","delay-in-days":16,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100000002","name":"National Institutes of Health","doi-asserted-by":"publisher","award":["T32-HG02536"],"award-info":[{"award-number":["T32-HG02536"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2023,4,3]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:sec>\n                  <jats:title>Motivation<\/jats:title>\n                  <jats:p>In a genome-wide association study, analyzing multiple correlated traits simultaneously is potentially superior to analyzing the traits one by one. Standard methods for multivariate genome-wide association study operate marker-by-marker and are computationally intensive.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Results<\/jats:title>\n                  <jats:p>We present a sparsity constrained regression algorithm for multivariate genome-wide association study based on iterative hard thresholding and implement it in a convenient Julia package MendelIHT.jl. In simulation studies with up to 100 quantitative traits, iterative hard thresholding exhibits similar true positive rates, smaller false positive rates, and faster execution times than GEMMA\u2019s linear mixed models and mv-PLINK\u2019s canonical correlation analysis. On UK Biobank data with 470\u00a0228 variants, MendelIHT completed a three-trait joint analysis (n=185\u2009656) in 20\u2009h and an 18-trait joint analysis (n=104\u2009264) in 53\u2009h with an 80\u00a0GB memory footprint. In short, MendelIHT enables geneticists to fit a single regression model that simultaneously considers the effect of all SNPs and dozens of traits.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Availability and implementation<\/jats:title>\n                  <jats:p>Software, documentation, and scripts to reproduce our results are available from https:\/\/github.com\/OpenMendel\/MendelIHT.jl.<\/jats:p>\n               <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btad193","type":"journal-article","created":{"date-parts":[[2023,4,17]],"date-time":"2023-04-17T14:24:57Z","timestamp":1681741497000},"source":"Crossref","is-referenced-by-count":0,"title":["Multivariate genome-wide association analysis by iterative hard thresholding"],"prefix":"10.1093","volume":"39","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-8342-2361","authenticated-orcid":false,"given":"Benjamin B","family":"Chu","sequence":"first","affiliation":[{"name":"Department of Computational Medicine, David Geffen School of Medicine at UCLA , Los Angeles, CA 90095-1554, United States"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Seyoon","family":"Ko","sequence":"additional","affiliation":[{"name":"Department of Computational Medicine, David Geffen School of Medicine at UCLA , Los Angeles, CA 90095-1554, United States"},{"name":"Department of Biostatistics, Fielding School of Public Health at UCLA , Los Angeles, CA 90095-1554, United States"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Jin J","family":"Zhou","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, Fielding School of Public Health at UCLA , Los Angeles, CA 90095-1554, United States"},{"name":"Department of Medicine, David Geffen School of Medicine at UCLA , Los Angeles, CA 90095-1554, United States"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Aubrey","family":"Jensen","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, Fielding School of Public Health at UCLA , Los Angeles, CA 90095-1554, United States"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Hua","family":"Zhou","sequence":"additional","affiliation":[{"name":"Department of Computational Medicine, David Geffen School of Medicine at UCLA , Los Angeles, CA 90095-1554, United States"},{"name":"Department of Biostatistics, Fielding School of Public Health at UCLA , Los Angeles, CA 90095-1554, United States"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Janet S","family":"Sinsheimer","sequence":"additional","affiliation":[{"name":"Department of Computational Medicine, David Geffen School of Medicine at UCLA , Los Angeles, CA 90095-1554, United States"},{"name":"Department of Biostatistics, Fielding School of Public Health at UCLA , Los Angeles, CA 90095-1554, United States"},{"name":"Department of Human Genetics, David Geffen School of Medicine at UCLA , Los Angeles, CA 90095-1554, United States"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-1313-5030","authenticated-orcid":false,"given":"Kenneth","family":"Lange","sequence":"additional","affiliation":[{"name":"Department of Computational Medicine, David Geffen School of Medicine at UCLA , Los Angeles, CA 90095-1554, United States"},{"name":"Department of Human Genetics, David Geffen School of Medicine at UCLA , Los Angeles, CA 90095-1554, United States"},{"name":"Department of Statistics at UCLA , Los Angeles, CA 90095-1554, United States"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2023,4,17]]},"reference":[{"key":"2023042704414742500_btad193-B1","doi-asserted-by":"crossref","first-page":"2776","DOI":"10.1093\/bioinformatics\/btx299","article-title":"FlashPCA2: principal component analysis of Biobank-scale genotype datasets","volume":"33","author":"Abraham","year":"2017","journal-title":"Bioinformatics"},{"key":"2023042704414742500_btad193-B2","doi-asserted-by":"crossref","first-page":"e1008773","DOI":"10.1371\/journal.pgen.1008773","article-title":"Scalable probabilistic PCA for large-scale genetic variation data","volume":"16","author":"Agrawal","year":"2020","journal-title":"PLoS Genet"},{"key":"2023042704414742500_btad193-B3","doi-asserted-by":"crossref","first-page":"722","DOI":"10.1002\/gepi.20623","article-title":"Stability selection for genome-wide association","volume":"35","author":"Alexander","year":"2011","journal-title":"Genet Epidemiol"},{"key":"2023042704414742500_btad193-B4","doi-asserted-by":"crossref","first-page":"2055","DOI":"10.1214\/15-AOS1337","article-title":"Controlling the false discovery rate via knockoffs","volume":"43","author":"Barber","year":"2015","journal-title":"Ann Statist"},{"key":"2023042704414742500_btad193-B5","doi-asserted-by":"crossref","first-page":"65","DOI":"10.1137\/141000671","article-title":"Julia: a fresh approach to numerical computing","volume":"59","author":"Bezanson","year":"2017","journal-title":"SIAM Rev"},{"key":"2023042704414742500_btad193-B6","doi-asserted-by":"crossref","first-page":"265","DOI":"10.1016\/j.acha.2009.04.002","article-title":"Iterative hard thresholding for compressed sensing","volume":"27","author":"Blumensath","year":"2009","journal-title":"Appl Comput Harmon Anal"},{"key":"2023042704414742500_btad193-B7","doi-asserted-by":"crossref","first-page":"232","DOI":"10.1214\/10-AOAS388","article-title":"Coordinate descent algorithms for nonconvex penalized regression, with applications to biological feature selection","volume":"5","author":"Breheny","year":"2011","journal-title":"Ann Appl Stat"},{"key":"2023042704414742500_btad193-B8","doi-asserted-by":"crossref","first-page":"giaa044","DOI":"10.1093\/gigascience\/giaa044","article-title":"Iterative hard thresholding in genome-wide association studies: generalized linear models, prior weights, and double sparsity","volume":"9","author":"Chu","year":"2020","journal-title":"GigaScience"},{"key":"2023042704414742500_btad193-B9","doi-asserted-by":"crossref","first-page":"265","DOI":"10.1093\/biomet\/68.1.265","article-title":"Some matrix-variate distribution theory: notational considerations and a Bayesian application","volume":"68","author":"Dawid","year":"1981","journal-title":"Biometrika"},{"key":"2023042704414742500_btad193-B10","doi-asserted-by":"crossref","first-page":"132","DOI":"10.1093\/bioinformatics\/btn563","article-title":"A multivariate test of association","volume":"25","author":"Ferreira","year":"2009","journal-title":"Bioinformatics"},{"key":"2023042704414742500_btad193-B11","doi-asserted-by":"crossref","first-page":"59","DOI":"10.1534\/genetics.114.171447","article-title":"Efficient multiple-trait association and estimation of genetic correlation using the matrix-variate linear mixed model","volume":"200","author":"Furlotte","year":"2015","journal-title":"Genetics"},{"key":"2023042704414742500_btad193-B12","doi-asserted-by":"crossref","first-page":"e95923","DOI":"10.1371\/journal.pone.0095923","article-title":"A comparison of multivariate genome-wide association methods","volume":"9","author":"Galesloot","year":"2014","journal-title":"PLoS One"},{"key":"2023042704414742500_btad193-B13","doi-asserted-by":"crossref","first-page":"248","DOI":"10.1002\/gepi.22276","article-title":"Ordered multinomial regression for genetic association analysis of ordinal phenotypes at Biobank scale","volume":"44","author":"German","year":"2020","journal-title":"Genet Epidemiol"},{"key":"2023042704414742500_btad193-B14","doi-asserted-by":"crossref","first-page":"1780","DOI":"10.1214\/11-AOAS455","article-title":"Bayesian variable selection regression for genome-wide association studies and other large-scale problems","volume":"5","author":"Guan","year":"2011","journal-title":"Ann Appl Stat"},{"key":"2023042704414742500_btad193-B15","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/s12859-021-04086-8","article-title":"Modern simulation utilities for genetic analysis","volume":"22","author":"Ji","year":"2021","journal-title":"BMC Bioinformatics"},{"key":"2023042704414742500_btad193-B16","doi-asserted-by":"crossref","first-page":"1749","DOI":"10.1038\/s41588-019-0530-8","article-title":"A resource-efficient tool for mixed model association analysis of large-scale data","volume":"51","author":"Jiang","year":"2019","journal-title":"Nat Genet"},{"key":"2023042704414742500_btad193-B17","doi-asserted-by":"crossref","first-page":"756","DOI":"10.1002\/gepi.22068","article-title":"Iterative hard thresholding for model selection in genome-wide association studies","volume":"41","author":"Keys","year":"2017","journal-title":"Genet Epidemiol"},{"key":"2023042704414742500_btad193-B18","author":"Ko","year":"2020"},{"key":"2023042704414742500_btad193-B19","doi-asserted-by":"crossref","first-page":"433","DOI":"10.1016\/j.ajhg.2022.01.018","article-title":"GWAS of longitudinal trajectories at biobank scale","volume":"109","author":"Ko","year":"2022","journal-title":"Am J Hum Genet"},{"key":"2023042704414742500_btad193-B20","doi-asserted-by":"crossref","first-page":"bbab256","DOI":"10.1093\/bib\/bbab256","article-title":"Computationally scalable regression modeling for ultrahigh-dimensional omics data with ParProx","volume":"22","author":"Ko","year":"2021","journal-title":"Brief Bioinform"},{"key":"2023042704414742500_btad193-B21","doi-asserted-by":"crossref","first-page":"494","DOI":"10.1214\/21-STS835","article-title":"High-performance statistical computing in the computing environments of the 2020s","volume":"37","author":"Ko","year":"2022","journal-title":"Statist Sci"},{"key":"2023042704414742500_btad193-B22","doi-asserted-by":"crossref","first-page":"308","DOI":"10.1145\/355841.355847","article-title":"Basic linear algebra subprograms for Fortran usage","volume":"5","author":"Lawson","year":"1979","journal-title":"ACM Trans Math Softw"},{"key":"2023042704414742500_btad193-B23","doi-asserted-by":"crossref","first-page":"906","DOI":"10.1038\/s41588-018-0144-6","article-title":"Mixed-model association for biobank-scale datasets","volume":"50","author":"Loh","year":"2018","journal-title":"Nat Genet"},{"key":"2023042704414742500_btad193-B24","doi-asserted-by":"crossref","first-page":"D896","DOI":"10.1093\/nar\/gkw1133","article-title":"The new NHGRI-EBI Catalog of published genome-wide association studies (GWAS Catalog)","volume":"45","author":"MacArthur","year":"2017","journal-title":"Nucleic Acids Res"},{"key":"2023042704414742500_btad193-B25","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/s12859-016-1437-3","article-title":"Scopa and META-SCOPA: software for the analysis and aggregation of genome-wide association studies of multiple correlated phenotypes","volume":"18","author":"M\u00e4gi","year":"2017","journal-title":"BMC Bioinformatics"},{"key":"2023042704414742500_btad193-B26","doi-asserted-by":"crossref","first-page":"649","DOI":"10.1093\/bioinformatics\/btz605","article-title":"gwasrapidd: an R package to query, download and wrangle GWAS catalog data","volume":"36","author":"Magno","year":"2020","journal-title":"Bioinformatics"},{"key":"2023042704414742500_btad193-B27","doi-asserted-by":"crossref","first-page":"e34861","DOI":"10.1371\/journal.pone.0034861","article-title":"MultiPhen: joint model of multiple phenotypes can increase discovery in GWAS","volume":"7","author":"O\u2019Reilly","year":"2012","journal-title":"PLoS One"},{"key":"2023042704414742500_btad193-B28","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1038\/srep38837","article-title":"Multivariate simulation framework reveals performance of multi-trait GWAS methods","volume":"7","author":"Porter","year":"2017","journal-title":"Sci Rep"},{"key":"2023042704414742500_btad193-B29","doi-asserted-by":"crossref","first-page":"459","DOI":"10.1038\/nrg2813","article-title":"New approaches to population stratification in genome-wide association studies","volume":"11","author":"Price","year":"2010","journal-title":"Nat Rev Genet"},{"key":"2023042704414742500_btad193-B30","doi-asserted-by":"crossref","first-page":"559","DOI":"10.1086\/519795","article-title":"Plink: a tool set for whole-genome association and population-based linkage analyses","volume":"81","author":"Purcell","year":"2007","journal-title":"Am J Hum Genet"},{"key":"2023042704414742500_btad193-B31","doi-asserted-by":"crossref","first-page":"e1009141","DOI":"10.1371\/journal.pgen.1009141","article-title":"A fast and scalable framework for large-scale and ultrahigh-dimensional sparse regression with application to the UK Biobank","volume":"16","author":"Qian","year":"2020","journal-title":"PLoS Genet"},{"key":"2023042704414742500_btad193-B32","doi-asserted-by":"crossref","first-page":"35","DOI":"10.1038\/ng.271","article-title":"Genome-wide association analysis of metabolic traits in a birth cohort from a founder population","volume":"41","author":"Sabatti","year":"2009","journal-title":"Nat Genet"},{"key":"2023042704414742500_btad193-B33","doi-asserted-by":"crossref","first-page":"e2105841118","DOI":"10.1073\/pnas.2105841118","article-title":"False discovery rate control in genome-wide association studies with population structure","volume":"118","author":"Sesia","year":"2021","journal-title":"Proc Natl Acad Sci USA"},{"key":"2023042704414742500_btad193-B34","doi-asserted-by":"crossref","first-page":"R111","DOI":"10.1093\/hmg\/ddv260","article-title":"Strategies for fine-mapping complex traits","volume":"24","author":"Spain","year":"2015","journal-title":"Hum Mol Genet"},{"key":"2023042704414742500_btad193-B35","doi-asserted-by":"crossref","first-page":"e65245","DOI":"10.1371\/journal.pone.0065245","article-title":"A unified framework for association analysis with multiple related phenotypes","volume":"8","author":"Stephens","year":"2013","journal-title":"PLoS One"},{"key":"2023042704414742500_btad193-B36","doi-asserted-by":"crossref","first-page":"e1001779","DOI":"10.1371\/journal.pmed.1001779","article-title":"UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age","volume":"12","author":"Sudlow","year":"2015","journal-title":"PLoS Med"},{"key":"2023042704414742500_btad193-B37","doi-asserted-by":"crossref","first-page":"e1008431","DOI":"10.1371\/journal.pgen.1008431","article-title":"Bayesian multivariate reanalysis of large genetic studies identifies many new associations","volume":"15","author":"Turchin","year":"2019","journal-title":"PLoS Genet"},{"key":"2023042704414742500_btad193-B38","doi-asserted-by":"crossref","first-page":"714","DOI":"10.1093\/bioinformatics\/btp041","article-title":"Genome-wide association analysis by lasso penalized logistic regression","volume":"25","author":"Wu","year":"2009","journal-title":"Bioinformatics"},{"key":"2023042704414742500_btad193-B39","doi-asserted-by":"crossref","first-page":"76","DOI":"10.1016\/j.ajhg.2010.11.011","article-title":"GCTA: a tool for genome-wide complex trait analysis","volume":"88","author":"Yang","year":"2011","journal-title":"Am J Hum Genet"},{"key":"2023042704414742500_btad193-B40","doi-asserted-by":"crossref","first-page":"119","DOI":"10.1016\/j.jmva.2012.01.005","article-title":"Model selection and estimation in the matrix normal graphical model","volume":"107","author":"Yin","year":"2012","journal-title":"J Multivar Anal"},{"key":"2023042704414742500_btad193-B41","doi-asserted-by":"crossref","first-page":"894","DOI":"10.1214\/09-AOS729","article-title":"Nearly unbiased variable selection under minimax concave penalty","volume":"38","author":"Zhang","year":"2010","journal-title":"Ann Statist"},{"key":"2023042704414742500_btad193-B42","first-page":"106","article-title":"Penalized regression for genome-wide association screening of sequence data","volume":"2011","author":"Zhou","year":"2011","journal-title":"Pac Symp Biocomput"},{"key":"2023042704414742500_btad193-B43","doi-asserted-by":"crossref","first-page":"311","DOI":"10.1214\/10-STS336","article-title":"Graphical processing units and high-dimensional optimization","volume":"25","author":"Zhou","year":"2010","journal-title":"Stat Sci"},{"key":"2023042704414742500_btad193-B44","doi-asserted-by":"crossref","first-page":"2375","DOI":"10.1093\/bioinformatics\/btq448","article-title":"Association screening of common and rare genetic variants by penalized regression","volume":"26","author":"Zhou","year":"2010","journal-title":"Bioinformatics"},{"key":"2023042704414742500_btad193-B45","doi-asserted-by":"crossref","first-page":"61","DOI":"10.1007\/s00439-019-02001-z","article-title":"OpenMendel: a cooperative programming project for statistical genetics","volume":"139","author":"Zhou","year":"2020","journal-title":"Hum Genet"},{"key":"2023042704414742500_btad193-B46","doi-asserted-by":"crossref","first-page":"1335","DOI":"10.1038\/s41588-018-0184-y","article-title":"Efficiently controlling for case-control imbalance and sample relatedness in large-scale genetic association studies","volume":"50","author":"Zhou","year":"2018","journal-title":"Nat Genet"},{"key":"2023042704414742500_btad193-B47","doi-asserted-by":"crossref","first-page":"407","DOI":"10.1038\/nmeth.2848","article-title":"Efficient multivariate linear mixed model algorithms for genome-wide association studies","volume":"11","author":"Zhou","year":"2014","journal-title":"Nat Methods"},{"key":"2023042704414742500_btad193-B48","doi-asserted-by":"crossref","first-page":"821","DOI":"10.1038\/ng.2310","article-title":"Genome-wide efficient mixed-model analysis for association studies","volume":"44","author":"Zhou","year":"2012","journal-title":"Nat Genet"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btad193\/49989160\/btad193.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/4\/btad193\/50105314\/btad193.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/4\/btad193\/50105314\/btad193.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,4,27]],"date-time":"2023-04-27T04:42:38Z","timestamp":1682570558000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btad193\/7126408"}},"subtitle":[],"editor":[{"given":"Tobias","family":"Marschall","sequence":"additional","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2023,4,1]]},"references-count":48,"journal-issue":{"issue":"4","published-print":{"date-parts":[[2023,4,3]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btad193","relation":{},"ISSN":["1367-4811"],"issn-type":[{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2023,4,1]]},"published":{"date-parts":[[2023,4,1]]},"article-number":"btad193"}}