Evolution of the mitochondrial genome: protist connections to animals, fungi and plants

The past decade has seen the determination of complete mitochondrial genome sequences from a taxonomically diverse set of organisms. These data have allowed an unprecedented understanding of the evolution of the mitochondrial genome in terms of gene content and order, as well as genome size and structure. In addition, phylogenetic reconstructions based on mitochondrial DNA (mtDNA)-encoded protein sequences have firmly established the identities of protistan relatives of the animal, fungal and plant lineages. Analysis of the mtDNAs of these protists has provided insight into the structure of the mitochondrial genome at the origin of these three, mainly multicellular, eukaryotic groups. Further research into mtDNAs of taxa ancestral and intermediate to currently characterized organisms will help to refine pathways and modes of mtDNA evolution, as well as provide valuable phylogenetic characters to assist in unraveling the deep branching order of all eukaryotes.     

The complete mitochondrial DNA sequences of Nephroselmis olivacea and Pedinomonas minor. Two radically different evolutionary patterns within green algae.

Green plants appear to comprise two sister lineages, Chlorophyta (classes Chlorophyceae, Ulvophyceae, Trebouxiophyceae, and Prasinophyceae) and Streptophyta (Charophyceae and Embryophyta, or land plants). To gain insight into the nature of the ancestral green plant mitochondrial genome, we have sequenced the mitochondrial DNAs (mtDNAs) of Nephroselmis olivacea and Pedinomonas minor. These two green algae are presumptive members of the Prasinophyceae. This class is thought to include descendants of the earliest diverging green algae. We find that Nephroselmis and Pedinomonas mtDNAs differ markedly in size, gene content, and gene organization. Of the green algal mtDNAs sequenced so far, that of Nephroselmis (45,223 bp) is the most ancestral (minimally diverged) and occupies the phylogenetically most basal position within the Chlorophyta. Its repertoire of 69 genes closely resembles that in the mtDNA of Prototheca wickerhamii, a later diverging trebouxiophycean green alga. Three of the Nephroselmis genes (nad10, rpl14, and rnpB) have not been identified in previously sequenced mtDNAs of green algae and land plants. In contrast, the 25,137-bp Pedinomonas mtDNA contains only 22 genes and retains few recognizably ancestral features. In several respects, including gene content and rate of sequence divergence, Pedinomonas mtDNA resembles the reduced mtDNAs of chlamydomonad algae, with which it is robustly affiliated in phylogenetic analyses. Our results confirm the existence of two radically different patterns of mitochondrial genome evolution within the green algae.     

The complete mitochondrial DNA sequence of Mesostigma viride identifies this green alga as the earliest green plant divergence and predicts a highly compact mitochondrial genome in the ancestor of all green plants.

To gain insights into the nature of the mitochondrial genome in the common ancestor of all green plants, we have completely sequenced the mitochondrial DNA (mtDNA) of Mesostigma viride. This green alga belongs to a morphologically heterogeneous class (Prasinophyceae) that includes descendants of the earliest diverging green plants. Recent phylogenetic analyses of ribosomal RNAs (rRNAs) and concatenated proteins encoded by the chloroplast genome identified Mesostigma as a basal branch relative to the Streptophyta and the Chlorophyta, the two phyla that were previously thought to contain all extant green plants. The circular mitochondrial genome of Mesostigma resembles the mtDNAs of green algae occupying a basal position within the Chlorophyta in displaying a small size (42,424 bp) and a high gene density (86.6% coding sequences). It contains 65 genes that are conserved in other mtDNAs. Although none of these genes represents a novel coding sequence among green plant mtDNAs, four of them (rps1, sdh3, sdh4, and trnL[caa]) have not been reported previously in chlorophyte mtDNAs, and two others (rpl14 and trnI[gau]) have not been identified in the streptophyte mtDNAs examined so far (land-plant mtDNAs). Phylogenetic analyses of 19 concatenated mtDNA-encoded proteins favor the hypothesis that Mesostigma represents the earliest branch of green plant evolution. Four group I introns (two in rnl and two in cox1) and three group II introns (two in nad3 and one in cox2), two of which are trans-spliced at the RNA level, reside in Mesostigma mtDNA. The insertion sites of the three group II introns are unique to this mtDNA, suggesting that trans-splicing arose independently in the Mesostigma lineage and in the Streptophyta. The few structural features that can be regarded as ancestral in Mesostigma mtDNA predict that the common ancestor of all green plants had a compact mtDNA containing a minimum of 75 genes and perhaps two group I introns. Considering that the mitochondrial genome is much larger in size in land plants than in Mesostigma, we infer that mtDNA size began to increase dramatically in the Streptophyta either during the evolution of charophyte green algae or during the transition from charophytes to land plants.     

