U12-type spliceosomal introns of Insecta

Most of eukaryotic genes are interrupted by introns that need to be removed from pre-mRNAs before they can perform their function. This is done by complex machinery called spliceosome. Many eukaryotes possess two separate spliceosomal systems that process separate sets of introns. The major (U2) spliceosome removes majority of introns, while minute fraction of intron repertoire is processed by the minor (U12) spliceosome. These two populations of introns are called U2-type and U12-type, respectively. The latter fall into two subtypes based on the terminal dinucleotides. The minor spliceosomal system has been lost independently in some lineages, while in some others few U12-type introns persist. We investigated twenty insect genomes in order to better understand the evolutionary dynamics of U12-type introns. Our work confirms dramatic drop of U12-type introns in Diptera, leaving these genomes just with a handful cases. This is mostly the result of intron deletion, but in a number of dipteral cases, minor type introns were switched to a major type, as well. Insect genes that harbor U12-type introns belong to several functional categories among which proteins binding ions and nucleic acids are enriched and these few categories are also overrepresented among these genes that preserved minor type introns in Diptera.     

Arabidopsis consensus intron sequences

We have analysed 998 Arabidopsis intron sequences in the EMBL database. All Arabidopsis introns to adhere to the :GU...AG: rule with the exception of 1% of introns with :GC at their 5 ends. Virtually all of the introns contained a putative branchpoint sequence (YUNAN) 18 to 60 nt upstream of the 3 splice site. Although a polypyrimidine tract was much less apparent than in vertebrate introns, the most common nucleotide in the region upstream of the 3 splice site was uridine. Consensus sequences for 5 and 3 splice sites and branchpoint sequences for Arabidopsis introns are presented.     

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.     

Comprehensive splice-site analysis using comparative genomics

We have collected over half a million splice sites from five species-Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans and Arabidopsis thaliana-and classified them into four subtypes: U2-type GT-AG and GC-AG and U12-type GT-AG and AT-AC. We have also found new examples of rare splice-site categories, such as U12-type introns without canonical borders, and U2-dependent AT-AC introns. The splice-site sequences and several tools to explore them are available on a public website (SpliceRack). For the U12-type introns, we find several features conserved across species, as well as a clustering of these introns on genes. Using the information content of the splice-site motifs, and the phylogenetic distance between them, we identify: (i) a higher degree of conservation in the exonic portion of the U2-type splice sites in more complex organisms; (ii) conservation of exonic nucleotides for U12-type splice sites; (iii) divergent evolution of C.elegans 3' splice sites (3'ss) and (iv) distinct evolutionary histories of 5' and 3'ss. Our study proves that the identification of broad patterns in naturally-occurring splice sites, through the analysis of genomic datasets, provides mechanistic and evolutionary insights into pre-mRNA splicing.     

Human GC-AG alternative intron isoforms with weak donor sites show enhanced consensus at acceptor exon positions

It has been previously observed that the intrinsically weak variant GC donor sites, in order to be recognized by the U2-type spliceosome, possess strong consensus sequences maximized for base pair formation with U1 and U5/U6 snRNAs. However, variability in signal strength is a fundamental mechanism for splice site selection in alternative splicing. Here we report human alternative GC-AG introns (for the first time from any species), and show that while constitutive GC-AG introns do possess strong signals at their donor sites, a large subset of alternative GC-AG introns possess weak consensus sequences at their donor sites. Surprisingly, this subset of alternative isoforms shows strong consensus at acceptor exon positions 1 and 2. The improved consensus at the acceptor exon can facilitate a strong interaction with U5 snRNA, which tethers the two exons for ligation during the second step of splicing. Further, these isoforms nearly always possess alternative acceptor sites and exhibit particularly weak polypyrimidine tracts characteristic of AG-dependent introns. The acceptor exon nucleotides are part of the consensus required for the U2AF³⁵-mediated recognition of AG in such introns. Such improved consensus at acceptor exons is not found in either normal or alternative GT-AG introns having weak donor sites or weak polypyrimidine tracts. The changes probably reflect mechanisms that allow GC-AG alternative intron isoforms to cope with two conflicting requirements, namely an apparent need for differential splice strength to direct the choice of alternative sites and a need for improved donor signals to compensate for the central mismatch base pair (C-A) in the RNA duplex of U1 snRNA and the pre-mRNA. The other important findings include (i) one in every twenty alternative introns is a GC-AG intron, and (ii) three of every five observed GC-AG introns are alternative isoforms.     

