Whole-genome sequencing of a laboratory-evolved yeast strain

Background

Two categories of introns are known, a common U2 type and a rare U12 type. These two types of introns are removed by distinct spliceosomes. The phylogenetic distribution of spliceosomal RNAs that are characteristic of the U12 spliceosome, i.e. the U11, U12, U4atac and U6atac RNAs, suggest that U12 spliceosomes were lost in many phylogenetic groups. We have now examined the distribution of U2 and U12 introns in many of these groups.

Results

U2 and U12 introns were predicted by making use of available EST and genomic sequences. The results show that in species or branches where U12 spliceosomal components are missing, also U12 type of introns are lacking. Examples are the choanoflagellate Monosiga brevicollis, Entamoeba histolytica, green algae, diatoms, and the fungal lineage Basidiomycota. Furthermore, whereas U12 splicing does not occur in Caenorhabditis elegans, U12 introns as well as U12 snRNAs are present in Trichinella spiralis, which is deeply branching in the nematode tree. A comparison of homologous genes in T. spiralis and C. elegans revealed different mechanisms whereby U12 introns were lost.

Conclusions

The phylogenetic distribution of U12 introns and spliceosomal RNAs give further support to an early origin of U12 dependent splicing. In addition, this distribution identifies a large number of instances during eukaryotic evolution where such splicing was lost.     

Evolutionary dynamics of U12-type spliceosomal introns

Background  

Many multicellular eukaryotes have two types of spliceosomes for the removal of introns from messenger RNA precursors. The major (U2) spliceosome processes the vast majority of introns, referred to as U2-type introns, while the minor (U12) spliceosome removes a small fraction (less than 0.5%) of introns, referred to as U12-type introns. U12-type introns have distinct sequence elements and usually occur together in genes with U2-type introns. A phylogenetic distribution of U12-type introns shows that the minor splicing pathway appeared very early in eukaryotic evolution and has been lost repeatedly.     

The abundance of the spliceosomal snRNPs is not limiting the splicing of U12-type introns

The rate of excision of U12-type introns has been reported to be slower than that of U2-type introns, suggesting a rate-limiting bottleneck that could down-regulate genes containing U12-type introns. The mechanistic reasons for this slower rate of intron excision are not known, but lower abundance of the U12-type snRNPs and slower rate of assembly or catalytic activity have been suggested. To investigate snRNP abundance we concentrated on the U4atac snRNA, which is the least abundant of the U12-type snRNAs and is limiting the formation of U4atac/U6atac complex. We identified mouse NIH-3T3 cell line isolates in which the level of both U4atac snRNA and U4atac/U6atac complexes is reduced to 10%-20% of the normal level. We used these cell lines to investigate splicing efficiency by transient transfection of a reporter gene containing a U12-type intron and by quantitative PCR analysis of endogenous genes. The splicing of the reporter U12-type intron was very inefficient, but the activity could be restored by overexpression of U4atac snRNA. Using these U4atac-deficient NIH-3T3 cells, we confirmed the results of previous studies showing that U12-type introns of endogenous genes are, indeed, excised more slowly than U2-type introns, but we found that the rate did not differ from that measured in cells displaying normal levels of U4atac snRNA. Thus our results suggest that the cellular abundance of the snRNPs does not limit U12-type intron splicing under normal conditions.     

Computational screen for spliceosomal RNA genes aids in defining the phylogenetic distribution of major and minor spliceosomal components

The RNA molecules of the spliceosome are critical for specificity and catalysis during splicing of eukaryotic pre-mRNA. In order to examine the evolution and phylogenetic distribution of these RNAs, we analyzed 149 eukaryotic genomes representing a broad range of phylogenetic groups. RNAs were predicted using high-sensitivity local alignment methods and profile HMMs in combination with covariance models. The results provide the most comprehensive view so far of the phylogenetic distribution of spliceosomal RNAs. RNAs were predicted in many phylogenetic groups where these RNA were not previously reported. Examples are RNAs of the major (U2-type) spliceosome in all fungal lineages, in lower metazoa and many protozoa. We also identified the minor (U12-type) spliceosomal U11 and U6atac RNAs in Acanthamoeba castellanii, where U12 spliceosomal RNA as well as minor introns were reported recently. In addition, minor-spliceosome-specific RNAs were identified in a number of phylogenetic groups where previously such RNAs were not observed, including the nematode Trichinella spiralis, the slime mold Physarum polycephalum and the fungal lineages Zygomycota and Chytridiomycota. The detailed map of the distribution of the U12-type RNA genes supports an early origin of the minor spliceosome and points to a number of occasions during evolution where it was lost.     

