Transcriptional profile of the rat cardiovascular system at single cell resolution.

Background: Despite the critical role of the cardiovascular system, our understanding of its cellular and transcriptional diversity remains limited. We therefore sought to characterize the cellular composition, phenotypes, molecular pathways, and communication networks between cell types at the tissue and sub-tissue level across the cardiovascular system of the healthy Wistar rat, an important model in preclinical cardiovascular research. We obtained high quality tissue samples under controlled conditions that reveal a level of cellular detail so far inaccessible in human studies. Methods and Results: We performed single nucleus RNA-sequencing in 78 samples in 10 distinct regions including the four chambers of the heart, ventricular septum, sinoatrial node, atrioventricular node, aorta, pulmonary artery, and pulmonary veins (PV), which produced an aggregate map of 505,835 nuclei. We identified 26 distinct cell types and additional subtypes, including a number of rare cell types such as PV cardiomyocytes and non-myelinating Schwann cells (NMSCs), and unique groups of vascular smooth muscle cells (VSMCs), endothelial cells (ECs) and fibroblasts (FBs), which gave rise to a detailed cell type distribution across tissues. We demonstrated differences in the cellular composition across different cardiac regions and tissue-specific differences in transcription for each cell type, highlighting the molecular diversity and complex tissue architecture of the cardiovascular system. Specifically, we observed great transcriptional heterogeneities among ECs and FBs. Importantly, several cell subtypes had a unique regional localization such as a subtype of VSMCs enriched in the large vasculature. We found the cellular makeup of PV tissue is closer to heart tissue than to the large arteries. We further explored the ligand-receptor repertoire across cell clusters and tissues, and observed tissue-enriched cellular communication networks, including heightened Nppa - Npr1/2/3 signaling in the sinoatrial node. Conclusions: Through a large single nucleus sequencing effort encompassing over 500,000 nuclei, we broadened our understanding of cellular transcription in the healthy cardiovascular system. The existence of tissue-restricted cellular phenotypes suggests regional regulation of cardiovascular physiology. The overall conservation in gene expression and molecular pathways across rat and human cell types, together with our detailed transcriptional characterization of each cell type, offers the potential to identify novel therapeutic targets and improve preclinical models of cardiovascular disease.

Single-nucleus RNA sequencing in ischemic cardiomyopathy reveals common transcriptional profile underlying end-stage heart failure.

Ischemic cardiomyopathy (ICM) is the leading cause of heart failure worldwide, yet the cellular and molecular signature of this disease is largely unclear. Using single-nucleus RNA sequencing (snRNA-seq) and integrated computational analyses, we profile the transcriptomes of over 99,000 human cardiac nuclei from the non-infarct region of the left ventricle of 7 ICM transplant recipients and 8 non-failing (NF) controls. We find the cellular composition of the ischemic heart is significantly altered, with decreased cardiomyocytes and increased proportions of lymphatic, angiogenic, and arterial endothelial cells in patients with ICM. We show that there is increased LAMININ signaling from endothelial cells to other cell types in ICM compared with NF. Finally, we find that the transcriptional changes that occur in ICM are similar to those in hypertrophic and dilated cardiomyopathies and that the mining of these combined datasets can identify druggable genes that could be used to target end-stage heart failure.

Aortic Cellular Diversity and Quantitative Genome-Wide Association Study Trait Prioritization Through Single-Nuclear RNA Sequencing of the Aneurysmal Human Aorta.

BACKGROUND: Mural cells in ascending aortic aneurysms undergo phenotypic changes that promote extracellular matrix destruction and structural weakening. To explore this biology, we analyzed the transcriptional features of thoracic aortic tissue. METHODS: Single-nuclear RNA sequencing was performed on 13 samples from human donors, 6 with thoracic aortic aneurysm, and 7 without aneurysm. Individual transcriptomes were then clustered based on transcriptional profiles. Clusters were used for between-disease differential gene expression analyses, subcluster analysis, and analyzed for intersection with genetic aortic trait data. RESULTS: We sequenced 71 689 nuclei from human thoracic aortas and identified 14 clusters, aligning with 11 cell types, predominantly vascular smooth muscle cells (VSMCs) consistent with aortic histology. With unbiased methodology, we found 7 vascular smooth muscle cell and 6 fibroblast subclusters. Differentially expressed genes analysis revealed a vascular smooth muscle cell group accounting for the majority of differential gene expression. Fibroblast populations in aneurysm exhibit distinct behavior with almost complete disappearance of quiescent fibroblasts. Differentially expressed genes were used to prioritize genes at aortic diameter and distensibility genome-wide association study loci highlighting the genes JUN, LTBP4 (latent transforming growth factor beta-binding protein 1), and IL34 (interleukin 34) in fibroblasts, ENTPD1, PDLIM5 (PDZ and LIM domain 5), ACTN4 (alpha-actinin-4), and GLRX in vascular smooth muscle cells, as well as LRP1 in macrophage populations. CONCLUSIONS: Using nuclear RNA sequencing, we describe the cellular diversity of healthy and aneurysmal human ascending aorta. Sporadic aortic aneurysm is characterized by differential gene expression within known cellular classes rather than by the appearance of novel cellular forms. Single-nuclear RNA sequencing of aortic tissue can be used to prioritize genes at aortic trait loci.

