DNA sequence and secondary structures of the large subunit rRNA coding regions and its two class I introns of mitochondrial DNA fromPodospora anserina

Summary DNA sequence analysis has shown that the gene coding for the mitochondrial (mt) large subunit ribosomal RNA (rRNA) fromPodospora anserina is interrupted by two class I introns. The coding region for the large subunit rRNA itself is 3715 bp and the two introns are 1544 (r1) and 2404 (r2) bp in length. Secondary structure models for the large subunit rRNA were constructed and compared with the equivalent structure fromEscherichia coli 23S rRNA. The two structures were remarkably similar despite an 800-base difference in length. The additional bases in theP. anserina rRNA appear to be mostly in unstructured regions in the 3 part of the RNA. Secondary structure models for the two introns show striking similarities with each other as well as with the intron models from the equivalent introns inSaccharomyces cerevisiae, Neurospora crassa, andAspergillus nidulans. The long open reading frames in each intron are different from each other, however, and the nucleotide sequence similarity diverges as it proceeds away from the core structure. Each intron is located within regions of the large subunit rRNA gene that are highly conserved in both sequence and structure. Computer analysis showed that the open reading frame for intron r1 contained a common maturase-like polypeptide. The open reading frames of intron r2 apeared to be chimeric, displaying high sequence similarity with the open reading frames in the r1 and ATPase 6 introns ofN. crassa.     

The nucleotide sequence of the entire ribosomal DNA operon and the structure of the large subunit rRNA of Giardia muris

Summary The total nucleotide sequence of the rDNA of Giardia muris, an intestinal protozoan parasite of rodents, has been determined. The repeat unit is 7668 basepairs (bp) in size and consists of a spacer of 3314 bp, a small-subunit rRNA (SSU-rRNA) gene of 1429, and a large-subunit rRNA (LSU-rRNA) gene of 2698 bp. The spacer contains long direct repeats and is heterogeneous in size. The LSU-rRNA of G. muris was compared to that of the human intestinal parasite Giardia duodenalis, to the bird parasite Giardia ardeae, and to that of Escherichia coli. The LSU-rRNA has a size comparable to the 23S rRNA of E. coli but shows structural features typical for eukaryotes. Some variable regions are typically small and account for the overall smaller size of this rRNA. The structure of the G. muris LSU-rRNA is similar to that of the other Giardia rRNA, but each rRNA has characteristic features residing in a number of variable regions.Offprint requests to: H. van Keulen     

DNA sequence,structure, and phylogenetic relationship of the mitochondrial small-subunit rRNA from the red alga Chondrus crispus (Gigartinales,Rhodophytes)

The entire nucleotide sequence containing the small-subunit ribosomal RNA gene (SSU rRNA) from the mitochondrial genome of Chondrus crispus was determined. To our knowledge, this is the first sequence of a mitochondrial 16S-like rRNA from a red alga. The length of this gene is 1,376 nucleotides. Its secondary structure was constructed and compared with other known secondary structures from eubacteria and from mitochondria of land plants, green and brown algae, and fungi. Phylogenetic trees were built upon SSU rRNA sequence alignment from mitochondria and eubacteria. The results show that rhodophytes and chromophytes provide additional links in the evolution of mitochondria between the green plant lineage and the nonplant lineages.Correspondence to: C. Boyen     

Nucleotide sequence of nine protein-coding genes and 22 tRNAs in the mitochondrial DNA of the sea starPisaster ochraceus

Summary We have cloned and sequenced over 9 kb of the mitochondrial genome from the sea starPisaster ochraceus. Within a continuous 8.0-kb fragment are located the genes for NADH dehydrogenase subunits 1, 2, 3, and 4L (ND1, ND2, ND3, and ND4L), cytochrome oxidase subunits I, II, and III (COI, COII, and COIII), and adenosine triphosphatase subunits 6 and 8 (ATPase 6 and ATPase 8). This large fragment also contains a cluster of 13 tRNA genes between ND1 and COI as well as the genes for isoleucine tRNA between ND1 and ND2, arginine tRNA between COI and ND4L, lysine tRNA between COII and ATPase 8, and the serine (UCN) tRNA between COIII and ND3. The genes for the other five tRNAs lie outside this fragment. The gene for phenylalanine tRNA is located between cytochrome b and the 12S ribosomal genes. The genes for tRNA^glu and tRNA^thr are 3 to the 12S ribosomal gene. The tRNAs for histidine and serine (AGN) are adjacent to each other and lie between ND4 and ND5. These data confirm the novel gene order in mitochondrial DNA (mtDNA) of sea stars and delineate additional distinctions between the sea star and other mtDNA molecules.     

