Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species

The mitochondrial genome (mtDNA) of Metazoa is a good model system for evolutionary genomic studies and the availability of more than 1000 sequences provides an almost unique opportunity to decode the mechanisms of genome evolution over a large phylogenetic range. In this paper, we review several structural features of the metazoan mtDNA, such as gene content, genome size, genome architecture and the new parameter of gene strand asymmetry in a phylogenetic framework. The data reviewed here show that: (1) the plasticity of Metazoa mtDNA is higher than previously thought and mainly due to variation in number and location of tRNA genes; (2) an exceptional trend towards stabilization of genomic features occurred in deuterostomes and was exacerbated in vertebrates, where gene content, genome architecture and gene strand asymmetry are almost invariant. Only tunicates exhibit a very high degree of genome variability comparable to that found outside deuterostomes. In order to analyse the genomic evolutionary process at short evolutionary distances, we have also compared mtDNAs of species belonging to the same genus: the variability observed in congeneric species significantly recapitulates the evolutionary dynamics observed at higher taxonomic ranks, especially for taxa showing high levels of genome plasticity and/or fast nucleotide substitution rates. Thus, the analysis of congeneric species promises to be a valuable approach for the assessment of the mtDNA evolutionary trend in poorly or not yet sampled metazoan groups.     

Characterization of a mitochondrial ORF from the gender-associated mtDNAs of Mytilus spp. (Bivalvia: Mytilidae): identification of the "missing" ATPase 8 gene

Bivalve species are characterized by extraordinary variability in terms of mitochondrial (mt) genome size, gene arrangement and tRNA gene number. Many species are thought to lack the mitochondrial protein-coding gene atp8. Of these species, the Mytilidae appears to be the only known taxon with doubly uniparental inheritance of mtDNA that does not possess the atp8 gene. This raises the question as to whether mytilids have completely lost the ATP8 protein, whether the gene has been transferred to the nucleus or whether they possess a highly modified version of the gene/protein that has led to its lack of annotation. In the present study, we re-investigated all complete (or nearly complete) F and M mytilid mt genomes previously sequenced for the presence of conserved open reading frames (ORFs) that might code for ATP8 and/or have other functional importance in these bivalves. We also revised the annotations of all available complete mitochondrial genomes of bivalves and nematodes that are thought to lack atp8 in an attempt to detect it. Our results indicate that a novel mytilid ORF of significant length (i.e., the ORF is >85 amino acids in length), with complete start and stop codons, is a candidate for the atp8 gene: (1) it possesses a pattern of evolution expected for a protein-coding gene evolving under purifying selection (i.e., the 3rd>1st>2nd codon pattern of evolution), (2) it is actively transcribed in Mytilus species, (3) it has one predicted transmembrane helix (as do other metazoan ATP8 proteins), (4) it has conserved functional motifs and (5), comparisons of its amino acid sequence with ATP8 sequences of other molluscan or bivalve species reveal similar hydropathy profiles. Furthermore, our revised annotations also confirmed the mt presence of atp8 in almost all bivalve species and in one nematode species. Our results thus support recognizing the presence of ATPase 8 in most bivalves mt genomes (if not all) rather than the continued characterization of these genomes as lacking this gene.     

Mitochondrial genomes and divergence times of crocodile newts: inter-islands distribution of Echinotriton andersoni and the origin of a unique repetitive sequence found in Tylototriton mt genomes

