Sequences upstream of the branch site are required to form helix II between U2 and U6 snRNA in a trans-splicing reaction

A nuclear function for an oncogenic microRNA as a modulator of snRNA and splicing

BACKGROUND

miRNAs are regulatory transcripts established as repressors of mRNA stability and translation that have been functionally implicated in carcinogenesis. miR-10b is one of the key onco-miRs associated with multiple forms of cancer. Malignant gliomas exhibit particularly striking dependence on miR-10b. However, despite the therapeutic potential of miR-10b targeting, this miRNA's poorly investigated and largely unconventional properties hamper the clinical translation.

METHODS

We utilized Covalent Ligation of Endogenous Argonaute-bound RNAs and their high-throughput RNA sequencing to identify miR-10b interactome and a combination of biochemical and imaging approaches for target validation. They included Crosslinking and RNA immunoprecipitation with spliceosomal proteins, a combination of miRNA FISH with protein immunofluorescence in glioma cells and patient-derived tumors, native Northern blotting, and the transcriptome-wide analysis of alternative splicing.

RESULTS

We demonstrate that miR-10b binds to U6 snRNA, a core component of the spliceosomal machinery. We provide evidence of the direct binding between miR-10b and U6, in situ imaging of miR-10b and U6 co-localization in glioma cells and tumors, and biochemical co-isolation of miR-10b with the components of the spliceosome. We further demonstrate that miR-10b modulates U6 N-6-adenosine methylation and pseudouridylation, U6 binding to splicing factors SART3 and PRPF8, and regulates U6 stability, conformation, and levels. These effects on U6 result in global splicing alterations, exemplified by the altered ratio of the isoforms of a small GTPase CDC42, reduced overall CDC42 levels, and downstream CDC42 -mediated effects on cell viability.

CONCLUSIONS

We identified U6 snRNA, the key RNA component of the spliceosome, as the top miR-10b target in glioblastoma. We, therefore, present an unexpected intersection of the miRNA and splicing machineries and a new nuclear function for a major cancer-associated miRNA.

A statistical approach for 5' splice site prediction using short sequence motifs and without encoding sequence data

Background

Most of the approaches for splice site prediction are based on machine learning techniques. Though, these approaches provide high prediction accuracy, the window lengths used are longer in size. Hence, these approaches may not be suitable to predict the novel splice variants using the short sequence reads generated from next generation sequencing technologies. Further, machine learning techniques require numerically encoded data and produce different accuracy with different encoding procedures. Therefore, splice site prediction with short sequence motifs and without encoding sequence data became a motivation for the present study.

Results

An approach for finding association among nucleotide bases in the splice site motifs is developed and used further to determine the appropriate window size. Besides, an approach for prediction of donor splice sites using sum of absolute error criterion has also been proposed. The proposed approach has been compared with commonly used approaches i.e., Maximum Entropy Modeling (MEM), Maximal Dependency Decomposition (MDD), Weighted Matrix Method (WMM) and Markov Model of first order (MM1) and was found to perform equally with MEM and MDD and better than WMM and MM1 in terms of prediction accuracy.

Conclusions

The proposed prediction approach can be used in the prediction of donor splice sites with higher accuracy using short sequence motifs and hence can be used as a complementary method to the existing approaches. Based on the proposed methodology, a web server was also developed for easy prediction of donor splice sites by users and is available at http://cabgrid.res.in:8080/sspred.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0362-6) contains supplementary material, which is available to authorized users.

Exon definition complexes contain the tri-snRNP and can be directly converted into B-like precatalytic splicing complexes

The first step in splicing of pre-mRNAs with long introns is exon definition, where U1 and U2 snRNPs bind at opposite ends of an exon. After exon definition, these snRNPs must form a complex across the upstream intron to allow splicing catalysis. Exon definition and conversion of cross-exon to cross-intron spliceosomal complexes are poorly understood. Here we demonstrate that, in addition to U1 and U2 snRNPs, cross-exon complexes contain U4, U5, and U6 (which form the tri-snRNP). Tri-snRNP docking involves the formation of U2/U6 helix II. This interaction is stabilized by a 5' splice site (SS)-containing oligonucleotide, which can bind the tri-snRNP and convert the cross-exon complex into a cross-intron, B-like complex. Our data suggest that the switch from cross-exon to cross-intron complexes can occur directly when an exon-bound tri-snRNP interacts with an upstream 5'SS, without prior formation of a cross-intron A complex, revealing an alternative spliceosome assembly pathway.

