Search

Mutations in the spliceosomal RNA binding protein RBM10 cause TARP syndrome and are frequently observed in lung adenocarcinoma (LUAD). We have previously shown that RBM10 enhances exon skipping of its target genes, including its paralog RBM5. Here, we report that RBM10 negatively regulates its own mRNA and protein expression and that of RBM5 by promoting alternative splicing-coupled nonsense-mediated mRNA decay (AS-NMD). Through computational analysis and experimental validation, we identified RBM10-promoted skipping of exon 6 or 12 in RBM10 and exon 6 or 16 in RBM5 as the underlying AS-NMD events. Importantly, we showed that LUAD-associated mutations affecting splice sites of RBM10 exons 6 or 12 abolished exon inclusion and correlated with reduced expression of RBM10 RNA. Together, our investigations have revealed novel molecular mechanisms underlying RBM10 autoregulation and cross-regulation of RBM5, thereby providing insights concerning the functions of RBM10 under various physiological and pathological conditions. Our combined computational and experimental approach should be useful for elucidating the role of AS-NMD in auto- and cross-regulation by other splicing regulators.

Objective— Cellular metabolism was recently shown to regulate endothelial cell phenotype profoundly. Whether the atheroprotective biomechanical stimulus elicited by laminar shear stress modulates endothelial cell metabolism is not known. Approach and Results— Here, we show that laminar flow exposure reduced glucose uptake and mitochondrial content in endothelium. Shear stress–mediated reduction of endothelial metabolism was reversed by silencing the flow-sensitive transcription factor Krüppel-like factor 2 (KLF2). Endothelial-specific deletion of KLF2 in mice induced glucose uptake in endothelial cells of perfused hearts. KLF2 overexpression recapitulates the inhibitory effects on endothelial glycolysis elicited by laminar flow, as measured by Seahorse flux analysis and glucose uptake measurements. RNA sequencing showed that shear stress reduced the expression of key glycolytic enzymes, such as 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase-3 (PFKFB3), phosphofructokinase-1, and hexokinase 2 in a KLF2-dependent manner. Moreover, KLF2 represses PFKFB3 promoter activity. PFKFB3 knockdown reduced glycolysis, and overexpression increased glycolysis and partially reversed the KLF2-mediated reduction in glycolysis. Furthermore, PFKFB3 overexpression reversed KLF2-mediated reduction in angiogenic sprouting and network formation. Conclusions— Our data demonstrate that shear stress–mediated repression of endothelial cell metabolism via KLF2 and PFKFB3 controls endothelial cell phenotype.

Rationale: The human genome harbors a large number of sequences encoding for RNAs that are not translated but control cellular functions by distinct mechanisms. The expression and function of the longer transcripts namely the long noncoding RNAs in the vasculature are largely unknown. Objective: Here, we characterized the expression of long noncoding RNAs in human endothelial cells and elucidated the function of the highly expressed metastasis-associated lung adenocarcinoma transcript 1 (MALAT1). Methods and Results: Endothelial cells of different origin express relative high levels of the conserved long noncoding RNAs MALAT1, taurine upregulated gene 1 (TUG1), maternally expressed 3 (MEG3), linc00657, and linc00493. MALAT1 was significantly increased by hypoxia and controls a phenotypic switch in endothelial cells. Silencing of MALAT1 by small interfering RNAs or GapmeRs induced a promigratory response and increased basal sprouting and migration, whereas proliferation of endothelial cells was inhibited. When angiogenesis was further stimulated by vascular endothelial growth factor, MALAT1 small interfering RNAs induced discontinuous sprouts indicative of defective proliferation of stalk cells. In vivo studies confirmed that genetic ablation of MALAT1 inhibited proliferation of endothelial cells and reduced neonatal retina vascularization. Pharmacological inhibition of MALAT1 by GapmeRs reduced blood flow recovery and capillary density after hindlimb ischemia. Gene expression profiling followed by confirmatory quantitative reverse transcriptase-polymerase chain reaction demonstrated that silencing of MALAT1 impaired the expression of various cell cycle regulators. Conclusions: Silencing of MALAT1 tips the balance from a proliferative to a migratory endothelial cell phenotype in vitro, and its genetic deletion or pharmacological inhibition reduces vascular growth in vivo.

The Drosophila melanogaster gene Dscam (Down syndrome cell adhesion molecule) can generate thousands of different ectodomains via mutual exclusive splicing of three large exon clusters. The isoform diversity plays a profound role in both neuronal wiring and pathogen recognition. However, the isoform expression pattern at the global level remained unexplored. Here, we developed a novel method that allows for direct quantification of the alternatively spliced exon combinations from over hundreds of millions of Dscam transcripts in one sequencing run. With unprecedented sequencing depth, we detected a total of 18 496 isoforms, out of 19 008 theoretically possible combinations. Importantly, we demonstrated that alternative splicing between different clusters is independent. Moreover, the isoforms were expressed across a broad dynamic range, with significant bias in cell/tissue and developmental stage-specific patterns. Hitherto underappreciated, such bias can dramatically reduce the ability of neurons to display unique surface receptor codes. Therefore, the seemingly excessive diversity encoded in the Dscam locus might nevertheless be essential for a robust self and non-self discrimination in neurons.