Mutational analysis of the U12-dependent branch site consensus sequence

Highly conserved sequences at the 5′ splice site and branch site of U12-dependent introns are important determinants for splicing by U12-dependent spliceosomes. This study investigates the in vivo splicing phenotypes of mutations in the branch site consensus sequence of the U12-dependent intron F from a human NOL1 (P120) minigene. Intron F contains a fully consensus branch site sequence (UUCCUUAAC). Mutations at each position were analyzed for their effects on U12-dependent splicing in vivo. Mutations at most positions resulted in a significant reduction of correct U12-dependent splicing. Defects observed included increased unspliced RNA levels, the activation of cryptic U2-dependent 5′ and 3′ splice sites, and the activation of cryptic U12-dependent branch/3′ splice sites. A strong correlation was observed between the predicted thermodynamic stability of the branch site: U12 snRNA interaction and correct U12-dependent splicing. The lack of a polypyrimidine tract between the branch site and 3′ splice site of U12-dependent introns and the observed reliance on base-pairing interactions for correct U12-dependent splicing emphasize the importance of RNA/RNA interactions during U12-dependent intron recognition and proper splice site selection.     

Evidence that U2/U6 helix I promotes both catalytic steps of pre-mRNA splicing and rearranges in between these steps

During pre-mRNA splicing, the spliceosome must configure the substrate, catalyze 5′ splice site cleavage, reposition the substrate, and catalyze exon ligation. The highly conserved U2/U6 helix I, which adjoins sequences that define the reactive sites, has been proposed to configure the substrate for 5′ splice site cleavage and promote catalysis. However, a role for this helix at either catalytic step has not been tested rigorously and previous observations question its role at the catalytic steps. Through a comprehensive molecular genetic study of U2/U6 helix I, we found that weakening U2/U6 helix I, but not mutually exclusive structures, compromised splicing of a substrate limited at the catalytic step of 5′ splice site cleavage, providing the first compelling evidence that this helix indeed configures the substrate during 5′ splice site cleavage. Further, mutations that we proved weaken only U2/U6 helix I suppressed a mutation in PRP16, a DEAH-box ATPase required after 5′ splice site cleavage, providing persuasive evidence that helix I is destabilized by Prp16p and suggesting that this structure is unwound between the catalytic steps. Lastly, weakening U2/U6 helix I also compromised splicing of a substrate limited at the catalytic step of exon ligation, providing evidence that U2/U6 helix I reforms and functions during exon ligation. Thus, our data provide evidence for a fundamental and apparently dynamic role for U2/U6 helix I during the catalytic stages of splicing.     

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.     

RNA induces conformational changes in the SF1/U2AF65 splicing factor complex

Spliceosomes assemble on pre-mRNA splice sites through a series of dynamic ribonucleoprotein complexes, yet the nature of the conformational changes remains unclear. Splicing factor 1 (SF1) and U2 auxiliary factor (U2AF⁶⁵) cooperatively recognize the 3′ splice site during the initial stages of pre-mRNA splicing. Here, we used small-angle X-ray scattering to compare the molecular dimensions and ab initio shape restorations of SF1 and U2AF⁶⁵ splicing factors, as well as the SF1/U2AF⁶⁵ complex in the absence and presence of AdML (adenovirus major late) splice site RNAs. The molecular dimensions of the SF1/U2AF⁶⁵/RNA complex substantially contracted by 15 Å in the maximum dimension, relative to the SF1/U2AF⁶⁵ complex in the absence of RNA ligand. In contrast, no detectable changes were observed for the isolated SF1 and U2AF⁶⁵ splicing factors or their individual complexes with RNA, although slight differences in the shapes of their molecular envelopes were apparent. We propose that the conformational changes that are induced by assembly of the SF1/U2AF⁶⁵/RNA complex serve to position the pre-mRNA splice site optimally for subsequent stages of splicing.     

