Discovery of a group I intron in the SSU rDNA of Ploeotia costata (Euglenozoa)

The gene coding for the small ribosomal subunit RNA of Ploeotia costata contains an actively splicing group I intron (Pco.S516) which is unique among euglenozoans. Secondary structure predictions indicate that paired segments P1-P10 as well as several conserved elements typical of group I introns and of subclass IC1 in particular are present. Phylogenetic analyses of SSU rDNA sequences demonstrate a well-supported placement of Ploeotia costata within the Euglenozoa; whereas, analyses of intron data sets uncover a close phylogenetic relation of Pco.S516 to S-516 introns from Acanthamoeba, Aureoumbra lagunensis (Stramenopila) and red algae of the order Bangiales. Discrepancies between SSU rDNA and intron phylogenies suggest horizontal spread of the group I intron. Monophyly of IC1 516 introns from Ploeotia costata, A. lagunensis and rhodophytes is supported by a unique secondary structure element: helix P5b possesses an insertion of 19 nt length with a highly conserved tetraloop which is supposed to take part in tertiary interactions. Neither functional nor degenerated ORFs coding for homing endonucleases can be identified in Pco.S516. Nevertheless, degenerated ORFs with His-Cys box motifs in closely related intron sequences indicate that homing may have occurred during evolution of the investigated intron group.     

Distinct sites of phosphorothioate substitution interfere with folding and splicing of the Anabaena group I intron

Although the active site of group I introns is phylogenetically conserved, subclasses of introns have evolved different mechanisms of stabilizing the catalytic core. Large introns contain weakly conserved 'peripheral' domains that buttress the core through predicted interhelical contacts, while smaller introns use loop-helix interactions for stability. In all cases, specific and non-specific magnesium ion binding accompanies folding into the active structure. Whether similar RNA-RNA and RNA-magnesium ion contacts play related functional roles in different introns is not clear, particularly since it can be difficult to distinguish interactions directly involved in catalysis from those important for RNA folding. Using phosphorothioate interference with RNA activity and structure in the small (249 nt) group I intron from Anabaena, we used two independent assays to detect backbone phosphates important for catalysis and those involved in intron folding. Comparison of the interference sites identified in each assay shows that positions affecting catalysis cluster primarily in the conserved core of the intron, consistent with conservation of functionally important phosphates, many of which are magnesium ion binding sites, in diverse group I introns, including those from Azoarcus and Tetrahymena. However, unique sites of folding interference located outside the catalytic core imply that different group I introns, even within the same subclass, use distinct sets of tertiary interactions to stabilize the structure of the catalytic core.     

Activity and thermostability of the small self-splicing group I intron in the pre-tRNA(lle) of the purple bacterium Azoarcus.

The 205-nt group I intron located in the pre-tRNA(lle) from the bacterium Azoarcus sp.BH72 is the smallest self-splicing group I intron identified to date. Comparative sequence analysis has placed this intron and the Anabaena pre-tRNA(Leu) intron into the same subgroup, IC3; we now compare their activity and stability. Unlike the Anabaena intron, the Azoarcus intron has two transitions in the kinetics of the first step of splicing. The faster transition occurs with a larger k(cat)/K(m) than that of the Anabaena or other group I introns, due to a rapid K(cat) (5 min(-1) at 32 degrees C) and a low K(m) for guanosine (17 microM). The excised intron circularizes by releasing a trinucleotide from the 5' end of the intron, another property unlike the Anabaena intron. Although it is smaller in size, the Azoarcus intron retains activity at higher temperatures, higher concentrations of urea, and higher pH than the Anabaena intron. Melting curves show that tertiary structure is disrupted at a lower temperature in the Anabaena intron. Some structural features that may explain the unusual stability of the Azoarcus intron include a G-C rich secondary structure and the presence of two 11-nt motifs, which are known to interact strongly with GAAA loops in group I and group II introns. The disruption of one of these interactions by substituting the Anabaena structural element in fact lowered the thermal stability of the Azoarcus intron. Thus, even superficially similar group I introns from the same structural subgroup can differ significantly in activity and stability.     

Structural metals in the group I intron: a ribozyme with a multiple metal ion core

