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.

Functional noncoding SNPs in human endothelial cells fine-map vascular trait associations

Comparative and Functional Genomic Resource for Mechanistic Studies of Human Blood Pressure–Associated Single Nucleotide Polymorphisms

Abstract 040: Comparative Genomic Analysis and Functional Insights for Blood Pressure-Associated Single Nucleotide Polymorphisms

Abstract 076: the Blood Pressure-associated Non-coding Single Nucleotide Polymorphism Rs28451064 Influences Gene Expression in Human Endothelial Cells

Chromatin Interactions Between Regulatory Elements Influence Gene Expression in Human Ipsc-Derived Endothelial and Vascular Smooth Muscle Cells

Genome-wide Maps of Human Kidney Proximal Tubules and Thick Ascending Limbs of Henle Reveal an Enrichment for Human Blood Pressure Associated Polymorphisms

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

Blood Pressure-Associated Non-Coding Single Nucleotide Polymorphism Rs217727 Modulates Gene Expression in Human Endothelial Cells

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

Abstract 067: Genome-wide Chromatin Conformation Map of Human Induced Pluripotent Stem Cell-derived Vascular Endothelial Cells Reveals New Regulatory Insights

Physiological role and mechanisms of action for a long noncoding haplotype region

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

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

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

DNA-binding factor footprints and enhancer RNAs identify functional non-coding genetic variants

Unveiling blood pressure-associated genes in aortic cells through integrative analysis of GWAS and RNA modification-associated variants

A Genetic Variant Associated with Five Vascular Diseases Is a Distal Regulator of Endothelin-1 Gene Expression

Integrative analysis of loss-of-function variants in clinical and genomic data reveals novel genes associated with cardiovascular traits

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

Variant-to-function dissection of rare non-coding GWAS loci with high impact on blood traits

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