Short report | Open | Published:
The mitochondrial genome of the egg-laying flatworm Aglaiogyrodactylus forficulatus (Platyhelminthes: Monogenoidea)
Parasites & Vectorsvolume 9, Article number: 285 (2016)
The rather species-poor oviparous gyrodactylids are restricted to South America. It was suggested that they have a basal position within the otherwise viviparous Gyrodactylidae. Accordingly, it was proposed that the species-rich viviparous gyrodactylids diversified and dispersed from there.
The mitochondrial genome of Aglaiogyrodactylus forficulatus was bioinformatically assembled from next-generation illumina MiSeq sequencing reads, annotated, and compared to previously published mitochondrial genomes of other monogenoidean flatworm species.
The mitochondrial genome of A. forficulatus consists of 14,371 bp with an average A + T content of 75.12 %. All expected 12 protein coding, 22 tRNA, and 2 rRNA genes were identified. Furthermore, there were two repetitive non-coding regions essentially consisting of 88 bp and 233 bp repeats, respectively. Maximum Likelihood analyses placed the mitochondrial genome of A. forficulatus in a well-supported clade together with the viviparous Gyrodactylidae species. The gene order differs in comparison to that of other monogenoidean species, with rearrangements mainly affecting tRNA genes. In comparison to Paragyrodactylus variegatus, four gene order rearrangements, i.e. three transpositions and one complex tandem-duplication-random-loss event, were detected.
Mitochondrial genome sequence analyses support a basal position of the oviparous A. forficulatus within Gyrodactylidae, and a sister group relationship of the oviparous and viviparous forms.
Within flatworms (Platyhelminthes) the parasitic Neodermata represent the most derived forms and are, at least when compared to the free-living flatworm lineages, particularly species rich. They include the ectoparasitic Monogenoidea, and the endoparasitic flukes (Trematoda) and tapeworms (Cestoda).
Among the Monogenoidea, Gyrodactylidae have attracted particular attention as some species have been of great concern for humans. For example, Gyrodactylus salaris Malmberg, 1957 has caused significant ecological and economic damage to wild stocks of Atlantic salmon (Salmo salar) in Norway and Russia as well as in aquaculture . The majority of the species within Gyrodactylidae are hyperviviparous ectoparasites on actinopterygian fish hosts, but there are also some oviparous gyrodactylid lineages. These lineages were originally placed within the Oogyrodactylidae , which was later rejected as paraphyletic by Boeger et al. , and Oogyrodactylidae were included in Gyrodactylidae based on morphological synapomorphies.
Oviparous gyrodactylids are restricted to South America and occur on freshwater catfishes, mainly on species of the Loricariidae (Siluriformes). Boeger et al.  considered the oviparous species a basal lineage within Gyrodactylidae. Taking into account their geographical restriction to continental South America they argued that the species-rich viviparous gyrodactylids diversified and dispersed from there . Accordingly, hyperviviparity and the loss of ‘sticky’ eggs were interpreted as synapomorphic key innovations of the viviparous Gyrodactylidae.
In comparison to the enormous diversity of the viviparous gyrodactylids the oviparous lineages appear rather species-poor; until today only 23 species in seven genera have been described . Among them is Aglaiogyrodactylus forficulatus Kritsky, Vianna & Boeger, 2007, the type-species of the genus, that was originally described from the loricarid host Kronichthys lacerta .
In the current study we present the first mitochondrial genome of an oviparous gyrodactylid. The molecular characteristics of the mitochondrial genome of A. forficulatus are compared with those of other previously described monogenoidean mitochondrial genomes, and the phylogenetic position of oviparous gyrodactylids is addressed.
The oviparous A. forficulatus was collected from its type-host, the catfish Kronichthys lacerta Nichols, 1919, from Rio Morato, basin of the Garaqueçaba, Paraná, Brazil (25°12'48''S, 48°17'52''W) on 30 October 2013.
About 200 individual parasites were pooled for DNA extraction with the E.Z.N.A. Tissue DNA kit (Omega Bio-Tek) following the Tissue DNA-Spin Protocol provided with the kit. For high-throughput next-generation sequencing (NGS) of the genomic DNA paired-end libraries were prepared, tagged, and analyzed (29,107,020 paired 300 bp reads) on an illumina MiSeq (outsourced to GENterprise GENOMICS, Mainz, Germany).
The mitochondrial genome of A. forficulatus was reconstructed by assembling the NGS reads using MITObim 1.8  using essentially the default settings of the program, and, in addition, the program’s quality trimming option. The COII sequence of A. ctenistus (GenBank accession number KF751723; ) was used as seed sequence for the assembly.
Annotation of the mitochondrial genes was done using MITOS  and DOGMA . In addition, tRNA genes were also identified using tRNAscan-SE 1.21 . For most protein coding genes (PCGs) only the conserved domains were identified by the two annotation programs. Phylogenetic comparisons with already published monogenoidean mitochondrial genomes (Table 1) were, therefore, used to manually complete these genes’ annotation, which was assisted by blast searches of GenBank. Codon usage was analyzed using the Sequence Manipulation Suite v2 . The dot-plot approach of YASS  was used to identify and visualize the repeat regions. The mitochondrial gene map was drawn using SnapGene 3.0.
The assembled mitochondrial genome of A. forficulatus was compared to previously published mitochondrial genomes of other monogenoidean species [13–25] that are listed in Table 1. The mitochondrial genomes of the cestodes Echinococcus oligarthrus  and Hymenolepis diminuta  served as outgroups. All genes were aligned individually using the online version of MAFFT Alignment v7.245 . Subsequently the individual tRNA, rRNA, and PCG alignments were concatenated into one extended alignment, which did not include the non-coding regions of the mitochondrial genes. GBlock v0.91b  and GUIDANCE2  were applied to remove ambiguous and unreliable sections from the concatenated MAFFT alignments.
PartitionFinder v.1.1.1_Mac  was used to select the best-fit model of molecular evolution for the concatenated alignments using as recommended the program’s “raxml” search option and the Bayesian Information Criterion (BIC) for model selection. Maximum Likelihood analyses were carried out in PHYML 3.0  applying the GTR + I + G model and 1,000 bootstrap replicates for all analyses except those for the Cytb and ND3 genes that were performed applying the GTR + G model.
Comparisons of gene orders of the mitochondrial genomes of A. forficulatus and the other monogenoidean species included in this study were conducted in CREx , a program that infers most parsimonious gene rearrangements based on common intervals.
All sampling of parasites (on fish) was done under license number 10007 (Instituto Chico Mendes de Conservação da Biodiversidade - ICMBio, Brazil).
Results and discussion
Next-generation sequencing of the genomic DNA of A. forficulatus delivered 29,107,020 million paired 300 bp reads. The initial assembly with the COII sequence of A. ctenistus as seed sequence as well as subsequent control assemblies using various A. forficulatus seeds delivered identical mitochondrial DNA (mtDNA) consensus sequences. During the assembly of the mtDNA of A. forficulatus particular attention was paid to the two repetitive repeat regions (see below) in order to avoid erroneous base-calling due to misassembled short reads; the final sequence was determined only from those reads that also in part covered the respective flanking regions.
The mtDNA of A. forficulatus is 14,371 bp long, and was assembled from 43,334 quality-trimmed reads with an average coverage of 765 x. The overall A + T content was 75.1 % and the nucleotide composition was A (30.2 %) C (9.7 %), G (15.1 %), and T (44.9 %). The base composition of the tRNA and rRNA genes, PCGs, and repeat regions are listed in Additional file 1. The nucleotide sequence of the assembled mtDNA of A. forficulatus was deposited in GenBank with the accession number KU679421.
All genes usually described for the mitochondrial genomes of other Monogenoidea (12 protein coding genes, two rRNAs, and 22 tRNAs) were also identified for A. forficulatus (Fig. 1, Table 2). Accordingly, ATP8 lacks also in A. forficulatus. All genes are coded on the same strand.
Protein coding genes (PCGs)
For most PCGs the applied annotation programs only detected the conserved domains. This worked best for the COI and Cytb genes, whereas ND2 was not detected at all. Thus, most PCGs were annotated manually through phylogenetic comparisons with reference to other previously published platyhelminth mitochondrial genomes (Table 1). Translation into amino acid sequences was straightforward using the flatworm mitochondrial code . All PCG except for ND2 and ND4 use ATG as start codon. For ND2 and ND4, however, it proved difficult to identify the beginning of the gene, as there was no canonical start codon of flatworm mitochondria (reviewed in ) in frame in the respective region. As a working hypothesis, we suggest that both ND2 and ND4 start with ATT, which was reported earlier as a rarely used alternative start codon in other flatworm mitochondria (e.g. [13, 34]). All PCGs terminate with one of the canonical flatworm stop codons TAG or TAA (Table 2). As reported for other monogenoidean mitochondrial genomes, A. forficulatus shows a strong codon bias for PCGs. The most commonly used codons for a particular amino acid in A. forficulatus were, however, in all cases the same as for P. variegatus (Additional file 2).
The 12S and 16S rRNA genes were identified of being 729 bp and 928 bp long, respectively. The A + T – content was 77.4 % for the 12S gene and for the 75.6 % 16S rRNA gene, and thus very similar to the average of the whole mitochondrial genome.
All 22 tRNA genes were identified by MITOS  and DOGMA  that, however, failed to detect tRNA P for prolin. tRNAscan-SE  confirmed eight tRNAs but could not predict the remaining ones (Table 2). The secondary cloverleaf structures of the predicted tRNAs are compiled in Additional file 3. The tRNAs C, S1, and S2 lack the DHU arm.
Genetic diversity and phylogenetic analyses
Based on the concatenated MAFFT alignments (15,261 positions), the mitochondrial genome of the egg laying A. forficulatus was found in a sister-group relationship to the mitochondrial genomes of the other viviparous Gyrodactylidea, i.e. G. salaris, G. thymalli, G. derjavinoides, and P. variegatus (Fig. 2) with high statistical support. However, it has to be taken into account that the depicted ML tree does not provide a comprehensive phylogenetic hypothesis for the Monogenoidea since it is only based on the limited number of available mitochondrial genomes. All included Capsalidea, Mazocraeidea, respectively, also clustered together with high statistical support. Despite the limitations of the analyses the observed phylogenetic affinity of A. forficulatus to the included viviparous Gyrodactylidea species makes sense, as the grouping is congruent with the results of earlier morphological analyses .
Similar tree topologies with respect to the Gyrodactylidea clade were also obtained when analyzing the concatenated tRNA (bootstrap support 50) genes as well as the COI (100), COII (86), COIII (65), Cytb (77), ND1 (100), ND3 (66), ND4L (93), ND4 (99), and ND5 (100) genes individually. When analyzing the rRNA genes as well as the ATP6, ND2, and ND6 genes individually somewhat different tree topologies with affinities of A. forficulatus to other species or clades were obtained. Similar tree topologies were also obtained when using the concatenated GBlock and GUIDANCE2 alignments, respectively (Additional file 4). With the unreliable parts of the original MAFFT alignment removed, the clade consisting of species of Gyrodactylidae received a bootstrap support of 100.
The clustering of the mitochondrial genome of A. forficulatus with the remaining Gyrodactylidae species is nevertheless not strongly reflected in a matrix of pairwise genetic distances. Although the genetic distance between A. forficulatus and the closest relative P. variegatus was lowest (0.402; p-distances), the obtained values for all other pairwise comparisons were only slightly higher ranging from 0.409 with N. melleni to 0.477 with M. sebastis (Additional file 5). The respective values for the GBlock alignment were 0.396 for the comparison with P. variegatus and a range from 0.402 with N. melleni to 0.471 with M. sebastis (Additional file 5).
Non-coding regions (NCRs) including repeat regions I and II
Two extended non-coding regions of 486 bp and 733 bp, respectively, were detected in the mitochondrial genome of A. forficulatus. They are located relatively close to each other, separated only by the tRNA V and tRNA D genes, i.e. 123 bp, within a cluster of five tRNA genes that represents in comparison to other previously published monogenoidean mitochondrial genomes a substantially rearranged part of the mtDNA (Additional file 6). Both non-coding regions include a repetitive region consisting of three repeats. In repeat region I the repeats are 88 bp long while they span 233 bp in repeat region II. For both repetitive sequence motifs blast searches revealed no significant sequence similarity to other GenBank entries. The average sequence similarity between the repeats was 81.1 % in repeat region I and 96.5 % in repeat region II, with somewhat higher sequence divergence in the beginning and the end of the repeated arrays (Additional file 7). The AT-content of repeat region I is 75.48 %, and repeat region II is particularly AT rich (83.38 %) and may contain important functional domains for replication and transcription. However, such functional domains have not been described yet in monogenoidean flatworms in more detail than e.g. “Poly T stretch” or “A + T – rich segment” , and similar domains can be certainly found in the repeat region of A. forficulatus as well. Nevertheless, without a comprehensive functional analysis the identification of such functional domains would remain highly speculative.
The structure of the non-coding regions of A. forficulatus is strikingly similar to that of P. variegatus  although there is only one non-coding region in the latter species. The non-coding region in P. variegatus also consists essentially of two repetitive regions, one consisting of two 394 bp repeats (termed part I and II), and one consisting of three 81 bp repeats (termed part III). In contrast, the two non-coding regions observed in the mitochondrial genomes of the other included Gyroadactylidae species G. salaris , G. thymalli , and G. derjavinoides  do not consist of internal repeats. However, in these species the sequences of the two repeat regions are highly similar to each other indicating on the one hand that they originated from a duplication and on the other hand that they bear some functional domains.
As proposed by the CREx program  P. variegatus and T. nebulosi (both species share identical mitochondrial gene orders except for the repetitive non-coding region in P. variegatus) have the highest similarity values in gene order based on inferred common intervals (Additional file 8). The program suggested a recombination scenario consisting of three transpositions and one complex tandem-duplication-random-loss (tdrl) event. The proposed transpositions (inversions; green boxes in Fig. 3) affect (a) tRNA-F and tRNA-M, (b) tRNA-A and tRNA-D, and (c) tRNA-R and tRNA-L2, and the proposed tdrl affect two larger gene blocks (red and blue boxes in Fig. 3). However, the proposed tdrl event can alternatively also be interpreted as a series of transpositions. Four larger gene blocks are conserved over all Gyrodactylidae species included in this comparison (Additional file 6). Rearrangements of gene order have already been described for several mitochondrial genomes of Neodermata species (e.g. ). The four rearrangements in gene order of A. forficulatus in comparison to P. variegatus are thus not very surprising. As reported for other Monogenoidea mitochondrial gene order rearrangements mainly affect the specific position of some tRNA genes and non-coding regions. Also in A. forficulatus individual protein coding genes are not rearranged. For the mitochondrial genomes of the Gyrodactylidea species included in this study there are essentially four conserved regions; these are (i) ND5, tRNA G, COIII, tRNA H, Cytb, ND4L, ND4, (ii) ATP6, ND2, (iii) tRNA S1, tRNA W, COI, tRNA T, 16S, tRNA C, 12S, COII, tRNA E, ND6, tRNA Y, and (iv) ND1, tRNA N, tRNA P, tRNA I, tRNA K, ND3 (Additional file 6).
The mitochondrial genome of A. forficulatus shows a hitherto unique gene order within the monogenoiden Gyrodactylidea with four rearrangements in comparison to P. variegatus. The previously proposed sister group relationship of the oviparous and viviparous Gyrodactylidae is corroborated. However, more comprehensive sampling is required to further test the proposed phylogenetic hypothesis. All Gyrodactylidea mitochondrial genomes sequenced so far include repetitive regions, although the structure of two regions consisting of short tandemly arranged repeats that was found in the basal Gyrodactylidea lineages represented by A. forficulatus and P. variegatus differs substantially from the structure of two longer dispersed repeats in Gyrodactylus spp. The biological function of these repetitive regions is yet unknown but the sequencing of mitochondrial genomes of further Gyrodactylidae species may shed some light on the evolution of these regions.
The authors declare that they have no competing interests.
LB and WAB designed the study. LPA and WAB collected the parasite specimens. LB performed the experimental work. LB, BF, and WAB analyzed the data. All authors co-drafted the manuscript and approved the final version.
Harris PD, Bachmann L, Bakke TA. The Parasites and Pathogens of the Atlantic Salmon: Lessons from Gyrodactylus salaris. In: Aas Ø, Einum S, Klemetsen A, Skurdal J, editors. Atlantic salmon ecology. Chichester (United Kingdom): Blackwell Publishing; 2011. p. 221–52.
Harris PD. The morphology and life cycle of the oviparous Oogyrodactylus farlowellae gen. et sp. n. (Monogenea, Gyrodactylidea). Parasitology. 1983; doi:10.1017/S0031182000082937.
Boeger WA, Kritsky DC, Belmont-Jégu E. Neotropical Monogenoidea. 20. Two new species of oviparous Gyrodactylidea (Polyonchoinea) from loricariid catfishes (Siluriformes) in Brazil and the phylogenetic status of Ooegyrodactylidae Harris, 1983. J Helminthol Soc Wash. 1994;61:34–44.
Boeger WA, Kritsky DC, Pie MR. Context of diversification of the viviparous Gyrodactylidae (Platyhelminthes, Monogenoidea). Zool Scr. 2003; doi:10.1046/j.1463-6409.2003.00130.x.
Kritsky DC, Vianna RT, Boeger WA. Neotropical Monogenoidea. 50. Oviparous gyrodactylids from loricariid and pimelodid catfishes in Brazil, with the proposal of Phanerothecioides n. g., Onychogyrodactylus n. g. and Aglaiogyrodactylus n. g. (Polyonchoinea: Gyrodactylidea). Syst Parasitol. 2007; doi:10.1007/s11230-006-9053-7.
Hahn C, Bachmann L, Chevreux B. Reconstructing mitochondrial genomes directly from genomic next generation sequencing reads - a baiting and iterative mapping approach. Nucl Acids Res. 2013; doi:10.1093/nar/gkt371.
Bueno-Silva M, Boeger WA. Neotropical Monogenoidea. 58. Three new species of Gyrodactylus (Gyrodactylidae) from Scleromystax spp. (Callichthyidae) and the proposal of COII gene as an additional fragment for barcoding gyrodactylids. Folia Parasitol. 2014; doi:10.14411/fp.2014.028.
Bernt M, Donath A, Jühling F, Externbrink F, Florentz C, Fritzsch G, Pütz J, Middendorf M, Stadler PF. MITOS: Improved de novo Metazoan Mitochondrial Genome Annotation. Mol Phylogenet Evol. 2013; doi:10.1016/j.ympev.2012.08.023. http://mitos.bioinf.uni-leipzig.de/index.py.
Ye F, King SD, Cone DK, You P. The mitochondrial genome of Paragyrodactylus variegatus (Platyhelminthes: Monogenea): differences in major non-coding region and gene order to Gyrodactylus. Parasit Vectors. 2014; doi:10.1186/1756-3305-7-377.
Huyse T, Plaisance L, Webster BL, Mo TA, Bakke TA, Bachmann L, Littlewood DT. The mitochondrial genome of Gyrodactylus salaris (Platyhelminthes: Monogenea), a pathogen of Atlantic salmon (Salmo salar). Parasitology. 2007; doi:10.1017/S003118200600201.
Plaisance L, Huyse T, Littlewood DT, Bakke TA, Bachmann L. The complete mitochondrial DNA sequence of the monogenean Gyrodactylus thymalli (Platyhelminthes: Monogenea), a parasite of grayling (Thymallus thymallus). Mol Biochem Parasitol. 2007; doi:10.1016/j.molbiopara.2007.04.012.
Huyse T, Buchmann K, Littlewood DT. The mitochondrial genome of Gyrodactylus derjavinoides (Platyhelminthes: Monogenea) - a mitogenomic approach for Gyrodactylus species and strain identification. Gene. 2008; doi:10.1016/j.gene.2008.03.008.
Perkins EM, Donnellan SC, Bertozzi T, Whittington ID. Closing the mitochondrial circle on paraphyly of the Monogenea (Platyhelminthes) infers evolution in the diet of parasitic flatworms. Int J Parasitol. 2010; doi:10.1016/j.ijpara.2010.02.017.
Kang S, Kim J, Lee J, Kim S, Min GS Park JK. The complete mitochondrial genome of an ectoparasitic monopisthocotylean fluke Benedenia hoshinai (Monogenea: Platyhelminthes), Mitochondrial DNA. 2012; doi: 10.3109/19401736.2012.668900.
Zhang J, Wu X, Li Y, Zhao M, Xie M, Li A. The complete mitochondrial genome of Neobenedenia melleni (Platyhelminthes: Monogenea): mitochondrial gene content, arrangement and composition compared with two Benedenia species. Mol Biol Rep. 2014; doi: 10.1007/s11033-014-3542-6.
Zhang J, Wu X, Li Y, Xie M, Li A. The complete mitochondrial genome of Tetrancistrum nebulosi (Monogenea: Ancyrocephalidae). Mitochondrial DNA. 2016; doi:10.3109/19401736.2013.867441.
Zhang J, Wu X, Xie M, Xu X, Li A. The mitochondrial genome of Polylabris halichoeres (Monogenea: Microcotylidae). Mitochondrial DNA. 2011; doi:10.3109/19401736.2011.588223.
Park JK, Kim KH, Kang S, Kim W, Eom KS, Littlewood DT. A common origin of complex life cycles in parasitic flatworms: evidence from the complete mitochondrial genome of Microcotyle sebastis (Monogenea: Platyhelminthes). BMC Evol Biol. 2007; doi:10.1186/1471-2148-7-11.
Zhang J, Wu X, Xie M, Li A. The complete mitochondrial genome of Pseudochauhanea macrorchis (Monogenea: Chauhaneidae) revealed a highly repetitive region and a gene rearrangement hot spot in Polyopisthocotylea. Mol Biol Rep. 2012; doi:10.1007/s11033-012-1659-z.
Nakao M, McManus DP, Schantz PM, Craig PS, Ito A. A molecular phylogeny of the genus Echinococcus inferred from complete mitochondrial genomes. Parasitology. 2007; doi:10.1017/S003118200600193.
von Nickisch-Rosenegk M, Brown WM, Boore JL. Complete sequence of the mitochondrial genome of the tapeworm Hymenolepis diminuta: gene arrangements indicate that Platyhelminths are Eutrochozoans. Mol Biol Evol. 2001;18:721–30.
MAFFT version 7. Multiple alignment program for amino acid or nucleotide sequences. http://mafft.cbrc.jp/alignment/server/.
Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17:540–52. http://molevol.cmima.csic.es/castresana/Gblocks_server.html.
Sela I, Ashkenazy H, Katoh K, Pupko T. GUIDANCE2: accurate detection of unreliable alignment regions accounting for the uncertainty of multiple parameters. Nucl Acids Res. 2015; doi:10.1093/nar/gkq443. http://guidance.tau.ac.il/ver2/.
Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012; doi:10.1093/molbev/mss020
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst Biol. 2010; doi: 10.1093/sysbio/syq010. http://www.atgc-montpellier.fr/phyml/.
Bernt M, Merkle D, Ramsch K, Fritzsch G, Perseke M, Bernhard D, Schlegel M, Stadler P, Middendorf M. CREx: Inferring Genomic Rearrangements Based on Common Intervals. Bioinformatics. 2007; doi:10.1093/bioinformatics/btm468.
Telford MJ, Herniou EA, Russell RB, Littlewood DT. Changes in mitochondrial genetic codes as phylogenetic characters: two examples from the flatworms. Proc Natl Acad Sci USA. 2000; doi: 10.1073/pnas.97.21.11359.
Wey-Fabrizius AR, Podsiadlowski L, Herlyn H, Hankeln T. Platyzoan mitochondrial genomes. Mol Phylogen Evol. 2013; doi:10.1016/j.ympev.2012.12.015.
Steinauer ML, Nickol BB, Broughton R, Orti G. First sequenced mitochondrial genome from the phylum Acanthocephala (Leptorhynchoides thecatus) and its phylogenetic position within Metazoa. J Mol Evol. 2005; doi:10.1007/s00239-004-0159-8.
Fromm B, Burow S, Hahn C, Bachmann L. MicroRNA loci support conspecificity of the parasites Gyrodactylus salaris and G. thymalli (Platyhelminthes: Monogenea). Int J Parasitol. 2014; doi:10.1016/j.ijpara.2014.05.010.
We are grateful for having access to the Abel high computing facility of the University of Oslo (UiO) through the Norwegian Metacenter for Computational Science (NOTUR; grant k9201nn) to LB. Article processing charges were kindly covered by the UiO publishing fund. Additional funding was provided from a grant of the Conselho de Desenvolvimento Científico e Tecnológico (CNPq), Brazil to WAB.
Base composition of the mtDNA. (DOCX 53 kb)
Codon usage of the protein coding genes (PCGs). (DOCX 107 kb)
The secondary structures of the tRNAs. (PDF 184 kb)
Maximum Likelihood trees after removal of ambiguous sites. (PDF 119 kb)
P-distances of concatenated mitochondrial genes. (DOCX 17 kb)
Schematic order of the mitochondrial genes. (PDF 113 kb)
Nucleotide sequence alignment of the repeat regions I and II. (DOCX 82 kb)
CREx distance matrix for mitochondrial gene order rearrangements. (DOCX 16 kb)