Rare Variant Association Tests Research Articles

Rare variants association testing for a binary outcome when pooling individual level data from heterogeneous studies.

Whole genome sequencing (WGS) and whole exome sequencing studies are used to test the association of rare genetic variants with health traits. Many existing WGS efforts now aggregate data from heterogeneous groups, for example, combining sets of individuals of European and African ancestries. We here investigate the statistical implications on rare variant association testing with a binary trait when combining together heterogeneous studies, defined as studies with potentially different disease proportion and different frequency of variant carriers. We study and compare in simulations the Type 1 error control and power of the naïve score test, the saddlepoint approximation to the score test, and the BinomiRare test in a range of settings, focusing on low numbers of variant carriers. We show that Type 1 error control and power patterns depend on both the number of carriers of the rare allele and on disease prevalence in each of the studies. We develop recommendations for association analysis of rare genetic variants. (1) The Score test is preferred when the case proportion in the sample is 50%. (2) Do not down-sample controls to balance case-control ratio, because it reduces power. Rather, use a test that controls the Type 1 error. (3) Conduct stratified analysis in parallel with combined analysis. Aggregated testing may have lower power when the variant effect size differs between strata.

MF-TOWmuT: Testing an optimally weighted combination of common and rare variants with multiple traits using family data.

With rapid advancements of sequencing technologies and accumulations of electronic health records, a large number of genetic variants and multiple correlated human complex traits have become available in many genetic association studies. Thus, it becomes necessary and important to develop new methods that can jointly analyze the association between multiple genetic variants and multiple traits. Compared with methods that only use a single marker or trait, the joint analysis of multiple genetic variants and multiple traits is more powerful since such an analysis can fully incorporate the correlation structure of genetic variants and/or traits and their mutual dependence patterns. However, most of existing methods that simultaneously analyze multiple genetic variants and multiple traits are only applicable to unrelated samples. We develop a new method called MF-TOWmuT to detect association of multiple phenotypes and multiple genetic variants in a genomic region with family samples. MF-TOWmuT is based on an optimally weighted combination of variants. Our method can be applied to both rare and common variants and both qualitative and quantitative traits. Our simulation results show that (1) the type I error of MF-TOWmuT is preserved; (2) MF-TOWmuT outperforms two existing methods such as Multiple Family-based Quasi-Likelihood Score Test and Multivariate Family-based Rare Variant Association Testin terms of power. We also illustrate the usefulness of MF-TOWmuT by analyzing genotypic and phenotipic data from the Genetics of Kidneys in Diabetes study. R program is available at https://github.com/gaochengPRC/MF-TOWmuT.

Increased Burden of Rare Sequence Variants in GnRH-Associated Genes in Women With Hypothalamic Amenorrhea.

Functional hypothalamic amenorrhea (HA) is a common, acquired form of hypogonadotropic hypogonadism that occurs in the setting of energy deficits and/or stress. Variability in individual susceptibility to these stressors, HA heritability, and previous identification of several rare sequence variants (RSVs) in genes associated with the rare disorder, isolated hypogonadotropic hypogonadism (IHH), in individuals with HA suggest a possible genetic contribution to HA susceptibility. We sought to determine whether the burden of RSVs in IHH-related genes is greater in women with HA than controls. We compared patients with HA to control women. The study was conducted at secondary referral centers. Women with HA (n = 106) and control women (ClinSeq study; n = 468). We performed exome sequencing in all patients and controls. The frequency of RSVs in 53 IHH-associated genes was determined using rare variant burden and association tests. RSVs were overrepresented in women with HA compared with controls (P = .007). Seventy-eight heterozygous RSVs in 33 genes were identified in 58 women with HA (36.8% of alleles) compared to 255 RSVs in 41 genes among 200 control women (27.2%). Women with HA are enriched for RSVs in genes that cause IHH, suggesting that variation in genes associated with gonadotropin-releasing hormone neuronal ontogeny and function may be a major determinant of individual susceptibility to developing HA in the face of diet, exercise, and/or stress.

A rank-based normalization method with the fully adjusted full-stage procedure in genetic association studies.

