Mitochondrial gene order is not conserved in arthropods: prostriate and metastriate tick mitochondrial genomes

The entire mitochondrial genome was sequenced in a prostriate tick, Ixodes hexagonus, and a metastriate tick, Rhipicephalus sanguineus. Both genomes encode 22 tRNAs, 13 proteins, and two ribosomal RNAs. Prostriate ticks are basal members of Ixodidae and have the same gene order as Limulus polyphemus. In contrast, in R. sanguineus, a block of genes encoding NADH dehydrogenase subunit 1 (ND1), tRNA(Leu)(UUR), tRNA(Leu)(CUN), 16S rDNA, tRNA(Val), 12S rDNA, the control region, and the tRNA(Ile) and tRNA(Gln) have translocated to a position between the tRNA(Glu) and tRNA(Phe) genes. The tRNA(Cys) gene has translocated between the control region and the tRNA(Met) gene, and the tRNA(Leu)(CUN) gene has translocated between the tRNA(Ser)(UCN) gene and the control region. Furthermore, the control region is duplicated, and both copies undergo concerted evolution. Primers that flank these rearrangements confirm that this gene order is conserved in all metastriate ticks examined. Correspondence analysis of amino acid and codon use in the two ticks and in nine other arthropod mitochondrial genomes indicate a strong bias in R. sanguineus towards amino acids encoded by AT-rich codons.      

Complete nucleotide sequence and gene organization of the mitochondrial genome of Paa spinosa (Anura: Ranoidae)

The mt genome of Paa spinosa (Anura: Ranoidae) is a circular molecule of 18,012 bp in length, containing 38 genes (including an extra copy of tRNA-Met gene). This mt genome is characterized by three distinctive features: a cluster of rearranged tRNA genes (LTPF tRNA gene cluster), a tandem duplication of tRNA-Met gene (Met1 and Met2), and distinct repeat regions at both 5′ and 3′-sides in the control region. Comparing the locations and the sequences of all tRNA-Met genes among Ranoidae, and constructing NJ tree of the nucleotide of those tRNA-Met genes, we suggested a tandem duplication of tRNA-Met gene can be regarded as a synapomorphy of Dicroglossinae. To further investigate the phylogenetic relationships of anurans, phylogenetic analyses (BI, ML and MP) based on the nucleotide dataset and the corresponding amino acid dataset of 11 protein-coding genes (except ND5 and ATP8) arrived at the similar topology.     

Complete mitochondrial genome sequence of the Korean water deer Hydropotes inermis argyropus (Cervidae, Capreolinae)

The complete mitochondrial genome sequence of Hydropotes inermis argyropus consists of 13 protein-coding, 22 tRNA, and two rRNA genes, and 1 control region (CR). Three overlaps among the 13 protein-coding genes were found: ATP8/ATP6, ND4L/ND4, and ND5/ND6. The CR was located between the tRNA-Pro and tRNA-Phe genes and is 928 bp in length. The typical conserved domains, such as TAS and CSB, were identified in the CR.     

Complete nucleotide sequences of mitochondrial genomes of two solitary entoprocts, Loxocorone allax and Loxosomella aloxiata: implications for lophotrochozoan phylogeny

The complete nucleotide sequences of the mitochondrial (mt) genomes of the entoprocts Loxocorone allax and Loxosomella aloxiata were determined. Both species carry the typical gene set of metazoan mt genomes and have similar organizations of their mt genes. However, they show differences in the positions of two tRNA(Leu) genes. Additionally, the tRNA(Val) gene, and half of the long non-coding region, is duplicated and inverted in the Loxos. aloxiata mt genome. The initiation codon of the Loxos. aloxiata cytochrome oxidase subunit I gene is expected to be ACG rather than AUG. The mt gene organizations in these two entoproct species most closely resemble those of mollusks such as Katharina tunicata and Octopus vulgaris, which have the most evolutionarily conserved mt gene organization reported to date in mollusks. Analyses of the mt gene organization in the lophotrochozoan phyla (Annelida, Brachiopoda, Echiura, Entoprocta, Mollusca, Nemertea, and Phoronida) suggested a close phylogenetic relationship between Brachiopoda, Annelida, and Echiura. However, Phoronida was excluded from this grouping. Molecular phylogenetic analyses based on the sequences of mt protein-coding genes suggested a possible close relationship between Entoprocta and Phoronida, and a close relationship among Brachiopoda, Annelida, and Echiura.     

Long-term conservation of six duplicated structural genes in cephalopod mitochondrial genomes

