Paramecium mitochondrial genes. I. Small subunit rRNA gene sequence and microevolution

The Oxytricha trifallax mitochondrial genome

The Oxytricha trifallax mitochondrial genome contains the largest sequenced ciliate mitochondrial chromosome (~70 kb) plus a ~5-kb linear plasmid bearing mitochondrial telomeres. We identify two new ciliate split genes (rps3 and nad2) as well as four new mitochondrial genes (ribosomal small subunit protein genes: rps- 2, 7, 8, 10), previously undetected in ciliates due to their extreme divergence. The increased size of the Oxytricha mitochondrial genome relative to other ciliates is primarily a consequence of terminal expansions, rather than the retention of ancestral mitochondrial genes. Successive segmental duplications, visible in one of the two Oxytricha mitochondrial subterminal regions, appear to have contributed to the genome expansion. Consistent with pseudogene formation and decay, the subtermini possess shorter, more loosely packed open reading frames than the remainder of the genome. The mitochondrial plasmid shares a 251-bp region with 82% identity to the mitochondrial chromosome, suggesting that it most likely integrated into the chromosome at least once. This region on the chromosome is also close to the end of the most terminal member of a series of duplications, hinting at a possible association between the plasmid and the duplications. The presence of mitochondrial telomeres on the mitochondrial plasmid suggests that such plasmids may be a vehicle for lateral transfer of telomeric sequences between mitochondrial genomes. We conjecture that the extreme divergence observed in ciliate mitochondrial genomes may be due, in part, to repeated invasions by relatively error-prone DNA polymerase-bearing mobile elements.

Intraspecific genetic variation in Paramecium revealed by mitochondrial cytochrome C oxidase I sequences

Studies of intraspecific genetic diversity of ciliates, such as population genetics and biogeography, are particularly hampered by the lack of suitable DNA markers. For example, sequences of the non-coding ribosomal internal transcribed spacer (ITS) regions are often too conserved for intraspecific analyses. We have therefore identified primers for the mitochondrial cytochrome c oxidase I (COI) gene and applied them for intraspecific investigations in Paramecium caudatum and Paramecium multimicronucleatum. Furthermore, we obtained sequences of the ITS regions from the same strains and carried out comparative sequence analyses of both data sets. The mitochondrial sequences revealed substantially higher variation in both Paramecium species, with intraspecific divergences up to 7% in P. caudatum and 9.5% in P. multimicronucleatum. Moreover, an initial survey of the population structure discovered different mitochondrial haplotypes of P. caudatum in one pond, thereby demonstrating the potential of this genetic marker for population genetic analyses. Our primers successfully amplified the COI gene of other Paramecium. This is the first report of intraspecific variation in free-living protozoans based on mitochondrial sequence data. Our results show that the high variation in mitochondrial DNA makes it a suitable marker for intraspecific and population genetic studies.

Mitochondria of Protists

Over the past several decades, our knowledge of the origin and evolution of mitochondria has been greatly advanced by determination of complete mitochondrial genome sequences. Among the most informative mitochondrial genomes have been those of protists (primarily unicellular eukaryotes), some of which harbor the most gene-rich and most eubacteria-like mitochondrial DNAs (mtDNAs) known. Comparison of mtDNA sequence data has provided insights into the radically diverse trends in mitochondrial genome evolution exhibited by different phylogenetically coherent groupings of eukaryotes, and has allowed us to pinpoint specific protist relatives of the multicellular eukaryotic lineages (animals, plants, and fungi). This comparative genomics approach has also revealed unique and fascinating aspects of mitochondrial gene expression, highlighting the mitochondrion as an evolutionary playground par excellence.

Identifying and distinguishing sibling species in the Tetrahymena pyriformis complex (Ciliophora, Oligohymenophorea) using PCR/RFLP analysis of nuclear ribosomal DNA

We describe a riboprinting strategy for identifying and distinguishing among sibling species in the Tetrahymena pyriformis complex. It involves use of the polymerase chain reaction to amplify a large segment of the nuclear ribosomal DNA and internal transcribed spacers, and digestion of this DNA with restriction enzymes. Unique restriction fragment length patterns or haplotypes were then used to distinguish species into: (1) six taxa that were identifiable to the species level, (2) eight taxa that were separated into four pairs, and (3) a group of eight taxa that were identical to each other. The latter result indicates that a more variable molecule is needed to distinguish the most closely related species in the complex. There was no intraspecific variation between two strains from one species (Tetrahymena thermophila) nor among multiple isolates from another species (Tetrahymena empidokyrea). This approach provides an alternative to traditional techniques for identifying T. pyriformis species that require living reference specimens and/or that reveal high levels of intraspecific variation.