In the area of genetic epidemiology, studies of the genotype-phenotype associations have made significant contributions to human complicated trait genetics. These studies depend on specialized statistical methods for uncover the association between traits and genetic variants, both common and rare variants. Often, in analyzing such studies, potentially confounding factors, such as social and environmental conditions, are required to be involved. Multiple linear regression is the most widely used type of regression analysis when the outcome of interest is quantitative traits. Many statistical tests for identifying genotype-phenotype associations using linear regression rely on the assumption that the traits (or the residuals) of the regression follow a normal distribution. In genomic research, the rank-based inverse normal transformation (INT) is one of the most popular approaches to reach normally distributed traits (or normally distributed residuals). Many researchers believe that applying the INT to the non-normality of the traits (or the non-normality of the residuals) is required for valid inference, because the phenotypic (or residual) outliers and non-normality have the significant influence on both the type I error rate control and statistical power, especially under the situation in rare-variant association testing procedures. Here we propose a test for exploring the association of the rare variant with the quantitative trait by using a fully adjusted full-stage INT. Using simulations we show that the fully adjusted full-stage INT is more appropriate than the existing INT methods, such as the fully adjusted two-stage INT and the INT-based omnibus test, in testing genotype-phenotype associations with rare variants, especially when genotypes are uncorrelated with covariates. The fully adjusted full-stage INT retains the advantages of the fully adjusted two-stage INT and ameliorates the problems of the fully adjusted two-stage INT for analysis of rare variants under non-normality of the trait. We also present theoretical results on these desirable properties. In addition, the two available methods with non-normal traits, the quantile/median regression method and the Yeo-Johnson power transformation, are also included in simulations for comparison with these desirable properties.

Novel Mutations and Decreased Expression of the Epigenetic Regulator TET2 in Pulmonary Arterial Hypertension.

Pulmonary arterial hypertension (PAH) is a lethal vasculopathy. Hereditary cases are associated with germline mutations in BMPR2 and 16 other genes; however, these mutations occur in <25% of patients with idiopathic PAH and are rare in PAH associated with connective tissue diseases. Preclinical studies suggest epigenetic dysregulation, including altered DNA methylation, promotes PAH. Somatic mutations of Tet-methylcytosine-dioxygenase-2 (TET2), a key enzyme in DNA demethylation, occur in cardiovascular disease and are associated with clonal hematopoiesis, inflammation, and adverse vascular remodeling. The role of TET2 in PAH is unknown. To test for a role of TET2, we used a cohort of 2572 cases from the PAH Biobank. Within this cohort, gene-specific rare variant association tests were performed using 1832 unrelated European patients with PAH and 7509 non-Finnish European subjects from the Genome Aggregation Database (gnomAD) as control subjects. In an independent cohort of 140 patients, we quantified TET2 expression in peripheral blood mononuclear cells. To assess causality, we investigated hemodynamic and histological evidence of PAH in hematopoietic Tet2-knockout mice. We observed an increased burden of rare, predicted deleterious germline variants in TET2 in PAH patients of European ancestry (9/1832) compared with control subjects (6/7509; relative risk=6; P=0.00067). Assessing the whole cohort, 0.39% of patients (10/2572) had 12 TET2 mutations (75% predicted germline and 25% somatic). These patients had no mutations in other PAH-related genes. Patients with TET2 mutations were older (71±7 years versus 48±19 years; P<0.0001), were more unresponsive to vasodilator challenge (0/7 versus 140/1055 [13.2%]), had lower pulmonary vascular resistance (5.2±3.1 versus 10.5±7.0 Wood units; P=0.02), and had increased inflammation (including elevation of interleukin-1β). Circulating TET2 expression did not correlate with age and was decreased in >86% of PAH patients. Tet2-knockout mice spontaneously developed PAH, adverse pulmonary vascular remodeling, and inflammation, with elevated levels of cytokines, including interleukin-1β. Long-term therapy with an antibody targeting interleukin-1β blockade resulted in regression of PAH. PAH is the first human disease related to potential TET2 germline mutations. Inherited and acquired abnormalities of TET2 occur in 0.39% of PAH cases. Decreased TET2 expression is ubiquitous and has potential as a PAH biomarker.

A novel association test for rare variants based on algebraic statistics

