Atomic level architecture of group I introns revealed

Twenty-two years after their discovery as ribozymes, the self-splicing group I introns are finally disclosing their architecture at the atomic level. The crystal structures of three group I introns solved at moderately high resolution (3.1-3.8A) reveal a remarkably conserved catalytic core bound to the metal ions required for activity. The structure of the core is stabilized by an intron-specific set of long-range interactions that involves peripheral elements. Group I intron structures thus provide much awaited and extremely valuable snapshots of how these ribozymes coordinate substrate binding and catalysis.     

Functional identification of ligands for a catalytic metal ion in group I introns

Many enzymes use metal ions within their active sites to achieve enormous rate acceleration. Understanding how metal ions mediate catalysis requires elucidation of metal ion interactions with both the enzyme and the substrate(s). The three-dimensional arrangement determined by X-ray crystallography provides a powerful starting point for identifying ground state interactions, but only functional studies can establish and interrogate transition state interactions. The Tetrahymena group I ribozyme is a paradigm for the study of RNA catalysis, and previous work using atomic mutagenesis and quantitative analysis of metal ion rescue behavior identified catalytic metal ions making five contacts with the substrate atoms. Here, we have combined atomic mutagenesis with site-specific phosphorothioate substitutions in the ribozyme backbone to establish transition state ligands on the ribozyme for one of the catalytic metal ions, referred to as M A. We identified the pro-S P oxygen atoms at nucleotides C208, A304, and A306 as ground state ligands for M A, verifying interactions suggested by the Azoarcus crystal structures. We further established that these interactions are present in the chemical transition state, a conclusion that requires functional studies, such as those carried out herein. Elucidating these active site connections is a crucial step toward an in-depth understanding of how specific structural features of the group I intron lead to catalysis.     

Mutational analysis of the yeast U2 snRNA suggests a structural similarity to the catalytic core of group I introns.

We have used an in vitro reconstitution system to determine the effects of a large number of mutations in the highly conserved 5' terminal domain of the yeast U2 snRNA on pre-mRNA splicing. Whereas many mutations have little or no functional consequence, base substitutions in two regions were found to have drastic effects on pre-mRNA splicing. A previously unrecognized function for the U2 snRNA in the second step of splicing was found by alteration of the absolutely conserved sequence AGA upstream of the branch point recognition sequence. The effects of these mutations suggest the formation of a structure involving the U2 snRNA similar to the guanosine-binding site found in the catalytic core of group I introns.     

The antiquity of group I introns.

The recent discovery of self-splicing introns in cyanobacteria has given renewed interest to the question of whether introns may have been present in the ancestor of all living things. The properties of introns in genes of bacteria and bacteriophages are discussed in the context of their possible origin and biological function.     

Structure and assembly of group I introns

Self-splicing group I introns have served as a model for RNA catalysis and folding for over two decades. New three-dimensional structures now bring the details into view. Revelations include an unanticipated turn in the RNA backbone around the guanosine-binding pocket. Two metal ions in the active site coordinate the substrate and phosphates from all three helical domains.     

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.     

Variations in sequence and occurrence of SSU rDNA group I introns in Monilinia fructicola isolates

A group I intron of 418 base pairs in the Monilinia fructicola ribosomal small-subunit sequence was characterized. The absence of such an intron in M. laxa and M. fructigena led to a PCR test for M. fructicola identification based on the presence of this intron. The failure to amplify a PCR fragment for some isolates of M. fructicola recently lead to speculation that the intron might not be present always in M. fructicola. In this study, we analyzed 13 isolates of M. fructicola and found that the intron was absent in four isolates and we determined from sequence analysis that there are several nucleotide variations that allow the M. fructicola ribosomal SSU intron to be grouped into 6 polymorphic types.     

Identification of a family of group II introns encoding LAGLIDADG ORFs typical of group I introns

Group I and group II introns are unrelated classes of introns that each encode proteins that facilitate intron splicing and intron mobility. Here we describe a new subfamily of nine introns in fungi that are group II introns but encode LAGLIDADG ORFs typical of group I introns. The introns have fairly standard group IIB1 RNA structures and are inserted into three different sites in SSU and LSU rRNA genes. Therefore, introns should not be assumed to be group I introns based solely on the presence of a LAGLIDADG ORF.     

Distribution, sequence homology, and homing of group I introns among T-even-like bacteriophages: evidence for recent transfer of old introns

