Search

The spliceosome catalyzes the removal of introns from pre-messenger RNA (mRNA) and subsequent pairing of exons with remarkable fidelity. Some exons are known to be skipped or included in the mature mRNA in a cell type- or context-dependent manner (cassette exons), thereby contributing to the diversification of the human proteome. Interestingly, splicing is initiated (and sometimes completed) co-transcriptionally. Here, we develop a kinetic mathematical modeling framework to investigate alternative co-transcriptional splicing (CTS) and, specifically, the control of cassette exons’ inclusion. We show that when splicing is co-transcriptional, default splice patterns of exon inclusion are more likely than when splicing is post-transcriptional, and that certain exons are more likely to be regulatable (i.e. cassette exons) than others, based on the exon–intron structure context. For such regulatable exons, transcriptional elongation rates may affect splicing outcomes. Within the CTS paradigm, we examine previously described hypotheses of co-operativity between splice sites of short introns (i.e. ‘intron definition’) or across short exons (i.e. ‘exon definition’), and find that models encoding these faithfully recapitulate observations in the fly and human genomes, respectively.

Alternative splicing (AS), in higher eukaryotes, is one of the mechanisms of post-transcriptional regulation that generate multiple transcripts from the same gene. One particular mode of AS is the skipping event where an exon may be alternatively excluded or constitutively included in the resulting mature mRNA. Both transcript isoforms from this skipping event site, i.e. in which the exon is either included (inclusion isoform) or excluded (skipping isoform), are typically present in one cell, and maintain a subtle balance that is vital to cellular function and dynamics. However, how the prevailing conditions dictate which isoform is expressed and what biological factors might influence the regulation of this process remain areas requiring further exploration. In this study, we have developed a novel computational method, graph-based exon-skipping scanner (GESS), for de novo detection of skipping event sites from raw RNA-seq reads without prior knowledge of gene annotations, as well as for determining the dominant isoform generated from such sites. We have applied our method to publicly available RNA-seq data in GM12878 and K562 cells from the ENCODE consortium and experimentally validated several skipping site predictions by RT-PCR. Furthermore, we integrated other sequencing-based genomic data to investigate the impact of splicing activities, transcription factors (TFs) and epigenetic histone modifications on splicing outcomes. Our computational analysis found that splice sites within the skipping-isoform-dominated group (SIDG) tended to exhibit weaker MaxEntScan-calculated splice site strength around middle, ‘skipping’, exons compared to those in the inclusion-isoform-dominated group (IIDG). We further showed the positional preference pattern of splicing factors, characterized by enrichment in the intronic splice sites immediately bordering middle exons. Finally, our analysis suggested that different epigenetic factors may introduce a variable obstacle in the process of exon–intron boundary establishment leading to skipping events.

Transcripts of the human tumor susceptibility gene 101 (TSG101) are aberrantly spliced in many cancers. A major aberrant splicing event on the TSG101 pre-mRNA involves joining of distant alternative 5′ and 3′ splice sites within exon 2 and exon 9, respectively, resulting in the extensive elimination of the mRNA. The estimated strengths of the alternative splice sites are much lower than those of authentic splice sites. We observed that the equivalent aberrant mRNA could be generated from an intron-less TSG101 gene expressed ectopically in breast cancer cells. Remarkably, we identified a pathway-specific endogenous lariat RNA consisting solely of exonic sequences, predicted to be generated by a re-splicing between exon 2 and exon 9 on the spliced mRNA. Our results provide evidence for a two-step splicing pathway in which the initial constitutive splicing removes all 14 authentic splice sites, thereby bringing the weak alternative splice sites into close proximity. We also demonstrate that aberrant multiple-exon skipping of the fragile histidine triad (FHIT) pre-mRNA in cancer cells occurs via re-splicing of spliced FHIT mRNA. The re-splicing of mature mRNA can potentially generate mutation-independent diversity in cancer transcriptomes. Conversely, a mechanism may exist in normal cells to prevent potentially deleterious mRNA re-splicing events.