An organelle-like small subunit ribosomal RNA gene from Babesia bovis: nucleotide sequence, secondary structure of the transcript and preliminary phylogenetic analysis

Investigations aimed at identifying the mitochondrial genome of Babesia bovis using the polymerase chain reaction (PCR) have established the existence of an organelle-like small subunit ribosomal RNA (SSU rRNA) gene in the parasite. The sequence, compiled from three main PCR products, was 1448 bp in length (including the primer regions), had a 73% A+T content and showed significant similarity (68% sequence identity) to the "organellar" SSU rRNA gene from Plasmodium falciparum. The proposed secondary structure of the transcript showed several features which were consistent with a eubacterial origin for the organelle-like gene. The presence of putative binding sites for streptomycin and tetracycline also supported an "organellar" location for the gene and suggested that the SSU rRNA transcript is functional in protein synthesis because tetracycline has anti-babesial activity. Phylogenetic analyses based on the conserved regions of the SSU-like rRNA genes from a wide variety of organisms showed only a weak association of the babesial sequence with its mitochondrial homologues and an even weaker association with the corresponding genes of plastid origin. The origin of this organelle-like gene in B. bovis therefore remains unresolved, as is the case for its homologue from P. falciparum.

The mitochondrial ribosomal RNA genes of the nematodes Caenorhabditis elegans and Ascaris suum: consensus secondary-structure models and conserved nucleotide sets for phylogenetic analysis

The small- and large-subunit mitochondrial ribosomal RNA genes (mt-s-rRNA and mt-l-rRNA) of the nematode worms Caenorhabditis elegans and Ascaris suum encode the smallest rRNAs so far reported for metazoa. These size reductions correlate with the previously described, smaller, structurally anomalous mt-tRNAs of C. elegans and A. suum. Using primer extension analysis, the 5' end nucleotides of the mt-s-rRNA and mt-l-rRNA genes were determined to be adjacent to the 3' end nucleotides of the tRNA(Glu) and tRNA(His) genes, respectively. Detailed, consensus secondary-structure models were constructed for the mt-s-rRNA genes and the 3' 64% of mt-l-rRNA genes of the two nematodes. The mt-s-rRNA secondary-structure model bears a remarkable resemblance to the previously defined universal core structure of E. coli 16S rRNA: most of the nucleotides that have been classified as variable or semiconserved in the E. coli model appear to have been eliminated from the C. elegans and A. suum sequences. Also, the secondary structure model constructed for the 3' 64% of the mt-l-rRNA is similar to the corresponding portion of the previously defined E. coli 23S rRNA core secondary structure. The proposed C. elegans/A. suum mt-s-rRNA and mt-l-rRNA models include all of the secondary-structure element-forming sequences that in E. coli rRNAs contain nucleotides important for A-site and P-site (but not E-site) interactions with tRNAs. Sets of apparently homologous sequences within the mt-s-rRNA and mt-l-rRNA core structures, derived by alignment of the C. elegans and A. suum mt-rRNAs to the corresponding mt-rRNAs of other eukaryotes, and E. coli rRNAs were used in maximum-likelihood analyses. The patterns of divergence of metazoan phyla obtained show considerable agreement with the most prevalent metazoan divergence patterns derived from more classical, morphological, and developmental data.

Mitochondrial DNA sequence analysis of the cytochrome oxidase subunit II gene from Podospora anserina. A group IA intron with a putative alternative splice site

A 5 kb region of the 95 kb mitochondrial genome of Podospora anserina race s has been mapped and sequenced (1 kb = 10(3) base-pairs). This DNA region is continuous with the sequence for the ND4L and ND5 gene complex in the accompanying paper. We show that this sequence contains the gene for cytochrome oxidase subunit II (COII). This gene is 4 kb in length and is interrupted by a subgroup IB intron (1267 base-pairs (bp) in length) and a subgroup IA intron (1992 bp in length). This group IA intron has a long open reading frame (ORF; 472 amino acid residues) discontinuous with the upstream exon sequence. A putative alternative splice site is present, which brings the ORF into phase with the 5' exon sequence. The 5'- and 3'-flanking regions of the COII gene contain G + C-rich palindromic sequences that resemble similar sequences flanking many Neurospora crassa mitochondrial genes.