Assessment of a model for intron RNA secondary structure relevant to RNA self-splicing--a review

A widespread class of introns is characterized by a particular RNA secondary structure, based upon four conserved nucleotide sequences. Among such "class I" introns are found the majority of introns in fungal mitochondrial genes and the self-splicing intron of the large ribosomal RNA of several species of Tetrahymena. A model of the RNA secondary structure, which must underlie the self-splicing activity, is here evaluated in the light of data on 16 further introns. The main body or "core structure" of the intron always consists of the base-paired regions P3 to P9 with the associated single-stranded loops, with P2 present also in most cases. Two minority sub-classes of core structure occur, one of which is typical of introns in fungal ribosomal RNA. Introns in which the core structure is close to the 5' splice site all have an internal guide sequence (IGS) which can pair with exon sequences adjacent to the 5' and 3' splice sites to align them precisely, as proposed by Davies et al. [Nature 300 (1982) 719-724]. In these cases, the internal guide model allows us to predict correctly the exact location of splice sites. All other introns probably use other mechanisms of alignment. This analysis provides strong support for the RNA splicing model which we have developed.     

Invasion of protein coding genes by green algal ribosomal group I introns

The spread of group I introns depends on their association with intron-encoded homing endonucleases. Introns that encode functional homing endonuclease genes (HEGs) are highly invasive, whereas introns that only encode the group I ribozyme responsible for self-splicing are generally stably inherited (i.e., vertical inheritance). A number of recent case studies have provided new knowledge on the evolution of group I introns, however, there are still large gaps in understanding of their distribution on the tree of life, and how they have spread into new hosts and genic sites. During a larger phylogenetic survey of chlorophyceaen green algae, we found that 23 isolates contain at least one group I intron in the rbcL chloroplast gene. Structural analyses show that the introns belong to one of two intron lineages, group IA2 intron-HEG (GIY-YIG family) elements inserted after position 462 in the rbcL gene, and group IA1 introns inserted after position 699. The latter intron type sometimes encodes HNH homing endonucleases. The distribution of introns was analyzed on an exon phylogeny and patterns were recovered that are consistent with vertical inheritance and possible horizontal transfer. The rbcL 462 introns are thus far reported only within the Volvocales, Hydrodictyaceae and Bracteacoccus, and closely related isolates of algae differ in the presence of rbcL introns. Phylogenetic analysis of the intron conserved regions indicates that the rbcL699 and rbcL462 introns have distinct evolutionary origins. The rbcL699 introns were likely derived from ribosomal RNA L2449 introns, whereas the rbcL462 introns form a close relationship with psbA introns.     

A novel spliceosome-mediated trans-splicing can change our view on genome complexity of the divergent eukaryote Giardia intestinalis

Although spliceosomal introns are an abundant landmark in eukaryotic genomes, the nuclear genome of the divergent eukaryote Giardia intestinalis, the causative agent of giardiasis, has been considered as “intron-poor” with only five canonical (cis-spliced) introns. However, three research groups (including ours) have independently reported a novel class of spliceosomal introns in the G. intestinalis genome. Three protein-coding genes are split into pieces in the G. intestinalis genome, and each of the partial coding regions was independently transcribed into polyadenylated premature mRNAs (pre-mRNAs). The two pre-mRNAs directly interact with each other by an intermolecular-stem structure formed between their non-coding portions, and are then processed into mature mRNAs by spliceosome-mediated trans-splicing. Here, we summarize the recently published works on split introns (“splintrons”) in the G. intestinalis genome, and then provide our speculation on the functional property of the Giardia spliceosomes based on the putative ratio of splintrons to canonical introns. Finally, we discuss a scenario for the transition from typical GT-AG boundaries to non-typical AT-AC boundaries in a particular splintron of Giardia.     

