Investigation of multi-trait associations using pathway-based analysis of GWAS summary statistics

Abstract

BackgroundGenome-wide association studies (GWAS) have been successful in identifying disease-associated genetic variants. Recently, an increasing number of GWAS summary statistics have been made available to the research community, providing extensive repositories for studies of human complex diseases. In particular, cross-trait associations at the genetic level can be beneficial from large-scale GWAS summary statistics by using genetic variants that are associated with multiple traits. However, direct assessment of cross-trait associations using susceptibility loci has been challenging due to the complex genetic architectures in most diseases, calling for advantageous methods that could integrate functional interpretation and imply biological mechanisms.ResultsWe developed an analytical framework for systematic integration of cross-trait associations. It incorporates two different approaches to detect enriched pathways and requires only summary statistics. We demonstrated the framework using 25 traits belonging to four phenotype groups. Our results revealed an average of 54 significantly associated pathways (ranged between 18 and 175) per trait. We further proved that pathway-based analysis provided increased power to estimate cross-trait associations compared to gene-level analysis. Based on Fisher’s Exact Test (FET), we identified a total of 24 (53) pairs of trait-trait association at adjusted pFET < 1 × 10− 3 (pFET < 0.01) among the 25 traits. Our trait-trait association network revealed not only many relationships among the traits within the same group but also novel relationships among traits from different groups, which warrants further investigation in future.ConclusionsOur study revealed that risk variants for 25 different traits aggregated in particular biological pathways and that these pathways were frequently shared among traits. Our results confirmed known mechanisms and also suggested several novel insights into the etiology of multi-traits.Electronic supplementary materialThe online version of this article (10.1186/s12864-018-5373-7) contains supplementary material, which is available to authorized users.

Mapping QTL Associated with Resistance to Avian Oncogenic Marek’s Disease Virus (MDV) Reveals Major Candidate Genes and Variants

“…We would hypothesize that many of our candidate genes are working in conjunction with each other to confer the resistance phenotype, whether that be through their proximity in the genome (e.g., CSTA, LAG3 and C1S on Chr5 and HAVCR1 and TIMD4 on Chr13) or their interaction in biological networks (e.g., SOCS1, TLR4 and TREML2 or RELT, LAG3, HAVCR1 and TIMD4 ). This would support the idea that some traits may be associated with pathway-level interactions as opposed to discrete gene functions [ 58 , 59 ].…”

Characterization of genome-wide association study data reveals spatiotemporal heterogeneity of mental disorders

“…This remained true whether we used all modules or parts of the modules (e.g., the most 25% or the most 50% variable modules across all spatiotemporal sites) for the clustering analysis. This is consistent with previous studies that SCZ and BIP shared common polygenic variations [4,14,[34][35][36]. More interestingly, we also observed ASD and ADHD formed in the same cluster away from the other three disorders in 7 out of 12 points and clustered together in 11 out of 12 (92%) points, indicating ASD and ADHD share more genetic background than the other three adult-onset disorders.…”

A Convergent Study of Genetic Variants Associated With Crohn’s Disease: Evidence From GWAS, Gene Expression, Methylation, eQTL and TWAS

“…DEGs and DMGs were obtained using blood samples. Hence, only Pascal genes from GWAS data were suitable for the determination of tissues (Pei et al, 2019b). We performed TSEA using Pascal genes defined at different threshold ( p < 0.05, p < 0.01, p < 5 × 10 -3 , p < 1 × 10 -3 , p < 5 × 10 -4 , p < 1 × 10 -4 , p < 5 × 10 -5 , p < 1 × 10 -5 , and p < 5 × 10 -6 , Figure 2).…”