The mitochondrial genome of the moss Physcomitrella patens sheds new light on mitochondrial evolution in land plants

The phylogenetic positions of bryophytes and charophytes, together with their genome features, are important for understanding early land plant evolution. Here we report the complete nucleotide sequence (105,340 bp) of the circular-mapping mitochondrial DNA of the moss Physcomitrella patens. Available evidence suggests that the multipartite structure of the mitochondrial genome in flowering plants does not occur in Physcomitrella. It contains genes for 3 rRNAs (rnl, rns, and rrn5), 24 tRNAs, and 42 conserved mitochondrial proteins (14 ribosomal proteins, 4 ccm proteins, 9 nicotinamide adenine dinucleotide dehydrogenase subunits, 5 ATPase subunits, 2 succinate dehydrogenase subunits, apocytochrome b, 3 cytochrome oxidase subunits, and 4 other proteins). We estimate that 5 tRNA genes are missing that might be encoded by the nuclear genome. The overall mitochondrial genome structure is similar in Physcomitrella, Chara vulgaris, Chaetosphaeridium globosum, and Marchantia polymorpha, with easily identifiable inversions and translocations. Significant synteny with angiosperm and chlorophyte mitochondrial genomes was not detected. Phylogenetic analysis of 18 conserved proteins suggests that the moss-liverwort clade is sister to angiosperms, which is consistent with a previous analysis of chloroplast genes but is not consistent with some analyses using mitochondrial sequences. In Physcomitrella, 27 introns are present within 16 genes. Nine of its intron positions are shared with angiosperms and 4 with Marchantia, which in turn shares only one intron position with angiosperms. The phylogenetic analysis as well as the syntenic structure suggest that the mitochondrial genomes of Physcomitrella and Marchantia retain prototype features among land plant mitochondrial genomes.     

The mitochondrial genome of the house centipede scutigera and the monophyly versus paraphyly of myriapods

Recent advances in molecular phylogenetics are continuously changing our perception of the phylogenetic relationships among the main arthropod lineages: crustaceans, hexapods, chelicerates, and myriapods. Besides the intrinsic interest in unraveling the evolution of the largest animal phylum, these studies are basic to an understanding of one of the major transitions in animal evolution-i.e., the conquest of land with all its associated structural and functional adaptations. Myriapods have been traditionally considered the closest relatives of hexapods, thus implying only one origin of terrestriality for the tracheate lineage, but this view is now challenged by molecular evidence. Sequence data available to date for centipedes and millipedes are very limited, and the taxon sampling is strongly biased. The most critical gap was the scutigeromorph centipedes, which are the sister group to all remaining Chilopoda from which they probably diverged in the Silurian if not earlier. We obtained the first complete mitochondrial sequence for a representative of this clade, the house centipede. In our phylogenetic analyses of the protein-coding genes in this mitochondrial genome, along with 16 further ones representing the other major arthropod clades plus two outgroups, the myriapods formed a clade with the chelicerates. This implies that water-to-land transition occurred at least three times (hexapods, myriapods, arachnids) during the evolution of the Arthropoda. In addition, in contrast to all previous studies, our best supported topologies favor paraphyly of the myriapods with respect to the chelicerates. This would increase to four the main events of land colonization in arthropods (once for centipedes, once for millipedes).     

Phylogeny, recombination, and mechanisms of stepwise mitochondrial genome reorganization in mantellid frogs from Madagascar

In Malagasy frogs of the family Mantellidae, the genus Mantellais known to possess highly reorganized mitochondrial (mt) genomeswith the following characteristics: 1) some rearranged genepositions, 2) 2 distinct genes and a pseudogene correspondingto the transfer RNA gene for methionine (trnM), and 3) 2 controlregions (CRs) with almost identical nucleotide sequences. Theseunique genomic features were observed concentrated between theduplicated CRs surrounding cytochrome b (cob) and nicotinamideadenine dinucleotide dehydrogenase subunit 2 (cnad2) genes.To elucidate the mechanisms and evolutionary pathway that yieldedthe derived genome condition, we surveyed the reorganized genomicportion for all 12 mantellid genera. Our results show that themt genomes of 7 genera retain the ancestral condition. In contrast,adding to Mantella, 4 genera of the subfamily Mantellinae, Blommersia,Guibemantis, Wakea, and Spinomantis, share several derived genomiccharacters. Furthermore, mt genomes of these mantellines showedadditional structural divergences, resulting in different genomeconditions between them. The high frequency of genomic reorganizationdoes not correlate with nucleotide substitution rate. The encounteredmt genomic conditions also suggest the occurrences of stepwisegene duplication and deletion events during the evolution ofmantellines. Simultaneously, the majority of duplication eventsseems to be mediated by general (homologous) or illegitimaterecombination, and general recombination also plays a role inconcerted sequence evolution between multiple CRs. Consideringour observations and recent conditional evidences, the followingoutlines can be expected for recombination processes in mt genomereorganization. 1) The CR is the "hot spot" of recombination;2) highly frequent recombination between CRs may be mediatedby a replication fork barrier lying in the CR; 3) general recombinationhas a potential to cause gene rearrangement in upstream regionsof multiple CRs as the results of gene conversion and unequalcrossing over processes. Our results also suggest that recombinationactivity is not a direct cause of convergent gene rearrangement;rather, homoplasious gene rearrangement seems to be mediatedby persistence of a copied genomic condition through severallineage splits and subsequent parallel deletions.     

