EPAS1 Attenuates Atherosclerosis Initiation at Disturbed Flow Sites Through Endothelial Fatty Acid Uptake.

BACKGROUND: Atherosclerotic plaques form unevenly due to disturbed blood flow, causing localized endothelial cell (EC) dysfunction. Obesity exacerbates this process, but the underlying molecular mechanisms are unclear. The transcription factor EPAS1 (HIF2A) has regulatory roles in endothelium, but its involvement in atherosclerosis remains unexplored. This study investigates the potential interplay between EPAS1, obesity, and atherosclerosis. METHODS: Responses to shear stress were analyzed using cultured porcine aortic EC exposed to flow in vitro coupled with metabolic and molecular analyses and by en face immunostaining of murine aortic EC exposed to disturbed flow in vivo. Obesity and dyslipidemia were induced in mice via exposure to a high-fat diet or through Leptin gene deletion. The role of Epas1 in atherosclerosis was evaluated by inducible endothelial Epas1 deletion, followed by hypercholesterolemia induction (adeno-associated virus-PCSK9 [proprotein convertase subtilisin/kexin type 9]; high-fat diet). RESULTS: En face staining revealed EPAS1 enrichment at sites of disturbed blood flow that are prone to atherosclerosis initiation. Obese mice exhibited substantial reduction in endothelial EPAS1 expression. Sulforaphane, a compound with known atheroprotective effects, restored EPAS1 expression and concurrently reduced plasma triglyceride levels in obese mice. Consistently, triglyceride derivatives (free fatty acids) suppressed EPAS1 in cultured EC by upregulating the negative regulator PHD2. Clinical observations revealed that reduced serum EPAS1 correlated with increased endothelial PHD2 and PHD3 in obese individuals. Functionally, endothelial EPAS1 deletion increased lesion formation in hypercholesterolemic mice, indicating an atheroprotective function. Mechanistic insights revealed that EPAS1 protects arteries by maintaining endothelial proliferation by positively regulating the expression of the fatty acid-handling molecules CD36 (cluster of differentiation 36) and LIPG (endothelial type lipase G) to increase fatty acid beta-oxidation. CONCLUSIONS: Endothelial EPAS1 attenuates atherosclerosis at sites of disturbed flow by maintaining EC proliferation via fatty acid uptake and metabolism. This endothelial repair pathway is inhibited in obesity, suggesting a novel triglyceride-PHD2 modulation pathway suppressing EPAS1 expression. These findings have implications for therapeutic strategies addressing vascular dysfunction in obesity.

CD163+ Macrophages Induce Endothelial-to-Mesenchymal Transition in Atheroma.

BACKGROUND: Cell phenotype switching is increasingly being recognized in atherosclerosis. However, our understanding of the exact stimuli for such cellular transformations and their significance for human atherosclerosis is still evolving. Intraplaque hemorrhage is thought to be a major contributor to plaque progression in part by stimulating the influx of CD163+ macrophages. Here, we explored the hypothesis that CD163+ macrophages cause plaque progression through the induction of proapoptotic endothelial-to-mesenchymal transition (EndMT) within the fibrous cap. METHODS: Human coronary artery sections from CVPath's autopsy registry were selected for pathological analysis. Athero-prone ApoE-/- and ApoE-/-/CD163-/- mice were used for in vivo studies. Human peripheral blood mononuclear cell-induced macrophages and human aortic endothelial cells were used for in vitro experiments. RESULTS: In 107 lesions with acute coronary plaque rupture, 55% had pathological evidence of intraplaque hemorrhage in nonculprit vessels/lesions. Thinner fibrous cap, greater CD163+ macrophage accumulation, and a larger number of CD31/FSP-1 (fibroblast specific protein-1) double-positive cells and TUNEL (terminal deoxynucleotidyl transferase-dUTP nick end labeling) positive cells in the fibrous cap were observed in nonculprit intraplaque hemorrhage lesions, as well as in culprit rupture sections versus nonculprit fibroatheroma sections. Human aortic endothelial cells cultured with supernatants from hemoglobin/haptoglobin-exposed macrophages showed that increased mesenchymal marker proteins (transgelin and FSP-1) while endothelial markers (VE-cadherin and CD31) were reduced, suggesting EndMT induction. Activation of NF-κB (nuclear factor kappa ß) signaling by proinflammatory cytokines released from CD163+ macrophages directly regulated the expression of Snail, a critical transcription factor during EndMT induction. Western blot analysis for cleaved caspase-3 and microarray analysis of human aortic endothelial cells indicated that apoptosis was stimulated during CD163+ macrophage-induced EndMT. Additionally, CD163 deletion in athero-prone mice suggested that CD163 is required for EndMT and plaque progression. Using single-cell RNA sequencing from human carotid endarterectomy lesions, a population of EndMT was detected, which demonstrated significant upregulation of apoptosis-related genes. CONCLUSIONS: CD163+ macrophages provoke EndMT, which may promote plaque progression through fibrous cap thinning.

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.

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.

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.

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.

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.

