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.     

Intron sequence and structure requirements for tRNA splicing in Saccharomyces cerevisiae

Predicted single-stranded structure at the 3' splice site is a conserved feature among intervening sequences (IVSs) in eukaryotic nuclear tRNA precursors. The role of 3' splice site structure in splicing was examined through hexanucleotide insertions at a central intron position in the Saccharomyces cerevisiae tRNA gene. These insertions were designed to alter the structure at the splice site without changing its sequence. Endonuclease cleavage of pre-tRNA substrates was then measured in vitro, and suppressor activity was examined in vivo. A precursor with fully double-stranded structure at the 3' splice site was not cleaved by endonuclease. The introduction of one unpaired nucleotide at the 3' splice site was sufficient to restore cleavage, although at a reduced rate. We have also observed that guanosine at the antepenultimate position provides a second consensus feature among IVSs in tRNA precursors. Point mutations at this position were found to affect splicing although there was no specific requirement for guanosine. These and previous results suggest that elements of secondary and/or tertiary structure at the 3' end of IVSs are primary determinants in pre-tRNA splice site utilization whereas specific sequence requirements are limited.     

Phylogenetic distribution of intron positions in alpha-amylase genes of bilateria suggests numerous gains and losses

Most eukaryotes have at least some genes interrupted by introns. While it is well accepted that introns were already present at moderate density in the last eukaryote common ancestor, the conspicuous diversity of intron density among genomes suggests a complex evolutionary history, with marked differences between phyla. The question of the rates of intron gains and loss in the course of evolution and factors influencing them remains controversial. We have investigated a single gene family, alpha-amylase, in 55 species covering a variety of animal phyla. Comparison of intron positions across phyla suggests a complex history, with a likely ancestral intronless gene undergoing frequent intron loss and gain, leading to extant intron/exon structures that are highly variable, even among species from the same phylum. Because introns are known to play no regulatory role in this gene and there is no alternative splicing, the structural differences may be interpreted more easily: intron positions, sizes, losses or gains may be more likely related to factors linked to splicing mechanisms and requirements, and to recognition of introns and exons, or to more extrinsic factors, such as life cycle and population size. We have shown that intron losses outnumbered gains in recent periods, but that "resets" of intron positions occurred at the origin of several phyla, including vertebrates. Rates of gain and loss appear to be positively correlated. No phase preference was found. We also found evidence for parallel gains and for intron sliding. Presence of introns at given positions was correlated to a strong protosplice consensus sequence AG/G, which was much weaker in the absence of intron. In contrast, recent intron insertions were not associated with a specific sequence. In animal Amy genes, population size and generation time seem to have played only minor roles in shaping gene structures.     

Identification of functional single nucleotide polymorphisms in the branchpoint site

Background

The human genome contains millions of single nucleotide polymorphisms (SNPs); many of these SNPs are intronic and have unknown functional significance. SNPs occurring within intron branchpoint sites, especially at the adenine (A), would presumably affect splicing; however, this has not been systematically studied. We employed a splicing prediction tool to identify human intron branchpoint sites and screened dbSNP for identifying SNPs located in the predicted sites to generate a genome-wide branchpoint site SNP database.

Results

We identified 600 SNPs located within branchpoint sites; among which, 216 showed a change in A. After scoring the SNPs by counting the As in the ±?10 nucleotide region, only four SNPs were identified without additional As (rs13296170, rs12769205, rs75434223, and rs67785924). Using minigene constructs, we examined the effects of these SNPs on splicing. The three SNPs (rs13296170, rs12769205, and rs75434223) with nucleotide substitution at the A position resulted in abnormal splicing (exon skipping and/or intron inclusion). However, rs67785924, a 5-bp deletion that abolished the branchpoint A nucleotide, exhibited normal RNA splicing pattern, presumably using two of the downstream As as alternative branchpoints. The influence of additional As on splicing was further confirmed by studying rs2733532, which contains three additional As in the ±?10 nucleotide region.

Conclusions

We generated a high-confidence genome-wide branchpoint site SNP database, experimentally verified the importance of A in the branchpoint, and suggested that other nearby As can protect branchpoint A substitution from abnormal splicing.

     

Yeast mRNA splicing in vitro