With the rapid growth of next-generation sequencing technology, more and more rare variants are available in the human genome. In recent years, the point of study has already changed direction to rare variants in genome-wide association studies (GWAS). Although a variety of approaches have been proposed to test associations between rare variants and phenotypes of interest, it is far from the end of this problem, and it is worth exploring new statistical methods based on special features of rare variants. As we all know, the most direct way is to evaluate the association in a two-way contingency table if the phenotype is a discrete variable. The numbers of observations are very close or equal to 0s for most of cells in the contingency table due to the extremely low mutation rates of rare variants. In this paper, we propose a novel association test for rare variants based on a generalization of Fisher’s exact test, and the p-value of this exact test can be computed under the multivariate hypergeometric distribution in the framework of algebraic statistics. Simulation results show that our proposed method outperforms the existing methods, despite there is heterogeneity among causal variants. We also successfully apply our method into the genetic association study of coronary artery disease and hypertension from the Wellcome Trust Case Control Consortium.

Leveraging Family History in Case-Control Analyses of Rare Variation.

Standard methods for case-control association studies of rare variation often treat disease outcome as a dichotomous phenotype. However, both theoretical and experimental studies have demonstrated that subjects with a family history of disease can be enriched for risk variation relative to subjects without such history. Assuming family history information is available, this observation motivates the idea of replacing the standard dichotomous outcome variable used in case-control studies with a more informative ordinal outcome variable that distinguishes controls (0), sporadic cases (1), and cases with a family history (2), with the expectation that we should observe increasing number of risk variants with increasing category of the ordinal variable. To leverage this expectation, we propose a novel rare-variant association test that incorporates family history information based on our previous GAMuT framework for rare-variant association testing of multivariate phenotypes. We use simulated data to show that, when family history information is available, our new method outperforms standard rare-variant association methods, like burden and SKAT tests, that ignore family history. We further illustrate our method using a rare-variant study of cleft lip and palate.

Privacy-Preserving Genome-Wide Association Study for Rare Mutations - A Secure FrameWork for Externalized Statistical Analysis

This paper proposes a new privacy-preserving framework to perform rare variant case-control association tests with information provided by two parties: a Genomic Research Unit (GRU) with sequencing data from individuals affected by a disease D (cases); a Genomic Research Center (GRC) with sequencing data from healthy individuals (controls). To identify genes with rare variants involved in D, GRU needs to compare cases against controls using association tests (genome-wide association study). The main originality of our proposal is twofold. First, it positions GRC as a proxy between GRU and the server. Doing so makes it possible to use classical cryptographic tools to securely conduct association tests with no computation complexity increase, contrarily to actual state of the art proposals which are of very high complexity being based on homomorphic encryption, for instance. In particular, we show how sensitive data confidentiality can be ensured with secret key based cryptographic hashing with no need to modify statistical algorithms. In our protocol the server simply conducts statistical analyses on partially hashed data. Secondly, we introduce a novel privacy constraint: GRU's identity should remain unknown to the server as this knowledge can give it clues about GRU's data (e.g., diseases and genes of interest). We exhibit how Pretty Good Privacy (PGP) can be used to solve this problem. We illustrate our protocol in the case of one rare variant association test, the Weighted-Sum Statistic (WSS) algorithm, carried out on real genetic data. This secure WSS achieves the same accuracy as its nonsecure version with no increase of complexity. Furthermore, we establish that our protocol can be extended to the different rare variant association tests available in the literature.

Identification of Influential Variants in Significant Aggregate Rare Variant Tests

Introduction: Studies that examine the role of rare variants in both simple and complex disease are increasingly common. Though the usual approach of testing rare variants in aggregate sets is more powerful than testing individual variants, it is of interest to identify the variants that are plausible drivers of the association. We present a novel method for prioritization of rare variants after a significant aggregate test by quantifying the influence of the variant on the aggregate test of association. Methods: In addition to providing a measure used to rank variants, we use outlier detection methods to present the computationally efficient Rare Variant Influential Filtering Tool (RIFT) to identify a subset of variants that influence the disease association. We evaluated several outlier detection methods that vary based on the underlying variance measure: interquartile range (Tukey fences), median absolute deviation, and SD. We performed 1,000 simulations for 50 regions of size 3 kb and compared the true and false positive rates. We compared RIFT using the Inner Tukey to 2 existing methods: adaptive combination of p values (ADA) and a Bayesian hierarchical model (BeviMed). Finally, we applied this method to data from our targeted resequencing study in idiopathic pulmonary fibrosis (IPF). Results: All outlier detection methods observed higher sensitivity to detect uncommon variants (0.001 < minor allele frequency, MAF > 0.03) compared to very rare variants (MAF <0.001). For uncommon variants, RIFT had a lower median false positive rate compared to the ADA. ADA and RIFT had significantly higher true positive rates than that observed for BeviMed. When applied to 2 regions found previously associated with IPF including 100 rare variants, we identified 6 polymorphisms with the greatest evidence for influencing the association with IPF. Discussion: In summary, RIFT has a high true positive rate while maintaining a low false positive rate for identifying polymorphisms influencing rare variant association tests. This work provides an approach to obtain greater resolution of the rare variant signals within significant aggregate sets; this information can provide an objective measure to prioritize variants for follow-up experimental studies and insight into the biological pathways involved.