Argonaute2 (Ago2) is an established component of the microRNA-induced silencing complex. Similar to miR-375 loss-of-function studies, inhibition of Ago2 in the pancreatic β-cell resulted in enhanced insulin release underlining the relationship between these two genes. Moreover, as the most abundant microRNA in pancreatic endocrine cells, miR-375 was also observed to be enriched in Ago2-associated complexes. Both Ago2 and miR-375 regulate the pancreatic β-cell secretome, and by using quantitative mass spectrometry, we identified the enhanced release of a set of proteins or secretion "signatures " in response to a glucose stimulus using the murine β-cell line MIN6. In addition, the loss of Ago2 resulted in the increased expression of miR-375 target genes, including gephyrin and ywhaz. These targets positively contribute to exocytosis indicating they may mediate the functional role of both miR-375 and Ago proteins in the pancreatic β-cell by influencing the secretory pathway. This study specifically addresses the role of Ago2 in the systemic release of proteins from β-cells and highlights the contribution of the microRNA pathway to the function of this cell type.

MicroRNAs (miRNAs) constitute an important class of small regulatory RNAs that are derived from distinct hairpin precursors (pre-miRNAs). In contrast to mature miRNAs, which have been characterized in numerous genome-wide studies of different organisms, research on global profiling of pre-miRNAs is limited. Here, using massive parallel sequencing, we have performed global characterization of both mouse mature and precursor miRNAs. In total, 87 369 704 and 252 003 sequencing reads derived from 887 mature and 281 precursor miRNAs were obtained, respectively. Our analysis revealed new aspects of miRNA/pre-miRNA processing and modification, including eight Ago2-cleaved pre-miRNAs, eight new instances of miRNA editing and exclusively 5' tailed mirtrons. Furthermore, based on the sequences of both mature and precursor miRNAs, we developed a miRNA discovery pipeline, miRGrep, which does not rely on the availability of genome reference sequences. In addition to 239 known mouse pre-miRNAs, miRGrep predicted 41 novel ones with high confidence. Similar as known ones, the mature miRNAs derived from most of these novel loci showed both reduced abundance following Dicer knockdown and the binding with Argonaute2. Evaluation on data sets obtained from Caenorhabditis elegans and Caenorhabditis sp.11 demonstrated that miRGrep could be widely used for miRNA discovery in metazoans, especially in those without genome reference sequences.

Background: Adenosine to inosine (A-to-I) RNA editing is catalysed by ADARs (adenosine deaminases acting on RNA) and is an important posttranscriptional regulator of RNA metabolism. Although RNA editing is essential for life, its role in cardiovascular disease is unknown. Methods and Results: RNA-sequencing (RNA-seq) of human endothelial cells revealed that ADAR1 is the main RNA editor in HUVECs. The vast majority of editing events are detected in Alu regions due to their ability to form double-stranded structures, a prerequisite for RNA editing. A-to-I RNA editing of the primate-specific Alu elements is observed in 25% of poly(A)+ transcripts and in 44% of ribosomal RNA-depleted transcripts. Cathepsin S (CTSS), an extracellular matrix degrading enzyme, has 3 Alu elements in its 3’-untranslated region and is among the most extensively edited mRNAs. Interestingly, ADAR1 silencing downregulates CTSS mRNA expression, whereas ADAR1 overexpression results in an increase of both CTSS RNA editing rate and CTSS mRNA expression. Mechanistically, RNA immunoprecipitation (RIP), luciferase reporter and iCLIP assays showed that A-to-I RNA editing of the Alu elements in the 3’UTR of CTSS mRNA regulates the recruitment of the stabilizing RNA-binding protein HuR to CTSS 3’UTR and, thus, controls CTSS mRNA stability and expression. CTSS Alu RNA editing is increased under hypoxic or pro-inflammatory conditions in vitro , as well as in patients with coronary or carotid atherosclerotic vascular disease. Further, the extent of CTSS Alu RNA editing rate correlated with various markers of subclinical atherosclerosis including intima-media thickness and number of atherosclerotic plaques. Importantly, ADAR1 and the extent of CTSS RNA editing correlated with CTSS expression levels in patients’ samples from 6 non-overlapping cohorts of patients with inflammatory vascular diseases (total n=366), including peripheral blood mononuclear cells, carotid atherosclerotic plaques and thoracic aortic aneurysms. Conclusion: This study shows for the first time that Alu A-to-I RNA editing is a critical modulator of inflammatory gene expression in all stages of atherosclerotic vascular disease development.

