Sequence comparison of the rDNA introns from six different species of Tetrahymena

Intraspecies Variation in Tetrahymena rostrata

Two distinct isolates of the facultative parasite, Tetrahymena rostrata were compared, identifying and utilising markers that are useful for studying clonal variation within the species were identified and utilised. The sequences of mitochondrial genomes and several nuclear genes were determined using Illumina short read sequencing. The two T. rostrata isolates had similar morphology. The linear mitogenomes had the gene content and organisation typical of the Tetrahymena genus, comprising 8 tRNA genes, 6 ribosomal RNA genes and 45 protein coding sequences (CDS), twenty-two of which had known function. The two isolates had nucleotide identity within common nuclear markers encoded within the histone H3 and H4 and small subunit ribosomal RNA genes and differed by only 2–4 nucleotides in a region of the characterised actin genes. Variation was observed in several mitochondrial genes and was used to determine intraspecies variation and may reflect the natural history of T. rostrata from different hosts or the geographic origins of the isolates.

Phylogenetic analyses suggest reverse splicing spread of group I introns in fungal ribosomal DNA

BACKGROUND

Group I introns have spread into over 90 different sites in nuclear ribosomal DNA (rDNA) with greater than 1700 introns reported in these genes. These ribozymes generally spread through endonuclease-mediated intron homing. Another putative pathway is reverse splicing whereby a free group I intron inserts into a homologous or heterologous RNA through complementary base-pairing between the intron and exon RNA. Reverse-transcription of the RNA followed by general recombination results in intron spread. Here we used phylogenetics to test for reverse splicing spread in a taxonomically broadly sampled data set of fungal group I introns including 9 putatively ancient group I introns in the rDNA of the yeast-like symbiont Symbiotaphrina buchneri.

RESULTS

Our analyses reveal a complex evolutionary history of the fungal introns with many cases of vertical inheritance (putatively for the 9 introns in S. buchneri) and intron lateral transfer. There are several examples in which introns, many of which are still present in S. buchneri, may have spread through reverse splicing into heterologous rDNA sites. If the S. buchneri introns are ancient as we postulate, then group I intron loss was widespread in fungal rDNA evolution.

CONCLUSION

On the basis of these results, we suggest that the extensive distribution of fungal group I introns is at least partially explained by the reverse splicing movement of existing introns into ectopic rDNA sites.

Multiple group I introns detected in the nuclear small subunit rDNA of the autosporic green alga Selenastrum capricornutum

A phylogenetic investigation of the autosporic chlorophycean alga species Selenastrum capricornutum using the small subunit (SSU) rRNA gene revealed the unusual presence of six group IC1 introns. Previous studies showed that numerous green algal taxa contain group IC1 introns. However, whereas some algal species harbor multiple introns in a single ribosomal gene, none have contained as many as S. capricornutum. Three of the S. capricornutum introns are located at conserved algal intron sites and the remaining three are located at novel eukaryotic positions. The SSU rRNA genes and their introns have been sequenced and putative secondary structures are proposed for the introns. Also, their similarity to other group IC1 introns from algal, fungal, and viral sources is investigated. Results suggest the initial presence of introns at conserved locations, followed by duplication and insertion to novel positions within the SSU rRNA gene.

Inter- and intra-strain variation in the 5.8S ribosomal RNA and internal transcribed spacer sequences of Entamoeba histolytica and comparison with Entamoeba dispar, Entamoeba moshkovskii and Entamoeba invadens

The ribosomal RNA genes in Entamoeba histolytica are located on circular DNA molecules in about 200 copies per genome equivalent. Nucleotide sequence analysis of the 5.8S rRNA gene and the flanking internal transcribed spacers was carried out to determine the degree of sequence divergence in the multiple rRNA gene copies of a given strain; amongst three different E. histolytica strains (HM-1:IMSS, Rahman and HK-9); and amongst four species of Entamoeba (Entamoeba histolytica, Entamoeba dispar, Entamoeba moshkovskii and Entamoeba invadens). The results show that all rRNA gene copies of a given strain are identical. Few nucleotide positions varied between strains of a species but the differences were very pronounced amongst species. In general, the internal transcribed spacer 2 sequence was more variable and may be useful for strain- and species-identification. The 5.8S rRNA gene and the internal transcribed spacer 2 of E. invadens were unusually small in size.

