Search

Short RNA expression was analyzed from Dicer -positive and Dicer -knockout mouse embroyonic stem (ES) cells, using high-throughput pyrosequencing. A correlation of miRNA quantification with sequencing frequency estimates that there are 110,000 miRNAs per ES cell, the majority of which can be accounted for by six distinct miRNA loci. Four of these miRNA loci or their human homologues have demonstrated roles in cell cycle regulation or oncogenesis, suggesting that a major function of the miRNA pathway in ES cells may be to shape their distinct cell cycle. Forty-six previously uncharacterized miRNAs were identified, most of which are expressed at low levels and are less conserved than the set of known miRNAs. Low-abundance short RNAs matching all classes of repetitive elements were present in cells lacking Dicer , although the production of some SINE- and simple repeat-associated short RNAs appeared to be Dicer-dependent. These and other Dicer-dependent sequences resembled miRNAs. At a depth of sequencing that approaches the total number of 5′ phosphorylated short RNAs per cell, miRNAs appeared to be Dicer's only substrate. The results presented suggest a model in which repeat-associated miRNAs serve as host defenses against repetitive elements, a function canonically ascribed to other classes of short RNA.

Knowledge of the functional cis-regulatory elements that regulate constitutive and alternative pre-mRNA splicing is fundamental for biology and medicine. Here we undertook a genome-wide comparative genomics approach using available mammalian genomes to identify conserved intronic splicing regulatory elements (ISREs). Our approach yielded 314 ISREs, and insertions of ~70 ISREs between competing splice sites demonstrated that 84% of ISREs altered 5′ and 94% altered 3′ splice site choice in human cells. Consistent with our experiments, comparisons of ISREs to known splicing regulatory elements revealed that 40%–45% of ISREs might have dual roles as exonic splicing silencers. Supporting a role for ISREs in alternative splicing, we found that 30%–50% of ISREs were enriched near alternatively spliced (AS) exons, and included almost all known binding sites of tissue-specific alternative splicing factors. Further, we observed that genes harboring ISRE-proximal exons have biases for tissue expression and molecular functions that are ISRE-specific. Finally, we discovered that for Nova1, neuronal PTB, hnRNP C, and FOX1, the most frequently occurring ISRE proximal to an alternative conserved exon in the splicing factor strongly resembled its own known RNA binding site, suggesting a novel application of ISRE density and the propensity for splicing factors to auto-regulate to associate RNA binding sites to splicing factors. Our results demonstrate that ISREs are crucial building blocks in understanding general and tissue-specific AS regulation and the biological pathways and functions regulated by these AS events., Author Summary During RNA splicing, sequences (introns) in a pre-mRNA are excised and discarded, and the remaining sequences (exons) are joined to form the mature RNA. Splicing is regulated not only by the binding of the basic splicing machinery to splice sites located at the exon–intron boundaries, but also by the combined effects of various other splicing factors that bind to a multitude of sequence elements located both in the exons as well as the flanking introns. Instances of alternative splicing, where usage of splice site(s) is incomplete or different between tissues, cell types, or lineages, can be created by the interaction of sequence elements and tissue, cell type, and stage-specific splicing factors. To better understand constitutive and alternative pre-mRNA splicing, the authors describe a comparative genomics approach, using available mammalian genomes, to systematically identify splicing regulatory elements located in the introns proximal to exons. A quarter of the elements were tested experimentally, and most of them altered splicing in human cells. The authors also showed that that the intronic elements are close to tissue-specific alternative exons and are more likely to be located in specific positions in the introns, suggestive of potential regulatory function. These elements are also frequently found in tissue-specific genes, suggesting a coupling between expression and alternative splicing of these genes. Finally, the authors propose a strategy using the elements to identify the binding sites of several splicing factors.

The central nervous system (CNS) is arguably one of the most complex systems in the universe. To understand the CNS, scientists have investigated a variety of molecules, including proteins, lipids, and various small molecules. However, one large class of molecules, noncoding RNAs (ncRNAs), has been relatively unexplored. ncRNAs function directly as structural, catalytic, or regulatory molecules rather than serving as templates for protein synthesis. The increasing variety of ncRNAs being identified in the CNS suggests a strong connection between the biogenesis, dynamics of action, and combinatorial regulatory potential of ncRNAs and the complexity of the CNS. In this review, we give an overview of the diversity and abundance of ncRNAs before delving into specific examples that illustrate their importance in the CNS. In particular, we cover recent evidence for the roles of microRNAs, small nucleolar RNAs, retrotransposons, the NRSE small modulatory RNA, and BC1/BC200 in the CNS. Finally, we speculate why ncRNAs are well adapted to improving organism-environment interactions.