Nuclear counterparts of the cytoplasmic mitochondrial 12S rRNA gene: A problem of ancient DNA and molecular phylogenies

Monkey mummy bones and teeth originating from the North Saqqara Baboon Galleries (Egypt), soft tissue from a mummified baboon in a museum collection, and nineteenth/twentieth-century skin fragments from mangabeys were used for DNA extraction and PCR amplification of part of the mitochondrial 12S rRNA gene. Sequences aligning with the 12S rRNA gene were recovered but were only distantly related to contemporary monkey mitochondrial 12S rRNA sequences. However, many of these sequences were identical or closely related to human nuclear DNA sequences resembling mitochondrial 12S rRNA (isolated from a cell line depleted in mitochondria) and therefore have to be considered contamination. Subsequently in a separate study we were able to recover genuine mitochondrial 12S rRNA sequences from many extant species of nonhuman Old World primates and sequences closely resembling the human nuclear integrations. Analysis of all sequences by the neighbor-joining (NJ) method indicated that mitochondrial DNA sequences and their nuclear counterparts can be divided into two distinct clusters. One cluster contained all temporary cytoplasmic mitochondrial DNA sequences and approximately half of the monkey nuclear mitochondriallike sequences. A second cluster contained most human nuclear sequences and the other half of monkey nuclear sequences with a separate branch leading to human and gorilla mitochondrial and nuclear sequences. Sequences recovered from ancient materials were equally divided between the two clusters. These results constitute a warning for when working with ancient DNA or performing phylogenetic analysis using mitochondrial DNA as a target sequence: Nuclear counterparts of mitochondrial genes may lead to faulty interpretation of results.Correspondence to: A.C. van der Kuyl     

Mitochondrial DNA coding region sequences support the phylogenetic distinction of two Indian wolf species

Two small endangered populations of Indian wolves were recently shown to be distant from other wolf and dog mtDNA lineages characterized so far. None of the inner branches in the tree of canid species based on partial hypervariable D-loop sequences were, however, statistically supported by the data raising the question whether the two Indian wolf lineages represent two new species, occupying an intermediate position between Canis latrans and C. lupus or have diverged from the sub-species of C. lupus due to isolation and drift. Here we report complete D-loop, cytochrome b, and 16S rRNA sequences data for 23 additional wolves from India analysed in the context of other canid species. Extended analyses of D-loop data and partial sequences of 16S rRNA showed highly reticulated pattern and were unable to resolve unambiguously the phylogenetic relationship of Indian wolves among other canid species. The phylogenetic reconstructions of cytochrome b sequences, however gave significant statistical support for the inner branches supporting genetic distinction of the two Indian wolf lineages within themselves as well as from all other wolves of the world, including individuals belonging to subspecies C. lupus chanco and C. lupus pallipes to which the two Indian wolf populations have been traditionally assigned. Their genetic differentiation relative to worldwide variation of wolves supports the suggestion to treat them as separate wolf species, C. himalayensis and C. indica .     

Phylogenetic estimation of Mytilidae in the East China Sea inferred from mitochondrial genes and nuclear DNA sequence variation

We performed a phylogenetic estimation of the family Mytilidae in the East China Sea based on nuclear internal transcribed spacer (ITS) genes and two mitochondrial genes (COI and 16S RNA). Analysis of five mytilid species based on each of the three genes resulted in mostly congruent trees, although there were some discrepancies in the classification of these species. We combine the results obtained from the three separate analyses to provide a phylogenetic estimation of Mytilidae. We found that the Mytilidae was divided into two major lineages: in one clade, Mytilus galloprovincialis was grouped with Mytilus coruscus; in the second clade, Septifer bilocularis was placed at the basal position in an individual clade, and Perna viridis and Musculista senhousia were recovered as a monophyletic group. Although these finding provide important insights into the taxonomic relationships among the Mytilidae, many aspects of Mytilidae phylogeny remain unresolved. Further analysis based on more molecular information and extensive taxon sampling is necessary to elucidate the phylogenetic relationships among the major lineages within the Mytilidae.     

The landscape of mitochondrial DNA variation in human colorectal cancer on the background of phylogenetic knowledge

