Search

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.

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

Objective: To delineate temporal and spatial dynamics of vascular smooth muscle cell (SMC) transcriptomic changes during aortic aneurysm development in Marfan syndrome (MFS). Approach and Results: We performed single-cell RNA sequencing to study aortic root/ascending aneurysm tissue from Fbn1 C1041G/ + (MFS) mice and healthy controls, identifying all aortic cell types. A distinct cluster of transcriptomically modulated SMCs (modSMCs) was identified in adult Fbn1 C1041G/ + mouse aortic aneurysm tissue only. Comparison with atherosclerotic aortic data (ApoE −/− mice) revealed similar patterns of SMC modulation but identified an MFS-specific gene signature, including plasminogen activator inhibitor-1 ( Serpine1 ) and Kruppel-like factor 4 ( Klf4 ). We identified 481 differentially expressed genes between modSMC and SMC subsets; functional annotation highlighted extracellular matrix modulation, collagen synthesis, adhesion, and proliferation. Pseudotime trajectory analysis of Fbn1 C1041G/ + SMC/modSMC transcriptomes identified genes activated differentially throughout the course of phenotype modulation. While modSMCs were not present in young Fbn1 C1041G/ + mouse aortas despite small aortic aneurysm, multiple early modSMCs marker genes were enriched, suggesting activation of phenotype modulation. modSMCs were not found in nondilated adult Fbn1 C1041G/ + descending thoracic aortas. Single-cell RNA sequencing from human MFS aortic root aneurysm tissue confirmed analogous SMC modulation in clinical disease. Enhanced expression of TGF-β (transforming growth factor beta)-responsive genes correlated with SMC modulation in mouse and human data sets. Conclusions: Dynamic SMC phenotype modulation promotes extracellular matrix substrate modulation and aortic aneurysm progression in MFS. We characterize the disease-specific signature of modSMCs and provide temporal, transcriptomic context to the current understanding of the role TGF-β plays in MFS aortopathy. Collectively, single-cell RNA sequencing implicates TGF-β signaling and Klf4 overexpression as potential upstream drivers of SMC modulation.

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.

Atherosclerosis is the process that underlies heart attack and stroke. A characteristic feature of the atherosclerotic plaque is the accumulation of apoptotic cells in the necrotic core. Prophagocytic antibody-based therapies are currently being explored to stimulate the phagocytic clearance of apoptotic cells; however, these therapies can cause off-target clearance of healthy tissues, which leads to toxicities such as anaemia. Here we developed a macrophage-specific nanotherapy based on single-walled carbon nanotubes loaded with a chemical inhibitor of the antiphagocytic CD47-SIRPα signalling axis. We demonstrate that these single-walled carbon nanotubes accumulate within the atherosclerotic plaque, reactivate lesional phagocytosis and reduce the plaque burden in atheroprone apolipoprotein-E-deficient mice without compromising safety, and thereby overcome a key translational barrier for this class of drugs. Single-cell RNA sequencing analysis reveals that prophagocytic single-walled carbon nanotubes decrease the expression of inflammatory genes linked to cytokine and chemokine pathways in lesional macrophages, which demonstrates the potential of ‘Trojan horse’ nanoparticles to prevent atherosclerotic cardiovascular disease. Single-walled carbon nanotubes can restore phagocytotic properties of macrophages in artherosclerotic plaques to promote plaque clearance and combat artherosclerosis.

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.

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.

