Nuclear-encoded rDNA group I introns: origin and phylogenetic relationships of insertion site lineages in the green algae

Group I introns are widespread in eukaryotic organelles and nuclear- encoded ribosomal DNAs (rDNAs). The green algae are particularly rich in rDNA group I introns. To better understand the origins and phylogenetic relationships of green algal nuclear-encoded small subunit rDNA group I introns, a secondary structure-based alignment was constructed with available intron sequences and 11 new subgroup ICI and three new subgroup IB3 intron sequences determined from members of the Trebouxiophyceae (common phycobiont components of lichen) and the Ulvophyceae. Phylogenetic analyses using a weighted maximum-parsimony method showed that most group I introns form distinct lineages defined by insertion sites within the SSU rDNA. The comparison of topologies defining the phylogenetic relationships of 12 members of the 1512 group I intron insertion site lineage (position relative to the E. coli SSU rDNA coding region) with that of the host cells (i.e., SSU rDNAs) that contain these introns provided insights into the possible origin, stability, loss, and lateral transfer of ICI group I introns. The phylogenetic data were consistent with a viral origin of the 1512 group I intron in the green algae. This intron appears to have originated, minimally, within the SSU rDNA of the common ancestor of the trebouxiophytes and has subsequently been vertically inherited within this algal lineage with loss of the intron in some taxa. The phylogenetic analyses also suggested that the 1512 intron was laterally transferred among later-diverging trebouxiophytes; these algal taxa may have coexisted in a developing lichen thallus, thus facilitating cell- to-cell contact and the lateral transfer. Comparison of available group I intron sequences from the nuclear-encoded SSU rDNA of phycobiont and mycobiont components of lichens demonstrated that these sequences have independent origins and are not the result of lateral transfer from one component to the other.      

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.     

The evolution of homing endonuclease genes and group I introns in nuclear rDNA

Group I introns are autonomous genetic elements that can catalyze their own excision from pre-RNA. Understanding how group I introns move in nuclear ribosomal (r)DNA remains an important question in evolutionary biology. Two models are invoked to explain group I intron movement. The first is termed homing and results from the action of an intron-encoded homing endonuclease that recognizes and cleaves an intronless allele at or near the intron insertion site. Alternatively, introns can be inserted into RNA through reverse splicing. Here, we present the sequences of two large group I introns from fungal nuclear rDNA, which both encode putative full-length homing endonuclease genes (HEGs). Five remnant HEGs in different fungal species are also reported. This brings the total number of known nuclear HEGs from 15 to 22. We determined the phylogeny of all known nuclear HEGs and their associated introns. We found evidence for intron-independent HEG invasion into both homologous and heterologous introns in often distantly related lineages, as well as the "switching" of HEGs between different intron peripheral loops and between sense and antisense strands of intron DNA. These results suggest that nuclear HEGs are frequently mobilized. HEG invasion appears, however, to be limited to existing introns in the same or neighboring sites. To study the intron-HEG relationship in more detail, the S943 group I intron in fungal small-subunit rDNA was used as a model system. The S943 HEG is shown to be widely distributed as functional, inactivated, or remnant ORFs in S943 introns.     

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.     

Separate origins of group I introns in two mitochondrial genes of the katablepharid Leucocryptos marina

Mitochondria are descendants of the endosymbiotic α-proteobacterium most likely engulfed by the ancestral eukaryotic cells, and the proto-mitochondrial genome should have been severely streamlined in terms of both genome size and gene repertoire. In addition, mitochondrial (mt) sequence data indicated that frequent intron gain/loss events contributed to shaping the modern mt genome organizations, resulting in the homologous introns being shared between two distantly related mt genomes. Unfortunately, the bulk of mt sequence data currently available are of phylogenetically restricted lineages, i.e., metazoans, fungi, and land plants, and are insufficient to elucidate the entire picture of intron evolution in mt genomes. In this work, we sequenced a 12 kbp-fragment of the mt genome of the katablepharid Leucocryptos marina. Among nine protein-coding genes included in the mt genome fragment, the genes encoding cytochrome b and cytochrome c oxidase subunit I (cob and cox1) were interrupted by group I introns. We further identified that the cob and cox1 introns host open reading frames for homing endonucleases (HEs) belonging to distantly related superfamilies. Phylogenetic analyses recovered an affinity between the HE in the Leucocryptos cob intron and two green algal HEs, and that between the HE in the Leucocryptos cox1 intron and a fungal HE, suggesting that the Leucocryptos cob and cox1 introns possess distinct evolutionary origins. Although the current intron (and intronic HE) data are insufficient to infer how the homologous introns were distributed to distantly related mt genomes, the results presented here successfully expanded the evolutionary dynamism of group I introns in mt genomes.     