Crocodile newts, which constitute the genera Echinotriton and Tylototriton, are known as living fossils, and these genera comprise many endangered species. To identify mitochondrial (mt) genes suitable for future population genetic analyses for endangered taxa, we determined the complete nucleotide sequences of the mt genomes of the Japanese crocodile newt Echinotriton andersoni and Himalayan crocodile newt Tylototriton verrucosus. Although the control region (CR) is known as the most variable mtDNA region in many animal taxa, the CRs of crocodile newts are highly conservative. Rather, the genes of NADH dehydrogenase subunits and ATPase subunit 6 were found to have high sequence divergences and to be usable for population genetics studies. To estimate the inter-population divergence ages of E. andersoni endemic to the Ryukyu Islands, we performed molecular dating analysis using whole and partial mt genomic data. The estimated divergence ages of the inter-island individuals are older than the paleogeographic segmentation ages of the islands, suggesting that the lineage splits of E. andersoni populations were not caused by vicariant events. Our phylogenetic analysis with partial mt sequence data also suggests the existence of at least two more undescribed species in the genus Tylototriton. We also found unusual repeat sequences containing the 3' region of cytochrome apoenzyme b gene, whole tRNA-Thr gene, and a noncoding region (the T-P noncoding region characteristic in caudate mtDNAs) from T. verrucosus mtDNA. Similar repeat sequences were found in two other Tylototriton species. The Tylototriton taxa with the repeats become a monophyletic group, indicating a single origin of the repeat sequences. The intra-and inter-specific comparisons of the repeat sequences suggest the occurrences of homologous recombination-based concerted evolution among the repeat sequences.     

Analysis of complete mitochondrial genome sequences increases phylogenetic resolution of bears (Ursidae), a mammalian family that experienced rapid speciation

Background  

Despite the small number of ursid species, bear phylogeny has long been a focus of study due to their conservation value, as all bear genera have been classified as endangered at either the species or subspecies level. The Ursidae family represents a typical example of rapid evolutionary radiation. Previous analyses with a single mitochondrial (mt) gene or a small number of mt genes either provide weak support or a large unresolved polytomy for ursids. We revisit the contentious relationships within Ursidae by analyzing complete mt genome sequences and evaluating the performance of both entire mt genomes and constituent mtDNA genes in recovering a phylogeny of extremely recent speciation events.     

A comprehensive and validated molecular taxonomy of beaked whales, family Ziphiidae

DNA sequences from orthologous loci can provide universal characters for taxonomic identification. Molecular taxonomy is of particular value for groups in which distinctive morphological features are difficult to observe or compare. To assist in species identification for the little known family Ziphiidae (beaked whales), we compiled a reference database of mitochondrial DNA (mtDNA) control region (437 bp) and cytochrome b (384 bp) sequences for all 21 described species in this group. This mtDNA database is complemented by a nuclear database of actin intron sequences (925 bp) for 17 of the 21 species. All reference sequences were derived from specimens validated by diagnostic skeletal material or other documentation, and included four holotypes. Phylogenetic analyses of mtDNA sequences confirmed the genetic distinctiveness of all beaked whale species currently recognized. Both mitochondrial loci were well suited for species identification, with reference sequences for all known ziphiids forming robust species-specific clades in phylogenetic reconstructions. The majority of species were also distinguished by nuclear alleles. Phylogenetic comparison of sequence data from "test" specimens to these reference databases resulted in three major taxonomic discoveries involving animals previously misclassified from morphology. Based on our experience with this family and the order Cetacea as a whole, we suggest that a molecular taxonomy should consider the following components: comprehensiveness, validation, locus sensitivity, genetic distinctiveness and exclusivity, concordance, and universal accessibility and curation.     

The complete mitochondrial genome of Fusarium oxysporum: insights into fungal mitochondrial evolution

The complete mitochondrial DNA (mtDNA) sequence was determined for the phytopathogenic fungus Fusarium oxysporum. It is 34,477 bp long, maps circularly, and encodes for 14 protein-coding, 25 tRNA and 2 rRNA genes. The nucleotide and amino acid data sets from its 14 concatenated protein-coding mitochondrial (mt) genes were used along with gene order comparisons for an extensive phylogenetic study of the Subphylum Pezizomycotina. Our results are in agreement with current taxonomic treatments and additionally provide better statistical support for all relationships within Pezizomycotina when compared to analyses based on single or few gene data sets. The gene order of F. oxysporum was consistent with that established in the order Hypocreales (Class: Sordariomycetes) and enhanced previous suppositions on the ancestral state of Sordariomycetes. In comparison with mt genomes of the other orders it added further insights to the evolution of Pezizomycotina.     