Three different group I introns in the nuclear large subunit ribosomal DNA of the amoeboflagellate Naegleria

We have amplified the large subunit ribosomal DNA (LSUrDNA) of the 12 described Naegleria spp. and of 34 other Naegleria lineages that might be distinct species. Two strains yielded a product that is longer than 3 kb, which is the length of the LSUrDNA of all described Naegleria spp. Sequencing data revealed that the insert in one of these strains is a group I intron without an open reading frame (ORF), while the other strain contains two different group I introns, of which the second intron has an ORF of 175 amino acids. In the latter ORF there is a conserved His-Cys box, as in the homing endonucleases present in group I introns in the small subunit ribosomal DNA (SSUrDNA) of Naegleria spp. Although the group I introns in the LSUrDNA differ in sequence, they are more related to each other than they are to the group I introns in the SSUrDNA of Naegleria spp. The three group I introns in the LSUrDNA in Naegleria are at different locations and are probably acquired by horizontal transfer, contrary to the SSUrDNA group I introns in this genus which are of ancestral origin and are transmitted vertically.

Comparison of sequence differences in a variable 23S rRNA domain among sets of cryptic species of ciliated protozoa

Studies were undertaken to discover the relative molecular distances separating some familiar forms of ciliated protozoa, and the genetic species they include. Sequences of 190 bases of the D2 domain of the large ribosomal nucleic acid molecule were obtained by polymerase chain reaction from protists of three distinctive groups of ciliated protozoa-Colpoda, Paramecium and Tetrahymena. Evolutionary trees were constructed for each set of sequences using the PHYLOGEN 1.0 string programs. All three groups of ciliates manifested large molecular diversity among strains difficult or impossible to distinguish morphologically. The largest single evolutionary distance within a group was the 75 differences separating Tetrahymena paravorax from the other tetrahymenids. The largest mean distance for a group was the 21.2 for the colpodids. In all the protist groups the large molecular diversity is obscured by morphological conservatism associated with constraints of ancient designs. The molecular diversity within morphotypes argues for long evolutionary coexistence of species differentiated from each other in significant physiological, ecological, or nutritional ways.

Characterization of a flatworm ribosomal RNA-encoding gene: promoter sequence and small subunit rRNA secondary structure

The transcription start point (tsp) of a ribosomal RNA (rRNA)-encoding gene (rDNA) from Echinococcus granulosus has been mapped at a position located 1.1 kb upstream from the small subunit (SSU) rRNA coding sequence. As expected from the analysis of the putative promoter sequence (-200 to +50), no homology was found with rDNA promoters from other organisms. Nevertheless, some interesting motifs were found. There is a d(T)11 track flanked by an inverted repeat (IR) centered at position -32, which is analogous to the position of the TATA box in promoters transcribed by RNA polymerase II. Two other IR were observed that are also present in the Trypanosoma cruzi rDNA promoter. We have also determined the SSU rDNA sequence and established a secondary structure model. The analysis of the secondary structure strongly suggests that this gene encodes a functional SSU rRNA. The fact that both the promoter and the rRNA coding sequence are derived from the same rDNA repeat indicates that the promoter is also functional.

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.

Evidence for the ancestral origin of group I introns in the SSUrDNA of Naegleria spp

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

Variation and in vitro splicing of group I introns in rRNA genes of Pneumocystis carinii

The sequences of the rRNA genes of Pneumocystis carinii from rat and human sources demonstrate three distinct genotypes based on the group I introns present in these genes. One rat isolate (Pc1) contains such introns in its 16S and 26S rRNA genes, while another rat isolate (Pc2) and a human isolate (Pc3) only contain an intron in the 26S rRNA gene. The four introns all catalyze their own excision from RNA transcripts, and this reaction is inhibited by the anti-pneumocystis drug pentamidine and its analogues. Although they differ in sequence, they are more similar to one another than to group I introns found in other eukaryotic microbes.