Human genetic variation recognizes functional elements in noncoding sequence

Noncoding DNA, particularly intronic DNA, harbors important functional elements that affect gene expression and RNA splicing. Yet, it is unclear which specific noncoding sites are essential for gene function and regulation. To identify functional elements in noncoding DNA, we characterized genetic variation within introns using ethnically diverse human polymorphism data from three public databases-PMT, NIEHS, and SeattleSNPs. We demonstrate that positions within introns corresponding to known functional elements involved in pre-mRNA splicing, including the branch site, splice sites, and polypyrimidine tract show reduced levels of genetic variation. Additionally, we observed regions of reduced genetic variation that are candidates for distance-dependent localization sites of functional elements, possibly intronic splicing enhancers (ISEs). Using several bioinformatics approaches, we provide additional evidence that supports our hypotheses that these regions correspond to ISEs. We conclude that studies of genetic variation can successfully discriminate and identify functional elements in noncoding regions. As more noncoding sequence data become available, the methods employed here can be utilized to identify additional functional elements in the human genome and provide possible explanations for phenotypic associations.

Large-scale comparative analysis of splicing signals and their corresponding splicing factors in eukaryotes

Introns are among the hallmarks of eukaryotic genes. Splicing of introns is directed by three main splicing signals: the 5' splice site (5'ss), the branch site (BS), and the polypyrimdine tract/3'splice site (PPT-3'ss). To study the evolution of these splicing signals, we have conducted a systematic comparative analysis of these signals in over 1.2 million introns from 22 eukaryotes. Our analyses suggest that all these signals have dramatically evolved: The PPT is weak among most fungi, intermediate in plants and protozoans, and strongest in metazoans. Within metazoans it shows a gradual strengthening from Caenorhabditis elegans to human. The 5'ss and the BS were found to be degenerate among most organisms, but highly conserved among some fungi. A maximum parsimony-based algorithm for reconstructing ancestral position-specific scoring matrices suggested that the ancestral 5'ss and BS were degenerate, as in metazoans. To shed light on the evolutionary variation in splicing signals, we have analyzed the evolutionary changes in the factors that bind these signals. Our analysis reveals coevolution of splicing signals and their corresponding splicing factors: The strength of the PPT is correlated to changes in key residues in its corresponding splicing factor U2AF2; limited correlation was found between changes in the 5'ss and U1 snRNA that binds it; but not between the BS and U2 snRNA. Thus, although the basic ability of eukaryotes to splice introns has remained conserved throughout evolution, the splicing signals and their corresponding splicing factors have considerably evolved, uniquely shaping the splicing mechanisms of different organisms.

DMS footprinting of structured RNAs and RNA-protein complexes

We describe a protocol in which dimethyl sulfate (DMS) modification of the base-pairing faces of unpaired adenosine and cytidine nucleotides is used for structural analysis of RNAs and RNA-protein complexes (RNPs). The protocol is optimized for RNAs of small to moderate size (< or = 500 nt). The RNA or RNP is first exposed to DMS under conditions that promote formation of the folded structure or complex, as well as 'control' conditions that do not allow folding or complex formation. The positions and extents of modification are then determined by primer extension, polyacrylamide gel electrophoresis and quantitative analysis. From changes in the extent of modification upon folding or protein binding (appearance of a 'footprint'), it is possible to detect local changes in the secondary and tertiary structure of RNA, as well as the formation of RNA-protein contacts. This protocol takes 1.5-3 d to complete, depending on the type of analysis used.

Dichotomous splicing signals in exon flanks

