Project description:BackgroundHundreds of candidate genes have been associated with coronary artery disease (CAD) through genome-wide association studies. However, a systematic way to understand the causal mechanism(s) of these genes, and a means to prioritize them for further study, has been lacking. This represents a major roadblock for developing novel disease- and gene-specific therapies for patients with CAD. Recently, powerful integrative genomics analyses pipelines have emerged to identify and prioritize candidate causal genes by integrating tissue/cell-specific gene expression data with genome-wide association study data sets.MethodsWe aimed to develop a comprehensive integrative genomics analyses pipeline for CAD and to provide a prioritized list of causal CAD genes. To this end, we leveraged several complimentary informatics approaches to integrate summary statistics from CAD genome-wide association studies (from UK Biobank and CARDIoGRAMplusC4D) with transcriptomic and expression quantitative trait loci data from 9 cardiometabolic tissue/cell types in the STARNET study (Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task).ResultsWe identified 162 unique candidate causal CAD genes, which exerted their effect from between one and up to 7 disease-relevant tissues/cell types, including the arterial wall, blood, liver, skeletal muscle, adipose, foam cells, and macrophages. When their causal effect was ranked, the top candidate causal CAD genes were CDKN2B (associated with the 9p21.3 risk locus) and PHACTR1; both exerting their causal effect in the arterial wall. A majority of candidate causal genes were represented in cross-tissue gene regulatory co-expression networks that are involved with CAD, with 22/162 being key drivers in those networks.ConclusionsWe identified and prioritized candidate causal CAD genes, also localizing their tissue(s) of causal effect. These results should serve as a resource and facilitate targeted studies to identify the functional impact of top causal CAD genes.

Project description:Coronary artery disease (CAD) is the leading cause of death worldwide. A recent meta-analysis of genome-wide association studies (GWAS) identified 175 loci associated with CAD. The majority of these loci are in non-coding regions and are predicted to regulate gene expression. We hypothesized that a subset of the CAD GWAS loci was associated with the regulation of transcription in vascular smooth muscle cells (SMCs), which play critical roles in the development of CAD. We measured gene expression in SMCs isolated from the ascending aortas of 151 ethnically diverse heart transplant donors in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single nucleotide polymorphism (SNP) markers across the genome. We identified 4,910 expression and 4,412 splice quantitative trait loci (sQTL) that represent regions of the genome associated with transcript abundance and splicing. 3,660 of the expression quantitative trait loci (eQTL) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 of the eQTLs were SMC- and sex-specific, respectively. To identify the effector transcripts regulated by the CAD GWAS loci, we used four distinct colocalization approaches and identified 84 eQTL and 164 sQTLs that colocalized with CAD loci, suggesting a more significant role for the genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. Two of the CAD loci colocolized with SMC eQTL showed sex-specific (AL160313.1 and TERF2IP) and one of the loci colocalized with SMC-specific eQTLs (ALKBH8). 27% and 37% of the sQTLs were unique to quiescent or proliferative SMCs, respectively. The most significantly associated CAD locus, 9p21, was an sQTL for the long non-coding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs. Collectively, these results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in vascular smooth muscle cells.'

Project description:Coronary artery disease (CAD) is the leading cause of mortality and morbidity driven by both genetic and environmental risk factors. Meta-analyses of genome-wide association studies (GWAS) have identified multiple single nucleotide polymorphisms (SNPs) associated with CAD and myocardial infarction (MI) susceptibility in multi-ethnic populations. The majority of these variants reside in non-coding regulatory regions and are co-inherited with hundreds of candidate regulatory SNPs. Herein, we use integrative genomic, epigenomic, and transcriptomic fine-mapping in human coronary artery smooth muscle cells (HCASMC) and tissues to identify causal regulatory variation and mechanisms responsible for CAD associations. Using these genome-wide maps we prioritize 65 candidate variants and perform allele-specific binding and expression analyses on 7 top candidates. We validate our findings in two independent cohorts of diseased human arterial expression quantitative trait loci (eQTL), which together demonstrate fundamental links between CAD associations and regulatory function in the appropriate disease context. We performed ATAC-seq, ChIP-seq, and RNA-seq on human coronary artery smooth muscle cells grown in SmGM-2 Smooth Muscle Growth Medium-2 including hEGF, insulin, hFGF-B and FBS, but without antibiotics (Lonza, #CC-3182). For ATAC-seq and RNA-seq we performed stimulations with growth factors (TGF-B1, PDGF-BB, PDGF-DD) versus serum-free control. We conducted two biological replicates for each condition using independent donors. For ATAC-seq experiments, sequencing was completed on an Illumina Hiseq 2500, paired-end 50bp reads. For ChIP-seq we performed immunoprecipitations using H3K27ac (Abcam ab4729). We conducted two biological replicates using HCASMC from independent donors, and also did an IgG control for these studies. For RNA-seq we also conducted two replicates using HCASMC from independent donors. For both ChIP-seq and RNA-seq experiments, sequencing was completed on an Illumina HiSeq 2500, paired-end 100bp reads. We also performed ex-vivo ATAC-seq on frozen tissues (isolated media) from normal and atherosclerotic human coronary arteries, using three independent donors for each. Sequencing was also completed on an Illumina HiSeq 2500, paired end 50bp reads.

Dataset Information

Prediction of Causal Candidate Genes in Coronary Artery Disease Loci.

Publications

Prediction of Causal Candidate Genes in Coronary Artery Disease Loci.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets