A Large-Scale Genome-Wide Study of Gene-Sleep Duration Interactions for Blood Pressure in 811,405 Individuals from Diverse Populations

Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to genes involved in neurological, thyroidal, bone metabolism, and hematopoietic pathways. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausibility of distinct influences of both sleep duration extremes in cardiovascular health. With several of our loci reflecting specificity towards population background or sex, our discovery sheds light on the importance of embracing granularity when addressing heterogeneity entangled in gene-environment interactions, and in therapeutic design approaches for blood pressure management.

Association between bone mineral density and coronary artery calcification: an observational and Mendelian randomization study

Observational studies have reported inconsistent associations between bone mineral density (BMD) and coronary artery calcification (CAC). We examined the observational association of BMD with CAC in two large population-based studies and evaluated the evidence for a potential causal relation between BMD and CAC using polygenic risk scores (PRS), 1- and 2-sample Mendelian randomization (MR) approaches. Our study populations comprised 1414 individuals (mean age 69.9 years, 52.0% women) from the Rotterdam Study and 2233 individuals (mean age 56.5 years, 50.9% women) from the Framingham Heart Study with complete information on CAC and BMD measurements at the total body (TB-), lumbar spine (LS-), and femoral neck (FN-). We used linear regression models to evaluate the observational association between BMD and CAC. Subsequently, we compared the mean CAC across PRSBMD quintile groups at different skeletal sites. In addition, we used the 2-stage least squares regression (2SLS) and the inverse variance weighted (IVW) model as primary methods for 1- and 2-sample MR to test evidence for a potentially causal association. We did not observe robust associations between measured BMD levels and CAC. These results were consistent with a uniform random distribution of mean CAC across PRSBMD quintile groups (p-value >0.05). Moreover, neither 1- nor 2-sample MR supported the possible causal association between BMD and CAC. Our results do not support the contention that lower BMD is (causally) associated with an increased CAC risk. These findings suggest that previously reported epidemiological associations of BMD with CAC are likely explained by unmeasured confounders or shared etiology, rather than by causal pathways underlying both osteoporosis and vascular calcification processes.

Multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci

Genome-wide association studies (GWASs) have identified hundreds of risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWASs, and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotypes to identify quantitative trait loci for expression (eQTLs) and splicing (sQTLs) in coronary arteries from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary artery; 19% exhibited cell-type-specific expression. Colocalization revealed subgroups of eGenes unique to CAD and blood pressure GWAS. Fine-mapping highlighted additional eGenes, including TBX20 and IL5. We also identified sQTLs for 1,690 genes, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing to accurately identify disease-relevant isoform expression. Our work provides a patient-derived coronary artery eQTL resource and exemplifies the need for diverse study populations and multifaceted approaches to characterize gene regulation in disease processes.

Polygenic Risk Score Associates With Atherosclerotic Plaque Characteristics at Autopsy

BACKGROUND

Polygenic risk scores (PRSs) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically.

METHODS

From 4327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner for sudden death between 1994 and 2015, 2455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas, and thrombotic CAD.

RESULTS

After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age, 48.8±14.7 years; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared with subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% versus 50.4%±38.7%; adjusted P<0.001) and a higher frequency of calcification (69.6% versus 35.8%; adjusted P=0.004) and thin-cap fibroatheroma (26.7% versus 9.5%; adjusted P=0.007). Even after adjustment for traditional CAD risk factors, subjects within the highest PRS quintile had higher odds of severe atherosclerosis (ie, ≥75% stenosis; adjusted odds ratio, 3.77 [95% CI, 2.10-6.78]; P<0.001) and plaque rupture (adjusted odds ratio, 4.05 [95% CI, 2.26-7.24]; P<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged ≤50 years (adjusted odds ratio, 4.08 [95% CI, 2.01-8.30]; P<0.001). No statistically significant associations were observed with plaque erosion after adjusting for covariates.

CONCLUSIONS

This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects.

Integrative single-cell meta-analysis reveals disease-relevant vascular cell states and markers in human atherosclerosis

Coronary artery disease (CAD) is characterized by atherosclerotic plaque formation in the arterial wall. CAD progression involves complex interactions and phenotypic plasticity among vascular and immune cell lineages. Single-cell RNA-seq (scRNA-seq) studies have highlighted lineage-specific transcriptomic signatures, but human cell phenotypes remain controversial. Here, we perform an integrated meta-analysis of 22 scRNA-seq libraries to generate a comprehensive map of human atherosclerosis with 118,578 cells. Besides characterizing granular cell-type diversity and communication, we leverage this atlas to provide insights into smooth muscle cell (SMC) modulation. We integrate genome-wide association study data and uncover a critical role for modulated SMC phenotypes in CAD, myocardial infarction, and coronary calcification. Finally, we identify fibromyocyte/fibrochondrogenic SMC markers (LTBP1 and CRTAC1) as proxies of atherosclerosis progression and validate these through omics and spatial imaging analyses. Altogether, we create a unified atlas of human atherosclerosis informing cell state-specific mechanistic and translational studies of cardiovascular diseases.

Treatment of calcific arterial disease via enhancement of autophagy using GSK343

Vascular calcification is a hallmark of atherosclerotic disease and serves as a strong predictor and risk factor for cardiovascular events. Growing evidence suggests that autophagy may play a protective role in early atherosclerosis. The precise effects of autophagy on VSMC-mediated calcification remain unknown. In this study, we utilized multi-omic profiling to investigate impaired autophagy at the transcriptional level as a key driver of VSMC calcification. Our findings revealed that impaired autophagy is an essential determinant of VSMC calcification. We observed that an osteogenic environment affects the open chromatin status of VSMCs, compromising the transcriptional activation of autophagy initiation genes. In vivo experiments involve pharmacological and genetic activation of autophagy using mouse models of spontaneous large (Mgp-/-) and small (Abcc6-/-) artery calcification. Taken together, these data advance our mechanistic understanding of vascular calcification and provide important insights for a broad range of cardiovascular diseases involving VSMC phenotype switch.

Gene-educational attainment interactions in a multi-population genome-wide meta-analysis identify novel lipid loci

Introduction: Educational attainment, widely used in epidemiologic studies as a surrogate for socioeconomic status, is a predictor of cardiovascular health outcomes. Methods: A two-stage genome-wide meta-analysis of low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), and triglyceride (TG) levels was performed while accounting for gene-educational attainment interactions in up to 226,315 individuals from five population groups. We considered two educational attainment variables: "Some College" (yes/no, for any education beyond high school) and "Graduated College" (yes/no, for completing a 4-year college degree). Genome-wide significant (p < 5 × 10-8) and suggestive (p < 1 × 10-6) variants were identified in Stage 1 (in up to 108,784 individuals) through genome-wide analysis, and those variants were followed up in Stage 2 studies (in up to 117,531 individuals). Results: In combined analysis of Stages 1 and 2, we identified 18 novel lipid loci (nine for LDL, seven for HDL, and two for TG) by two degree-of-freedom (2 DF) joint tests of main and interaction effects. Four loci showed significant interaction with educational attainment. Two loci were significant only in cross-population analyses. Several loci include genes with known or suggested roles in adipose (FOXP1, MBOAT4, SKP2, STIM1, STX4), brain (BRI3, FILIP1, FOXP1, LINC00290, LMTK2, MBOAT4, MYO6, SENP6, SRGAP3, STIM1, TMEM167A, TMEM30A), and liver (BRI3, FOXP1) biology, highlighting the potential importance of brain-adipose-liver communication in the regulation of lipid metabolism. An investigation of the potential druggability of genes in identified loci resulted in five gene targets shown to interact with drugs approved by the Food and Drug Administration, including genes with roles in adipose and brain tissue. Discussion: Genome-wide interaction analysis of educational attainment identified novel lipid loci not previously detected by analyses limited to main genetic effects.

The accumulation of erythrocytes quantified and visualized by Glycophorin C in carotid atherosclerotic plaque reflects intraplaque hemorrhage and pre-procedural neurological symptoms

The accumulation of erythrocyte membranes within an atherosclerotic plaque may contribute to the deposition of free cholesterol and thereby the enlargement of the necrotic core. Erythrocyte membranes can be visualized and quantified in the plaque by immunostaining for the erythrocyte marker glycophorin C. Hence, we theorized that the accumulation of erythrocytes quantified by glycophorin C could function as a marker for plaque vulnerability, possibly reflecting intraplaque hemorrhage (IPH), and offering predictive value for pre-procedural neurological symptoms. We employed the CellProfiler-integrated slideToolKit workflow to visualize and quantify glycophorin C, defined as the total plaque area that is positive for glycophorin C, in single slides of culprit lesions obtained from the Athero-Express Biobank of 1819 consecutive asymptomatic and symptomatic patients who underwent carotid endarterectomy. Our assessment included the evaluation of various parameters such as lipid core, calcifications, collagen content, SMC content, and macrophage burden. These parameters were evaluated using a semi-quantitative scoring method, and the resulting data was dichotomized as predefined criteria into categories of no/minor or moderate/heavy staining. In addition, the presence or absence of IPH was also scored. The prevalence of IPH and pre-procedural neurological symptoms were 62.4% and 87.1%, respectively. The amount of glycophorin staining was significantly higher in samples from men compared to samples of women (median 7.15 (IQR:3.37, 13.41) versus median 4.06 (IQR:1.98, 8.32), p < 0.001). Glycophorin C was associated with IPH adjusted for clinical confounders (OR 1.90; 95% CI 1.63, 2.21; p =  < 0.001). Glycophorin C was significantly associated with ipsilateral pre-procedural neurological symptoms (OR:1.27, 95%CI:1.06-1.41, p = 0.005). Sex-stratified analysis, showed that this was also the case for men (OR 1.37; 95%CI 1.12, 1.69; p = 0.003), but not for women (OR 1.15; 95%CI 0.77, 1.73; p = 0.27). Glycophorin C was associated with classical features of a vulnerable plaque, such as a larger lipid core, a higher macrophage burden, less calcifications, a lower collagen and SMC content. There were marked sex differences, in men, glycophorin C was associated with calcifications and collagen while these associations were not found in women. To conclude, the accumulation of erythrocytes in atherosclerotic plaque quantified and visualized by glycophorin C was independently associated with the presence of IPH, preprocedural symptoms in men, and with a more vulnerable plaque composition in both men and women. These results strengthen the notion that the accumulation of erythrocytes quantified by glycophorin C can be used as a marker for plaque vulnerability.

Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

Lowering of Circulating Sclerostin May Increase Risk of Atherosclerosis and Its Risk Factors: Evidence From a Genome-Wide Association Meta-Analysis Followed by Mendelian Randomization

OBJECTIVE

In this study, we aimed to establish the causal effects of lowering sclerostin, target of the antiosteoporosis drug romosozumab, on atherosclerosis and its risk factors.

METHODS

A genome-wide association study meta-analysis was performed of circulating sclerostin levels in 33,961 European individuals. Mendelian randomization (MR) was used to predict the causal effects of sclerostin lowering on 15 atherosclerosis-related diseases and risk factors.

RESULTS

We found that 18 conditionally independent variants were associated with circulating sclerostin. Of these, 1 cis signal in SOST and 3 trans signals in B4GALNT3, RIN3, and SERPINA1 regions showed directionally opposite signals for sclerostin levels and estimated bone mineral density. Variants with these 4 regions were selected as genetic instruments. MR using 5 correlated cis-SNPs suggested that lower sclerostin increased the risk of type 2 diabetes mellitus (DM) (odds ratio [OR] 1.32 [95% confidence interval (95% CI) 1.03-1.69]) and myocardial infarction (MI) (OR 1.35 [95% CI 1.01-1.79]); sclerostin lowering was also suggested to increase the extent of coronary artery calcification (CAC) (β = 0.24 [95% CI 0.02-0.45]). MR using both cis and trans instruments suggested that lower sclerostin increased hypertension risk (OR 1.09 [95% CI 1.04-1.15]), but otherwise had attenuated effects.

CONCLUSION

This study provides genetic evidence to suggest that lower levels of sclerostin may increase the risk of hypertension, type 2 DM, MI, and the extent of CAC. Taken together, these findings underscore the requirement for strategies to mitigate potential adverse effects of romosozumab treatment on atherosclerosis and its related risk factors.

Leukocytes carrying Clonal Hematopoiesis of Indeterminate Potential (CHIP) Mutations invade Human Atherosclerotic Plaques

Background

Leukocyte progenitors derived from clonal hematopoiesis of undetermined potential (CHIP) are associated with increased cardiovascular events. However, the prevalence and functional relevance of CHIP in coronary artery disease (CAD) are unclear, and cells affected by CHIP have not been detected in human atherosclerotic plaques.

Methods

CHIP mutations in blood and tissues were identified by targeted deep-DNA-sequencing (DNAseq: coverage >3,000) and whole-genome-sequencing (WGS: coverage >35). CHIP-mutated leukocytes were visualized in human atherosclerotic plaques by mutaFISH™. Functional relevance of CHIP mutations was studied by RNAseq.

Results

DNAseq of whole blood from 540 deceased CAD patients of the Munich cardIovaScular StudIes biObaNk (MISSION) identified 253 (46.9%) CHIP mutation carriers (mean age 78.3 years). DNAseq on myocardium, atherosclerotic coronary and carotid arteries detected identical CHIP mutations in 18 out of 25 mutation carriers in tissue DNA. MutaFISH™ visualized individual macrophages carrying DNMT3A CHIP mutations in human atherosclerotic plaques. Studying monocyte-derived macrophages from Stockholm-Tartu Atherosclerosis Reverse Networks Engineering Task (STARNET; n=941) by WGS revealed CHIP mutations in 14.2% (mean age 67.1 years). RNAseq of these macrophages revealed that expression patterns in CHIP mutation carriers differed substantially from those of non-carriers. Moreover, patterns were different depending on the underlying mutations, e.g. those carrying TET2 mutations predominantly displayed upregulated inflammatory signaling whereas ASXL1 mutations showed stronger effects on metabolic pathways.

Conclusions

Deep-DNA-sequencing reveals a high prevalence of CHIP mutations in whole blood of CAD patients. CHIP-affected leukocytes invade plaques in human coronary arteries. RNAseq data obtained from macrophages of CHIP-affected patients suggest that pro-atherosclerotic signaling differs depending on the underlying mutations. Further studies are necessary to understand whether specific pathways affected by CHIP mutations may be targeted for personalized treatment.

Polygenic Risk Score Associates with Atherosclerotic Plaque Characteristics at Autopsy

Background

Polygenic risk scores (PRS) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically.

Methods

From 4,327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner (OCME) for sudden death between 1994 and 2015, 2,455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas (TCFA), and thrombotic CAD.

Results