The complete nucleotide sequences of the mitochondrial (mt) genomes of three cephalopods, Octopus vulgaris (Octopodiformes, Octopoda, Incirrata), Todarodes pacificus (Decapodiformes, Oegopsida, Ommastrephidae), and Watasenia scintillans (Decapodiformes, Oegopsida, Enoploteuthidae), were determined. These three mt genomes encode the standard set of metazoan mt genes. However, W. scintillans and T. pacificus mt genomes share duplications of the longest noncoding region, three cytochrome oxidase subunit genes and two ATP synthase subunit genes, and the tRNA(Asp) gene. Southern hybridization analysis of the W. scintillans mt genome shows that this single genome carries both duplicated regions. The near-identical sequence of the duplicates suggests that there are certain concerted evolutionary mechanisms, at least in cephalopod mitochondria. Molecular phylogenetic analyses of mt protein genes are suggestive, although not statistically significantly so, of a monophyletic relationship between W. scintillans and T. pacificus.     

Complete nucleotide sequence of the mitochondrial genome of a Malagasy poison frog Mantella madagascariensis: evolutionary implications on mitochondrial genomes of higher anuran groups

We determined the complete nucleotide sequence of the mitochondrial (mt) genome of a Malagasy poison frog, Mantella madagascariensis (family Mantellidae), and partial sequences of two Mantella (M. baroni and M. bernhardi) and two additional mantellid species (Boophis madagascariensis and Mantidactylus cf. ulcerosus). The M. madagascariensis genome was shown to be the largest (23kbp) of all vertebrate mtDNAs investigated so far. Furthermore, the following unique features were revealed: (1) the positions of some genes and gene regions were rearranged compared to mitochondrial genomes typical for vertebrates and other anuran groups, (2) two distinct genes and a pseudogene corresponding to transfer RNA gene for methionine (tRNA-Met) were encoded, and (3) two control regions with very high sequence homology were present. These features were shared by the two other Mantella species but not the other mantellid species, indicating dynamic genome reorganization in a common ancestor linage before divergence of the Mantella genus. The reorganization pathway could be explained by a model of gene duplication and deletion. Duplication and deletion events also seem to have been responsible for concerted sequence evolution of the control regions in Mantella mt genomes. It is also suggested that the pseudo tRNA-Met gene sustained for a long time in Mantella mt genomes possibly functions as a punctuation marker for NADH dehydrogenase subunit (ND) 2 mRNA processing. Phylogenetic analyses employing a large sequence data set of mt genes supported the monophyly of Mantellidae and Rhacophoridae and other recent phylogenetic views for ranoid frogs. The resultant phylogenetic relationship also suggested parallel occurrence of two tRNA-Met genes, duplicated control regions, and ND5 gene translocation in independent ranoid lineages.     

The Australian fresh water isopod (Phreatoicidea: Isopoda) allows insights into the early mitogenomic evolution of isopods

The complete mitochondrial (mt) genome sequence of the Australian fresh water isopod Eophreatoicus sp.-14 has been determined. The new species is a member of the taxon Phreatoicidea, a clade of particular interest, as it is often regarded as the sister group to all other Isopoda. Although the overall genome organization of Eophreatoicus sp.-14 conforms to the typical state of Metazoa—it is a circular ring of DNA hosting the usual 37 genes and one major non-coding region—it bears a number of derived characters that fall within the scope of “genome morphology”. Earlier studies have indicated that the isopod mitochondrial gene order is not as conserved as that of other crustaceans. Indeed, the mt genome of Eophreatoicus sp.-14 shows an inversion of seven genes (including cox1), which is as far as we know unique. Even more interesting is the derived arrangement of nad1, trnL(CUN), rrnS, control region, cob, trnT, nad5 and trnF that is shared by nearly all available isopod mt genomes. A striking feature is the close proximity of the rearranged genes to the mt control region. Inferable gene translocation events are, however, more suitable to trace the evolution of mt genomes. Genes like nad1/trnL(CUN) and nad5/trnF, which retained their adjacent position after being rearranged, were most likely translocated together. A very good example for the need to understand the mechanisms of translocations is the remolding of trnL(UUR) to trnL(CUN). Both tRNA genes are adjacent and have a high sequence similarity, probably the result of a gene duplication and subsequent anticodon mutation. Modified secondary structures were found in three tRNAs of Eophreatoicus sp.-14, which are all characterized by the loss of the DHU-arm. This is common to crustaceans for tRNA Serine(AGY), while the arm-loss in tRNA Cysteine within Malacostraca is only shared by other isopods. Modification of the third tRNA, Isoleucine, is not known from any other related species. Nucleotide frequencies of genes have been found to be indirectly correlated to the orientation of the mitochondrial replication process. In Eophreatoicus sp.-14 and in other Isopoda the associated nucleotide bias is inversed to the state of other Malacostraca. This is a strong indication for an inversion of the control region that most likely evolved in the isopod ancestor.     