Complete mtDNA sequence of the North American freshwater mussel,Lampsilis ornata (Unionidae): an examination of the evolution and phylogenetic utility of mitochondrial genome organization in Bivalvia (Mollusca)

Molluscs in general, and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. Two factors inhibiting our understanding the evolution of gene rearrangement in bivalves are inadequate taxonomic sampling and failure to examine gene order in a phylogenetic framework. Here, we report the first complete nucleotide sequence (16,060 bp) of the mitochondrial (mt) genome of a North American freshwater bivalve, Lampsilis ornata (Mollusca: Paleoheterodonta: Unionidae). Gene order and mt genome content is examined in a comparative phylogenetic framework for Lampsilis and five other bivalves, representing five families. Mitochondrial genome content is shown to vary by gene duplication and loss among taxa and between male and female mitotypes within a species. Although mt gene arrangement is highly variable among bivalves, when optimized on an independently derived phylogenetic hypothesis, it allows for the reconstruction of ancestral gene order states and indicates the potential phylogenetic utility of the data. However, the interpretation of reconstructed ancestral gene order states must take in to account both the accuracy of the phylogenetic estimation and the probability of character state change across the topology, such as the presence/absence of atp8 in bivalve lineages. We discuss what role, if any, doubly uniparental inheritance (DUI) and recombination between sexual mitotypes may play in influencing gene rearrangement of the mt genome in some bivalve lineages.     

Transfer of chloroplast genomic DNA to mitochondrial genome occurred at least 300 MYA

With the completion of the first gymnosperm mitochondrial genome (mtDNA) from Cycas taitungensis and the availability of more mtDNA taxa in the past 5 years, we have conducted a systematic analysis of DNA transfer from chloroplast genomes (cpDNAs) to mtDNAs (mtpts) in 11 plants, including 2 algae, 1 liverwort, 1 moss, 1 gymnosperm, 3 monocots, and 3 eudicots. By using shared gene order and boundaries between different mtpts as the criterion, the timing of cpDNA transfer during plant evolution was estimated from the phylogenetic tree reconstructed independently from concatenated protein-coding genes of 11 available mtDNAs. Several interesting findings emerged. First, frequent DNA transfer from cpDNA to mtDNA occurred at least as far back as the common ancestor of extant gymnosperms and angiosperms, about 300 MYA. The oldest mtpt is trnV(uac)-trnM(cau)-atpE-atpB-rbcL. Three other mtpts--psaA-psaB, rps19-trnH(gug)-rpl2-rpl23, and psbE-psbF--were dated to the common ancestor of extant angiosperms, at least 150 MYA. However, all protein-coding genes of mtpts have degenerated since their first transfer. Therefore, mtpts contribute nothing to the functioning of mtDNA but junk sequences. We discovered that the cpDNA transfers have occurred randomly at any positions of the cpDNAs. We provide strong evidence that the cp-derived tRNA-trnM(cau) is the only mtpt (1 out of 3 cp-derived tRNA shared by seed plants) truly transferred from cpDNA to mtDNA since the time of the common ancestor of extant gymnosperms and angiosperms. Our observations support the proposition of Richly and Leister (2004) that "primary insertions of organellar DNAs are large and then diverge and fragment over evolutionary time."     

The mitochondrial genome sequence and molecular phylogeny of the turkey, Meleagris gallopavo