Background: Adenosine to inosine (A-to-I) RNA editing is catalysed by ADARs (adenosine deaminases acting on RNA) and is an important posttranscriptional regulator of RNA metabolism. Although RNA editing is essential for life, its role in cardiovascular disease is unknown. Methods and Results: Next generation RNA-sequencing (RNA-seq) of human endothelial cells revealed that ADAR1 edits almost 25% of endothelial transcriptome, regulates gene expression (135 transcripts were at least 2-fold regulated after ADAR1 knockdown) and controls endothelial cell function. Among the highest edited transcripts is cathepsin S (CTSS), an extracellular matrix degradation enzyme with an established role in cardiovascular and inflammatory diseases. The number of edited nucleotides in the 3’UTR of CTSS mRNA was increased under hypoxic conditions and was strongly associated with increased CTSS expression (r=0.971, P=0.029). Most importantly, absence of RNA editing resulted in a 3-fold decrease of CTSS mRNA expression under hypoxic conditions. RNA-seq analysis of 8 patients with coronary heart disease (CHD) and 4 healthy subjects revealed that A-to-I RNA editing of CTSS 3’UTR is significantly increased in patients with CHD compared to control healthy subjects (P0.7, P Conclusion: This study suggests a novel primate-specific mechanism that may account for the increased cathepsin S expression in ischemic/inflammatory diseases and it may serve as a prototypical example for the evaluation of RNA-based mechanisms in health and disease.

Circular RNAs (circRNAs) have re-emerged as an interesting RNA species. Here, by deep RNA profiling in different mouse tissues, we observed that circRNAs are significantly enriched in brain.and a disproportionate fraction of them is derived from host genes that code for synaptic proteins. Moreover, based on separate profiling of the RNAs localized in neuronal cell bodies and neuropil, on average, circRNAs are more enriched in the neuropil than their host gene mRNA isoforms. Using high resolution in situ hybridization we, for the first time, visualized circRNA punctae in the dendrites of neurons. Consistent with the idea that circRNAs might regulate synaptic function, during development, many circRNAs change their abundance abruptly at a time corresponding to synaptogenesis. In addition, following a homeostatic downscaling of neuronal activity many circRNAs exhibit significant up or down-regulation. Together, our data indicate that brain circRNAs are positioned to respond to and regulate synaptic function.

Background Exome sequencing is increasingly used to search for phenotypically-relevant sequence variants in the mouse genome. All of the current hybridization-based mouse exome capture systems are designed based on the genome reference sequences of the C57BL/6 J strain. Given that the substantial sequence divergence exists between C57BL/6 J and other distantly-related strains, the impact of sequence divergence on the efficiency of such capture systems needs to be systematically evaluated before they can be widely applied to the study of those strains. Results Using the Agilent SureSelect mouse exome capture system, we performed exome sequencing on F1 generation hybrid mice that were derived by crossing two divergent strains, C57BL/6 J and SPRET/EiJ. Our results showed that the C57BL/6 J-based probes captured the sequences derived from C57BL/6 J alleles more efficiently and that the bias was higher for the target regions with greater sequence divergence. At low sequencing depths, the bias also affected the efficiency of variant detection. However, the effects became negligible when sufficient sequencing depth was achieved. Conclusion Sufficient sequence depth needs to be planned to match the sequence divergence between C57BL/6 J and the strain to be studied, when the C57BL/6 J–based Agilent SureSelect exome capture system is to be used.

Freshwater planaria are a very attractive model system for stem cell biology, tissue homeostasis, and regeneration. The genome of the planarian Schmidtea mediterranea has recently been sequenced and is estimated to contain >20,000 protein-encoding genes. However, the characterization of its transcriptome is far from complete. Furthermore, not a single proteome of the entire phylum has been assayed on a genome-wide level. We devised an efficient sequencing strategy that allowed us to de novo assemble a major fraction of the S. mediterranea transcriptome. We then used independent assays and massive shotgun proteomics to validate the authenticity of transcripts. In total, our de novo assembly yielded 18,619 candidate transcripts with a mean length of 1118 nt after filtering. A total of 17,564 candidate transcripts could be mapped to 15,284 distinct loci on the current genome reference sequence. RACE confirmed complete or almost complete 5′ and 3′ ends for 22/24 transcripts. The frequencies of frame shifts, fusion, and fission events in the assembled transcripts were computationally estimated to be 4.2%–13%, 0%–3.7%, and 2.6%, respectively. Our shotgun proteomics produced 16,135 distinct peptides that validated 4200 transcripts (FDR ≤1%). The catalog of transcripts assembled in this study, together with the identified peptides, dramatically expands and refines planarian gene annotation, demonstrated by validation of several previously unknown transcripts with stem cell-dependent expression patterns. In addition, our robust transcriptome characterization pipeline could be applied to other organisms without genome assembly. All of our data, including homology annotation, are freely available at SmedGD, the S. mediterranea genome database.