Recently, an increasing number of studies indicate that mutations in mitochondrial genome may contribute to cancer development or metastasis. Hence, it is important to determine whether the mitochondrial DNA might be a good, clinically applicable marker of cancer. This review describes hereditary as well as somatic mutations reported in mitochondrial DNA of colorectal cancer cells. We showed here that the entire mitochondrial genome mutational spectra are different in colorectal cancer and non-tumor cells. We also placed the described mutations on the phylogenetic context, which highlighted the recurrent problem of data quality. Therefore, the most important rules for adequately assessing the quality of mitochondrial DNA sequence analysis in cancer have been summarized. As follows from this review, neither the reliable spectrum of mtDNA somatic mutations nor the association between hereditary mutations and colorectal cancer risk have been resolved. This indicates that only high resolution studies on mtDNA variability, followed by a proper data interpretation employing phylogenetic knowledge may finally verify the utility of mtDNA sequence (if any) in clinical practice.     

DNA sequence and secondary structure of the mitochondrial small subunit ribosomal RNA coding region including a group-IC2 intron from the cultivated basidiomycete Agrocybe aegerita

Due to their structural complexity and their evolutionary dimension, rRNAs are the most investigated nucleic acids in prokaryotes, eukaryotes and organelles. However, no complete sequence of a mitochondrial small subunit (SSU) rRNA was available in the basidiomycotina subdivision. The mitochondrial gene encoding the SSU rRNA of the cultivated basidiomycete Agrocybe aegerita was cloned and its complete nucleotide sequence achieved; the 5′- and 3′-ends were localized by nuclease S1 mapping, leading to a size of 3277 nt. The secondary structure of the SSU rRNA (1906 nt in size) possessed all the helices and loops of the prokaryotic model; a unique modification was found in a conserved nucleotide predicted by the model: the nt 487 was A instead of C. The same modification, has been found in all the partial basidiomycete mitochondrial sequences available in databases. The Agrocybe aegerita SSU rRNA was characterized by large and unusual extensions leading to additional helices in the variable domains V4, V6 and V9, which were the longest of the known prokaryotic or mitochondrial SSU rRNAs. Nucleotide sequence analysis indicated a 1371-bp intron, belonging to subgroup-IC2, located in a conserved loop in the 3′-part of the SSU rRNA. This intron, which is the second example reported in a fungal mitochondrial SSU rDNA, encoded a putative protein (407 aa) sharing homologies with endonucleases involved in group-I intron mobility. This report constitutes the first complete mitochondrial SSU rRNA sequence and secondary structure of any member of the basidiomycotina subdivision.     

Characterization of a novel small RNA encoded by Dictyostelium discoideum mitochondrial DNA

In this study, we analyzed a mitochondrial small (ms) RNA in Dictyostelium discoideum, which is 129 nucleotides long and has a GC content of only 22.5%. In the mitochondrial DNA, a single-copy gene (msr) for the ms RNA was located downstream of the gene for large-subunit rRNA. The location of msr was similar to that of the 5S rRNA gene in prokaryotes and chloroplasts, but clearly different from that in mitochondria of plants, liverwort and the chlorophycean alga Prototheca wikerhamii, in which small-subunit rRNA and 5S rRNA genes are closely linked. The primary sequence of ms RNA showed low homology with mitochondrial 5S rRNA from plants, liverwort and the chlorophycean alga, but the proposed secondary structure of ms RNA was similar to that of cytoplasmic 5S rRNA. In addition, ms RNA showed a highly conserved GAAC sequence in the same loop as in common 5S rRNA. However, ms RNA was detected mainly in the mitochondrial 25 000 × g supernatant fraction which was devoid of ribosomes. It is possible that ms RNA is an evolutionary derivative of mitochondrial 5S rRNA. Received: 17 May 1997 / Accepted: 26 August 1997     

The complete nucleotide sequence of the mitochondrial DNA of the fin whale,Balaenoptera physalus

Summary The composition of the mitochondrial DNA (mtDNA) of the fin whale,Balaenoptera physalus, was determined. The length of the molecule is 16,398 bp, and its organization conforms with that of other mammals. The general similarity between the mtDNA of the fin whale and the cow is greater than the similarity between the fin whale and other species (human, mouse, rat) in which the composition of the entire molecule has been described. The D-loop region of the mtDNA of the fin whale is 81% identical to the D-loop of dolphin DNA, and the central portion of the D-loop is similar to the bovine D-loop. The accumulation of transversions and gaps in the 12S and 16S rRNA genes was assessed by comparing the fin whale, cow, and human. The sequence difference between human and the whale and human and the cow was at the same level, indicating that the rate of evolution of the mtDNA rRNA genes is about the same in artiodactyls and cetaceans. In the 12S rRNA gene an accumulation rate of 0.05% per million years places the separation of cetaceans and artiodactyls at about 55 million years ago. The corresponding figure for human and either the whale or the cow is about 80 million years. In the 16S rRNA gene a 0.08% accumulation rate of transversions and gaps per million years yields concurring figures. A comparison between the cytochromeb gene of the fin whale and cytochromeb sequences in the literature, including dolphin (Stenella) sequences, identified the cetaceans as monophyletic and the artiodactyls as their closest relatives. The comparison between the cytochromeb sequences of the fin whale andStenella showed that differences in codon positions one or two were frequently associated with a change in another codon position.     

