Nested Evolution of a tRNALeu(UAA) Group I Intron by both Horizontal Intron Transfer and Recombination of the Entire tRNA Locus

The origin and evolution of bacterial introns are still controversial issues. Here we present data on the distribution and evolution of a recently discovered divergent tRNA(Leu)(UAA) intron. The intron shows a higher sequence affiliation with introns in tRNA(Ile)(CAU) and tRNA(Arg)(CCU) genes in alpha- and beta-proteobacteria, respectively, than with other cyanobacterial tRNA(Leu)(UAA) group I introns. The divergent tRNA(Leu)(UAA) intron is sporadically distributed both within the Nostoc and the Microcystis radiations. The complete tRNA gene, including flanking regions and intron from Microcystis aeruginosa strain NIVA-CYA 57, was sequenced in order to elucidate the evolutionary pattern of this intron. Phylogenetic reconstruction gave statistical evidence for different phylogenies for the intron and exon sequences, supporting an evolutionary model involving horizontal intron transfer. The distribution of the tRNA gene, its flanking regions, and the introns were addressed by Southern hybridization and PCR amplification. The tRNA gene, including the flanking regions, were absent in the intronless stains but present in the intron-containing strains. This suggests that the sporadic distribution of this intron within the Microcystis genus cannot be attributed to intron mobility but rather to an instability of the entire tRNA(Leu)(UAA) intron-containing genome region. Taken together, the complete data set for the evolution of this intron can best be explained by a model involving a nested evolution of the intron, i.e., wherein the intron has been transferred horizontally (probably through a single or a few events) to a tRNA(Leu)(UAA) gene which is located within a unstable genome region.     

Origin and evolution of group I introns in cyanobacterial tRNA genes.

Many tRNA(Leu)UAA genes from plastids contain a group I intron. An intron is also inserted in the same gene at the same position in cyanobacteria, the bacterial progenitors of plastids, suggesting an ancient bacterial origin for this intron. A group I intron has also been found in the tRNA(fMet) gene of some cyanobacteria but not in plastids, suggesting a more recent origin for this intron. In this study, we investigate the phylogenetic distributions of the two introns among cyanobacteria, from the earliest branching to the more derived species. The phylogenetic distribution of the tRNA(Leu)UAA intron follows the clustering of rRNA sequences, being either absent or present in clades of closely related species, with only one exception in the Pseudanabaena group. Our data support the notion that the tRNA(Leu)UAA intron was inherited by cyanobacteria and plastids through a common ancestor. Conversely, the tRNA(fMet) intron has a sporadic distribution, implying that many gains and losses occurred during cyanobacterial evolution. Interestingly, a phylogenetic tree inferred from intronic sequences clearly separates the different tRNA introns, suggesting that each family has its own evolutionary history.     

Self-splicing group I intron in cyanobacterial initiator methionine tRNA: evidence for lateral transfer of introns in bacteria.

A group I self-splicing intron has been found in the anticodon loop of tRNA(fMet) genes in three cyanobacterial genera: Dermocarpa, Scytonema and Synechocystis; it is absent in nine others. The Synechocystis intron is also interrupted by an open reading frame (ORF) of 150 codons. Of these three bacteria, only Scytonema also contains the group I intron that has previously been reported in tRNA(Leu) (UAA) genes in both cyanobacteria and chloroplasts. The presence of an ORF in the tRNA(fMet) intron, the sporadic distribution of the intron among cyanobacteria and the lack of correlation between relatedness of the intron sequences and the bacteria in which they reside, are all consistent with recent introduction of this intron by lateral transfer.     

The natural history of group I introns

There are four major classes of introns: self-splicing group I and group II introns, tRNA and/or archaeal introns and spliceosomal introns in nuclear pre-mRNA. Group I introns are widely distributed in protists, bacteria and bacteriophages. Group II introns are found in fungal and land plant mitochondria, algal plastids, bacteria and Archaea. Group II and spliceosomal introns share a common splicing pathway and might be related to each other. The tRNA and/or archaeal introns are found in the nuclear tRNA of eukaryotes and in archaeal tRNA, rRNA and mRNA. The mechanisms underlying the self-splicing and mobility of a few model group I introns are well understood. By contrast, the role of these highly distinct processes in the evolution of the 1500 group I introns found thus far in nature (e.g. in algae and fungi) has only recently been clarified. The explosion of new sequence data has facilitated the use of comparative methods to understand group I intron evolution in a broader context and to generate hypotheses about intron insertion, splicing and spread that can be tested experimentally.     

The anticodon of the maize chloroplast gene for tRNA Leu UAA is split by a large intron.