Metal ions play key roles in the folding and function for many structured RNAs, including group I introns. We determined the X-ray crystal structure of the Azoarcus bacterial group I intron in complex with its 5' and 3' exons. In addition to 222 nucleotides of RNA, the model includes 18 Mg(2+) and K(+) ions. Five of the metals bind within 12 A of the scissile phosphate and coordinate the majority of the oxygen atoms biochemically implicated in conserved metal-RNA interactions. The metals are buried deep within the structure and form a multiple metal ion core that is critical to group I intron structure and function. Eight metal ions bind in other conserved regions of the intron structure, and the remaining five interact with peripheral structural elements. Each of the 18 metals mediates tertiary interactions, facilitates local bends in the sugar-phosphate backbone or binds in the major groove of helices. The group I intron has a rich history of biochemical efforts aimed to identify RNA-metal ion interactions. The structural data are correlated to the biochemical results to further understand the role of metal ions in group I intron structure and function.     

Leaving group stabilization by metal ion coordination and hydrogen bond donation is an evolutionarily conserved feature of group I introns.

To understand the behavior of group I introns on a biologically fundamental level, we must distinguish those traits that arise as the products of natural selection (selected traits) from those that arise as the products of neutral drift (non-selected traits). In practice, this distinction relies on comparing the similarities and differences among widely divergent introns to identify conserved traits. Here we address whether the strategies used by the eukaryotic group I intron from the Tetrahymena ciliate to stabilize the leaving group during splicing are maintained in the group I intron from the widely divergent Azoarcus bacterium. A substrate analogue containing a 3'-phosphorothiolate linkage, in which a sulfur atom replaces the bridging 3'-oxygen atom of the scissile phosphate, reacts 20-fold slower in the Azoarcus reaction than the corresponding unmodified substrate in the presence of Mg(II) as the only divalent cation. However, Mn(II) relieves this negative effect such that the 3'-S-P bond cleaves 21-fold faster than does the 3'O-P bond. Other thiophilic divalent metal ions such as Co(II), Cd(II), and Zn(II) similarly support cleavage of the S-P bond. These results indicate that a metal ion directly coordinates to the leaving group in the transition state of the Azoarcus ribozyme reaction. Additionally, the 3'-sulfur substitution eliminates the approximately 10(3)-fold contribution of the adjacent 2'-OH to transition state stabilization. Considering that sulfur accepts hydrogen bonds weakly compared to oxygen, this result suggests that the 2'-OH contributes to catalysis by donating a hydrogen bond to the 3'-oxygen leaving group in the transition state, presumably acting in conjunction with the metal ion to stabilize the developing negative charge. These same catalytic strategies of metal ion coordination and hydrogen bond donation operate in the Tetrahymena ribozyme reaction, suggesting that these features of catalysis have been conserved during evolution and thus extend to all group I introns. The two ribozymes also exhibit quantitative differences in their response to 3'-sulfur substitution. The Azoarcus ribozyme binds and cleaves the phosphorothiolate substrate more efficiently relative to the natural substrate than the Tetrahymena ribozyme under the same conditions, suggesting that the Azoarcus ribozyme better accommodates the phosphorothiolate at the active site both in the ground state and in the transition state. These differences may reflect either a less tightly knit Azoarcus structure and/or spatial deviations between backbone atoms in the two ribozymes that arise during divergent evolution, analogous to the well-documented relationship between protein sequence and structure.     

The P5 activator of a group IC ribozyme can replace the P7.1/7.2 activator of a group IA ribozyme

The P5 or P7 extensions in the group I intron ribozyme serve as "modular activator units" by stabilizing the conserved core of the ribozyme. The P5 extension of a group IC1 intron was introduced to a barely active group IA2 intron lacking its original P7 extension. The inserted P5 extension significantly activated the chimeric construct. Because the CYT-18 protein factor is also known to activate mutant group IA2 and IC1 introns lacking their P7 and P5 extensions, respectively, the RNA and protein activator units function in an analogous manner.     

Linking the group II intron catalytic domains: tertiary contacts and structural features of domain 3