Intronic elements flanking the splice-site consensus sequences are thought to play a role in pre-mRNA splicing. However, the generality of this role, the catalog of effective sequences, and the mechanisms involved are still lacking. Using molecular genetic tests, we first showed that the approximately 50-nt intronic flanking sequences of exons beyond the splice-site consensus are generally important for splicing. We then went on to characterize exon flank sequences on a genomic scale. The G+C content of flanks displayed a bimodal distribution reflecting an exaggeration of this base composition in flanks relative to the gene as a whole. We divided all exons into two classes according to their flank G+C content and used computational and statistical methods to define pentamers of high relative abundance and phylogenetic conservation in exon flanks. Upstream pentamers were often common to the two classes, whereas downstream pentamers were totally different. Upstream and downstream pentamers were often identical around low G+C exons, and in contrast, were often complementary around high G+C exons. In agreement with this complementarity, predicted base pairing was more frequent between the flanks of high G+C exons. Pseudo exons did not exhibit this behavior, but rather tended to form base pairs between flanks and exon bodies. We conclude that most exons require signals in their immediate flanks for efficient splicing. G+C content is a sequence feature correlated with many genetic and genomic attributes. We speculate that there may be different mechanisms for splice site recognition depending on G+C content.

Allosteric cascade of spliceosome activation

Introns are removed from precursor messenger RNAs in the cell nucleus by a large ribonucleoprotein complex called the spliceosome. The spliceosome contains five subcomplexes called snRNPs, each with one RNA and several protein components. Interactions of the snRNPs with each other and the intron are highly dynamic, changing in an ordered progression throughout the splicing process. This allosteric cascade of interactions is programmed into the RNA and protein components of the spliceosome, and is driven by a family of DExD/H-box RNA-dependent ATPases. The dependence of cascade progression on multiple intron-recognition events likely serves to enforce the accuracy of splicing. Here, the progression of the allosteric cascade from the first recognition event to the first catalytic step of splicing is reviewed.

The U1 snRNP base pairs with the 5' splice site within a penta-snRNP complex

Recognition of the 5' splice site is an important step in mRNA splicing. To examine whether U1 approaches the 5' splice site as a solitary snRNP or as part of a multi-snRNP complex, we used a simplified in vitro system in which a short RNA containing the 5' splice site sequence served as a substrate in a binding reaction. This system allowed us to study the interactions of the snRNPs with the 5' splice site without the effect of other cis-regulatory elements of precursor mRNA. We found that in HeLa cell nuclear extracts, five spliceosomal snRNPs form a complex that specifically binds the 5' splice site through base pairing with the 5' end of U1. This system can accommodate RNA-RNA rearrangements in which U5 replaces U1 binding to the 5' splice site, a process that occurs naturally during the splicing reaction. The complex in which U1 and the 5' splice site are base paired sediments in the 200S fraction of a glycerol gradient together with all five spliceosomal snRNPs. This fraction is functional in mRNA spliceosome assembly when supplemented with soluble nuclear proteins. The results argue that U1 can bind the 5' splice site in a mammalian preassembled penta-snRNP complex.

The U5/U6 snRNA genomic repeat of Taenia solium

The U6 and U5 snRNA (small nuclear ribonucleic acid) genes were identified in Taenia solium with the aim of characterizing their sequence and genomic structures. They are contained within a shared 1,009-nt tandem genomic repeat and present at approximately 3 copies per haploid genome. The U6 snRNA gene shares 92 and 95% sequence similarity with the U6 homologs from humans and Schistosoma mansoni, respectively. The U5 snRNA gene of T. solium is 70% similar to the human U5 sequence in the 5' stem and loop 1 domains. The U6 and U5 snRNA genes are on complementary genomic strands and separated by 458 nt at their "heads" and 306 nt at their "tails." The nucleotides upstream of the U6 gene lack a recognizable TATA box and proximal sequence elements (PSEs), and the putative gene promoter for U5 snRNA does not resemble vertebrate examples. There are short blocks of similarity between the sequences upstream of the U5 and U6 snRNA genes, and these may be sites of shared transcription factor binding at the respective RNA polymerase II and III promoters. It is possible that this unusual allied U5/U6 snRNA genomic repeat may help mediate coordinated regulation of expression of the 2 snRNAs.