The maize chloroplast gene encoding tRNA Leu UAA has been sequenced. It contains a 458 base pair intron between the first and second bases of the anticodon. The tRNA is 88 nucleotides long (the 3'-terminal CCA sequence included which, however, is not encoded by the gene) and differs in only four nucleotides (modified nucleotides are not considered) from the corresponding isoacceptor from bean chloroplasts. The unusual position of the intron in this maize chloroplast tRNA gene suggests a splicing model different from that generally accepted for eukaryotic split tRNA genes.     

A group I plant intron accumulates as circular RNA forms with extensive 5' deletions in vivo.

Analysis by a PAGE approach for detecting small circular RNAs showed the existence of one such molecular species (RNA 1) accumulating at high levels in cherimoya. Sequencing of cDNA clones of RNA 1 revealed a size of 281 nt and a sequence identical to the 3'-terminal region of the 494-nt tRNALeu(UAA) group I intron from cherimoya. Northern blot hybridizations with a probe complementary to RNA 1 showed that this RNA coexists in vivo with its corresponding linear form, with the presumed full-length intron, and with minor amounts of two additional small circular species (RNAs 2 and 3). RNAs 2 and 3 had sizes of 216 and 156 nt, respectively, and sequences identical to different moieties of the 3'-terminal region of the tRNALeu(UAA) intron. The three cyclization sites giving rise to RNAs 1, 2, and 3, located within loop 8, are preceded by CUU or UUU trinucleotides and followed by sequences capable of forming base pairing interactions with the internal guide sequence characteristic of group I introns. The good correlation observed between the stabilities of these interactions and the in vivo accumulation levels of the corresponding cherimoya circular RNAs support the hypothesis that they emerge through a common mechanism similar to that advanced previously for the generation of circular RNAs derived from other group I introns. The lack of interactions of similar stabilities in tobacco, in which no circular RNAs derived from the tRNALeu(UAA) intron were detected, is consistent with this proposal, although other factors are also probably important in the synthesis and accumulation of the small circular RNAs in cherimoya.     

The restriction fold turns to the dark side: a bacterial homing endonuclease with a PD-(D/E)-XK motif

The homing endonuclease I-Ssp6803I causes the insertion of a group I intron into a bacterial tRNA gene-the only example of an invasive mobile intron within a bacterial genome. Using a computational fold prediction, mutagenic screen and crystal structure determination, we demonstrate that this protein is a tetrameric PD-(D/E)-XK endonuclease - a fold normally used to protect a bacterial genome from invading DNA through the action of restriction endonucleases. I-Ssp6803I uses its tetrameric assembly to promote recognition of a single long target site, whereas restriction endonuclease tetramers facilitate cooperative binding and cleavage of two short sites. The limited use of the PD-(D/E)-XK nucleases by mobile introns stands in contrast to their frequent use of LAGLIDADG and HNH endonucleases - which in turn, are rarely incorporated into restriction/modification systems.     

A chemical phylogeny of group I introns based upon interference mapping of a bacterial ribozyme

Despite its small size, the 205 nt group I intron from Azoarcus tRNA(Ile) is an exceptionally stable self-splicing RNA. This IC3 class intron retains the conserved secondary structural elements common to group I ribozymes, but lacks several peripheral helices. These features make it an ideal system to establish the conserved chemical basis of group I intron activity. We collected nucleotide analog interference mapping (NAIM) data of the Azoarcus intron using 14 analogs that modified the phosphate backbone, the ribose sugar, or the purine base functional groups. In conjunction with a complete interference set collected on the Tetrahymena group I intron (IC1 class), these data define a "chemical phylogeny" of functional groups that are important for the activity of both introns and that may be common chemical features of group I intron catalysts. The data identify the functional moieties most likely to play a conserved role as ligands for catalytic metal ions, the substrate helix, and the guanosine cofactor. These include backbone functional groups whose nucleotide identity is not conserved, and hence are difficult to identify by standard phylogenetic sequence comparisons. The data suggest that both introns utilize an equivalent set of long range tertiary interactions for 5'-splice site selection between the P1 substrate helix and its receptor in the J4/5 asymmetric bulge, as well as an equivalent set of 2'-OH groups for P1 helix docking into most of the single stranded segment J8/7. However, the Azoarcus intron appears to make an alternative set of interactions at the base of the P1 helix and at the 5'-end of the J8/7. Extensive differences were observed within the intron peripheral domains, particularly in P2 and P8 where the Azoarcus data strongly support the proposed formation of a tetraloop-tetraloop receptor interaction. This chemical phylogeny for group I intron catalysis helps to refine structural models of the RNA active site and identifies functional groups that should be carefully investigated for their role in transition state stabilization.     

The cyanobacterial tRNA(Leu) (UAA) intron: evolutionary patterns in a genetic marker

Cyanobacterial tRNA(Leu) (UAA) intron sequences from natural populations of Nostoc and other cyanobacteria were compared. Variation between the different introns was not randomly distributed but strongly restricted by the secondary and tertiary structure of the intron. Although all Nostoc sequences examined shared high similarity, differences were observed in one stem-loop. This stem-loop could be divided into two classes, both built up from two base pairing heptanucleotide repeats. Size variation was primarily caused by different numbers of repeats, but some strains also contained additional sequences in this stem-loop not following the heptanucleotide repeat motif. Several sequences showing similarity with these additional sequences were identified in the Nostoc punctiforme genome. Furthermore, the regions flanking these sequences contained the same, or similar, heptanucleotide repeats as those flanking the corresponding sequences in the intron. It is proposed that both slipped strand mispairing during replication and homologous recombination among different loci in the genome are important processes causing variation between introns.     

NMR Structure of the C-terminal domain of a tyrosyl-tRNA synthetase that functions in group I intron splicing

The mitochondrial tyrosyl-tRNA synthetases (mt TyrRSs) of Pezizomycotina fungi are bifunctional proteins that aminoacylate mitochondrial tRNA(Tyr) and are structure-stabilizing splicing cofactors for group I introns. Studies with the Neurospora crassa synthetase (CYT-18 protein) showed that splicing activity is dependent upon Pezizomycotina-specific structural adaptations that form a distinct group I intron-binding site in the N-terminal catalytic domain. Although CYT-18's C-terminal domain also binds group I introns, it has been intractable to X-ray crystallography in the full-length protein. Here, we determined an NMR structure of the isolated C-terminal domain of the Aspergillus nidulans mt TyrRS, which is closely related to but smaller than CYT-18's. The structure shows an S4 fold like that of bacterial TyrRSs, but with novel features, including three Pezizomycontia-specific insertions. (15)N-(1)H two-dimensional NMR showed that C-terminal domains of the full-length A. nidulans and Geobacillus stearothermophilus synthetases do not tumble independently in solution, suggesting restricted orientations. Modeling onto a CYT-18/group I intron cocrystal structure indicates that the C-terminal domains of both subunits of the homodimeric protein bind different ends of the intron RNA, with one C-terminal domain having to undergo a large shift on its flexible linker to bind tRNA(Tyr) or the intron RNA on either side of the catalytic domain. The modeling suggests that the C-terminal domain acts together with the N-terminal domain to clamp parts of the intron's catalytic core, that at least one C-terminal domain insertion functions in group I intron binding, and that some C-terminal domain regions bind both tRNA(Tyr) and group I intron RNAs.     

Nucleolar introns from Physarum flavicomum contain insertion elements that may explain how mobile group I introns gained their open reading frames.

Comparison of two group I intron sequences in the nucleolar genome of the myxomycete Physarum flavicomum to their homologs in the closely related Physarum polycephalum revealed insertion-like elements. One of the insertion-like elements consists of two repetitive sequence motifs of 11 and 101 bp in five and three copies, respectively. The smaller motif, which flanks the larger, resembles a target duplication and indicates a relationship to transposons or retroelements. The insertion-like elements are found in the peripheral loops of the RNA structure; the positions occupied by the ORFs of mobile nucleolar group I introns. The P. flavicomum introns are 1184 and 637 bp in size, located in the large subunit ribosomal RNA gene, and can be folded into group I intron structures at the RNA level. However, the intron 2s from both P. flavicomum and P. polycephalum contain an unusual core region that lacks the P8 segment. None of the introns are able to self-splice in vitro. Southern analysis of different isolates indicates that the introns are not optional in myxomycetes.     

Novel system for analysis of group I 3' splice site reactions based on functional trans-interaction of the P1/P10 reaction helix with the ribozyme's catalytic core.

A group I intron from a bacterial tRNA precursor has been converted into an RNA enzyme that catalyzes the efficient polymerization of oligoribonucleotide analogs of tRNA exons using a reaction scheme consisting of multiple cycles of reverse and forward exon ligation reactions. Here, we present results showing that this system represents a novel and useful tool for the analysis of 3' splice site reactions of group I ribozymes. First, analysis of variant substrates containing base substitutions in group I secondary structure elements P1, P9.0 and P10 confirms that exon polymerization is dependent on these structures, and therefore constitutes an appropriate and relevant model system for studying the exon ligation step of splicing. Second, to probe interactions between the intron's catalytic core and the bases and backbone of the P1/P10 reaction helix, two successful strategies for separating the internal guide sequence from the intron core were devised. One such strategy uses a construct in which the reaction helix interacts functionally with the catalytic core using only tertiary contacts. Further stabilization of this interaction through the inclusion of a 7 bp intermolecular P2 helix generates increased reaction efficiency. Third, when provided with two reaction helices, the ribozyme synthesizes mixed polymers through a mechanism that involves sequential binding and release of the duplexes. Fourth, in these reactions, turnover of the external guide sequence requires unwinding and annealing of the P2 helix, suggesting that P2 unwinding may occur during group I splicing. These results provide novel experimental tools to probe the relatively poorly understood 3' splice site reactions of group I introns, and may be relevant to ribozyme-catalyzed assembly and recombination of oligomers in prebiotic scenarios.     

The IStron CdISt1 of Clostridium difficile: molecular symbiosis of a group I intron and an insertion element

The IStron CdISt1 was first discovered as an insertion into the tcdA gene of the clinical isolate C34. It combines structural and functional properties of a group I intron at its 5'-end with those of an insertion element at its 3'-end. Up to date four different types could be found, mainly differing in their IS-element portions. Contrasting classical group I introns, CdISt1 is always integrated in ORFs encoding bacterial protein. In case CdISt1 had only the IS-element function such insertion would inactivate the protein encoded by the host gene. It is only due to the self-splicing activity of the group I intron parts that CdISt1 integration does not abolish protein function. Both elements seem to exist in molecular symbiosis and CdISt1 could thus be a prototype of a novel class of genetic elements. Moreover, integration of the CdISt1 into the genome could be advantageous for the bacterium, a motor function for evolution of bacterial proteins is discussed. In clinical practice CdISt1 might well serve as a tool for epidemiological studies of C. difficile infections.     

The recent transfer of a homing endonuclease gene

The myxomycete Didymium iridis (isolate Panama 2) contains a mobile group I intron named Dir.S956-1 after position 956 in the nuclear small subunit (SSU) rRNA gene. The intron is efficiently spread through homing by the intron-encoded homing endonuclease I-DirI. Homing endonuclease genes (HEGs) usually spread with their associated introns as a unit, but infrequently also spread independent of introns (or inteins). Clear examples of HEG mobility are however sparse. Here, we provide evidence for the transfer of a HEG into a group I intron named Dir.S956-2 that is inserted into the SSU rDNA of the Costa Rica 8 isolate of D.iridis. Similarities between intron sequences that flank the HEG and rDNA sequences that flank the intron (the homing endonuclease recognition sequence) suggest that the HEG invaded the intron during the recent evolution in a homing-like event. Dir.S956-2 is inserted into the same SSU site as Dir.S956-1. Remarkably, the two group I introns encode distantly related splicing ribozymes with phylogenetically related HEGs inserted on the opposite strands of different peripheral loop regions. The HEGs are both interrupted by small spliceosomal introns that must be removed during RNA maturation.     

Complex Evolutionary Patterns of tRNAUAALeu Group I Introns in Cyanobacterial Radiation

Based on the findings that plastids and cyanobacteria have similar group I introns inserted into tRNAUAALeu genes, these introns have been suggested to be immobile and of ancient origin. In contrast, recent evidence suggests lateral transfer of cyanobacterial group I introns located in tRNAUAALeu genes. In light of these new findings, we have readdressed the evolution and lateral transfer of tRNAUAALeu group I introns in cyanobacteral radiation. We determined the presence of introns in 38 different strains, representing the major cyanobacterial lineages, and characterized the introns in 22 of the strains. Notably, two of these strains have two tRNAUAALeu genes, with each of these genes interrupted by introns, while three of the strains have both interrupted and uninterrupted genes. Two evolutionary distinct clusters of tRNA genes, with the genes interrupted by introns belonging to two distinct intron clusters, were identified. We also compared 16S rDNA and intron evolution for both closely and distantly related strains. The distribution of the introns in the clustered groups, as defined from 16S rDNA analysis, indicates relatively recent gain and/or loss of the introns in some of these lineages. The comparative analysis also suggests differences in the phylogenetic trees for 16S rDNA and the tRNAUAALeu group I introns. Taken together, our results show that the evolution of the intron is considerably more complex than previous studies found to be the case. We discuss, based on our results, evolutionary models involving lateral intron transfer and models involving differential loss of the intron.     

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.     

Intron phylogeny: a new hypothesis

The three major classes of intron are clearly of unequal antiquity. Structured (often self-splicing and sometimes mobile) introns are the most ancient, probably dating (at least for group I) from the ancestral (eubacterial) cell 3500 million years ago, and were originally restricted to tRNA. Protein-spliced introns (usually in tRNA) probably evolved from them by a radical change in splicing mechanism in the common ancestor of eukaryotes and archaebacteria, perhaps only about 1700 million years ago. Spliceosomal introns probably evolved from group-II-like self-splicing introns after the origin of the nucleus between 1700 and 1000 million years ago, and were probably mostly inserted into previously unsplit protein-coding genes after the origin of mitochondria 1000 million years ago.     

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.