The complete mitochondrial DNA sequence of the green alga Pseudendoclonium akinetum (Ulvophyceae) highlights distinctive evolutionary trends in the chlorophyta and suggests a sister-group relationship between the Ulvophyceae and Chlorophyceae

The mitochondrial genome has undergone radical changes in both the Chlorophyta and Streptophyta, yet little is known about the dynamics of mtDNA evolution in either of these lineages. In the Chlorophyta, which comprises four of the five recognized classes of green algae (Prasinophyceae, Trebouxiophyceae, Ulvophyceae, and Chlorophyceae), the mitochondrial genome varies from 16 to 55 kb. This genome has retained a compact gene organization and a relatively complex gene repertoire ("ancestral" pattern) in the basal lineages represented by the Trebouxiophyceae and Prasinophyceae, whereas it has been reduced in size and gene complement and tends to evolve much more rapidly at the sequence level ("reduced-derived" pattern of evolution) in the Chlorophyceae and the lineage leading to the enigmatic chlorophyte Pedinomonas. To gain information about the evolutionary trends of mtDNA in the Ulvophyceae and also to gain insights into the phylogenetic relationships between ulvophytes and other chlorophytes, we have determined the mtDNA sequence of Pseudendoclonium akinetum. At 95,880 bp, Pseudendoclonium mtDNA is the largest green-algal mitochondrial genome sequenced to date and has the lowest gene density. These derived features are reminiscent of the "expanded" pattern exhibited by embryophyte mtDNAs, indicating that convergent evolution towards genome expansion has occurred independently in the Chlorophyta and Streptophyta. With 57 conserved genes, the gene repertoire of Pseudendoclonium mtDNA is slightly smaller than those of the prasinophyte Nephroselmis olivacea and the trebouxiophyte Prototheca wickerhamii. This ulvophyte mtDNA contains seven group I introns, four of which have homologs in green-algal mtDNAs displaying an "ancestral" or a "reduced-derived" pattern of evolution. Like its counterpart in the chlorophycean green alga Scenedesmus obliquus, it features numerous small, dispersed repeats in intergenic regions and introns. Its overall rate of sequence evolution appears to be accelerated to an intermediary level as compared with the rates observed in "ancestral" and "reduced-derived" mtDNAs. In agreement with the finding that Pseudendoclonium mtDNA exhibits features typical of both the "ancestral" and "reduced-derived" patterns of evolution, phylogenetic analyses of seven mtDNA-encoded proteins revealed a sister-group relationship between this ulvophyte and chlorophytes displaying "reduced-derived" mtDNAs.     

Seventeen new complete mtDNA sequences reveal extensive mitochondrial genome evolution within the Demospongiae

Two major transitions in animal evolution--the origins of multicellularity and bilaterality--correlate with major changes in mitochondrial DNA (mtDNA) organization. Demosponges, the largest class in the phylum Porifera, underwent only the first of these transitions and their mitochondrial genomes display a peculiar combination of ancestral and animal-specific features. To get an insight into the evolution of mitochondrial genomes within the Demospongiae, we determined 17 new mtDNA sequences from this group and analyzing them with five previously published sequences. Our analysis revealed that all demosponge mtDNAs are 16- to 25-kbp circular molecules, containing 13-15 protein genes, 2 rRNA genes, and 2-27 tRNA genes. All but four pairs of sampled genomes had unique gene orders, with the number of shared gene boundaries ranging from 1 to 41. Although most demosponge species displayed low rates of mitochondrial sequence evolution, a significant acceleration in evolutionary rates occurred in the G1 group (orders Dendroceratida, Dictyoceratida, and Verticillitida). Large variation in mtDNA organization was also observed within the G0 group (order Homosclerophorida) including gene rearrangements, loss of tRNA genes, and the presence of two introns in Plakortis angulospiculatus. While introns are rare in modern-day demosponge mtDNA, we inferred that at least one intron was present in cox1 of the common ancestor of all demosponges. Our study uncovered an extensive mitochondrial genomic diversity within the Demospongiae. Although all sampled mitochondrial genomes retained some ancestral features, including a minimally modified genetic code, conserved structures of tRNA genes, and presence of multiple non-coding regions, they vary considerably in their size, gene content, gene order, and the rates of sequence evolution. Some of the changes in demosponge mtDNA, such as the loss of tRNA genes and the appearance of hairpin-containing repetitive elements, occurred in parallel in several lineages and suggest general trends in demosponge mtDNA evolution.     