Synthetic actin and CYH2 pre-mRNAs containing a single intron are accurately spliced in a soluble whole cell extract of yeast. Splicing in vitro requires ATP. The excised intron is released as a lariat in which an RNA branch connects the 5' end of the molecule to the last A in the "intron conserved sequence" UACUAAC. Two other discrete RNA species produced during splicing in vitro may represent reaction intermediates: free, linear exon 1 and a form of the intron lariat extending beyond the 3' splice site to include exon 2. Both lariat forms correspond to molecules previously shown to be produced during yeast pre-mRNA splicing in vivo.     

Evidence for the Ancestral Origin of Group I Introns in the SSUrDNA of Naegleria spp.

ABSTRACT The sequence variation within the group I intron in five Naegleria spp. was studied and compared with the sequence variation within the flanking small subunit ribosomal DNA. Considerable sequence divergence was observed in the introns as well as in the rDNA. In the intron deletions and insertions are only detected in the sequence contributing to the secondary structure, not in the open reading frame. Most of the sequence variation is detected in the unpaired loops. In the case of nucleotide substitution in helices, compensating base pair changes were observed. The sequence variation does not induce variation in the secondary structure model. The phylogenetic tree based on the intron sequences is similar to the tree based on the flanking rDNA sequences. This observation indicates that the intron might have been acquired at an early stage in evolution, and lost in the majority of Naegleria spp.     

Intron open reading frames as mobile elements and evolution of a group I intron

Group I introns are proposed to have become mobile following the acquisition of open reading frames (ORFs) that encode highly specific DNA endonucleases. This proposal implies that intron ORFs could behave as autonomously mobile entities. This was supported by abundant circumstantial evidence but no experiment of ORF transfer from an ORF- containing intron to its ORF-less counterpart has been described. In this paper we present such experiments, which demonstrate the efficient mobility of the mitochondrial nad1-i4-orf1 between two Podospora strains. The homing of this mobile ORF was accompanied by a bidirectional co-conversion that did not systematically involve the whole intron sequence. Orf1 acquisition would be the most recent step in the evolution of the nad1-i4 intron, which has resulted in many strains of Podospora having an intron with two ORFs (biorfic) and four splicing pathways. We show that two of the splicing events that operate in this biorfic intron, as evidenced by PCR experiments, are generated by a 5'-alternative splice site, which is most probably a remnant of the monoorfic ancestral form of the intron. We propose a sequential evolution model that is consistent with the four organizations of the corresponding nad1 locus that we found among various species of the Pyrenomycete family; these organizations consist of no intron, an intron alone, a monoorfic intron, and a biorfic intron.      

Mutational analysis of conserved nucleotides in a self-splicing group I intron.

We have constructed all single base substitutions in almost all of the highly conserved residues of the Tetrahymena self-splicing intron. Mutation of highly conserved residues almost invariably leads to loss of enzymatic activity. In many cases, activity could be regained by making additional mutations that restored predicted base-pairings; these second site suppressors in general confirm the secondary structure derived from phylogenetic data. At several positions, our suppression data can be most readily explained by assuming non-Watson-Crick base-pairings. In addition to the requirements imposed by the secondary structure, the sequence of the intron is constrained by "negative interactions", the exclusion of particular nucleotide sequences that would form undesirable secondary structures. A comparison of genetic and phylogenetic data suggests sites that may be involved in tertiary structural interactions.     

裸子植物5S rRNA基因序列变异及二级结构特征

在高等植物中,5SrRNA基因一级结构是高度保守的,二级结构也相当一致。通过比较18种裸子植物5SrRNA基因序列和二级结构变异,发现55％的核苷酸位点是可变的,这种变异有68％发生在干区（双链区）,其中一些变异,如双链的互补性核苷酸替代,GU配对等能够维系5SrRNA二级结构的稳定性。环区相对保守,这与5SrRNA三级结构折叠或在转录翻译过程中蛋白质、RNA的结合相关。另外,首次报道了松属环E区核苷酸的变异性,这可能与其他区域的变异一样,是假基因造成的结果。5SrRNA基因信息可反映大分类群的系统进化关系,但由于基因长度短,信息量小,其在近缘种系统分类的应用受到限制。     

Molecular evolution of a plastid tandem repeat locus in an orchid lineage