核小体定位与RNA剪接

根据核小体定位序列和缺失序列的碱基分布特征,应用多样性增量二次判别方法(IDQD)构建模型对这两类序列进行了区分,受试者操作特性曲线下的面积达到了0.958．应用这一模型研究了核小体在人类基因组剪接位点(GT/AG)邻近序列中的分布方式,发现外显子所对应的DNA序列通常倾向参与核小体的形成,并且由它所转录的RNA统计上具有较强的刚性,而剪接位点及其邻近的内含子对应的DNA序列则避免参与核小体的形成,所转录的RNA统计上具有较强的柔性．进一步还发现,DNA序列的核小体定位/缺失和RNA的刚性/柔性具有统计相关性,为从机制上解释为何前体RNA剪接事件与DNA序列中的核小体定位信息有关提供了依据．     

Primordial spliceosomal introns were probably U2-type

The two types of eukaryotic spliceosomal introns, U2 and U12, possess different splice signals and are excised by distinct spliceosomes. The nature of the primordial introns remains uncertain. A comparison of the amino acid distributions at insertion sites of introns that retained their positions throughout eukaryotic evolution with the distributions for human and Arabidopsis thaliana U2 and U12 introns reveals close similarity with U2 but not U12. Thus, the primordial spliceosomal introns were, most likely, U2-type.     

Gene expression profiling of U12-type spliceosome mutant Drosophila reveals widespread changes in metabolic pathways

Background

The U12-type spliceosome is responsible for the removal of a subset of introns from eukaryotic mRNAs. U12-type introns are spliced less efficiently than normal U2-type introns, which suggests a rate-limiting role in gene expression. The Drosophila genome contains about 20 U12-type introns, many of them in essential genes, and the U12-type spliceosome has previously been shown to be essential in the fly.

Methodology/Principal Findings

We have used a Drosophila line with a P-element insertion in U6atac snRNA, an essential component of the U12-type spliceosome, to investigate the impact of U12-type introns on gene expression at the organismal level during fly development. This line exhibits progressive accumulation of unspliced U12-type introns during larval development and the death of larvae at the third instar stage. Surprisingly, microarray and RT-PCR analyses revealed that most genes containing U12-type introns showed only mild perturbations in the splicing of U12-type introns. In contrast, we detected widespread downstream effects on genes that do not contain U12-type introns, with genes related to various metabolic pathways constituting the largest group.

Conclusions/Significance

U12-type intron-containing genes exhibited variable gene-specific responses to the splicing defect, with some genes showing up- or downregulation, while most did not change significantly. The observed residual U12-type splicing activity could be explained with the mutant U6atac allele having a low level of catalytic activity. Detailed analysis of all genes suggested that a defect in the splicing of the U12-type intron of the mitochondrial prohibitin gene may be the primary cause of the various downstream effects detected in the microarray analysis.     

Gene Expression,Intron Density,and Splice Site Strength in <Emphasis Type="Italic">Drosophila</Emphasis> and <Emphasis Type="Italic">Caenorhabditis</Emphasis>

In this paper we investigate the relationships among intron density (number of introns per kilobase of coding sequence), gene expression level, and strength of splicing signals in two species: Drosophila melanogaster and Caenorhabditis elegans. We report a negative correlation between intron density and gene expression levels, opposite to the effect previously observed in human. An increase in splice site strength has been observed in long introns in D. melanogaster. We show this is also true of C. elegans. We also examine the relationship between intron density and splice site strength. There is an increase in splice site strength as the intron structure becomes less dense. This could suggest that introns are not recognized in isolation but could function in a cooperative manner to ensure proper splicing. This effect remains if we control for the effects of alternative splicing on splice site strength. Reviewing Editor: Dr. Nicolas Galtier