After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age 48.8±14.7; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared to subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% vs. 50.4%±38.7%; adjusted p<0.001) and a higher frequency of calcification (69.6% vs. 35.8%; adjusted p=0.004) and TCFAs (26.7% vs. 9.5%; adjusted p=0.007). Even after adjustment for traditional CAD risk factors subjects within the highest PRS quintile had higher odds of severe atherosclerosis (i.e., ≥75% stenosis; adjusted OR 3.77; 95%CI 2.10-6.78; p<0.001) and plaque rupture (adjusted OR 4.05; 95%CI 2.26-7.24; p<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged 50 years and younger (adjusted OR 4.08; 95%CI 2.01-8.30; p<0.001). No associations were observed with plaque erosion.

Conclusions

This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects.

Highlights

In this autopsy study including 954 subjects within the CVPath Sudden Death Registry, high PRS correlated with plaque burden and atherosclerosis severity.The PRS showed differential associations with plaque rupture and plaque erosion, suggesting different etiologies to these two causes of thrombotic CAD.PRS may be useful for risk stratification, particularly in the young. Further examination of individual risk loci and their association with plaque morphology may help understand molecular mechanisms of atherosclerosis, potentially revealing new therapy targets of CAD.

Graphic Abstract

A polygenic risk score, generated from 291 known CAD risk loci, was assessed in 954 subjects within the CVPath Sudden Death Registry. Histopathologic examination of the coronary arteries was performed in all subjects. Subjects in the highest PRS quintile exhibited more severe atherosclerosis as compared to subjects in the lowest quintile, with a greater plaque burden, more calcification, and a higher frequency of plaque rupture.

INKILN is a Novel Long Noncoding RNA Promoting Vascular Smooth Muscle Inflammation via Scaffolding MKL1 and USP10

BACKGROUND

Activation of vascular smooth muscle cell (VSMC) inflammation is vital to initiate vascular disease. The role of human-specific long noncoding RNAs in VSMC inflammation is poorly understood.

METHODS

Bulk RNA sequencing in differentiated human VSMCs revealed a novel human-specific long noncoding RNA called inflammatory MKL1 (megakaryoblastic leukemia 1) interacting long noncoding RNA (INKILN). INKILN expression was assessed in multiple in vitro and ex vivo models of VSMC phenotypic modulation as well as human atherosclerosis and abdominal aortic aneurysm. The transcriptional regulation of INKILN was verified through luciferase reporter and chromatin immunoprecipitation assays. Loss-of-function and gain-of-function studies and multiple RNA-protein and protein-protein interaction assays were used to uncover a mechanistic role of INKILN in the VSMC proinflammatory gene program. Bacterial artificial chromosome transgenic mice were used to study INKILN expression and function in ligation injury-induced neointimal formation.

RESULTS

INKILN expression is downregulated in contractile VSMCs and induced in human atherosclerosis and abdominal aortic aneurysm. INKILN is transcriptionally activated by the p65 pathway, partially through a predicted NF-κB (nuclear factor kappa B) site within its proximal promoter. INKILN activates proinflammatory gene expression in cultured human VSMCs and ex vivo cultured vessels. INKILN physically interacts with and stabilizes MKL1, a key activator of VSMC inflammation through the p65/NF-κB pathway. INKILN depletion blocks interleukin-1β-induced nuclear localization of both p65 and MKL1. Knockdown of INKILN abolishes the physical interaction between p65 and MKL1 and the luciferase activity of an NF-κB reporter. Furthermore, INKILN knockdown enhances MKL1 ubiquitination through reduced physical interaction with the deubiquitinating enzyme USP10 (ubiquitin-specific peptidase 10). INKILN is induced in injured carotid arteries and exacerbates ligation injury-induced neointimal formation in bacterial artificial chromosome transgenic mice.

CONCLUSIONS

These findings elucidate an important pathway of VSMC inflammation involving an INKILN/MKL1/USP10 regulatory axis. Human bacterial artificial chromosome transgenic mice offer a novel and physiologically relevant approach for investigating human-specific long noncoding RNAs under vascular disease conditions.

Cardiovascular utility of single cell RNA-Seq

PURPOSE OF REVIEW

Cardiovascular diseases remain the leading causes of morbidity and mortality globally. Single-cell RNA sequencing has the potential to improve diagnostics, risk stratification, and provide novel therapeutic targets that have the potential to improve patient outcomes.

RECENT FINDINGS

Here, we provide an overview of the basic processes underlying single-cell RNA sequencing, including library preparation, data processing, and downstream analyses. We briefly discuss how the technique has been adapted to related medical disciplines, including hematology and oncology, with short term translational impact. We discuss potential applications of this technology within cardiology as well as recent innovative research within the field. We also discuss future directions to translate this technology to other high impact clinical areas.

SUMMARY

The use of single-cell RNA sequencing technology has made significant advancements in the field of cardiology, with ongoing growth in terms of applications and uptake. Most of the current research has focused on structural or atherosclerotic heart disease. Future areas that stand to benefit from this technology include cardiac electrophysiology and cardio-oncology.

FHL5 Controls Vascular Disease-Associated Gene Programs in Smooth Muscle Cells

BACKGROUND

Genome-wide association studies have identified hundreds of loci associated with common vascular diseases, such as coronary artery disease, myocardial infarction, and hypertension. However, the lack of mechanistic insights for many GWAS loci limits their translation into the clinic. Among these loci with unknown functions is UFL1-four-and-a-half LIM (LIN-11, Isl-1, MEC-3) domain 5 (FHL5; chr6q16.1), which reached genome-wide significance in a recent coronary artery disease/ myocardial infarction GWAS meta-analysis. UFL1-FHL5 is also associated with several vascular diseases, consistent with the widespread pleiotropy observed for GWAS loci.

METHODS

We apply a multimodal approach leveraging statistical fine-mapping, epigenomic profiling, and ex vivo analysis of human coronary artery tissues to implicate FHL5 as the top candidate causal gene. We unravel the molecular mechanisms of the cross-phenotype genetic associations through in vitro functional analyses and epigenomic profiling experiments in coronary artery smooth muscle cells.

RESULTS

We prioritized FHL5 as the top candidate causal gene at the UFL1-FHL5 locus through expression quantitative trait locus colocalization methods. FHL5 gene expression was enriched in the smooth muscle cells and pericyte population in human artery tissues with coexpression network analyses supporting a functional role in regulating smooth muscle cell contraction. Unexpectedly, under procalcifying conditions, FHL5 overexpression promoted vascular calcification and dysregulated processes related to extracellular matrix organization and calcium handling. Lastly, by mapping FHL5 binding sites and inferring FHL5 target gene function using artery tissue gene regulatory network analyses, we highlight regulatory interactions between FHL5 and downstream coronary artery disease/myocardial infarction loci, such as FOXL1 and FN1 that have roles in vascular remodeling.

CONCLUSIONS

Taken together, these studies provide mechanistic insights into the pleiotropic genetic associations of UFL1-FHL5. We show that FHL5 mediates vascular disease risk through transcriptional regulation of downstream vascular remodeling gene programs. These transacting mechanisms may explain a portion of the heritable risk for complex vascular diseases.

Gene-repressing epigenetic reader EED unexpectedly enhances cyclinD1 gene activation

Epigenetically switched, proliferative vascular smooth muscle cells (SMCs) form neointima, engendering stenotic diseases. Histone-3 lysine-27 trimethylation (H3K27me3) and acetylation (H3K27ac) marks are associated with gene repression and activation, respectively. The polycomb protein embryonic ectoderm development (EED) reads H3K27me3 and also enhances its deposition, hence is a canonical gene repressor. However, herein we found an unexpected role for EED in activating the bona fide pro-proliferative gene Ccnd1 (cyclinD1). EED overexpression in SMCs increased Ccnd1 mRNA, seemingly contradicting its gene-repressing function. However, consistently, EED co-immunoprecipitated with gene-activating H3K27ac reader BRD4, and they co-occupied at both mitogen-activated Ccnd1 and mitogen-repressed P57 (bona fide anti-proliferative gene), as indicated by chromatin immunoprecipitation qPCR. These results were abolished by an inhibitor of either the EED/H3K27me3 or BRD4/H3K27ac reader function. In accordance, elevating BRD4 increased H3K27me3. In vivo, while EED was upregulated in rat and human neointimal lesions, selective EED inhibition abated angioplasty-induced neointima and reduced cyclinD1 in rat carotid arteries. Thus, results uncover a previously unknown role for EED in Ccnd1 activation, likely via its cooperativity with BRD4 that enhances each other's reader function; i.e., activating pro-proliferative Ccnd1 while repressing anti-proliferative P57. As such, this study confers mechanistic implications for the epigenetic intervention of neointimal pathology.

CD163+ macrophages restrain vascular calcification, promoting the development of high-risk plaque

Vascular calcification (VC) is concomitant with atherosclerosis, yet it remains uncertain why rupture-prone high-risk plaques do not typically show extensive calcification. Intraplaque hemorrhage (IPH) deposits erythrocyte-derived cholesterol, enlarging the necrotic core and promoting high-risk plaque development. Pro-atherogenic CD163+ alternative macrophages engulf hemoglobin:haptoglobin (HH) complexes at IPH sites. However, their role in VC has never been examined to our knowledge. Here we show, in human arteries, the distribution of CD163+ macrophages correlated inversely with VC. In vitro experiments using vascular smooth muscle cells (VSMCs) cultured with HH-exposed human macrophage - M(Hb) - supernatant reduced calcification, while arteries from ApoE-/- CD163-/- mice showed greater VC. M(Hb) supernatant-exposed VSMCs showed activated NF-κB, while blocking NF-κB attenuated the anticalcific effect of M(Hb) on VSMCs. CD163+ macrophages altered VC through NF-κB-induced transcription of hyaluronan synthase (HAS), an enzyme that catalyzes the formation of the extracellular matrix glycosaminoglycan, hyaluronan, within VSMCs. M(Hb) supernatants enhanced HAS production in VSMCs, while knocking down HAS attenuated its anticalcific effect. NF-κB blockade in ApoE-/- mice reduced hyaluronan and increased VC. In human arteries, hyaluronan and HAS were increased in areas of CD163+ macrophage presence. Our findings highlight an important mechanism by which CD163+ macrophages inhibit VC through NF-κB-induced HAS augmentation and thus promote the high-risk plaque development.

Integrative multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci

Genome-wide association studies (GWAS) have identified hundreds of genetic risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWAS and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotype information to identify quantitative trait loci (QTL) for gene expression and splicing in coronary arteries obtained from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary arteries and 19% exhibited cell-type-specific expression. Colocalization analysis with GWAS identified subgroups of eGenes unique to CAD and blood pressure. Fine-mapping highlighted additional eGenes of interest, including TBX20 and IL5 . Splicing (s)QTLs for 1,690 genes were also identified, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing events to accurately identify disease-relevant gene expression. Our work provides the first human coronary artery eQTL resource from a patient sample and exemplifies the necessity of diverse study populations and multi-omic approaches to characterize gene regulation in critical disease processes.

Study Design Overview

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes

BACKGROUND

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes.

METHODS

We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome.

RESULTS

We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs.

CONCLUSIONS

Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

INKILN is a novel long noncoding RNA promoting vascular smooth muscle inflammation via scaffolding MKL1 and USP10

Background

Activation of vascular smooth muscle cells (VSMCs) inflammation is vital to initiate vascular disease. However, the role of human-specific long noncoding RNAs (lncRNAs) in VSMC inflammation is poorly understood.

Methods

Bulk RNA-seq in differentiated human VSMCs revealed a novel human-specific lncRNA called IN flammatory M K L1 I nteracting L ong N oncoding RNA ( INKILN ). INKILN expression was assessed in multiple in vitro and ex vivo models of VSMC phenotypic modulation and human atherosclerosis and abdominal aortic aneurysm (AAA) samples. The transcriptional regulation of INKILN was determined through luciferase reporter system and chromatin immunoprecipitation assay. Both loss- and gain-of-function approaches and multiple RNA-protein and protein-protein interaction assays were utilized to uncover the role of INKILN in VSMC proinflammatory gene program and underlying mechanisms. Bacterial Artificial Chromosome (BAC) transgenic (Tg) mice were utilized to study INKLIN expression and function in ligation injury-induced neointimal formation.

Results

INKILN expression is downregulated in contractile VSMCs and induced by human atherosclerosis and abdominal aortic aneurysm. INKILN is transcriptionally activated by the p65 pathway, partially through a predicted NF-κB site within its proximal promoter. INKILN activates the proinflammatory gene expression in cultured human VSMCs and ex vivo cultured vessels. Mechanistically, INKILN physically interacts with and stabilizes MKL1, a key activator of VSMC inflammation through the p65/NF-κB pathway. INKILN depletion blocks ILIβ-induced nuclear localization of both p65 and MKL1. Knockdown of INKILN abolishes the physical interaction between p65 and MKL1, and the luciferase activity of an NF-κB reporter. Further, INKILN knockdown enhances MKL1 ubiquitination, likely through the reduced physical interaction with the deubiquitinating enzyme, USP10. INKILN is induced in injured carotid arteries and exacerbates ligation injury-induced neointimal formation in BAC Tg mice.

Conclusions

These findings elucidate an important pathway of VSMC inflammation involving an INKILN /MKL1/USP10 regulatory axis. Human BAC Tg mice offer a novel and physiologically relevant approach for investigating human-specific lncRNAs under vascular disease conditions.

Transcriptomic-based clustering of human atherosclerotic plaques identifies subgroups with different underlying biology and clinical presentation

Histopathological studies have revealed key processes of atherosclerotic plaque thrombosis. However, the diversity and complexity of lesion types highlight the need for improved sub-phenotyping. Here we analyze the gene expression profiles of 654 advanced human carotid plaques. The unsupervised, transcriptome-driven clustering revealed five dominant plaque types. These plaque phenotypes were associated with clinical presentation and showed differences in cellular compositions. Validation in coronary segments showed that the molecular signature of these plaques was linked to coronary ischemia. One of the plaque types with the most severe clinical symptoms pointed to both inflammatory and fibrotic cell lineages. Further, we did a preliminary analysis of potential circulating biomarkers that mark the different plaques phenotypes. In conclusion, the definition of the plaque at risk for a thrombotic event can be fine-tuned by in-depth transcriptomic-based phenotyping. These differential plaque phenotypes prove clinically relevant for both carotid and coronary artery plaques and point to distinct underlying biology of symptomatic lesions.

Polygenic risk score associates with atherosclerosis severity at autopsy

Background

Polygenic risk scores (PRS) for coronary artery disease (CAD) are emerging as a potential method to improve cardiovascular risk prediction. Questions remain about their applicability to diverse populations as well as their correlation with coronary histopathology.

Purpose

To assess whether high genetic risk associates with histopathologic coronary plaque morphology.

Methods

We assessed 122 known CAD risk loci in 954 Black and White subjects within our sudden death registry to generate a PRS. The cohort was divided into quintiles according to z-score-standardized PRS, both in a race-stratified fashion and in the pooled sample. Detailed histopathologic examination of the coronary arteries was performed in all subjects.

Results

Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared to subjects in the lowest quintile, with greater mean cross-sectional luminal narrowing (71.5% (95% CI, 66.6%-76.5%) vs. 56.6% (95% CI, 51.1%-62.1%); adjusted p&lt;0.001; Figure 1) and a higher frequency of calcification (adjusted OR 2.19; 95% CI 1.31–3.68; p=0.003) after adjustment for the first 10 principal components, age, sex, and race. Higher z-score-standardized PRS was predictive for the finding of severe atherosclerosis (i.e., ≥75% cross-sectional luminal narrowing) even after additional controlling for traditional CAD risk factors including hypertension, smoking, diabetes mellitus, and hyperlipidemia (adjusted OR 1.39; 95% CI 1.19–1.63; p&lt;0.001; Figure 2). Among Black subjects, higher PRS was associated with higher odds of plaque rupture (adjusted OR 1.31; 95% CI 1.03–1.66; p=0.03) and predicted CAD-associated cause of death among subjects younger than 50 years old (adjusted OR 1.26; 95% CI 1.01–1.58; p=0.04).

Conclusions

This is the first autopsy study investigating associations between PRS and atherosclerotic plaque morphology at a histopathologic level. Our pathological analysis suggests PRS correlates with plaque burden and coronary artery calcification and may be useful as a method for CAD risk stratification, especially in younger subjects.

Funding Acknowledgement

Type of funding sources: Foundation. Main funding source(s): R01 HL141425 Leducq Foundation Grant

PlaqView 2.0: A comprehensive web portal for cardiovascular single-cell genomics

Single-cell RNA-seq (scRNA-seq) is a powerful genomics technology to interrogate the cellular composition and behaviors of complex systems. While the number of scRNA-seq datasets and available computational analysis tools have grown exponentially, there are limited systematic data sharing strategies to allow rapid exploration and re-analysis of single-cell datasets, particularly in the cardiovascular field. We previously introduced PlaqView, an open-source web portal for the exploration and analysis of published atherosclerosis single-cell datasets. Now, we introduce PlaqView 2.0 (www.plaqview.com), which provides expanded features and functionalities as well as additional cardiovascular single-cell datasets. We showcase improved PlaqView functionality, backend data processing, user-interface, and capacity. PlaqView brings new or improved tools to explore scRNA-seq data, including gene query, metadata browser, cell identity prediction, ad hoc RNA-trajectory analysis, and drug-gene interaction prediction. PlaqView serves as one of the largest central repositories for cardiovascular single-cell datasets, which now includes data from human aortic aneurysm, gene-specific mouse knockouts, and healthy references. PlaqView 2.0 brings advanced tools and high-performance computing directly to users without the need for any programming knowledge. Lastly, we outline steps to generalize and repurpose PlaqView's framework for single-cell datasets from other fields.

The HDAC9-associated risk locus promotes coronary artery disease by governing TWIST1

Genome wide association studies (GWAS) have identified thousands of single nucleotide polymorphisms (SNPs) associated with the risk of common disorders. However, since the large majority of these risk SNPs reside outside gene-coding regions, GWAS generally provide no information about causal mechanisms regarding the specific gene(s) that are affected or the tissue(s) in which these candidate gene(s) exert their effect. The ‘gold standard’ method for understanding causal genes and their mechanisms of action are laborious basic science studies often involving sophisticated knockin or knockout mouse lines, however, these types of studies are impractical as a high-throughput means to understand the many risk variants that cause complex diseases like coronary artery disease (CAD). As a solution, we developed a streamlined, data-driven informatics pipeline to gain mechanistic insights on complex genetic loci. The pipeline begins by understanding the SNPs in a given locus in terms of their relative location and linkage disequilibrium relationships, and then identifies nearby expression quantitative trait loci (eQTLs) to determine their relative independence and the likely tissues that mediate their disease-causal effects. The pipeline then seeks to understand associations with other disease-relevant genes, disease sub-phenotypes, potential causality (Mendelian randomization), and the regulatory and functional involvement of these genes in gene regulatory co-expression networks (GRNs). Here, we applied this pipeline to understand a cluster of SNPs associated with CAD within and immediately adjacent to the gene encoding HDAC9. Our pipeline demonstrated, and validated, that this locus is causal for CAD by modulation of TWIST1 expression levels in the arterial wall, and by also governing a GRN related to metabolic function in skeletal muscle. Our results reconciled numerous prior studies, and also provided clear evidence that this locus does not govern HDAC9 expression, structure or function. This pipeline should be considered as a powerful and efficient way to understand GWAS risk loci in a manner that better reflects the highly complex nature of genetic risk associated with common disorders.

Single-nucleus chromatin accessibility profiling highlights regulatory mechanisms of coronary artery disease risk

Coronary artery disease (CAD) is a complex inflammatory disease involving genetic influences across cell types. Genome-wide association studies (GWAS) have identified over 200 loci associated with CAD, where the majority of risk variants reside in noncoding DNA sequences impacting cis-regulatory elements (CREs). Here, we applied single-nucleus ATAC-seq to profile 28,316 nuclei across coronary artery segments from 41 patients with varying stages of CAD, which revealed 14 distinct cellular clusters. We mapped ~320,000 accessible sites across all cells, identified cell type-specific elements, transcription factors, and prioritized functional CAD risk variants. . We identified elements in smooth muscle cell (SMC) transition states (e.g. fibromyocytes) and functional variants predicted to alter SMC and macrophage-specific regulation of MRAS (3q22) and LIPA (10q23), respectively. We further nominated key driver transcription factors such as PRDM16 and TBX2. Together, this single nucleus atlas provides a critical step towards interpreting regulatory mechanisms across the continuum of CAD risk.

Integrative Prioritization of Causal Genes for Coronary Artery Disease

BACKGROUND

Hundreds 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.

METHODS

We 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).