The mitochondrial genome of Chara vulgaris: insights into the mitochondrial DNA architecture of the last common ancestor of green algae and land plants

Mitochondrial DNA (mtDNA) has undergone radical changes during the evolution of green plants, yet little is known about the dynamics of mtDNA evolution in this phylum. Land plant mtDNAs differ from the few green algal mtDNAs that have been analyzed to date by their expanded size, long spacers, and diversity of introns. We have determined the mtDNA sequence of Chara vulgaris (Charophyceae), a green alga belonging to the charophycean order (Charales) that is thought to be the most closely related alga to land plants. This 67,737-bp mtDNA sequence, displaying 68 conserved genes and 27 introns, was compared with those of three angiosperms, the bryophyte Marchantia polymorpha, the charophycean alga Chaetosphaeridium globosum (Coleochaetales), and the green alga Mesostigma viride. Despite important differences in size and intron composition, Chara mtDNA strikingly resembles Marchantia mtDNA; for instance, all except 9 of 68 conserved genes lie within blocks of colinear sequences. Overall, our genome comparisons and phylogenetic analyses provide unequivocal support for a sister-group relationship between the Charales and the land plants. Only four introns in land plant mtDNAs appear to have been inherited vertically from a charalean algar ancestor. We infer that the common ancestor of green algae and land plants harbored a tightly packed, gene-rich, and relatively intron-poor mitochondrial genome. The group II introns in this ancestral genome appear to have spread to new mtDNA sites during the evolution of bryophytes and charalean green algae, accounting for part of the intron diversity found in Chara and land plant mitochondria.     

The Complete Mitochondrial Genome Sequence of the Hornwort <Emphasis Type="Italic">Megaceros</Emphasis><Emphasis Type="Italic">aenigmaticus</Emphasis> Shows a Mixed Mode of Conservative Yet Dynamic Evolution in Early Land Plant Mitochondrial Genomes

Land plants possess some of the most unusual mitochondrial genomes among eukaryotes. However, in early land plants these genomes resemble those of green and red algae or early eukaryotes. The question of when during land plant evolution the dramatic change in mtDNAs occurred remains unanswered. Here we report the first completely sequenced mitochondrial genome of the hornwort, Megaceros aenigmaticus, a member of the sister group of vascular plants. It is a circular molecule of 184,908 base pairs, with 32 protein genes, 3 rRNA genes, 17 tRNA genes, and 30 group II introns. The genome contains many genes arranged in the same order as in those of a liverwort, a moss, several green and red algae, and Reclinomonas americana, an early-branching eukaryote with the most ancestral form of mtDNA. In particular, the gene order between mtDNAs of the hornwort and Physcomitrella patens (moss) differs by only 8 inversions and translocations. However, the hornwort mtDNA possesses 4 derived features relative to green alga mtDNAs—increased genome size, RNA editing, intron gains, and gene losses—which were all likely acquired during the origin and early evolution of land plants. Overall, this genome and those of other 2 bryophytes show that mitochondrial genomes in early land plants, unlike their seed plant counterparts, exhibit a mixed mode of conservative yet dynamic evolution. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Libo Li and Bin Wang contributed equally to this work.     

Complete sequence of the mitochondrial DNA of the red alga Porphyra purpurea. Cyanobacterial introns and shared ancestry of red and green algae.