Structure and evolution of myxomycete nuclear group I introns: a model for horizontal transfer by intron homing

We have examined five nuclear group I introns, located at three different positions in the large subunit ribosomal RNA (LSU rRNA) gene of the two myxomycete species, Didymium iridis and Physarum polycephalum. Structural models of intron RNAs, including secondary and tertiary interactions, are proposed. This analysis revealed that the Physarum intron 2 contains an unusual core region that lacks the P8 segment, as well as several of the base-triples known to be conserved among group I introns. Structural and evolutionary comparisons suggest that the corresponding introns 1 and 2 were present in a common ancestor of Didymium and Physarum, and that the five introns in LSU rRNA genes of these myxomycetes were acquired in three different events. Evolutionary relationships, inferred from the sequence analysis of several different nuclear group I introns and the ribosomal RNA genes of the intron-harbouring organisms, strongly support horizontal transfer of introns in the course of evolution. We propose a model that may explain how myxomycetes in natural environments obtained their nuclear group I introns.

Sequence and variability of the 5.8S and 26S rRNA genes of Pneumocystis carinii

The sequence of the coding region of the rRNA operon of rat-derived Pneumocystis carinii has been completed, including the genes for 5.8S and 26S rRNA. These genes show homology to the rRNA genes of yeast, and an apparent group I self-splicing intron is present in the 26S rRNA gene. Like a similar intron in the 16S rRNA gene, this intron is in a phylogenetically conserved region. Variation in the 26S rRNA sequence was noted between P. carinii organisms isolated from two different sources.

The gene coding for small ribosomal subunit RNA in the basidiomycete Ustilago maydis contains a group I intron

The nucleotide sequence of the gene coding for small ribosomal subunit RNA in the basidiomycete Ustilago maydis was determined. It revealed the presence of a group I intron with a length of 411 nucleotides. This is the third occurrence of such an intron discovered in a small subunit rRNA gene encoded by a eukaryotic nuclear genome. The other two occurrences are in Pneumocystis carinii, a fungus of uncertain taxonomic status, and Ankistrodesmus stipitatus, a green alga. The nucleotides of the conserved core structure of 101 group I intron sequences present in different genes and genome types were aligned and their evolutionary relatedness was examined. This revealed a cluster including all group I introns hitherto found in eukaryotic nuclear genes coding for small and large subunit rRNAs. A secondary structure model was designed for the area of the Ustilago maydis small ribosomal subunit RNA precursor where the intron is situated. It shows that the internal guide sequence pairing with the intron boundaries fits between two helices of the small subunit rRNA, and that minimal rearrangement of base pairs suffices to achieve the definitive secondary structure of the 18S rRNA upon splicing.

Inheritance of the group I rDNA intron in Tetrahymena pigmentosa

We have previously argued from phylogenetic sequence data that the group I intron in the rRNA genes of Tetrahymena was acquired by different Tetrahymena species at different times during evolution. We have now approached the question of intron mobility experimentally by crossing intron+ and intron- strains looking for a strong polarity in the inheritance of the intron (intron homing). Based on the genetic analysis we find that the intron in T. pigmentosa is inherited as a neutral character and that intron+ and intron- alleles segregate in a Mendelian fashion with no sign of intron homing. In an analysis of vegetatively growing cells containing intron+ and intron- rDNA, initially in the same macronucleus, we similarly find no evidence of intron homing. During the course of this work, we observed to our surprise that progeny clones from some crosses contained three types of rDNA. One possible explanation is that T. pigmentosa has two rdn loci in contrast to the single locus found in T. thermophila. Some of the progeny clones from the genetic analysis were expanded for several hundred generations, and allelic assortment of the rDNA was demonstrated by subcloning analysis.

Group I introns as mobile genetic elements: facts and mechanistic speculations--a review