Self-splicing group I introns are being found in an increasing number of bacteriophages. Most introns contain an open reading frame coding for a homing endo-nuclease that confers mobility to both the intron and the homing endonuclease gene (HEG). The frequent occurrence of intron/HEG has raised questions whether group I introns are spread via horizontal transfer between phage populations. We have determined complete sequences for the known group I introns among T-even-like bacteriophages together with sequences of the intron-containing genes td, nrdB, and nrdD from phages with and without introns. A previously uncharacterized phage isolate, U5, is shown to contain all three introns, the only phage besides T4 found with a "full set" of these introns. Sequence analysis of td and nrdB genes from intron-containing and intronless phages provides evidence that recent horizontal transmission of introns has occurred among the phages. The fact that several of the HEGs have suffered deletions rendering them non-functional implies that the homing endonucleases are of no selective advantage to the phage and are rapidly degenerating and probably dependent upon frequent horizontal transmissions for maintenance within the phage populations. Several of the introns can home to closely related intronless phages during mixed infections. However, the efficiency of homing varies and is dependent on homology in regions flanking the intron insertion site. The occurrence of optional genes flanking the respective intron-containing gene can strongly affect the efficiency of homing. These findings give further insight into the mechanisms of propagation and evolution of group I introns among the T-even-like bacteriophages.     

Mutagenesis and comparative sequence analysis of a base triple joining the two domains of group I ribozymes.

Tertiary interactions are important in the higher-order folding of catalytic RNAs. Recently, a base triple, joining the two major domains of the catalytic core, was determined in group I introns from the cyanobacterium Anabaena PCC7120 and the eukaryote Tetrahymena thermophila. This base triple involves the fifth base pair of P4 and the fifth base of the single-stranded region J8/7. We made base pair and single-nucleotide substitutions in the fifth base pair of P4, a G-C in the wild-type Anabaena intron, and tested them for self-splicing activity. The results suggest a hydrogen bonding model in which only the C of the base pair interacts directly with the fifth base of J8/7. Comparative sequence analysis was used to determine the different combinations of base triples that occur in approximately 450 natural group I introns identified to date. About 94% of the base triples analyzed are compatible with the proposed hydrogen bonding model. Disrupting this base triple in the Tetrahymena intron resulted in the disappearance of splicing intermediates (intron 3' exon and 5' exon), even though the first step of splicing was not affected. Restoration of the base triple by a compensatory mutation reverted the intermediates to wild-type levels. These results suggest that disruption of the base triple increases the rate of the second step of splicing or of a conformational change preceding the second step. Repositioning of the base triple to form a new set of interactions may be required for the second step of splicing.     

Modeling a minimal ribosome based on comparative sequence analysis

We have determined the three-dimensional organization of ribosomal RNAs and proteins essential for minimal ribosome function. Comparative sequence analysis identifies regions of the ribosome that have been evolutionarily conserved, and the spatial organization of conserved domains is determined by mapping these onto structures of the 30S and 50S subunits determined by X-ray crystallography. Several functional domains of the ribosome are conserved in their three-dimensional organization in the Archaea, Bacteria, Eucaryotic nuclear, mitochondria and chloroplast ribosomes. In contrast, other regions from both subunits have shifted their position in three-dimensional space during evolution, including the L11 binding domain and the alpha-sarcin-ricin loop (SRL). We examined conserved bridge interactions between the two ribosomal subunits, giving an indication of which contacts are more significant. The tRNA contacts that are conserved were also determined, highlighting functional interactions as the tRNA moves through the ribosome during protein synthesis. To augment these studies of a large collection of comparative structural models sampled from all major branches on the phylogenetic tree, Caenorhabditis elegans mitochondrial rRNA is considered individually because it is among the smallest rRNA sequences known. The C.elegans model supports the large collection of comparative structure models while providing insight into the evolution of mitochondrial ribosomes.     

Strain-dependent sequence heterogeneity in the nuclear group I introns and large subunit ribosomal RNA in Physarum polycephalum.

The DNA sequence of an internal EcoRI fragment in the large subunit ribosomal RNA gene from the myxomycete Physarum polycephalum, strain M3C, has been determined. This ribosomal DNA fragment contain two group I introns. Sequence heterogeneity in the two introns and in the coding region of the ribosomal RNA were identified upon comparison to the strain PPO-1 sequence. The nature of the sequence variations is discussed.     

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 three-dimensional architecture of the class I ligase ribozyme

The class I ligase ribozyme catalyzes a Mg(++)-dependent RNA-ligation reaction that is chemically analogous to a single step of RNA polymerization. Indeed, this ribozyme constitutes the catalytic domain of an accurate and general RNA polymerase ribozyme. The ligation reaction is also very rapid in both single- and multiple-turnover contexts and thus is informative for the study of RNA catalysis as well as RNA self-replication. Here we report the initial characterization of the three-dimensional architecture of the ligase. When the ligase folds, several segments become protected from hydroxyl-radical cleavage, indicating that the RNA adopts a compact tertiary structure. Ribozyme folding was largely, though not completely, Mg(++) dependent, with a K(1/2[Mg]) < 1 mM, and was observed over a broad temperature range (20 degrees C -50 degrees C). The hydroxyl-radical mapping, together with comparative sequence analyses and analogy to a region within 23S ribosomal RNA, were used to generate a three-dimensional model of the ribozyme. The predictive value of the model was tested and supported by a photo-cross-linking experiment.     