The mitochondrial genome (mtGenome) has been little studied in the turkey ( Meleagris gallopavo ), a species for which there is no publicly available mtGenome sequence. Here, we used PCR-based methods with 19 pairs of primers designed from the chicken and other species to develop a complete turkey mtGenome sequence. The entire sequence (16 717 bp) of the turkey mtGenome was obtained, and it exhibited 85% similarity to the chicken mtGenome sequence. Thirteen genes and 24 RNAs (22 tRNAs and 2 rRNAs) were annotated. An mtGenome-based phylogenetic analysis indicated that the turkey is most closely related to the chicken, Gallus gallus , and quail, Corturnix japonica . Given the importance of the mtGenome, the present work adds to the growing genomic resources needed to define the genetic mechanisms that underlie some economically significant traits in the turkey.     

The complete mitochondrial genome of the onychophoran Epiperipatus biolleyi reveals a unique transfer RNA set and provides further support for the ecdysozoa hypothesis

Onychophora (velvet worms) play a crucial role in current discussions on position of arthropods. The ongoing Articulata/Ecdysozoa debate is in need of additional ground pattern characters for Panarthropoda (Arthropoda, Tardigrada, and Onychophora). Hence, Onychophora is an important outgroup taxon in resolving the relationships among arthropods, irrespective of whether morphological or molecular data are used. To date, there has been a noticeable lack of mitochondrial genome data from onychophorans. Here, we present the first complete mitochondrial genome sequence of an onychophoran, Epiperipatus biolleyi (Peripatidae), which shows several characteristic features. Specifically, the gene order is considerably different from that in other arthropods and other bilaterians. In addition, there is a lack of 9 tRNA genes usually present in bilaterian mitochondrial genomes. All these missing tRNAs have anticodon sequences corresponding to 4-fold degenerate codons, whereas the persisting 13 tRNAs all have anticodons pairing with 2-fold degenerate codons. Sequence-based phylogenetic analysis of the mitochondrial protein-coding genes provides a robust support for a clade consisting of Onychophora, Priapulida, and Arthropoda, which confirms the Ecdysozoa hypothesis. However, resolution of the internal ecdysozoan relationships suffers from a cluster of long-branching taxa (including Nematoda and Platyhelminthes) and a lack of data from Tardigrada and further nemathelminth taxa in addition to nematodes and priapulids.     

The mitochondrial genome of the lizard Calotes versicolor and a novel gene inversion in South Asian draconine agamids

A complete mitochondrial DNA (mtDNA) sequence was determinedfor the lizard Calotes versicolor (Reptilia; Agamidae). The16,670-bp genome with notable shorter genes for some protein-codingand tRNA genes had the same gene content as that found in othervertebrates. However, a novel gene arrangement was found inwhich the proline tRNA (trnP) gene is located in the light strandinstead of its typical heavy-strand position, providing thefirst known example of gene inversion in vertebrate mtDNAs.A segment of mtDNA encompassing the trnP gene and its flankinggenes and the control region was amplified and sequenced forvarious agamid taxa to investigate timing and mechanism of thegene inversion. The inverted trnP gene organization was sharedby all South Asian draconine agamids examined but by none ofthe other Asian and African agamids. Phylogenetic analyses includingclock-free Bayesian analyses for divergence time estimationsuggested a single occurrence of the gene inversion on a lineageleading to the draconine agamids during the Paleogene period.This gene inversion could not be explained by the tandem duplication/randomloss model for mitochondrial gene rearrangements. Our availablesequence data did not provide evidence for remolding of thetrnP gene by an anticodon switch in a duplicated tRNA gene.Based on results of sequence comparisons and other circumstantialevidence, we hypothesize that inversion of the trnP gene wasoriginally mediated by a homologous DNA recombination and thatthe de novo gene organization that does not disrupt expressionof mitochondrial genes has been maintained in draconine mtDNAsfor such a long period of time.     

A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.     

Physical map and gene organization of the mitochondrial genome from the unicellular green alga Platymonas (Tetraselmis) subcordiformis (Prasinophyceae)

the entire mitochondrial genome (mt genome) of the unicellular green alga Platymonas subcordiformis (synonym Tetraselmis subcordiformis; Prasinophyceae) was cloned and a physical map for the four restriction enzymes Hind III, Eco RI, Bgl II and Xba I was constructed. The mt genome of P. subcordiformis is a 42.8 kb circular molecule, coding for at least 23 genes. Hybridization and sequence analysis revealed the presence of a ca. 1.5 kb inverted repeat on the mt genome of P. subcordiformis. Phylogenetic analyses based on sequences of several coxI genes were carried out. Our data indicate that mitochondria from P. subcordiformis and from land plants form a natural, monophyletic group.     

The linear 20 kb mitochondrial genome of Pandorina morum (Volvocaceae,Chlorophyta)