The mitochondrial genome of the predatory mite Metaseiulus occidentalis (Arthropoda: Chelicerata: Acari: Phytoseiidae) is unexpectedly large and contains several novel features

The complete mitochondrial genome of the phytoseiid Metaseiulus occidentalis (Arthropoda: Chelicerata: Acari: Phytoseiidae) has been sequenced. It is 24,961 bp in length and contains a 14,695-bp unique region, a 345-bp triplicated region, and a 9921-bp duplicated region, in that order. The A+T content of the unique region is 76.9% and contains 11 protein coding (COI-III; ATP6-8; CytB; ND1, 2, 4, 5 and 4L), two ribosomal RNA (srRNA and lrRNA), 22 transfer RNA (tRNA) genes, and two copies of D-loop control sequence. Two genes (ND3 and 6) appear to be missing, but there is a large intergenic spacer (390 bp) present, which could contain ND3 if a different codon usage is employed. The gene order is completely different from the pattern in all other known chelicerates, including the horseshoe crab Limulus polyphemus [Lavrov et al., Mol. Biol. Evol., 2000; 17:813-824]. All the inferred tRNA genes are missing the TPsiC arm, but this arm has fused with the variable arm to generate a TV replacement loop. The duplicated region (9921 bp) contains 18 genes in the same order as in the unique region from CytB to tRNA-His, plus one copy of D-loop control sequence (311 bp) and a partial tRNA-Leu2 sequence (34 bp). The small triplicated region (345 bp) contains a D-loop control sequence (311 bp) and a partial tRNA-Leu2 sequence (34 bp). Because of these anomalies, amplifying sequences posed technical difficulties, but were accomplished by using a primer-walking strategy and increasing the AT content to 75% in the high-fidelity PCR dNTP mix. This is the first phytoseiid mitochondrial genome to be completely sequenced and the largest (25 kb) detected from the Chelicerata.     

The complete mitochondrial genome of the mackerel icefish,Champsocephalus gunnari (Actinopterygii: Channichthyidae), with reference to the evolution of mitochondrial genomes in Antarctic notothenioids

The mackerel icefish (Champsocephalus gunnari Lönnberg, 1905) is a ray‐finned fish living in the Southern Ocean around Antarctica. We sequenced the complete mitochondrial (mt) genome of the mackerel icefish and a segment from cytochrome b to the control region (CR) in 32 individuals. The mt genome of the mackerel icefish was rearranged, containing two nicotinamide adenine dinucleotide (reduced form) dehydrogenase subunit 6 (ND6), two tRNA^Glu, and two CRs. However, variations in numbers of ND6 and tRNA^Glu were observed amongst individuals. These variations included type 1 (containing two ND6 and two tRNA^Glu), type 2 (containing one ND6, one incomplete ND6, and one tRNA^Glu), and type 3 (containing one ND6 and one tRNA^Glu). The gene orders of types 1 and 2, and variations in numbers of ND6 and tRNA^Glu were not previously found in any Antarctic notothenioids, whereas type 3 is the same as that of Racovitzia glacialis. Phylogenetic analyses of CR DNA sequences showed that duplicated CRs of the same species formed a monophyletic group, suggesting that duplication of CRs occurred in each species. The frequent duplication of mt genomes in Antarctic notothenioids is an unusual feature in vertebrates. We propose that interspecific hybridization and impairment of mismatch repair might account for the high frequency of gene duplications and rearrangement of mt genomes in Antarctic notothenioids.     

A unique tRNA gene family and a novel,highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS^− AGY and trnQ genes, a trnM isomerism, but it lacks trnS^− CUN. We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP → trnS^− AGY) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, “Oxidored_q1_N”, which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.     

Complete mitochondiral genome of the Korean red fox Vulpes vulpes (Carnivora, Canidae)

The complete mitochondrial genome sequence of Vulpes vulpes consists of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 control region (CR). CR is located between the tRNA-Pro and tRNA-Phe genes and is 1173?base pairs (bp) in length. It consists of a short non-repetitive sequence followed by 8-bp 5'-ACACACGT-3' tandem repeat between conserved sequence black I and conserved sequence black II.     

