Complete Mitochondrial DNA Sequences of Six Snakes: Phylogenetic Relationships and Molecular Evolution of Genomic Features

Complete mitochondrial DNA (mtDNA) sequences were determined for representative species from six snake families: the acrochordid little file snake, the bold boa constrictor, the cylindrophiid red pipe snake, the viperid himehabu, the pythonid ball python, and the xenopeltid sunbeam snake. Thirteen protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 2 control regions were identified in these mtDNAs. Duplication of the control region and translocation of the tRNA^Leu gene were two notable features of the snake mtDNAs. The duplicate control regions had nearly identical nucleotide sequences within species but they were divergent among species, suggesting concerted sequence evolution of the two control regions. In addition, the duplicate control regions appear to have facilitated an interchange of some flanking tRNA genes in the viperid lineage. Phylogenetic analyses were conducted using a large number of sites (9570 sites in total) derived from the complete mtDNA sequences. Our data strongly suggested a new phylogenetic relationship among the major families of snakes: ((((Viperidae, Colubridae), Acrochordidae), (((Pythonidae, Xenopeltidae), Cylindrophiidae), Boidae)), Leptotyphlopidae). This conclusion was distinct from a widely accepted view based on morphological characters in denying the sister-group relationship of boids and pythonids, as well as the basal divergence of nonmacrostomatan cylindrophiids. These results imply the significance to reconstruct the snake phylogeny with ample molecular data, such as those from complete mtDNA sequences.[Reviewing Editor: Dr. Bill Ballard]     

Mitochondrial DNA sequences of five squamates: phylogenetic affiliation of snakes.