A physical restriction map of the mitochondrial genome from one clone (TCC 854) of the sexually isolated populations (syngens) of the morphologically uniform species Pandorina morum Bory has been constructed using restriction endonucleases Ava I, Bam HI, Bgl II, Eco RI, Kpn I, and Pst I. The 20 kb linear genome can easily be separated from plastid DNA, nuclear satellite rDNA, and main band (nuclear) DNA on a Hoechst/CsCl buoyant density gradient. The Pandorina mitochondrial DNA shows sufficient similarity to the 16 kb mitochondrial genome of Chlamydomonas reinhardtii to cross-hybridize, and also hybridizes with a probe containing maize mitochondrial 18S rRNA genes. Double digests, self-probing, and Bal31 exonuclease experiments suggest that 1.8 to 3.3 kb of sequence is repeated at each end of the genome as an inverted repeat. Mitochondrial genome sizes of other P. morum syngens were found to range from ca. 20 to ca. 38 kb. The mitochondrial genome should be valuable for taxonomic studies; it can be used for comparative organellar studies; and it should be of interest to compare with that of other plant and animal mitochondrial genomes.     

The mitochondrial genome of the spinyhead croaker Collichthys lucida: genome organization and phylogenetic consideration

The complete mitochondrial genome of the spiny head croaker Collichthys lucida was determined in the present study. The mitochondrial DNA was 16,442 base pairs in length, and contained 13 protein coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and one major non-coding control region, with the content and order of genes being similar to those in typical teleosts. Most of the genes of C. lucida were encoded on the H-strand, while the ND6 and eight tRNA (Gln, Ala, Asn, Cys, Tyr, Ser (UCN), Glu and Pro) genes were encoded on the L-strand. The reading frames of two pairs of genes overlapped: ATPase 8 and 6 and ND4L and ND4 by ten and seven nucleotides, respectively. The control region was unusually short at only 768bp, and absence of typical conserved blocks (CSB-D, CSB-E, and CSB-F). Phylogenetic analyses indicated that C. lucida was located in the cluster of fish species from the family Sciaenidae, supporting the traditional taxonomic classification of fish, and in the cluster of Serranidae, the divergence time in Plectropomus leopardus is longer than that among its coordinal species. On the other hand, phylogenetic analyses do not support the monophyletic of family Centracanthidae and genera Larimichthys and Collichthys, which is against the morphological results.     

A Mitochondrial Genome Sequence of the Tibetan Antelope （ Pantholops hodgsonii ）

To investigate genetic mechanisms of high altitude adaptations of native mammals on the Tibetan Plateau, we compared mitochondrial sequences of the endangered Pantholops hodgsonii with its lowland distant relatives Ovis ames and Capra hircus, as well as other mammals. The complete mitochondrial genome of P. hodgsonii （16,498 bp） revealed a similar gene order as of other mammals. Because of tandem duplications, the control region of P. hodgsonii mitochondrial genome is shorter than those of O. ames and C. hircus, but longer than those of Bos species. Phylogenetic analysis based on alignments of the entire cytochrome b genes suggested that P. hodgsonii is more closely related to O. ames and C. hircus, rather than to species of the Antilopinae subfamily. The estimated divergence time between P. hodgsonii and O. ames is about 2.25 million years ago. Eutther analysis on natural selection indicated that the COXI （cytochrome c oxidase subunit I） gene was under positive selection in P. hodgsonii and Bos grunniens. Considering the same climates and environments shared by these two mammalian species, we proposed that the mitochondrial COXI gene is probably relevant for these native mammals to adapt the high altitude environment unique to the Tibetan Plateau.     

Bird mitochondrial gene order: insight from 3 warbler mitochondrial genomes

Two main gene orders exist in birds: the ancestral gene order and the remnant control region (CR) 2 gene order. These gene orders differ by the presence of 1 or 2 copies of the CR, respectively. Among songbirds, Oscines were thought to follow the ancestral gene order, with the exception of the lyrebird and Phylloscopus warblers. Here, we determined the complete mitochondrial genome sequence of 3 non-Phylloscopus warblers species and found that the blackcap (Sylvia atricapilla) and the reed warbler (Acrocephalus scirpaceus) have 2 almost identical copies of the CR, whereas the eastern orphean warbler (Sylvia crassirostris) follows the remnant CR 2 gene order. Our results contradict previous studies suggesting that Acrocephalus and most sylvioid warblers exhibit the ancestral gene order. We were able to trace this contradiction to a misidentification of gene order from polymerase chain reaction length determination. We thus suggest that passerine gene order evolution needs to be revised.