1
|
Boutry S, Helaers R, Lenaerts T, Vikkula M. Rare variant association on unrelated individuals in case-control studies using aggregation tests: existing methods and current limitations. Brief Bioinform 2023; 24:bbad412. [PMID: 37974506 DOI: 10.1093/bib/bbad412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 10/14/2023] [Accepted: 10/28/2023] [Indexed: 11/19/2023] Open
Abstract
Over the past years, progress made in next-generation sequencing technologies and bioinformatics have sparked a surge in association studies. Especially, genome-wide association studies (GWASs) have demonstrated their effectiveness in identifying disease associations with common genetic variants. Yet, rare variants can contribute to additional disease risk or trait heterogeneity. Because GWASs are underpowered for detecting association with such variants, numerous statistical methods have been recently proposed. Aggregation tests collapse multiple rare variants within a genetic region (e.g. gene, gene set, genomic loci) to test for association. An increasing number of studies using such methods successfully identified trait-associated rare variants and led to a better understanding of the underlying disease mechanism. In this review, we compare existing aggregation tests, their statistical features and scope of application, splitting them into the five classical classes: burden, adaptive burden, variance-component, omnibus and other. Finally, we describe some limitations of current aggregation tests, highlighting potential direction for further investigations.
Collapse
Affiliation(s)
- Simon Boutry
- Human Molecular Genetics, de Duve Institute, University of Louvain, Avenue Hippocrate 74 (+5) bte B1.74.06, 1200 Brussels, Belgium
- Interuniversity Institute of Bioinformatics in Brussels, Université Libre de Bruxelles-Vrije Universiteit Brussels, 1050 Brussels, Belgium
| | - Raphaël Helaers
- Human Molecular Genetics, de Duve Institute, University of Louvain, Avenue Hippocrate 74 (+5) bte B1.74.06, 1200 Brussels, Belgium
| | - Tom Lenaerts
- Interuniversity Institute of Bioinformatics in Brussels, Université Libre de Bruxelles-Vrije Universiteit Brussels, 1050 Brussels, Belgium
- Machine Learning Group, Université Libre de Bruxelles, 1050 Brussels, Belgium
- Artificial Intelligence laboratory, Vrije Universiteit Brussel, 1050 Brussels, Belgium
| | - Miikka Vikkula
- Human Molecular Genetics, de Duve Institute, University of Louvain, Avenue Hippocrate 74 (+5) bte B1.74.06, 1200 Brussels, Belgium
- WELBIO department, WEL Research Institute, avenue Pasteur, 6, 1300 Wavre, Belgium
| |
Collapse
|
2
|
A combination test for detection of gene-environment interaction in cohort studies. Genet Epidemiol 2017; 41:396-412. [DOI: 10.1002/gepi.22043] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2016] [Revised: 02/06/2017] [Accepted: 02/06/2017] [Indexed: 12/24/2022]
|
3
|
Callegaro A, Spiessens B, Dizier B, Montoya FU, van Houwelingen HC. Testing interaction between treatment and high-dimensional covariates in randomized clinical trials. Biom J 2016; 59:672-684. [PMID: 27763683 DOI: 10.1002/bimj.201500194] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2015] [Revised: 05/18/2016] [Accepted: 08/16/2016] [Indexed: 11/05/2022]
Abstract
In this paper, we considered different methods to test the interaction between treatment and a potentially large number (p) of covariates in randomized clinical trials. The simplest approach was to fit univariate (marginal) models and to combine the univariate statistics or p-values (e.g., minimum p-value). Another possibility was to reduce the dimension of the covariates using the principal components (PCs) and to test the interaction between treatment and PCs. Finally, we considered the Goeman global test applied to the high-dimensional interaction matrix, adjusted for the main (treatment and covariates) effects. These tests can be used for personalized medicine to test if a large set of biomarkers can be useful to identify a subset of patients who may be more responsive to treatment. We evaluated the performance of these methods on simulated data and we applied them on data from two early phases oncology clinical trials.
Collapse
Affiliation(s)
| | - Bart Spiessens
- GSK Vaccines, Rue de l'Institut 89, 1330, Rixensart, Belgium
| | - Benjamin Dizier
- GSK Vaccines, Rue de l'Institut 89, 1330, Rixensart, Belgium
| | | | - Hans C van Houwelingen
- Department of Medical Statistics, Leiden University Medical Center, P.O. Box 9604, 2300 RC Leiden, The Netherlands
| |
Collapse
|
4
|
Wang Y, Li D, Wei P. Powerful Tukey's One Degree-of-Freedom Test for Detecting Gene-Gene and Gene-Environment Interactions. Cancer Inform 2015; 14:209-18. [PMID: 26064040 PMCID: PMC4459566 DOI: 10.4137/cin.s17305] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2014] [Revised: 04/20/2015] [Accepted: 04/28/2015] [Indexed: 12/17/2022] Open
Abstract
Genome-wide association studies (GWASs) have identified thousands of single nucleotide polymorphisms (SNPs) robustly associated with hundreds of complex human diseases including cancers. However, the large number of GWAS-identified genetic loci only explains a small proportion of the disease heritability. This “missing heritability” problem has been partly attributed to the yet-to-be-identified gene–gene (G × G) and gene–environment (G × E) interactions. In spite of the important roles of G × G and G × E interactions in understanding disease mechanisms and filling in the missing heritability, straightforward GWAS scanning for such interactions has very limited statistical power, leading to few successes. Here we propose a two-step statistical approach to test G × G/G × E interactions: the first step is to perform principal component analysis (PCA) on the multiple SNPs within a gene region, and the second step is to perform Tukey’s one degree-of-freedom (1-df) test on the leading PCs. We derive a score test that is computationally fast and numerically stable for the proposed Tukey’s 1-df interaction test. Using extensive simulations we show that the proposed approach, which combines the two parsimonious models, namely, the PCA and Tukey’s 1-df form of interaction, outperforms other state-of-the-art methods. We also demonstrate the utility and efficiency gains of the proposed method with applications to testing G × G interactions for Crohn’s disease using the Wellcome Trust Case Control Consortium (WTCCC) GWAS data and testing G × E interaction using data from a case–control study of pancreatic cancer.
Collapse
Affiliation(s)
- Yaping Wang
- Department of Biostatistics, School of Public Health, University of Texas Health Science Center
| | - Donghui Li
- Department of Gastrointestinal Medical Oncology, The University of Texas MD Anderson Cancer Center
| | - Peng Wei
- Department of Biostatistics, School of Public Health, University of Texas Health Science Center ; Human Genetics Center, School of Public Health, University of Texas Health Science Center, Houston, TX, USA
| |
Collapse
|
5
|
Fan Y, Kong Y, Li D, Zheng Z. Innovated interaction screening for high-dimensional nonlinear classification. Ann Stat 2015. [DOI: 10.1214/14-aos1308] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
6
|
Kruzliak P, Haley AP, Starcevic JN, Gaspar L, Petrovic D. Polymorphisms of the peroxisome proliferator-activated receptor-γ (rs1801282) and its coactivator-1 (rs8192673) are associated with obesity indexes in subjects with type 2 diabetes mellitus. Cardiovasc Diabetol 2015; 14:42. [PMID: 25928419 PMCID: PMC4450508 DOI: 10.1186/s12933-015-0197-0] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/08/2015] [Accepted: 03/20/2015] [Indexed: 12/19/2022] Open
Abstract
ᅟ The aim of this study was to clarify whether common single nucleotide polymorphisms (SNPs) of the Peroxisome Proliferator-Activated Receptor-γ (PPAR-γ) gene (rs1801282) and the Peroxisome Proliferator-Activated Receptor-γ Coactivator-1 (PGC-1α) gene (rs8192673) are associated with obesity indexes (BMI, waist circumference) in subjects with type 2 diabetes mellitus (T2DM) in Caucasian population. The second aim was to find an association of both polymorphisms with T2DM. Methods Two exonic SNPs of both genes rs1801282 of the PPAR-γ gene and rs8192673 of the PGC-1α gene) were genotyped in 881 unrelated Slovene subjects (Caucasians) with T2DM and in 348 subjects without T2DM (control subjects). Results Female homozygotes with the CC genotype of the rs8192673 had higher waist circumference in comparison with subjects with other genotypes. Homozygotes (females, males) with wild allele (Pro) of the rs1801282 (Pro12Ala polymorphism) had higher waist circumference in comparison with subjects with other genotypes. In the study, there were no differences in the distributions of the rs8192673 and the rs1801282 genotypes between patients with T2DM and controls. Linear regression analyses for both polymorphisms were performed and demonstrated an independent effect of the rs1801282 of the PPAR-γ on waist circumference in subjects with T2DM, whereas an independent effect on waist circumference was not demonstrated for the rs8192673 of the PGC-1α gene. Conclusions In a large sample of the Caucasians the rs8192673 of the PGC-1α gene and the rs1801282 of the PPAR-γ gene were associated with waist circumference in subjects with T2DM.
Collapse
Affiliation(s)
- Peter Kruzliak
- Department of Cardiovascular Diseases, International Clinical Research Center, St Anne's University Hospital and Masaryk University, Brno, Czech Republic.
| | - Andreana P Haley
- Department of Psychology, The University of Texas, Austin, TX, USA. .,University of Texas Imaging Research Center, Austin, TX, USA.
| | | | - Ludovit Gaspar
- 2nd Department of Internal Medicine, University Hospital and Comenius University, Bratislava, Slovak Republic.
| | - Daniel Petrovic
- Institute of Histology and Embryology, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia.
| |
Collapse
|
7
|
Wei P, Tang H, Li D. Functional logistic regression approach to detecting gene by longitudinal environmental exposure interaction in a case-control study. Genet Epidemiol 2014; 38:638-51. [PMID: 25219575 DOI: 10.1002/gepi.21852] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2014] [Revised: 05/29/2014] [Accepted: 07/29/2014] [Indexed: 12/26/2022]
Abstract
Most complex human diseases are likely the consequence of the joint actions of genetic and environmental factors. Identification of gene-environment (G × E) interactions not only contributes to a better understanding of the disease mechanisms, but also improves disease risk prediction and targeted intervention. In contrast to the large number of genetic susceptibility loci discovered by genome-wide association studies, there have been very few successes in identifying G × E interactions, which may be partly due to limited statistical power and inaccurately measured exposures. Although existing statistical methods only consider interactions between genes and static environmental exposures, many environmental/lifestyle factors, such as air pollution and diet, change over time, and cannot be accurately captured at one measurement time point or by simply categorizing into static exposure categories. There is a dearth of statistical methods for detecting gene by time-varying environmental exposure interactions. Here, we propose a powerful functional logistic regression (FLR) approach to model the time-varying effect of longitudinal environmental exposure and its interaction with genetic factors on disease risk. Capitalizing on the powerful functional data analysis framework, our proposed FLR model is capable of accommodating longitudinal exposures measured at irregular time points and contaminated by measurement errors, commonly encountered in observational studies. We use extensive simulations to show that the proposed method can control the Type I error and is more powerful than alternative ad hoc methods. We demonstrate the utility of this new method using data from a case-control study of pancreatic cancer to identify the windows of vulnerability of lifetime body mass index on the risk of pancreatic cancer as well as genes that may modify this association.
Collapse
Affiliation(s)
- Peng Wei
- Division of Biostatistics and Human Genetics Center, The University of Texas School of Public Health, Houston, Texas, United States of America
| | | | | |
Collapse
|
8
|
Abstract
This article focuses on conducting global testing for association between a binary trait and a set of rare variants (RVs), although its application can be much broader to other types of traits, common variants (CVs), and gene set or pathway analysis. We show that many of the existing tests have deteriorating performance in the presence of many nonassociated RVs: their power can dramatically drop as the proportion of nonassociated RVs in the group to be tested increases. We propose a class of so-called sum of powered score (SPU) tests, each of which is based on the score vector from a general regression model and hence can deal with different types of traits and adjust for covariates, e.g., principal components accounting for population stratification. The SPU tests generalize the sum test, a representative burden test based on pooling or collapsing genotypes of RVs, and a sum of squared score (SSU) test that is closely related to several other powerful variance component tests; a previous study (Basu and Pan 2011) has demonstrated good performance of one, but not both, of the Sum and SSU tests in many situations. The SPU tests are versatile in the sense that one of them is often powerful, although its identity varies with the unknown true association parameters. We propose an adaptive SPU (aSPU) test to approximate the most powerful SPU test for a given scenario, consequently maintaining high power and being highly adaptive across various scenarios. We conducted extensive simulations to show superior performance of the aSPU test over several state-of-the-art association tests in the presence of many nonassociated RVs. Finally we applied the SPU and aSPU tests to the GAW17 mini-exome sequence data to compare its practical performance with some existing tests, demonstrating their potential usefulness.
Collapse
|
9
|
Shigemizu D, Abe T, Morizono T, Johnson TA, Boroevich KA, Hirakawa Y, Ninomiya T, Kiyohara Y, Kubo M, Nakamura Y, Maeda S, Tsunoda T. The construction of risk prediction models using GWAS data and its application to a type 2 diabetes prospective cohort. PLoS One 2014; 9:e92549. [PMID: 24651836 PMCID: PMC3961382 DOI: 10.1371/journal.pone.0092549] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2013] [Accepted: 02/24/2014] [Indexed: 02/07/2023] Open
Abstract
Recent genome-wide association studies (GWAS) have identified several novel single nucleotide polymorphisms (SNPs) associated with type 2 diabetes (T2D). Various models using clinical and/or genetic risk factors have been developed for T2D risk prediction. However, analysis considering algorithms for genetic risk factor detection and regression methods for model construction in combination with interactions of risk factors has not been investigated. Here, using genotype data of 7,360 Japanese individuals, we investigated risk prediction models, considering the algorithms, regression methods and interactions. The best model identified was based on a Bayes factor approach and the lasso method. Using nine SNPs and clinical factors, this method achieved an area under a receiver operating characteristic curve (AUC) of 0.8057 on an independent test set. With the addition of a pair of interaction factors, the model was further improved (p-value 0.0011, AUC 0.8085). Application of our model to prospective cohort data showed significantly better outcome in disease-free survival, according to the log-rank trend test comparing Kaplan-Meier survival curves (p--value 2:09 x 10(-11)). While the major contribution was from clinical factors rather than the genetic factors, consideration of genetic risk factors contributed to an observable, though small, increase in predictive ability. This is the first report to apply risk prediction models constructed from GWAS data to a T2D prospective cohort. Our study shows our model to be effective in prospective prediction and has the potential to contribute to practical clinical use in T2D.
Collapse
Affiliation(s)
- Daichi Shigemizu
- Laboratory for Medical Science Mathematics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Testuo Abe
- Laboratory for Medical Science Mathematics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Takashi Morizono
- Laboratory for Medical Science Mathematics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Todd A. Johnson
- Laboratory for Medical Science Mathematics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Keith A. Boroevich
- Laboratory for Medical Science Mathematics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Yoichiro Hirakawa
- Department of Environmental Medicine, Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan
| | - Toshiharu Ninomiya
- Department of Environmental Medicine, Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan
- Department of Medicine and Clinical Science, Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan
| | - Yutaka Kiyohara
- Department of Environmental Medicine, Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan
| | - Michiaki Kubo
- Laboratory for Genotyping Development, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Yusuke Nakamura
- Human Genome Center, Institute of Medical Science, The University of Tokyo, Tokyo, Japan
| | - Shiro Maeda
- Laboratory for Endocrinology, Metabolism and Kidney Diseases, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Tatsuhiko Tsunoda
- Laboratory for Medical Science Mathematics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| |
Collapse
|
10
|
Simino J, Sung YJ, Kume R, Schwander K, Rao DC. Gene-alcohol interactions identify several novel blood pressure loci including a promising locus near SLC16A9. Front Genet 2013; 4:277. [PMID: 24376456 PMCID: PMC3860258 DOI: 10.3389/fgene.2013.00277] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2013] [Accepted: 11/22/2013] [Indexed: 01/11/2023] Open
Abstract
Alcohol consumption is a known risk factor for hypertension, with recent candidate studies implicating gene-alcohol interactions in blood pressure (BP) regulation. We used 6882 (predominantly) Caucasian participants aged 20-80 years from the Framingham SNP Health Association Resource (SHARe) to perform a genome-wide analysis of SNP-alcohol interactions on BP traits. We used a two-step approach in the ABEL suite to examine genetic interactions with three alcohol measures (ounces of alcohol consumed per week, drinks consumed per week, and the number of days drinking alcohol per week) on four BP traits [systolic (SBP), diastolic (DBP), mean arterial (MAP), and pulse (PP) pressure]. In the first step, we fit a linear mixed model of each BP trait onto age, sex, BMI, and antihypertensive medication while accounting for the phenotypic correlation among relatives. In the second step, we conducted 1 degree-of-freedom (df) score tests of the SNP main effect, alcohol main effect, and SNP-alcohol interaction using the maximum likelihood estimates (MLE) of the parameters from the first step. We then calculated the joint 2 df score test of the SNP main effect and SNP-alcohol interaction using MixABEL. The effect of SNP rs10826334 (near SLC16A9) on SBP was significantly modulated by both the number of alcoholic drinks and the ounces of alcohol consumed per week (p-values of 1.27E-08 and 3.92E-08, respectively). Each copy of the G-allele decreased SBP by 3.79 mmHg in those consuming 14 drinks per week vs. a 0.461 mmHg decrease in non-drinkers. Index SNPs in 20 other loci exhibited suggestive (p-value ≤ 1E-06) associations with BP traits by the 1 df interaction test or joint 2 df test, including 3 rare variants, one low-frequency variant, and SNPs near/in genes ESRRG, FAM179A, CRIPT-SOCS5, KAT2B, ADCY2, GLI3, ZNF716, SLIT1, PDE3A, KERA-LUM, RNF219-AS1, CLEC3A, FBXO15, and IGSF5. SNP-alcohol interactions may enhance discovery of novel variants with large effects that can be targeted with lifestyle modifications.
Collapse
Affiliation(s)
- Jeannette Simino
- Division of Biostatistics, Washington University School of MedicineSt. Louis, MO, USA
| | | | | | | | | |
Collapse
|
11
|
Lohmueller KE, Sparsø T, Li Q, Andersson E, Korneliussen T, Albrechtsen A, Banasik K, Grarup N, Hallgrimsdottir I, Kiil K, Kilpeläinen TO, Krarup NT, Pers TH, Sanchez G, Hu Y, Degiorgio M, Jørgensen T, Sandbæk A, Lauritzen T, Brunak S, Kristiansen K, Li Y, Hansen T, Wang J, Nielsen R, Pedersen O. Whole-exome sequencing of 2,000 Danish individuals and the role of rare coding variants in type 2 diabetes. Am J Hum Genet 2013; 93:1072-86. [PMID: 24290377 DOI: 10.1016/j.ajhg.2013.11.005] [Citation(s) in RCA: 116] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2013] [Revised: 10/16/2013] [Accepted: 11/04/2013] [Indexed: 12/15/2022] Open
Abstract
It has been hypothesized that, in aggregate, rare variants in coding regions of genes explain a substantial fraction of the heritability of common diseases. We sequenced the exomes of 1,000 Danish cases with common forms of type 2 diabetes (including body mass index > 27.5 kg/m(2) and hypertension) and 1,000 healthy controls to an average depth of 56×. Our simulations suggest that our study had the statistical power to detect at least one causal gene (a gene containing causal mutations) if the heritability of these common diseases was explained by rare variants in the coding regions of a limited number of genes. We applied a series of gene-based tests to detect such susceptibility genes. However, no gene showed a significant association with disease risk after we corrected for the number of genes analyzed. Thus, we could reject a model for the genetic architecture of type 2 diabetes where rare nonsynonymous variants clustered in a modest number of genes (fewer than 20) are responsible for the majority of disease risk.
Collapse
Affiliation(s)
- Kirk E Lohmueller
- Department of Integrative Biology, University of California, Berkeley, Berkeley, CA 94720, USA
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|