The Role of Rare Coding Variants in Parkinson's Disease GWAS Loci.

Introduction: Genome-wide association studies (GWAS) have identified multiple loci associated with Parkinson's disease (PD) risk. The presence of rare variants within these loci that may account for the increased susceptibility requires further investigation.Methods: Using exome sequencing, we performed a comprehensive rare variant screen of genes located within 56 novel PD loci. We first analyzed exomes from 109 subjects in the discovery cohort (85 diagnosed with PD and 24 healthy controls) and filtered for rare coding variants with minor allele frequency <0.01 and combined annotation-dependent depletion > 15. Further investigation of exome data from a replication cohort of 2,859 European patients with PD (International Parkinson's Disease Genomics Consortium) and 24,146 non-Finnish European controls from gnomAD were used for association testing of specific rare variants found in the discovery cohort.Results: Our genetic screening identified 54 potential disease-relevant variants in 71 genes in 109 subjects. Six out of 54 variants were found in two or more patients and were not observed in healthy controls: DNAH1 p.A3639T, STAB1 p.S1089G, ANK2 p.V3634D, ANK2 p.R3906W, SH3GL2 p.G276V, and NOD2 p.G908R. Replication in the International Parkinson's Disease Genomics Consortium (IPDGC) confirmed the association with PD risk for three out of the six identified variants (STAB1 p.S1089G, SH3GL2 p.G276V, and NOD2 p.G908R) (p < 10−3).Conclusion: Our study suggests that some of the associations identified in PD risk loci can be ascribed to rare variants with likely functional effects that modify PD risk.

Powerful rare variant association testing in a copula-based joint analysis of multiple phenotypes.

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.

No genetic evidence for involvement of alcohol dehydrogenase genes in risk for Parkinson's disease

Multiple genes have been implicated in Parkinson's disease (PD), including causal gene variants and risk variants typically identified using genome-wide association studies. Variants in the alcohol dehydrogenase genes ADH1C and ADH1B are among the genes that have been associated with PD, suggesting that this family of genes may be important in PD. As part of the International Parkinson's Disease Genomics Consortium's efforts to scrutinize previously reported risk factors for PD, we explored genetic variation in the alcohol dehydrogenase genes ADH1A, ADH1B, ADH1C, ADH4, ADH5, ADH6, and ADH7 using imputed genome-wide association study data from 15,097 cases and 17,337 healthy controls. Rare-variant association tests and single-variant score tests did not show any statistically significant association of alcohol dehydrogenase genetic variation with the risk for PD.

Effects of Common and Rare Chromosome 4 GABAergic Gene Variation on Alcohol Use and Antisocial Behavior

Epidemiological estimates suggest that nearly half of individuals diagnosed with alcohol use disorder will be diagnosed with another mental health disorder, with strong associations involving other externalizing disorders. Molecular genetic studies investigating the relation between alcohol use disorder and externalizing behaviors (e.g., antisocial behavior) have focused on a cluster of chromosome 4 γ-aminobutyric acid (GABA) receptor genes (GABRG1-A2-A4-B1) but have generated varying results. The current study examined associations between common and rare variation in this region with alcohol use disorder and antisocial behavior using genetic sequencing data. Specifically, the University of California at San Francisco Family Alcoholism Sample (n = 1,610; 62% female) was used to conduct common and rare variant association tests in the GABRG1-A2-A4-B1 cluster with DSM-5 alcohol use disorder symptom counts, antisocial behavior, and a product term representing their interaction. Gene-based analyses of rare variation resulted in a significant association between rare GABRA2 variation and the interaction term. Single-variant analysis yielded only nominally significant associations. The strongest association for alcohol use disorder (rs3756007) was located in GABRA2, the strongest association for antisocial behavior (rs11941860) was located in GABRG1, and the interaction term yielded top associations in GABRA2 (rs2119183) and the intergenic region between GABRA2 and GABRG1 (rs536599). Common and rare variant associations for the interaction remained similar when covarying for the effects of the other type of variation, suggesting that the significant rare variant signal is independent of common variant contributions. The present study suggests that both rare and common variant associations in GABRA2 confer risk for alcohol use disorder and antisocial behaviors, indicating a potential liability toward externalizing behavior more broadly.