RESULTS

We 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.

CONCLUSIONS

We 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.

Intersecting single-cell transcriptomics and genome-wide association studies identifies crucial cell populations and candidate genes for atherosclerosis

Aims

Genome-wide association studies (GWASs) have discovered hundreds of common genetic variants for atherosclerotic disease and cardiovascular risk factors. The translation of susceptibility loci into biological mechanisms and targets for drug discovery remains challenging. Intersecting genetic and gene expression data has led to the identification of candidate genes. However, previously studied tissues are often non-diseased and heterogeneous in cell composition, hindering accurate candidate prioritization. Therefore, we analysed single-cell transcriptomics from atherosclerotic plaques for cell-type-specific expression to identify atherosclerosis-associated candidate gene–cell pairs.

Methods and results

We applied gene-based analyses using GWAS summary statistics from 46 atherosclerotic and cardiovascular disease, risk factors, and other traits. We then intersected these candidates with single-cell RNA sequencing (scRNA-seq) data to identify genes specific for individual cell (sub)populations in atherosclerotic plaques. The coronary artery disease (CAD) loci demonstrated a prominent signal in plaque smooth muscle cells (SMCs) (SKI, KANK2, and SORT1) P-adj. = 0.0012, and endothelial cells (ECs) (SLC44A1, ATP2B1) P-adj. = 0.0011. Finally, we used liver-derived scRNA-seq data and showed hepatocyte-specific enrichment of genes involved in serum lipid levels.