Although core mechanisms and machinery of premRNA splicing are conserved from yeast to human, the details of intron recognition often differ, even between closely related organisms. For example, genes from the pufferfish Fugu rubripes generally contain one or more introns that are not properly spliced in mouse cells. Exploiting available genome sequence data, a battery of sequence analysis techniques was used to reach several conclusions about the organization and evolution of splicing regulatory elements in vertebrate genes. The classical splice site and putative branch site signals are completely conserved across the vertebrates studied (human, mouse, pufferfish, and zebrafish), and exonic splicing enhancers also appear broadly conserved in vertebrates. However, another class of splicing regulatory elements, the intronic splicing enhancers, appears to differ substantially between mammals and fish, with G triples (GGG) very abundant in mammalian introns but comparatively rare in fish. Conversely, short repeats of AC and GT are predicted to function as intronic splicing enhancers in fish but are not enriched in mammalian introns. Consistent with this pattern, exonic splicing enhancer-binding SR proteins are highly conserved across all vertebrates, whereas heterogeneous nuclear ribonucleoproteins, which bind many intronic sequences, vary in domain structure and even presence/absence between mammals and fish. Exploiting differences in intronic sequence composition, a statistical model was developed to predict the splicing phenotype of Fugu introns in mammalian systems and was used to engineer the spliceability of a Fugu intron in human cells by insertion of specific sequences, thereby rescuing splicing in human cells.

Alternative pre-mRNA splicing is widely used to regulate gene expression by tuning the levels of tissue-specific mRNA isoforms. Few regulatory mechanisms are understood at the level of combinatorial control despite numerous sequences, distinct from splice sites, that have been shown to play roles in splicing enhancement or silencing. Here we use molecular approaches to identify a ternary combination of exonic UAGG and 5′-splice-site-proximal GGGG motifs that functions cooperatively to silence the brain-region-specific CI cassette exon (exon 19) of the glutamate NMDA R1 receptor (GRIN1) transcript. Disruption of three components of the motif pattern converted the CI cassette into a constitutive exon, while predominant skipping was conferred when the same components were introduced, de novo, into a heterologous constitutive exon. Predominant exon silencing was directed by the motif pattern in the presence of six competing exonic splicing enhancers, and this effect was retained after systematically repositioning the two exonic UAGGs within the CI cassette. In this system, hnRNP A1 was shown to mediate silencing while hnRNP H antagonized silencing. Genome-wide computational analysis combined with RT-PCR testing showed that a class of skipped human and mouse exons can be identified by searches that preserve the sequence and spatial configuration of the UAGG and GGGG motifs. This analysis suggests that the multi-component silencing code may play an important role in the tissue-specific regulation of the CI cassette exon, and that it may serve more generally as a molecular language to allow for intricate adjustments and the coordination of splicing patterns from different genes., Many genes are alternatively spliced, but the signals that regulate the process are unclear. These authors have found a sequence motif that appears to function at many alternatively spliced genes

We propose a framework for modeling sequence motifs based on the maximum entropy principle (MEP). We recommend approximating short sequence motif distributions with the maximum entropy distribution (MED) consistent with low-order marginal constraints estimated from available data, which may include dependencies between nonadjacent as well as adjacent positions. Many maximum entropy models (MEMs) are specified by simply changing the set of constraints. Such models can be utilized to discriminate between signals and decoys. Classification performance using different MEMs gives insight into the relative importance of dependencies between different positions. We apply our framework to large datasets of RNA splicing signals. Our best models out-perform previous probabilistic models in the discrimination of human 5' (donor) and 3' (acceptor) splice sites from decoys. Finally, we discuss mechanistically motivated ways of comparing models.

