Integrating Genome-wide Association Studies Data Research Articles

Abstract 11993: Integrative Analysis of Single-Cell Transcriptomics and Polygenic Risk Scores Uncovers Vascular Smooth Muscle Cell Heterogeneity in Human Aortic Aneurisms

Background: Atherosclerosis is accompanied by distinct vascular phenotypes that underly its initiation and progression, though their contribution to its pathogenesis is currently unknown. By integrating polygenetic risk profiles from genome-wide association studies (GWAS) of vascular disease with single-cell transcriptional analyses, it is possible to implicate causal cellular subpopulations and functional states. Hypothesis: We sought to build a more precise understanding of the cellular mechanisms underlying vascular disease in humans by combining single-nuclear RNA-sequencing (snRNA-seq) data from vascular tissue with known polygenic risk loci associated with vascular disease. Methods: We performed snRNA-seq on five flash-frozen human aortas, implementing a novel dissociation method to isolate and sequence 14,228 quality-controlled nuclei from three thoracic aortic aneurysms and two samples without known aortopathy, obtained from patients aged 37-70 years. Cell types were defined using canonical gene markers, and non-linear dimensional reduction via tSNE was employed to identify transcriptional sub-phenotypes. The identified pathological programs were validated using polygenic risk and by integrating GWAS data from other vascular single-nuclear transcriptomic datasets. Results: Our analysis revealed ten distinct clusters representing all major vascular cell populations, with vascular smooth muscle cells (VSMCs) and myofibroblasts emerging as the most relevant cell types in relation to vascular disease based on polygenic risk. Within the VSMC population, two distinct subpopulations exhibited differential expression of genes related to the extracellular matrix, contractility, and proliferation. RNA-trajectory analysis further uncovered a continuum of gene expression linking VSMC1, VSMC2, and myofibroblast transcriptional clusters. Conclusions: This study validates novel methods for identifying pathologic VSMC states associated with vascular disease and implicates phenotypic switching as a pathogenic mechanism. We highlight genes and pathways associated with endothelial cell and VSMC state transitions that are enriched in GWAS of coronary artery disease (CAD).

Integration of epigenetic and genetic profiles identifies multiple sclerosis disease-critical cell types and genes

Genome-wide association studies (GWAS) successfully identified multiple sclerosis (MS) susceptibility variants. Despite this notable progress, understanding the biological context of these associations remains challenging, due in part to the complexity of linking GWAS results to causative genes and cell types. Here, we aimed to address this gap by integrating GWAS data with single-cell and bulk chromatin accessibility data and histone modification profiles from immune and nervous systems. MS-GWAS associations are significantly enriched in regulatory regions of microglia and peripheral immune cell subtypes, especially B cells and monocytes. Cell-specific polygenic risk scores were developed to examine the cumulative impact of the susceptibility genes on MS risk and clinical phenotypes, showing significant associations with risk and brain white matter volume. The findings reveal enrichment of GWAS signals in B cell and monocyte/microglial cell-types, consistent with the known pathology and presumed targets of effective MS therapeutics.

A transcriptome-wide association study of uterine fibroids to identify potential genetic markers and toxic chemicals.

Uterine fibroid is one of the most prevalent benign tumors in women, with high socioeconomic costs. Although genome-wide association studies (GWAS) have identified several loci associated with uterine fibroid risks, they could not successfully interpret the biological effects of genomic variants at the gene expression levels. To prioritize uterine fibroid susceptibility genes that are biologically interpretable, we conducted a transcriptome-wide association study (TWAS) by integrating GWAS data of uterine fibroid and expression quantitative loci data. We identified nine significant TWAS genes including two novel genes, RP11-282O18.3 and KBTBD7, which may be causal genes for uterine fibroid. We conducted functional enrichment network analyses using the TWAS results to investigate the biological pathways in which the overall TWAS genes were involved. The results demonstrated the immune system process to be a key pathway in uterine fibroid pathogenesis. Finally, we carried out chemical-gene interaction analyses using the TWAS results and the comparative toxicogenomics database to determine the potential risk chemicals for uterine fibroid. We identified five toxic chemicals that were significantly associated with uterine fibroid TWAS genes, suggesting that they may be implicated in the pathogenesis of uterine fibroid. In this study, we performed an integrative analysis covering the broad application of bioinformatics approaches. Our study may provide a deeper understanding of uterine fibroid etiologies and informative notifications about potential risk chemicals for uterine fibroid.

