Molecular mechanisms of extensive mitochondrial gene rearrangement in plethodontid salamanders

Extensive gene rearrangement is reported in the mitochondrial genomes of lungless salamanders (Plethodontidae). In each genome with a novel gene order, there is evidence that the rearrangement was mediated by duplication of part of the mitochondrial genome, including the presence of both pseudogenes and additional, presumably functional, copies of duplicated genes. All rearrangement-mediating duplications include either the origin of light-strand replication and the nearby tRNA genes or the regions flanking the origin of heavy-strand replication. The latter regions comprise nad6, trnE, cob, trnT, an intergenic spacer between trnT and trnP and, in some genomes, trnP, the control region, trnF, rrnS, trnV, rrnL, trnL1, and nad1. In some cases, two copies of duplicated genes, presumptive regulatory regions, and/or sequences with no assignable function have been retained in the genome following the initial duplication; in other genomes, only one of the duplicated copies has been retained. Both tandem and nontandem duplications are present in these genomes, suggesting different duplication mechanisms. In some of these mitochondrial DNAs, up to 25% of the total length is composed of tandem duplications of noncoding sequence that includes putative regulatory regions and/or pseudogenes of tRNAs and protein-coding genes along with the otherwise unassignable sequences. These data indicate that imprecise initiation and termination of replication, slipped-strand mispairing, and intramolecular recombination may all have played a role in generating repeats during the evolutionary history of plethodontid mitochondrial genomes.     

Complete sequence of the mitochondrial genome of Taenia saginata: comparison with T. solium and T. asiatica

The complete sequence of the Taenia saginata mitochondrial genome was determined, and its organization and structure were compared to other human-tropic Taenia tapeworms for which complete mitochondrial sequence data were available. The mitochondrial genome was 13,670 bp long, contained 12 protein-coding genes, two ribosomal RNAs (rRNAs, a small and a large subunit), and 22 transfer RNAs (tRNAs). It did not encode the atp8 gene. Overlapping regions were found between nad4L and nad4, nad1 and trnN, and cox1 and trnT. The ATG initiation codon was used for 10 protein-coding genes, and the GTG initiation codon was used for the remaining 2 genes (nad4 and atp6). The size of the protein-coding genes of the three human Taenia tapeworms did not vary, except for Taenia solium nad1 (891 aa) and nad4 (1212 aa) and Taenia asiatica cox2 (576 aa). The tRNA genes were 57-75 bp long, and the predicted secondary structures of 18 of these genes had typical clover-leaf shapes with paired dihydrouridine (DHU) arms. The genes in all human Taenia tapeworms for the two mitochondrial rRNA subunits rrnL and rrnS are separated by trnC. The putative T. saginata rrnL and rrnS are 972 and 732 bp long, respectively. The non-coding regions of the mt genome of T. saginata consisted of 2 regions: a short non-coding region (SNR, 66 nucleotides) and a long non-coding region (LNR, 159 nucleotides). The overall sequence difference in the full mitochondrial genome between T. saginata and T. asiatica was 4.6%, while T. solium differed by 11%. In conclusion, the complete sequence of the T. saginata mitochondrial genome will serve as a resource for comparative mitochondrial genomics and systematic studies of the parasitic cestodes.     

Complex repeat structures and novel features in the mitochondrial genomes of the diatoms Phaeodactylum tricornutum and Thalassiosira pseudonana

The mitochondrial genome of the raphid pennate diatom Phaeodactylum tricornutum has several novel features compared with the mitochondrial genomes of the centric diatom Thalassiosira pseudonana and the araphid pennate diatom Synedra acus. It is almost double the size (77,356 bp) due to a 35,454 bp sequence block consisting of an elaborate combination of direct repeats, making it the largest stramenopile (heterokont) mitochondrial genome known. In addition, the cox1 gene has a +1 translational frameshift involving Pro codons CCC and CCT, the first translational frameshift to be detected in an algal mitochondrial genome. The nad9 and rps14 genes are fused by the insertion of an in-frame sequence and cotranscribed. The nad11 gene is split into two parts corresponding to the FeS and molybdate-binding domains, but both parts are still on the mitochondrial genome, in contrast to the brown algae where the second domain appears to have been transferred to the nucleus. In contrast to P. tricornutum, the repeat region of T. pseudonana consists of a much smaller 4790 bp string of almost identical double-hairpin elements, evidence of slipped-strand mispairing and active gene conversion. The diatom mitochondrial genomes have undergone considerable gene rearrangement since the three lineages of diatoms diverged, but all three have kept their repeat regions segregated from their relatively compact coding regions.     

