ZEB2 Shapes the Epigenetic Landscape of Atherosclerosis.

BACKGROUND: Smooth muscle cells (SMCs) transition into a number of different phenotypes during atherosclerosis, including those that resemble fibroblasts and chondrocytes, and make up the majority of cells in the atherosclerotic plaque. To better understand the epigenetic and transcriptional mechanisms that mediate these cell state changes, and how they relate to risk for coronary artery disease (CAD), we have investigated the causality and function of transcription factors at genome-wide associated loci. METHODS: We used CRISPR-Cas 9 genome and epigenome editing to identify the causal gene and cells for a complex CAD genome-wide association study signal at 2q22.3. Single-cell epigenetic and transcriptomic profiling in murine models and human coronary artery smooth muscle cells were used to understand the cellular and molecular mechanism by which this CAD risk gene exerts its function. RESULTS: CRISPR-Cas 9 genome and epigenome editing showed that the complex CAD genetic signals within a genomic region at 2q22.3 lie within smooth muscle long-distance enhancers for ZEB2, a transcription factor extensively studied in the context of epithelial mesenchymal transition in development of cancer. Zeb2 regulates SMC phenotypic transition through chromatin remodeling that obviates accessibility and disrupts both Notch and transforming growth factor ß signaling, thus altering the epigenetic trajectory of SMC transitions. SMC-specific loss of Zeb2 resulted in an inability of transitioning SMCs to turn off contractile programing and take on a fibroblast-like phenotype, but accelerated the formation of chondromyocytes, mirroring features of high-risk atherosclerotic plaques in human coronary arteries. CONCLUSIONS: These studies identify ZEB2 as a new CAD genome-wide association study gene that affects features of plaque vulnerability through direct effects on the epigenome, providing a new therapeutic approach to target vascular disease.

Embryologic Origin Influences Smooth Muscle Cell Phenotypic Modulation Signatures in Murine Marfan Syndrome Aortic Aneurysm.

BACKGROUND: Aortic root smooth muscle cells (SMC) develop from both the second heart field (SHF) and neural crest. Disparate responses to disease-causing Fbn1 variants by these lineages are proposed to promote focal aortic root aneurysm formation in Marfan syndrome (MFS), but lineage-stratified SMC analysis in vivo is lacking. METHODS: We generated SHF lineage-traced MFS mice and performed integrated multiomic (single-cell RNA and assay for transposase-accessible chromatin sequencing) analysis stratified by embryological origin. SMC subtypes were spatially identified via RNA in situ hybridization. Response to TWIST1 overexpression was determined via lentiviral transduction in human aortic SMCs. RESULTS: Lineage stratification enabled nuanced characterization of aortic root cells. We identified heightened SHF-derived SMC heterogeneity including a subset of Tnnt2 (cardiac troponin T)-expressing cells distinguished by altered proteoglycan expression. MFS aneurysm-associated SMC phenotypic modulation was identified in both SHF-traced and nontraced (neural crest-derived) SMCs; however, transcriptomic responses were distinct between lineages. SHF-derived modulated SMCs overexpressed collagen synthetic genes and small leucine-rich proteoglycans while nontraced SMCs activated chondrogenic genes. These modulated SMCs clustered focally in the aneurysmal aortic root at the region of SHF/neural crest lineage overlap. Integrated RNA-assay for transposase-accessible chromatin analysis identified enriched Twist1 and Smad2/3/4 complex binding motifs in SHF-derived modulated SMCs. TWIST1 overexpression promoted collagen and SLRP gene expression in vitro, suggesting TWIST1 may drive SHF-enriched collagen synthesis in MFS aneurysm. CONCLUSIONS: SMCs derived from both SHF and neural crest lineages undergo phenotypic modulation in MFS aneurysm but are defined by subtly distinct transcriptional responses. Enhanced TWIST1 transcription factor activity may contribute to enriched collagen synthetic pathways SHF-derived SMCs in MFS.

Genomic profiling of human vascular cells identifies TWIST1 as a causal gene for common vascular diseases.