Summary Background Blood pressure lowering is an established strategy for preventing microvascular and macrovascular complications of diabetes, but its role in the prevention of diabetes itself is unclear. We aimed to examine this question using individual participant data from major randomised controlled trials. Methods We performed a one-stage individual participant data meta-analysis, in which data were pooled to investigate the effect of blood pressure lowering per se on the risk of new-onset type 2 diabetes. An individual participant data network meta-analysis was used to investigate the differential effects of five major classes of antihypertensive drugs on the risk of new-onset type 2 diabetes. Overall, data from 22 studies conducted between 1973 and 2008, were obtained by the Blood Pressure Lowering Treatment Trialists’ Collaboration (Oxford University, Oxford, UK). We included all primary and secondary prevention trials that used a specific class or classes of antihypertensive drugs versus placebo or other classes of blood pressure lowering medications that had at least 1000 persons-years of follow-up in each randomly allocated arm. Participants with a known diagnosis of diabetes at baseline and trials conducted in patients with prevalent diabetes were excluded. For the one-stage individual participant data meta-analysis we used stratified Cox proportional hazards model and for the individual participant data network meta-analysis we used logistic regression models to calculate the relative risk (RR) for drug class comparisons. Findings 145 939 participants (88 500 [60·6%] men and 57 429 [39·4%] women) from 19 randomised controlled trials were included in the one-stage individual participant data meta-analysis. 22 trials were included in the individual participant data network meta-analysis. After a median follow-up of 4·5 years (IQR 2·0), 9883 participants were diagnosed with new-onset type 2 diabetes. Systolic blood pressure reduction by 5 mm Hg reduced the risk of type 2 diabetes across all trials by 11% (hazard ratio 0·89 [95% CI 0·84–0·95]). Investigation of the effects of five major classes of antihypertensive drugs showed that in comparison to placebo, angiotensin-converting enzyme inhibitors (RR 0·84 [95% 0·76–0·93]) and angiotensin II receptor blockers (RR 0·84 [0·76–0·92]) reduced the risk of new-onset type 2 diabetes; however, the use of β blockers (RR 1·48 [1·27–1·72]) and thiazide diuretics (RR 1·20 [1·07–1·35]) increased this risk, and no material effect was found for calcium channel blockers (RR 1·02 [0·92–1·13]). Interpretation Blood pressure lowering is an effective strategy for the prevention of new-onset type 2 diabetes. Established pharmacological interventions, however, have qualitatively and quantitively different effects on diabetes, likely due to their differing off-target effects, with angiotensin-converting enzyme inhibitors and angiotensin II receptor blockers having the most favourable outcomes. This evidence supports the indication for selected classes of antihypertensive drugs for the prevention of diabetes, which could further refine the selection of drug choice according to an individual's clinical risk of diabetes. Funding British Heart Foundation, National Institute for Health Research, and Oxford Martin School.

Atherosclerotic plaques consist mostly of smooth muscle cells (SMC), and genes that influence SMC biology 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 expression of SMAD3 in SMC, but the mechanism by which this gene modifies CAD risk remains poorly understood. SMC-specific deletion of Smad3 in a murine atherosclerosis model resulted in greater plaque burden, positive remodeling, and increased vascular calcification. Single-cell transcriptomic analyses revealed that loss of Smad3 altered SMC progeny phenotype toward the previously described chondromyocyte fate, but importantly also promoted transition to a novel cell-state that governs remodeling and recruitment of inflammatory cells. This new remodeling population was marked by uniquely high Mmp3 and Cxcl12 expression, and its appearance correlated with higher-risk plaque features such as increased positive remodeling and macrophage content. Further, investigation of transcriptional mechanisms by which Smad3 alters SMC cell-fate revealed novel roles for Hox and Sox transcription factors whose direct interaction with Smad3 regulate an extensive transcriptional program balancing remodeling and vascular ECM with significant implications for human Mendelian aortic aneurysmal diseases. Together, these data suggest that Smad3 expression in SMC inhibits the emergence of specific SMC phenotypic transition cells that mediate adverse plaque features, including positive remodeling, monocyte recruitment, and vascular calcification.