Both U2 snRNA and U12 snRNA are required for accurate splicing of exon 5 of the rat calcitonin/CGRP gene

Two classes of spliceosome are present in eukaryotic cells. Most introns in nuclear pre-mRNAs are removed by a spliceosome that requires U1, U2, U4, U5, and U6 small nuclear ribonucleoprotein particles (snRNPs). A minor class of introns are removed by a spliceosome containing U11, U12, U5, U4atac, and U6 atac snRNPs. We describe experiments that demonstrate that splicing of exon 5 of the rat calcitonin/CGRP gene requires both U2 snRNA and U12 snRNA. In vitro, splicing to calcitonin/ CGRP exon 5 RNA was dependent on U2 snRNA, as preincubation of nuclear extract with an oligonucleotide complementary to U2 snRNA abolished exon 5 splicing. Addition of an oligonucleotide complementary to U12 snRNA increased splicing at a cryptic splice site in exon 5 from <5% to 50% of total spliced RNA. Point mutations in a candidate U12 branch sequence in calcitonin/CGRP intron 4, predicted to decrease U12-pre-mRNA base-pairing, also significantly increased cryptic splicing in vitro. Calcitonin/CGRP genes containing base changes disrupting the U12 branch sequence expressed significantly decreased CGRP mRNA levels when expressed in cultured cells. Coexpression of U12 snRNAs containing base changes predicted to restore U12-pre-mRNA base pairing increased CGRP mRNA synthesis to the level of the wild-type gene. These observations indicate that accurate, efficient splicing of calcitonin/CGRP exon 5 is dependent upon both U2 and U12 snRNAs.     

Identification,characterization and molecular phylogeny of U12-dependent introns in the Arabidopsis thaliana genome

U12-dependent introns are spliced by the minor U12-type spliceosome and occur in a variety of eukaryotic organisms, including Arabidopsis. In this study, a set of putative U12-dependent introns was compiled from a large collection of cDNA/EST- confirmed introns in the Arabidopsis thaliana genome by means of high-throughput bioinformatic analysis combined with manual scrutiny. A total of 165 U12-type introns were identified based upon stringent criteria. This number of sequences well exceeds the total number of U12-type introns previously reported for plants and allows a more thorough statistical analysis of U12-type signals. Of particular note is the discovery that the distance between the branch site adenosine and the acceptor site ranges from 10 to 39 nt, significantly longer than the previously postulated limit of 21 bp. Further analysis indicates that, in addition to the spacing constraint, the sequence context of the potential acceptor site may have an important role in 3′ splice site selection. Several alternative splicing events involving U12-type introns were also captured in this study, providing evidence that U12-dependent acceptor sites can also be recognized by the U2-type spliceosome. Furthermore, phylogenetic analysis suggests that both U12-type AT-AC and U12-type GT-AG introns occurred in Na⁺/H⁺ antiporters in a progenitor of animals and plants.     

Base pairing with U6atac snRNA is required for 5' splice site activation of U12-dependent introns in vivo.

The minor U12-dependent class of eukaryotic nuclear pre-mRNA introns is spliced by a distinct spliceosomal mechanism that requires the function of U11, U12, U5, U4atac, and U6atac snRNAs. Previous work has shown that U11 snRNA plays a role similar to U1 snRNA in the major class spliceosome by base pairing to the conserved 5'' splice site sequence. Here we show that U6atac snRNA also base pairs to the 5'' splice site in a manner analogous to that of U6 snRNA in the major class spliceosome. We show that splicing defective mutants of the 5'' splice site can be activated for splicing in vivo by the coexpression of compensatory U6atac snRNA mutants. In some cases, maximal restoration of splicing required the coexpression of compensatory U11 snRNA mutants. The allelic specificity of mutant phenotype suppression is consistent with Watson-Crick base pairing between the pre-mRNA and the snRNAs. These results provide support for a model of the RNA-RNA interactions at the core of the U12-dependent spliceosome that is strikingly similar to that of the major class U2-dependent spliceosome.     