Group I introns form a structural and functional group of introns with widespread but irregular distribution among very diverse organisms and genetic systems. Evidence is now accumulating that several group I introns are mobile genetic elements with properties similar to those originally described for the omega system of Saccharomyces cerevisiae: mobile group I introns encode sequence-specific double-strand (ds) endoDNases, which recognize and cleave intronless genes to insert a copy of the intron by a ds-break repair mechanism. This mechanism results in: the efficient propagation of group I introns into their cognate sites; their maintenance at the site against spontaneous loss; and, perhaps, their transposition to different sites. The spontaneous loss of group I introns occurs with low frequency by an RNA-mediated mechanism. This mechanism eliminates introns defective for mobility and/or for RNA splicing. Mechanisms of intron acquisition and intron loss must create an equilibrium, which explains the irregular distribution of group I introns in various genetic systems. Furthermore, the observed distribution also predicts that horizontal transfer of intron sequences must occur between unrelated species, using vectors yet to be discovered.

On finding all suboptimal foldings of an RNA molecule

An algorithm and a computer program have been prepared for determining RNA secondary structures within any prescribed increment of the computed global minimum free energy. The mathematical problem of determining how well defined a minimum energy folding is can now be solved. All predicted base pairs that can participate in suboptimal structures may be displayed and analyzed graphically. Representative suboptimal foldings are generated by selecting these base pairs one at a time and computing the best foldings that contain them. A distance criterion that ensures that no two structures are "too close" is used to avoid multiple generation of similar structures. Thermodynamic parameters, including free-energy increments for single-base stacking at the ends of helices and for terminal mismatched pairs in interior and hairpin loops, are incorporated into the underlying folding model of the above algorithm.

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

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

Three-dimensional model of the active site of the self-splicing rRNA precursor of Tetrahymena

The rRNA intervening sequence of Tetrahymena is a catalytic RNA molecule, or "ribozyme." A tertiary-structure model of the active site of this ribozyme has been constructed based on comparative sequence analysis of related group I intervening sequences, data on the accessibility of each nucleotide to chemical and enzymatic probes, and principles of RNA folding derived from a consideration of the structure of tRNA determined by x-ray crystallography. In the model, the catalytic center has a two-helix structural framework composed of the base-paired segments of the group I conserved sequence elements. The structural framework supports and orients the conserved nucleotides that are adjacent to the base-paired sequence elements; these conserved nucleotides are proposed to form the active site and to bind the 5' splice-site duplex and the guanine nucleotide substrate. Tests of the model are proposed.

Variable occurrence of the nrdB intron in the T-even phages suggests intron mobility

The bacteriophage T4 nrdB gene, encoding nucleoside diphosphate reductase subunit B, contains a self-splicing group I intervening sequence. The nrdB intron was shown to be absent from the genomes of the closely related T-even phages T2 and T6. Evidence for variable intron distribution was provided by autocatalytic 32P-guanosine 5'-triphosphate labeling of T-even RNAs, DNA and RNA hybridization analyses, and DNA sequencing studies. The results indicate the nonessential nature of the intron in nrdB expression and phage viability. Furthermore, they suggest that either precise intron loss from T2 and T6 or lateral intron acquisition by T4 occurred since the evolution of these phages from a common ancestor. Intron movement in the course of T-even phage divergence raises provocative questions about the origin of these self-splicing elements in prokaryotes.

Phylogenetic evidence for the acquisition of ribosomal RNA introns subsequent to the divergence of some of the major Tetrahymena groups

Previous work has demonstrated the presence of a self-splicing intron in the large subunit ribosomal RNA coding region in some strains of the ciliate protozoan Tetrahymena. Sequence comparisons of the intron regions from six Tetrahymena species showed these to fall into three homology groups. In an attempt to evaluate the evolutionary origins of the intervening sequences, we have now determined complete small subunit ribosomal RNA gene sequences from 13 species of Tetrahymena and the absolute number of nucleotide differences between the sequences was used to construct a phylogenetic tree. This phylogeny was consistent with the groupings suggested by comparisons of other biochemical characters including cytoskeletal proteins, isozyme analyses, and restriction maps of complete rRNA transcription units. The homology groupings that were based upon the intron sequence data do not agree with the relationships inferred from the small subunit rRNA sequence data. These observations are taken to indicate that the intron character has been acquired independently in different species at a stage later than the branching out of the species.