Background To investigate the epigenetic and transcriptional mechanisms of coronary artery disease (CAD) risk, as well as the functional regulation of chromatin structure and function, we create a catalog of genetic variants associated with three stages of transcriptional cis-regulation in primary human coronary artery vascular smooth muscle cells (HCASMCs). Results We use a pooling approach with HCASMC lines to map regulatory variants that mediate binding of the CAD-associated transcription factor TCF21 with ChIPseq studies (bQTLs), variants that regulate chromatin accessibility with ATACseq studies (caQTLs), and chromosomal looping with Hi-C methods (clQTLs). We examine the overlap of these QTLs and their relationship to smooth muscle-specific genes and transcription factors. Further, we use multiple analyses to show that these QTLs are highly associated with CAD GWAS loci and correlate to lead SNPs where they show allelic effects. By utilizing genome editing, we verify that identified functional variants can regulate both chromatin accessibility and chromosomal looping, providing new insights into functional mechanisms regulating chromatin state and chromosomal structure. Finally, we directly link the disease-associated TGFB1-SMAD3 pathway to the CAD-associated FN1 gene through a response QTL that modulates both chromatin accessibility and chromosomal looping. Conclusions Together, these studies represent the most thorough mapping of multiple QTL types in a highly disease-relevant primary cultured cell type and provide novel insights into their functional overlap and mechanisms that underlie these genomic features and their relationship to disease risk.

Background: Smooth muscle cells (SMC) play a critical role in atherosclerosis. The Aryl hydrocarbon receptor (AHR) is an environment-sensing transcription factor that contributes to vascular development, and has been implicated in coronary artery disease risk. We hypothesized that AHR can affect atherosclerosis by regulating phenotypic modulation of SMC. Methods: We combined RNA-sequencing, chromatin immunoprecipitation followed by sequencing, assay for transposase-accessible chromatin using sequencing, and in vitro assays in human coronary artery SMCs, with single-cell RNA-sequencing, histology, and RNAscope in an SMC-specific lineage-tracing Ahr knockout mouse model of atherosclerosis to better understand the role of AHR in vascular disease. Results: Genomic studies coupled with functional assays in cultured human coronary artery SMCs revealed that AHR modulates the human coronary artery SMC phenotype and suppresses ossification in these cells. Lineage-tracing and activity-tracing studies in the mouse aortic sinus showed that the Ahr pathway is active in modulated SMCs in the atherosclerotic lesion cap. Furthermore, single-cell RNA-sequencing studies of the SMC-specific Ahr knockout mice showed a significant increase in the proportion of modulated SMCs expressing chondrocyte markers such as Col2a1 and Alpl , which localized to the lesion neointima. These cells, which we term “chondromyocytes,” were also identified in the neointima of human coronary arteries. In histological analyses, these changes manifested as larger lesion size, increased lineage-traced SMC participation in the lesion, decreased lineage-traced SMCs in the lesion cap, and increased alkaline phosphatase activity in lesions in the Ahr knockout in comparison with wild-type mice. We propose that AHR is likely protective based on these data and inference from human genetic analyses. Conclusions: Overall, we conclude that AHR promotes the maintenance of lesion cap integrity and diminishes the disease-related SMC-to-chondromyocyte transition in atherosclerotic tissues.

IntroductionSmooth muscle cells (SMC) play a critical role in atherosclerosis. The Aryl hydrocarbon receptor (AHR) is an environment-sensing transcription factor that contributes to vascular development, and has been implicated in coronary artery disease (CAD) risk. We hypothesized that AHR can affect atherosclerosis by regulating phenotypic modulation of SMC.MethodsWe combined RNA-Seq, ChIP-Seq, ATAC-Seq and in-vitro assays in human coronary artery SMC (HCASMC), with single-cell RNA-Seq (scRNA-Seq), histology, and RNAscope in an SMC-specific lineage-tracing Ahr knockout mouse model of atherosclerosis to better understand the role of AHR in vascular disease.ResultsGenomic studies coupled with functional assays in cultured HCASMC revealed that AHR modulates HCASMC phenotype and suppresses ossification in these cells. Lineage tracing and activity tracing studies in the mouse aortic sinus showed that the Ahr pathway is active in modulated SMC in the atherosclerotic lesion cap. Furthermore, scRNA-Seq studies of the SMC-specific Ahr knockout mice showed a significant increase in the proportion of modulated SMC expressing chondrocyte markers such as Col2a1 and Alpl, which localized to the lesion neointima. These cells, which we term “chondromyocytes” (CMC), were also identified in the neointima of human coronary arteries. In histological analyses, these changes manifested as larger lesion size, increased lineage-traced SMC participation in the lesion, decreased lineage-traced SMC in the lesion cap, and increased alkaline phosphatase activity in lesions in the Ahr knockout compared to wild-type mice. We propose that AHR is likely protective based on these data and inference from human genetic analyses.ConclusionOverall, we conclude that AHR promotes maintenance of lesion cap integrity and diminishes the disease related SMC-to-CMC transition in atherosclerotic tissues.

