Gene-set association tests for next-generation sequencing data

The association between single-nucleotide polymorphisms within type 1 interferon pathway genes and human immunodeficiency virus type 1 viral load in antiretroviral-naïve participants

BACKGROUND

Human genetic contribution to HIV progression remains inadequately explained. The type 1 interferon (IFN) pathway is important for host control of HIV and variation in type 1 IFN genes may contribute to disease progression. This study assessed the impact of variations at the gene and pathway level of type 1 IFN on HIV-1 viral load (VL).

METHODS

Two cohorts of antiretroviral (ART) naïve participants living with HIV (PLWH) with either early (START) or advanced infection (FIRST) were analysed separately. Type 1 IFN genes (n = 17) and receptor subunits (IFNAR1, IFNAR2) were examined for both cumulated type 1 IFN pathway analysis and individual gene analysis. SKAT-O was applied to detect associations between the genotype and HIV-1 study entry viral load (log10 transformed) as a proxy for set point VL; P-values were corrected using Bonferroni (P < 0.0025).

RESULTS

The analyses among those with early infection included 2429 individuals from five continents. The median study entry HIV VL was 14,623 (IQR 3460-45100) copies/mL. Across 673 SNPs within 19 type 1 IFN genes, no significant association with study entry VL was detected. Conversely, examining individual genes in START showed a borderline significant association between IFNW1, and study entry VL (P = 0.0025). This significance remained after separate adjustments for age, CD4+ T-cell count, CD4+/CD8+ T-cell ratio and recent infection. When controlling for population structure using linear mixed effects models (LME), in addition to principal components used in the main model, this was no longer significant (p = 0.0244). In subgroup analyses stratified by geographical region, the association between IFNW1 and study entry VL was only observed among African participants, although, the association was not significant when controlling for population structure using LME. Of the 17 SNPs within the IFNW1 region, only rs79876898 (A > G) was associated with study entry VL (p = 0.0020, beta = 0.32; G associated with higher study entry VL than A) in single SNP association analyses. The findings were not reproduced in FIRST participants.

CONCLUSION

Across 19 type 1 IFN genes, only IFNW1 was associated with HIV-1 study entry VL in a cohort of ART-naïve individuals in early stages of their infection, however, this was no longer significant in sensitivity analyses that controlled for population structures using LME.

Matrisome and Immune Pathways Contribute to Extreme Vascular Outcomes in Williams-Beuren Syndrome

BACKGROUND

Supravalvar aortic stenosis (SVAS) is a characteristic feature of Williams-Beuren syndrome (WBS). Its severity varies: ~20% of people with Williams-Beuren syndrome have SVAS requiring surgical intervention, whereas ~35% have no appreciable SVAS. The remaining individuals have SVAS of intermediate severity. Little is known about genetic modifiers that contribute to this variability.

METHODS AND RESULTS

We performed genome sequencing on 473 individuals with Williams-Beuren syndrome and developed strategies for modifier discovery in this rare disease population. Approaches include extreme phenotyping and nonsynonymous variant prioritization, followed by gene set enrichment and pathway-level association tests. We next used GTEx v8 and proteomic data sets to verify expression of candidate modifiers in relevant tissues. Finally, we evaluated overlap between the genes/pathways identified here and those ascertained through larger aortic disease/trait genome-wide association studies. We show that SVAS severity in Williams-Beuren syndrome is associated with increased frequency of common and rarer variants in matrisome and immune pathways. Two implicated matrisome genes (ACAN and LTBP4) were uniquely expressed in the aorta. Many genes in the identified pathways were previously reported in genome-wide association studies for aneurysm, bicuspid aortic valve, or aortic size.

CONCLUSIONS

Smaller sample sizes in rare disease studies necessitate new approaches to detect modifiers. Our strategies identified variation in matrisome and immune pathways that are associated with SVAS severity. These findings suggest that, like other aortopathies, SVAS may be influenced by the balance of synthesis and degradation of matrisome proteins. Leveraging multiomic data and results from larger aorta-focused genome-wide association studies may accelerate modifier discovery for rare aortopathies like SVAS.

From SNP to pathway-based GWAS meta-analysis: do current meta-analysis approaches resolve power and replication in genetic association studies?

Genome-wide association studies (GWAS) have benefited greatly from enhanced high-throughput technology in recent decades. GWAS meta-analysis has become increasingly popular to highlight the genetic architecture of complex traits, informing about the replicability and variability of effect estimations across human ancestries. A wealth of GWAS meta-analysis methodologies have been developed depending on the input data and the outcome information of interest. We present a survey of current approaches from SNP to pathway-based meta-analysis by acknowledging the range of resources and methodologies in the field, and we provide a comprehensive review of different categories of Genome-Wide Meta-analysis methods employed. These methods highlight different levels at which GWAS meta-analysis may be done, including Single Nucleotide Polymorphisms, Genes and Pathways, for which we describe their framework outline. We also discuss the strengths and pitfalls of each approach and make suggestions regarding each of them.

Whole exome sequencing in patients with Williams-Beuren syndrome followed by disease modeling in mice points to four novel pathways that may modify stenosis risk