Rhabdomyosarcoma is the most commonly diagnosed pediatric soft tissue cancer. Tumors may develop in several locations, but all tumors express markers of undifferentiated skeletal muscle tissue. Of the two main rhabdomyosarcoma subtypes, the alveolar subtype leads to a much poorer prognosis than the embryonal subtype. Alveolar rhabdomyosarcoma is characterized by a chromosomal translocation between either chromosomes 1 and 13, or more commonly 2 and 13, which creates the fusion genes PAX7-FKHR or PAX3-FKHR, respectively. The fusion genes encode a protein that contains the PAX DNA binding domain and the FKHR transactivation domain, but the mechanism of disease of the fusion protein has not yet been fully elucidated. Here we use gain of function and loss of function techniques to study the effect of PAX3-FKHR on miRNA expression. Using small RNA-seq, we show that PAX3-FKHR expression leads to an upregulation of microRNAs mir-495 and miR-543. Examination of the mir-495 and mir-543 genes revealed conserved PAX3 binding sites. We further investigate the effect of the microRNAs on rhabdomyosarcoma cell differentiation and tumorigenicity. Citation Format: Shannon Muir, Jason Nathanson, Melissa Wilbert, Gene Yeo, Frank Furnari, Karen Arden, Webster Cavenee. The role of miRNA in PAX3-FKHR positive rhabdomyosarcoma. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3558. doi:10.1158/1538-7445.AM2014-3558

Human embryonic stem cells (hESCs) and neural progenitor (NP) cells are excellent models for recapitulating early neuronal development in vitro, and are key to establishing strategies for the treatment of degenerative disorders. While much effort had been undertaken to analyze transcriptional and epigenetic differences during the transition of hESC to NP, very little work has been performed to understand post-transcriptional changes during neuronal differentiation. Alternative RNA splicing (AS), a major form of post-transcriptional gene regulation, is important in mammalian development and neuronal function. Human ESC, hESC-derived NP, and human central nervous system stem cells were compared using Affymetrix exon arrays. We introduced an outlier detection approach, REAP (Regression-based Exon Array Protocol), to identify 1,737 internal exons that are predicted to undergo AS in NP compared to hESC. Experimental validation of REAP-predicted AS events indicated a threshold-dependent sensitivity ranging from 56% to 69%, at a specificity of 77% to 96%. REAP predictions significantly overlapped sets of alternative events identified using expressed sequence tags and evolutionarily conserved AS events. Our results also reveal that focusing on differentially expressed genes between hESC and NP will overlook 14% of potential AS genes. In addition, we found that REAP predictions are enriched in genes encoding serine/threonine kinase and helicase activities. An example is a REAP-predicted alternative exon in the SLK (serine/threonine kinase 2) gene that is differentially included in hESC, but skipped in NP as well as in other differentiated tissues. Lastly, comparative sequence analysis revealed conserved intronic cis-regulatory elements such as the FOX1/2 binding site GCAUG as being proximal to candidate AS exons, suggesting that FOX1/2 may participate in the regulation of AS in NP and hESC. In summary, a new methodology for exon array analysis was introduced, leading to new insights into the complexity of AS in human embryonic stem cells and their transition to neural stem cells., Author Summary Deriving neural progenitors (NP) from human embryonic stem cells (hESC) is the first step in creating homogeneous populations of cells that will differentiate into myriad neuronal subtypes necessary to form a human brain. During alternative RNA splicing (AS), noncoding sequences (introns) in a pre-mRNA are differentially removed in different cell types and tissues, and the remaining sequences (exons) are joined to form multiple forms of mature RNA, playing an important role in cellular diversity. The authors utilized Affymetrix exon arrays with probes targeting hundreds of thousands of exons to study AS comparing human ES to NP. To accomplish this, a novel computational method, REAP (Regression-based Exon Array Protocol), is introduced to analyze the exon array data. The authors showed that REAP candidates are consistent with other types of methods for discovering alternative exons. In addition, REAP candidate alternative exons are enriched in genes encoding serine/theronine kinases and helicase activities. An example is the alternative exon in the SLK (serine/threonine kinase 2) gene that is included in hESC, but excluded in NP as well as in other differentiated tissues. Finally, by comparing genomic sequences across multiple mammals, the authors identified dozens of conserved candidate binding sites that were enriched proximal to REAP candidate exons.

In PLoS Genetics, volume 3, issue 5: 10.1371/journal.pgen.0030085 Mr. Eric L. Van Nostrand's name was incorrectly listed in the online citation as Nostrand ELV. The correct citation is: Yeo GW, Van Nostrand EL, Liang TY (2007) Discovery and Analysis of Evolutionarily Conserved Intronic Splicing Regulatory Elements. PLoS Genet 3(5): e85 doi:10.1371/journal.pgen.0030085.

Silencing of splicing regulators by RNA interference, combined with splicing-specific microarrays, has revealed a complex network of distinct alternative splicing events in Drosophila, while a high-throughput screen of more than 6,000 compounds has identified drugs that interfere specifically and directly with one class of splicing regulators in human cells.