In conventionally-expressed eukaryotic genes, transcription start sites (TSSs) can be identified by mapping the mature mRNA 5'-terminal sequence onto the genome. However, this approach is not applicable to genes that undergo pre-mRNA 5'-leader trans-splicing (SL trans-splicing) because the original 5'-segment of the primary transcript is replaced by the spliced leader sequence during the trans-splicing reaction and is discarded. Thus TSS mapping for trans-spliced genes requires different approaches. We describe two such approaches and show that they generate precisely agreeing results for an SL trans-spliced gene encoding the muscle protein troponin I in the ascidian tunicate chordate Ciona intestinalis. One method is based on experimental deletion of trans-splice acceptor sites and the other is based on high-throughput mRNA 5'-RACE sequence analysis of natural RNA populations in order to detect minor transcripts containing the pre-mRNA's original 5'-end. Both methods identified a single major troponin I TSS located ∼460 nt upstream of the trans-splice acceptor site. Further experimental analysis identified a functionally important TATA element 31 nt upstream of the start site. The two methods employed have complementary strengths and are broadly applicable to mapping promoters/TSSs for trans-spliced genes in tunicates and in trans-splicing organisms from other phyla.

Ultraconservation, defined as perfect human-to-rodent sequence identity at least 200-bp long, is a strong indicator of evolutionary and functional importance and has been explored extensively at the genome level. However, it has not been investigated at the transcript level, where such extreme conservation might highlight loci with important post-transcriptional regulatory roles. We present 96 ultraconserved cDNA segments (UCSs), stretches of human mature mRNAs that match identically with orthologous regions in the mouse and rat genomes. UCSs can span multiple exons, a feature we leverage here to elucidate the role of ultraconservation in post-transcriptional regulation. UCS sites are implicated in functions at essentially every post-transcriptional stage: pre-mRNA splicing and degradation through alternative splicing and nonsense-mediated decay (AS-NMD), mature mRNA silencing by miRNA, fast mRNA decay rate and translational repression by upstream AUGs. We also found UCSs to exhibit resistance to formation of RNA secondary structure. These multiple layers of regulation underscore the importance of the UCS-containing genes as key global RNA processing regulators, including members of the serine/arginine-rich protein and heterogeneous nuclear ribonucleoprotein (hnRNP) families of essential splicing regulators. The discovery of UCSs shed new light on the multifaceted, fine-tuned and tight post-transcriptional regulation of gene families as conserved through the majority of the mammalian lineage.

Messenger RNA (mRNA) 3' end formation is a nuclear process through which all eukaryotic primary transcripts are endonucleolytically cleaved and most of them acquire a poly(A) tail. This process, which consists in the recognition of defined poly(A) signals of the pre-mRNAs by a large cleavage/polyadenylation machinery, plays a critical role in gene expression. Indeed, the poly(A) tail of a mature mRNA is essential for its functions, including stability, translocation to the cytoplasm and translation. In addition, this process serves as a bridge in the network connecting the different transcription, capping, splicing and export machineries. It also participates in the quantitative and qualitative regulation of gene expression in a variety of biological processes through the selection of single or alternative poly(A) signals in transcription units. A large number of protein factors associates with this machinery to regulate the efficiency and specificity of this process and to mediate its interaction with other nuclear events. Here, we review the eukaryotic 3' end processing machineries as well as the comprehensive set of regulatory factors and discuss the different molecular mechanisms of 3' end processing regulation by proposing several overlapping models of regulation.