Oesophagostomum dentatum and Oesophagostomum quadrispinulatum: characterization of the complete mitochondrial genome sequences of the two pig nodule worms

In the present study, the complete mitochondrial DNA (mtDNA) sequences of the pig nodule worm Oesophagostomum quadrispinulatum were determined for the first time, and the mt genome of Oesophagostomum dentatum from China was also sequenced for comparative analysis of their gene contents and genome organizations. The mtDNA sequences of O. dentatum China isolate and O. quadrispinulatum were 13,752 and 13,681 bp in size, respectively. Each of the two mt genomes comprises 36 genes, including 12 protein-coding genes, two ribosomal RNA and 22 transfer RNA genes, but lacks the ATP synthetase subunit 8 gene. All genes are transcribed in the same direction and have a nucleotide composition high in A and T. The contents of A+T are 75.79% and 77.52% for the mt genomes of O. dentatum and O. quadrispinulatum, respectively. Phylogenetic analyses using concatenated amino acid sequences of the 12 protein-coding genes, with three different computational algorithms (maximum likelihood, maximum parsimony and Bayesian inference), all revealed that O. dentatum and O. quadrispinulatum represent distinct but closely-related species. These data provide novel and useful markers for studying the systematics, population genetics and molecular diagnosis of the two pig nodule worms.     

Characterization of the complete mitochondrial genomes of five Eimeria species from domestic chickens

Chicken coccidiosis caused by members of the genus Eimeria causes significant economic losses worldwide. In the present study we sequenced the complete mitochondrial DNA (mtDNA) sequences of six Eimeria species and analyzed features of their gene contents and genome organizations. The complete mt genomes of E. acervulina, E. brunetti, E. maxima, E. necatrix, E. tenella and E. praecox were 6179bp, 6148bp, 6169bp, 6214bp, 6213bp and 6174bp in size, respectively. All of the mt genomes consist of 3 genes for proteins (cox1, cox3, and cytb), 12 gene fragments for the large subunit (LSU) rRNA, and 7 gene fragments for the small subunit (SSU) rRNA, but no transfer RNA genes. The organization of the mt genomes is similar to that of Plasmodium, but distinct from Babesia and Theileria. The putative direction of translation for 3 genes (cox1, cox3, and cytb) was the same in all six Eimeria species. The contents of A+T of the mt genomes were 65.35% for E. acervulina, 65.43% for E. brunetti, 64.53% for E. maxima, 65.04% for E. necatrix, 64.98% for E. tenella and 65.59% for E. praecox. The AT bias has a significant effect on both the codon usage pattern and amino acid composition of proteins. Phylogenetic analyses using concatenated nucleotide sequences of the 2 protein-coding genes (cytb and cox1), with three different computational algorithms (Bayesian analysis, maximum parsimony and maximum likelihood), all revealed distinct groups with high statistical support, indicating that the six Eimeria spp. represent six distinct but closely-related species. These data provide novel mtDNA markers for studying the molecular epidemiology and population genetics of the six Eimeria spp., and should have implications for the molecular diagnosis, prevention and control of coccidiosis in domestic chickens.     

The Complete Mitochondrial Genome of Galba pervia (Gastropoda: Mollusca), an Intermediate Host Snail of Fasciola spp