Polymorphic sequence in the D-loop region of equine mitochondrial DNA

The D-loop regions in equine mitochondrial DNA were cloned from three thoroughbred horses by polymerase chain reaction (PCR). The total number of bases in the D-loop region were 1114bp, 1115bp and 1146bp. The equine D-loop region is A/T rich like many other mammalian D-loops. The large central conserved sequence block and small conserved sequence blocks 1, 2 and 3, that are common to other mammals, were observed. Between conserved sequence blocks 1 and 2 there were tandem repeats of an 8bp equine-specific sequence TGTGCACC, and the number of tandem repeats differed among individual horses. The base composition in the unit of these repeats is G/C rich as are the short repeats in the D-loops of rabbit and pig. Comparing DNA sequences between horse and other mammals, the difference in the D-loop region length is mostly due to the difference in the number of DNA sequences at both extremities. The similarities of the DNA sequences are in the middle part of the D-loop. In comparison of the sequences among three thoroughbred horses, it was determined that the region between tRNA^Pro and the large central conserved sequence block was the richest in variation. PCR primers in the D-loop region were designed and the expected maternal inheritance was confirmed by PCR-RFLP (restriction fragment length polymorphism).     

Identification of the most informative regions of the mitochondrial genome for phylogenetic and coalescent analyses

Analysis of complete mitochondrial genome sequences is becoming increasingly common in genetic studies. The availability of full genome datasets enables an analysis of the information content distributed throughout the mitochondrial genome in order to optimize the research design of future evolutionary studies. The goal of our study was to identify informative regions of the human mitochondrial genome using two criteria: (1) accurate reconstruction of a phylogeny and (2) consistent estimates of time to most recent common ancestor (TMRCA). We created two series of datasets by deleting individual genes of varied length and by deleting 10 equal-size fragments throughout the coding region. Phylogenies were statistically compared to the full-coding-region tree, while coalescent methods were used to estimate the TMRCA and associated credible intervals. Individual fragments important for maintaining a phylogeny similar to the full-coding-region tree encompassed bp 577-2122 and 11,399-16,023, including all or part of 12S rRNA, 16S rRNA, ND4, ND5, ND6, and cytb. The control region only tree was the most poorly resolved with the majority of the tree manifest as an unresolved polytomy. Coalescent estimates of TMRCA were less sensitive to removal of any particular fragment(s) than reconstruction of a consistent phylogeny. Overall, we discovered that half the genome, i.e., bp 3669-11,398, could be removed with no significant change in the phylogeny (p(AU)=0.077) while still maintaining overlap of TMRCA 95% credible intervals. Thus, sequencing a contiguous fragment from bp 11,399 through the control region to bp 3668 would create a dataset that optimizes the information necessary for phylogenetic and coalescent analyses and also takes advantage of the wealth of data already available on the control region.     

Molecular phylogeny of parabasalids inferred from small subunit rRNA sequences, with emphasis on the Hypermastigea

Small subunit rRNA gene sequences were identified without cultivation from parabasalid symbionts of termites belonging to the hypermastigid orders Trichonymphida (the genera Hoplonympha, Staurojoenina, Teranympha, and Eucomonympha) and Spirotrichonymphida (Spirotrichonymphella), and from four yet-unidentified parabasalid symbionts of the termite Incisitermes minor. All these new sequences were analyzed by Bayesian, likelihood, and parsimony methods in a broad phylogeny including all identified parabasalid sequences available in databases and some as yet unidentified sequences probably derived from hypermastigids. A salient point of our study focused on hypermastigids was the polyphyly of this class. We also noted a clear dichotomy between Trichonymphida and the other parabasalid taxa. However, this hypermastigid order was apparently polyphyletic, probably reflecting its morphological diversity. Among Trichonymphida, Teranympha (Teranymphidae) grouped together with the members of the family Eucomonymphidae, suggesting that its family status is ambiguous. The monophyletic lineage composed by Spirotrichonymphida exhibited a narrower branching pattern than Trichonymphida. The root of parabasalids was examined but could not be discerned accurately.     

A novel mutation in the mitochondrial 16S rRNA gene in a patient with MELAS syndrome, diabetes mellitus, hyperthyroidism and cardiomyopathy