Despite its importance for group II intron catalytic activity, structural information on conserved domain 3 (D3) is extremely limited. This domain is known to specifically stimulate the chemical rate of catalysis and to function as a 'catalytic effector'. Of all the long-range tertiary contacts that have been identified within group II introns, none has included D3 residues. Furthermore, little is known about the atoms and functional groups in D3 that contribute to catalysis. Using a nucleotide analog interference mapping assay with an extended repertoire of nucleotide analogs, we have identified functional groups in D3 that are critical for ribozyme activity. These data, together with mutational analysis, suggest the formation of noncanonical base pairs within the phylogenetically conserved internal loop at the base of D3. Finally, a related nucleotide analog interference suppression study resulted in the identification of a direct tertiary interaction between D3 and catalytic domain 5, which sheds new light on D3 function in the group II intron structure and mechanism.     

RNA splicing: group I intron crystal structures reveal the basis of splice site selection and metal ion catalysis

The group I intron has served as a model for RNA catalysis since its discovery 25 years ago. Four recently determined high-resolution crystal structures complement extensive biochemical studies on this system. Structures of the Azoarcus, Tetrahymena and bacteriophage Twort group I introns mimic different states of the splicing or ribozyme reaction pathway and provide information on splice site selection and metal ion catalysis. The 5'-splice site is selected by formation of a conserved G.U wobble pair between the 5'-exon terminus and the intron. The 3'-splice site is identified through stacking of three base triples, in which the middle triple contains the conserved terminal nucleotide of the intron, OmegaG. The structures support a two-metal-ion mechanism for group I intron splicing that might have corollaries to group II intron and pre-mRNA splicing by the spliceosome.     

Characterization of the Azoarcus ribozyme: tight binding to guanosine and substrate by an unusually small group I ribozyme

We report novel chemical properties of the ribozyme derived from the smallest group I intron (subgroup IC3) that comes from the pre-tRNA(Ile) of the bacterium Azoarcus sp. BH72. Despite the small size of the Azoarcus ribozyme (195 nucleotides (nt)), it binds tightly to the guanosine nucleophile (Kd = 15 +/- 3 microM) and exhibits activity at high temperatures (approximately 60-70 degrees C). These features may be due to the two GA3 tetraloop interactions postulated in the intron and the high GC content of the secondary structure. The second order rate constant for the Azoarcus ribozyme, ((k(cat)/Km)S = 8.4 +/- 2.1 x 10(-5) M(-1) min(-1)) is close to that found for the related ribozyme derived from the pre-tRNA(Ile) of the cyanobacterium Anabaena PCC7120. pH dependence studies and kinetic analyses of deoxy-substituted substrates suggest that the chemical cleavage step is the rate-determining process in the Azoarcus ribozyme. This may be due to the short 3-nt guide sequence-substrate pairing present in the Azoarcus ribozyme. Finally, the Azoarcus ribozyme shares features conserved in other group I ribozymes including the pH profile, the stereospecificity for the Rp-phosphorothioate at the cleavage site and the 1000-fold decrease in cleavage rate with a deoxyribonucleoside leaving group.     

Evolution of rbcL group IA introns and intron open reading frames within the colonial Volvocales (Chlorophyceae)