Molecular characterization of a cloned dolphin mitochondrial genome

Summary DNA clones have been isolated that span the complete mitochondrial (mt) genome of the dolphin,Cephalorhynchus commersonii. Hybridization experiments with purified primate mtDNA probes have established that there is close resemblance in the general organization of the dolphin mt genome and the terrestrial mammalian mt genomes. Sequences covering 2381 bp of the dolphin mt genome from the major noncoding region, three tRNA genes, and parts of the genes encoding cytochrome b, NADH dehydrogenase subunit 3 (ND3), and 16S rRNA have been compared with corresponding regions from other mammalian genomes. There is a general tendency throughout the sequenced regions for greater similarity between dolphin and bovine mt genomes than between dolphin and rodent or human mt genomes.     

Mitochondrial Genome of <Emphasis Type="Italic">Ciona savignyi</Emphasis> (Urochordata,Ascidiacea, Enterogona): Comparison of Gene Arrangement and tRNA Genes with <Emphasis Type="Italic">Halocynthia roretzi</Emphasis> Mitochondrial Genome

The complete nucleotide sequence of the urochordate Ciona savignyi (Ascidiacea, Enterogona) mitochondrial (mt) genome (14,737 bp) was determined. The Ciona mt genome does not encode a gene for ATP synthetase subunit 8 but encodes an additional tRNA^Gly gene (anticodon UCU), as is the case in another urochordate, Halocynthia roretzi (Ascidiacea, Pleurogona), mt genome. In addition, the Ciona mt genome encodes two tRNA^Met genes; anticodon CAT and anticodon TAT. The tRNA^Cys gene is thought to lack base pairs at the D-stem. Thus, the Ciona mt genome encodes 12 protein, 2 rRNA, and 24 tRNA genes. The gene arrangement of the Ciona mt genome differs greatly from those of any other metazoan mt genomes reported to date. Only three gene boundaries are shared between the Halocynthia and the Ciona mt genomes. Molecular phylogenetic analyses based on amino acid sequences of mt protein genes failed to demonstrate the monophyly of the chordates.     

The Complete Mitochondrial Genome of Gossypium hirsutum and Evolutionary Analysis of Higher Plant Mitochondrial Genomes

Background

Mitochondria are the main manufacturers of cellular ATP in eukaryotes. The plant mitochondrial genome contains large number of foreign DNA and repeated sequences undergone frequently intramolecular recombination. Upland Cotton (Gossypium hirsutum L.) is one of the main natural fiber crops and also an important oil-producing plant in the world. Sequencing of the cotton mitochondrial (mt) genome could be helpful for the evolution research of plant mt genomes.

Methodology/Principal Findings

We utilized 454 technology for sequencing and combined with Fosmid library of the Gossypium hirsutum mt genome screening and positive clones sequencing and conducted a series of evolutionary analysis on Cycas taitungensis and 24 angiosperms mt genomes. After data assembling and contigs joining, the complete mitochondrial genome sequence of G. hirsutum was obtained. The completed G.hirsutum mt genome is 621,884 bp in length, and contained 68 genes, including 35 protein genes, four rRNA genes and 29 tRNA genes. Five gene clusters are found conserved in all plant mt genomes; one and four clusters are specifically conserved in monocots and dicots, respectively. Homologous sequences are distributed along the plant mt genomes and species closely related share the most homologous sequences. For species that have both mt and chloroplast genome sequences available, we checked the location of cp-like migration and found several fragments closely linked with mitochondrial genes.

Conclusion

The G. hirsutum mt genome possesses most of the common characters of higher plant mt genomes. The existence of syntenic gene clusters, as well as the conservation of some intergenic sequences and genic content among the plant mt genomes suggest that evolution of mt genomes is consistent with plant taxonomy but independent among different species.     

Complete nucleotide sequence and gene arrangement of the mitochondrial genome of the crab-eating frog Fejervarya cancrivora and evolutionary implications