Recycling of the U12-type spliceosome requires p110, a component of the U6atac snRNP

U12-dependent introns are spliced by the so-called minor spliceosome, requiring the U11, U12, and U4atac/U6atac snRNPs in addition to the U5 snRNP. We have recently identified U6-p110 (SART3) as a novel human recycling factor that is related to the yeast splicing factor Prp24. U6-p110 transiently associates with the U6 and U4/U6 snRNPs during the spliceosome cycle, regenerating functional U4/U6 snRNPs from singular U4 and U6 snRNPs. Here we investigated the involvement of U6-p110 in recycling of the U4atac/U6atac snRNP. In contrast to the major U6 and U4/U6 snRNPs, p110 is primarily associated with the U6atac snRNP but is almost undetectable in the U4atac/U6atac snRNP. Since p110 does not occur in U5 snRNA-containing complexes, it appears to be transiently associated with U6atac during the cycle of the minor spliceosome. The p110 binding site was mapped to U6 nucleotides 38 to 57 and U6atac nucleotides 10 to 30, which are highly conserved between these two functionally related snRNAs. With a U12-dependent in vitro splicing system, we demonstrate that p110 is required for recycling of the U4atac/U6atac snRNP.     

Evolutionary conservation of minor U12-type spliceosome between plants and humans

Splicing of rare, U12-type or AT-AC introns is mediated by a distinct spliceosome that assembles from U11, U12, U4atac, U6atac, and U5 snRNPs. Although in human cells the protein composition of minor and major snRNPs is similar, differences, particularly in U11 and U12 snRNPs, have been recently described. We have identified an Arabidopsis U11 snRNP-specific 35K protein as an interacting partner of an RS-domain-containing cyclophilin. By using a transient expression system in Arabidopsis protoplasts, we show that the 35K protein incorporates into snRNP. Oligo affinity selection and glycerol gradient centrifugation revealed that the Arabidopsis 35K protein is present in monomeric U11 snRNP and in U11/U12-di snRNP. The interaction of the 35K protein with Arabidopsis SR proteins together with its strong sequence similarity to U1-70K suggests that its function in splicing of minor introns is analogous to that of U1-70K. Analysis of Arabidopsis and Oryza sativa genome sequences revealed that all U11/U12-di-snRNP-specific proteins are conserved in dicot and monocot plants. In addition, we have identified an Arabidopsis gene encoding the homolog of U4atac snRNA and a second Arabidopsis gene encoding U6atac snRNA. Secondary structure predictions indicate that the Arabidopsis U4atac is able to form dimeric complexes with both Arabidopsis U6atac snRNAs. As revealed by RNaseA/T1 protection assay, the U4atac snRNA gene is expressed as an ~160-nt RNA, whereas the second U6atac snRNA gene seems to be a pseudogene. Taken together, our data indicate that recognition and splicing of minor, AT-AC introns in plants is highly similar to that in humans.     

The intramolecular stem-loop structure of U6 snRNA can functionally replace the U6atac snRNA stem-loop

The U6 spliceosomal snRNA forms an intramolecular stem-loop structure during spliceosome assembly that is required for splicing and is proposed to be at or near the catalytic center of the spliceosome. U6atac snRNA, the analog of U6 snRNA used in the U12-dependent splicing of the minor class of spliceosomal introns, contains a similar stem-loop whose structure but not sequence is conserved between humans and plants. To determine if the U6 and U6atac stem-loops are functionally analogous, the stem-loops from human and budding yeast U6 snRNAs were substituted for the U6atac snRNA structure and tested in an in vivo genetic suppression assay. Both chimeric U6/U6atac snRNA constructs were active for splicing in vivo. In contrast, several mutations of the native U6atac stem-loop that either delete putatively unpaired residues or disrupt the putative stem regions were inactive for splicing. Compensatory mutations that are expected to restore base pairing within the stem regions restored splicing activity. However, other mutants that retained base pairing potential were inactive, suggesting that functional groups within the stem regions may contribute to function. These results show that the U6atac snRNA stem-loop structure is required for in vivo splicing within the U12-dependent spliceosome and that its role is likely to be similar to that of the U6 snRNA intramolecular stem-loop.     