BackgroundTo investigate the epigenetic and transcriptional mechanisms of coronary artery disease (CAD) risk, as well as the functional regulation of chromatin structure and function, we have created a catalog of genetic variants associated with three stages of transcriptionalcis-regulation in primary human coronary artery vascular smooth muscle cells (HCASMC).ResultsTo this end, we have used a pooling approach with HCASMC lines to map regulatory variation that mediates binding of the CAD associated transcription factor TCF21 with ChIPseq studies (bQTLs), variation that regulates chromatin accessibility with ATACseq studies (caQTLs), and chromosomal looping with HiC methods (clQTLs). We show significant overlap of the QTLs, and their relationship to smooth muscle specific genes and the binding of smooth muscle transcription factors. Further, we use multiple analyses to show that these QTLs are highly associated with CAD GWAS loci and correlated to lead SNPs in these loci where they show allelic effects. We have verified with genome editing that identified functional variants can regulate both chromatin accessibility and chromosomal looping, providing new insights into functional mechanisms regulating chromatin state and chromosomal structure. Finally, we directly link the disease associatedTGFβ1-SMAD3pathway to the CAD associatedFN1gene through a response QTL that modulates both chromatin accessibility and chromosomal looping.ConclusionsTogether, these studies represent the most thorough mapping of multiple QTL types in a highly disease relevant primary cultured cell type, and provide novel insights into their functional overlap and mechanisms that underlie these genomic features and their relationship to disease risk.

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.

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., Author summary Genome-wide association studies (GWAS) have identified hundreds of genetic variants that are associated with human vascular disease including coronary artery disease. These are predominantly common single nucleotide polymorphisms (SNPs) in non-coding regions, which makes the identification of the causal genes and their underlying connection to pathophysiology challenging. Mapping of expression quantitative trait loci (eQTLs) has been performed to associate GWAS SNPs with risk genes in vascular cells and tissues. However, atherosclerotic vascular tissues contain multiple cell types. We perform deep transcriptomic profiling of genotyped human-derived vascular cells–endothelial cells and smooth muscle cells–and use splicing quantitative trait locus, allele-specific expression, and colocalization analyses to annotate genetic variants associated with vascular diseases and gain insight into their potential function in a cell-type specific manner. Based on these analyses, we identified computationally and then validate experimentally an association between the CAD risk locus rs2107595 and the gene TWIST1. We propose that the minor allele for this locus can affect transcription factor binding and provide data supporting a role for TWIST1 in modulating smooth muscle cell phenotype.

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.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.Modulation ofThese data indicate that

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.