Origin and Inheritance of Group I Introns in 26S rRNA Genes of Gaeumannomyces graminis

Studies of the distribution of the three group I introns (intron A, intron T, and intron AT) in the 26S rDNA of Gaeumannomyces graminis had suggested that they were transferred to a common ancestor of G. graminis var. avenae and var. tritici after it had branched off from var. graminis. Intron AT and intron A exhibited vertical inheritance and coevolved in concert with their hosts. Intron loss could occur after its acquisition. Loss of any one of the three introns could occur in var. tritici whereas only loss of intron T had been found in the majority of var. avenae isolates. The existence of isolates of var. tritici and var. avenae with three introns suggested that intron loss could be reversed by intron acquisition and that the whole process is a dynamic one. This process of intron acquisition and intron loss reached different equilibrium points for different varieties and subgroups, which explained the irregular distribution of these introns in G. graminis. Each of the three group I introns was more closely related to other intron sequences that share the same insertion point in the 26S rDNA than to each other. These introns in distantly related organisms appeared to have a common ancestry. This system had provided a good model for studies on both the lateral transfer and common ancestry of group I introns in the 26S rRNA genes. Received: 17 May 1996 / Accepted: 14 January 1997     

Evidence of the presence of tRNAfMet group I intron in the marine cyanobacterium Synechococcus elongatus

Self-splicing group I introns in tRNA anticodon loops have been found in diverse groups of bacteria. In this work, we identified tRNAfMet group I introns in six strains of marine Synechococcus elongatus. Introns with sizes around 280 bp were consistently obtained in all the strains tested. In a phylogenetic analysis using the nucleotide sequence determined in this study with other cyanobacterial tRNAfMet and tRNALeu intron sequences, the Synechococcus sequence was grouped together with the sequences from other unicellular cyanobacterial strains. Interestingly, the phylogenetic tree inferred from the intronic sequences clearly separates the different tRNA introns, suggesting that each family has its own evolutionary history.     

Group I Introns from Zygomycota: Evolutionary Implications for the Fungal IC1 Intron Subgroup

The origins of fungal group I introns within nuclear small-subunit (nSSU) rDNA are enigmatic. This is partly because they have never been reported in basal fungal phyla (Zygomycota and Chytridiomycota), which are hypothesized to be ancestral to derived phyla (Ascomycota and Basidiomycota). Here we report group I introns from the nSSU rDNA of two zygomycete fungi, Zoophagus insidians (Zoopagales) and Coemansia mojavensis (Kickxellales). Secondary structure analyses predicted that both introns belong to the IC1 subgroup and that they are distantly related to each other, which is also suggested by different insertion sites. Molecular phylogenetic analyses indicated that the IC1 intron of Z. insidians is closely related to the IC1 intron inserted in the LSU rDNA of the basidiomycete fungus Clavicorona taxophila, which strongly suggests interphylum horizontal transfer. The IC1 intron of C. mojavensis has a low phylogenetic affinity to other fungal IC1 introns inserted into site 943 of nSSU rDNA (relative to E. coli 16S rDNA). It is noteworthy that this intron contains a putative ORF containing a His–Cys box motif in the antisense strand, a hallmark for nuclear-encoded homing endonucleases. Overall, molecular phylogenetic analyses do not support the placement of these two introns in basal fungal IC1 intron lineages. This result leads to the suggestion that fungal IC1 introns might have invaded or been transferred laterally after the divergence of the four major fungal phyla. Received: 8 February 2001 / Accepted: 1 November 2001     

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.     

Divergent Histories of rDNA Group I Introns in the Lichen Family Physciaceae

The wide but sporadic distribution of group I introns in protists, plants, and fungi, as well as in eubacteria, likely resulted from extensive lateral transfer followed by differential loss. The extent of horizontal transfer of group I introns can potentially be determined by examining closely related species or genera. We used a phylogenetic approach with a large data set (including 62 novel large subunit [LSU] rRNA group I introns) to study intron movement within the monophyletic lichen family Physciaceae. Our results show five cases of horizontal transfer into homologous sites between species but do not support transposition into ectopic sites. This is in contrast to previous work with Physciaceae small subunit (SSU) rDNA group I introns where strong support was found for multiple ectopic transpositions. This difference in the apparent number of ectopic intron movements between SSU and LSU rDNA genes may in part be explained by a larger number of positions in the SSU rRNA, which can support the insertion and/or retention of group I introns. In contrast, we suggest that the LSU rRNA may have fewer acceptable positions and therefore intron spread is limited in this gene. Reviewing Editor: Dr. W. Ford Doolittle     