DNA sequence analysis of the mitochondrial ND4L-ND5 gene complex from Podospora anserina. Duplication of the ND4L gene within its intron

A 15 kb region of the 100 kb mitochondrial genome of Podospora anserina has been mapped and sequenced (1 kb = 10(3) base-pairs). The genes for ND4L and ND5 are identified as contiguous genes with overlapping termination and initiation codons. In race A (101 kb) the gene for ND4L (4.3 kb) has a gene duplication within an intron including a second subgroup IC intron. Race s (95 kb) lacks this second gene complex. Each intron has the identical 5' exon boundary. Secondary structure analysis showed that the closest relative of the second intron is the first intron itself. The open reading frames of the two introns are also closely related to each other as well as to their counterpart in the ND4L gene of Neurospora crassa. The 9.9 kb ND5 gene starts immediately at the termination codon of ND4L and is split by two group IB introns, one group IC intron and one group II intron. The group II intron is closely related to other group II introns although its open reading frame sequence similarity with retroviral reverse transcriptase appears to be more divergent. The similarities in secondary structure and open reading frames for these six introns are discussed.

DNA sequence analysis of the 24.5 kilobase pair cytochrome oxidase subunit I mitochondrial gene from Podospora anserina: a gene with sixteen introns

The DNA sequence of a 26.7 Kilobase pair (10(3) base pairs = 1 Kb) region of the mitochondrial genomes of races s and A from Podospora anserina was determined. Within this region, the 24.5 Kb cytochrome oxidase subunit I gene was located and its exon sequences determined by computer analysis comparisons with other fungal genes. The Podospora COI gene was interrupted by two group II introns (one in race s) and fourteen group I introns ranging in size from about 2.2 Kb to 404 bp. Earlier studies on secondary structure analysis, as well as comparison of their open reading frames (ORFs), showed that the two group II introns were closely related. The fourteen group I introns were representatives of three subgroupings (IB, C and a new category, subgroup ID). Two of these group I introns were separated by just a single exon codon. The analysis of all these introns is discussed in comparison with other fungal introns as well as with the known Podospora anserina introns.

An unusual region of Paramecium mitochondrial DNA containing chloroplast-like genes

Based on DNA and amino acid comparisons with known genes and their products, a region of the Paramecium aurelia mitochondrial (mt) genome has been found to encode the following gene products: (1) photosystem II protein G (psbG); (2) a large open reading frame (ORF400) which is also found encoded in the chloroplast (cp) DNA of tobacco (as ORF393) and liverwort (as ORF392), and in the kinetoplast maxicircle DNA of Leishmania tarentolae (as ORFs 3 and 4); (3) ribosomal protein L2 (rpl2); (4) ribosomal protein S12 (rps12); (5) ribosomal protein S14 (rps14); and (6) NADH dehydrogenase subunit 2 (ndh2). All of these genes have been found in cp DNA, but the psbG gene has never been identified in a mt genome, and ribosomal protein genes have never been located in an animal or protozoan mitochondrion. The ndh2 gene has been found in both mitochondria and plastids. The Paramecium genes are among the most divergent of those sequenced to date. Two of the genes are encoded on the strand of DNA complementary to that encoding all other known Paramecium mt genes. No gene contains an identifiable intron. The rps12 and psbG genes are probably overlapping. It is not yet known whether these genes are transcribed or have functional gene products. The presence of these genes in the mt genome raises interesting questions concerning their evolutionary origin.

DNA sequence, structure, and phylogenetic relationship of the small subunit rRNA coding region of mitochondrial DNA from Podospora anserina

DNA sequence analysis and the localization of the 5' and 3' termini by S1 mapping have shown that the mitochondrial (mt) small subunit rRNA coding region from Podospora anserina is 1980 bp in length. The analogous coding region for mt rRNA is 1962 bp in maize, 1686 bp in Saccharomyces cerevisiae, and 956 bp in mammals, whereas its counterpart in Escherichia coli is 1542 bp. The P. anserina mt 16S-like rRNA is 400 bases longer than that from E. coli, but can be folded into a similar secondary structure. The additional bases appear to be clustered at specific locations, including extensions at the 5' and 3' termini. Comparison with secondary structure diagrams of 16S-like RNAs from several organisms allowed us to specify highly conserved and variable regions of this gene. Phylogenetic tree construction indicated that this gene is grouped with other mitochondrial genes, but most closely, as expected, with the fungal mitochondrial genes.