Conclusion

We discovered novel and known gene–cell pairs pointing to new biological mechanisms of atherosclerotic disease. We highlight that loci associated with CAD reveal prominent association levels in mainly plaque SMC and EC populations. We present an intuitive single-cell transcriptomics-driven workflow rooted in human large-scale genetic studies to identify putative candidate genes and affected cells associated with cardiovascular traits. Collectively, our workflow allows for the identification of cell-specific targets relevant for atherosclerosis and can be universally applied to other complex genetic diseases and traits.

Enhanced single-cell RNA-seq workflow reveals coronary artery disease cellular cross-talk and candidate drug targets

BACKGROUND AND AIMS

The atherosclerotic plaque microenvironment is highly complex, and selective agents that modulate plaque stability are not yet available. We sought to develop a scRNA-seq analysis workflow to investigate this environment and uncover potential therapeutic approaches. We designed a user-friendly, reproducible workflow that will be applicable to other disease-specific scRNA-seq datasets.

METHODS

Here we incorporated automated cell labeling, pseudotemporal ordering, ligand-receptor evaluation, and drug-gene interaction analysis into a ready-to-deploy workflow. We applied this pipeline to further investigate a previously published human coronary single-cell dataset by Wirka et al. Notably, we developed an interactive web application to enable further exploration and analysis of this and other cardiovascular single-cell datasets.

RESULTS

We revealed distinct derivations of fibroblast-like cells from smooth muscle cells (SMCs), and showed the key changes in gene expression along their de-differentiation path. We highlighted several key ligand-receptor interactions within the atherosclerotic environment through functional expression profiling and revealed several avenues for future pharmacological development for precision medicine. Further, our interactive web application, PlaqView (www.plaqview.com), allows lay scientists to explore this and other datasets and compare scRNA-seq tools without prior coding knowledge.

CONCLUSIONS

This publicly available workflow and application will allow for more systematic and user-friendly analysis of scRNA datasets in other disease and developmental systems. Our analysis pipeline provides many hypothesis-generating tools to unravel the etiology of coronary artery disease. We also highlight potential mechanisms for several drugs in the atherosclerotic cellular environment. Future releases of PlaqView will feature more scRNA-seq and scATAC-seq atherosclerosis-related datasets to provide a critical resource for the field, and to promote data harmonization and biological interpretation.

Dynamic changes in chromatin accessibility are associated with the atherogenic transitioning of vascular smooth muscle cells

AIMS

De-differentiation and activation of pro-inflammatory pathways are key transitions vascular smooth muscle cells (SMCs) make during atherogenesis. Here, we explored the upstream regulators of this 'atherogenic transition'.

METHODS AND RESULTS

Genome-wide sequencing studies, including ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) and RNA-seq, were performed on cells isolated from both murine SMC-lineage tracing models of atherosclerosis and human atherosclerotic lesions. At the bulk level, alterations in chromatin accessibility were associated with the atherogenic transitioning of lesional SMCs, especially in relation to genes that govern differentiation status and complement-dependent inflammation. Using computational biology, we observed that a transcription factor previously related to coronary artery disease, ATF3 (Activating transcription factor 3), was predicted to be an upstream regulator of genes altered during the transition. At the single-cell level, our results indicated that ATF3 is a key repressor of SMC transitioning towards the subset of cells that promote vascular inflammation by activating the complement cascade. The expression of ATF3 and complement component C3 were negatively correlated in SMCs from human atherosclerotic lesions, suggesting translational relevance. Phenome-wide association studies indicated that genetic variation that results in reduced expression of ATF3 is correlated with an increased risk for atherosclerosis, and the expression of ATF3 was significantly downregulated in humans with advanced vascular disease.

CONCLUSION

Our study indicates that the plasticity of atherosclerotic SMCs may in part be explained by dynamic changes in their chromatin architecture, which in turn may contribute to their maladaptive response to inflammation-induced stress.

TRANSLATIONAL PERSPECTIVE

The recent CANTOS and COLCOT trials have shown that targeting inflammatory pathways lowers the risk of major adverse cardiovascular events. However, more specific targets are needed to avoid immunosuppressive side effects. Our data identify an upstream regulator of pro-inflammatory SMCs, ATF3, which is involved in the initial atherogenic transitioning of lesional SMCs. Restoring ATF3 activity may prevent the de-differentiation of SMCs and offer a novel translational approach for the suppression of complement-dependent inflammation in atherosclerotic lesions.

Vascular smooth muscle cell dysfunction contribute to neuroinflammation and Tau hyperphosphorylation in Alzheimer disease

Despite the emerging evidence implying early vascular contributions to neurodegenerative syndromes, the role of vascular smooth muscle cells (VSMCs) in the pathogenesis of Alzheimer disease (AD) is still not well understood. Herein, we show that VSMCs in brains of patients with AD and animal models of the disease are deficient in multiple VSMC contractile markers which correlated with Tau accumulation in brain arterioles. Ex vivo and in vitro experiments demonstrated that VSMCs undergo dramatic phenotypic transitions under AD-like conditions, adopting pro-inflammatory phenotypes. Notably, these changes coincided with Tau hyperphosphorylation at residues Y18, T205, and S262. We also observed that VSMC dysfunction occurred in an age-dependent manner and that expression of Sm22α protein was inversely correlated with CD68 and Tau expression in brain arterioles of the 3xTg-AD and 5xFAD mice. Together, these findings further support the contribution of dysfunctional VSMCs in AD pathogenesis and nominate VSMCs as a potential therapeutic target in AD.

•

Loss of VSMC contractile phenotypes correlates with Tau accumulation in brain arterioles

•

VSMC dysfunction promotes the hyperphosphorylation of Tau protein at multiple residues

•

VSMC dysfunction occurs in an age-dependent manner in brain arterioles of patients with AD

•

Vascular smooth muscle cell is a promising therapeutic target in AD

Abstract MP39: FHL5 , A Novel Cofactor Associated With Coronary Artery Disease, Regulates Smooth Muscle Cell Function Through A Transcriptional Network Linking Multiple Risk Loci

Introduction: Coronary artery disease (CAD) is the leading cause of death worldwide with an estimated heritability of ~50%. The most recent genome wide association study (GWAS) for CAD and myocardial infarction (MI) identified over 200 loci. One such novel CAD/MI locus, UFL1-FHL5 (P = 1.1E-8), is associated with additional vascular pathologies, including hypertension, migraines, and coronary calcification. We previously showed through statistical fine-mapping that FHL5 is the top candidate causal gene underlying each of these vascular trait associations and highly enriched in the contractile mural cell populations in the artery. These preliminary studies motivate our hypothesis that FHL5 functions as a transcriptional regulator of SMC contractility to affect vascular disease risk.