Genome-wide association studies have identified multiple novel genomic loci associated with vascular diseases. Many of these loci are common non-coding variants that affect the expression of disease-relevant genes within coronary vascular cells. To identify such genes on a genome-wide level, we performed deep transcriptomic analysis of genotyped primary human coronary artery smooth muscle cells (HCASMCs) and coronary endothelial cells (HCAECs) from the same subjects, including splicing Quantitative Trait Loci (sQTL), allele-specific expression (ASE), and colocalization analyses. We identified sQTLs for TARS2, YAP1, CFDP1, and STAT6 in HCASMCs and HCAECs, and 233 ASE genes, a subset of which are also GTEx eGenes in arterial tissues. Colocalization of GWAS association signals for coronary artery disease (CAD), migraine, stroke and abdominal aortic aneurysm with GTEx eGenes in aorta, coronary artery and tibial artery discovered novel candidate risk genes for these diseases. At the CAD and stroke locus tagged by rs2107595 we demonstrate colocalization with expression of the proximal gene TWIST1. We show that disrupting the rs2107595 locus alters TWIST1 expression and that the risk allele has increased binding of the NOTCH signaling protein RBPJ. Finally, we provide data that TWIST1 expression influences vascular SMC phenotypes, including proliferation and calcification, as a potential mechanism supporting a role for TWIST1 in CAD.

Coronary Disease-Associated Gene TCF21 Inhibits Smooth Muscle Cell Differentiation by Blocking the Myocardin-Serum Response Factor Pathway.

RATIONALE: The gene encoding TCF21 (transcription factor 21) has been linked to coronary artery disease risk by human genome-wide association studies in multiple racial ethnic groups. In murine models, Tcf21 is required for phenotypic modulation of smooth muscle cells (SMCs) in atherosclerotic tissues and promotes a fibroblast phenotype in these cells. In humans, TCF21 expression inhibits risk for coronary artery disease. The molecular mechanism by which TCF21 regulates SMC phenotype is not known. OBJECTIVE: To better understand how TCF21 affects the SMC phenotype, we sought to investigate the possible mechanisms by which it regulates the lineage determining MYOCD (myocardin)-SRF (serum response factor) pathway. METHODS AND RESULTS: Modulation of TCF21 expression in human coronary artery SMC revealed that TCF21 suppresses a broad range of SMC markers, as well as key SMC transcription factors MYOCD and SRF, at the RNA and protein level. We conducted chromatin immunoprecipitation-sequencing to map SRF-binding sites in human coronary artery SMC, showing that binding is colocalized in the genome with TCF21, including at a novel enhancer in the SRF gene, and at the MYOCD gene promoter. In vitro genome editing indicated that the SRF enhancer CArG box regulates transcription of the SRF gene, and mutation of this conserved motif in the orthologous mouse SRF enhancer revealed decreased SRF expression in aorta and heart tissues. Direct TCF21 binding and transcriptional inhibition at colocalized sites were established by reporter gene transfection assays. Chromatin immunoprecipitation and protein coimmunoprecipitation studies provided evidence that TCF21 blocks MYOCD and SRF association by direct TCF21-MYOCD interaction. CONCLUSIONS: These data indicate that TCF21 antagonizes the MYOCD-SRF pathway through multiple mechanisms, further establishing a role for this coronary artery disease-associated gene in fundamental SMC processes and indicating the importance of smooth muscle response to vascular stress and phenotypic modulation of this cell type in coronary artery disease risk.

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.

The ESCRT-III pathway facilitates cardiomyocyte release of cBIN1-containing microparticles.

