Discovery of Transacting Long Noncoding RNAs That Regulate Smooth Muscle Cell Phenotype.

BACKGROUND: Smooth muscle cells (SMC), the major cell type in atherosclerotic plaques, are vital in coronary artery diseases (CADs). SMC phenotypic transition, which leads to the formation of various cell types in atherosclerotic plaques, is regulated by a network of genetic and epigenetic mechanisms and governs the risk of disease. The involvement of long noncoding RNAs (lncRNAs) has been increasingly identified in cardiovascular disease. However, SMC lncRNAs have not been comprehensively characterized, and their regulatory role in SMC state transition remains unknown. METHODS: A discovery pipeline was constructed and applied to deeply strand-specific RNA sequencing from perturbed human coronary artery SMC with different disease-related stimuli, to allow for the detection of novel lncRNAs. The functional relevance of a select few novel lncRNAs were verified in vitro. RESULTS: We identified 4579 known and 13 655 de novo lncRNAs in human coronary artery SMC. Consistent with previous long noncoding RNA studies, these lncRNAs overall have fewer exons, are shorter in length than protein-coding genes (pcGenes), and have relatively low expression level. Genomic location of these long noncoding RNA is disproportionately enriched near CAD-related TFs (transcription factors), genetic loci, and gene regulators of SMC identity, suggesting the importance of their function in disease. Two de novo lncRNAs, ZIPPOR (ZEB-interacting suppressor) and TNS1-AS2 (TNS1-antisense 2), were identified by our screen. Combining transcriptional data and in silico modeling along with in vitro validation, we identified CAD gene ZEB2 as a target through which these lncRNAs exert their function in SMC phenotypic transition. CONCLUSIONS: Expression of a large and diverse set of lncRNAs in human coronary artery SMC are highly dynamic in response to CAD-related stimuli. The dynamic changes in expression of these lncRNAs correspond to alterations in transcriptional programs that are relevant to CAD, suggesting a critical role for lncRNAs in SMC phenotypic transition and human atherosclerotic disease.

Comprehensive Integration of Multiple Single-Cell Transcriptomic Data Sets Defines Distinct Cell Populations and Their Phenotypic Changes in Murine Atherosclerosis.

BACKGROUND: The application of single-cell transcriptomic (single-cell RNA sequencing) analysis to the study of atherosclerosis has provided unique insights into the molecular and genetic mechanisms that mediate disease risk and pathophysiology. However, nonstandardized methodologies and relatively high costs associated with the technique have limited the size and replication of existing data sets and created disparate or contradictory findings that have fostered misunderstanding and controversy. METHODS: To address these uncertainties, we have performed a conservative integration of multiple published single-cell RNA sequencing data sets into a single meta-analysis, performed extended analysis of native resident vascular cells, and used in situ hybridization to map the disease anatomic location of the identified cluster cells. To investigate the transdifferentiation of smooth muscle cells to macrophage phenotype, we have developed a classifying algorithm based on the quantification of reporter transgene expression. RESULTS: The reporter gene expression tool indicates that within the experimental limits of the examined studies, transdifferentiation of smooth muscle cell to the macrophage lineage is extremely rare. Validated transition smooth muscle cell phenotypes were defined by clustering, and the location of these cells was mapped to lesion anatomy with in situ hybridization. We have also characterized 5 endothelial cell phenotypes and linked these cellular species to different vascular structures and functions. Finally, we have identified a transcriptomically unique cellular phenotype that constitutes the aortic valve. CONCLUSIONS: Taken together, these analyses resolve a number of outstanding issues related to differing results reported with vascular disease single-cell RNA sequencing studies, and significantly extend our understanding of the role of resident vascular cells in anatomy and disease.

Autophagy Is Differentially Regulated in Leukocyte and Nonleukocyte Foam Cells During Atherosclerosis.

