Abstract:Genetics is estimated to account for 50% of Coronary Artery Disease (CAD) risk. Genome-Wide Association Studies (GWAS) have revealed 100s of genomic loci which are associated with CAD (1). However, moving from GWAS signals to causal genes has proved challenging, slowing progress in developing interventions for CAD and other diseases. The greatest challenge in interpreting GWAS data is that only 10% of variants are in protein-coding regions and 90% are in "non-coding" genomic regions. Many of these non-coding variants are found in regulatory elements such as enhancers. Enhancers are distal regulatory elements that bind transcription factors and regulate genes up to 1 million base pairs away. This is achieved by DNA looping bringing enhancers into physical proximity with target genes. However, it is not possible to predict which gene(s) a given enhancer will regulate based on genomic location alone. We aimed to develop a generalisable workflow to solve CAD GWAS signals accurately and sensitively, by identifying causal enhancer SNPs and linking them mechanistically to output genes in a cell-type specific manner. We combined cell-specific epigenomic mapping with machine learning to prioritise 21,959 SNPs from 861 candidate CAD loci, in endothelial cells (HUVECs) and smooth muscle cells (HCASMCs). Micro-Capture-C (MCC) was used to resolve enhancer-promoter interactions at base-pair resolution (2). MNase footprinting determined differential transcription factor binding at specific enhancer SNPs. In HCASMCs, 20 putative enhancers containing CAD SNPs were found to interact with 32 genes, and in HUVECs 25 enhancers interacted with 48 genes. These genes were enriched for pathways including TGFb signalling, RUNX3 activity, and MAPK6/MAPK4 signalling. To prove causality at specific loci, heterozygous SNPs were identified in the individual primary cell lines. The interaction frequency at these elements was statistically compared between risk and non-risk alleles. For example, at the 9p21.3 locus, an enhancer containing a non-risk T allele was found to interact more strongly with the CDKN2B gene, compared to the risk G allele. Interestingly, this imbalance between allele interaction was potentiated by stimulation with TGFb, a key mediator in atherosclerosis. Furthermore, the transcription factor footprint at the SNP location was altered by the risk allele, indicating the mechanistic cause of the CAD-associated change in CDKN2B expression. We demonstrate a combined computational and experimental approach to identify causal non-coding variants in CAD risk loci, at scale. Among the key outputs of this study is the identification of the risk variant (rs1537373), output gene (CDKN2B), and specific cell type (HCASMC) underlying the 9p21.3 locus association with CAD risk. Using this approach to translate GWAS signals into cell-type specific outputs genes, will lead to a better understanding of the pathogenesis of CAD, and new therapeutic opportunities.

The Integrated Landscape of Biological Candidate Causal Genes in Coronary Artery Disease

Integrative Prioritization of Causal Genes for Coronary Artery Disease

Prediction of Causal Candidate Genes in Coronary Artery Disease Loci.

Identification of novel proteins for coronary artery disease by integrating GWAS data and human plasma proteomes

Integrative genomic analysis identified common regulatory networks underlying the correlation between coronary artery disease and plasma lipid levels

Integration of Summary Data from GWAS and Eqtl Studies Identified Novel Risk Genes for Coronary Artery Disease

Candidate pathway-based genome-wide association studies identify novel associations of genomic variants in the complement system associated with coronary artery disease.

Transcriptome-wide association study of coronary artery disease identifies novel susceptibility genes

Comprehensive identification of coronary artery disease-associated variants regulating vascular smooth muscle cell gene expression

Convergence of coronary artery disease genes onto endothelial cell programs

Novel Common Variants Associated with Body Mass Index and Coronary Artery Disease Detected Using a Pleiotropic Cfdr Method.

Weighted gene co-expression network analysis identifies specific modules and hub genes related to coronary artery disease

Micro-capture-C in vessel wall cell types identifies mechanisms underlying coronary artery disease susceptibility loci

Identification Of Key Genes And Crucial Modules Associated With Coronary Artery Disease By Bioinformatics Analysis

Integrative analysis prioritizes the relevant genes and risk factors for chronic venous disease

Integrative analysis of liver-specific non-coding regulatory SNPs associated with the risk of coronary artery disease

Comprehensive Analysis of Alzheimer's Disease Biologically Candidate Causal Genes Revealed by Function Association Study with GWAS.

Protein interaction networks in the vasculature prioritize genes and pathways underlying coronary artery disease

Genome-wide linkage scan identifies two novel genetic loci for coronary artery disease: in GeneQuest families.

Landscape of the Relationship Between Type 2 Diabetes and Coronary Heart Disease Through an Integrated Gene Network Analysis.

Coronary Artery Disease Associated Specific Modules and Feature Genes Revealed by Integrative Methods of WGCNA, MetaDE and Machine Learning.