Using site-specific incorporation of the photo-chemical cross-linking reagent 4-thiouridine, we demonstrate the previously unknown association of two proteins with yeast 3' splice sites. One of these is an unidentified approximately 122 kDa protein that cross-links to 3' splice sites during formation of the pre--spliceosome. The other factor is the DExH-box RNA helicase, Prp22p. With substrates functional in the second step of splicing, only very weak cross-linking of Prp22p to intron sequences at the 3' splice site is observed. In contrast, substrates blocked at the second step exhibit strong cross-linking of Prp22 to intron sequences at the 3' splice site, but not to adjacent exon sequences. In vitro reconstitution experiments also show that the association of Prp22p with intron sequences at the 3' splice site is dependent on Prp16p and does not persist when release of mature mRNA from the spliceosome is blocked. Taken together, these results suggest that the 3' splice site of yeast introns is contacted much earlier than previously envisioned by a protein of approximately 120 kDa, and that a transient association of Prp22p with the 3' splice site occurs between the first and second catalytic steps.

Group II introns are self-splicing RNAs and retroelements found in bacteria and lower eukaryotic organelles. During the past several years, they have been uncovered in surprising numbers in bacteria due to the genome sequencing projects; however, most of the newly sequenced introns are not correctly identified. We have initiated an ongoing web site database for mobile group II introns in order to provide correct information on the introns, particularly in bacteria. Information in the web site includes: (1) introductory information on group II introns; (2) detailed information on subfamilies of intron RNA structures and intron-encoded proteins; (3) a listing of identified introns with correct boundaries, RNA secondary structures and other detailed information; and (4) phylogenetic and evolutionary information. The comparative data should facilitate study of the function, spread and evolution of group II introns. The database can be accessed at http:// www.fp.ucalgary.ca/group2introns/.

It has been widely acknowledged that non-coding RNAs are master-regulators of genomic functions. However, the significance of the presence of ncRNA within introns has not received proper attention. ncRNA within introns are commonly produced through the post-splicing process and are specific signals of gene transcription events, impacting many other genes and modulating their expression. This study, along with the following discussion, details the association of thousands of ncRNAs--snoRNA, miRNA, siRNA, piRNA and long ncRNA--within human introns. We propose that such an association between human introns and ncRNAs has a pronounced synergistic effect with important implications for fine-tuning gene expression patterns across the entire genome.

The cell line G403.4.7, isolated as a spontaneous variant of the MPC-11 derived myeloma G403.4, produces a truncated gamma 2b HC protein, but no light chain (LC), and a single gamma 2b specific transcript of 2.4kb. This gamma 2b transcript consists of the VDJ and CH1 exons, the CH1 to Hinge (Hi) intervening sequence (IVS) and HI exon, part of the IVS between the two membrane exons M1 and M2, and most of the membrane 3' untranslated (UT) region. Even though the mature mRNA contains intronic sequences, it is abundant in the cytoplasm. Analysis of the gamma 2b genomic organization reveals that this unusual transcript results in part from two genomic deletions of 2.5kb and 588bp and in part from an altered splicing pattern. This altered splicing pattern is probably a consequence of the sequence alterations resulting from the genomic deletions. Analysis of these events provides some interesting insights into the mechanism of splice site selection and the evolution of introns and exons.

The intronic splicing silencer (ISS) of CFTR exon 9 promotes exclusion of this exon from the mature mRNA. This negative influence has important consequences with regards to human pathologic events, as lack of exon 9 correlates well with the occurrence of monosymptomatic and full forms of CF disease. We have previously shown that the ISS element interacts with members of the SR protein family. In this work, we now provide the identification of SF2/ASF and SRp40 as the specific SR proteins binding to this element and map their precise binding sites in IVS9. We have also performed a functional analysis of the ISS element using a variety of unrelated SR-binding sequences and different splicing systems. Our results suggest that SR proteins mediate CFTR exon 9 exclusion by providing a 'decoy' sequence in the vicinity of its suboptimal donor site. The results of this study give an insight on intron 'exonization' mechanisms and provide useful indications for the development of novel therapeutic strategies aimed at the recovery of exon inclusion.