The mitochondrial DNA (mtDNA) of Porphyra purpurea, a circular-mapping genome of 36,753 bp, has been completely sequenced. A total of 57 densely packed genes has been identified, including the basic set typically found in animals and fungi, as well as seven genes characteristic of protist and plant mtDNAs and specifying ribosomal proteins and subunits of succinate:ubiquinone oxidoreductase. The mitochondrial large subunit rRNA gene contains two group II introns that are extraordinarily similar to those found in the cyanobacterium Calothrix sp, suggesting a recent lateral intron transfer between a bacterial and a mitochondrial genome. Notable features of P. purpurea mtDNA include the presence of two 291-bp inverted repeats that likely mediate homologous recombination, resulting in genome rearrangement, and of numerous sequence polymorphisms in the coding and intergenic regions. Comparative analysis of red algal mitochondrial genomes from five different, evolutionarily distant orders reveals that rhodophyte mtDNAs are unusually uniform in size and gene order. Finally, phylogenetic analyses provide strong evidence that red algae share a common ancestry with green algae and plants.     

Genome digging: insight into the mitochondrial genome of Homo

Background

A fraction of the Neanderthal mitochondrial genome sequence has a similarity with a 5,839-bp nuclear DNA sequence of mitochondrial origin (numt) on the human chromosome 1. This fact has never been interpreted. Although this phenomenon may be attributed to contamination and mosaic assembly of Neanderthal mtDNA from short sequencing reads, we explain the mysterious similarity by integration of this numt (mtAncestor-1) into the nuclear genome of the common ancestor of Neanderthals and modern humans not long before their reproductive split.

Principal Findings

Exploiting bioinformatics, we uncovered an additional numt (mtAncestor-2) with a high similarity to the Neanderthal mtDNA and indicated that both numts represent almost identical replicas of the mtDNA sequences ancestral to the mitochondrial genomes of Neanderthals and modern humans. In the proteins, encoded by mtDNA, the majority of amino acids distinguishing chimpanzees from humans and Neanderthals were acquired by the ancestral hominins. The overall rate of nonsynonymous evolution in Neanderthal mitochondrial protein-coding genes is not higher than in other lineages. The model incorporating the ancestral hominin mtDNA sequences estimates the average divergence age of the mtDNAs of Neanderthals and modern humans to be 450,000–485,000 years. The mtAncestor-1 and mtAncestor-2 sequences were incorporated into the nuclear genome approximately 620,000 years and 2,885,000 years ago, respectively.

Conclusions

This study provides the first insight into the evolution of the mitochondrial DNA in hominins ancestral to Neanderthals and humans. We hypothesize that mtAncestor-1 and mtAncestor-2 are likely to be molecular fossils of the mtDNAs of Homo heidelbergensis and a stem Homo lineage. The d_N/d_S dynamics suggests that the effective population size of extinct hominins was low. However, the hominin lineage ancestral to humans, Neanderthals and H. heidelbergensis, had a larger effective population size and possessed genetic diversity comparable with those of chimpanzee and gorilla.     

The mitochondrial genome of the stramenopile alga Chrysodidymus synuroideus. Complete sequence, gene content and genome organization

This is the first report of a complete mitochondrial genome sequence from a photosynthetic member of the stramenopiles, the chrysophyte alga Chrysodidymus synuroideus. The circular-mapping mitochondrial DNA (mtDNA) of 34 119 bp contains 58 densely packed genes (all without introns) and five unique open reading frames (ORFs). Protein genes code for components of respiratory chain complexes, ATP synthase and the mitoribosome, as well as one product of unknown function, encoded in many other protist mtDNAs (YMF16). In addition to small and large subunit ribosomal RNAs, 23 tRNAs are mtDNA-encoded, permitting translation of all codons present in protein-coding genes except ACN (Thr) and CGN (Arg). The missing tRNAs are assumed to be imported from the cytosol. Comparison of the C.synuroideus mtDNA with that of other stramenopiles allowed us to draw conclusions about mitochondrial genome organization, expression and evolution. First, we provide evidence that mitochondrial ORFs code for highly derived, unrecognizable versions of ribosomal or respiratory genes otherwise ‘missing’ in a particular mtDNA. Secondly, the observed constraints in mitochondrial genome rearrangements suggest operon-based, co-ordinated expression of genes functioning in common biological processes. Finally, stramenopile mtDNAs reveal an unexpectedly low variability in genome size and gene complement, testifying to substantial differences in the tempo of mtDNA evolution between major eukaryotic lineages.     