Long noncoding RNAs (lncRNA) comprise an underlying regulatory network in the human genome that remains largely unexplored and could represent a molecular mechanism connecting GWAS risk variants with the disease causal genes. Human coronary artery smooth muscle cells (HCASMC) are one of the most important cells in the regulation of atherosclerotic disease progression. In this study, we aimed to discover the HCASMC-specific lncRNA collection and to explore the mechanistic relationship between HCAMSC lncRNAs with increased disease risk for coronary artery disease. We applied 6 different stimuli relevant to athero-condition, including stimulations with TGFbeta, TNFalpha, PDGFD and serum, and two transcription factor knock-downs TCF21 and SMAD3, and performed deep RNA sequencing on 48 HCASMC samples. We designed a lncRNA discovery pipeline to maximize the detection of tissue specific and low expressed lncRNAs. This generated a set of 53076 lncRNAs, that showed increased association with CAD GWAS variants, various HCASMC eQTL data sets and GTEx eQTLs for tissues enriched in smooth muscle. Module analysis revealed lncRNAs highly associated with TGFbeta and general pro-differentiation traits located in the CAD GWAS loci, such as, FES/FURIN, COL4A1/A2, CDKN2A/B, TGFB1, and FN1. Transcription factor motif analysis revealed the presence of HCASMC relevant factors, such as, TCF21, ZEB1, ZEB2, JUN, JUND, and an SRF cofactor ELK4, near the boundaries of HCASMC lncRNA. We defined lncRNA QTLs using the GTEx Coronary Artery dataset, and showed colocalization with other HCASMC QTLs, such as expression QTLs, chromatin looping QTLs, chromatin accessibility QTLs and TCF21 binding QTLs. We show several examples of CAD GWAS loci where lncRNAs show regulatory function, such as FES/FURIN, FN1, COL4A1/A2 and TGFB1. We unraveled a complex network of regulatory interactions at FES/FURIN locus, involving TCF21 and lncRNA 43779.9, and present a model in which TCF21 inhibited lncRNA 43779.9 regulates FES gene and indirectly the pro-differentiation FURIN gene by creating a looping interaction with the intron 1 of FES transcript. We define 5 regulatory variants at the FES/FURIN locus and propose a causal variant, rs35346340, that disrupts TCF21 binding and subsequently influences 43779.9 expression. Finally, we show the presence of lncRNAs in deeply sequenced TCF21 pooled ChIPSeq and discover TCF21 binding lncRNAs, 3938.1 and 3852.1 as ACTA2-regulating transcripts. This study defines lncRNAs as essential regulators of GWAS loci and shows the importance of using deep sequencing approach to further explore the genomic regulatory landscape of lncRNAs in human tissues.

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.

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., Author summary Coronary artery disease (CAD) is the worldwide leading cause of death. The majority of risk for CAD is genetic in nature, i.e., a feature of the genetic information that is transmitted to each individual from both parents, and primarily affects the disease processes in the blood vessel wall that regulate the disease molecular pathways. Modern genetic approaches have allowed mapping of the regions of the human genome that encode information that mediates this risk. The SMAD3 gene has been identified through these studies, a known master regulatory of other genes and molecular pathways, and we have investigated the functions of this gene that are important for disease risk. SMAD3 affects basic functions of a cellular component of the vessel wall, the smooth muscle cell (SMC), that is responsible for responding to vascular stresses to heal the lesions that are produced in conjunction with elevated lipids and other classic risk factors. Studies reported here show that SMAD3 actually inhibits the cellular processes that allow SMC to repair the vascular lesions, and its expression is promoted by the disease related variable sequences in the disease associated regions of the genome. SMAD3 is opposed by another CAD gene, TCF21, that functions to block the effects of SMAD3 expression, and these studies identify genetic mechanisms by which this is done. Thus, these studies identify an interactive pathway that directly contributes to disease risk, and the ability to block SMAD3 or promote TCF21 function could be exploited to inhibit vascular events such as myocardial infarction.

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.

Coronary artery disease (or coronary heart disease), is the leading cause of mortality in many of the developing as well as the developed countries of the world. Cholesterol-enriched plaques in the heart's blood vessels combined with inflammation lead to the lesion expansion, narrowing of blood vessels, reduced blood flow, and may subsequently cause lesion rupture and a heart attack. Even though several environmental risk factors have been established, such as high LDL-cholesterol, diabetes, and high blood pressure, the underlying genetic composition may substantially modify the disease risk; hence, genome composition and gene-environment interactions may be critical for disease progression. Ongoing scientific efforts have seen substantial advancements related to the fields of genetics and genomics, with the major breakthroughs yet to come. As genomics is the most rapidly advancing field in the life sciences, it is important to present a comprehensive overview of current efforts. Here, we present a summary of various genetic and genomics assays and approaches applied to coronary artery disease research.