The alternatively spliced RNA species of tumor suppressor gene p53, containing an additional 96 bases derived from intron 10, is present at approximately 25 to 30% the level of regularly spliced p53 RNA in both normal epidermal and carcinoma cells. The presence of this alternatively spliced RNA in 10T1/2 fibroblast cells, mouse liver and testis suggests that this alternative splicing may be universal. The level of alternatively spliced p53 RNA was increased coordinately with that of regularly spliced p53 in 10T1/2 cells in response to epidermal growth factor. Immunoprecipitation analysis of epidermal cells using monoclonal antibodies which recognize different epitopes of p53 suggested that distinct p53 proteins may be translated from both RNA species. Considering previous observations on the potential importance of carboxyl terminal sequences in p53 function, knowledge of the ubiquitous presence of alternatively spliced p53 is important for future studies of p53 function in normal cells and in oncogenesis.

We have isolated and characterised a differentially-regulated gene family in the protozoan parasite Leishmania major. The family contains 5 genes linked within a 10Kb region of the genome: three of the genes are closely related in DNA sequence, the other two have only limited homology. Post-transcriptional control of the differential expression pattern is suggested by detection of precursor RNA molecules containing intergenic sequences and evidence that mature mRNA molecules contain a 35nt spliced leader sequence at their 5' ends. These features support a model of polycistronic transcription in which the stability and differential processing of precursor RNA molecules determine the steady state levels of mature mRNA. We have identified several DNA sequence motifs within the gene family that have potential roles in differential processing and/or RNA stability: an alternative 5' splice acceptor site for trans-splicing; a putative polyadenylation site; and a region of potential secondary structure within 3' flanking sequences. The 3' sequence elements are conserved in those genes that share the same pattern of differential regulation. To our knowledge, this is the first example of coordinated differential-regulation of a non-identical gene cluster in Leishmania.

The Yeast Intron DataBase (YIDB) contains currently available information about all introns encoded in the nuclear and mitochondrial genomes of the yeast Saccharomyces cerevisiae. Introns are divided according to their mechanism of excision: group I and group II introns, pre-mRNA introns, tRNA introns and the HAC1 intron. Information about the host genome, the type of RNA in which they are inserted and their primary structure are provided together with references. For nuclear pre-mRNA introns, transcription frequencies, as determined by microarray experiments, have also been included. This updated database is accessible at: http://www.embl-heidelberg. de/ExternalInfo/seraphin/yidb.html

We have previously cloned a gene for a developmentally regulated transport protein from the trypanosomatid protozoan Leishmania enriettii. We demonstrate here that this transporter is encoded by a single family of tandemly clustered genes containing approximately 8 copies of the 3.6 kilobase repeat unit. Transcriptional mapping defines a contiguous 3.3 kilobase region of the repeat unit that encodes the mRNA. The 5' end of the mature mRNA contains the spliced leader or mini-exon previously identified in kinetoplastid protozoa, while the 3' ends of the mRNA are heterogeneous in sequence and in location of the polyadenylation site. We have identified genomic restriction fragments that flank the tandem repeat on the 5' and 3' sides and which may be linked to sequences required for expression of the gene family. Other species of Leishmania also contain sequences that hybridize to the cloned L. enriettii gene at high stringency.

The secondary structure of the 9S RNA precursor to ribosomal 5S RNA in Escherichia coli has been determined using chemical reagents and ribonucleases in combination with a reverse transcription procedure. The 9S RNA precursor was generated in vitro by T7 RNA polymerase, and the rrnB operon terminator, T1, was able to terminate the in vitro transcript. The secondary structure model exhibits three structural domains corresponding to a 5' region, a mature region and a terminator region. The mature domain is structurally identical to 5S RNA, and the ribosomal proteins L18 and L25 are able to bind to the precursor. The processing endoribonuclease RNase E cleaves between the structural domains. Moreover, an intramolecular refolding of the nascent transcript must take place if the current view of RNase III processing stems is correct.

Thermodynamically stable secondary structures have been found for adjacent splice junctions of 17 introns from 10 eukaryotic and viral mRNA precursors. Complementary regions in these structures position the exons to be ligated in a favourable orientation for splicing. The possible role of these structures in a general splicing mechanism in described.