Microparticles (MPs) are cell-cell communication vesicles derived from the cell surface plasma membrane, although they are not known to originate from cardiac ventricular muscle. In ventricular cardiomyocytes, the membrane deformation protein cardiac bridging integrator 1 (cBIN1 or BIN1+13+17) creates transverse-tubule (t-tubule) membrane microfolds, which facilitate ion channel trafficking and modulate local ionic concentrations. The microfold-generated microdomains continuously reorganize, adapting in response to stress to modulate the calcium signaling apparatus. We explored the possibility that cBIN1-microfolds are externally released from cardiomyocytes. Using electron microscopy imaging with immunogold labeling, we found in mouse plasma that cBIN1 exists in membrane vesicles about 200 nm in size, which is consistent with the size of MPs. In mice with cardiac-specific heterozygous Bin1 deletion, flow cytometry identified 47% less cBIN1-MPs in plasma, supporting cardiac origin. Cardiac release was also evidenced by the detection of cBIN1-MPs in medium bathing a pure population of isolated adult mouse cardiomyocytes. In human plasma, osmotic shock increased cBIN1 detection by enzyme-linked immunosorbent assay (ELISA), and cBIN1 level decreased in humans with heart failure, a condition with reduced cardiac muscle cBIN1, both of which support cBIN1 release in MPs from human hearts. Exploring putative mechanisms of MP release, we found that the membrane fission complex endosomal sorting complexes required for transport (ESCRT)-III subunit charged multivesicular body protein 4B (CHMP4B) colocalizes and coimmunoprecipitates with cBIN1, an interaction enhanced by actin stabilization. In HeLa cells with cBIN1 overexpression, knockdown of CHMP4B reduced the release of cBIN1-MPs. Using truncation mutants, we identified that the N-terminal BAR (N-BAR) domain in cBIN1 is required for CHMP4B binding and MP release. This study links the BAR protein superfamily to the ESCRT pathway for MP biogenesis in mammalian cardiac ventricular cells, identifying elements of a pathway by which cytoplasmic cBIN1 is released into blood.

Advances in Transcriptomics: Investigating Cardiovascular Disease at Unprecedented Resolution.

Whole-genome transcriptional profiling has become a standard genomic approach to investigate biological processes. RNA sequencing (RNAseq) in particular has witnessed myriad applications in genetics and various biomedical fields. RNAseq involves a relatively simple experimental protocol of RNA extraction and cDNA library preparation and, because of decreasing next-generation sequencing cost and lower computational burden for data processing, has obtained a central role in the modern biology. The recent application of RNAseq methodology to single-cell transcriptional profiling has enabled the more precise characterization of cell lineage and cell state genetic profiles. The development of bioinformatic and statistical tools has provided for differential gene expression analysis, RNA isoform analysis, haplotype-specific analysis of gene expression (allele-specific expression), and analysis of expression quantitative trait loci. We give an overview of these and recent developments in RNAseq methodology with emphasis on quality control, read mapping, feature counting, differential gene expression, allele-specific expression and expression quantitative trait loci analysis, and fusion transcript detection. We describe utilization of RNAseq as a diagnostic tool in Mendelian diseases, complex phenotypes, and cancer and give an overview of long read RNAseq technology. Furthermore, we discuss in detail the recent revolution in single-cell transcriptomics that is reshaping modern biology.

Large-Scale Single-Cell RNA-Seq Reveals Molecular Signatures of Heterogeneous Populations of Human Induced Pluripotent Stem Cell-Derived Endothelial Cells.

RATIONALE: Human-induced pluripotent stem cell-derived endothelial cells (iPSC-ECs) have risen as a useful tool in cardiovascular research, offering a wide gamut of translational and clinical applications. However, inefficiency of the currently available iPSC-EC differentiation protocol and underlying heterogeneity of derived iPSC-ECs remain as major limitations of iPSC-EC technology. OBJECTIVE: Here, we performed droplet-based single-cell RNA sequencing (scRNA-seq) of the human iPSCs after iPSC-EC differentiation. Droplet-based scRNA-seq enables analysis of thousands of cells in parallel, allowing comprehensive analysis of transcriptional heterogeneity. METHODS AND RESULTS: Bona fide iPSC-EC cluster was identified by scRNA-seq, which expressed high levels of endothelial-specific genes. iPSC-ECs, sorted by CD144 antibody-conjugated magnetic sorting, exhibited standard endothelial morphology and function including tube formation, response to inflammatory signals, and production of NO. Nonendothelial cell populations resulting from the differentiation protocol were identified, which included immature cardiomyocytes, hepatic-like cells, and vascular smooth muscle cells. Furthermore, scRNA-seq analysis of purified iPSC-ECs revealed transcriptional heterogeneity with 4 major subpopulations, marked by robust enrichment of CLDN5, APLNR, GJA5, and ESM1 genes, respectively. CONCLUSIONS: Massively parallel, droplet-based scRNA-seq allowed meticulous analysis of thousands of human iPSCs subjected to iPSC-EC differentiation. Results showed inefficiency of the differentiation technique, which can be improved with further studies based on identification of molecular signatures that inhibit expansion of nonendothelial cell types. Subtypes of bona fide human iPSC-ECs were also identified, allowing us to sort for iPSC-ECs with specific biological function and identity.