Phylogenetic evidence for horizontal transmission of group I introns in the nuclear ribosomal DNA of mushroom-forming fungi

Group I introns were discovered inserted at the same position in the nuclear small-subunit ribosomal DNA (nuc-ssu-rDNA) in several species of homobasidiomycetes (mushroom-forming fungi). Based on conserved intron sequences, a pair of intron-specific primers was designed for PCR amplification and sequencing of intron-containing rDNA repeats. Using the intron-specific primers together with flanking rDNA primers, a PCR assay was conducted to determine presence or absence of introns in 39 species of homobasidiomycetes. Introns were confined to the genera Panellus, Clavicorona, and Lentinellus. Phylogenetic analyses of nuc-ssu-rDNA and mitochondrial ssu-rDNA sequences suggest that Clavicorona and Lentinellus are closely related, but that Panellus is not closely related to these. The simplest explanation for the distribution of the introns is that they have been twice independently gained via horizontal transmission, once on the lineage leading to Panellus, and once on the lineage leading to Lentinellus and Clavicorona. BLAST searches using the introns from Panellus and Lentinellus as query sequences retrieved 16 other similar group I introns of nuc-ssu-rDNA and nuclear large-subunit rDNA (nuc-lsu-rDNA) from fungal and green algal hosts. Phylogenetic analyses of intron sequences suggest that the mushroom introns are monophyletic, and are nested within a clade that contains four other introns that insert at the same position as the mushroom introns, two from different groups of fungi and two from green algae. The distribution of host lineages and insertion sites among the introns suggests that horizontal and vertical transmission, homing, and transposition have been factors in intron evolution. As distinctive, heritable features of nuclear rDNAs in certain lineages, group I introns have promise as phylogenetic markers. Nevertheless, the possibility of horizontal transmission and homing also suggest that their use poses certain pitfalls.      

Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution

Group I and group II introns are different catalytic self-splicing and mobile RNA elements that contribute to genome dynamics. In this study, we have analyzed their distribution and evolution in 29 sequenced genomes from the Bacillus cereus group of bacteria. Introns were of different structural classes and evolutionary origins, and a large number of nearly identical elements are shared between multiple strains of different sources, suggesting recent lateral transfers and/or that introns are under a strong selection pressure. Altogether, 73 group I introns were identified, inserted in essential genes from the chromosome or newly described prophages, including the first elements found within phages in bacterial plasmids. Notably, bacteriophages are an important source for spreading group I introns between strains. Furthermore, 77 group II introns were found within a diverse set of chromosomal and plasmidic genes. Unusual findings include elements located within conserved DNA metabolism and repair genes and one intron inserted within a novel retroelement. Group II introns are mainly disseminated via plasmids and can subsequently invade the host genome, in particular by coupling mobility with host cell replication. This study reveals a very high diversity and variability of mobile introns in B. cereus group strains.     

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.     

Introns and their flanking sequences of Bombyx mori rDNA.

We obtained two different clones (16 kb and 13 kb) of B. mori rDNA with intron sequence within the 28S-rRNA coding region. The sequence surrounding the intron was found to be highly conserved as indicated in several eukaryotes (Tetrahymena, Drosophila and Xenopus). The 28S rRNA-coding sequence of 16 kb and 13 kb clone was interrupted at precisely the same sites as those where the D. melanogaster rDNA interrupted by the type I and type II intron, respectively. The intron sequences of B. mori were different from those of D. melanogaster. In 16 kb clone, the intron was flanked by 14 bp duplication of the junction sequence, which was also present once within the 28S rRNA-coding region of rDNA without intron. This 14 bp sequence was identical with those surrounding the introns of Dipteran rDNAs.     

Detection of Seven Major Evolutionary Lineages in Cyanobacteria Based on the 16S rRNA Gene Sequence Analysis with New Sequences of Five Marine Synechococcus Strains