Population genetic simulation study of power in association testing across genetic architectures and study designs.

While it is well established that genetics can be a major contributor to population variation of complex traits, the relative contributions of rare and common variants to phenotypic variation remains a matter of considerable debate. Here, we simulate genetic and phenotypic data across different case/control panel sampling strategies, sequencing methods, and genetic architecture models based on evolutionary forces to determine the statistical performance of rare variant association tests (RVATs) widely in use. We find that the highest statistical power of RVATs is achieved by sampling case/control individuals from the extremes of an underlying quantitative trait distribution. We also demonstrate that the use of genotyping arrays, in conjunction with imputation from a whole‐genome sequenced (WGS) reference panel, recovers the vast majority (90%) of the power that could be achieved by sequencing the case/control panel using current tools. Finally, we show that for dichotomous traits, the statistical performance of RVATs decreases as rare variants become more important in the trait architecture. Our results extend previous work to show that RVATs are insufficiently powered to make generalizable conclusions about the role of rare variants in dichotomous complex traits.

Detecting Clustered Independent Rare Variant Associations Using Genetic Algorithms.

The availability of an increasing collection of sequencing data provides the opportunity to study genetic variation with an unprecedented level of detail. There is much interest in uncovering the role of rare variants and their contribution to disease. However, detecting associations of rare variants with small minor allele frequencies (MAF) and modest effects remains a challenge for rare variant association methods. Due to this low signal-to-noise ratio, most methods are underpowered to detect associations even when conducting rare variant association tests at the gene level. We present a new method for detecting rare variant associations. The algorithm consists of two steps. In the first step, a genetic algorithm searches for a promising genomic region containing a collection of genes with causal rare variants. In the second step, a genetic algorithm aims at removing false positives from the located genomic region. We tested the proposed method with a collection of datasets obtained from real exome data. The proposed method possesses sufficient power for detecting associations of rare variants with complex phenotypes. This method can be used for studying the contribution of rare variants with complex disease, particularly in cases where single-variant or gene-based tests are underpowered.

FunSPU: A versatile and adaptive multiple functional annotation-based association test of whole-genome sequencing data.

Despite ongoing large-scale population-based whole-genome sequencing (WGS) projects such as the NIH NHLBI TOPMed program and the NHGRI Genome Sequencing Program, WGS-based association analysis of complex traits remains a tremendous challenge due to the large number of rare variants, many of which are non-trait-associated neutral variants. External biological knowledge, such as functional annotations based on the ENCODE, Epigenomics Roadmap and GTEx projects, may be helpful in distinguishing causal rare variants from neutral ones; however, each functional annotation can only provide certain aspects of the biological functions. Our knowledge for selecting informative annotations a priori is limited, and incorporating non-informative annotations will introduce noise and lose power. We propose FunSPU, a versatile and adaptive test that incorporates multiple biological annotations and is adaptive at both the annotation and variant levels and thus maintains high power even in the presence of noninformative annotations. In addition to extensive simulations, we illustrate our proposed test using the TWINSUK cohort (n = 1,752) of UK10K WGS data based on six functional annotations: CADD, RegulomeDB, FunSeq, Funseq2, GERP++, and GenoSkyline. We identified genome-wide significant genetic loci on chromosome 19 near gene TOMM40 and APOC4-APOC2 associated with low-density lipoprotein (LDL), which are replicated in the UK10K ALSPAC cohort (n = 1,497). These replicated LDL-associated loci were missed by existing rare variant association tests that either ignore external biological information or rely on a single source of biological knowledge. We have implemented the proposed test in an R package “FunSPU”.

A modified association test for rare and common variants based on affected sib-pair design