Central region of the human splicing factor Hprp3p interacts with Hprp4p

Human splicing factors Hprp3p and Hprp4p are associated with the U4/U6 small nuclear ribonucleoprotein particle, which is essential for the assembly of an active spliceosome. Currently, little is known about the specific roles of these factors in splicing. In this study, we characterized the molecular interaction between Hprp3p and Hprp4p. Constructs were created for expression of Hprp3p or its mutants in bacterial or mammalian cells. We showed that antibodies against either Hprp3p or Hprp4p were able to pull-down the Hprp3p-Hprp4p complex formed in Escherichia coli lysates. By co-immunoprecipitation and isothermal titration calorimetry, we demonstrated that purified Hprp3p and its mutants containing the central region, but lacking either the N-terminal 194 amino acids or the C-terminal 240 amino acids, were able to interact with Hprp4p. Conversely, Hprp3p mutants containing only the N- or C-terminal region did not interact with Hprp4p. In addition, by co-immunoprecipitation, we showed that intact Hprp3p and its mutants containing the central region interacted with Hprp4p in HeLa cell nuclear extracts. Primer extension analysis illustrated that the central region of Hprp3p is required to maintain the association of Hprp3p-Hprp4p with U4/U6 small nuclear RNAs, suggesting that this Hprp3p/Hprp4p interaction allows the recruitment of Hprp4p, and perhaps other protein(s), to the U4/U6 small nuclear ribonucleoprotein particle.

Probing interactions between the U2 small nuclear ribonucleoprotein and the DEAD-box protein, Prp5

Pre-mRNA binding to the yeast U2 small nuclear ribonucleoprotein (snRNP) during prespliceosome formation requires ATP hydrolysis, the highly conserved UACUAAC box of the branch point region of the pre-mRNA, and several factors. Here we analyzed the binding of a radiolabeled 2'-O-methyl oligonucleotide complementary to U2 small nuclear RNA to study interactions between the UACUAAC box, U2 snRNP, and Prp5p, a DEAD box protein necessary for prespliceosome formation. Binding of the 2'-O-methyl oligonucleotide to the U2 snRNP in yeast cell extract was assayed by gel electrophoresis. Binding was rapid, enhanced by ATP, and dependent on the integrity and conformation of the U2 snRNP. It was also stimulated by Prp5p that was found to associate physically with U2 snRNP. In vitro heat inactivation of the temperature-sensitive prp5-1 mutant extract decreased oligonucleotide binding to U2 and the ATP enhancement of binding by 3-fold. Furthermore, the temperature-sensitive prp5-1 mutation maps to the ATP-binding motif I within the helicase-like domain. Thus the catalytic activity of Prp5p likely promotes a conformational change in the U2 snRNP.

The ATP requirement for U2 snRNP addition is linked to the pre-mRNA region 5' to the branch site

Association of U2 snRNP with the pre-mRNA branch region is the first ATP-dependent step in spliceosome assembly. The basis of this energy dependence is not known. Previously, we identified minimal intron-derived substrates that form complexes with U2 independent of ATP. Here, we identify the intron region linked to the ATP dependence of this step by comparing these substrates to longer RNAs that recapitulate the ATP requirement. This region needed to impose ATP dependence lies immediately 5' to the branch site. Sequences ranging from 6 to 14 nt yield a near linear inhibitory effect on efficiency of complex formation with U2 snRNP, with 18 nt yielding near maximal ATP dependence. This region is not protected prior to U2 addition, and RNase H targeting of the region within nuclear extract converts an ATP-dependent substrate into an ATP-independent one. Within this region, there is no sequence specificity linked with the ATP requirement, as neither a specific sequence is needed, nor even nucleobases. These data and the results of other modifications suggest models in which the 18-nt region is a target for interactions with U2 snRNP in an ATP-bound or -activated conformation.