Rabbit polyribosomal globin messenger RNA (mRNA) and messenger ribonucleoprotein (mRNP) were labelled at the 3′ poly(A) tail to high specific activity with T4 RNA ligase and [5′-32P]pCp without consequent loss of functional activity. Labelled message was translated in both micrococcal nuclease treated and untreated rabbit reticulocyte lysates, as shown by the formation of labelled polyribosomes. The utilisation of labelled messenger was abolished by T2 toxin or sodium fluoride which are known to inhibit protein synthesis.

The alternative RNA processing pathways in human calcitonin gene (CALC-I gene) expression were investigated using steady state RNA isolated from human medullary thyroid carcinoma (MTC) and from a culture line derived from this tumor. On Northern blots the mature 1.0 kilobases (Kb) calcitonin (CT) - and 1.1 Kb calcitonin gene-related peptide (CGRP) mRNAs were detected with CALCI gene specific probes as well as high molecular weight poly (A) containing RNAs of 2.1, 2.3, 3.3, 4.2, 5.0 and 5.7 Kb. The 5.7 Kb RNA was identified as the poly(A) tailed primary transcript containing sequences corresponding to all 6 exons and 5 introns of the CALC-I gene. From the composition of the other RNAs the splicing order of the different introns could be deduced. The results suggest the following model. First all introns not involved in alternative processing (introns 1, 2 and 5) are spliced from the 5.7 Kb RNA in rapid successive reactions yielding a 3.3 Kb RNA, which accumulates. From this 3.3 Kb RNA, the last common intermediate in the alternative processing pathway, CT mRNA is formed by splicing of intron 3 and poly(A) addition at exon 4, in this order or the reverse order via 2.3 Kb or 2.1 Kb RNA intermediates respectively. Alternatively, the whole intron 3-exon 4-intron 4 region is spliced from the 3.3 Kb RNA yielding CGRP mRNA. The temporal sequence of poly(A) addition at exons 4 and 6 may relate to the observed structural differences between the poly(A) addition signals at these sites. The ratio of CT- to CGRP mRNA may relate also to the differences in the primary structures of the intron 3- and intron 4 splice acceptor sites.