Sequence analysis of mitochondrial DNA from Podospora anserina. Pervasiveness of a class I intron in three separate genes

A 48 kb region of the 95 kb mitochondrial genome of Podospora anserina has been mapped and sequenced (1 kb = 10(3) base-pairs). The DNA sequence of the genes for ND2, 3, 4, ATPase 6 and URFC are presented here. As in Neurospora crassa, the ND2 and 3 genes consist of a unit separated by one TAA stop codon. ND3, 4 and ATPase 6 are interrupted by class I introns. All three introns are remarkably similar in the C-domain of their secondary structure, sufficient enough to designate them as new subgroup, class IC introns. The open reading frames of the ND3 and 4 introns bear a high sequence similarity to the open reading frame of the class IB introns of ATPase 6 from N. crassa and ND1 from Neurospora intermedia Varkud. We also show that the tRNA Met-2 gene is duplicated and is involved in a recombinational event. The 5' region of URFC is also duplicated but no involvement of this gene with recombination or formation of plasmids is known. The evolutionary significance of the similarities of intron secondary structures and open reading frames of the ND3, 4 and ATPase 6 genes is discussed, including the possible separate evolution of structural and coding sequences.

DNA sequence and organization of the mitochondrial ND1 gene from Podospora anserina: analysis of alternate splice sites

Earlier, we reported that the ND1 mitochondrial gene of Podospora anserina is mosaic, containing at least three class I introns. We have now completed the sequence of the ND1 gene and have determined that it contains four class I introns of 1,820, 2,631, 2,256 and 2,597 bp with the entire gene complex containing 10,505 bp, only 1,101 of which are exon sequences. Introns 1 and 3 appear to be related in that their open reading frames (ORFS) exhibit extensive amino acid sequence similarity and like the URFN sequence from Neurospora crassa have multiple sequence repetitions. Introns 2 and 4 are similar in that both appear to be mosaic introns. Where intron 2 has many short ORFS, intron 4 has two, 391 and 262 aa respectively. The first ORF has some patch work sequence similarity with one of the intron 2 ORFs but the second ORF is strikingly similar to the single intron ORF in the ND1 gene of N. crassa. Just upstream of the sequences necessary to form the central core of the P. anserina intron 4 secondary structure, there is a 17 bp sequence which is an exact replica of the exon sequence abutting the 5' flank of the 1,118 bp N. crassa ND1 intron. Secondary structure analysis suggests that the 2,597 bp intron 4 can fold as an entity but a similar structure can be constructed just from an 1,130 bp portion by utilizing the 17 bp element as an alternate splice site. Detailed structural analysis suggests that intron 4 (as well as the single ND1 intron from N. crassa) can utilize helical configurations which bring the downstream open reading frame into juxtaposition with the exon sequences.

Unexpectedly long 14S ribosomal RNA gene in Tetrahymena mitochondria

Extraction of RNA from Tetrahymena mitochondrial ribosomes yields several RNA species, including a "large" 21S molecule, a "small" 14S molecule, a 7S molecule, and other smaller RNAs. The molecular weight of the 14S rRNA indicates that it is about 1,300 bases in length. We have sequenced the 14S rRNA gene and, by aligning our sequence with that of the corresponding small rRNA from E. coli, find that the 14S rDNA is at least 1,635 bases in length. We propose, based on the results of hybridization studies, that this unexpected length is due to the presence of 7S RNA sequence within the 14S gene sequence. The 7S region is apparently lost from the 14S rRNA, yet is still a component of the ribosome.

Paramecium mitochondrial DNA sequences and RNA transcripts for cytochrome oxidase subunit I, URF1, and three ORFs adjacent to the replication origin

A 2-kb region adjacent to the replication origin (ori) and a 3-kb region located between the small and large ribosomal RNAs of Paramecium mitochondrial (mt) DNA have been sequenced and the locations of their transcripts determined. The ori segment contains four transcripts, some of which are overlapping, which encode a known protein and two other open reading frames. The other segment encodes, on separate transcripts, the cytochrome c oxidase subunit one gene (COI) and the URF1 gene (ND1) common to most mt genomes. All these genes have the same orientation and do not contain introns. The COI gene is the most divergent of those known and has an internal 108 amino acid 'insert' not found in COI genes from other organisms. With these data it is possible to define a probable Paramecium mt genetic code. With the exception that TGA codes for tryptophan and the use of different start codons, Paramecium mtDNA appears to follow the universal code. GTA possibly can be used as a start codon.