RATIONALE: Atherosclerosis is characterized by an accumulation of foam cells within the arterial wall, resulting from excess cholesterol uptake and buildup of cytosolic lipid droplets (LDs). Autophagy promotes LD clearance by freeing stored cholesterol for efflux, a process that has been shown to be atheroprotective. While the role of autophagy in LD catabolism has been studied in macrophage-derived foam cells, this has remained unexplored in vascular smooth muscle cell (VSMC)-derived foam cells that constitute a large fraction of foam cells within atherosclerotic lesions. OBJECTIVE: We performed a comparative analysis of autophagy flux in lipid-rich aortic intimal populations to determine whether VSMC-derived foam cells metabolize LDs similarly to their macrophage counterparts. METHODS AND RESULTS: Atherosclerosis was induced in GFP-LC3 (microtubule-associated proteins 1A/1B light chain 3) transgenic mice by PCSK9 (proprotein convertase subtilisin/kexin type 9)-adeno-associated viral injection and Western diet feeding. Using flow cytometry of aortic digests, we observed a significant increase in dysfunctional autophagy of VSMC-derived foam cells during atherogenesis relative to macrophage-derived foam cells. Using cell culture models of lipid-loaded VSMCs and macrophages, we show that autophagy-mediated cholesterol efflux from VSMC foam cells was poor relative to macrophage foam cells, and largely occurs when HDL (high-density lipoprotein) was used as a cholesterol acceptor, as opposed to apoA-1 (apolipoproteinA-1). This was associated with the predominant expression of ABCG1 in VSMC foam cells. Using metformin, an autophagy activator, cholesterol efflux to HDL was significantly increased in VSMC, but not in macrophage, foam cells. CONCLUSIONS: These data demonstrate that VSMC and macrophage foam cells perform cholesterol efflux by distinct mechanisms, and that autophagy flux is highly impaired in VSMC foam cells, but can be induced by pharmacological means. Further investigation is warranted into targeting autophagy specifically in VSMC foam cells, the predominant foam cell subtype of advanced atherosclerotic plaques, to promote reverse cholesterol transport and resolution of the atherosclerotic plaque.

Human Coronary Plaque T Cells Are Clonal and Cross-React to Virus and Self.

BACKGROUND: Coronary artery disease is an incurable, life-threatening disease that was once considered primarily a disorder of lipid deposition. Coronary artery disease is now also characterized by chronic inflammation' notable for the buildup of atherosclerotic plaques containing immune cells in various states of activation and differentiation. Understanding how these immune cells contribute to disease progression may lead to the development of novel therapeutic strategies. METHODS: We used single-cell technology and in vitro assays to interrogate the immune microenvironment of human coronary atherosclerotic plaque at different stages of maturity. RESULTS: In addition to macrophages, we found a high proportion of αß T cells in the coronary plaques. Most of these T cells lack high expression of CCR7 and L-selectin, indicating that they are primarily antigen-experienced memory cells. Notably, nearly one-third of these cells express the HLA-DRA surface marker, signifying activation through their TCRs (T-cell receptors). Consistent with this, TCR repertoire analysis confirmed the presence of activated αß T cells (CD4

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.

Environment-Sensing Aryl Hydrocarbon Receptor Inhibits the Chondrogenic Fate of Modulated Smooth Muscle Cells in Atherosclerotic Lesions.

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.

Single-Cell Transcriptomic Profiling of Vascular Smooth Muscle Cell Phenotype Modulation in Marfan Syndrome Aortic Aneurysm.

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 Fbn1C1041G/+ (MFS) mice and healthy controls, identifying all aortic cell types. A distinct cluster of transcriptomically modulated SMCs (modSMCs) was identified in adult Fbn1C1041G/+ 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 Fbn1C1041G/+ SMC/modSMC transcriptomes identified genes activated differentially throughout the course of phenotype modulation. While modSMCs were not present in young Fbn1C1041G/+ 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 Fbn1C1041G/+ 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.

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.

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.

Genetics and Genomics of Coronary Artery Disease.

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.

Local tissue mechanics control cardiac pacemaker cell embryonic patterning.

Cardiac pacemaker cells (CPCs) initiate the electric impulses that drive the rhythmic beating of the heart. CPCs reside in a heterogeneous, ECM-rich microenvironment termed the sinoatrial node (SAN). Surprisingly, little is known regarding the biochemical composition or mechanical properties of the SAN, and how the unique structural characteristics present in this region of the heart influence CPC function remains poorly understood. Here, we have identified that SAN development involves the construction of a "soft" macromolecular ECM that specifically encapsulates CPCs. In addition, we demonstrate that subjecting embryonic CPCs to substrate stiffnesses higher than those measured in vivo results in loss of coherent electrical oscillation and dysregulation of the HCN4 and NCX1 ion channels required for CPC automaticity. Collectively, these data indicate that local mechanics play a critical role in maintaining the embryonic CPC function while also quantitatively defining the range of material properties that are optimal for embryonic CPC maturation.

Association of VEGF and VEGFR2 single nucleotide polymorphisms with hypertension and clinical outcome in metastatic clear cell renal cell carcinoma patients treated with sunitinib.

PURPOSE: 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. PATIENT AND METHODS: 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. RESULTS: 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). CONCLUSIONS: 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.

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.