Braveheart, a long noncoding RNA required for cardiovascular lineage commitment.

Long noncoding RNAs (lncRNAs) are often expressed in a development-specific manner, yet little is known about their roles in lineage commitment. Here, we identified Braveheart (Bvht), a heart-associated lncRNA in mouse. Using multiple embryonic stem cell (ESC) differentiation strategies, we show that Bvht is required for progression of nascent mesoderm toward a cardiac fate. We find that Bvht is necessary for activation of a core cardiovascular gene network and functions upstream of mesoderm posterior 1 (MesP1), a master regulator of a common multipotent cardiovascular progenitor. We also show that Bvht interacts with SUZ12, a component of polycomb-repressive complex 2 (PRC2), during cardiomyocyte differentiation, suggesting that Bvht mediates epigenetic regulation of cardiac commitment. Finally, we demonstrate a role for Bvht in maintaining cardiac fate in neonatal cardiomyocytes. Together, our work provides evidence for a long noncoding RNA with critical roles in the establishment of the cardiovascular lineage during mammalian development.

The Drosophila HP1 homolog Rhino is required for transposon silencing and piRNA production by dual-strand clusters.

Piwi-interacting RNAs (piRNAs) silence transposons and maintain genome integrity during germline development. In Drosophila, transposon-rich heterochromatic clusters encode piRNAs either on both genomic strands (dual-strand clusters) or predominantly one genomic strand (uni-strand clusters). Primary piRNAs derived from these clusters are proposed to drive a ping-pong amplification cycle catalyzed by proteins that localize to the perinuclear nuage. We show that the HP1 homolog Rhino is required for nuage organization, transposon silencing, and ping-pong amplification of piRNAs. rhi mutations virtually eliminate piRNAs from the dual-strand clusters and block production of putative precursor RNAs from both strands of the major 42AB dual-strand cluster, but not of transcripts or piRNAs from the uni-strand clusters. Furthermore, Rhino protein associates with the 42AB dual-strand cluster,but does not bind to uni-strand cluster 2 or flamenco. Rhino thus appears to promote transcription of dual-strand clusters, leading to production of piRNAs that drive the ping-pong amplification cycle.

Screening for novel regulators of embryonic stem cell identity.

Two recent studies, including one in this issue of Cell Stem Cell, have identified novel regulators of embryonic stem cell (ESC) self-renewal using genome-wide RNAi screens in mouse ESCs (Ding et al., 2009; Hu et al., 2009) and have further expanded the repertoire of factors that regulate ESC identity.

Collapse of germline piRNAs in the absence of Argonaute3 reveals somatic piRNAs in flies.

Piwi-interacting RNAs (piRNAs) silence transposons in animal germ cells. piRNAs are thought to derive from long transcripts spanning transposon-rich genomic loci and to direct an autoamplification loop in which an antisense piRNA, bound to Aubergine or Piwi protein, triggers production of a sense piRNA bound to the PIWI protein Argonaute3 (Ago3). In turn, the new piRNA is envisioned to produce a second antisense piRNA. Here, we describe strong loss-of-function mutations in ago3, allowing a direct genetic test of this model. We find that Ago3 acts to amplify piRNA pools and to enforce on them an antisense bias, increasing the number of piRNAs that can act to silence transposons. We also detect a second, Ago3-independent piRNA pathway centered on Piwi. Transposons targeted by this second pathway often reside in the flamenco locus, which is expressed in somatic ovarian follicle cells, suggesting a role for piRNAs beyond the germline.

Biogenesis and germline functions of piRNAs.

Small interfering RNAs and microRNAs are generated from double-stranded RNA precursors by the Dicer endonucleases, and function with Argonaute-family proteins to target transcript destruction or to silence translation. A distinct class of 24- to 30-nucleotide-long RNAs, produced by a Dicer-independent mechanism, associates with Piwi-class Argonaute proteins. Studies in flies, fish and mice implicate these Piwi-associated RNAs (piRNAs) in germline development, silencing of selfish DNA elements, and in maintaining germline DNA integrity. However, whether piRNAs primarily control chromatin organization, gene transcription, RNA stability or RNA translation is not well understood, neither is piRNA biogenesis. Here, we review recent studies of piRNA production and function, and discuss unanswered questions about this intriguing new class of small RNAs.

Drosophila rasiRNA pathway mutations disrupt embryonic axis specification through activation of an ATR/Chk2 DNA damage response.