Using RNase protection analysis, we found a novel C to G mutation at nucleotide position 3093 of mitochondrial DNA (mtDNA) in a previously reported 35-year-old woman exhibiting clinical features of mitochondrial myopathy, encephalopathy, lactic acidosis and stroke-like episodes (MELAS) syndrome together with diabetes mellitus, hyperthyroidism and cardiomyopathy. The patient also had an A3243G mutation in the tRNA(Leu(UUR)) gene and a 260-base pair duplication in the D-loop of mtDNA. The fibroblasts of the patient were cultured and used for the construction of cybrids using cytoplasmic transfer of the patient's mtDNA to the mtDNA-less rho(0) cells. RNA isolated from the cybrids was subjected to RNase protection analysis, and a C3093G transversion at the 16S rRNA gene and a MELAS-associated A3243G mutation of mtDNA were detected. The novel C3093G mutation together with the A3243G transition were found in muscle biopsies, hair follicles and blood cells of this patient and also in her skin fibroblasts and cybrids. The proportion of the C3093G mutant mtDNA in muscle biopsies of the patient was 51%. In contrast, the mutation was not detected in three sons of the proband. To characterize the impact of the mtDNA mutation-associated defects on mitochondrial function, we determined the respiratory enzyme activities of the primary culture of fibroblasts established from the proband, her mother and her three sons. The proportions of mtDNA with the C3093G transversion and the A3243G transition in the fibroblasts of the proband were 45 and 58%, respectively. However, the fibroblasts of the proband's mother and children harbored lower levels of mtDNA with the A3243G mutation but did not contain the C3093G mutation. The complex I activity in the proband's fibroblasts was decreased to 47% of the control but those of the fibroblasts of the mother and three sons of the proband were not significantly changed. These findings suggest that the C3093G transversion together with the A3243G transition of mtDNA impaired the respiratory function of mitochondria and caused the atypical MELAS syndrome associated with diabetes mellitus, hyperthyroidism and cardiomyopathy in this patient.     

Molecular phylogeny of parabasalids based on small subunit rRNA sequences, with emphasis on the Trichomonadinae subfamily

We determined small subunit ribosomal DNA sequences from three parabasalid species, Trichomitus batrachorum strain R105, Tetratrichomonas gallinarum, and Pentatrichomonas hominis belonging to the Trichomonadinae subfamily. Unrooted molecular phylogenetic trees inferred by distance, parsimony, and likelihood methods reveal four discrete clades among the parabasalids. The Trichomonadinae form a robust monophyletic group. Within this subfamily T. gallinarum is closely related to Trichomonas species as supported by morphological data, with P. hominis and Pseudotrypanosoma giganteum occupying basal positions. Our analysis does not place T. batrachorum within the Trichomonadinae. Trichomitus batrachorum (strains R105 and BUB) and Hypotrichomonas acosta form a well-separated cluster, suggesting the genus Trichomitus is polyphyletic. The emergence of T. batrachorum precedes the Trichomonadinae-Tritrichomonadinae dichotomy, emphasizing its pivotal evolutionary position among the Trichomonadidae. A third cluster unites the Devescovinidae and the Calonymphidae. The fourth clade contains the three hypermastigid sequences from the genus Trichonympha, which exhibit the earliest emergence among the parabasalids. The addition of these three new parabasalid species did not however resolve ambiguities regarding the relative branching order of the parabasalid clades. The phylogenetic positions of Tritrichomonas faetus, Monocercomonas sp., Dientamoeba fragilis, and the unidentified Reticulitermes flavipes gut symbiont 1 remain unclear.     

Phylogenetic position of the trichomonad parasite of turkeys, Histomonas meleagridis (Smith) Tyzzer, inferred from small subunit rRNA sequence.

The phylogenetic position of the trichomonad, Histomonas meleagridis was determined by analysis of small subunit rRNAs. Molecular trees including all identified parabasalid sequences available in data bases were inferred by distance, parsimony, and likelihood methods. All reveal a close relationship between H. meleagridis, and Dientamoeba fragilis. Moreover, small subunit rRNAs of both amoeboid species have a reduced G + C content and increased chain length relative to other parabasalids. Finally, the rRNA genes from H. meleagridis and D. fragilis share a recent common ancestor with Tritrichomonasfoetus, which exhibits a more developed cytoskeleton. This indicates that Histomonas and Dientamoeba secondarily lost most of the typical trichomonad cytoskeletal structures and hence, do not represent primitive morphologies. A global phylogeny of parabasalids revealed significant discrepancies with morphology-based classifications, such as the polyphyly of most of the parabasalid families and classes included in our study.