Low prevalence of connexin-40 gene variants in atrial tissues and blood from atrial fibrillation subjects.

BACKGROUND: The atrial gap junction protein connexin-40 (Cx40) has been implicated to play an important role in atrial conduction and development of atrial fibrillation (AF). However, the frequency of Cx40 mutations in AF populations and their impact on Cx40 expression remains unclear. In this study, we sought to identify polymorphisms in the Cx40 gene GJA5, investigate the potential functional role of these polymorphisms, and determine their allelic frequencies. The prevalence of nonsynonymous Cx40 mutations in blood and atrial tissue was also compared to mutation frequencies reported in prior studies. METHODS: We conducted direct sequencing of the GJA5 coding and 3' UTR regions in blood samples from 91 lone AF subjects and 67 atrial tissue-derived samples from a lone cohort, a mixed AF cohort, and several transplant donors. Reporter gene transfection and tissue allelic expression imbalance assays were used to assess the effects of a common insertion/deletion polymorphism on Cx40 mRNA stability and expression. RESULTS: We identified one novel synonymous SNP in blood-derived DNA from a lone AF subject. In atrial tissue-derived DNA from lone and mixed AF subjects, we observed one novel nonsynonymous SNP, one rare previously reported synonymous SNP, and one novel 3' UTR SNP. A previously noted 25 bp insertion/deletion polymorphism in the 3' UTR was found to be common (minor allele frequency = 0.45) but had no effect on Cx40 mRNA stability and expression. The observed prevalence of nonsynonymous Cx40 mutations in atrial tissues derived from lone AF subjects differed significantly (p = 0.03) from a prior atrial tissue study reporting a high mutation frequency in a group of highly selected young lone AF subjects. CONCLUSIONS: Our results suggest that Cx40 coding SNPs are uncommon in AF populations, although rare mutations in this gene may certainly lead to AF pathogenesis. Furthermore, a common insertion/deletion polymorphism in the Cx40 3' UTR does not appear to play a role in modulating Cx40 mRNA levels.

Smad3 regulates smooth muscle cell fate and mediates adverse remodeling and calcification of the atherosclerotic plaque.

Atherosclerotic plaques consist mostly of smooth muscle cells (SMC), and genes that influence SMC phenotype can modulate coronary artery disease (CAD) risk. Allelic variation at 15q22.33 has been identified by genome-wide association studies to modify the risk of CAD and is associated with the expression of SMAD3 in SMC. However, the mechanism by which this gene modifies CAD risk remains poorly understood. Here we show that SMC-specific deletion of Smad3 in a murine atherosclerosis model resulted in greater plaque burden, more outward remodelling and increased vascular calcification. Single-cell transcriptomic analyses revealed that loss of Smad3 altered SMC transition cell state toward two fates: a SMC phenotype that governs both vascular remodelling and recruitment of inflammatory cells, as well as a chondromyocyte fate. Together, the findings reveal that Smad3 expression in SMC inhibits the emergence of specific SMC phenotypic transition cells that mediate adverse plaque features, including outward remodelling, monocyte recruitment, and vascular calcification.

AMPA-Type Glutamate Receptors Associated With Vascular Smooth Muscle Cell Subpopulations in Atherosclerosis and Vascular Injury.