Supravalvular aortic stenosis (SVAS) is a narrowing of the aorta caused by elastin (ELN) haploinsufficiency. SVAS severity varies among patients with Williams-Beuren syndrome (WBS), a rare disorder that removes one copy of ELN and 25-27 other genes. Twenty percent of children with WBS require one or more invasive and often risky procedures to correct the defect while 30% have no appreciable stenosis, despite sharing the same basic genetic lesion. There is no known medical therapy. Consequently, identifying genes that modify SVAS offers the potential for novel modifier-based therapeutics. To improve statistical power in our rare-disease cohort (N = 104 exomes), we utilized extreme-phenotype cohorting, functional variant filtration and pathway-based analysis. Gene set enrichment analysis of exome-wide association data identified increased adaptive immune system variant burden among genes associated with SVAS severity. Additional enrichment, using only potentially pathogenic variants known to differ in frequency between the extreme phenotype subsets, identified significant association of SVAS severity with not only immune pathway genes, but also genes involved with the extracellular matrix, G protein-coupled receptor signaling and lipid metabolism using both SKAT-O and RQTest. Complementary studies in Eln+/-; Rag1-/- mice, which lack a functional adaptive immune system, showed improvement in cardiovascular features of ELN insufficiency. Similarly, studies in mixed background Eln+/- mice confirmed that variations in genes that increase elastic fiber deposition also had positive impact on aortic caliber. By using tools to improve statistical power in combination with orthogonal analyses in mice, we detected four main pathways that contribute to SVAS risk.

Meta-Qtest: meta-analysis of quadratic test for rare variants

Background

In genome-wide association studies (GWASs), meta-analysis has been widely used to improve statistical power by combining the results of different studies. Meta-analysis can detect phenotype associated variants that are failed to be detected in single studies. Especially, in biomedical sciences, meta-analysis is often necessary not only for improving statistical power, but also for reducing unavoidable limitation in data collection. As next-generation sequencing (NGS) technology has been developed, meta-analysis of rare variants is proceeding briskly along with meta-analysis of common variants in GWASs. However, meta-analysis on a single variant that is commonly used in common variant association test is improper for rare variants. A sparse signal of rare variant undermines the association signal and its large number causes multiple testing problem. To over-come these problems, we propose a meta-analysis method at the gene-level rather than variant level.

Results

Among many methods that have been developed, we used the unified quadratic tests (Q-tests); Q-test is more powerful than or as powerful as other tests such as Sequence Kernel Association Tests (SKAT). Since there are three different versions of Q-test (QTest1, QTest2, QTest3), each assumes different relationships among multiple rare variants, we extended them into meta-study accordingly. For meta-analysis, we consider two types of approaches, the one is to combine regression coefficients and the other is to combine test statistics from each single study. We extend the Q-test for meta-analysis, proposing Meta Quadratic Test (Meta-Qtest). Meta Q-test avoids the limitations of MetaSKAT. It does not only consider genetic heterogeneity among studies as MetaSKAT but also reflects diverse real situations; since we extend three different Q-tests into meta-analysis respectively, flexible Meta Q-test suggests way to deal with gene-level rare variant meta-analysis efficiently From the results of real data analysis of blood pressure trait, our meta-analysis could successfully discovered genes, KCNA5 and CABIN1 that are already well known for relevance with hypertension disease and they are not detected in MetaSKAT.

Conclusion

As exemplified by an application to T2D Genes projects data set, Meta-Qtest more effectively identified genes associated with hypertension disease than MetaSKAT did.

rqt: an R package for gene-level meta-analysis

Motivation

Despite recent advances of modern GWAS methods, it is still remains an important problem of addressing calculation an effect size and corresponding p-value for the whole gene rather than for single variant.

Results

We developed an R package rqt, which offers gene-level GWAS meta-analysis. The package can be easily included into bioinformatics pipeline or used stand-alone. We applied this tool to the analysis of Alzheimer's disease data from three datasets CHS, FHS and LOADFS. Test results from meta-analysis of three Alzheimer studies show its applicability for association testing.

Availability and implementation

The package rqt is freely available under the following link: https://github.com/izhbannikov/rqt.

Contact

ilya.zhbannikov@duke.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

WISARD: workbench for integrated superfast association studies for related datasets

BACKGROUND

A Mendelian transmission produces phenotypic and genetic relatedness between family members, giving family-based analytical methods an important role in genetic epidemiological studies-from heritability estimations to genetic association analyses. With the advance in genotyping technologies, whole-genome sequence data can be utilized for genetic epidemiological studies, and family-based samples may become more useful for detecting de novo mutations. However, genetic analyses employing family-based samples usually suffer from the complexity of the computational/statistical algorithms, and certain types of family designs, such as incorporating data from extended families, have rarely been used.

RESULTS

We present a Workbench for Integrated Superfast Association studies for Related Data (WISARD) programmed in C/C++. WISARD enables the fast and a comprehensive analysis of SNP-chip and next-generation sequencing data on extended families, with applications from designing genetic studies to summarizing analysis results. In addition, WISARD can automatically be run in a fully multithreaded manner, and the integration of R software for visualization makes it more accessible to non-experts.

CONCLUSIONS

Comparison with existing toolsets showed that WISARD is computationally suitable for integrated analysis of related subjects, and demonstrated that WISARD outperforms existing toolsets. WISARD has also been successfully utilized to analyze the large-scale massive sequencing dataset of chronic obstructive pulmonary disease data (COPD), and we identified multiple genes associated with COPD, which demonstrates its practical value.