Alternative pre-mRNA splicing (AS) is widely used by higher eukaryotes to generate different protein isoforms in specific cell or tissue types. To compare AS events across human tissues, we analyzed the splicing patterns of genomically aligned expressed sequence tags (ESTs) derived from libraries of cDNAs from different tissues. Controlling for differences in EST coverage among tissues, we found that the brain and testis had the highest levels of exon skipping. The most pronounced differences between tissues were seen for the frequencies of alternative 3' splice site and alternative 5' splice site usage, which were about 50 to 100% higher in the liver than in any other human tissue studied. Quantifying differences in splice junction usage, the brain, pancreas, liver and the peripheral nervous system had the most distinctive patterns of AS. Analysis of available microarray expression data showed that the liver had the most divergent pattern of expression of serine-arginine protein and heterogeneous ribonucleoprotein genes compared to the other human tissues studied, possibly contributing to the unusually high frequency of alternative splice site usage seen in liver. Sequence motifs enriched in alternative exons in genes expressed in the brain, testis and liver suggest specific splicing factors that may be important in AS regulation in these tissues. This study distinguishes the human brain, testis and liver as having unusually high levels of AS, highlights differences in the types of AS occurring commonly in different tissues, and identifies candidate cis-regulatory elements and trans-acting factors likely to have important roles in tissue-specific AS in human cells.

SummaryAlternative splicing (AS) defects that adversely affect gene expression and function have been identified in diabetic hearts; however, the mechanisms responsible are largely unknown. Here, we show that the RNA-binding protein RBFOX2 contributes to transcriptome changes under diabetic conditions. RBFOX2 controls AS of genes with important roles in heart function relevant to diabetic cardiomyopathy. RBFOX2 protein levels are elevated in diabetic hearts despite low RBFOX2 AS activity. A dominant-negative (DN) isoform of RBFOX2 that blocks RBFOX2-mediated AS is generated in diabetic hearts. DN RBFOX2 interacts with wild-type (WT) RBFOX2, and ectopic expression of DN RBFOX2 inhibits AS of RBFOX2 targets. Notably, DN RBFOX2 expression is specific to diabetes and occurs at early stages before cardiomyopathy symptoms appear. Importantly, DN RBFOX2 expression impairs intracellular calcium release in cardiomyocytes. Our results demonstrate that RBFOX2 dysregulation by DN RBFOX2 is an early pathogenic event in diabetic hearts.

SummaryThe generation of pancreas, liver, and intestine from a common pool of progenitors in the foregut endoderm requires the establishment of organ boundaries. How dorsal foregut progenitors activate pancreatic genes and evade the intestinal lineage choice remains unclear. Here, we identify Pdx1 and Sox9 as cooperative inducers of a gene regulatory network that distinguishes the pancreatic from the intestinal lineage. Genetic studies demonstrate dual and cooperative functions for Pdx1 and Sox9 in pancreatic lineage induction and repression of the intestinal lineage choice. Pdx1 and Sox9 bind to regulatory sequences near pancreatic and intestinal differentiation genes and jointly regulate their expression, revealing direct cooperative roles for Pdx1 and Sox9 in gene activation and repression. Our study identifies Pdx1 and Sox9 as important regulators of a transcription factor network that initiates pancreatic fate and sheds light on the gene regulatory circuitry that governs the development of distinct organs from multi-lineage-competent foregut progenitors.

SummaryUnderstanding how RNA binding proteins control the splicing code is fundamental to human biology and disease. Here, we present a comprehensive study to elucidate how heterogeneous nuclear ribonucleoparticle (hnRNP) proteins, among the most abundant RNA binding proteins, coordinate to regulate alternative pre-mRNA splicing (AS) in human cells. Using splicing-sensitive microarrays, crosslinking and immunoprecipitation coupled with high-throughput sequencing (CLIP-seq), and cDNA sequencing, we find that more than half of all AS events are regulated by multiple hnRNP proteins and that some combinations of hnRNP proteins exhibit significant synergy, whereas others act antagonistically. Our analyses reveal position-dependent RNA splicing maps, in vivo consensus binding sites, a surprising level of cross- and autoregulation among hnRNP proteins, and the coordinated regulation by hnRNP proteins of dozens of other RNA binding proteins and genes associated with cancer. Our findings define an unprecedented degree of complexity and compensatory relationships among hnRNP proteins and their splicing targets that likely confer robustness to cells.