Complete mitochondrial (mt) genomes and the gene rearrangements are increasingly used as molecular markers for investigating phylogenetic relationships. Contributing to the complete mt genomes of Gastropoda, especially Pulmonata, we determined the mt genome of the freshwater snail Galba pervia, which is an important intermediate host for Fasciola spp. in China. The complete mt genome of G. pervia is 13,768 bp in length. Its genome is circular, and consists of 37 genes, including 13 genes for proteins, 2 genes for rRNA, 22 genes for tRNA. The mt gene order of G. pervia showed novel arrangement (tRNA-His, tRNA-Gly and tRNA-Tyr change positions and directions) when compared with mt genomes of Pulmonata species sequenced to date, indicating divergence among different species within the Pulmonata. A total of 3655 amino acids were deduced to encode 13 protein genes. The most frequently used amino acid is Leu (15.05%), followed by Phe (11.24%), Ser (10.76%) and IIe (8.346%). Phylogenetic analyses using the concatenated amino acid sequences of the 13 protein-coding genes, with three different computational algorithms (maximum parsimony, maximum likelihood and Bayesian analysis), all revealed that the families Lymnaeidae and Planorbidae are closely related two snail families, consistent with previous classifications based on morphological and molecular studies. The complete mt genome sequence of G. pervia showed a novel gene arrangement and it represents the first sequenced high quality mt genome of the family Lymnaeidae. These novel mtDNA data provide additional genetic markers for studying the epidemiology, population genetics and phylogeographics of freshwater snails, as well as for understanding interplay between the intermediate snail hosts and the intra-mollusca stages of Fasciola spp..     

Characterization of the complete mitochondrial genome sequence of Spirometra erinaceieuropaei (Cestoda: Diphyllobothriidae) from China

Sparganosis, caused by the plerocercoid larvae of members of the genus Spirometra, can cause significant public health problem and considerable economic losses. In the present study, the complete mitochondrial DNA (mtDNA) sequence of Spirometra erinaceieuropaei from China was determined, characterized and compared with that of S. erinaceieuropaei from Japan. The gene arrangement in the mt genome sequences of S. erinaceieuropaei from China and Japan is identical. The identity of the mt genomes was 99.1% between S. erinaceieuropaei from China and Japan, and the complete mtDNA sequence of S. erinaceieuropaei from China is slightly shorter (2 bp) than that from Japan. Phylogenetic analysis of S. erinaceieuropaei with other representative cestodes using two different computational algorithms [Bayesian inference (BI) and maximum likelihood (ML)] based on concatenated amino acid sequences of 12 protein-coding genes, revealed that S. erinaceieuropaei is closely related to Diphyllobothrium spp., supporting classification based on morphological features. The present study determined the complete mtDNA sequences of S. erinaceieuropaei from China that provides novel genetic markers for studying the population genetics and molecular epidemiology of S. erinaceieuropaei in humans and animals.     

The complete mitochondrial genome of the American lobster, Homarus americanus (Crustacea, Decapoda)

Although relatively a large number of the complete mitochondrial genome sequences have been determined from various decapod species (29 mtDNA sequences reported so far), the information for the infraorder Astacidea (including lobsters, crayfishes, and their relatives) is very limited and represented by only one complete sequence from the Australian freshwater crayfish species Cherax destructor. In this study, we determined the complete mitochondrial DNA sequence of Homarus americanus, the first representative of the family Nephropidae to be fully characterized. Comparison of the gene arrangement reveals that H. americanus mtDNA is identical to those of other pancrustacean species but differs from the other astacidean species (C. destructor). Based on these data, it can be assumed that an idiosyncratic gene order discovered in C. destructor mtDNA may be secondarily acquired from the ancestral lineage of the Astacidea.     

A gene coding for an RPS2 protein is present in the mitochondrial genome of several cereals, but not in dicotyledons

A gene coding for a protein that shows homologies to prokaryotic ribosomal protein S2 is present in the mitochondrial (mt) genome of wheat (Triticum aestivum). The wheat gene is transcribed as a single mRNA which is edited by C-to-U conversions at seven positions, all resulting in alteration of the encoded amino acid. Homologous gene sequences are also present in the mt genomes of rice and maize, but we failed to identify the corresponding sequences in the mtDNA of all dicotyledonous species tested; in these species the mitochondrial RPS2 is probably encoded in the nucleus. The protein sequence deduced from the wheat rps2 gene sequence has a long C-terminal extension when compared to other prokaryotic RPS2 sequences. This extension presents no similarity with any known sequence and is not conserved in the maize or rice mitochondrial rps2 gene. Most probably, after translation, this peptide extension is processed by a specific peptidase to give rise to the mature wheat mitochondrial RPS2.     