Preeclampsia (PE) is a devastating complication of pregnancy that puts both mother and child at risk for clinical complications. In severe cases, PE can progress to multi-organ dysfunction and death. Characterized by hypertension (high blood pressure) and proteinuria (protein in the urine), the syndrome occurs in up to 8% of pregnancies. Although it is widely accepted that the pathophysiology of PE involves disturbance of normal blood vessel function in the placenta and the development and function of the placenta, the exact mechanisms underlying PE have not been defined. Discovery of the causal mechanisms in PE is urgently needed so that early assessment of risk and therapeutic interventions can be developed. On page [707][1] of this issue, Ho et al. ( 1 ) report breakthrough insights into the mechanisms of PE. They reveal a central role for ELABELA (ELA), a recently discovered circulating peptide ligand, in placental vascular development and PE in mice. [1]: /lookup/doi/10.1126/science.aam6607

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., Author summary Microparticles are small vesicles generated from the cell surface membrane and externally released for communication with other cells. We now show that heart ventricular muscle cells, which form the main pumping chambers of the heart, release microparticles in both mouse and human. Ventricular microparticles arise from surface membrane microdomains organized by cardiac bridging integrator 1 (cBIN1), a membrane deformation protein that has been shown to be reduced during human heart failure. Here we identify microparticles containing cBIN1 in blood, which were reduced in mutant mice lacking heart cBIN1 expression. Furthermore, the process leading to microparticle release involves the recruitment of CHMP4B protein to snip the cBIN1 membrane. In humans, cBIN1 is present in blood and within microparticles. Upon osmotic shock, human microparticles burst, allowing for the quantification of cBIN1 in plasma by enzyme-linked immunosorbent assay (ELISA). The measured cBIN1 level was greatly reduced in patients with heart failure. Thus, we introduce the biology for a new blood-based diagnostic tool that can assess cardiac muscle health and identify failing heart in human patients.

Biomarkers that predict response or toxicity to antiangiogenic therapy are sought to favorably inform the risk/benefit ratio. This study evaluated the association of vascular endothelial growth factor (VEGF) and VEGF receptor 2 (VEGFR2) genetic polymorphisms with the development of hypertension (HTN) and clinical outcome in metastatic clear cell renal cell carcinoma (MCCRCC) patients treated with sunitinib.Sixty-three MCCRCC patients receiving sunitinib (50 mg 4/2) with available blood pressure (BP) data and germline DNA were retrospectively identified. A panel of candidate VEGF and VEGFR2 single nucleotide polymorphisms (SNPs) were evaluated for associations with the development of hypertension and clinical outcome.VEGF SNP -634 genotype was associated with the prevalence and duration of sunitinib-induced hypertension (as defined by systolic pressure ≥150 mmHg and/or diastolic pressure ≥90 mmHg) in both univariable analysis (P = .03 and .01, respectively) and multivariable analysis, which adjusted for baseline BP and use of antihypertension medication (P = .05 and .02, respectively). Patients with the GG genotype were estimated to have a greater likelihood of being hypertensive during treatment compared with patients with the CC genotype (odds ratio of 13.62, 95% confidence interval [CI] 3.71-50.04). No single VEGF or VEGFR SNPs were found to correlate with clinical outcome. However, the combination of VEGF SNP 936 and VEGFR2 SNP 889 were associated with overall survival after adjustment for prognostic risk group (P = .03).In MCCRCC patients treated with sunitinib, VEGF SNP -634 is associated with hypertension and a combination of VEGF SNP 936 and VEGFR2 SNP 889 genotypes is associated with overall survival.

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