Identification of novel proteins associated with movement-related adverse antipsychotic effects by integrating GWAS data and human brain proteomes

Genome-wide association studies (GWAS) have identified some variants for movement-related adverse antipsychotic effects (MAAE), while how these variants confer MAAE remains unclear. We used the probabilistic Mendelian randomization (PMR) method to identify candidate proteins for MAAE by integrating MAAE GWASs and protein quantitative trait loci (pQTL) data. An independent pQTL data from the Banner project and brain-derived eQTL data were used to perform confirmatory PMR. A total of 56 proteins were identified as candidate targets for MAAE after false discovery rates (FDR) correction, such as GRIN2B, ADRA1A, and PED4B. 12 genes were replicated in the confirmatory PMR, and 18 genes had consistent evidence at the transcript level. Furthermore, we investigated the associations between candidate proteins and the motor symptoms of Parkinson's disease (PD). There were 24, 38, and 10 candidate proteins that were significantly associated with PD, PD motor subtypes, and PD motor progression, respectively. Enrichment analysis identified 34 GO terms and 17 pathways that may be involved in MAAE, such as glutamatergic synapse, glutamate receptor complex, and GABAergic synapse. Our study identified multiple candidate genes and pathways that were associated with MAAE, providing new insights into the biological mechanism of MAAE and targets for further mechanistic and therapeutic studies.

Genetic Mechanism Revealed of Age-Related Macular Degeneration Based on Fusion of Statistics and Machine Learning Method.

Age-related macular degeneration (AMD) is the most common cause of irreversible vision loss in the developed world which affects the quality of life for millions of elderly individuals worldwide. Genome-wide association studies (GWAS) have identified genetic variants at 34 loci contributing to AMD. To better understand the disease pathogenesis and identify causal genes for AMD, we applied random walk (RW) and support vector machine (SVM) to identify AMD-related genes based on gene interaction relationship and significance of genes. Our model achieved 0.927 of area under the curve (AUC), and 65 novel genes have been identified as AMD-related genes. To verify our results, a statistics method called summary data-based Mendelian randomization (SMR) has been implemented to integrate GWAS data and transcriptome data to verify AMD susceptibility-related genes. We found 45 genes are related to AMD by SMR. Among these genes, 37 genes overlap with those found by SVM-RW. Finally, we revealed the biological process of genetic mutations leading to changes in gene expression leading to AMD. Our results reveal the genetic pathogenic factors and related mechanisms of AMD.

Identification of Critical Core Genes of Sarcoma Based on Centrality Analysis of Networks Nodes

Genome-wide association studies (GWAS) are powerful tools for identifying pathogenic genes of complex diseases and revealing genetic structure of diseases. However, due to gene-to-gene interactions, only a part of the hereditary factors can be revealed. The meta-analysis based on GWAS can integrate gene expression data at multiple levels and reveal the complex relationship between genes. Therefore, we used meta-analysis to integrate GWAS data of sarcoma to establish complex networks and discuss their significant genes. Firstly, we established gene interaction networks based on the data of different subtypes of sarcoma to analyze the node centralities of genes. Secondly, we calculated the significant score of each gene according to the Staged Significant Gene Network Algorithm (SSGNA). Then, we obtained the critical gene set HYC of sarcoma by ranking the scores, and then combined Gene Ontology enrichment analysis and protein network analysis to further screen it. Finally, the critical core gene set Hcore containing 47 genes was obtained and validated by GEPIA analysis. Our method has certain generalization performance to the study of complex diseases with prior knowledge and it is a useful supplement to genome-wide association studies.

Integrative analyses of major histocompatibility complex loci in the genome-wide association studies of major depressive disorder.

Recent European genome-wide association studies (GWAS) have revealed strong statistical correlations between MDD and numerous zero-to-high linked variants in the genomic region containing major histocompatibility complex (MHC) genes (MHC region), but the underlying biological mechanisms are still unclear. To better understand the roles of this genomic region in the neurobiology of MDD, we applied a convergent functional genomics approach to integrate GWAS data of MDD relevant biological phenotypes, gene-expression analyses results obtained from brain samples, and genetic analyses of independent Chinese MDD samples. We observed that independent MDD risk variants in the MHC region were also significantly associated with the relevant biological phenotypes in the predicted directions, including the emotional and cognitive-related phenotypes. Gene-expression analyses further revealed that mRNA expression levels of several MHC region genes in the human brain were associated with MDD risk SNPs and diagnostic status. For instance, a brain-enriched gene ZNF603P consistently showed lower mRNA levels in the individuals carrying MDD risk alleles and in MDD patients. Remarkably, we further found that independent MDD risk SNPs in the MHC region likely converged to affect the mRNA level(s) of the same gene(s), and Europeans and Han Chinese populations have a substantial shared genetic and molecular basis underlying MDD risk associations in the MHC region. These results highlighted several potential pivotal genes at the MHC region in the pathogenesis of MDD. Their common impacts on multiple psychiatric relevant phenotypes also implicated the neurological processes shared by different psychological processes, such as mood and/or cognition, shedding lights on their potential biological mechanisms.