We have used a modification of the Berk-Sharp technique to determine that the 5' termini of the mouse 15 S beta-globin precursor and the mature mRNA have identical map coordinates. The modification involves the use of 5' (or 3') terminally labeled probes; it allows the detection of the precursor in the presence of excess mature mRNA.

The initiation of transcription of the yeast mitochondrial genes coding for subunit I of cytochrome c oxidase (COX1) and for tRNA1Thr has been examined. COX1 messenger RNA synthesis is initiated in a conserved nonanucleotide sequence (ATATAAGTA) which we have previously found immediately upstream of ribosomal RNA genes at positions at which RNA synthesis starts. The 5'-end of the precursor of tRNA1Thr is located in a variant nonanucleotide motif (TTATAAGTA), which may be characteristic for tRNA genes. Using a partially purified fraction of mtRNA polymerase, we demonstrate that RNA synthesis is precisely initiated in vitro in nonanucleotide sequences preceding both ribosomal RNA-, tRNA- and messenger RNA-encoding genes and origins of replication.

Two independent systems and several analytical procedures have been used to establish that isolated mammalian nuclei selectively transport mature RNA polymerase I and II products. Murine myeloma nuclei retain physiologic restriction in our transport assay as assessed by the transport of the immunoglobulin kappa light chain mRNA and 18S and 28S rRNAs. Nearly 50% of the total kappa exons are transported as structurally intact mature mRNA molecules while less than 8% of either pulse-labeled or steady state kappa intron sequences are detected in the transported fraction. Ribosomal external transcribed spacer sequences also are absent in transported RNA. Release of cytoplasmic RNA from the outer nuclear membrane during the transport assay accounts for less than 10% of transported RNA. Nuclei isolated from adenovirus-infected HeLa cells at 20 hours post infection retain cellular actin mRNA and transport viral poly A+RNA. Ribosomal RNA is transported from infected nuclei although at a reduced rate compared to transport from mock-infected nuclei. Inhibition of transport of host mRNA is paralleled by the absence of pulse-labeled actin mRNA in the cytoplasm of infected cells. The implications of our transport data in relationship to intranuclear RNA trafficking are discussed.

We have analyzed SV40 specific RNA sequences which are associated with purified viral chromatin. Three discrete size classes of large unprocessed transcripts were identified: 4.8, 3.6 and 3.0 kb. Furthermore, two polyadenylated species of 2.6 and 1.6 kb were found, the latter of which representing mature 16 S mRNA. S1 analysis revealed both the presence of spliced and unspliced RNA sequences in transcriptionally active SV40 chromatin.

Using a cDNA for an abundant Trypanosoma cruzi mRNA as probe, we have cloned and sequenced a gene which is organized in at least 20 nearly perfect tandem repeats of 940 base pairs. The 5' end of the mRNA has been sequenced by primer extension and found to contain a 35 nucleotide mini-exon (or spliced-leader) sequence that is ubiquitous in trypanosome mRNAs. This sequence, however, is not present in the tandem genomic repeats which encode the exon containing the major portion of the mRNA. Previous studies have shown that the 35-nucleotide sequence is encoded by a separate tandem gene family. One model to explain the formation of a segmented mRNA invokes priming of transcription by a small RNA which contains the leader sequence at its 5' end. However, northern blot analysis of total trypanosome RNA reveals a ladder of molecules larger than the mature mRNA, which appear to be faithful multimeric copies of the tandem gene. The discrete sizes of these RNAs correspond to those expected for partially processed precursors. These observations lend credence to the possibility of an alternative model where segmented mRNAs are generated by inter-molecular splicing.

RNA expression in plant mitochondria implies a large number of post-transcriptional events in which transcript processing and stabilization are essential. In this study, we analyzed the function of the Arabidopsis mitochondrial stability factor 2 gene (MTSF2) and show that the encoded pentatricopeptide repeat protein is essential for the accumulation of stable nad1 mRNA. The production of mature nad1 requires the assembly of three independent RNA precursors via two trans-splicing reactions. Genetic analyses revealed that the lack of nad1 in mtsf2 mutants results from the specific destabilization of the nad1 exons 2-3 precursor transcript. We further demonstrated that MTSF2 binds to its 3΄ extremity with high affinity, suggesting a protective action by blocking exoribonuclease progression. By defining the 3΄ end of nad1 exons 2-3 precursor, MTSF2 concomitantly determines the 3΄ extremity of the first half of the trans-intron found at the end of the transcript. Therefore, binding of the MTSF2 protein to nad1 exons 2-3 precursor evolved both to stabilize the transcript and to define a 3΄ extremity compatible with the trans-splicing reaction needed to reconstitute mature nad1. We thus reveal that the range of transcripts stabilized by association with protective protein on their 3΄ end concerns also mitochondrial precursor transcripts.

Transcript splicing in plant mitochondria involves numerous nucleus-encoded factors, most of which are of eukaryotic origin. Some of these belong to protein families initially characterised to perform unrelated functions. The RAD52-like ODB1 protein has been reported to have roles in homologous recombination-dependent DNA repair in the nuclear and mitochondrial compartments in Arabidopsis thaliana. We show that it is additionally involved in splicing and facilitates the excision of two cis-spliced group II introns, nad1 intron 2 and nad2 intron 1, in Arabidopsis mitochondria. odb1 mutants lacking detectable amounts of ODB1 protein over-accumulated incompletely spliced nad1 and nad2 transcripts. The two ODB1-dependent introns were both found to splice via first-step hydrolysis and to be released as linear or circular molecules instead of lariats. Our systematic analysis of the structures of excised introns in Arabidopsis mitochondria revealed several other hydrolytically spliced group II introns in addition to nad1 intron 2 and nad2 intron 1, indicating that ODB1 is not a general determinant of the hydrolytic splicing pathway.