Joint association testing of common and rare genetic variants using hierarchical modeling

Identification of grouped rare and common variants via penalized logistic regression

In spite of the success of genome-wide association studies in finding many common variants associated with disease, these variants seem to explain only a small proportion of the estimated heritability. Data collection has turned toward exome and whole genome sequencing, but it is well known that single marker methods frequently used for common variants have low power to detect rare variants associated with disease, even with very large sample sizes. In response, a variety of methods have been developed that attempt to cluster rare variants so that they may gather strength from one another under the premise that there may be multiple causal variants within a gene. Most of these methods group variants by gene or proximity, and test one gene or marker window at a time. We propose a penalized regression method (PeRC) that analyzes all genes at once, allowing grouping of all (rare and common) variants within a gene, along with subgrouping of the rare variants, thus borrowing strength from both rare and common variants within the same gene. The method can incorporate either a burden-based weighting of the rare variants or one in which the weights are data driven. In simulations, our method performs favorably when compared to many previously proposed approaches, including its predecessor, the sparse group lasso [Friedman et al., 2010].

Sequence kernel association tests for the combined effect of rare and common variants

Recent developments in sequencing technologies have made it possible to uncover both rare and common genetic variants. Genome-wide association studies (GWASs) can test for the effect of common variants, whereas sequence-based association studies can evaluate the cumulative effect of both rare and common variants on disease risk. Many groupwise association tests, including burden tests and variance-component tests, have been proposed for this purpose. Although such tests do not exclude common variants from their evaluation, they focus mostly on testing the effect of rare variants by upweighting rare-variant effects and downweighting common-variant effects and can therefore lose substantial power when both rare and common genetic variants in a region influence trait susceptibility. There is increasing evidence that the allelic spectrum of risk variants at a given locus might include novel, rare, low-frequency, and common genetic variants. Here, we introduce several sequence kernel association tests to evaluate the cumulative effect of rare and common variants. The proposed tests are computationally efficient and are applicable to both binary and continuous traits. Furthermore, they can readily combine GWAS and whole-exome-sequencing data on the same individuals, when available, and are also applicable to deep-resequencing data of GWAS loci. We evaluate these tests on data simulated under comprehensive scenarios and show that compared with the most commonly used tests, including the burden and variance-component tests, they can achieve substantial increases in power. We next show applications to sequencing studies for Crohn disease and autism spectrum disorders. The proposed tests have been incorporated into the software package SKAT.

Does accounting for gene-environment interactions help uncover association between rare variants and complex diseases?

OBJECTIVE

To determine whether accounting for gene-environment (G×E) interactions improves the power to detect associations between rare variants and a disease, we have extended three statistical methods and compared their power under various simulated disease models.

METHODS

To test for association of a group of rare variants with a disease, Min-P uses the lowest p value within the group of variants, CAST (Cohort Allelic Sums Test) uses an indicator variable to quantify the rare alleles within the group of variants, and SKAT (Sequence Kernel Association Test) uses a logistic regression based on kernel machine. For each method, we incorporate a term for the G×E interaction and test for association and interaction jointly.

RESULTS

When testing for disease association with a set of rare variants, accounting for G×E interactions can improve power in specific situations (pure interaction or high proportion of causal variants interacting with the environment). However, the power of this approach can decrease, in particular in the presence of main genetic or environmental effects. Among the methods compared, the optimized and weighted SKAT performed best, whether to test for genetic association or to test it jointly with G×E interactions.

CONCLUSION

This approach can be used in specific situations but is not appropriate for a primary analysis.

Incorporating prior biologic information for high-dimensional rare variant association studies

BACKGROUND

Given the increasing scale of rare variant association studies, we introduce a method for high-dimensional studies that integrates multiple sources of data as well as allows for multiple region-specific risk indices.

METHODS

Our method builds upon the previous Bayesian risk index by integrating external biological variant-specific covariates to help guide the selection of associated variants and regions. Our extension also incorporates a second level of uncertainty as to which regions are associated with the outcome of interest.

RESULTS

Using a set of study-based simulations, we show that our approach leads to an increase in power to detect true associations in comparison to several commonly used alternatives. Additionally, the method provides multi-level inference at the pathway, region and variant levels.

CONCLUSION

To demonstrate the flexibility of the method to incorporate various types of information and the applicability to high-dimensional data, we apply our method to a single region within a candidate gene study of second primary breast cancer and to multiple regions within a candidate pathway study of colon cancer.