Complete Mitochondrial Genome of Haplorchis taichui and Comparative Analysis with Other Trematodes

Mitochondrial genomes have been extensively studied for phylogenetic purposes and to investigate intra- and interspecific genetic variations. In recent years, numerous groups have undertaken sequencing of platyhelminth mitochondrial genomes. Haplorchis taichui (family Heterophyidae) is a trematode that infects humans and animals mainly in Asia, including the Mekong River basin. We sequenced and determined the organization of the complete mitochondrial genome of H. taichui. The mitochondrial genome is 15,130 bp long, containing 12 protein-coding genes, 2 ribosomal RNAs (rRNAs, a small and a large subunit), and 22 transfer RNAs (tRNAs). Like other trematodes, it does not encode the atp8 gene. All genes are transcribed from the same strand. The ATG initiation codon is used for 9 protein-coding genes, and GTG for the remaining 3 (nad1, nad4, and nad5). The mitochondrial genome of H. taichui has a single long non-coding region between trnE and trnG. H. taichui has evolved as being more closely related to Opisthorchiidae than other trematode groups with maximal support in the phylogenetic analysis. Our results could provide a resource for the comparative mitochondrial genome analysis of trematodes, and may yield genetic markers for molecular epidemiological investigations into intestinal flukes.     

The complete mitochondrial genome of the grand jackknife clam, Solen grandis (Bivalvia: Solenidae): a novel gene order and unusual non-coding region

Molluscs in general, and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. The complete mitochondrial genome of Solen grandis (Bivalvia: Solenidae) was determined using long-PCR and genome walking techniques. The entire mitochondrial genome sequence of S. grandis is 16,784 bp in length, and contains 36 genes including 12 protein-coding genes (atp8 is absent), 2 ribosomal RNAs, and 22 tRNAs. All genes are encoded on the same strand. Compared with other species, it bears a novel gene order. Besides these, we find a peculiar non-coding region of 435 bp with a microsatellite-like (TA)₁₂ element, poly-structures and many hairpin structures. In contrast to the available heterodont mitochondrial genomes from GenBank, the complete mtDNA of S. grandis has the shortest cox3 gene, and the longest atp6, nad4, nad5 genes.     

The mitochondrial genome of Polylabris halichoeres (Monogenea: Microcotylidae)

This study presents the complete mitochondrial (mt) genome of Polylabris halichoeres, which is the largest mt genome sequenced of monogeneans so far and the second complete sequence after Microcotyle sebastis from the Microcotylidae. It is basically similar to that of M. sebastis, with the exception of a high level of gene rearrangement located between trnC and trnL((UUR)), a translocation of trnM and trnH, as well as a highly repetitive region (HRR) in the large non-coding region (NCR). We also find a series of trnI pseudogenes (ΨI) and one unknown short open reading frame (ORF) in the large NCR. Although the ORF cannot be unambiguously regarded as an atp8 gene, we cannot rule out the possibility that it has other functional importance, but it need further study in the future.     

Multiple rearrangements in mitochondrial genomes of Isopoda and phylogenetic implications

In this study, we analyse the evolutionary dynamics and phylogenetic implications of gene order rearrangements in five newly sequenced mitochondrial (mt) genomes and four published mt genomes of isopod crustaceans. The sequence coverage is nearly complete for four of the five newly sequenced species, with only the control region and some tRNA genes missing, while in Janira maculosa only two thirds of the genome could be determined. Mitochondrial gene order in isopods seems to be more plastic than that in other crustacean lineages, making all nine known mt gene orders different. Especially the asellote Janira is characterized by many autapomorphies. The following inferred ancestral isopod mt gene order exists slightly modified in modern isopods: nad1, tnrL1, rrnS, control region, trnS1, cob, trnT, nad5, trnF. We consider the inferred gene translocation events leading to gene rearrangements as valuable characters in phylogenetic analyses. In this first study covering major isopod lineages, potential apomorphies were identified, e.g., a shared relative position of trnR in Valvifera. We also report one of the first findings of homoplasy in mitochondrial gene order, namely a shared relative position of trnV in unrelated isopod lineages. In addition to increased taxon sampling secondary structure, modification in tRNAs and GC-skew inversion may be potentially fruitful subjects for future mt genome studies in a phylogenetic context.     

