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.

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

Single cell variant to enhancer to gene map for coronary artery disease

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

Abstract 218: Integration Of Coronary Artery Disease GWAS With Bulk And Single-cell Transcriptomics From Atherosclerotic Plaques By Deconvolution, Reveals Novel Smooth Muscle Cell Genes

Genetic Regulation of Atherosclerosis-Relevant Phenotypes in Human Vascular Smooth Muscle Cells

Convergence of coronary artery disease genes onto endothelial cell programs

Molecular Mechanisms of Coronary Disease Revealed Using Quantitative Trait Loci for TCF21 Binding, Chromatin Accessibility, and Chromosomal Looping

The 9p21.3 coronary artery disease risk locus drives vascular smooth muscle cells to an osteochondrogenic state

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

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

High-Resolution Regulatory Maps Connect Vascular Risk Variants to Disease-Related Pathways

Multimodal CRISPR perturbations of GWAS loci associated with coronary artery disease in vascular endothelial cells

CRISPR perturbations at many coronary artery disease loci impair vascular endothelial cell functions

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

Statistical and Functional Studies Identify Epistasis of Cardiovascular Risk Genomic Variants From Genome-Wide Association Studies

A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease.

Abstract 151: Single Cell Transcriptional and Epigenetic Map of Human Coronary Artery Disease in A Multi-ethnic Population

Abstract 2072: The 9p21.3 Coronary Artery Disease Risk Locus Induces Cell State Transitions In Vascular Smooth Muscle Cells

Mapping atherogenesis mechanisms in smooth muscle cells by targeting genes linked to coronary artery disease

Integrative Functional Analysis of Super Enhancer SNPs for Coronary Artery Disease.

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