{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,27]],"date-time":"2026-02-27T17:59:00Z","timestamp":1772215140887,"version":"3.50.1"},"reference-count":49,"publisher":"Public Library of Science (PLoS)","issue":"11","license":[{"start":{"date-parts":[[2022,11,9]],"date-time":"2022-11-09T00:00:00Z","timestamp":1667952000000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["www.plosone.org"],"crossmark-restriction":false},"short-container-title":["PLoS ONE"],"abstract":"<jats:p>\n                    The emergence of genetic data coupled to longitudinal electronic medical records (EMRs) offers the possibility of phenome-wide association studies (PheWAS). In PheWAS, the whole phenome can be divided into numerous phenotypic categories according to the genetic architecture across phenotypes. Currently, statistical analyses for PheWAS are mainly univariate analyses, which test the association between one genetic variant and one phenotype at a time. In this article, we derived a novel and powerful multivariate method for PheWAS. The proposed method involves three steps. In the first step, we apply the bottom-up hierarchical clustering method to partition a large number of phenotypes into disjoint clusters within each phenotypic category. In the second step, the clustering linear combination method is used to combine test statistics within each category based on the phenotypic clusters and obtain p-values from each phenotypic category. In the third step, we propose a new false discovery rate (FDR) control approach. We perform extensive simulation studies to compare the performance of our method with that of other existing methods. The results show that our proposed method controls FDR very well and outperforms other methods we compared with. We also apply the proposed approach to a set of EMR-based phenotypes across more than 300,000 samples from the UK Biobank. We find that the proposed approach not only can well-control FDR at a nominal level but also successfully identify 1,244 significant SNPs that are reported to be associated with some phenotypes in the GWAS catalog. Our open-access tools and instructions on how to implement HCLC-FC are available at\n                    <jats:ext-link xmlns:xlink=\"http:\/\/www.w3.org\/1999\/xlink\" ext-link-type=\"uri\" xlink:href=\"https:\/\/github.com\/XiaoyuLiang\/HCLCFC\" xlink:type=\"simple\">https:\/\/github.com\/XiaoyuLiang\/HCLCFC<\/jats:ext-link>\n                    .\n                  <\/jats:p>","DOI":"10.1371\/journal.pone.0276646","type":"journal-article","created":{"date-parts":[[2022,11,9]],"date-time":"2022-11-09T13:48:26Z","timestamp":1668001706000},"page":"e0276646","update-policy":"https:\/\/doi.org\/10.1371\/journal.pone.corrections_policy","source":"Crossref","is-referenced-by-count":10,"title":["HCLC-FC: A novel statistical method for phenome-wide association studies"],"prefix":"10.1371","volume":"17","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-7796-2441","authenticated-orcid":true,"given":"Xiaoyu","family":"Liang","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2136-0964","authenticated-orcid":true,"given":"Xuewei","family":"Cao","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-9342-3269","authenticated-orcid":true,"given":"Qiuying","family":"Sha","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-9478-1199","authenticated-orcid":true,"given":"Shuanglin","family":"Zhang","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"340","published-online":{"date-parts":[[2022,11,9]]},"reference":[{"issue":"12","key":"pone.0276646.ref001","doi-asserted-by":"crossref","first-page":"e1002822","DOI":"10.1371\/journal.pcbi.1002822","article-title":"Chapter 11: Genome-wide association studies","volume":"8","author":"WS Bush","year":"2012","journal-title":"PLoS Comput Biol"},{"issue":"10","key":"pone.0276646.ref002","doi-asserted-by":"crossref","first-page":"e0140348","DOI":"10.1371\/journal.pone.0140348","article-title":"Association Tests of Multiple Phenotypes: ATeMP","volume":"10","author":"X Guo","year":"2015","journal-title":"PLoS One"},{"issue":"10","key":"pone.0276646.ref003","doi-asserted-by":"crossref","first-page":"1294","DOI":"10.1093\/bioinformatics\/btm108","article-title":"GenABEL: an R library for genome-wide association analysis","volume":"23","author":"YS Aulchenko","year":"2007","journal-title":"Bioinformatics"},{"key":"pone.0276646.ref004","doi-asserted-by":"crossref","first-page":"134","DOI":"10.1186\/1471-2105-11-134","article-title":"ProbABEL package for genome-wide association analysis of imputed data","volume":"11","author":"YS Aulchenko","year":"2010","journal-title":"BMC Bioinformatics"},{"key":"pone.0276646.ref005","doi-asserted-by":"crossref","first-page":"387","DOI":"10.1146\/annurev.genom.9.081307.164242","article-title":"Genotype imputation","volume":"10","author":"Y Li","year":"2009","journal-title":"Annu Rev Genomics Hum Genet"},{"issue":"8","key":"pone.0276646.ref006","doi-asserted-by":"crossref","first-page":"816","DOI":"10.1002\/gepi.20533","article-title":"MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes","volume":"34","author":"Y Li","year":"2010","journal-title":"Genet Epidemiol"},{"issue":"10","key":"pone.0276646.ref007","doi-asserted-by":"crossref","first-page":"833","DOI":"10.1038\/nmeth.1681","article-title":"FaST linear mixed models for genome-wide association studies","volume":"8","author":"C Lippert","year":"2011","journal-title":"Nat Methods"},{"issue":"7","key":"pone.0276646.ref008","doi-asserted-by":"crossref","first-page":"906","DOI":"10.1038\/ng2088","article-title":"A new multipoint method for genome-wide association studies by imputation of genotypes","volume":"39","author":"J Marchini","year":"2007","journal-title":"Nat Genet"},{"issue":"3","key":"pone.0276646.ref009","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":"S Purcell","year":"2007","journal-title":"Am J Hum Genet"},{"issue":"2","key":"pone.0276646.ref010","first-page":"983","article-title":"Large-Scale Simultaneous Testing of Cross-Covariance Matrices with Applications to PheWAS","volume":"29","author":"T Cai","year":"2019","journal-title":"Statistica Sinica"},{"issue":"5","key":"pone.0276646.ref011","doi-asserted-by":"crossref","first-page":"410","DOI":"10.1002\/gepi.20589","article-title":"The use of phenome-wide association studies (PheWAS) for exploration of novel genotype-phenotype relationships and pleiotropy discovery","volume":"35","author":"SA Pendergrass","year":"2011","journal-title":"Genet Epidemiol"},{"issue":"2","key":"pone.0276646.ref012","doi-asserted-by":"crossref","first-page":"157","DOI":"10.1111\/imm.12195","article-title":"The challenges, advantages and future of phenome-wide association studies","volume":"141","author":"SJ Hebbring","year":"2014","journal-title":"Immunology"},{"issue":"4","key":"pone.0276646.ref013","doi-asserted-by":"crossref","first-page":"444","DOI":"10.1007\/s11684-012-0225-3","article-title":"Genome-wide association studies: inherent limitations and future challenges","volume":"6","author":"Y Du","year":"2012","journal-title":"Front Med"},{"issue":"4","key":"pone.0276646.ref014","doi-asserted-by":"crossref","first-page":"529","DOI":"10.1016\/j.ajhg.2011.09.008","article-title":"Variants near FOXE1 are associated with hypothyroidism and other thyroid conditions: using electronic medical records for genome- and phenome-wide studies","volume":"89","author":"JC Denny","year":"2011","journal-title":"Am J Hum Genet"},{"issue":"9","key":"pone.0276646.ref015","doi-asserted-by":"crossref","first-page":"1205","DOI":"10.1093\/bioinformatics\/btq126","article-title":"PheWAS: demonstrating the feasibility of a phenome-wide scan to discover gene-disease associations","volume":"26","author":"JC Denny","year":"2010","journal-title":"Bioinformatics"},{"issue":"4","key":"pone.0276646.ref016","doi-asserted-by":"crossref","first-page":"560","DOI":"10.1016\/j.ajhg.2010.03.003","article-title":"Robust replication of genotype-phenotype associations across multiple diseases in an electronic medical record","volume":"86","author":"MD Ritchie","year":"2010","journal-title":"Am J Hum Genet"},{"issue":"1","key":"pone.0276646.ref017","doi-asserted-by":"crossref","first-page":"e1003087","DOI":"10.1371\/journal.pgen.1003087","article-title":"Phenome-wide association study (PheWAS) for detection of pleiotropy within the Population Architecture using Genomics and Epidemiology (PAGE) Network","volume":"9","author":"SA Pendergrass","year":"2013","journal-title":"PLoS Genet"},{"issue":"1","key":"pone.0276646.ref018","doi-asserted-by":"crossref","first-page":"55","DOI":"10.1016\/j.ajhg.2018.11.006","article-title":"Human-Disease Phenotype Map Derived from PheWAS across 38,682 Individuals","volume":"104","author":"A Verma","year":"2019","journal-title":"Am J Hum Genet"},{"issue":"2","key":"pone.0276646.ref019","first-page":"746","article-title":"An introduction to multiplicity issues in clinical trials: the what, why, when and how","volume":"46","author":"G Li","year":"2017","journal-title":"Int J Epidemiol"},{"issue":"3","key":"pone.0276646.ref020","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":"DA Cole","year":"1994","journal-title":"Psychological bulletin"},{"issue":"5","key":"pone.0276646.ref021","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":"PF O\u2019Reilly","year":"2012","journal-title":"PLoS One"},{"issue":"1","key":"pone.0276646.ref022","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":"S van der Sluis","year":"2013","journal-title":"PLoS Genet"},{"key":"pone.0276646.ref023","doi-asserted-by":"crossref","first-page":"34323","DOI":"10.1038\/srep34323","article-title":"An Adaptive Fisher\u2019s Combination Method for Joint Analysis of Multiple Phenotypes in Association Studies","volume":"6","author":"X Liang","year":"2016","journal-title":"Sci Rep"},{"issue":"8","key":"pone.0276646.ref024","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":"Q Sha","year":"2019","journal-title":"Bioinformatics"},{"issue":"1","key":"pone.0276646.ref025","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":"X Li","year":"2020","journal-title":"Genet Epidemiol"},{"issue":"4","key":"pone.0276646.ref026","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":"X Liang","year":"2018","journal-title":"Genet Epidemiol"},{"issue":"8","key":"pone.0276646.ref027","doi-asserted-by":"crossref","first-page":"904","DOI":"10.1038\/ng1847","article-title":"Principal components analysis corrects for stratification in genome-wide association studies","volume":"38","author":"AL Price","year":"2006","journal-title":"Nat Genet"},{"issue":"6","key":"pone.0276646.ref028","doi-asserted-by":"crossref","first-page":"561","DOI":"10.1002\/gepi.21649","article-title":"Detecting association of rare and common variants by testing an optimally weighted combination of variants","volume":"36","author":"Q Sha","year":"2012","journal-title":"Genet Epidemiol"},{"issue":"11","key":"pone.0276646.ref029","doi-asserted-by":"crossref","first-page":"1835","DOI":"10.1016\/j.jspi.2013.05.019","article-title":"Correlated variables in regression: clustering and sparse estimation","volume":"143","author":"P B\u00fchlmann","year":"2013","journal-title":"Journal of Statistical Planning and Inference"},{"issue":"6","key":"pone.0276646.ref030","doi-asserted-by":"crossref","first-page":"2048","DOI":"10.1038\/s41380-020-0670-3","article-title":"GWAS significance thresholds for deep phenotyping studies can depend upon minor allele frequencies and sample size","volume":"26","author":"H Asif","year":"2021","journal-title":"Mol Psychiatry"},{"issue":"4","key":"pone.0276646.ref031","doi-asserted-by":"crossref","first-page":"765","DOI":"10.1086\/383251","article-title":"A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other","volume":"74","author":"DR Nyholt","year":"2004","journal-title":"Am J Hum Genet"},{"issue":"3","key":"pone.0276646.ref032","doi-asserted-by":"crossref","first-page":"221","DOI":"10.1038\/sj.hdy.6800717","article-title":"Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix","volume":"95","author":"J Li","year":"2005","journal-title":"Heredity (Edinb)"},{"issue":"318","key":"pone.0276646.ref033","first-page":"626","article-title":"Rectangular confidence regions for the means of multivariate normal distributions","volume":"62","author":"Z. \u0160id\u00e1k","year":"1967","journal-title":"Journal of the American Statistical Association"},{"issue":"7726","key":"pone.0276646.ref034","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":"C Bycroft","year":"2018","journal-title":"Nature"},{"issue":"3","key":"pone.0276646.ref035","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":"C Sudlow","year":"2015","journal-title":"PLoS Med"},{"key":"pone.0276646.ref036","first-page":"565903","article-title":"Hierarchical clustering of gene-level association statistics reveals shared and differential genetic architecture among traits in the UK Biobank","author":"MR McGuirl","year":"2019","journal-title":"bioRxiv"},{"issue":"1","key":"pone.0276646.ref037","doi-asserted-by":"crossref","DOI":"10.1186\/s13742-015-0047-8","article-title":"Second-generation PLINK: rising to the challenge of larger and richer datasets","volume":"4","author":"CC Chang","year":"2015","journal-title":"Gigascience"},{"key":"pone.0276646.ref038","doi-asserted-by":"crossref","first-page":"353","DOI":"10.1146\/annurev-genom-090314-024956","article-title":"Phenome-wide association studies as a tool to advance precision medicine","volume":"17","author":"JC Denny","year":"2016","journal-title":"Annual review of genomics and human genetics"},{"issue":"1","key":"pone.0276646.ref039","doi-asserted-by":"crossref","first-page":"37","DOI":"10.1016\/j.ajhg.2017.05.014","article-title":"A Fast and Accurate Algorithm to Test for Binary Phenotypes and Its Application to PheWAS","volume":"101","author":"R Dey","year":"2017","journal-title":"Am J Hum Genet"},{"key":"pone.0276646.ref040","doi-asserted-by":"crossref","DOI":"10.1007\/978-1-4419-7338-2","volume-title":"The fundamentals of modern statistical genetics","author":"NM Laird","year":"2011"},{"issue":"8","key":"pone.0276646.ref041","doi-asserted-by":"crossref","first-page":"467","DOI":"10.1038\/s41576-019-0127-1","article-title":"Benefits and limitations of genome-wide association studies","volume":"20","author":"V Tam","year":"2019","journal-title":"Nature Reviews Genetics"},{"key":"pone.0276646.ref042","doi-asserted-by":"crossref","first-page":"343","DOI":"10.1146\/annurev.genom.7.080505.115721","article-title":"The uneasy ethical and legal underpinnings of large-scale genomic biobanks","volume":"8","author":"HT Greely","year":"2007","journal-title":"Annu Rev Genomics Hum Genet"},{"issue":"9","key":"pone.0276646.ref043","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":"W Zhou","year":"2018","journal-title":"Nat Genet"},{"issue":"9","key":"pone.0276646.ref044","doi-asserted-by":"crossref","first-page":"1311","DOI":"10.1038\/ng.3926","article-title":"Bayesian analysis of genetic association across tree-structured routine healthcare data in the UK Biobank","volume":"49","author":"A Cortes","year":"2017","journal-title":"Nat Genet"},{"issue":"9","key":"pone.0276646.ref045","doi-asserted-by":"crossref","first-page":"1295","DOI":"10.1038\/ng.3946","article-title":"Reaching for the next branch on the biobank tree of knowledge","volume":"49","author":"NJ Cox","year":"2017","journal-title":"Nat Genet"},{"issue":"4","key":"pone.0276646.ref046","doi-asserted-by":"crossref","first-page":"321","DOI":"10.1007\/s40471-017-0127-7","article-title":"Current Scope and Challenges in Phenome-Wide Association Studies","volume":"4","author":"A Verma","year":"2017","journal-title":"Curr Epidemiol Rep"},{"issue":"7","key":"pone.0276646.ref047","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":"Z Wang","year":"2018","journal-title":"PLoS One"},{"issue":"3","key":"pone.0276646.ref048","doi-asserted-by":"crossref","first-page":"291","DOI":"10.1038\/ng.3211","article-title":"LD Score regression distinguishes confounding from polygenicity in genome-wide association studies","volume":"47","author":"Schizophrenia Working Group of the Psychiatric Genomics C","year":"2015","journal-title":"Nat Genet"},{"issue":"11","key":"pone.0276646.ref049","doi-asserted-by":"crossref","first-page":"1236","DOI":"10.1038\/ng.3406","article-title":"An atlas of genetic correlations across human diseases and traits","volume":"47","author":"B Bulik-Sullivan","year":"2015","journal-title":"Nat Genet"}],"container-title":["PLOS ONE"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pone.0276646","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2022,11,9]],"date-time":"2022-11-09T13:49:15Z","timestamp":1668001755000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pone.0276646"}},"subtitle":[],"editor":[{"given":"Chunyu","family":"Liu","sequence":"first","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2022,11,9]]},"references-count":49,"journal-issue":{"issue":"11","published-online":{"date-parts":[[2022,11,9]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pone.0276646","relation":{"has-preprint":[{"id-type":"doi","id":"10.1101\/2022.03.14.484203","asserted-by":"object"}]},"ISSN":["1932-6203"],"issn-type":[{"value":"1932-6203","type":"electronic"}],"subject":[],"published":{"date-parts":[[2022,11,9]]}}}