Excision-amplification of mitochondrial DNA during senescence in Podospora anserina. DNA sequence analysis of three unique "plasmids"

During senescence in the filamentous fungus Podospora anserina, specific regions of the mitochondrial genome, termed senDNA are excised, ligated and amplified. We have cloned in their entirety three such autonomously replicating plasmids, alpha, beta and epsilon senDNA. None of these plasmids displayed cross-hybridization nor did we detect any significant DNA homology by computer analysis. The complete DNA sequence of the 2.5 kb alpha, the 5.5 kb epsilon and about 3.4 kb of the 9.8 kb beta senDNA is presented (kb = 10(3) base-pairs). These sequences were analyzed for the presence of consensus sequences common to introns, and it was found that alpha senDNA has the characteristics of a group II intron, epsilon senDNA contains three group I introns, and beta senDNA did not show relevant sequences in the 3.4 kb examined. Comparison of the 5' and 3'-flanking sequences of alpha senDNA with oxi 3 (Co I) amino acid sequences from Neurospora crassa and Saccharomyces cerevisiae revealed significant homology and provided strong support that the excised alpha senDNA itself consists entirely of an intron. Upstream from the oxi 3 gene a transfer RNA cysteine sequence was detected. beta senDNA contained four tRNA sequences, aspartic acid, serine, valine and tryptophan, and sequences homologous to URFC (untranslated reading frame C) as well as two new URFs. epsilon senDNA contained sequences homologous to ATPase 8 and URFl; URFl was interrupted by three group I introns. The excision site sequences, as located by S1 nuclease mapping were unique for each senDNA. Analysis for repeated units showed that each plasmid contained elements which could be involved in secondary structure required for the alignment of distal ends preparatory to excision. These results are interpreted in terms of the structural requirements of mobile elements including the possible involvement of reverse transcriptase in the excision-ligation-amplification process.

A fine restriction map of the linear mitochondrial DNA of Tetrahymena pyriformis: genome size, map locations of rRNA and tRNA genes, terminal inversion repeat, and restriction site polymorphism

A fine restriction map of the linear mitochondrial DNA of Tetrahymena pyriformis strain ST is presented. 1. Based on agarose gel electrophoresis data together with limited nucleotide sequences available on some restriction fragments, we estimate the actual size of this genome to be about 55,000 base pairs. 2. Seven tRNA gene locations have been assigned, which are scattered along the genome length. Six of these locations encode the genes for tRNA(phe), tRNA(his), tRNA(trp), and tRNA(glu), and the duplicate tRNA(tyr) genes which are located at the inverted terminal repeat segments. The tRNA gene(s) encoded in one location has not been identified. We have not yet found the tRNA(leu) and tRNA(met) genes, which were previously shown to be encoded in the genome (Chiu et al. 1974; Suyama 1982). 3. We have mapped the 14S rRNA gene by sequencing the 170 bp segment of EcoRI fragment 8 and by aligning its sequence with E. coli 16S rRNA. From our recent complete sequence data the gene size was found to be about 1,650 bp, which is unexpectedly large for the 14S rRNA which has an estimated size of 1,300 bp. The 14S rRNA is probably a cleavage product of the larger primary transcript of which 200-300 bases of the 5' end are missing. 4. The duplicate copies of the 21S rRNA gene at the terminal duplication inversion segments were analyzed. ClaI fragment 7 (1,500 bp) corresponds in sequence from base position 850 to 2,390 of the 20S rRNA gene of Paramecium mitochondrial DNA (Seilhamer et al. 1984b). The 21S gene is approximately 2,500 bp long. The presence of some restriction site polymorphism is apparent in this segment. 5. Each of the 21S gene copies precedes the tRNA(tyr) gene, but the space flanking one tRNA(tyr) gene differs in size and restriction sites from the space flanking another tRNA(tyr) gene. Thus, this space corresponds to the segment of an imperfect match in the terminal duplication inversion of Goldbach et al. (1978a). 6. Saccharomyces cerevisiae mitochondrial probes including Cob, ATPase VI and IX, and cytochrome oxidase I gene sequences, 21S and 15S rRNAs, and mouse mitochondrial DNA showed no significant hybridization with any restriction fragments of Tetrahymena mitochondrial DNA. The results are in accordance with an extensive sequence divergence previously found in the Tetrahymena mitochondrial genome (Goldbach et al. 1977).