The molecular evolution of a chloroplast minisatellite locus in the Anacamptis palustris (Orchidaceae) lineage and haplotype variation in two Italian A. palustris populations were investigated. A phylogenetic analyses of the chloroplast tRNA(LEU) intron, where the minisatellite locus is located, revealed that a deletion in the ancestor of the A. palustris lineage led to the formation of two noncontiguous, complementary sequence motifs. We propose a model to explain the initial formation of the minisatellite repeat motif, starting with the two noncontiguous, complementary sequence motifs. A survey of minisatellite variation in four species of the A. palustris lineage revealed several haplotypes that differed not only in repeat number, but also in repeat organization. A haplotype network suggests that three different minisatellite loci evolved independently at the same position in the tRNA(LEU) intron. A secondary structure model revealed that the A. palustris minisatellite repeat forms a stem region of the tRNA(LEU) intron, which allows its notable expansion without negatively affecting splicing. Minisatellite variation was high in the two examined A. palustris populations where 20 haplotypes were detected, whereas no length variation was detected in a neighboring poly (A) microsatellite locus. We estimated a chloroplast minisatellite mutation rate of 3.2 x 10(-3) mutations per generation. Southern blot analyses did not find evidence for chloroplast heteroplasmy. Based on the analysis of the largest known, extant A. palustris population, a stepwise mutation model (SMM) was inferred.     

Remarkable sequence conservation of the last intron in the PKD1 gene

The last intron of the PKD1 gene (intron 45) was found to have exceptionally high sequence conservation across four mammalian species: human, mouse, rat, and dog. This conservation did not extend to the comparable intron in pufferfish. Pairwise comparisons for intron 45 showed 91% identity (human vs. dog) to 100% identity (mouse vs. rat) for an average for all four species of 94% identity. In contrast, introns 43 and 44 of the PKD1 gene had average pairwise identities of 57% and 54%, and exons 43, 44, and 45 and the coding region of exon 46 had average pairwise identities of 80%, 84%, 82%, and 80%. Intron 45 is 90 to 95 bp in length, with the major region of sequence divergence being in a central 4-bp to 9-bp variable region. RNA secondary structure analysis of intron 45 predicts a branching stem-loop structure in which the central variable region lies in one loop and the putative branch point sequence lies in another loop, suggesting that the intron adopts a specific stem-loop structure that may be important for its removal. Although intron 45 appears to conform to the class of small, G-triplet-containing introns that are spliced by a mechanism utilizing intron definition, its high sequence conservation may be a reflection of constraints imposed by a unique mechanism that coordinates splicing of this last PKD1 intron with polyadenylation.     

Molecular consequences of specific intron mutations on yeast mRNA splicing in vivo and in vitro

We have altered the TACTAAC sequence in the yeast CYH2m gene intron to TACTACC. This mutation changes the nucleotide at the normal position of the branch in intron RNA lariats produced during pre-mRNA splicing, and it prevents splicing in vivo. In a yeast pre-mRNA splicing system, CYH2m pre-mRNA carrying the TACTACC mutation is not specifically cut or rearranged in any way. Substitution of an A for the first G of the CYH2m intron, converting the highly conserved GTATGT 5' splice site sequence to ATATGT, also blocks intron excision in vivo and in vitro: pre-mRNA carrying this mutation was still cut normally at the mutant 5' splice site in vitro, to give authentic exon 1 and an intron-exon 2 lariat RNA with an A-A 2'-5' phosphodiester linkage at the branch point. This lariat RNA is a dead-end product. The subsequent cleavage at the 3' splice site is therefore sensitive to the sequence of the 5' end of the intron attached at the branch point.     

A self-splicing group I intron in the nuclear pre-rRNA of the green alga, Ankistrodesmus stipitatus.

The nuclear small subunit ribosomal RNA gene of the unicellular green alga Ankistrodesmus stipitatus contains a group I intron, the first of its kind to be found in the nucleus of a member of the plant kingdom. The intron RNA closely resembles the group I intron found in the large subunit rRNA precursor of Tetrahymena thermophila, differing by only eight nucleotides of 48 in the catalytic core and having the same peripheral secondary structure elements. The Ankistrodesmus RNA self-splices in vitro, yielding the typical group I intron splicing intermediates and products. Unlike the Tetrahymena intron, however, splicing is accelerated by high concentrations of monovalent cations and is rate-limited by the exon ligation step. This system provides an opportunity to understand how limited changes in intron sequence and structure alter the properties of an RNA catalytic center.     