The complete mitochondrial genome of Evania appendigaster (Hymenoptera: Evaniidae) has low A+T content and a long intergenic spacer between atp8 and atp6

The apocritan Hymenoptera show extraordinary features in mitochondrial genomes, but no complete sequence has been reported for the basal lineage, Evanioidea. Here, we sequenced the complete mitochondrial genome of Evania appendigaster. This genome is 17,817 bp long; with low A+T content, 77.8%, compared with other hymenopteran species. Four tRNA genes were rearranged, among which remote inversion is the dominant gene rearrangement event. Gene shuffling is caused by tandem duplication-random loss while remote inversion is best explained by recombination. The start codon of nad1 was found as TTG, which might be common across Hymenoptera. trnS2 and trnK use abnormal anticodons TCT and TTT, respectively, and the D-stem pairings in trnS2 are absent. The secondary structure of two rRNA genes are predicted and compared with those in other insects. Five long intergenic spacers were present, including a long intergenic spacer between atp8 and atp6, where these two genes overlap in the previously reported animal genomes. A conserved motif was found between trnS1 and nad1, which is proposed to be associated with mtTERM. The A+T-rich region is 2,325 bp long, among the longest in insects, and contains a tandem repeat region.     

Rapid evolution of the compact and unusual mitochondrial genome in the ctenophore, Pleurobrachia bachei

Ctenophores are one of the most basally branching lineages of metazoans with the largest mitochondrial organelles in the animal kingdom. We sequenced the mitochondrial (mtDNA) genome from the Pacific cidipid ctenophore, Pleurobrachia bachei. The circular mitochondrial genome is 11,016 nts, with only 12 genes, and one of the smallest metazoan mtDNA genomes recorded. The protein coding genes are intronless cox1-3, cob, nad1, 3, 4, 4L and 5. The nad2 and 6 genes are represented as short fragments whereas the atp6 gene was found in the nuclear genome. Only the large ribosomal RNA subunit and two tRNAs were present with possibly the small subunit unidentifiable due to extensive fragmentation. The observed unique features of this mitochondrial genome suggest that nuclear and mitochondrial genomes have evolved at very different rates. This reduced mtDNA genome sharply contrasts with the very large sizes of mtDNA found in other basal metazoans including Porifera (sponges), and Placozoa (Trichoplax).     

Comparative analysis of mitochondrial genome data for Necator americanus from two endemic regions reveals substantial genetic variation

Necator americanus is a blood-sucking, intestinal nematode of major human health importance in many tropical and subtropical regions of the world. The aim of the present study was to compare the complete mitochondrial genome sequence from one N. americanus individual from Togo with another from China, in order to estimate the magnitude of genetic variability for different mitochondrial genes and non-coding regions. For the 12 protein genes, this comparison revealed sequence differences at both the nucleotide (3-7%) and amino acid (1-7%) levels. The most conserved of these was the nad4L gene, whereas the nad1 gene was least conserved at both the nucleotide and amino acid levels. Nucleotide differences were also detected in 14 of the 22 transfer RNAs (trns) (1-13%), the AT-rich region ( approximately 8%), non-coding regions (8-25%) and in the small (rrnS) and large (rrnL) subunits of mitochondrial ribosomal RNA (rrn) ( approximately 1%). Comparison of the rrnL sequences among multiple individual worms revealed nine unequivocal nucleotide differences between N. americanus from the two countries. Consistent with previous studies, these findings provide evidence for substantial genetic variation within N. americanus, which may have implications for the transmission and control of hookworm disease.     

The complete mitochondrial genome of the black coral Chrysopathes formosa (Cnidaria:Anthozoa:Antipatharia) supports classification of antipatharians within the subclass Hexacorallia

Black corals comprise a globally distributed shallow- and deep-water taxon whose phylogenetic position within the Anthozoa has been debated. We sequenced the complete mitochondrial genome of the antipatharian Chrysopathes formosa to further evaluate its phylogenetic relationships. The circular mitochondrial genome (18,398 bp) consists of 13 energy pathway protein-coding genes and two ribosomal RNAs, but only two transfer RNA genes (trnM and trnW), as well as a group I intron within the nad5 gene that contains the only copies of nad1 and nad3. No novel genes were found in the antipatharian mitochondrial genome. Gene order and genome content are most similar to those of the sea anemone Metridium senile (subclass Hexacorallia), with differences being the relative location of three contiguous genes (cox2-nad4-nad6) and absence (from the antipatharian) of a group I intron within the cox1 gene. Phylogenetic analyses of multiple protein-coding genes support classifying the Antipatharia within the subclass Hexacorallia and not the subclass Ceriantipatharia; however, the sister-taxon relationships of black corals within Hexacorallia remain inconclusive.     