The Mitochondrial Genome of an Aquatic Plant,Spirodela polyrhiza

Background

Spirodela polyrhiza is a species of the order Alismatales, which represent the basal lineage of monocots with more ancestral features than the Poales. Its complete sequence of the mitochondrial (mt) genome could provide clues for the understanding of the evolution of mt genomes in plant.

Methods

Spirodela polyrhiza mt genome was sequenced from total genomic DNA without physical separation of chloroplast and nuclear DNA using the SOLiD platform. Using a genome copy number sensitive assembly algorithm, the mt genome was successfully assembled. Gap closure and accuracy was determined with PCR products sequenced with the dideoxy method.

Conclusions

This is the most compact monocot mitochondrial genome with 228,493 bp. A total of 57 genes encode 35 known proteins, 3 ribosomal RNAs, and 19 tRNAs that recognize 15 amino acids. There are about 600 RNA editing sites predicted and three lineage specific protein-coding-gene losses. The mitochondrial genes, pseudogenes, and other hypothetical genes (ORFs) cover 71,783 bp (31.0%) of the genome. Imported plastid DNA accounts for an additional 9,295 bp (4.1%) of the mitochondrial DNA. Absence of transposable element sequences suggests that very few nuclear sequences have migrated into Spirodela mtDNA. Phylogenetic analysis of conserved protein-coding genes suggests that Spirodela shares the common ancestor with other monocots, but there is no obvious synteny between Spirodela and rice mtDNAs. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly four-fifths of the Spirodela mitochondrial genome is of unknown origin and function. Although it contains a similar chloroplast DNA content and range of RNA editing as other monocots, it is void of nuclear insertions, active gene loss, and comprises large regions of sequences of unknown origin in non-coding regions. Moreover, the lack of synteny with known mitochondrial genomic sequences shed new light on the early evolution of monocot mitochondrial genomes.     

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.     

The Mitochondrial Genomes of Two Nematodes, Caenorhabditis Elegans and Ascaris Suum

The nucleotide sequences of the mitochondrial DNA (mtDNA) molecules of two nematodes, Caenorhabditis elegans [13,794 nucleotide pairs (ntp)], and Ascaris suum (14,284 ntp) are presented and compared. Each molecule contains the genes for two ribosomal RNAs (s-rRNA and l-rRNA), 22 transfer RNAs (tRNAs) and 12 proteins, all of which are transcribed in the same direction. The protein genes are the same as 12 of the 13 protein genes found in other metazoan mtDNAs: Cyt b, cytochrome b; COI-III, cytochrome c oxidase subunits I-III; ATPase6, Fo ATPase subunit 6; ND1-6 and 4L, NADH dehydrogenase subunits 1-6 and 4L: a gene for ATPase subunit 8, common to other metazoan mtDNAs, has not been identified in nematode mtDNAs. The C. elegans and A. suum mtDNA molecules both include an apparently noncoding sequence that contains runs of AT dinucleotides, and direct and inverted repeats (the AT region: 466 and 886 ntp, respectively). A second, apparently noncoding sequence in the C. elegans and A. suum mtDNA molecules (109 and 117 ntp, respectively) includes a single, hairpin-forming structure. There are only 38 and 89 other intergenic nucleotides in the C. elegans and A. suum mtDNAs, and no introns. Gene arrangements are identical in the C. elegans and A. suum mtDNA molecules except that the AT regions have different relative locations. However, the arrangement of genes in the two nematode mtDNAs differs extensively from gene arrangements in all other sequenced metazoan mtDNAs. Unusual features regarding nematode mitochondrial tRNA genes and mitochondrial protein gene initiation codons, previously described by us, are reviewed. In the C. elegans and A. suum mt-genetic codes, AGA and AGG specify serine, TGA specifies tryptophan and ATA specifies methionine. From considerations of amino acid and nucleotide sequence similarities it appears likely that the C. elegans and A. suum ancestral lines diverged close to the time of divergence of the cow and human ancestral lines, about 80 million years ago.     

The mitochondrial genome of the gymnosperm Cycas taitungensis contains a novel family of short interspersed elements, Bpu sequences, and abundant RNA editing sites