Integrating Coexpression Networks with GWAS to Prioritize Causal Genes in Maize.

Genome-wide association studies (GWAS) have identified loci linked to hundreds of traits in many different species. Yet, because linkage equilibrium implicates a broad region surrounding each identified locus, the causal genes often remain unknown. This problem is especially pronounced in nonhuman, nonmodel species, where functional annotations are sparse and there is frequently little information available for prioritizing candidate genes. We developed a computational approach, Camoco, that integrates loci identified by GWAS with functional information derived from gene coexpression networks. Using Camoco, we prioritized candidate genes from a large-scale GWAS examining the accumulation of 17 different elements in maize (Zea mays) seeds. Strikingly, we observed a strong dependence in the performance of our approach based on the type of coexpression network used: expression variation across genetically diverse individuals in a relevant tissue context (in our case, roots that are the primary elemental uptake and delivery system) outperformed other alternative networks. Two candidate genes identified by our approach were validated using mutants. Our study demonstrates that coexpression networks provide a powerful basis for prioritizing candidate causal genes from GWAS loci but suggests that the success of such strategies can highly depend on the gene expression data context. Both the software and the lessons on integrating GWAS data with coexpression networks generalize to species beyond maize.

Simultaneous inference of phenotype-associated genes and relevant tissues from GWAS data via Bayesian integration of multiple tissue-specific gene networks.

Although genome-wide association studies (GWAS) have successfully identified thousands of genomic loci associated with hundreds of complex traits in the past decade, the debate about such problems as missing heritability and weak interpretability has been appealing for effective computational methods to facilitate the advanced analysis of the vast volume of existing and anticipated genetic data. Towards this goal, gene-level integrative GWAS analysis with the assumption that genes associated with a phenotype tend to be enriched in biological gene sets or gene networks has recently attracted much attention, due to such advantages as straightforward interpretation, less multiple testing burdens, and robustness across studies. However, existing methods in this category usually exploit non-tissue-specific gene networks and thus lack the ability to utilize informative tissue-specific characteristics. To overcome this limitation, we proposed a Bayesian approach called SIGNET (Simultaneously Inference of GeNEs and Tissues) to integrate GWAS data and multiple tissue-specific gene networks for the simultaneous inference of phenotype-associated genes and relevant tissues. Through extensive simulation studies, we showed the effectiveness of our method in finding both associated genes and relevant tissues for a phenotype. In applications to real GWAS data of 14 complex phenotypes, we demonstrated the power of our method in both deciphering genetic basis and discovering biological insights of a phenotype. With this understanding, we expect to see SIGNET as a valuable tool for integrative GWAS analysis, thereby boosting the prevention, diagnosis, and treatment of human inherited diseases and eventually facilitating precision medicine.

Genome-wide association studies in asthma; perhaps, the end of the beginning

A large number of genetic loci contribute towards an individual's susceptibility to asthma and other complex diseases. Genome-wide association studies (GWASs) have provided us with a wealth of loci associated with asthma susceptibility, asthma endotypes and responsiveness to asthma medications. The reproducibility of these genetic loci across different studies highlights the interplay of general and population-specific risk alleles in asthma. Although GWASs have been successful in identifying disease-associated loci, there is still large potential for such studies to provide further insights into asthma pathogenesis. GWASs over the past year have extended study design well beyond the simple case-control and continuous phenotype association formats, for example, including interactions with environmental factors, integrating GWAS data with epigenetic data and GWASs in animal models, incorporating pathway analyses and utilising emerging sequencing technologies. Moving beyond traditional GWAS formats is likely to significantly enhance our understanding of the genetic basis for asthma. This review discusses where we are after half a decade of asthma GWASs, and focuses on advances over the past year that show where the GWAS field is headed in the future.