The complete mitochondrial genome of the nudibranch Roboastra europaea (Mollusca: Gastropoda) supports the monophyly of opisthobranchs

The complete nucleotide sequence (14,472 bp) of the mitochondrial genome of the nudibranch Roboastra europaea (Gastropoda: Opisthobranchia) was determined. This highly compact mitochondrial genome is nearly identical in gene organization to that found in opisthobranchs and pulmonates (Euthyneura) but not to that in prosobranchs (a paraphyletic group including the most basal lineages of gastropods). The newly determined mitochondrial genome differs only in the relative position of the trnC gene when compared with the mitochondrial genome of Pupa strigosa, the only opisthobranch mitochondrial genome sequenced so far. Pupa and Roboastra represent the most basal and derived lineages of opisthobranchs, respectively, and their mitochondrial genomes are more similar in sequence when compared with those of pulmonates. All phylogenetic analyses (maximum parsimony, minimum evolution, maximum likelihood, and Bayesian) based on the deduced amino acid sequences of all mitochondrial protein-coding genes supported the monophyly of opisthobranchs. These results are in agreement with the classical view that recognizes Opisthobranchia as a natural group and contradict recent phylogenetic studies of the group based on shorter sequence data sets. The monophyly of opisthobranchs was further confirmed when a fragment of 2,500 nucleotides including the mitochondrial cox1, rrnL, nad6, and nad5 genes was analyzed in several species representing five different orders of opisthobranchs with all common methods of phylogenetic inference. Within opisthobranchs, the polyphyly of cephalaspideans and the monophyly of nudibranchs were recovered. The evolution of mitochondrial tRNA rearrangements was analyzed using the cox1+rrnL+nad6+nad5 gene phylogeny. The relative position of the trnP gene between the trnA and nad6 genes was found to be a synapomorphy of opisthobranchs that supports their monophyly.     

The mitochondrial genome of the thermal dimorphic fungus Penicillium marneffei is more closely related to those of molds than yeasts

We report the complete sequence of the mitochondrial genome of Penicillium marneffei, the first complete mitochondrial DNA sequence of a thermal dimorphic fungus. This 35 kb mitochondrial genome contains the genes encoding ATP synthase subunits 6, 8, and 9 (atp6, atp8, and atp9), cytochrome oxidase subunits I, II, and III (cox1, cox2, and cox3), apocytochrome b (cob), reduced nicotinamide adenine dinucleotide ubiquinone oxireductase subunits (nad1, nad2, nad3, nad4, nad4L, nad5, and nad6), ribosomal protein of the small ribosomal subunit (rps), 28 tRNAs, and small and large ribosomal RNAs. Analysis of gene contents, gene orders, and gene sequences revealed that the mitochondrial genome of P. marneffei is more closely related to those of molds than yeasts.     

Complete mitochondrial genomes of two Japanese precious corals, Paracorallium japonicum and Corallium konojoi (Cnidaria, Octocorallia, Coralliidae): notable differences in gene arrangement

Precious coral are taxonomically a group of corals that belong to the family Coralliidae within the order Alcyonacea, subclass Octocorallia, and class Anthozoa, whose skeletal axes are used for jewelry. They are distributed in the Mediterranean Sea and in waters adjacent to Japan, Taiwan, Midway Island and the Hawaiian Islands. The genus Corallium of the family Coralliidae was recently divided into two genera, Corallium and Paracorallium, based on morphological observations, but insufficient molecular evidence to support this classification has been presented to date. We determined for the first time the complete mitochondrial genome sequence of two precious corals P. japonicum and C. konojoi, in order to clarify their systematic positions. The circular mitochondrial genomes of P. japonicum and C. konojoi are 18,913bp and 18,969bp in length, respectively, and encode 13 typical energy pathway protein coding genes (nad1-6, nad4L, cox1-3, cob, atp6 and atp8), two ribosomal RNA genes (rns and rnl), a transfer RNA (trnM) and a mismatch repair gene homologue msh1. The two genomes have an overall nucleotide sequence identity of 97.5%, which is comparable to that between Acanella eburnea and Keratoisidinae sp. belonging to Octocorallia. Surprisingly, however, their gene arrangements were not identical. Phylogenetic analyses using seven complete mitochondrial genome sequences belonging to species in the subclass Octocorallia indicated that within the subclass, at least three gene order rearrangement events occurred during evolution. Our results support the validity of the morphological classification that separated the family Coralliidae into two genera, Corallium and Paracorallium.     