Introduction: Heritability studies of atrial fibrillation (AF) suggest a complex genetic basis for common AF. Recently a locus on chromosome 4q25 was associated with AF in a genome-wide association study (GWAS), but the mechanism by which this locus is associated with AF is unknown. Methods: To identify genetic variations associated with AF, we performed an analysis of AF as a secondary phenotype from a case-control GWAS of myocardial infarction. 210,178 autosomal SNPs of high quality were genotyped in 138 AF cases and 546 controls. Stratified score tests of the trend/additive model and general genotypic model were performed for each SNP. Genomic control and multidimensional scaling were used to control and assess population stratification. The SNP meeting genome-wide significance was genotyped in 46 left atrial appendage tissue samples for which we obtained expression levels for 22,184 transcripts. Cis expression traits within 1 megabase of the top SNP were tested for association with the SNP genotypes. Results: One SNP (rs4611994) met our per-SNP genome-wide significance thresholds (both a conservative Bonferroni threshold or a more liberal FDR-based threshold) for association with AF, which was in perfect linkage disequilibrium with the locus on chromosome 4q25 found in the deCODE genetics AF study (r 2 with SNP rs2200733=1.0). This SNP yielded an AF odds ratio of 2.28 under the additive model, with a p-value using genomic control/general genotype model of 5.0E-08, and a minor allele frequency 0.13 in controls and 0.25 in AF cases, and a genome-wide p-value of 0.046 by permutation analysis. In the 46 left atrial tissue samples, cis analysis demonstrated significant association between rs4611994 genotype and expression of LRIT3 (Leucine-rich repeat, immunoglobulin-like domain and transmembrane domain-containing protein 3) with q-value 2.22E-03 and the C isoform of PITX2 with q-value 7.50E-03. Conclusions: Our GWAS for common AF has replicated the AF-associated locus on chromosome 4q25. The genotype at this locus was correlated with cis regulation of the LRIT3 and PITX2C transcripts, suggesting these genes as candidates for mechanistic studies.

Objective Endovascular interventions cause arterial injury and induce a healing response to restore vessel wall homeostasis. Complications of defective or excessive healing are common and result in increased morbidity and repeated interventions. Experimental models of intimal hyperplasia are vital for understanding the vascular healing mechanisms and resolving the clinical problems of restenosis, vein graft stenosis, and dialysis access failure. Our aim was to systematically investigate the transcriptional, histologic, and systemic reaction to vascular injury during a prolonged time. Methods Balloon injury of the left common carotid artery was performed in male rats. Animals (n = 69) were euthanized before or after injury, either directly or after 2 hours, 20 hours, 2 days, 5 days, 2 weeks, 6 weeks, and 12 weeks. Both injured and contralateral arteries were subjected to microarray profiling, followed by bioinformatic exploration, histologic characterization of the biopsy specimens, and plasma lipid analyses. Results Immune activation and coagulation were key mechanisms in the early response, followed by cytokine release, tissue remodeling, and smooth muscle cell modulation several days after injury, with reacquisition of contractile features in later phases. Novel pathways related to clonal expansion, inflammatory transformation, and chondro-osteogenic differentiation were identified and immunolocalized to neointimal smooth muscle cells. Analysis of uninjured arteries revealed a systemic component of the reaction after local injury, underlined by altered endothelial signaling, changes in overall tissue bioenergy metabolism, and plasma high-density lipoprotein levels. Conclusions We demonstrate that vascular injury induces dynamic transcriptional landscape and metabolic changes identifiable as early, intermediate, and late response phases, reaching homeostasis after several weeks. This study provides a temporal “roadmap” of vascular healing as a publicly available resource for the research community., Clinical Relevance Endovascular intervention causes an injury to the arterial wall that subsequently induces a healing response to restore the vessel wall homeostasis. Complications after vascular interventions related to defective or excessive healing response, such as thrombosis or restenosis, are common and result in increased morbidity, suffering of the patient, need for repeated interventions, and possibly death. Thus, there is a need for better understanding of the underlying molecular mechanisms during vascular injury and healing response to identify and to assess the risk of complications in patients. Using an experimental model of vascular injury, this study demonstrates the full landscape of dynamic transcriptional changes in the resolution of vascular injury, accompanied also by systemic variations in plasma lipid levels and reaching homeostasis several weeks after injury. These results can guide the development of new strategies and molecular targets for modulation of the intimal response on endovascular interventions.