Human pluripotent stem cells (hPSCs) require precise control of post-transcriptional RNA networks to maintain proliferation and survival. Using enhanced UV crosslinking and immunoprecipitation (eCLIP), we identify RNA targets of the IMP/IGF2BP family of RNA-binding proteins in hPSCs. At the broad region and binding site levels, IMP1 and IMP2 show reproducible binding to a large and overlapping set of 3′ UTR-enriched targets. RNA Bind-N-seq applied to recombinant full-length IMP1 and IMP2 reveals CA-rich motifs that are enriched in eCLIP-defined binding sites. We observe that IMP1 loss in hPSCs recapitulates IMP1 phenotypes, including a reduction in cell adhesion and increase in cell death. For cell adhesion, we find IMP1 maintains levels of integrin mRNA specifically regulating RNA stability of ITGB5 in hPSCs. Additionally, we show that IMP1 can be linked to hPSC survival via direct target BCL2. Thus, transcriptome-wide binding profiles identify hPSC targets modulating well-characterized IMP1 roles., National Institutes of Health (U.S.) (Grant HG007005)

Summary Proper neuronal function and several forms of synaptic plasticity are highly dependent on precise control of mRNA translation, particularly in dendrites. We find that eIF4AIII, a core exon junction complex (EJC) component loaded onto mRNAs by pre-mRNA splicing, is associated with neuronal mRNA granules and dendritic mRNAs. eIF4AIII knockdown markedly increases both synaptic strength and GLUR1 AMPA receptor abundance at synapses. eIF4AIII depletion also increases ARC, a protein required for maintenance of long-term potentiation; arc mRNA, one of the most abundant in dendrites, is a natural target for nonsense-mediated decay (NMD). Numerous new NMD candidates, some with potential to affect synaptic activity, were also identified computationally. Two models are presented for how translation-dependent decay pathways such as NMD might advantageously function as critical brakes for protein synthesis in cells such as neurons that are highly dependent on spatially and temporally restricted protein expression.

SummaryAutism spectrum disorders (ASD) are complex neurodevelopmental diseases in which different combinations of genetic mutations may contribute to the phenotype. Using Rett syndrome (RTT) as an ASD genetic model, we developed a culture system using induced pluripotent stem cells (iPSCs) from RTT patients' fibroblasts. RTT patients' iPSCs are able to undergo X-inactivation and generate functional neurons. Neurons derived from RTT-iPSCs had fewer synapses, reduced spine density, smaller soma size, altered calcium signaling and electrophysiological defects when compared to controls. Our data uncovered early alterations in developing human RTT neurons. Finally, we used RTT neurons to test the effects of drugs in rescuing synaptic defects. Our data provide evidence of an unexplored developmental window, before disease onset, in RTT syndrome where potential therapies could be successfully employed. Our model recapitulates early stages of a human neurodevelopmental disease and represents a promising cellular tool for drug screening, diagnosis and personalized treatment.PaperFlick

The N6-methyladenosine (m6A) modification is the most prevalent post-transcriptional mRNA modification, regulating mRNA decay, translation and splicing. It plays a major role during normal development, differentiation, and disease progression. The modification is dynamically regulated by a set of writer, eraser and reader proteins. The YTH-domain family of proteins: Ythdf1, Ythdf2, and Ythdf3, are three homologous m6A binding proteins, which have different cellular functions. However, their sequence similarity and their tendency to bind the same targets suggest that they may have overlapping roles. We systematically knocked out (KO) the Mettl3 writer for each of the Ythdf readers and for the three readers together (triple-KO). We then estimated the effect in-vivo, in mouse gametogenesis and viability, and in-vitro, in mouse embryonic stem cells (mESCs). We show that in gametogenesis, Mettl3-KO severity is increased as the deletion occurs earlier in the process, and Ythdf2 has a dominant role that cannot be compensated by Ythdf1 or Ythdf3, possibly due to differences in readers’ expression, both in quantity and in spatial location. By knocking out the three readers together and systematically testing offspring genotypes, we have revealed a redundancy in the readers’ role during early development, a redundancy which is dosage-dependent. Additionally, we show that in mESCs there is compensation between the three readers, since the inability to differentiate and the significant effect on mRNA decay occur only in the triple-KO cells and not in the single KOs. Thus, we suggest a novel model for the Ythdf readers function. There is a dosage-dependent redundancy when all three readers are co-expressed in the same location in the cells.