The analysis of the complete mitochondrial genome of Lecanicillium muscarium (synonym Verticillium lecanii) suggests a minimum common gene organization in mtDNAs of Sordariomycetes: phylogenetic implications

The mitochondrial genome (mtDNA) of the entomopathogenic fungus Lecanicillium muscarium (synonym Verticillium lecanii) with a total size of 24,499-bp has been analyzed. So far, it is the smallest known mitochondrial genome among Pezizomycotina, with an extremely compact gene organization and only one group-I intron in its large ribosomal RNA (rnl) gene. It contains the 14 typical genes coding for proteins related to oxidative phosphorylation, the two rRNA genes, one intronic ORF coding for a possible ribosomal protein (rps), and a set of 25 tRNA genes which recognize codons for all amino acids, except alanine and cysteine. All genes are transcribed from the same DNA strand. Gene order comparison with all available complete fungal mtDNAs-representatives of all four Phyla are included-revealed some characteristic common features like uninterrupted gene pairs, overlapping genes, and extremely variable intergenic regions, that can all be exploited for the study of fungal mitochondrial genomes. Moreover, a minimum common mtDNA gene order could be detected, in two units, for all known Sordariomycetes namely nad1-nad4-atp8-atp6 and rns-cox3-rnl, which can be extended in Hypocreales, to nad4L-nad5-cob-cox1-nad1-nad4-atp8-atp6 and rns-cox3-rnl nad2-nad3, respectively. Phylogenetic analysis of all fungal mtDNA essential protein-coding genes as one unit, clearly demonstrated the superiority of small genome (mtDNA) over single gene comparisons.     

Complete sequences of mitochondria genomes of Aedes aegypti and Culex quinquefasciatus and comparative analysis of mitochondrial DNA fragments inserted in the nuclear genomes

We present complete sequences of the mitochondrial genomes for two important mosquitoes, Aedes aegypti and Culex quinquefasciatus, that are major vectors of dengue virus and lymphatic filariasis, respectively. The A. aegypti mitochondrial genome is 16,655 bp in length and that of C. quinquefasciatus is 15,587 bp, yet both contain 13 protein coding genes, 22 transfer RNA (tRNA) genes, one 12S ribosomal RNA (rRNA) gene, one 16S rRNA gene and a control region (CR) in the same order. The difference in the genome size is due to the difference in the length of the control region. We also analyzed insertions of nuclear copies of mtDNA-like sequences (NUMTs) in a comparative manner between the two mosquitoes. The NUMT sequences occupy ~0.008% of the A. aegypti genome and ~0.001% of the C. quinquefasciatus genome. Several NUMTs were found localized in the introns of predicted protein coding genes in both genomes (32 genes in A. aegypti but only four in C. quinquefasciatus). None of these NUMT-containing genes had an ortholog between the two species or had paralogous copies within a genome that was also NUMT-containing. It was further observed that the NUMT-containing genes were relatively longer but had lower GC content compared to the NUMT-less paralogous copies. Moreover, stretches of homologies are present among the genic and non-genic NUMTs that may play important roles in genomic rearrangement of NUMTs in these genomes. Our study provides new insights on understanding the roles of nuclear mtDNA sequences in genome complexities of these mosquitoes.     

The mitochondrial genomes of the human hookworms, Ancylostoma duodenale and Necator americanus (Nematoda: Secernentea).