Vertical evolution and intragenic spread of lichen-fungal group I introns

One family within the Euascomycetes (Ascomycota), the lichen-forming Physciaceae, is particularly rich in nuclear ribosomal [r]DNA group I introns. We used phylogenetic analyses of group I introns and lichen-fungal host cells to address four questions about group I intron evolution in lichens, and generally in all eukaryotes: 1) Is intron spread in the lichens associated with the intimate association of the fungal and photosynthetic cells that make up the lichen thallus? 2) Are the multiple group I introns in the lichen-fungi of independent origins, or have existing introns spread into novel sites in the rDNA? 3) If introns have moved to novel sites, then does the exon context of these sites provide insights into the mechanism of intron spread? and 4) What is the pattern of intron loss in the small subunit rDNA gene of lichen-fungi? Our analyses show that group I introns in the lichen-fungi and in the lichen-algae (and lichenized cyanobacteria) do not share a close evolutionary relationship, suggesting that these introns do not move between the symbionts. Many group I introns appear to have originated in the common ancestor of the Lecanorales, whereas others have spread within this lineage (particularly in the Physciaceae) putatively through reverse-splicing into novel rRNA sites. We suggest that the evolutionary history of most lichen-fungal group I introns is characterized by rare gains followed by extensive losses in descendants, resulting in a sporadic intron distribution. Detailed phylogenetic analyses of the introns and host cells are required, therefore, to distinguish this scenario from the alternative hypothesis of widespread and independent intron gains in the different lichen-fungal lineages.     

Molecular phylogeny of phyllopharyngean ciliates and their group I introns

We analyzed small subunit ribosomal DNA (ssu-rDNA) sequences to evaluate both the monophyly of the ciliate class Phyllopharyngea de Puytorac et al. (1974), and relationships among subclasses. Classifications based on morphology and ultrastructure divide the Phyllopharyngea into four subclasses, the Phyllopharyngia, Chonotrichia, Rhynchodia, and Suctoria. Our analyses of ssu-rDNA genealogies derived from sequence data collected from diverse members representing three of the four subclasses of Phyllopharyngea (Suctoria: Ephelota spp., Prodiscophyra collini, Acineta sp.; Phyllopharyngia: Chlamydodon exocellatus, Chlamydodon triquetrus, Dysteria sp.; and Chonotrichia: Isochona sp.) provide strong support for the monophyly of the Phyllopharyngea, and show that the Chonotrichia emerge from within the Phyllopharyngia. Based on this initial sampling, suctorian budding types are monophyletic, and exogenous budding appears to be basal to evaginative and endogenous budding. Further, we report the discovery of a group I intron at position 891 in the Suctoria Acineta sp. and Tokophrya lemnarum, and a second group I intron at position 1506 in T. lemnarum. These introns represent only the second examples of group I introns in a ciliate ribosomal gene, since the discovery of ribozymes in the LSU rRNA gene of Tetrahymena thermophila. Phylogenetic analyses of Group I introns suggest a complex evolutionary history involving either multiple loses or gains of introns within endogenously budding Suctoria.     

Toward predicting self-splicing and protein-facilitated splicing of group I introns

In the current era of massive discoveries of noncoding RNAs within genomes, being able to infer a function from a nucleotide sequence is of paramount interest. Although studies of individual group I introns have identified self-splicing and nonself-splicing examples, there is no overall understanding of the prevalence of self-splicing or the factors that determine it among the >2300 group I introns sequenced to date. Here, the self-splicing activities of 12 group I introns from various organisms were assayed under six reaction conditions that had been shown previously to promote RNA catalysis for different RNAs. Besides revealing that assessing self-splicing under only one condition can be misleading, this survey emphasizes that in vitro self-splicing efficiency is correlated with the GC content of the intron (>35% GC was generally conductive to self-splicing), and with the ability of the introns to form particular tertiary interactions. Addition of the Neurospora crassa CYT-18 protein activated splicing of two nonself-splicing introns, but inhibited the second step of self-splicing for two others. Together, correlations between sequence, predicted structure and splicing begin to establish rules that should facilitate our ability to predict the self-splicing activity of any group I intron from its sequence.     

Selection of the 3′-splice site in group I introns

A model for selection of 3′-splice sites in splicing of RNA precursors containing group I introns is presented. The key feature of this model is a newly identified tertiary interaction between the catalytic core of the intron and the 3′-splice site. This tertiary pairing would bring the 3′-splice site into the core of the intron, which is known to contain RNA sequences and structures essential for catalyzing the splicing reactions. The proposed tertiary interaction can coexist with P10, a pairing between 3′-exon sequences and the ‘internal guide sequence’ near the 5′-end of the intron. The model predicts that three RNA-RNA interactions are important in selection of 3′-splice sites: (i) binding of intron sequences with the core; (ii) pairing of exon sequences with the internal guide sequence; and (iii) binding of the terminal guanosine to an unknown site within the core.