1
|
Wang M, Zhang S, Sha Q. A computationally efficient clustering linear combination approach to jointly analyze multiple phenotypes for GWAS. PLoS One 2022; 17:e0260911. [PMID: 35482827 PMCID: PMC9049312 DOI: 10.1371/journal.pone.0260911] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2021] [Accepted: 04/13/2022] [Indexed: 11/18/2022] Open
Abstract
There has been an increasing interest in joint analysis of multiple phenotypes in genome-wide association studies (GWAS) because jointly analyzing multiple phenotypes may increase statistical power to detect genetic variants associated with complex diseases or traits. Recently, many statistical methods have been developed for joint analysis of multiple phenotypes in genetic association studies, including the Clustering Linear Combination (CLC) method. The CLC method works particularly well with phenotypes that have natural groupings, but due to the unknown number of clusters for a given data, the final test statistic of CLC method is the minimum p-value among all p-values of the CLC test statistics obtained from each possible number of clusters. Therefore, a simulation procedure needs to be used to evaluate the p-value of the final test statistic. This makes the CLC method computationally demanding. We develop a new method called computationally efficient CLC (ceCLC) to test the association between multiple phenotypes and a genetic variant. Instead of using the minimum p-value as the test statistic in the CLC method, ceCLC uses the Cauchy combination test to combine all p-values of the CLC test statistics obtained from each possible number of clusters. The test statistic of ceCLC approximately follows a standard Cauchy distribution, so the p-value can be obtained from the cumulative density function without the need for the simulation procedure. Through extensive simulation studies and application on the COPDGene data, the results demonstrate that the type I error rates of ceCLC are effectively controlled in different simulation settings and ceCLC either outperforms all other methods or has statistical power that is very close to the most powerful method with which it has been compared.
Collapse
Affiliation(s)
- Meida Wang
- Mathematical Sciences, Michigan Technological University, Houghton, MI, United States of America
| | - Shuanglin Zhang
- Mathematical Sciences, Michigan Technological University, Houghton, MI, United States of America
| | - Qiuying Sha
- Mathematical Sciences, Michigan Technological University, Houghton, MI, United States of America
| |
Collapse
|
2
|
Associating Multivariate Traits with Genetic Variants Using Collapsing and Kernel Methods with Pedigree- or Population-Based Studies. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2021; 2021:8812282. [PMID: 33628328 PMCID: PMC7889379 DOI: 10.1155/2021/8812282] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/06/2020] [Revised: 01/02/2021] [Accepted: 01/08/2021] [Indexed: 11/18/2022]
Abstract
In genetic association analysis, several relevant phenotypes or multivariate traits with different types of components are usually collected to study complex or multifactorial diseases. Over the past few years, jointly testing for association between multivariate traits and multiple genetic variants has become more popular because it can increase statistical power to identify causal genes in pedigree- or population-based studies. However, most of the existing methods mainly focus on testing genetic variants associated with multiple continuous phenotypes. In this investigation, we develop a framework for identifying the pleiotropic effects of genetic variants on multivariate traits by using collapsing and kernel methods with pedigree- or population-structured data. The proposed framework is applicable to the burden test, the kernel test, and the omnibus test for autosomes and the X chromosome. The proposed multivariate trait association methods can accommodate continuous phenotypes or binary phenotypes and further can adjust for covariates. Simulation studies show that the performance of our methods is satisfactory with respect to the empirical type I error rates and power rates in comparison with the existing methods.
Collapse
|
3
|
Chen L, Zhou Y. A fast and powerful aggregated Cauchy association test for joint analysis of multiple phenotypes. Genes Genomics 2021; 43:69-77. [PMID: 33432394 DOI: 10.1007/s13258-020-01034-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2020] [Accepted: 12/23/2020] [Indexed: 11/27/2022]
Abstract
BACKGROUND Pleiotropy is a widespread phenomenon in complex human diseases. Jointly analyzing multiple phenotypes can improve power performance of detecting genetic variants and uncover the underlying genetic mechanism. OBJECTIVE This study aims to detect the association between genetic variants in a genomic region and multiple phenotypes. METHODS We develop the aggregated Cauchy association test to detect the association between rare variants in a genomic region and multiple phenotypes (abbreviated as "Multi-ACAT"). Multi-ACAT first detects the association between each rare variant and multiple phenotypes based on reverse regression and obtains variant-level p-values, then takes linear combination of transformed p-values as the test statistic which approximately follows Cauchy distribution under the null hypothesis. RESULTS Extensive simulation studies show that when the proportion of causal variants in a genomic region is extremely small, Multi-ACAT is more powerful than the other several methods and is robust to bi-directional effects of causal variants. Finally, we illustrate our proposed method by analyzing two phenotypes [systolic blood pressure (SBP) and diastolic blood pressure (DBP)] from Genetic Analysis Workshop 19 (GAW19). CONCLUSION The Multi-ACAT computes extremely fast, does not consider complex distributions of multiple correlated phenotypes, and can be applied to the case with noise phenotypes.
Collapse
Affiliation(s)
- Lili Chen
- School of Mathematical Sciences, Heilongjiang University, No. 74 Xuefu Road, Nangang District, Harbin, 150080, People's Republic of China
| | - Yajing Zhou
- School of Mathematical Sciences, Heilongjiang University, No. 74 Xuefu Road, Nangang District, Harbin, 150080, People's Republic of China.
| |
Collapse
|
4
|
Konigorski S, Yilmaz YE, Janke J, Bergmann MM, Boeing H, Pischon T. Powerful rare variant association testing in a copula-based joint analysis of multiple phenotypes. Genet Epidemiol 2019; 44:26-40. [PMID: 31732979 DOI: 10.1002/gepi.22265] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2019] [Revised: 08/13/2019] [Accepted: 09/16/2019] [Indexed: 12/16/2022]
Abstract
In genetic association studies of rare variants, the low power of association tests is one of the main challenges. In this study, we propose a new single-marker association test called C-JAMP (Copula-based Joint Analysis of Multiple Phenotypes), which is based on a joint model of multiple phenotypes given genetic markers and other covariates. We evaluated its performance and compared its empirical type I error and power with existing univariate and multivariate single-marker and multi-marker rare-variant tests in extensive simulation studies. C-JAMP yielded unbiased genetic effect estimates and valid type I errors with an adjusted test statistic. When strongly dependent traits were jointly analyzed, C-JAMP had the highest power in all scenarios except when a high percentage of variants were causal with moderate/small effect sizes. When traits with weak or moderate dependence were analyzed, whether C-JAMP or competing approaches had higher power depended on the effect size. When C-JAMP was applied with a misspecified copula function, it still achieved high power in some of the scenarios considered. In a real-data application, we analyzed sequencing data using C-JAMP and performed the first genome-wide association studies of high-molecular-weight and medium-molecular-weight adiponectin plasma concentrations. C-JAMP identified 20 rare variants with p-values smaller than 10-5 , while all other tests resulted in the identification of fewer variants with higher p-values. In summary, the results indicate that C-JAMP is a powerful, flexible, and robust method for association studies, and we identified novel candidate markers for adiponectin. C-JAMP is implemented as an R package and freely available from https://cran.r-project.org/package=CJAMP.
Collapse
Affiliation(s)
- Stefan Konigorski
- Molecular Epidemiology Research Group, Max Delbrück Center (MDC) for Molecular Medicine in the Helmholtz Association, Berlin, Germany.,Digital Health and Machine Learning Research Group, Hasso Plattner Institute for Digital Engineering, Potsdam, Germany
| | - Yildiz E Yilmaz
- Department of Mathematics and Statistics, Memorial University of Newfoundland, St. John's, NL, Canada.,Discipline of Genetics, Faculty of Medicine, Memorial University of Newfoundland, St. John's, NL, Canada.,Discipline of Medicine, Faculty of Medicine, Memorial University of Newfoundland, St. John's, NL, Canada
| | - Jürgen Janke
- Molecular Epidemiology Research Group, Max Delbrück Center (MDC) for Molecular Medicine in the Helmholtz Association, Berlin, Germany
| | - Manuela M Bergmann
- Department of Epidemiology, German Institute of Human Nutrition Potsdam-Rehbrücke (DIfE), Nuthetal, Germany
| | - Heiner Boeing
- Department of Epidemiology, German Institute of Human Nutrition Potsdam-Rehbrücke (DIfE), Nuthetal, Germany
| | - Tobias Pischon
- Molecular Epidemiology Research Group, Max Delbrück Center (MDC) for Molecular Medicine in the Helmholtz Association, Berlin, Germany.,Charité-Universitätsmedizin Berlin, Berlin, Germany.,DZHK (German Center for Cardiovascular Research), partner site Berlin, Berlin, Germany
| |
Collapse
|
5
|
Onogi A. Comparison of F-tests for Univariate and Multivariate Mixed-Effect Models in Genome-Wide Association Mapping. Front Genet 2019; 10:30. [PMID: 30778369 PMCID: PMC6369166 DOI: 10.3389/fgene.2019.00030] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2018] [Accepted: 01/17/2019] [Indexed: 01/24/2023] Open
Abstract
Genome-wide association mapping (GWA) has been widely applied to a variety of species to identify genomic regions responsible for quantitative traits. The use of multivariate information could enhance the detection power of GWA. Although mixed-effect models are frequently used for GWA, the utility of F-tests for multivariate mixed-effect models is not well-recognized. Thus, we compared the F-tests for univariate and multivariate mixed-effect models with simulations. The superiority of the multivariate F-test over the univariate test varied depending on three parameters: phenotypic correlation between variates (r), relative size of quantitative trait locus effects between variates (ad), and missing proportion of phenotypic records (mprop). Simulation results showed that, when mprop was low, the multivariate F-test outperformed the univariate test as r and ad differ, and as mprop increased, the multivariate F-test outperformed as ad increased. These observations were consistent with results of the analytical evaluation of the F-value. When mprop was at the maximum, i.e., when no individual had phenotypic values for multiple variates, as in the case of meta-analysis, the multivariate F-test gained more detection power as ad increased. Although using multivariate information in mixed-effect model contexts did not always ensure more detection power than with univariate tests, the multivariate F-test will be a method applied when multivariate data are available because it does not show inflation of signals and could lead to new findings.
Collapse
Affiliation(s)
- Akio Onogi
- Institute of Crop Science, National Agriculture and Food Research Organization, Tsukuba, Japan.,Japan Science and Technology Agency PRESTO, Kawaguchi, Japan
| |
Collapse
|
6
|
Wang Z, Sha Q, Fang S, Zhang K, Zhang S. Testing an optimally weighted combination of common and/or rare variants with multiple traits. PLoS One 2018; 13:e0201186. [PMID: 30048520 PMCID: PMC6062080 DOI: 10.1371/journal.pone.0201186] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2018] [Accepted: 07/10/2018] [Indexed: 12/25/2022] Open
Abstract
Recently, joint analysis of multiple traits has become popular because it can increase statistical power to identify genetic variants associated with complex diseases. In addition, there is increasing evidence indicating that pleiotropy is a widespread phenomenon in complex diseases. Currently, most of existing methods test the association between multiple traits and a single genetic variant. However, these methods by analyzing one variant at a time may not be ideal for rare variant association studies because of the allelic heterogeneity as well as the extreme rarity of rare variants. In this article, we developed a statistical method by testing an optimally weighted combination of variants with multiple traits (TOWmuT) to test the association between multiple traits and a weighted combination of variants (rare and/or common) in a genomic region. TOWmuT is robust to the directions of effects of causal variants and is applicable to different types of traits. Using extensive simulation studies, we compared the performance of TOWmuT with the following five existing methods: gene association with multiple traits (GAMuT), multiple sequence kernel association test (MSKAT), adaptive weighting reverse regression (AWRR), single-TOW, and MANOVA. Our results showed that, in all of the simulation scenarios, TOWmuT has correct type I error rates and is consistently more powerful than the other five tests. We also illustrated the usefulness of TOWmuT by analyzing a whole-genome genotyping data from a lung function study.
Collapse
Affiliation(s)
- Zhenchuan Wang
- Department of Mathematical Sciences, Michigan Technological University, Houghton, Michigan, United States of America
| | - Qiuying Sha
- Department of Mathematical Sciences, Michigan Technological University, Houghton, Michigan, United States of America
| | - Shurong Fang
- Department of Mathematics and Computer Science, John Carroll University, University Heights, Ohio, United States of America
| | - Kui Zhang
- Department of Mathematical Sciences, Michigan Technological University, Houghton, Michigan, United States of America
| | - Shuanglin Zhang
- Department of Mathematical Sciences, Michigan Technological University, Houghton, Michigan, United States of America
| |
Collapse
|
7
|
Hackinger S, Zeggini E. Statistical methods to detect pleiotropy in human complex traits. Open Biol 2018; 7:rsob.170125. [PMID: 29093210 PMCID: PMC5717338 DOI: 10.1098/rsob.170125] [Citation(s) in RCA: 78] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2017] [Accepted: 09/29/2017] [Indexed: 12/13/2022] Open
Abstract
In recent years pleiotropy, the phenomenon of one genetic locus influencing several traits, has become a widely researched field in human genetics. With the increasing availability of genome-wide association study summary statistics, as well as the establishment of deeply phenotyped sample collections, it is now possible to systematically assess the genetic overlap between multiple traits and diseases. In addition to increasing power to detect associated variants, multi-trait methods can also aid our understanding of how different disorders are aetiologically linked by highlighting relevant biological pathways. A plethora of available tools to perform such analyses exists, each with their own advantages and limitations. In this review, we outline some of the currently available methods to conduct multi-trait analyses. First, we briefly introduce the concept of pleiotropy and outline the current landscape of pleiotropy research in human genetics; second, we describe analytical considerations and analysis methods; finally, we discuss future directions for the field.
Collapse
|
8
|
Chen L, Wang Y, Zhou Y. Association analysis of multiple traits by an approach of combining P values. J Genet 2018; 97:79-85. [PMID: 29666327] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Increasing evidence shows that one variant can affect multiple traits, which is a widespread phenomenon in complex diseases. Joint analysis of multiple traits can increase statistical power of association analysis and uncover the underlying genetic mechanism. Although there are many statistical methods to analyse multiple traits, most of these methods are usually suitable for detecting common variants associated with multiple traits. However, because of low minor allele frequency of rare variant, these methods are not optimal for rare variant association analysis. In this paper, we extend an adaptive combination of P values method (termed ADA) for single trait to test association between multiple traits and rare variants in the given region. For a given region, we use reverse regression model to test each rare variant associated with multiple traits and obtain the P value of single-variant test. Further, we take the weighted combination of these P values as the test statistic. Extensive simulation studies show that our approach is more powerful than several other comparison methods in most cases and is robust to the inclusion of a high proportion of neutral variants and the different directions of effects of causal variants.
Collapse
Affiliation(s)
- Lili Chen
- Department of Mathematics, School of Sciences, Harbin Institute of Technology, Harbin 150001, People's Republic of China.
| | | | | |
Collapse
|
9
|
Chen L, Wang Y, Zhou Y. Association analysis of multiple traits by an approach of combining
$$P$$
P
values. J Genet 2018. [DOI: 10.1007/s12041-018-0885-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
10
|
Association analysis of rare and common variants with multiple traits based on variable reduction method. Genet Res (Camb) 2018; 100:e2. [PMID: 29386084 DOI: 10.1017/s0016672317000052] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
Pleiotropy, the effect of one variant on multiple traits, is widespread in complex diseases. Joint analysis of multiple traits can improve statistical power to detect genetic variants and uncover the underlying genetic mechanism. Currently, a large number of existing methods target one common variant or only rare variants. Increasing evidence shows that complex diseases are caused by common and rare variants. Here we propose a region-based method to test both rare and common variant associated multiple traits based on variable reduction method (abbreviated as MULVR). However, in the presence of noise traits, the MULVR method may lose power, so we propose the MULVR-O method, which jointly analyses the optimal number of traits associated with genetic variants by the MULVR method, to guard against the effect of noise traits. Extensive simulation studies show that our proposed method (MULVR-O) is applied to not only multiple quantitative traits but also qualitative traits, and is more powerful than several other comparison methods in most scenarios. An application to the two genes (SHBG and CHRM3) and two phenotypes (systolic blood pressure and diastolic blood pressure) from the GAW19 dataset illustrates that our proposed methods (MULVR and MULVR-O) are feasible and efficient as a region-based method.
Collapse
|
11
|
Zhu H, Zhang S, Sha Q. A novel method to test associations between a weighted combination of phenotypes and genetic variants. PLoS One 2018; 13:e0190788. [PMID: 29329304 PMCID: PMC5766098 DOI: 10.1371/journal.pone.0190788] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2017] [Accepted: 12/20/2017] [Indexed: 11/18/2022] Open
Abstract
Many complex diseases like diabetes, hypertension, metabolic syndrome, et cetera, are measured by multiple correlated phenotypes. However, most genome-wide association studies (GWAS) focus on one phenotype of interest or study multiple phenotypes separately for identifying genetic variants associated with complex diseases. Analyzing one phenotype or the related phenotypes separately may lose power due to ignoring the information obtained by combining phenotypes, such as the correlation between phenotypes. In order to increase statistical power to detect genetic variants associated with complex diseases, we develop a novel method to test a weighted combination of multiple phenotypes (WCmulP). We perform extensive simulation studies as well as real data (COPDGene) analysis to evaluate the performance of the proposed method. Our simulation results show that WCmulP has correct type I error rates and is either the most powerful test or comparable to the most powerful test among the methods we compared. WCmulP also has an outstanding performance for identifying single-nucleotide polymorphisms (SNPs) associated with COPD-related phenotypes.
Collapse
Affiliation(s)
- Huanhuan Zhu
- Department of Mathematical Sciences, Michigan Technological University, Houghton, Michigan, United States of America
| | - Shuanglin Zhang
- Department of Mathematical Sciences, Michigan Technological University, Houghton, Michigan, United States of America
| | - Qiuying Sha
- Department of Mathematical Sciences, Michigan Technological University, Houghton, Michigan, United States of America
- * E-mail:
| |
Collapse
|
12
|
Surowy H, Varga D, Burwinkel B, Marmé F, Sohn C, Luedeke M, Rinckleb A, Maier C, Deissler H, Volcic M, Wiesmüller L, Hasenburg A, Klar M, Hoegel J, Vogel W. A low-frequency haplotype spanning SLX4/FANCP constitutes a new risk locus for early-onset breast cancer (<60 years) and is associated with reduced DNA repair capacity. Int J Cancer 2017; 142:757-768. [PMID: 29044504 DOI: 10.1002/ijc.31105] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Revised: 08/31/2017] [Accepted: 09/25/2017] [Indexed: 01/21/2023]
Abstract
Only a fraction of breast cancer (BC) cases can be yet explained by mutations in genes or genomic variants discovered in linkage, genome-wide association and sequencing studies. The known genes entailing medium or high risk for BC are strongly enriched for a function in DNA double strand repair. Thus, aiming at identifying low frequency variants conferring an intermediate risk, we here investigated 17 variants (MAF: 0.01-0.1) in 10 candidate genes involved in DNA repair or cell cycle control. In an exploration cohort of 437 cases and 1189 controls, we show the variant rs3810813 in the SLX4/FANCP gene to be significantly associated with both BC (≤60 years; OR = 2.6(1.6-3.9), p = 1.6E-05) and decreased DNA repair capacity (≤60 years; beta = 37.8(17.9-57.8), p = 5.3E-4). BC association was confirmed in a verification cohort (N = 2441). Both associations were absent from cases diagnosed >60 years and stronger the earlier the diagnosis. By imputation we show that rs3810813 tags a haplotype with 5 additional variants with the same allele frequency (R2 > 0.9), and a pattern of association very similar for both phenotypes (cases <60 years, p < 0.001, the Bonferroni threshold derived from unlinked variants in the region). In young cases (≤60 years) carrying the risk haplotype, micronucleus test results are predictive for BC (AUC > 0.9). Our findings propose a risk variant with high penetrance on the haplotype spanning SLX4/FANCP to be functionally associated to BC predisposition via decreased repair capacity and suggest this variant is carried by a fraction of these haplotypes that is enriched in early onset BC cases.
Collapse
Affiliation(s)
- Harald Surowy
- Institute of Human Genetics, University of Ulm, Albert-Einstein-Allee 11, Ulm, D-89081, Germany.,Division of Molecular Biology of Breast Cancer, Department of Obstetrics and Gynecology, University of Heidelberg, Im Neuenheimer Feld 440, Heidelberg, D-69120, Germany.,Molecular Epidemiology, C080, German Cancer Research Center, Im Neuenheimer Feld 581, Heidelberg, D-69120, Germany.,Harald Surowy's current address is: Institute of Human Genetics, University of Duesseldorf, Universitaetsstr. 1, Duesseldorf, 40225, Germany
| | - Dominic Varga
- Department of Obstetrics and Gynecology, Ulm University, Prittwitzstr. 43, Ulm, D-89075, Germany.,Dominic Varga's current address is: Department of Gynecology, Donauklinik, Krankenhausstr 11, Neu-Ulm, 89231, Germany
| | - Barbara Burwinkel
- Division of Molecular Biology of Breast Cancer, Department of Obstetrics and Gynecology, University of Heidelberg, Im Neuenheimer Feld 440, Heidelberg, D-69120, Germany.,Molecular Epidemiology, C080, German Cancer Research Center, Im Neuenheimer Feld 581, Heidelberg, D-69120, Germany
| | - Frederik Marmé
- Division of Molecular Biology of Breast Cancer, Department of Obstetrics and Gynecology, University of Heidelberg, Im Neuenheimer Feld 440, Heidelberg, D-69120, Germany.,National Centre for Tumor Diseases, Heidelberg, D-69120, Germany
| | - Christof Sohn
- Division of Molecular Biology of Breast Cancer, Department of Obstetrics and Gynecology, University of Heidelberg, Im Neuenheimer Feld 440, Heidelberg, D-69120, Germany
| | - Manuel Luedeke
- Institute of Human Genetics, University of Ulm, Albert-Einstein-Allee 11, Ulm, D-89081, Germany
| | - Antje Rinckleb
- Institute of Human Genetics, University of Ulm, Albert-Einstein-Allee 11, Ulm, D-89081, Germany
| | - Christiane Maier
- Institute of Human Genetics, University of Ulm, Albert-Einstein-Allee 11, Ulm, D-89081, Germany
| | - Helmut Deissler
- Department of Obstetrics and Gynecology, Ulm University, Prittwitzstr. 43, Ulm, D-89075, Germany.,Helmut Deissler's current address is: HD/U Biomedical Services, Im Wiblinger Hart 62, Ulm, 89079, Germany
| | - Meta Volcic
- Department of Obstetrics and Gynecology, Ulm University, Prittwitzstr. 43, Ulm, D-89075, Germany
| | - Lisa Wiesmüller
- Department of Obstetrics and Gynecology, Ulm University, Prittwitzstr. 43, Ulm, D-89075, Germany
| | - Annette Hasenburg
- Department of Obstetrics and Gynecology, University of Freiburg, Hugstetter Straße 55, Freiburg, D-79106, Germany.,Annete Hasenburg's current address is: Klinik und Poliklinik für Geburtshilfe und Frauengesundheit, University of Mainz, Langenbeckstr. 1, Mainz, 55131, Germany
| | - Maximilian Klar
- Department of Obstetrics and Gynecology, University of Freiburg, Hugstetter Straße 55, Freiburg, D-79106, Germany
| | - Josef Hoegel
- Institute of Human Genetics, University of Ulm, Albert-Einstein-Allee 11, Ulm, D-89081, Germany
| | - Walther Vogel
- Institute of Human Genetics, University of Ulm, Albert-Einstein-Allee 11, Ulm, D-89081, Germany
| |
Collapse
|
13
|
An Adaptive Fisher's Combination Method for Joint Analysis of Multiple Phenotypes in Association Studies. Sci Rep 2016; 6:34323. [PMID: 27694844 PMCID: PMC5046106 DOI: 10.1038/srep34323] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2016] [Accepted: 09/12/2016] [Indexed: 12/22/2022] Open
Abstract
Currently, the analyses of most genome-wide association studies (GWAS) have been performed on a single phenotype. There is increasing evidence showing that pleiotropy is a widespread phenomenon in complex diseases. Therefore, using only one single phenotype may lose statistical power to identify the underlying genetic mechanism. There is an increasing need to develop and apply powerful statistical tests to detect association between multiple phenotypes and a genetic variant. In this paper, we develop an Adaptive Fisher’s Combination (AFC) method for joint analysis of multiple phenotypes in association studies. The AFC method combines p-values obtained in standard univariate GWAS by using the optimal number of p-values which is determined by the data. We perform extensive simulations to evaluate the performance of the AFC method and compare the power of our method with the powers of TATES, Tippett’s method, Fisher’s combination test, MANOVA, MultiPhen, and SUMSCORE. Our simulation studies show that the proposed method has correct type I error rates and is either the most powerful test or comparable with the most powerful test. Finally, we illustrate our proposed methodology by analyzing whole-genome genotyping data from a lung function study.
Collapse
|