The complete mitochondrial genome sequences were determined for two species of human hookworms, Ancylostoma duodenale (13,721 bp) and Necator americanus (13,604 bp). The circular hookworm genomes are amongst the smallest reported to date for any metazoan organism. Their relatively small size relates mainly to a reduced length in the AT-rich region. Both hookworm genomes encode 12 protein, two ribosomal RNA and 22 transfer RNA genes, but lack the ATP synthetase subunit 8 gene, which is consistent with three other species of Secernentea studied to date. All genes are transcribed in the same direction and have a nucleotide composition high in A and T, but low in G and C. The AT bias had a significant effect on both the codon usage pattern and amino acid composition of proteins. For both hookworm species, genes were arranged in the same order as for Caenorhabditis elegans, except for the presence of a non-coding region between genes nad3 and nad5. In A. duodenale, this non-coding region is predicted to form a stem-and-loop structure which is not present in N. americanus. The mitochondrial genome structure for both hookworms differs from Ascaris suum only in the location of the AT-rich region, whereas there are substantial differences when compared with Onchocerca volvulus, including four gene or gene-block translocations and the positions of some transfer RNA genes and the AT-rich region. Based on genome organisation and amino acid sequence identity, A. duodenale and N. americanus were more closely related to C. elegans than to A. suum or O. volvulus (all secernentean nematodes), consistent with a previous phylogenetic study using ribosomal DNA sequence data. Determination of the complete mitochondrial genome sequences for two human hookworms (the first members of the order Strongylida ever sequenced) provides a foundation for studying the systematics, population genetics and ecology of these and other nematodes of socio-economic importance.     

Extensive variation in nuclear mitochondrial DNA content between the genomes of Phytophthora sojae and Phytophthora ramorum

Fragments of mitochondrial DNA (mtDNA) transferred to the nuclear genome are called nuclear mitochondrial DNAs (NUMTs). We report here a comparison of NUMT content between genomes from two species of the same genus. Analysis of the genomes of Phytophthora sojae and P. ramorum revealed large differences in the NUMT content of the two genomes: 16.27 x 10(-3) and 2.28 x 10(-3)% of each genome, respectively. Substantial differences also exist between the two species in the sizes of the NUMTs found in each genome, with ranges of 20 to 405 bp for P. sojae and 19 to 137 bp for P. ramorum. Furthermore, in P. sojae, fragments from the mitochondrial genes rns, rnl, coxl, and nad (various subunits) are found most frequently, whereas P. ramorum NUMTs most often originate from the cox3, rpsl4, nad4, and nad5 genes. The large differences in the presumptive mtDNA insertions suggest that the insertions occurred subsequent to the divergence of the two species, and this is supported by sequence comparisons among the NUMTs and the mtDNA sequences of the two species. P. sojae mtDNA sequences inserted in the nuclear genome appear to have been altered as a result of insertions, deletions, inversions, and translocations and provide insights into active mechanisms of sequence divergence in this plant pathogen. No clear examples were found of NUMTs forming functional nuclear genes or of NUMTs inserted into exons or introns of any nuclear gene.     

Making the most of mitochondrial genomes--markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea)

An increasing number of complete sequences of mitochondrial (mt) genomes provides the opportunity to optimise the choice of molecular markers for phylogenetic and ecological studies. This is particularly the case where mt genomes from closely related taxa have been sequenced; e.g., within Schistosoma. These blood flukes include species that are the causative agents of schistosomiasis, where there has been a need to optimise markers for species and strain recognition. For many phylogenetic and population genetic studies, the choice of nucleotide sequences depends primarily on suitable PCR primers. Complete mt genomes allow individual gene or other mt markers to be assessed relative to one another for potential information content, prior to broad-scale sampling. We assess the phylogenetic utility of individual genes and identify regions that contain the greatest interspecific variation for molecular ecological and diagnostic markers. We show that variable characters are not randomly distributed along the genome and there is a positive correlation between polymorphism and divergence. The mt genomes of African and Asian schistosomes were compared with the available intraspecific dataset of Schistosoma mansoni through sliding window analyses, in order to assess whether the observed polymorphism was at a level predicted from interspecific comparisons. We found a positive correlation except for the two genes (cox1 and nad1) adjoining the putative control region in S. mansoni. The genes nad1, nad4, nad5, cox1 and cox3 resolved phylogenies that were consistent with a benchmark phylogeny and in general, longer genes performed better in phylogenetic reconstruction. Considering the information content of entire mt genome sequences, partial cox1 would not be the ideal marker for either species identification (barcoding) or population studies with Schistosoma species. Instead, we suggest the use of cox3 and nad5 for both phylogenetic and population studies. Five primer pairs designed against Schistosoma mekongi and Schistosoma malayensis were tested successfully against Schistosoma japonicum. In combination, these fragments encompass 20-27% of the variation amongst the genomes (average total length approximately 14,000bp), thus providing an efficient means of encapsulating the greatest amount of variation within the shortest sequence. Comparative mitogenomics provides the basis of a rational approach to molecular marker selection and optimisation.