GSE193817
Colocalization of GWAS and molecular QTL signals in human vascular smooth muscle cells predict candidate causal genes for coronary artery disease and vascular disorders
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.' Overall design: RNA sequencing data across proliferative and quiescent smooth muscle cells (SMCs) in ~150 multiethnic donors.