Mobile group I introns sometimes contain an open reading frame (ORF) possibly encoding a site-specific DNA endonuclease. However, previous phylogenetic studies have not clearly deduced the evolutionary roles of the group I intron ORFs. In this paper, we examined the phylogeny of group IA2 introns inserted in the position identical to that of the chloroplast-encoded rbcL coding region (rbcL-462 introns) and their ORFs from 13 strains of five genera (Volvox, Pleodorina, Volvulina, Astrephomene, and Gonium) of the colonial Volvocales (Chlorophyceae) and a related unicellular green alga, Vitreochlamys. The rbcL-462 introns contained an intact or degenerate ORF of various sizes except for the Gonium multicoccum rbcL-462 intron. Partial amino acid sequences of some rbcL-462 intron ORFs exhibited possible homology to the endo/excinuclease amino acid terminal domain. The distribution of the rbcL-462 introns is sporadic in the phylogenetic trees of the colonial Volvocales based on the five chloroplast exon sequences (6021 bp). Phylogenetic analyses of the conserved intron sequences resolved that the G. multicoccum rbcL-462 intron had a phylogenetic position separate from those of other colonial volvocalean rbcL-462 introns, indicating the recent horizontal transmission of the intron in the G. multicoccum lineage. However, the combined data set from conserved intron sequences and ORFs from most of the rbcL-462 introns resolved robust phylogenetic relationships of the introns that were consistent with those of the host organisms. Therefore, most of the extant rbcL-462 introns may have been vertically inherited from the common ancestor of their host organisms, whereas such introns may have been lost in other lineages during evolution of the colonial Volvocales. In addition, apparently higher synonymous substitutions than nonsynonymous substitutions in the rbcL-462 intron ORFs indicated that the ORFs might evolve under functional constraint, which could result in homing of the rbcL-462 intron in cases of spontaneous intron loss. On the other hand, the presence of intact to largely degenerate ORFs of the rbcL-462 introns within the three isolates of Gonium viridistellatum and the rare occurrence of the ORF-lacking rbcL-462 intron suggested that the ORFs might degenerate to result in the spontaneous intron loss during a very short evolutionary time following the loss of the ORF function. Thus, the sporadic distribution of the rbcL-462 introns within the colonial Volvocales can be largely explained by an equilibrium between maintenance of the introns by the intron ORF and spontaneous loss of introns when the introns do not have a functional ORF.     

Long-term evolution of the S788 fungal nuclear small subunit rRNA group I introns

More than 1000 group I introns have been identified in fungal rDNA. Little is known, however, of the splicing and secondary structure evolution of these ribozymes. Here, we use a combination of comparative and biochemical methods to address the evolution and splicing of a vertically inherited group I intron found at position 788 in the fungal small subunit (S) rRNA. The ancestral state of the S788 intron contains a highly conserved core and an extended P5 domain typical of IC1 introns. In contrast, the more derived introns have lost most of P5, and have an accelerated divergence rate within the core region with three functionally important substitutions that unambiguously separate them from the ancestral pool. Of 14 S788 group I introns that were tested for splicing, five, all of the ancestral type, were able to self-splice and produced intron RNA circles in vitro. The more derived S788 introns did not self-splice, and potentially rely on fungal-specific factors to facilitate splicing. In summary, we demonstrate one possible fate of vertically inherited group I introns, the loss of secondary structure elements, lessened selective constraints in the intron core, and ultimately, dependence on host-mediated splicing.     

The Neurospora crassa CYT-18 protein C-terminal RNA-binding domain helps stabilize interdomain tertiary interactions in group I introns

The Neurospora crassa mitochondrial tyrosyl-tRNA synthetase (CYT-18 protein) promotes the splicing of group I introns by stabilizing the catalytically active RNA structure. To accomplish this, CYT-18 recognizes conserved structural features of group I intron RNAs using regions of the N-terminal nucleotide-binding fold, intermediate alpha-helical, and C-terminal RNA-binding domains that also function in binding tRNA(Tyr). Curiously, whereas the splicing of the N. crassa mitochondrial large subunit rRNA intron is completely dependent on CYT-18's C-terminal RNA-binding domain, all other group I introns tested thus far are spliced efficiently by a truncated protein lacking this domain. To investigate the function of the C-terminal domain, we used an Escherichia coli genetic assay to isolate mutants of the Saccharomyces cerevisiae mitochondrial large subunit rRNA and phage T4 td introns that can be spliced in vivo by the wild-type CYT-18 protein, but not by the C-terminally truncated protein. Mutations that result in dependence on CYT-18's C-terminal domain include those disrupting two long-range GNRA tetraloop/receptor interactions: L2-P8, which helps position the P1 helix containing the 5'-splice site, and L9-P5, which helps establish the correct relative orientation of the P4-P6 and P3-P9 domains of the group I intron catalytic core. Our results indicate that different structural mutations in group I intron RNAs can result in dependence on different regions of CYT-18 for RNA splicing.     

A structural and phylogenetic analysis of the group IC1 introns in the order Bangiales (Rhodophyta)