The mtDNA of Cycas taitungensis is a circular molecule of 414,903 bp, making it 2- to 6-fold larger than the known mtDNAs of charophytes and bryophytes, but similar to the average of 7 elucidated angiosperm mtDNAs. It is characterized by abundant RNA editing sites (1,084), more than twice the number found in the angiosperm mtDNAs. The A + T content of Cycas mtDNA is 53.1%, the lowest among known land plants. About 5% of the Cycas mtDNA is composed of a novel family of mobile elements, which we designated as "Bpu sequences." They share a consensus sequence of 36 bp with 2 terminal direct repeats (AAGG) and a recognition site for the Bpu 10I restriction endonuclease (CCTGAAGC). Comparison of the Cycas mtDNA with other plant mtDNAs revealed many new insights into the biology and evolution of land plant mtDNAs. For example, the noncoding sequences in mtDNAs have drastically expanded as land plants have evolved, with abrupt increases appearing in the bryophytes, and then in the seed plants. As a result, the genomic organizations of seed plant mtDNAs are much less compact than in other plants. Also, the Cycas mtDNA appears to have been exempted from the frequent gene loss observed in angiosperm mtDNAs. Similar to the angiosperms, the 3 Cycas genes nad1, nad2, and nad5 are disrupted by 5 group II intron squences, which have brought the genes into trans-splicing arrangements. The evolutionary origin and invasion/duplication mechanism of the Bpu sequences in Cycas mtDNA are hypothesized and discussed.     

Linear mitochondrial DNAs of yeasts: frequency of occurrence and general features.

In most yeast species, the mitochondrial DNA (mtDNA) has been reported to be a circular molecule. However, two cases of linear mtDNA with specific termini have previously been described. We examined the frequency of occurrence of linear forms of mtDNA among yeasts by pulsed-field gel electrophoresis. Among the 58 species from the genera Pichia and Williopsis that we examined, linear mtDNA was found with unexpectedly high frequency. Thirteen species contained a linear mtDNA, as confirmed by restriction mapping, and labeling, and electron microscopy. The mtDNAs from Pichia pijperi, Williopsis mrakii, and P. jadinii were studied in detail. In each case, the left and right terminal fragments shared homologous sequences. Between the terminal repeats, the order of mitochondrial genes was the same in all of the linear mtDNAs examined, despite a large variation of the genome size. This constancy of gene order is in contrast with the great variation of gene arrangement in circular mitochondrial genomes of yeasts. The coding sequences determined on several genes were highly homologous to those of the circular mtDNAs, suggesting that these two forms of mtDNA are not of distant origins.     

原生动物线粒体DNA的研究进展

目前已完成了多种原生动物mtDNA全序列的测定,对线粒体基因组有了较全面和深入的认识。本文总结了原生动物线粒体基因组结构、基因组成及基因表达等方面的研究进展,有助于进一步了解原始线粒体基因组的组成及其mtDNA进化,并为细胞的起源和真核生物进化的研究提供有价值的线索。     

The complete nucleotide sequence of a snake (Dinodon semicarinatus) mitochondrial genome with two identical control regions.

The 17,191-bp mitochondrial DNA (mtDNA) of a Japanese colubrid snake, akamata (Dinodon semicarinatus), was cloned and sequenced. The snake mtDNA has some peculiar features that were found in our previous study using polymerase chain reaction: duplicate control regions that have completely identical sequences over 1 kbp, translocation of tRNALeu(UUR) gene, shortened TpsiC arm for most tRNA genes, and a pseudogene for tRNAPro. Phylogenetic analysis of amino acid sequences of protein genes suggested an unusually high rate of molecular evolution in the snake compared to other vertebrates. Southern hybridization experiments using mtDNAs purified from multiple akamata individuals showed that the duplicate state of the control region is not a transient or unstable feature found in a particular individual, but that it stably occurs in mitochondrial genomes of the species. This may, therefore, be regarded as an unprecedented example of stable functional redundancy in animal mtDNA. However, some of the examined individuals contain a rather scanty proportion of heteroplasmic mtDNAs with an organization of genes distinct from that of the major mtDNA. The gene organization of the minor mtDNA is in agreement with one of models that we present to account for the concerted evolution of duplicate control regions.