Role of the 5' hairpin structure in the splicing accuracy of the fourth intron of the yeast cob-box gene

The splicing mechanism of the maturase-coding introns is poorly understood. We have systematically examined the phenotypes of a large number of revertants from the mitochondrial mutation G2457. This mutation results from a single base change near the 5' splicing site. We show here that this base change does not completely block the splicing of the intron but rather affects the specificity of the splicing process. We examine four classes of revertants which allow us to characterize the crucial role of a stem and loop structure in the accuracy of the intron excision process. An unexpected class of revertant suggests that other elements are involved in this mechanism. Reversion of G2457 can also occur via the excision in the mitochondrial genome of the intron coding sequence. These results are discussed in relation to the possible role fulfilled by the maturase in the control of intron splicing.     

The length of the downstream exon and the substitution of specific sequences affect pre-mRNA splicing in vitro.

We have shown previously that truncation of the human beta-globin pre-mRNA in the second exon, 14 nucleotides downstream from the 3' splice site, leads to inhibition of splicing but not cleavage at the 5' splice site. We now show that several nonglobin sequences substituted at this site can restore splicing and that the efficiency of splicing depends on the length of the second (downstream) exon and not a specific sequence. Deletions in the first exon have no effect on the efficiency of in vitro splicing. Surprisingly, an intron fragment from the 5' region of the human or rabbit beta-globin intron 2, when placed 14 nucleotides downstream from the 3' splice site, inhibited all the steps in splicing beginning with cleavage at the 5' splice site. This result suggests that the intron 2 fragment carries a "poison" sequence that can inhibit the splicing of an upstream intron.     

Molecular genetics of group I introns: RNA structures and protein factors required for splicing--a review

In vivo and in vitro genetic techniques have been widely used to investigate the structure-function relationships and requirements for splicing of group-I introns. Analyses of group-I introns from extremely diverse genetic systems, including fungal mitochondria, protozoan nuclei, and bacteriophages, have yielded results which are complementary and highly consistent. In vivo genetic studies of fungal mitochondrial systems have served to identify cis-acting sequences within mitochondrial introns, and trans-acting protein products of mitochondrial and nuclear genes which are important for splicing, and to show that some mitochondrial introns are mobile genetic elements. In vitro genetic studies of the self-splicing intron within the Tetrahymena thermophila nuclear large ribosomal RNA precursor (Tetrahymena LSU intron) have been used to examine essential and nonessential RNA sequences and structures in RNA-catalyzed splicing. In vivo and in vitro genetic analysis of the intron within the bacteriophage T4 td gene has permitted the detailed examination of mutant phenotypes by analyzing splicing in vivo and self-splicing in vitro. The genetic studies combined with phylogenetic analysis of intron structure based on comparative nucleotide sequence data [Cech 73 (1988) 259-271] and with biochemical data obtained from in vitro splicing experiments have resulted in significant advances in understanding the biology and chemistry of group-I introns.     

Three variant introns of the same general class in the mitochondrial gene for cytochrome oxidase subunit 1 in Aspergillus nidulans

The oxiA gene of Aspergillus nidulans, coding for cytochrome oxidase subunit 1, is shown by DNA sequencing to contain three introns. An AUG start codon is not present at the beginning of the sequence, suggesting that either another codon, possibly the four base codon AUGA, is used for initiation or there is a further short intron between the true start codon and the beginning of the recognisable coding region. The second and third introns have long open reading frames, which could code for maturase proteins. The lack of conservation of amino acid sequence in the putative region of proteolytic cleavage for maturase formation suggests that the first conserved decapeptide may act as the recognition signal for protein processing. The third intron is remarkably (70%) homologous to the second intron of the cytochrome oxidase subunit 1 gene of Schizosaccharomyces pombe and both are located in exactly the same position. The third Aspergillus intron has an in-frame insertion of a 37-bp GC-rich DNA sequence which is now flanked by a 5-bp repeat, a well-known feature of transposable elements. All three introns in the oxiA gene have a 'core' RNA secondary structure found in a class of introns fitting the RNA splicing model of Davies et al. (1982). This core RNA structure may play a catalytic as well as a structural role in intron splicing. A sequence within the intron could act as a guide to align the splice sites of two of the introns in accordance with the model of Davies et al.