Objectives and Aims: Vascular smooth muscle cells (VSMCs) are key constituents of both normal arteries and atherosclerotic plaques. They have an ability to adapt to changes in the local environment by undergoing phenotypic modulation. An improved understanding of the mechanisms that regulate VSMC phenotypic changes may provide insights that suggest new therapeutic targets in treatment of cardiovascular disease (CVD). The amino-acid glutamate has been associated with CVD risk and VSMCs metabolism in experimental models, and glutamate receptors regulate VSMC biology and promote pulmonary vascular remodeling. However, glutamate-signaling in human atherosclerosis has not been explored. Methods and Results: We identified glutamate receptors and glutamate metabolism-related enzymes in VSMCs from human atherosclerotic lesions, as determined by single cell RNA sequencing and microarray analysis. Expression of the receptor subunits glutamate receptor, ionotropic, α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic (AMPA)-type subunit 1 (GRIA1) and 2 (GRIA2) was restricted to cells of mesenchymal origin, primarily VSMCs, as confirmed by immunostaining. In a rat model of arterial injury and repair, changes of GRIA1 and GRIA2 mRNA level were most pronounced at time points associated with VSMC proliferation, migration, and phenotypic modulation. In vitro, human carotid artery SMCs expressed GRIA1, and selective AMPA-type receptor blocking inhibited expression of typical contractile markers and promoted pathways associated with VSMC phenotypic modulation. In our biobank of human carotid endarterectomies, low expression of AMPA-type receptor subunits was associated with higher content of inflammatory cells and a higher frequency of adverse clinical events such as stroke. Conclusion: AMPA-type glutamate receptors are expressed in VSMCs and are associated with phenotypic modulation. Patients suffering from adverse clinical events showed significantly lower mRNA level of GRIA1 and GRIA2 in their atherosclerotic lesions compared to asymptomatic patients. These results warrant further mapping of neurotransmitter signaling in the pathogenesis of human atherosclerosis.

Atheroprotective roles of smooth muscle cell phenotypic modulation and the TCF21 disease gene as revealed by single-cell analysis.

In response to various stimuli, vascular smooth muscle cells (SMCs) can de-differentiate, proliferate and migrate in a process known as phenotypic modulation. However, the phenotype of modulated SMCs in vivo during atherosclerosis and the influence of this process on coronary artery disease (CAD) risk have not been clearly established. Using single-cell RNA sequencing, we comprehensively characterized the transcriptomic phenotype of modulated SMCs in vivo in atherosclerotic lesions of both mouse and human arteries and found that these cells transform into unique fibroblast-like cells, termed 'fibromyocytes', rather than into a classical macrophage phenotype. SMC-specific knockout of TCF21-a causal CAD gene-markedly inhibited SMC phenotypic modulation in mice, leading to the presence of fewer fibromyocytes within lesions as well as within the protective fibrous cap of the lesions. Moreover, TCF21 expression was strongly associated with SMC phenotypic modulation in diseased human coronary arteries, and higher levels of TCF21 expression were associated with decreased CAD risk in human CAD-relevant tissues. These results establish a protective role for both TCF21 and SMC phenotypic modulation in this disease.

A common connexin-40 gene promoter variant affects connexin-40 expression in human atria and is associated with atrial fibrillation.

BACKGROUND: A common single-nucleotide polymorphism (SNP) in the promoter of the Connexin-40 (Cx40) gene GJA5 was suggested to affect Cx40 promoter activity and the risk of atrial fibrillation (AF), but the role of other common Cx40 polymorphisms is unknown. METHODS AND RESULTS: Eight SNPs within the Cx40 gene region were tested for association with Cx40 levels measured in atrial tissue from 61 individuals. The previously described Cx40 promoter SNP (rs35594137, -44G→A) was not associated with Cx40 mRNA levels. However, a common SNP (rs10465885) located in the TATA box of an alternative Cx40 promoter was strongly associated with Cx40 mRNA expression (P<0.0001) and displayed strong and consistent allelic expression imbalance in human atrial tissue. A promoter-luciferase assay in cultured murine cardiomyocytes demonstrated reduced activity of the promoter containing the minor allele of this SNP (P<0.0001). Both rs35594137 and rs10465885 were tested for association with early-onset lone AF (≤60 years of age) in 384 cases and 3010 population control subjects. rs10465885 was associated with the AF phenotype (odds ratio, 1.18; P=0.046). This result was confirmed in a meta-analysis including 2 additional early-onset lone AF case-control cohorts (odds ratio, 1.16, P=0.022). rs35594137 was not associated with the lone AF phenotype in any of the cohorts studied or in a combined analysis. CONCLUSIONS: A previously described Cx40 promoter SNP was not found to influence Cx40 expression or risk of AF. We describe an alternate promoter polymorphism that directly affects levels of Cx40 mRNA in vivo and is associated with early-onset lone AF.