Methods and Results: Given the reported role of FHL5 as a cofactor, we mapped 17,201 candidate FHL5 binding sites in coronary artery SMCs using the Cleavage Under Targets and Release Using Nuclease (CUT&RUN) method. Binding sites were enriched for AP-1 family motifs and strongly overlapped CREB binding sites and H3K27ac enhancer marks. FHL5 target genes were functionally enriched in well characterized CAD pathways, such as extracellular matrix organization (P=5.2E-15, OR=2.1) and TGF-beta signaling (P=8.0E-5, OR=2.0). Interestingly, FHL5 binding sites were also enriched for CAD (P=1.8E-7, OR=2.3) and blood pressure (BP) risk variants (P=0.04, OR=1.2), thereby linking FHL5 with the regulation of multiple downstream CAD/BP loci. Weighted gene co-expression network and key driver analyses of human coronary artery transcriptomic data further support this regulatory role. Lastly, we validated these findings by performing confocal based calcium imaging and collagen gel contraction assays. Consistent with our functional genomic analyses, overexpression of FHL5 elevated intracellular calcium levels 2.5X and increased SMC contractility by ~45%.

Conclusion: Taken together, our results provide evidence that FHL5 may impact CAD risk by regulating a network of disease-associated genes mediating SMC functions. These findings further contribute to our understanding of the heritable risk for multiple vascular diseases.

Abstract 113: Cell-specific Chromatin Landscape Of Human Coronary Artery Resolves Mechanisms Of Disease Risk

Coronary artery disease (CAD) is a complex inflammatory disease involving genetic influences across several cell types. Genome-wide association studies (GWAS) have identified over 170 loci associated with CAD, where the majority of risk variants reside in noncoding DNA sequences impacting cis -regulatory elements (CREs). Here, we applied single-cell ATAC-seq to profile 28,316 cells across coronary artery segments from 41 patients with varying stages of CAD, which revealed 14 distinct cellular clusters. We mapped over 320,000 accessible sites across all cells, identified cell type-specific elements, transcription factors, and prioritized functional CAD risk variants via quantitative trait locus and sequence-based predictive modeling. Using differential peak analyses we identified a number of candidate mechanisms for smooth muscle cell transition states (e.g. fibromyocytes). By integrating these profiles with GWAS meta-analysis summary data we resolved cell type-specific putative binding sites for the majority of CAD risk variants. In particular, we prioritized functional variants predicted to alter MEF2 binding in smooth muscle cells at the MRAS locus. We also identify variants predicted to alter macrophage-specific regulation of LIPA. We further employed DNA to gene linkage to nominate disease-associated key driver transcription factors such as PRDM16 and TBX2. Together, this single cell atlas provides a critical step towards interpreting cis -regulatory mechanisms in the vessel wall across the continuum of CAD risk.

Genome-wide analysis identifies novel susceptibility loci for myocardial infarction

AIMS

While most patients with myocardial infarction (MI) have underlying coronary atherosclerosis, not all patients with coronary artery disease (CAD) develop MI. We sought to address the hypothesis that some of the genetic factors which establish atherosclerosis may be distinct from those that predispose to vulnerable plaques and thrombus formation.

METHODS AND RESULTS

We carried out a genome-wide association study for MI in the UK Biobank (n∼472 000), followed by a meta-analysis with summary statistics from the CARDIoGRAMplusC4D Consortium (n∼167 000). Multiple independent replication analyses and functional approaches were used to prioritize loci and evaluate positional candidate genes. Eight novel regions were identified for MI at the genome wide significance level, of which effect sizes at six loci were more robust for MI than for CAD without the presence of MI. Confirmatory evidence for association of a locus on chromosome 1p21.3 harbouring choline-like transporter 3 (SLC44A3) with MI in the context of CAD, but not with coronary atherosclerosis itself, was obtained in Biobank Japan (n∼165 000) and 16 independent angiography-based cohorts (n∼27 000). Follow-up analyses did not reveal association of the SLC44A3 locus with CAD risk factors, biomarkers of coagulation, other thrombotic diseases, or plasma levels of a broad array of metabolites, including choline, trimethylamine N-oxide, and betaine. However, aortic expression of SLC44A3 was increased in carriers of the MI risk allele at chromosome 1p21.3, increased in ischaemic (vs. non-diseased) coronary arteries, up-regulated in human aortic endothelial cells treated with interleukin-1β (vs. vehicle), and associated with smooth muscle cell migration in vitro.

CONCLUSIONS

A large-scale analysis comprising ∼831 000 subjects revealed novel genetic determinants of MI and implicated SLC44A3 in the pathophysiology of vulnerable plaques.

APOL1 Genetic Variants Are Associated With Increased Risk of Coronary Atherosclerotic Plaque Rupture in the Black Population

[Figure: see text].

Comparison of perometry-based volumetric arm measurements and bioimpedance spectroscopy for early identification of lymphedema in a prospectively-screened cohort of breast cancer patients

Breast cancer-related lymphedema (BCRL) affects more than one in five women treated for breast cancer, and women remain at lifelong risk. Screening for BCRL is recommended by several national and international organizations for women at risk of BCRL, and multiple methods of objective screening measurement exist. The goal of this study was to compare the use of perometry and bioimpedance spectroscopy (BIS) for early identification of BCRL in a cohort of 138 prospectivelyscreened patients. At each screening visit, a patient's relative volume change (RVC) from perometer measurements and change in L-Dex from baseline (ΔL-Dex) using BIS was calculated. There was a negligible correlation between RVC and ΔL-Dex (r=0.195). Multiple thresholds of BCRL were examined: RVC ≥5% and ≥10% as well as and ΔL-Dex ≥6.5 and ≥10. While some patients developed an elevated RVC and ΔL-Dex, many demonstrated elevations in only one threshold category. Moreover, the majority of patients with RVC ≥5%, ΔL-Dex ≥6.5, or ΔL-Dex ≥10 regressed to non-elevated measurements without intervention. These findings suggest a role for combining multiple screening methods for early identification of BCRL; furthermore, BCRL diagnosis must incorporate patient symptoms and clinical evaluation with objective measurements obtained from techniques such as perometry and bioimpedance spectroscopy.

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

RATIONALE

Coronary artery disease (CAD) is a major cause of morbidity and mortality worldwide. Recent genome-wide association studies revealed 163 loci associated with CAD. However, the precise molecular mechanisms by which the majority of these loci increase CAD risk are not known. Vascular smooth muscle cells (VSMCs) are critical in the development of CAD. They can play either beneficial or detrimental roles in lesion pathogenesis, depending on the nature of their phenotypic changes.

OBJECTIVE

To identify genetic variants associated with atherosclerosis-relevant phenotypes in VSMCs.

METHODS AND RESULTS

We quantified 12 atherosclerosis-relevant phenotypes related to calcification, proliferation, and migration in VSMCs isolated from 151 multiethnic heart transplant donors. After genotyping and imputation, we performed association mapping using 6.3 million genetic variants. We demonstrated significant variations in calcification, proliferation, and migration. These phenotypes were not correlated with each other. We performed genome-wide association studies for 12 atherosclerosis-relevant phenotypes and identified 4 genome-wide significant loci associated with at least one VSMC phenotype. We overlapped the previously identified CAD loci with our data set and found nominally significant associations at 79 loci. One of them was the chromosome 1q41 locus, which harbors MIA3. The G allele of the lead risk single nucleotide polymorphism (SNP) rs67180937 was associated with lower VSMC MIA3 expression and lower proliferation. Lentivirus-mediated silencing of MIA3 (melanoma inhibitory activity protein 3) in VSMCs resulted in lower proliferation, consistent with human genetics findings. Furthermore, we observed a significant reduction of MIA3 protein in VSMCs in thin fibrous caps of late-stage atherosclerotic plaques compared to early fibroatheroma with thick and protective fibrous caps in mice and humans.

CONCLUSIONS

Our data demonstrate that genetic variants have significant influences on VSMC function relevant to the development of atherosclerosis. Furthermore, high MIA3 expression may promote atheroprotective VSMC phenotypic transitions, including increased proliferation, which is essential in the formation or maintenance of a protective fibrous cap.

Clonally expanding smooth muscle cells promote atherosclerosis by escaping efferocytosis and activating the complement cascade

Atherosclerosis is the process underlying heart attack and stroke. Despite decades of research, its pathogenesis remains unclear. Dogma suggests that atherosclerotic plaques expand primarily via the accumulation of cholesterol and inflammatory cells. However, recent evidence suggests that a substantial portion of the plaque may arise from a subset of "dedifferentiated" vascular smooth muscle cells (SMCs) which proliferate in a clonal fashion. Herein we use multicolor lineage-tracing models to confirm that the mature SMC can give rise to a hyperproliferative cell which appears to promote inflammation via elaboration of complement-dependent anaphylatoxins. Despite being extensively opsonized with prophagocytic complement fragments, we find that this cell also escapes immune surveillance by neighboring macrophages, thereby exacerbating its relative survival advantage. Mechanistic studies indicate this phenomenon results from a generalized opsonin-sensing defect acquired by macrophages during polarization. This defect coincides with the noncanonical up-regulation of so-called don't eat me molecules on inflamed phagocytes, which reduces their capacity for programmed cell removal (PrCR). Knockdown or knockout of the key antiphagocytic molecule CD47 restores the ability of macrophages to sense and clear opsonized targets in vitro, allowing for potent and targeted suppression of clonal SMC expansion in the plaque in vivo. Because integrated clinical and genomic analyses indicate that similar pathways are active in humans with cardiovascular disease, these studies suggest that the clonally expanding SMC may represent a translational target for treating atherosclerosis.