The 203 kbp Mitochondrial Genome of the Phytopathogenic Fungus Sclerotinia borealis Reveals Multiple Invasions of Introns and Genomic Duplications

Here we report the complete sequence of the mitochondrial (mt) genome of the necrotrophic phytopathogenic fungus Sclerotinia borealis, a member of the order Helotiales of Ascomycetes. The 203,051 bp long mtDNA of S. borealis represents one of the largest sequenced fungal mt genomes. The large size is mostly determined by the presence of mobile genetic elements, which include 61 introns. Introns contain a total of 125,394 bp, are scattered throughout the genome, and are found in 12 protein-coding genes and in the ribosomal RNA genes. Most introns contain complete or truncated ORFs that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Integrations of mobile elements are also evidenced by the presence of two regions similar to fragments of inverton-like plasmids. Although duplications of some short genome regions, resulting in the appearance of truncated extra copies of genes, did occur, we found no evidences of extensive accumulation of repeat sequences accounting for mitochondrial genome size expansion in some other fungi. Comparisons of mtDNA of S. borealis with other members of the order Helotiales reveal considerable gene order conservation and a dynamic pattern of intron acquisition and loss during evolution. Our data are consistent with the hypothesis that horizontal DNA transfer has played a significant role in the evolution and size expansion of the S. borealis mt genome.     

The Cephalopod Loligo bleekeri Mitochondrial Genome: Multiplied Noncoding Regions and Transposition of tRNA Genes

We previously reported the sequence of a 9260-bp fragment of mitochondrial (mt) DNA of the cephalopod Loligo bleekeri [J. Sasuga et al. (1999) J. Mol. Evol. 48:692–702]. To clarify further the characteristics of Loligo mtDNA, we have sequenced an 8148-bp fragment to reveal the complete mt genome sequence. Loligo mtDNA is 17,211 bp long and possesses a standard set of metazoan mt genes. Its gene arrangement is not identical to any other metazoan mt gene arrangement reported so far. Three of the 19 noncoding regions longer than 10 bp are 515, 507, and 509 bp long, and their sequences are nearly identical, suggesting that multiplication of these noncoding regions occurred in an ancestral Loligo mt genome. Comparison of the gene arrangements of Loligo, Katharina tunicata, and Littorina saxatilis mt genomes revealed that 17 tRNA genes of the Loligo mt genome are adjacent to noncoding regions. A majority (15 tRNA genes) of their counterparts is found in two tRNA gene clusters of the Katharina mt genome. Therefore, the Loligo mt genome (17 tRNA genes) may have spread over the genome, and this may have been coupled with the multiplication of the noncoding regions. Maximum likelihood analysis of mt protein genes supports the clade Mollusca + Annelida + Brachiopoda but fails to infer the relationships among Katharina, Loligo, and three gastropod species. Received: 9 May 2001 / Accepted: 3 October 2001     

The complete mitochondrial genome of the Japanese mud shrimp Upogebia major (Crustacea, Decapoda)

We determined a full-length sequence of mitochondrial (mt) genome from Upogebia major. This is the first complete mt genome report for infraorder Thalassinidea in Decapoda, Crustacea. Our result showed that U. major generally followed a typical pancrustacean gene order but some tRNA genes showed a very unique gene arrangement such as duplication or translocation. Since none of the complete mt genome sequences in the infraorder Thalassinidea are available yet, this report will provide additional information in relation to mt genome diversity and evolution of the decapods.     

Mitochondrial phylogeny of hedgehogs and monophyly of Eulipotyphla