The RNA-binding protein hnRNPLL induces a T cell alternative splicing program delineated by differential intron retention in polyadenylated RNA

Background

Retention of a subset of introns in spliced polyadenylated mRNA is emerging as a frequent, unexplained finding from RNA deep sequencing in mammalian cells.

Results

Here we analyze intron retention in T lymphocytes by deep sequencing polyadenylated RNA. We show a developmentally regulated RNA-binding protein, hnRNPLL, induces retention of specific introns by sequencing RNA from T cells with an inactivating Hnrpll mutation and from B lymphocytes that physiologically downregulate Hnrpll during their differentiation. In Ptprc mRNA encoding the tyrosine phosphatase CD45, hnRNPLL induces selective retention of introns flanking exons 4 to 6; these correspond to the cassette exons containing hnRNPLL binding sites that are skipped in cells with normal, but not mutant or low, hnRNPLL. We identify similar patterns of hnRNPLL-induced differential intron retention flanking alternative exons in 14 other genes, representing novel elements of the hnRNPLL-induced splicing program in T cells. Retroviral expression of a normally spliced cDNA for one of these targets, Senp2, partially corrects the survival defect of Hnrpll-mutant T cells. We find that integrating a number of computational methods to detect genes with differentially retained introns provides a strategy to enrich for alternatively spliced exons in mammalian RNA-seq data, when complemented by RNA-seq analysis of purified cells with experimentally perturbed RNA-binding proteins.

Conclusions

Our findings demonstrate that intron retention in mRNA is induced by specific RNA-binding proteins and suggest a biological significance for this process in marking exons that are poised for alternative splicing.     

Minor-class splicing occurs in the nucleus of the Xenopus oocyte

A small fraction of premessenger RNA introns in certain eukaryotes is excised by the minor spliceosome, which contains low-abundance small nuclear ribonucleoproteins (snRNPs). Recently, it was suggested that minor-class snRNPs are localized to and function in the cytoplasm of vertebrate cells. To test whether U12-type splicing occurs in the cytoplasm of Xenopus oocytes, we performed microinjections of the well-characterized P120 minor-class splicing substrate into the nucleus or into the cytoplasm. Our results demonstrate that accurate splicing of this U12-dependent intron occurs exclusively in the nuclear compartment of the oocyte, where U12 and U6atac snRNPs are primarily localized. We further demonstrate that splicing of both a major-class and a minor-class intron is inhibited after nuclear envelope breakdown during meiosis.     

Dynamic exchanges of RNA interactions leading to catalytic core formation in the U12-dependent spliceosome

Important general insights into the mechanism of pre-mRNA splicing have emerged from studies of the U12-dependent spliceosome. Here, photochemical cross-linking analyses during U12-dependent spliceosome assembly have surprisingly revealed that an upstream 5' exon region is required for establishing two essential catalytic core interactions, U12/U6atac helix Ib and U6atac/5' splice site contacts, but not for U5/5' exon interactions or partial unwinding of U4atac/U6atac. A novel intermediate, representing an alternative pathway for catalytic core formation, is a ternary snRNA complex containing U4atac/U6atac stem II and U12/U6atac helix Ia that forms even without U6atac replacing U11 at the 5' splice site. A powerful oligonucleotide displacement method suggests that the blocked complexes analyzed to deduce the interdependence of these multiple RNA exchanges are authentic intermediates in U12-dependent spliceosome assembly.     

Alternative splicing and bioinformatic analysis of human U12-type introns