Genetic Insights Into Smooth Muscle Cell Contributions to Coronary Artery Disease

Coronary artery disease is a complex cardiovascular disease involving an interplay of genetic and environmental influences over a lifetime. Although considerable progress has been made in understanding lifestyle risk factors, genetic factors identified from genome-wide association studies may capture additional hidden risk undetected by traditional clinical tests. These genetic discoveries have highlighted many candidate genes and pathways dysregulated in the vessel wall, including those involving smooth muscle cell phenotypic modulation and injury responses. Here, we summarize experimental evidence for a few genome-wide significant loci supporting their roles in smooth muscle cell biology and disease. We also discuss molecular quantitative trait locus mapping as a powerful discovery and fine-mapping approach applied to smooth muscle cell and coronary artery disease-relevant tissues. We emphasize the critical need for alternative genetic strategies, including cis/trans-regulatory network analysis, genome editing, and perturbations, as well as single-cell sequencing in smooth muscle cell tissues and model organisms, under both normal and disease states. By integrating multiple experimental and analytical modalities, these multidimensional datasets should improve the interpretation of coronary artery disease genome-wide association studies and molecular quantitative trait locus signals and inform candidate targets for therapeutic intervention or risk prediction.

PCSK6 Is a Key Protease in the Control of Smooth Muscle Cell Function in Vascular Remodeling

RATIONALE

PCSKs (Proprotein convertase subtilisins/kexins) are a protease family with unknown functions in vasculature. Previously, we demonstrated PCSK6 upregulation in human atherosclerotic plaques associated with smooth muscle cells (SMCs), inflammation, extracellular matrix remodeling, and mitogens.

OBJECTIVE

Here, we applied a systems biology approach to gain deeper insights into the PCSK6 role in normal and diseased vessel wall.

METHODS AND RESULTS

Genetic analyses revealed association of intronic PCSK6 variant rs1531817 with maximum internal carotid intima-media thickness progression in high-cardiovascular risk subjects. This variant was linked with PCSK6 mRNA expression in healthy aortas and plaques but also with overall plaque SMA+ cell content and pericyte fraction. Increased PCSK6 expression was found in several independent human cohorts comparing atherosclerotic lesions versus healthy arteries, using transcriptomic and proteomic datasets. By immunohistochemistry, PCSK6 was localized to fibrous cap SMA+ cells and neovessels in plaques. In human, rat, and mouse intimal hyperplasia, PCSK6 was expressed by proliferating SMA+ cells and upregulated after 5 days in rat carotid balloon injury model, with positive correlation to PDGFB (platelet-derived growth factor subunit B) and MMP (matrix metalloprotease) 2/MMP14. Here, PCSK6 was shown to colocalize and cointeract with MMP2/MMP14 by in situ proximity ligation assay. Microarrays of carotid arteries from Pcsk6^-/- versus control mice revealed suppression of contractile SMC markers, extracellular matrix remodeling enzymes, and cytokines/receptors. Pcsk6^-/- mice showed reduced intimal hyperplasia response upon carotid ligation in vivo, accompanied by decreased MMP14 activation and impaired SMC outgrowth from aortic rings ex vivo. PCSK6 silencing in human SMCs in vitro leads to downregulation of contractile markers and increase in MMP2 expression. Conversely, PCSK6 overexpression increased PDGFBB (platelet-derived growth factor BB)-induced cell proliferation and particularly migration.

CONCLUSIONS

PCSK6 is a novel protease that induces SMC migration in response to PDGFB, mechanistically via modulation of contractile markers and MMP14 activation. This study establishes PCSK6 as a key regulator of SMC function in vascular remodeling. Visual Overview: An online visual overview is available for this article.

BART: a transcription factor prediction tool with query gene sets or epigenomic profiles

Summary

Identification of functional transcription factors that regulate a given gene set is an important problem in gene regulation studies. Conventional approaches for identifying transcription factors, such as DNA sequence motif analysis, are unable to predict functional binding of specific factors and not sensitive enough to detect factors binding at distal enhancers. Here, we present binding analysis for regulation of transcription (BART), a novel computational method and software package for predicting functional transcription factors that regulate a query gene set or associate with a query genomic profile, based on more than 6000 existing ChIP-seq datasets for over 400 factors in human or mouse. This method demonstrates the advantage of utilizing publicly available data for functional genomics research.

Availability and implementation

BART is implemented in Python and available at http://faculty.virginia.edu/zanglab/bart.

Supplementary information

Supplementary data are available at Bioinformatics online.

The novel coronary artery disease risk gene JCAD/KIAA1462 promotes endothelial dysfunction and atherosclerosis

AIMS

Recent genome-wide association studies (GWAS) have identified that the JCAD locus is associated with risk of coronary artery disease (CAD) and myocardial infarction (MI). However, the mechanisms whereby candidate gene JCAD confers disease risk remain unclear. We addressed whether and how JCAD affects the development of atherosclerosis, the common cause of CAD.

METHODS AND RESULTS

By mining data in the Genotype-Tissue Expression (GTEx) database, we found that CAD-associated risk variants at the JCAD locus are linked to increased JCAD gene expression in human arteries, implicating JCAD as a candidate causal CAD gene. We therefore generated global and endothelial cell (EC) specific-JCAD knockout mice, and observed that JCAD deficiency attenuated high fat diet-induced atherosclerosis in ApoE-deficient mice. JCAD-deficiency in mice also improved endothelium-dependent relaxation. Genome-wide transcriptional profiling of JCAD-depleted human coronary artery ECs showed that JCAD depletion inhibited the activation of YAP/TAZ pathway, and the expression of downstream pro-atherogenic genes, including CTGF and Cyr61. As a result, JCAD-deficient ECs attracted fewer monocytes in response to lipopolysaccharide (LPS) stimulation. Moreover, JCAD expression in ECs was decreased under unidirectional laminar flow in vitro and in vivo. Proteomics studies suggest that JCAD regulates YAP/TAZ activation by interacting with actin-binding protein TRIOBP, thereby stabilizing stress fiber formation. Finally, we observed that endothelial JCAD expression was increased in mouse and human atherosclerotic plaques.

CONCLUSION

The present study demonstrates that the GWAS-identified CAD risk gene JCAD promotes endothelial dysfunction and atherosclerosis, thus highlighting the possibility of new therapeutic strategies for CAD by targeting JCAD.

Inducible Pluripotent Stem Cell-Derived Cardiomyocytes Reveal Aberrant Extracellular Regulated Kinase 5 and Mitogen-Activated Protein Kinase Kinase 1/2 Signaling Concomitantly Promote Hypertrophic Cardiomyopathy in RAF1-Associated Noonan Syndrome

BACKGROUND

More than 90% of individuals with Noonan syndrome (NS) with mutations clustered in the CR2 domain of RAF1 present with severe and often lethal hypertrophic cardiomyopathy (HCM). The signaling pathways by which NS RAF1 mutations promote HCM remain elusive, and so far, there is no known treatment for NS-associated HCM.

METHODS

We used patient-derived RAF1^S257L/+ and CRISPR-Cas9-generated isogenic control inducible pluripotent stem cell (iPSC)-derived cardiomyocytes to model NS RAF1-associated HCM and to further delineate the molecular mechanisms underlying the disease.

RESULTS

We show that mutant iPSC-derived cardiomyocytes phenocopy the pathology seen in hearts of patients with NS by exhibiting hypertrophy and structural defects. Through pharmacological and genetic targeting, we identify 2 perturbed concomitant pathways that, together, mediate HCM in RAF1 mutant iPSC-derived cardiomyocytes. Hyperactivation of mitogen-activated protein kinase kinase 1/2 (MEK1/2), but not extracellular regulated kinase 1/2, causes myofibrillar disarray, whereas the enlarged cardiomyocyte phenotype is a direct consequence of increased extracellular regulated kinase 5 (ERK5) signaling, a pathway not previously known to be involved in NS. RNA-sequencing reveals genes with abnormal expression in RAF1 mutant iPSC-derived cardiomyocytes and identifies subsets of genes dysregulated by aberrant MEK1/2 or ERK5 pathways that could contribute to the NS-associated HCM.

CONCLUSIONS

Taken together, the results of our study identify the molecular mechanisms by which NS RAF1 mutations cause HCM and reveal downstream effectors that could serve as therapeutic targets for treatment of NS and perhaps other, more common, congenital HCM disorders.

