{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,4,13]],"date-time":"2026-04-13T15:53:43Z","timestamp":1776095623276,"version":"3.50.1"},"reference-count":40,"publisher":"Oxford University Press (OUP)","issue":"12","license":[{"start":{"date-parts":[[2023,11,22]],"date-time":"2023-11-22T00:00:00Z","timestamp":1700611200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"Portage Health Foundation Graduate Assistantship"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2023,12,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:sec>\n                  <jats:title>Motivation<\/jats:title>\n                  <jats:p>Genome-wide association studies is an essential tool for analyzing associations between phenotypes and single nucleotide polymorphisms (SNPs). Most of binary phenotypes in large biobanks are extremely unbalanced, which leads to inflated type I error rates for many widely used association tests for joint analysis of multiple phenotypes. In this article, we first propose a novel method to construct a Multi-Layer Network (MLN) using individuals with at least one case status among all phenotypes. Then, we introduce a computationally efficient community detection method to group phenotypes into disjoint clusters based on the MLN. Finally, we propose a novel approach, MLN with Omnibus (MLN-O), to jointly analyse the association between phenotypes and a SNP. MLN-O uses the score test to test the association of each merged phenotype in a cluster and a SNP, then uses the Omnibus test to obtain an overall test statistic to test the association between all phenotypes and a SNP.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Results<\/jats:title>\n                  <jats:p>We conduct extensive simulation studies to reveal that the proposed approach can control type I error rates and is more powerful than some existing methods. Meanwhile, we apply the proposed method to a real data set in the UK Biobank. Using phenotypes in Chapter XIII (Diseases of the musculoskeletal system and connective tissue) in the UK Biobank, we find that MLN-O identifies more significant SNPs than other methods we compare with.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Availability and implementation<\/jats:title>\n                  <jats:p>https:\/\/github.com\/Hongjing-Xie\/Multi-Layer-Network-with-Omnibus-MLN-O.<\/jats:p>\n               <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btad707","type":"journal-article","created":{"date-parts":[[2023,11,22]],"date-time":"2023-11-22T20:02:42Z","timestamp":1700683362000},"source":"Crossref","is-referenced-by-count":4,"title":["Joint analysis of multiple phenotypes for extremely unbalanced case\u2013control association studies using multi-layer network"],"prefix":"10.1093","volume":"39","author":[{"given":"Hongjing","family":"Xie","sequence":"first","affiliation":[{"name":"Department of Mathematical Sciences, Michigan Technological University , Houghton, MI 49931, United States"}]},{"given":"Xuewei","family":"Cao","sequence":"additional","affiliation":[{"name":"Department of Mathematical Sciences, Michigan Technological University , Houghton, MI 49931, United States"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-9478-1199","authenticated-orcid":false,"given":"Shuanglin","family":"Zhang","sequence":"additional","affiliation":[{"name":"Department of Mathematical Sciences, Michigan Technological University , Houghton, MI 49931, United States"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-9342-3269","authenticated-orcid":false,"given":"Qiuying","family":"Sha","sequence":"additional","affiliation":[{"name":"Department of Mathematical Sciences, Michigan Technological University , Houghton, MI 49931, United States"}]}],"member":"286","published-online":{"date-parts":[[2023,11,22]]},"reference":[{"key":"2023120522522952300_btad707-B1","doi-asserted-by":"crossref","first-page":"662","DOI":"10.1016\/j.ajhg.2014.03.016","article-title":"Maximizing the power of principal-component analysis of correlated phenotypes in genome-wide association studies","volume":"94","author":"Aschard","year":"2014","journal-title":"Am J Hum Genet"},{"key":"2023120522522952300_btad707-B2","doi-asserted-by":"crossref","first-page":"203","DOI":"10.1038\/s41586-018-0579-z","article-title":"The UK biobank resource with deep phenotyping and genomic data","volume":"562","author":"Bycroft","year":"2018","journal-title":"Nature"},{"key":"2023120522522952300_btad707-B3","author":"Cao","year":"2023"},{"key":"2023120522522952300_btad707-B4","doi-asserted-by":"crossref","first-page":"e1001323","DOI":"10.1371\/journal.pgen.1001323","article-title":"Differential genetic associations for systemic lupus erythematosus based on anti-dsDNA autoantibody production","volume":"7","author":"Chung","year":"2011","journal-title":"PLoS Genet"},{"key":"2023120522522952300_btad707-B5","doi-asserted-by":"crossref","first-page":"465","DOI":"10.1037\/0033-2909.115.3.465","article-title":"How the power of MANOVA can both increase and decrease as a function of the intercorrelations among the dependent variables","volume":"115","author":"Cole","year":"1994","journal-title":"Psychol Bull"},{"key":"2023120522522952300_btad707-B6","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":"2023120522522952300_btad707-B7","doi-asserted-by":"crossref","first-page":"141","DOI":"10.1093\/biomet\/80.1.141","article-title":"A likelihood-based method for analysing longitudinal binary responses","volume":"80","author":"Fitzmaurice","year":"1993","journal-title":"Biometrika"},{"key":"2023120522522952300_btad707-B8","doi-asserted-by":"crossref","first-page":"75","DOI":"10.1016\/j.physrep.2009.11.002","article-title":"Community detection in graphs","volume":"486","author":"Fortunato","year":"2010","journal-title":"Phys Rep"},{"key":"2023120522522952300_btad707-B9","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":"2023120522522952300_btad707-B10","doi-asserted-by":"crossref","first-page":"e3000836","DOI":"10.1371\/journal.pbio.3000836","article-title":"Extent and context dependence of pleiotropy revealed by high-throughput single-cell phenotyping","volume":"18","author":"Geiler-Samerotte","year":"2020","journal-title":"PLoS Biol"},{"key":"2023120522522952300_btad707-B11","doi-asserted-by":"crossref","first-page":"591","DOI":"10.1016\/j.ajhg.2012.02.008","article-title":"On sharing quantitative trait GWAS results in an era of multiple-omics data and the limits of genomic privacy","volume":"90","author":"Im","year":"2012","journal-title":"Am J Hum Genet"},{"key":"2023120522522952300_btad707-B12","doi-asserted-by":"crossref","first-page":"S33","DOI":"10.1002\/gepi.20470","article-title":"Analysis of multiple phenotypes","volume":"33","author":"Kent","year":"2009","journal-title":"Genet Epidemiol"},{"key":"2023120522522952300_btad707-B13","doi-asserted-by":"crossref","first-page":"651","DOI":"10.1002\/gepi.21931","article-title":"An adaptive association test for multiple phenotypes with GWAS summary statistics","volume":"39","author":"Kim","year":"2015","journal-title":"Genet Epidemiol"},{"key":"2023120522522952300_btad707-B14","doi-asserted-by":"crossref","first-page":"67","DOI":"10.1002\/gepi.22263","article-title":"Joint analysis of multiple phenotypes using a clustering linear combination method based on hierarchical clustering","volume":"44","author":"Li","year":"2020","journal-title":"Genet Epidemiol"},{"key":"2023120522522952300_btad707-B15","doi-asserted-by":"crossref","first-page":"13","DOI":"10.1093\/biomet\/73.1.13","article-title":"Longitudinal data analysis using generalized linear models","volume":"73","author":"Liang","year":"1986","journal-title":"Biometrika"},{"key":"2023120522522952300_btad707-B16","doi-asserted-by":"crossref","first-page":"e0276646","DOI":"10.1371\/journal.pone.0276646","article-title":"HCLC-FC: a novel statistical method for phenome-wide association studies","volume":"17","author":"Liang","year":"2022","journal-title":"PLoS ONE"},{"key":"2023120522522952300_btad707-B17","doi-asserted-by":"crossref","first-page":"344","DOI":"10.1002\/gepi.22124","article-title":"A hierarchical clustering method for dimension reduction in joint analysis of multiple phenotypes","volume":"42","author":"Liang","year":"2018","journal-title":"Genet Epidemiol"},{"key":"2023120522522952300_btad707-B18","doi-asserted-by":"crossref","first-page":"648831","DOI":"10.3389\/fgene.2021.648831","article-title":"AGNEP: an agglomerative nesting clustering algorithm for phenotypic dimension reduction in joint analysis of multiple phenotypes","volume":"12","author":"Liu","year":"2021","journal-title":"Front Genet"},{"key":"2023120522522952300_btad707-B19","first-page":"10","article-title":"Pleiotropy: one gene can affect multiple traits","volume":"1","author":"Lobo","year":"2008","journal-title":"Nat Educ"},{"key":"2023120522522952300_btad707-B20","doi-asserted-by":"crossref","first-page":"95","DOI":"10.1016\/j.physrep.2013.08.002","article-title":"Clustering and community detection in directed networks: a survey","volume":"533","author":"Malliaros","year":"2013","journal-title":"Phys Rep"},{"key":"2023120522522952300_btad707-B21","volume-title":"Policy for Sharing of Data Obtained in NIH Supported or Conducted Genome-Wide Association Studies (GWAS)","author":"National Institutes of Health","year":"2007"},{"key":"2023120522522952300_btad707-B22","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'Reilly","year":"2012","journal-title":"PLoS ONE"},{"key":"2023120522522952300_btad707-B23","doi-asserted-by":"crossref","first-page":"1081","DOI":"10.1534\/genetics.114.165035","article-title":"A powerful and adaptive association test for rare variants","volume":"197","author":"Pan","year":"2014","journal-title":"Genetics"},{"key":"2023120522522952300_btad707-B24","doi-asserted-by":"crossref","first-page":"832","DOI":"10.1016\/j.ajhg.2010.04.005","article-title":"Pooled association tests for rare variants in exon-resequencing studies","volume":"86","author":"Price","year":"2010","journal-title":"Am J Hum Genet"},{"key":"2023120522522952300_btad707-B25","doi-asserted-by":"crossref","first-page":"20","DOI":"10.1002\/gepi.21937","article-title":"USAT: a unified score-based association test for multiple phenotype-genotype analysis","volume":"40","author":"Ray","year":"2016","journal-title":"Genet Epidemiol"},{"key":"2023120522522952300_btad707-B26","doi-asserted-by":"crossref","first-page":"1373","DOI":"10.1093\/bioinformatics\/bty810","article-title":"A clustering linear combination approach to jointly analyze multiple phenotypes for GWAS","volume":"35","author":"Sha","year":"2019","journal-title":"Bioinformatics"},{"key":"2023120522522952300_btad707-B27","doi-asserted-by":"crossref","first-page":"607","DOI":"10.1016\/j.ajhg.2011.10.004","article-title":"Abundant pleiotropy in human complex diseases and traits","volume":"89","author":"Sivakumaran","year":"2011","journal-title":"Am J Hum Genet"},{"key":"2023120522522952300_btad707-B28","doi-asserted-by":"crossref","first-page":"483","DOI":"10.1038\/nrg3461","article-title":"Pleiotropy in complex traits: challenges and strategies","volume":"14","author":"Solovieff","year":"2013","journal-title":"Nat Rev Genet"},{"key":"2023120522522952300_btad707-B29","doi-asserted-by":"crossref","first-page":"767","DOI":"10.1534\/genetics.110.122549","article-title":"One hundred years of pleiotropy: a retrospective","volume":"186","author":"Stearns","year":"2010","journal-title":"Genetics"},{"key":"2023120522522952300_btad707-B30","doi-asserted-by":"crossref","first-page":"2355","DOI":"10.1056\/NEJMoa0801197","article-title":"Multiple genetic loci for bone mineral density and fractures","volume":"358","author":"Styrkarsdottir","year":"2008","journal-title":"N Engl J Med"},{"key":"2023120522522952300_btad707-B31","doi-asserted-by":"crossref","first-page":"845","DOI":"10.1093\/bioinformatics\/bts051","article-title":"A gene-based test of association using canonical correlation analysis","volume":"28","author":"Tang","year":"2012","journal-title":"Bioinformatics"},{"key":"2023120522522952300_btad707-B32","doi-asserted-by":"crossref","first-page":"82","DOI":"10.1038\/nature14962","article-title":"The UK10K project identifies rare variants in health and disease","volume":"526","author":"UK10K Consortium","year":"2015","journal-title":"Nature"},{"key":"2023120522522952300_btad707-B33","doi-asserted-by":"crossref","first-page":"e1003235","DOI":"10.1371\/journal.pgen.1003235","article-title":"TATES: efficient multivariate genotype-phenotype analysis for genome-wide association studies","volume":"9","author":"Van der Sluis","year":"2013","journal-title":"PLoS Genet"},{"key":"2023120522522952300_btad707-B34","doi-asserted-by":"crossref","first-page":"e0260911","DOI":"10.1371\/journal.pone.0260911","article-title":"A computationally efficient clustering linear combination approach to jointly analyze multiple phenotypes for GWAS","volume":"17","author":"Wang","year":"2022","journal-title":"PLoS ONE"},{"key":"2023120522522952300_btad707-B35","doi-asserted-by":"crossref","first-page":"e0201186","DOI":"10.1371\/journal.pone.0201186","article-title":"Testing an optimally weighted combination of common and\/or rare variants with multiple traits","volume":"13","author":"Wang","year":"2018","journal-title":"PLoS ONE"},{"key":"2023120522522952300_btad707-B36","doi-asserted-by":"crossref","first-page":"82","DOI":"10.1016\/j.ajhg.2011.05.029","article-title":"Rare-variant association testing for sequencing data with the sequence kernel association test","volume":"89","author":"Wu","year":"2011","journal-title":"Am J Hum Genet"},{"key":"2023120522522952300_btad707-B37","doi-asserted-by":"crossref","first-page":"223","DOI":"10.1093\/biostatistics\/4.2.223","article-title":"Combining dependent tests for linkage or association across multiple phenotypic traits","volume":"4","author":"Xu","year":"2003","journal-title":"Biostatistics"},{"key":"2023120522522952300_btad707-B38","doi-asserted-by":"crossref","first-page":"652569","DOI":"10.1155\/2012\/652569","article-title":"Methods for analyzing multivariate phenotypes in genetic association studies","volume":"2012","author":"Yang","year":"2012","journal-title":"J Probab Stat"},{"key":"2023120522522952300_btad707-B39","doi-asserted-by":"crossref","first-page":"309","DOI":"10.1016\/j.neuroimage.2014.03.061","article-title":"Testing for association with multiple traits in generalized estimation equations, with application to neuroimaging data","volume":"96","author":"Zhang","year":"2014","journal-title":"NeuroImage"},{"key":"2023120522522952300_btad707-B40","doi-asserted-by":"crossref","first-page":"144","DOI":"10.1159\/000446239","article-title":"Power comparisons of methods for joint association analysis of multiple phenotypes","volume":"80","author":"Zhu","year":"2015","journal-title":"Hum Hered"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btad707\/53686138\/btad707.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/12\/btad707\/54025726\/btad707.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/12\/btad707\/54025726\/btad707.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,12,5]],"date-time":"2023-12-05T23:00:16Z","timestamp":1701817216000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btad707\/7441501"}},"subtitle":[],"editor":[{"given":"Russell","family":"Schwartz","sequence":"additional","affiliation":[]}],"short-title":[],"issued":{"date-parts":[[2023,11,22]]},"references-count":40,"journal-issue":{"issue":"12","published-print":{"date-parts":[[2023,12,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btad707","relation":{},"ISSN":["1367-4811"],"issn-type":[{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2023,12,1]]},"published":{"date-parts":[[2023,11,22]]},"article-number":"btad707"}}