Complete or nearly complete mitochondrial DNA sequences were determined from four lizards (Western fence lizard, Warren's spinytail lizard, Terrestrial arboreal alligator lizard, and Chinese crocodile lizard) and a snake (Texas blind snake). These genomes had a typical gene organization found in those of most mammals and fishes, except for a translocation of the glutamine tRNA gene in the blind snake and a tandem duplication of the threonine and proline tRNA genes in the spinytail lizard. Although previous work showed the existence of duplicate control regions in mitochondrial DNAs of several snakes, the blind snake did not have this characteristic. Phylogenetic analyses based on different tree-building methods consistently supported that the blind snake and a colubrid snake (akamata) make a sister clade relative to all the lizard taxa from six different families. An alternative hypothesis that snakes evolved from a lineage of varanoids was not favored and nearly statistically rejected by the Kishino-Hasegawa test. It is therefore likely that the apparent similarity of the tongue structure between snakes and varanoids independently evolved and that the duplication of the control region occurred on a snake lineage after divergence of the blind snake.     

Gene rearrangements in snake mitochondrial genomes: highly concerted evolution of control-region-like sequences duplicated and inserted into a tRNA gene cluster

Mitochondrial DNA (mtDNA) regions corresponding to two major tRNA gene clusters were amplified and sequenced for the Japanese pit viper, himehabu. In one of these clusters, which in most vertebrates characterized to date contains three tightly connected genes for tRNA(Ile), and tRNA(Gln), and tRNA(Met), a sequence of approximately 1.3 kb was found to be inserted between the genes for tRNA(Ile) and tRNA(Gln). The insert consists of a control-region-like sequence possessing some conserved sequence blocks, and short flanking sequences which may be folded into tRNA(Pro), tRNA(Phe), and tRNA(Leu) genes. Several other snakes belonging to different families were also found to possess a control-region-like sequence and tRNA(Leu) gene between the tRNA(Ile)and tRNA(Gln) genes. We also sequenced a region surrounded by genes for cytochrome b and 12S rRNA, where the control region and genes for tRNA(Pro) and tRNA(Phe) are normally located in the mtDNAs of most vertebrates. In this region of three examined snakes, a control-region- like sequence exists that is almost completely identical to the one found between the tRNA(Ile) and tRNA(Gln) genes. The mtDNAs of these snakes thus possess two nearly identical control-region-like sequences which are otherwise divergent to a large extent between the species. These results suggest that the duplicate state of the control-region- like sequences has long persisted in snake mtDNAs, possibly since the original insertion of the control-region-like sequence and tRNA(Leu) gene into the tRNA gene cluster, which occurred in the early stage of the divergence of snakes. It is also suggested that the duplicated control-region-like sequences at two distant locations of mtDNA have evolved concertedly by a mechanism such as frequent gene conversion. The secondary structures of the determined tRNA genes point to the operation of simplification pressure on the T psi C arm of snake mitochondrial tRNAs.      

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.     

Genome structure and gene content in protist mitochondrial DNAs.

Although the collection of completely sequenced mitochondrial genomes is expanding rapidly, only recently has a phylogenetically broad representation of mtDNA sequences from protists (mostly unicellular eukaryotes) become available. This review surveys the 23 complete protist mtDNA sequences that have been determined to date, commenting on such aspects as mitochondrial genome structure, gene content, ribosomal RNA, introns, transfer RNAs and the genetic code and phylogenetic implications. We also illustrate the utility of a comparative genomics approach to gene identification by providing evidence that orfB in plant and protist mtDNAs is the homolog of atp8 , the gene in animal and fungal mtDNA that encodes subunit 8 of the F0portion of mitochondrial ATP synthase. Although several protist mtDNAs, like those of animals and most fungi, are seen to be highly derived, others appear to be have retained a number of features of the ancestral, proto-mitochondrial genome. Some of these ancestral features are also shared with plant mtDNA, although the latter have evidently expanded considerably in size, if not in gene content, in the course of evolution. Comparative analysis of protist mtDNAs is providing a new perspective on mtDNA evolution: how the original mitochondrial genome was organized, what genes it contained, and in what ways it must have changed in different eukaryotic phyla.     

Diversity in organization and the origin of gene orders in the mitochondrial DNA molecules of the genus Saccharomyces

Sequencing of the Saccharomyces cerevisiae nuclear and mitochondrial genomes provided a new background for studies on the evolution of the genomes. In this study, mitochondrial genomes of a number of Saccharomyces yeasts were mapped by restriction enzyme analysis, the orders of the genes were determined, and two of the genes were sequenced. The genome organization, i.e., the size, presence of intergenic sequences, and gene order, as well as polymorphism within the coding regions, indicate that Saccharomyces mtDNA molecules are dynamic structures and have undergone numerous changes during their evolution. Since the separation and sexual isolation of different yeast lineages, the coding parts have been accumulating point mutations, presumably in a linear manner with the passage of time. However, the accumulation of other changes may not have been a simple function of time. Larger mtDNA molecules belonging to Saccharomyces sensu stricto yeasts have acquired extensive intergenic sequences, including guanosine-cytosine-rich clusters, and apparently have rearranged the gene order at higher rates than smaller mtDNAs belonging to the Saccharomyces sensu lato yeasts. While within the sensu stricto group transposition has been a predominant mechanism for the creation of novel gene orders, the sensu lato yeasts could have used both transposition- and inversion-based mechanisms.     

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.     

Mitochondrial genome of the Komodo dragon: efficient sequencing method with reptile-oriented primers and novel gene rearrangements.

The mitochondrial genome of the Komodo dragon (Varanus komodoensis) was nearly completely sequenced, except for two highly repetitive noncoding regions. An efficient sequencing method for squamate mitochondrial genomes was established by combining the long polymerase chain reaction (PCR) technology and a set of reptile-oriented primers designed for nested PCR amplifications. It was found that the mitochondrial genome had novel gene arrangements in which genes from NADH dehydrogenase subunit 6 to proline tRNA were extensively shuffled with duplicate control regions. These control regions had 99% sequence similarity over 700 bp. Although snake mitochondrial genomes are also known to possess duplicate control regions with nearly identical sequences, the location of the second control region suggested independent occurrence of the duplication on lineages leading to snakes and the Komodo dragon. Another feature of the mitochondrial genome of the Komodo dragon was the considerable number of tandem repeats, including sequences with a strong secondary structure, as a possible site for the slipped-strand mispairing in replication. These observations are consistent with hypotheses that tandem duplications via the slipped-strand mispairing may induce mitochondrial gene rearrangements and may serve to maintain similar copies of the control region.     

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.     

Mitochondrial genomes of two demosponges provide insights into an early stage of animal evolution

Mitochondrial DNA (mtDNA) of multicellular animals (Metazoa) is typically a small ( approximately 16 kbp), circular-mapping molecule that encodes 37 tightly packed genes. The structures of mtDNA-encoded transfer RNAs (tRNAs) and ribosomal RNAs (rRNAs) are usually highly unorthodox, and proteins are translated with multiple deviations from the standard genetic code. In contrast, mtDNA of the choanoflagellate Monosiga brevicollis, the closest unicellular relative of animals, is four times larger, contains 1.5 times as many genes, and lacks mentioned peculiarities of animal mtDNA. To investigate the evolutionary transition that led to the specific organization of metazoan mtDNA, we determined complete mitochondrial sequences from the demosponges Geodia neptuni and Tethya actinia, two representatives of the most basal animal phylum, the Porifera. We found that poriferan mtDNAs resemble those of other animals in their compact organization, lack of introns, and a well-conserved animal-like gene order. Yet, they contain several extra genes, encode bacterial-like rRNAs and tRNAs, and use a minimally derived genetic code. Our findings suggest that the evolution of the typical metazoan mtDNA has been a multistep process in which the compact genome organization and the reduced gene content were established prior to the reduction of tRNA and rRNA structures and the introduction of multiple changes of the translation code.     

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.     

Complete sequence of the mitochondrial genome of Tetrahymena pyriformis and comparison with Paramecium aurelia mitochondrial DNA

We report the complete nucleotide sequence of the Tetrahymena pyriformis mitochondrial genome and a comparison of its gene content and organization with that of Paramecium aurelia mtDNA. T. pyriformis mtDNA is a linear molecule of 47,172 bp (78.7 % A+T) excluding telomeric sequences (identical tandem repeats of 31 bp at each end of the genome). In addition to genes encoding the previously described bipartite small and large subunit rRNAs, the T. pyriformis mitochondrial genome contains 21 protein-coding genes that are clearly homologous to genes of defined function in other mtDNAs, including one (yejR) that specifies a component of a cytochrome c biogenesis pathway. As well, T. pyriformis mtDNA contains 22 open reading frames of unknown function larger than 60 codons, potentially specifying proteins ranging in size from 74 to 1386 amino acid residues. A total of 13 of these open reading frames ("ciliate-specific") are found in P. aurelia mtDNA, whereas the remaining nine appear to be unique to T. pyriformis; however, of the latter, five are positionally equivalent and of similar size in the two ciliate mitochondrial genomes, suggesting they may also be homologous, even though this is not evident from sequence comparisons. Only eight tRNA genes encoding seven distinct tRNAs are found in T. pyriformis mtDNA, formally confirming a long-standing proposal that most T. pyriformis mitochondrial tRNAs are nucleus-encoded species imported from the cytosol. Atypical features of mitochondrial gene organization and expression in T. pyriformis mtDNA include split and rearranged large subunit rRNA genes, as well as a split nad1 gene (encoding subunit 1 of NADH dehydrogenase of respiratory complex I) whose two segments are located on and transcribed from opposite strands, as is also the case in P. aurelia. Gene content and arrangement are very similar in T. pyriformis and P. aurelia mtDNAs, the two differing by a limited number of duplication, inversion and rearrangement events. Phylogenetic analyses using concatenated sequences of several mtDNA-encoded proteins provide high bootstrap support for the monophyly of alveolates (ciliates, dinoflagellates and apicomplexans) and slime molds.     

Nucleotide Sequence and Gene Organization of the Starfish Asterina Pectinifera Mitochondrial Genome

The 16,260-bp mitochondrial DNA (mtDNA) from the starfish Asterina pectinifera has been sequenced. The genes for 13 proteins, two rRNAs and 22 tRNAs are organized in an extremely economical fashion, similar to those of other animal mtDNAs, with some of the genes overlapping each other. The gene organization is the same as that for another echinoderm, sea urchin, except for the inversion of a 4.6-kb segment that contains genes for two proteins, 13 tRNAs and the 16S rRNA. Judging from the organization of the protein coding genes, mammalian mtDNAs resemble the sea urchin mtDNA more than that of the starfish. The region around the 3' end of the 12S rRNA gene of the starfish shows a high similarity with those for vertebrates. This region encodes a possible stem and loop structure; similar potential structures occur in this region of vertebrate mtDNAs and also in nonmitochondrial small subunit rRNA. A similar stem and loop structure is also found at the 3' end of the 16S rRNA genes in A. pectinifera, in another starfish Pisaster ochraceus, in vertebrates and in Drosophila, but not in sea urchins. The full sequence data confirm the presumption that AGA/AGG, AUA and AAA codons, respectively, code for serine, isoleucine, and asparagine in the starfish mitochondria, and that AGA/AGG codons are read by tRNA(GCU)(Ser), which possesses a truncated dihydrouridine arm, that was previously suggested from a partial mtDNA sequence. The structural characteristics of tRNAs and possible mechanisms for the change in the mitochondrial genetic code are also discussed.     

A comparative genomics approach to the evolution of eukaryotes and their mitochondria.

The Organelle Genome Megasequencing Program (OGMP) investigates mitochondrial genome diversity and evolution by systematically determining the complete mitochondrial DNA (mtDNA) sequences of a phylogenetically broad selection of protists. The mtDNAs of lower fungi and choanoflagellates are being analyzed by the Fungal Mitochondrial Genome Project (FMGP), a sister project to the OGMP. Some of the most interesting protists include the jakobid flagellates Reclinomonas americana, Malawimonas jakobiformis, and Jakoba libera, which share ultrastructural similarities with amitochondriate retortamonads, and harbor mitochondrial genes not seen before in mtDNAs of other organisms. In R. americana and J. libera, gene clusters are found that resemble, to an unprecedented degree, the contiguous ribosomal protein operons str, S10, spc, and alpha of eubacteria. In addition, their mtDNAs code for an RNase P RNA that displays all the elements of a bacterial minimum consensus structure. This structure has been instrumental in detecting the rnpB gene in additional protists. Gene repertoire and gene order comparisons as well as multiple-gene phylogenies support the view of a single endosymbiotic origin of mitochondria, whose closest extant relatives are Rickettsia-type alpha-Proteobacteria.     

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.