We sequenced the complete mitochondrial (mt) genomes of three insectivores: the long-eared hedgehog Hemiechinus auritus, the Japanese mole Mogera wogura, and the greater Japanese shrew-mole Urotrichus talpoides. These mtDNA data together with other previously sequenced mtDNAs were analyzed using a maximum likelihood method to infer their phylogenetic relationships among eutherians. Previous mitochondrial protein analyses used a simple model that did not consider site-heterogeneity, and Erinaceoidea (hedgehogs and moonrats) was placed at the basal eutherian position that is separated from Soricoidea (shrews) and Talpoidea (moles), suggesting the exclusion of the Erinaceoidea-Eulipotyphla tree. By including the new mtDNA sequences and introducing site-heterogeneity into the model, the Erinaceoidea-Eulipotyphla tree emerges as the best tree or as a tree with a log-likelihood score indistinguishable from that of the best tree. However, this conclusion depends on species sampling in Erinaceoidea, demonstrating the importance of both species sampling and use of an appropriate substitution model when inferring phylogenetic relationships.     

Complete mitochondrial genomes and novel gene rearrangements in two dicroglossid frogs, Hoplobatrachus tigerinus and Euphlyctis hexadactylus, from Bangladesh

We determined the complete nucleotide sequences of mitochondrial (mt) genomes from two dicroglossid frogs, Hoplobatrachus tigerinus (Indian Bullfrog) and Euphlyctis hexadactylus (Indian Green frog). The genome sizes are 20462 bp in H. tigerinus and 20280 bp in E. hexadactylus. Although both genomes encode the typical 37 mt genes, the following unique features are observed: 1) the ND5 genes are duplicated in H. tigerinus that have completely identical sequences, whereas duplicated ND5 genes in E. hexadactylus possessed dissimilar substitutions; 2) duplicated control region (CR) in H. tigerinus has almost identical sequences whereas single control region (CR) was found in E. hexadactylus; 3) the tRNA-Leu (CUN) gene is translocated from the LTPF tRNA cluster to downstream of ND5-1 in H. tigerinus, and the tRNA-Pro gene is translocated from the LTPF tRNA cluster to downstream of CR in E. hexadactylus; 4) pseudo tRNA-Leu (CUN) and tRNA-Pro genes are observed in E. hexadactylus; and 5) two tRNA-Met genes are encoded in both species, as observed in the previously reported dicroglossid mt genomes. Almost all observed gene rearrangements in H. tigerinus and E. hexadactylus can be explained by the tandem duplication and random loss model, except translocation of tRNA-Pro in E. hexadactylus. The novel mt genomic features found in this study may be useful for future phylogenetic studies in the dicroglossid taxa. However, the mt genome with interesting features found in the present study reveal a high level of variation of gene order and gene content, inspiring more research to understand the mechanisms behind gene and genome evolution in the dicroglossid and as well as in the amphibian taxa in future studies.     

A gene coding for an RPS2 protein is present in the mitochondrial genome of several cereals, but not in dicotyledons

A gene coding for a protein that shows homologies to prokaryotic ribosomal protein S2 is present in the mitochondrial (mt) genome of wheat (Triticum aestivum). The wheat gene is transcribed as a single mRNA which is edited by C-to-U conversions at seven positions, all resulting in alteration of the encoded amino acid. Homologous gene sequences are also present in the mt genomes of rice and maize, but we failed to identify the corresponding sequences in the mtDNA of all dicotyledonous species tested; in these species the mitochondrial RPS2 is probably encoded in the nucleus. The protein sequence deduced from the wheat rps2 gene sequence has a long C-terminal extension when compared to other prokaryotic RPS2 sequences. This extension presents no similarity with any known sequence and is not conserved in the maize or rice mitochondrial rps2 gene. Most probably, after translation, this peptide extension is processed by a specific peptidase to give rise to the mature wheat mitochondrial RPS2. Received: 20 November 1997 / Accepted: 29 January 1998