The complete nucleotide sequence of the mitochondrial (mt) genome of the crab-eating frog, Fejervarya cancrivora Gravenhorst (Amphibia: Anura: Ranidae), was determined. The mt genome is 17,843 bp long and contains 13 protein-coding (ATP6, ATP8, COI-III, ND1-6 and 4L, and Cyt b) and two ribosomal RNA (12S and 16SrRNA) genes. Although metazoan mt genomes typically encode 22 transfer RNA genes (tRNAs), the F. cancrivora mtDNA contains 23 tRNAs due to the presence of an extra copy of tRNA(Met). A major noncoding region and a prominent intergenic spacer corresponding to the control region and light-strand replication origin were also found. To confirm the phylogenetic position of F. cancrivora, we compared the gene arrangement with that of other anurans and performed phylogenetic analyses based on mt genomic data. The genome organization of F. cancrivora mtDNA differs from that of typical vertebrates and neobatrachian frogs but is identical with that of F. limnocharis, suggesting that the unique gene arrangement occurred in the common ancestor of the genus. Phylogenetic analyses supported the monophyly of the Fejervarya species used here as well as the dicroglossini clade. Although the family Ranidae as previously recognized (= Ranidae, Discoglossidae, and some other natatanuran families; sensu Frost et al., 2006) is shown as a clade in the maximum parsimony analysis, the maximum likelihood and the Bayesian analyses suggest the paraphyly of the Ranidae with respect to the families, Mantellidae and Rhacophoridae. Three-tandem duplications of gene regions followed by subsequent deletions of supernumerary genes were proposed to explain the evolution of the extra tRNA(Met) and translocation of ND5 from the original neobatrachian gene order.     

Mitochondrial genomes of Vanhornia eucnemidarum (Apocrita: Vanhorniidae) and Primeuchroeus spp. (Aculeata: Chrysididae): Evidence of rearranged mitochondrial genomes within the Apocrita (Insecta: Hymenoptera).

We sequenced most of the mitochondrial (mt) genomes of 2 apocritan taxa: Vanhornia eucnemidarum and Primeuchroeus spp. These mt genomes have similar nucleotide composition and codon usage to those of mt genomes reported for other Hymenoptera, with a total A + T content of 80.1% and 78.2%, respectively. Gene content corresponds to that of other metazoan mt genomes, but gene organization is not conserved. There are a total of 6 tRNA genes rearranged in V. eucnemidarum and 9 in Primeuchroeus spp. Additionally, several noncoding regions were found in the mt genome of V. eucnemidarum, as well as evidence of a sustained gene duplication involving 3 tRNA genes. We also report an inversion of the large and small ribosomal RNA genes in Primeuchroeus spp. mt genome. However, none of the rearrangements reported are phylogenetically informative with respect to the current taxon sample.     

The complete sequence of the mitochondrial genome of Buteo buteo (Aves, Accipitridae) indicates an early split in the phylogeny of raptors

The complete sequence of the mitochondrial (mt) genome of Buteo buteo was determined. Its gene content and nucleotide composition are typical for avian genomes. Due to expanded noncoding sequences, Buteo possesses the longest mt genome sequenced so far (18,674 bp). The gene order comprising the control region and neighboring genes is identical to that of Falco peregrinus, suggesting that the corresponding rearrangement occurred before the falconid/accipitrid split. Phylogenetic analyses performed with the mt sequence of Buteo and nine other mt genomes suggest that for investigations at higher taxonomic levels (e.g., avian orders), concatenated rRNA and tRNA gene sequences are more informative than protein gene sequences with respect to resolution and bootstrap support. Phylogenetic analyses indicate an early split between Accipitridae and Falconidae, which, according to molecular dating of other avian divergence times, can be assumed to have taken place in the late Cretaceous 65-83 MYA.     

Cloning of the mitochondrial genome of Rana catesbeiana and the nucleotide sequences of the ND2 and five tRNA genes

The entire mitochondrial genome of Rana catesbeiana was cloned into a plasmid vector pBR322 at the unique BamHI site and the nucleotide sequences of the ND2 gene and of its flanking genes were determined. The ND2 gene was encoded by 1,033 base pairs and, as deduced from the nucleotide sequence, the ND2 product consisted of 344 amino acids with a molecular weight of 37,561. This gene was flanked on the 5' side by the tRNA genes for isoleucine, glutamine, and methionine and on the 3' side by those for tryptophan and alanine. These genes were the same in their organization as those found in the mammalian and Xenopus laevis mitochondrial genomes. A comparison of the putative amino acid sequences of the ND2 proteins of different animal species revealed that six regions in the sequence were well conserved during evolution, suggesting that some of these conserved sequences are crucial for biological activity of the ND2 protein. The nucleotide sequence homologies between the five tRNA genes of R. catesbeiana and their counterparts of mammals and X. laevis were in the range of 55 to 85%, depending on the tRNA and animal species.