Mitochondria-localized AMPK responds to local energetics and contributes to exercise and energetic stress-induced mitophagy.

Mitochondria form a complex, interconnected reticulum that is maintained through coordination among biogenesis, dynamic fission, and fusion and mitophagy, which are initiated in response to various cues to maintain energetic homeostasis. These cellular events, which make up mitochondrial quality control, act with remarkable spatial precision, but what governs such spatial specificity is poorly understood. Herein, we demonstrate that specific isoforms of the cellular bioenergetic sensor, 5' AMP-activated protein kinase (AMPKα1/α2/ß2/γ1), are localized on the outer mitochondrial membrane, referred to as mitoAMPK, in various tissues in mice and humans. Activation of mitoAMPK varies across the reticulum in response to energetic stress, and inhibition of mitoAMPK activity attenuates exercise-induced mitophagy in skeletal muscle in vivo. Discovery of a mitochondrial pool of AMPK and its local importance for mitochondrial quality control underscores the complexity of sensing cellular energetics in vivo that has implications for targeting mitochondrial energetics for disease treatment.

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.

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

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.

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

[Figure: see text].

CD47-blocking antibodies restore phagocytosis and prevent atherosclerosis.

Atherosclerosis is the disease process that underlies heart attack and stroke. Advanced lesions at risk of rupture are characterized by the pathological accumulation of diseased vascular cells and apoptotic cellular debris. Why these cells are not cleared remains unknown. Here we show that atherogenesis is associated with upregulation of CD47, a key anti-phagocytic molecule that is known to render malignant cells resistant to programmed cell removal, or 'efferocytosis'. We find that administration of CD47-blocking antibodies reverses this defect in efferocytosis, normalizes the clearance of diseased vascular tissue, and ameliorates atherosclerosis in multiple mouse models. Mechanistic studies implicate the pro-atherosclerotic factor TNF-α as a fundamental driver of impaired programmed cell removal, explaining why this process is compromised in vascular disease. Similar to recent observations in cancer, impaired efferocytosis appears to play a pathogenic role in cardiovascular disease, but is not a fixed defect and may represent a novel therapeutic target.

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.

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.

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.

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 RAF1S257L/+ 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.

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.

TCF21 and the environmental sensor aryl-hydrocarbon receptor cooperate to activate a pro-inflammatory gene expression program in coronary artery smooth muscle cells.

Both environmental factors and genetic loci have been associated with coronary artery disease (CAD), however gene-gene and gene-environment interactions that might identify molecular mechanisms of risk are not easily studied by human genetic approaches. We have previously identified the transcription factor TCF21 as the causal CAD gene at 6q23.2 and characterized its downstream transcriptional network that is enriched for CAD GWAS genes. Here we investigate the hypothesis that TCF21 interacts with a downstream target gene, the aryl hydrocarbon receptor (AHR), a ligand-activated transcription factor that mediates the cellular response to environmental contaminants, including dioxin and polycyclic aromatic hydrocarbons (e.g., tobacco smoke). Perturbation of TCF21 expression in human coronary artery smooth muscle cells (HCASMC) revealed that TCF21 promotes expression of AHR, its heterodimerization partner ARNT, and cooperates with these factors to upregulate a number of inflammatory downstream disease related genes including IL1A, MMP1, and CYP1A1. TCF21 was shown to bind in AHR, ARNT and downstream target gene loci, and co-localization was noted for AHR-ARNT and TCF21 binding sites genome-wide in regions of HCASMC open chromatin. These regions of co-localization were found to be enriched for GWAS signals associated with cardio-metabolic as well as chronic inflammatory disease phenotypes. Finally, we show that similar to TCF21, AHR gene expression is increased in atherosclerotic lesions in mice in vivo using laser capture microdissection, and AHR protein is localized in human carotid atherosclerotic lesions where it is associated with protein kinases with a critical role in innate immune response. These data suggest that TCF21 can cooperate with AHR to activate an inflammatory gene expression program that is exacerbated by environmental stimuli, and may contribute to the overall risk for CAD.