Current genome-wide association analysis has identified a great number of rare and common variants associated with common complex traits, however, more effective approaches for detecting associations between rare and common variants with common diseases are still demanded. Approaches for detecting rare variant association analysis will compromise the power when detecting the effects of rare and common variants simultaneously. In this paper, we extend an existing method of testing for rare variant association based on affected sib pairs (TOW-sib) and propose a variable weight test for rare and common variants association based on affected sib pairs (abbreviated as VW-TOWsib). The VW-TOWsib can be used to achieve the purpose of detecting the association of rare and common variants with complex diseases. Simulation results in various scenarios show that our proposed method is more powerful than existing methods for detecting effects of rare and common variants. At the same time, the VW-TOWsib also performs well as a method for rare variant association analysis.

Family-based association tests for rare variants with censored traits.

We propose a set of family-based burden and kernel tests for censored traits (FamBAC and FamKAC). Here, censored traits refer to time-to-event outcomes, for instance, age-at-onset of a disease. To model censored traits in family-based designs, we used the frailty model, which incorporated not only fixed genetic effects of rare variants in a region of interest but also random polygenic effects shared within families. We first partitioned genotype scores of rare variants into orthogonal between- and within-family components, and then derived their corresponding efficient score statistics from the frailty model. Finally, FamBAC and FamKAC were constructed by aggregating the weighted efficient scores of the within-family components across rare variants and subjects. FamBAC collapsed rare variants within subject first to form a burden test that followed a chi-squared distribution; whereas FamKAC was a variant component test following a mixture of chi-squared distributions. For FamKAC, p-values can be computed by permutation tests or for computational efficiency by approximation methods. Through simulation studies, we showed that type I error was correctly controlled by FamBAC for various variant weighting schemes (0.0371 to 0.0527). However, FamKAC type I error rates based on approximation methods were deflated (max 0.0376) but improved by permutation tests. Our simulations also demonstrated that burden test FamBAC had higher power than kernel test FamKAC when high proportion (e.g. ≥ 80%) of causal variants had effects in the same direction. In contrast, when the effects of causal variants on the censored trait were in mixed directions, FamKAC outperformed FamBAC and had comparable or higher power than an existing method, RVFam. Our proposed framework has the flexibility to accommodate general nuclear families, and can be used to analyze sequence data for censored traits such as age-at-onset of a complex disease of interest.

Multivariate association test for rare variants controlling for cryptic and family relatedness

AbstractIn genetic studies of complex diseases, multiple measures of related phenotypes are often collected. Jointly analyzing these phenotypes may improve statistical power to detect sets of rare variants affecting multiple traits. In this work, we consider association testing between a set of rare variants and multiple phenotypes in family‐based designs. We use a mixed linear model to express the correlations among the phenotypes and between related individuals. Given the many sources of correlations in this situation, deriving an appropriate test statistic is not straightforward. We derive a vector of score statistics, whose joint distribution is approximated using a copula. This allows us to have closed‐form expressions for the p‐values of several test statistics. A comprehensive simulation study and an application to Genetic Analysis Workshop 18 data highlight the gains associated with joint testing over univariate approaches, especially in the presence of pleiotropy or highly correlated phenotypes. The Canadian Journal of Statistics 47: 90–107; 2019 © 2018 Statistical Society of Canada

A simple and accurate method to determine genomewide significance for association tests in sequencing studies.

Whole-exome sequencing (WES) and whole-genome sequencing (WGS) studies are underway to investigate the impact of genetic variants on complex diseases and traits. It is customary to perform single-variant association tests for common variants and region-based association tests for rare variants. The latter may target variants with similar or opposite effects, interrogate variants with different frequencies or different functional annotations, and examine a variety of regions. The large number of tests that are performed necessitates adjustment for multiple testing. The conventional Bonferroni correction is overly conservative as the test statistics are correlated. To address this challenge, we propose a simple and accurate method based on parametric bootstrap to assess genomewide significance. We show that the correlations of the test statistics are determined primarily by the genotypes, such that the same significance threshold can be used in different studies that share a common sequencing platform. We demonstrate the usefulness of the proposed method with WES data from the National Heart, Lung, and Blood Institute Exome Sequencing Project and WGS data from the 1000 Genomes Project. We recommend the p value of as the genomewide significance threshold for testing all common and low-frequency variants (MAFs 0.1%) in the human genome.