Small repeat-associated siRNAs (rasiRNAs) mediate silencing of retrotransposons and the Stellate locus. Mutations in the Drosophila rasiRNA pathway genes armitage and aubergine disrupt embryonic axis specification, triggering defects in microtubule polarization as well as asymmetric localization of mRNA and protein determinants in the developing oocyte. Mutations in the ATR/Chk2 DNA damage signal transduction pathway dramatically suppress these axis specification defects, but do not restore retrotransposon or Stellate silencing. Furthermore, rasiRNA pathway mutations lead to germline-specific accumulation of gamma-H2Av foci characteristic of DNA damage. We conclude that rasiRNA-based gene silencing is not required for axis specification, and that the critical developmental function for this pathway is to suppress DNA damage signaling in the germline.

Normal microRNA maturation and germ-line stem cell maintenance requires Loquacious, a double-stranded RNA-binding domain protein.

microRNAs (miRNAs) are single-stranded, 21- to 23-nucleotide cellular RNAs that control the expression of cognate target genes. Primary miRNA (pri-miRNA) transcripts are transformed to mature miRNA by the successive actions of two RNase III endonucleases. Drosha converts pri-miRNA transcripts to precursor miRNA (pre-miRNA); Dicer, in turn, converts pre-miRNA to mature miRNA. Here, we show that normal processing of Drosophila pre-miRNAs by Dicer-1 requires the double-stranded RNA-binding domain (dsRBD) protein Loquacious (Loqs), a homolog of human TRBP, a protein first identified as binding the HIV trans-activator RNA (TAR). Efficient miRNA-directed silencing of a reporter transgene, complete repression of white by a dsRNA trigger, and silencing of the endogenous Stellate locus by Suppressor of Stellate, all require Loqs. In loqs(f00791) mutant ovaries, germ-line stem cells are not appropriately maintained. Loqs associates with Dcr-1, the Drosophila RNase III enzyme that processes pre-miRNA into mature miRNA. Thus, every known Drosophila RNase-III endonuclease is paired with a dsRBD protein that facilitates its function in small RNA biogenesis.

Cholesterol depletion induces PKA-mediated basolateral-to-apical transcytosis of the scavenger receptor class B type I in MDCK cells.

Cholesterol-based membrane microdomains, or lipid rafts, are believed to play important, yet poorly defined, roles in protein trafficking and signal transduction. In polarized epithelial cells, the current view is that rafts are involved in apical but not in basolateral protein transport from the trans-Golgi network (TGN). We report here that cholesterol is required in a post-TGN mechanism of basolateral regionalization. Permanently transfected Madin-Darby canine kidney cells segregated the caveolae/raft-associated high-density lipoprotein scavenger receptor class B type I (SR-BI) predominantly to the basolateral domain where it was constitutively internalized and recycled basolaterally. Acute cholesterol depletion did not significantly alter SR-BI internalization, implying a cholesterol depletion-insensitive endocytic process but instead induced its transcytosis through a protein kinase A (PKA)- and microtubule-dependent mechanism. Forskolin also elicited SR-BI transcytosis. The basolateral distribution of endogenous epidermal growth factor receptor remained unaffected. Strikingly, cholesterol depletion induced PKA activity without increasing the cAMP levels. Thus, our results are consistent with a scenario in which cholesterol-based rafts promote internalization and basolateral recycling of internalized SR-BI whereas a PKA pool sensitive to cholesterol depletion mediates SR-BI transcytosis. Regulated transcytosis of SR-BI may provide an additional mechanism to control cholesterol homeostasis. These results disclose relationships between cholesterol-based rafts and PKA activity operating in a post-TGN mechanism of regulated apical-to-basolateral cell surface protein distribution.

Human brain synembryn interacts with Gsalpha and Gqalpha and is translocated to the plasma membrane in response to isoproterenol and carbachol.

Heterotrimeric G-proteins transduce signals from heptahelical transmembrane receptors to different effector systems, regulating diverse complex intracellular pathways and functions. In brain, facilitation of depolarization-induced neurotransmitter release for synaptic transmission is mediated by Gsalpha and Gqalpha. To identify effectors for Galpha-proteins, we performed a yeast two-hybrid screening of a human brain cDNA library, using the human Galphas protein as a bait. We identified a protein member of the synembryn family as one of the interacting proteins. Extending the study to other Galpha subunits, we found that Gqalpha also interacts with synembryn, and these interactions were confirmed by in vitro pull down studies and by in vivo confocal laser microscopy analysis. Furthermore, synembryn was shown to translocate to the plasma membrane in response to carbachol and isoproterenol. This study supports recent findings in C. elegans where, through genetic studies, synembryn was shown to act together with Gqalpha regulating neuronal transmitter release. Based on these observations, we propose that synembryn is playing a similar role in human neuronal cells.