Although molecular phylogenetic studies of cyanobacteria on the basis of the 16S rRNA gene sequence have been reported, the topologies were unstable, especially in the inner branchings. Our analysis of 16S rRNA gene phylogeny by the maximum-likelihood and neighbor-joining methods combined with rate homogeneous and heterogeneous models revealed seven major evolutionary lineages of the cyanobacteria, including prochlorophycean organisms. These seven lineages are always stable on any combination of these methods and models, fundamentally corresponding to phylogenetic relationships based on other genes, e.g., psbA, rbcL, rnpB, rpoC, and tufA. Moreover, although known genotypic and phenotypic characters sometimes appear paralleled in independent lineages, many characters are not contradictory within each group. Therefore we propose seven evolutionary groups as a working hypothesis for successive taxonomic reconstruction. New 16S rRNA sequences of five unicellular cyanobacterial strains, PCC 7001, PCC 7003, PCC 73109, PCC 7117, and PCC 7335 of Synechococcus sp., were determined in this study. Although all these strains have been assigned to ``marine clusters B and C,' they were separated into three lineages. This suggests that the organisms classified in the genus Synechococcus evolved diversely and should be reclassified in several independent taxonomic units. Moreover, Synechococcus strains and filamentous cyanobacteria make a monophyletic group supported by a comparatively high statistical confidence value (80 to 100%) in each of the two independent lineages; therefore, these monophylies probably reflect the convergent evolution of a multicellular organization. Received: 3 September 1998 / Accepted: 30 November 1998     

Remarkable interkingdom conservation of intron positions and massive,lineage-specific intron loss and gain in eukaryotic evolution

Sequencing of eukaryotic genomes allows one to address major evolutionary problems, such as the evolution of gene structure. We compared the intron positions in 684 orthologous gene sets from 8 complete genomes of animals, plants, fungi, and protists and constructed parsimonious scenarios of evolution of the exon-intron structure for the respective genes. Approximately one-third of the introns in the malaria parasite Plasmodium falciparum are shared with at least one crown group eukaryote; this number indicates that these introns have been conserved through >1.5 billion years of evolution that separate Plasmodium from the crown group. Paradoxically, humans share many more introns with the plant Arabidopsis thaliana than with the fly or nematode. The inferred evolutionary scenario holds that the common ancestor of Plasmodium and the crown group and, especially, the common ancestor of animals, plants, and fungi had numerous introns. Most of these ancestral introns, which are retained in the genomes of vertebrates and plants, have been lost in fungi, nematodes, arthropods, and probably Plasmodium. In addition, numerous introns have been inserted into vertebrate and plant genes, whereas, in other lineages, intron gain was much less prominent.     

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.     

Structural features and evolutionary considerations of group IB introns in SSU rDNA of the lichen fungus Teloschistes

Different species of the lichen-forming ascomycete fungus Teloschistes were found to contain group IB introns at position S1506 in the small subunit ribosomal RNA gene. We have characterized the structural organization and phylogeny of the Teloschistes introns Tco.S1506, Tla.S1506, and Tvi.S1506. Common features to all the introns are a small size, a compact RNA structure, and an atypical catalytic ribozyme core sequence motif. Variations in intron sizes, due to sequence extensions in the P1 and P8 loop segments, were observed in different species and isolates. Phylogenetic analyses based on the ITS1-5.8S-ITS2 region as well as the introns show that the Teloschistes S1506 introns represent a distinct evolutionary isolated cluster among the nuclear group I introns. Furthermore, introns from different lineages of Teloschistes villosus appear not strictly vertically inherited probably due to horizontal transfer in one of the lineages.     

USING RDNA GENES TO UNDERSTAND THE ORIGIN AND EVOLUTION OF SPLICEOSOMAL AND GROUP I INTRONS

Ribosomal DNA genes in lichen algae and lichen fungi are astonishingly rich in spliceosomal and group I introns. We use phylogenetic, secondary structure, and biochemical analyses to understand the evolution of these introns. Despite the widespread distribution of spliceosomal introns in nuclear pre-mRNA genes, their general mechanism of origin remains an open question because few proven cases of recent and pervasive intron origin have been documented. The lichen introns are valuable in this respect because they are undoubtedly of a "recent" origin and limited to the Euascomycetes. Our analyses suggest that rDNA spliceosomal introns have arisen through aberrant reverse-splicing (in trans) of free pre-mRNA introns into r RNAs. We propose that the spliceosome itself (and not an external agent; e.g. transposable elements, group II introns) has given rise to the introns. The rDNA introns are found most often between the flanking sequence G (78%) - intron-G (72%), and their clustered positions on secondary structures suggest that particular r RNA regions are preferred sites (i.e., proto-splice sites) for insertion. Mapping of intron positions on the newly available tertiary structures show that they are found most often in exposed regions of the ribosomes. This again is consistent with an intron origin through reverse-splicing. Remarkably, the distribution and phylogenetic relationships of most group I introns in nuclear rDNA genes are also consistent with a reverse-splicing origin. These data underline the value of lichens as a model system for understanding intron origin and stress the importance of RNA-level processes in the spread of these sequences in nuclear coding regions.