TCF21 and AP-1 interact through epigenetic modifications to regulate coronary artery disease gene expression

Background

Genome-wide association studies have identified over 160 loci that are associated with coronary artery disease. As with other complex human diseases, risk in coronary disease loci is determined primarily by altered expression of the causal gene, due to variation in binding of transcription factors and chromatin-modifying proteins that directly regulate the transcriptional apparatus. We have previously identified a coronary disease network downstream of the disease-associated transcription factor TCF21, and in work reported here extends these studies to investigate the mechanisms by which it interacts with the AP-1 transcription complex to regulate local epigenetic effects in these downstream coronary disease loci.

Methods

Genomic studies, including chromatin immunoprecipitation sequencing, RNA sequencing, and protein-protein interaction studies, were performed in human coronary artery smooth muscle cells.

Results

We show here that TCF21 and JUN regulate expression of two presumptive causal coronary disease genes, SMAD3 and CDKN2B-AS1, in part by interactions with histone deacetylases and acetyltransferases. Genome-wide TCF21 and JUN binding is jointly localized and particularly enriched in coronary disease loci where they broadly modulate H3K27Ac and chromatin state changes linked to disease-related processes in vascular cells. Heterozygosity at coronary disease causal variation, or genome editing of these variants, is associated with decreased binding of both JUN and TCF21 and loss of expression in cis, supporting a transcriptional mechanism for disease risk.

Conclusions

These data show that the known chromatin remodeling and pioneer functions of AP-1 are a pervasive aspect of epigenetic control of transcription, and thus, the risk in coronary disease-associated loci, and that interaction of AP-1 with TCF21 to control epigenetic features, contributes to the genetic risk in loci where they co-localize.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0635-9) contains supplementary material, which is available to authorized users.

Multi-Omics Approaches to Study Long Non-coding RNA Function in Atherosclerosis

Atherosclerosis is a complex inflammatory disease of the vessel wall involving the interplay of multiple cell types including vascular smooth muscle cells, endothelial cells, and macrophages. Large-scale genome-wide association studies (GWAS) and the advancement of next generation sequencing technologies have rapidly expanded the number of long non-coding RNA (lncRNA) transcripts predicted to play critical roles in the pathogenesis of the disease. In this review, we highlight several lncRNAs whose functional role in atherosclerosis is well-documented through traditional biochemical approaches as well as those identified through RNA-sequencing and other high-throughput assays. We describe novel genomics approaches to study both evolutionarily conserved and divergent lncRNA functions and interactions with DNA, RNA, and proteins. We also highlight assays to resolve the complex spatial and temporal regulation of lncRNAs. Finally, we summarize the latest suite of computational tools designed to improve genomic and functional annotation of these transcripts in the human genome. Deep characterization of lncRNAs is fundamental to unravel coronary atherosclerosis and other cardiovascular diseases, as these regulatory molecules represent a new class of potential therapeutic targets and/or diagnostic markers to mitigate both genetic and environmental risk factors.

Functional regulatory mechanism of smooth muscle cell-restricted LMOD1 coronary artery disease locus

Recent genome-wide association studies (GWAS) have identified multiple new loci which appear to alter coronary artery disease (CAD) risk via arterial wall-specific mechanisms. One of the annotated genes encodes LMOD1 (Leiomodin 1), a member of the actin filament nucleator family that is highly enriched in smooth muscle-containing tissues such as the artery wall. However, it is still unknown whether LMOD1 is the causal gene at this locus and also how the associated variants alter LMOD1 expression/function and CAD risk. Using epigenomic profiling we recently identified a non-coding regulatory variant, rs34091558, which is in tight linkage disequilibrium (LD) with the lead CAD GWAS variant, rs2820315. Herein we demonstrate through expression quantitative trait loci (eQTL) and statistical fine-mapping in GTEx, STARNET, and human coronary artery smooth muscle cell (HCASMC) datasets, rs34091558 is the top regulatory variant for LMOD1 in vascular tissues. Position weight matrix (PWM) analyses identify the protective allele rs34091558-TA to form a conserved Forkhead box O3 (FOXO3) binding motif, which is disrupted by the risk allele rs34091558-A. FOXO3 chromatin immunoprecipitation and reporter assays show reduced FOXO3 binding and LMOD1 transcriptional activity by the risk allele, consistent with effects of FOXO3 downregulation on LMOD1. LMOD1 knockdown results in increased proliferation and migration and decreased cell contraction in HCASMC, and immunostaining in atherosclerotic lesions in the SMC lineage tracing reporter mouse support a key role for LMOD1 in maintaining the differentiated SMC phenotype. These results provide compelling functional evidence that genetic variation is associated with dysregulated LMOD1 expression/function in SMCs, together contributing to the heritable risk for CAD.

Coronary artery disease genes SMAD3 and TCF21 promote opposing interactive genetic programs that regulate smooth muscle cell differentiation and disease risk

Although numerous genetic loci have been associated with coronary artery disease (CAD) with genome wide association studies, efforts are needed to identify the causal genes in these loci and link them into fundamental signaling pathways. Recent studies have investigated the disease mechanism of CAD associated gene SMAD3, a central transcription factor (TF) in the TGFβ pathway, investigating its role in smooth muscle biology. In vitro studies in human coronary artery smooth muscle cells (HCASMC) revealed that SMAD3 modulates cellular phenotype, promoting expression of differentiation marker genes while inhibiting proliferation. RNA sequencing and chromatin immunoprecipitation sequencing studies in HCASMC identified downstream genes that reside in pathways which mediate vascular development and atherosclerosis processes in this cell type. HCASMC phenotype, and gene expression patterns promoted by SMAD3 were noted to have opposing direction of effect compared to another CAD associated TF, TCF21. At sites of SMAD3 and TCF21 colocalization on DNA, SMAD3 binding was inversely correlated with TCF21 binding, due in part to TCF21 locally blocking chromatin accessibility at the SMAD3 binding site. Further, TCF21 was able to directly inhibit SMAD3 activation of gene expression in transfection reporter gene studies. In contrast to TCF21 which is protective toward CAD, SMAD3 expression in HCASMC was shown to be directly correlated with disease risk. We propose that the pro-differentiation action of SMAD3 inhibits dedifferentiation that is required for HCASMC to expand and stabilize disease plaque as they respond to vascular stresses, counteracting the protective dedifferentiating activity of TCF21 and promoting disease risk.

Genetic Regulatory Mechanisms of Smooth Muscle Cells Map to Coronary Artery Disease Risk Loci

Coronary artery disease (CAD) is the leading cause of death globally. Genome-wide association studies (GWASs) have identified more than 95 independent loci that influence CAD risk, most of which reside in non-coding regions of the genome. To interpret these loci, we generated transcriptome and whole-genome datasets using human coronary artery smooth muscle cells (HCASMCs) from 52 unrelated donors, as well as epigenomic datasets using ATAC-seq on a subset of 8 donors. Through systematic comparison with publicly available datasets from GTEx and ENCODE projects, we identified transcriptomic, epigenetic, and genetic regulatory mechanisms specific to HCASMCs. We assessed the relevance of HCASMCs to CAD risk using transcriptomic and epigenomic level analyses. By jointly modeling eQTL and GWAS datasets, we identified five genes (SIPA1, TCF21, SMAD3, FES, and PDGFRA) that may modulate CAD risk through HCASMCs, all of which have relevant functional roles in vascular remodeling. Comparison with GTEx data suggests that SIPA1 and PDGFRA influence CAD risk predominantly through HCASMCs, while other annotated genes may have multiple cell and tissue targets. Together, these results provide tissue-specific and mechanistic insights into the regulation of a critical vascular cell type associated with CAD in human populations.

GWAS Reveal Targets in Vessel Wall Pathways to Treat Coronary Artery Disease

Coronary artery disease (CAD) is the leading cause of mortality worldwide and poses a considerable public health burden. Recent genome-wide association studies (GWAS) have revealed >100 genetic loci associated with CAD susceptibility in humans. While a number of these loci harbor gene targets of currently approved therapies, such as statins and PCSK9 inhibitors, the majority of the annotated genes at these loci encode for proteins involved in vessel wall function with no known drugs available. Importantly many of the associated genes linked to vascular (smooth muscle, endothelial, and macrophage) cell processes are now organized into distinct functional pathways, e.g., vasodilation, growth factor responses, extracellular matrix and plaque remodeling, and inflammation. In this mini-review, we highlight the most recently identified loci that have predicted roles in the vessel wall and provide genetic context for pre-existing therapies as well as new drug targets informed from GWAS. With the development of new modalities to target these pathways, (e.g., antisense oligonucleotides, CRISPR/Cas9, and RNA interference) as well as the computational frameworks to prioritize or reposition therapeutics, there is great opportunity to close the gap from initial genetic discovery to clinical translation for many patients affected by this common disease.