U12-type introns exist, albeit rarely, in a variety of multicellular organisms. Splicing of U12 intron-containing precursor mRNAs takes place in the U12-type spliceosome that is distinct from the major U2-type spliceosome. Due to incompatibility of these two spliceosomes, alternative splicing involving a U12-type intron may give rise to a relatively complicated impact on gene expression. We studied alternative U12-type intron splicing in an attempt to gain more mechanistic insights. First, we characterized mutually exclusive exon selection of the human JNK2 gene, which involves an unusual intron possessing the U12-type 5′ splice site and the U2-type 3′ splice site. We demonstrated that the long and evolutionary conserved polypyrimidine tract of this hybrid intron provides important signals for inclusion of its downstream alternative exon. In addition, we examined the effects of single nucleotide polymorphisms in the human WDFY1 U12-type intron on pre-mRNA splicing. These results provide mechanistic implications on splice-site selection of U12-type intron splicing. We finally discuss the potential effects of splicing of a U12-type intron with genetic defects or within a set of genes encoding RNA processing factors on global gene expression.     

Functional interaction of a novel 15.5kD [U4/U6.U5] tri-snRNP protein with the 5' stem-loop of U4 snRNA.

Activation of the spliceosome for splicing catalysis requires the dissociation of U4 snRNA from the U4/U6 snRNA duplex prior to the first step of splicing. We characterize an evolutionarily conserved 15.5 kDa protein of the HeLa [U4/U6.U5] tri-snRNP that binds directly to the 5' stem-loop of U4 snRNA. This protein shares a novel RNA recognition motif with several RNP-associated proteins, which is essential, but not sufficient for RNA binding. The 15.5kD protein binding site on the U4 snRNA consists of an internal purine-rich loop flanked by the stem of the 5' stem-loop and a stem comprising two base pairs. Addition of an RNA oligonucleotide comprising the 5' stem-loop of U4 snRNA (U4SL) to an in vitro splicing reaction blocked the first step of pre-mRNA splicing. Interestingly, spliceosomal C complex formation was inhibited while B complexes accumulated. This indicates that the 15.5kD protein, and/or additional U4 snRNP proteins associated with it, play an important role in the late stage of spliceosome assembly, prior to step I of splicing catalysis. Our finding that the 15.5kD protein also efficiently binds to the 5' stem-loop of U4atac snRNA indicates that it may be shared by the [U4atac/U6atac.U5] tri-snRNP of the minor U12-type spliceosome.     

Evolutionarily emerged G tracts between the polypyrimidine tract and 3′ AG are splicing silencers enriched in genes involved in cancer

Background

The 3′ splice site (SS) at the end of pre-mRNA introns has a consensus sequence (Y)_nNYAG for constitutive splicing of mammalian genes. Deviation from this consensus could change or interrupt the usage of the splice site leading to alternative or aberrant splicing, which could affect normal cell function or even the development of diseases. We have shown that the position “N” can be replaced by a CA-rich RNA element called CaRRE1 to regulate the alternative splicing of a group of genes.

Results

Taking it a step further, we searched the human genome for purine-rich elements between the -3 and -10 positions of the 3′ splice sites of annotated introns. This identified several thousand such 3′SS; more than a thousand of them contain at least one copy of G tract. These sites deviate significantly from the consensus of constitutive splice sites and are highly associated with alterative splicing events, particularly alternative 3′ splice and intron retention. We show by mutagenesis analysis and RNA interference that the G tracts are splicing silencers and a group of the associated exons are controlled by the G tract binding proteins hnRNP H/F. Species comparison of a group of the 3′SS among vertebrates suggests that most (~87%) of the G tracts emerged in ancestors of mammals during evolution. Moreover, the host genes are most significantly associated with cancer.

Conclusion

We call these elements together with CaRRE1 regulatory RNA elements between the Py and 3′AG (REPA). The emergence of REPA in this highly constrained region indicates that this location has been remarkably permissive for the emergence of de novo regulatory RNA elements, even purine-rich motifs, in a large group of mammalian genes during evolution. This evolutionary change controls alternative splicing, likely to diversify proteomes for particular cellular functions.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1143) contains supplementary material, which is available to authorized users.