Our previous study of the North American biogeography of Bangia revealed the presence of two introns inserted at positions 516 and 1506 in the nuclear-encoded SSU rRNA gene. We subsequently sequenced nuclear SSU rRNA in additional representatives of this genus and the sister genus Porphyra in order to examine the distribution, phylogeny, and structural characteristics of these group I introns. The lengths of these introns varied considerably, ranging from 467 to 997 nt for intron 516 and from 509 to 1,082 nt for intron 1506. The larger introns contained large insertions in the P2 domain of intron 516 and the P1 domain of intron 1506 that correspond to open reading frames (ORFs) with His-Cys box homing endonuclease motifs. These ORFs were found on the complementary strand of the 1506 intron in Porphyra fucicola and P. umbilicalis (HG), unlike the 516 intron in P. abbottae, P. kanakaensis, P. tenera (SK), Bangia fuscopurpurea (Helgoland), and B. fuscopurpurea (MA). Frameshifts were noted in the ORFs of the 516 introns in P. kanakaensis and B. fuscopurpurea (HL), and all ORFs terminated prematurely relative to the amino acid sequence for the homing endonuclease I-Ppo I. This raises the possibility that these sequences are pseudogenes. Phylogenies generated using sequences of both introns and the 18S rRNA gene were congruent, which indicated long-term immobility and vertical inheritance of the introns followed by subsequent loss in more derived lineages. The introns within the florideophyte species Hildenbrandia rubra (position 1506) were included to determine relationships with those in the Bangiales. The two sequences of intron 1506 analyzed in Hildenbrandia were positioned on a well-supported branch associated with members of the Bangiales, indicating possible common ancestry. Structural analysis of the intron sequences revealed a signature structural feature in the P5b domain of intron 516 that is unique to all Bangialean introns in this position and not seen in intron 1506 or other group IC1 introns.     

Analysis of the CYT-18 Protein Binding Site at the Junction of Stacked Helices in a Group I Intron RNA by Quantitative Binding Assays andin vitroSelection

TheNeurospora crassamitochondrial tyrosyl-tRNA synthetase (CYT-18 protein) functions in splicing group I introns by promoting the formation of the catalytically active structure of the intron RNA. Previous studies showed that CYT-18 binds with high affinity to the P4-P6 domain of the catalytic core and that there is some additional contribution to binding from the P3-P9 domain. Here, quantitative binding assays with deletion derivatives of theN. crassamitochondrial large rRNA intron showed that at least 70% of the binding energy can be accounted for by the interaction of CYT-18 with the P4-P6 domain. Within this domain, P4 and P6 are required for high affinity CYT-18 binding, while the distal elements P5 and P6a may contribute indirectly by stabilizing the correct structure of the binding site in P4 and P6. CYT-18 binds to a small RNA corresponding to the isolated P4-P6 domain, but not to a permuted version of this RNA in which P4-P6 is a continuous rather than a stacked helix. Iterativein vitroselection experiments with the isolated P4-P6 domain showed a requirement for base-pairing to maintain helices P4, P6 and P6a, but indicate that P5 is subject to fewer constraints. The most strongly conserved nucleotides in the selections were clustered around the junction of the P4-P6 stacked helix, with ten nucleotides (J3/4-2,3, P4 bp -1 and 3, and P6 bp -1 and 2) found invariant in the context of the wild-type RNA structure.In vitromutagenesis confirmed that replacement of the wild-type nucleotides at J3/4-2 and 3 or P4 bp-3 markedly decreased CYT-18 binding, reflecting either base specific contacts or indirect readout of RNA structure by the protein. Our results suggest that a major function of CYT-18 is to promote assembly of the P4-P6 domain by stabilizing the correct geometry at the junction of the P4-P6 stacked helix. The relatively large number of conserved nucleotides at the binding site suggests that the interaction of CYT-18 with group I introns is unlikely to have arisen by chance and could reflect either an evolutionary relationship between group I introns and tRNAs or interaction with a common stacked-helical structural motif that evolved separately in these RNAs.     

Small structural costs for evolution from RNA to RNP-based catalysis

Typical RNA-based cellular catalysts achieve their active structures only as complexes with protein cofactors, implying that protein binding compensates for some structural deficiencies in the RNA. An unresolved question was the extent to which protein-facilitation imposes additional structural costs, by requiring that an RNA maintain structures required for protein binding, beyond those required for catalysis. We used nucleotide analog interference to identify initially 71 functional group substitutions at phosphate, 2'-ribose, and adenosine base positions that compromise RNA self-splicing in the bI5 group I intron. Protein-facilitated splicing by CBP2 suppresses 11 of 30 interfering substitutions at the RNA backbone and a greater fraction, 27 of 41, at the adenosine base, including at structures conserved among group I introns. Only one substitution directly interferes with protein binding but not with self-splicing. This substitution, plus three adenosine base modifications that interfere more strongly in CBP2-dependent splicing than in self-splicing, yield a cost for protein facilitation of only four functional groups, as approximated by this set of analogs. The small observed structural cost provides a strong physical rationale for the evolutionary drive from RNA to RNP-based function in biology. Remarkably, the four extra requirements do not appear to report disruption of direct protein-RNA contacts and instead likely reflect design against misfolding rather than for maintenance of a protein-binding site.     

Naegleria nucleolar introns contain two group I ribozymes with different functions in RNA splicing and processing.

We have characterized the structural organization and catalytic properties of the large nucleolar group I introns (NaSSU1) of the different Naegleria species N. jamiesoni, N. andersoni, N. italica, and N. gruberi. NaSSU1 consists of three distinct RNA domains: an open reading frame encoding a homing-type endonuclease, and a small group I ribozyme (NaGIR1) inserted into the P6 loop of a second group I ribozyme (NaGIR2). The two ribozymes have different functions in RNA splicing and processing. NaGIR1 is an unusual self-cleaving group I ribozyme responsible for intron processing at two internal sites (IPS1 and IPS2), both close to the 5' end of the open reading frame. This processing is hypothesized to lead to formation of a messenger RNA for the endonuclease. Structurally, NaGIR2 is a typical group IC1 ribozyme, catalyzing intron excision and exon ligation reactions. NaGIR2 is responsible for circularization of the excised intron, a reaction that generates full-length RNA circles of wild-type intron. Although it is only distantly related in primary sequence, NaSSU1 RNA has a predicted organization and function very similar to that of the mobile group I intron DiSSU1 of Didymium, the only other group I intron known to encode two ribozymes. We propose that these twin-ribozyme introns define a distinct category of group I introns with a conserved structural organization and function.     

Exclusive conservation of mitochondrial group II intron nad4i548 among liverworts and its use for phylogenetic studies in this ancient plant clade

Liverworts occupy a pivotal position in land plant (embryophyte) phylogeny as the presumed earliest-branching major clade, sister to all other land plants, including the mosses, hornworts, lycophytes, monilophytes and seed plants. Molecular support for this earliest dichotomy in land plant phylogeny comes from strikingly different occurrences of introns in mitochondrial genes distinguishing liverworts from all other embryophytes. Exceptionally, however, the nad5 gene--the mitochondrial locus hitherto used most widely to elucidate early land plant phylogeny--carries a group I type intron that is shared between liverworts and mosses. We here explored whether a group II intron, the other major type of organellar intron, would similarly be conserved in position across the entire diversity of extant liverworts and could be of use for phylogenetic analyses in this supposedly most ancient embryophyte clade. To this end, we investigated the nad4 gene as a candidate locus possibly featuring different introns in liverworts as opposed to the non-liverwort embryophyte (NLE) lineage. We indeed found group II intron nad4i548 universally conserved in a wide phylogenetic sampling of 55 liverwort taxa, confirming clade specificity and surprising evolutionary stability of plant mitochondrial introns. As expected, intron nad4i548g2 carries phylogenetic information in its variable sequences, which confirms and extends previous cladistic insights on liverwort evolution. We integrate the new nad4 data with those of the previously established mitochondrial nad5 and the chloroplast rbcL and rps4 genes and present a phylogeny based on the fused datasets. Notably, the phylogenetic analyses suggest a reconsideration of previous phylogenetic and taxonomic assignments for the genera Calycularia and Mylia and resolve a sister group relationship of Ptilidiales and Porellales.