Characterization of the first mitochondrial genome of a catenulid flatworm: Stenostomum leucops (Platyhelminthes)

The first complete mitochondrial genome of a catenulid, Stenostomum leucops, was characterized. Illumina sequencing and 90 813 reads were utilized in the assembly, producing one contig with an average coverage of 1118×. The length of this genome is 15 742 bp with 12 protein‐coding, two rRNA and 22 tRNA genes. Although the atp8 gene is absent in other Platyhelminthes, a highly divergent putative atp8 gene was found in S. leucops. In contrast to other Platyhelminthes, the mitochondrial genes of S. leucops are encoded on both strands. The gene order in the S. leucops mitogenome is very divergent from those observed in other Platyhelminthes, showing only small blocks of synteny. With AAA as the codon for lysine S. leucops shows the probable plesiomorphic condition, whereas Rhabditophora possess the derived GAA. This evolutionary transition is correlated with changes in the respective anticodons in trnK. It remains unclear whether the absence of the D arm and loop in trnS1 is a convergence in Catenulida and Neodermata.     

The complete mitochondrial genomes of two vent squat lobsters,Munidopsis lauensis and M. verrilli: Novel gene arrangements and phylogenetic implications

Hydrothermal vents are considered as one of the most extremely harsh environments on the Earth. In this study, the complete mitogenomes of hydrothermal vent squat lobsters, Munidopsis lauensis and M. verrilli, were determined through Illumina sequencing and compared with other available mitogenomes of anomurans. The mitogenomes of M. lauensis (17,483 bp) and M. verrilli (17,636 bp) are the largest among all Anomura mitogenomes, while the A+T contents of M. lauensis (62.40%) and M. verrilli (63.99%) are the lowest. The mitogenomes of M. lauensis and M. verrilli display novel gene arrangements, which might be the result of three tandem duplication–random loss (tdrl) events from the ancestral pancrustacean pattern. The mitochondrial gene orders of M. lauensis and M. verrilli shared the most similarities with S. crosnieri. The phylogenetic analyses based on both gene order data and nucleotide sequences (PCGs and rRNAs) revealed that the two species were closely related to Shinkaia crosnieri. Positive selection analysis revealed that eighteen residues in seven genes (atp8, Cytb, nad3, nad4, nad4l, nad5, and nad6) of the hydrothermal vent anomurans were positively selected sites.     

The complete mitochondrial genome of Tuta absoluta (Lepidoptera: Gelechiidae) and genetic variation in two newly invaded populations in China

Tuta absoluta (Meyrick) (Lepidoptera: Gelechiidae) is a devastating invasive pest worldwide, causing severe damage to tomatoes. Recently, it has been recorded in the northwestern and southwestern parts of China. Here, the mitogenomes and genetic variation of two newly invaded T. absoluta populations in Xinjiang and Yunnan, were determined. The results showed that the complete mitogenome size of T. absoluta is 15298 bp for the individual from Xinjiang and 15296 bp for the individual from Yunnan, which were both longer than the reported mitogenome from Spain (15290 bp). The mitogenome sequences of individuals collected from three locations showed high levels of sequence similarity, except for 8 polymorphic sites, which were in genes cox2 (1 site), cox3 (2 sites), cob (1 site), atp6 (1 site), nad1 (2 sites) and nad5 (1 site). Tuta absoluta mitogenomes share many features with other 6 Gelechiidae mitogenomes, except for several differences in the start and stop codons of protein-coding genes and the length of intergenic spacers. Seven partial mitochondrial genes (cox1, cox2, cox3, atp6, cob, nad1, and nad5) were used for genetic variation analysis, and significant population differentiation was found between the two populations based on cox2, atp6, nad1, and nad5. The complete mitogenomes and sensitive mitochondrial gene markers reported here provide useful data for further population genetics study of this pest.     

Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera)

The complete arrangement of genes in the mitochondrial (mt) genome is known for 12 species of insects, and part of the gene arrangement in the mt genome is known for over 300 other species of insects. The arrangement of genes in the mt genome is very conserved in insects studied, since all of the protein-coding and rRNA genes and most of the tRNA genes are arranged in the same way. We sequenced the entire mt genome of the wallaby louse, Heterodoxus macropus, which is 14,670 bp long and has the 37 genes typical of animals and some noncoding regions. The largest noncoding region is 73 bp long (93% A+T), and the second largest is 47 bp long (92% A+T). Both of these noncoding regions seem to be able to form stem-loop structures. The arrangement of genes in the mt genome of this louse is unlike that of any other animal studied. All tRNA genes have moved and/or inverted relative to the ancestral gene arrangement of insects, which is present in the fruit fly Drosophila yakuba. At least nine protein-coding genes (atp6, atp8, cox2, cob, nad1-nad3, nad5, and nad6) have moved; moreover, four of these genes (atp6, atp8, nad1, and nad3) have inverted. The large number of gene rearrangements in the mt genome of H. macropus is unprecedented for an arthropod.     

Mitochondrial genome sequence of <Emphasis Type="Italic">Unionicola foili</Emphasis> (Acari: Unionicolidae): a unique gene order with implications for phylogenetic inference

The mitochondrial genome of Unionicola foili is circular, 14,738 bp in length, and contains several notable features. The sequence and annotation revealed a unique gene order, continuing a pattern of highly rearranged mitochondrial genomes in the Trombidiformes. U. foili mitochondrial tRNA sequences predict non-canonical secondary structures for these molecules, and our annotation suggests an in-frame fusion between the nad4L and nad5 genes in this genome. The unique gene order and unusual tRNA structures could serve as idiosyncratic characters and have the potential to be phylogenetically informative.     

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.     

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.     

The mitochondrial genome of <Emphasis Type="Italic">Euphausia superba</Emphasis> (Prydz Bay) (Crustacea: Malacostraca: Euphausiacea) reveals a novel gene arrangement and potential molecular markers

Euphausiid krill are dominant organisms in the zooplankton population and play a central role in marine ecosystems. In this paper, we described the gene organization, gene rearrangement and codon usage in the mitochondrial genome of Euphausia superba Dana 1852 (sampling from Prydz Bay, PB). The mitochondrial genome of E. superba is more than 15,498 bp in length (partial non-coding region was not determined). Translocation of four tRNAs (trnL ₁ , trnL ₂ , trnW and trnI) and duplication of one tRNA (trnN) were founded in the mitochondrial genome of E. superba when comparing its genome with the pancrustacean ground pattern. To investigate the phylogenetic relationship within Malacostraca, phylogenetic trees based on currently available malacostracan mitochondrial genomes were built with the maximum likelihood and the Bayesian models. All analyses based on nucleotide and amino acid data strongly support the monophyly of Stomatopoda, Penaeidae, Caridea, and Brachyura, which is consistent with previous research. However, the taxonomic position of Euphausiacea within Malacostraca is unstable. From comparing the mitochondrial genome between E. superba (PB) and E. superba (sampling from Weddell Sea, WS), we found that nad2 gene contains maximal variation with 61 segregating sites, following by nad5 gene which has 12 segregating sites. Thus, nad2 and nad5 genes may be used as potential molecular markers to study the inherit diversity among different E. superba groups, which would be helpful to the exploitation and management of E. superba resources.     

The complete mitochondrial genome of the oriental fruit moth <Emphasis Type="Italic">Grapholita molesta</Emphasis> (Busck) (Lepidoptera: Tortricidae)

The oriental fruit moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae) currently is one of the economically most destructive pest species of stone and pome fruits worldwide. Here we sequenced the complete mitochondrial genome of this pest. This genome is 15,776 bp long, with an A + T content of 81.24%, containing 37 typical animal mitochondrial genes and an A + T-rich region. All gene are arranged as hypothesized ancestral gene order of insects except for trnM, which was shuffled from 3′ downstream of trnQ to 5′ upstream of trnI. cox1 gene uses unusual CGA start codon, as that in all other sequenced lepidopteran mitochondrial genome. The secondary structures for the two rRNA genes were predicted. All helices typically present in insect mitochondrial rRNA genes are generated. A microsatellite sequence was inserted into the region of H2347 in rrnL in G. molesta and two other sequenced tortricid mitochondrial genomes, indicating that the insertion event in this helix might occurred anciently in family Tortricidae. All of the 22 typical animal tRNA genes have a typical cloverleaf structure except for trnS2, in which the D-stem pairings in the DHU arm are absent. An intergenic sequence is present between trnQ and nad2 as well as in other sequenced lepidopteran mitochondrial genomes, which was presumed to be a remnant of trnM gene and its boundary sequences after the duplication of trnM to the upstream of trnI in Lepidoptera. The A + T-rich region is 836 bp, containing six repeat sequences of “TTATTATTATTATTAAATA(G)TTT.”     

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 mitochondrial genome of Priapulus caudatus Lamarck (Priapulida: Priapulidae)

We sequenced and annotated the complete mitochondrial (mt) genome of the priapulid Priapulus caudatus in order to provide a source of phylogenetic characters including an assessment of gene order arrangement. The genome was 14,919 bp in its entirety with few, short non-coding regions. A number of protein-coding and tRNA genes overlapped, making the genome relatively compact. The gene order was: cox1, cox2, trnK, trnD, atp8, atp6, cox3, trnG, nad3, trnA, trnR, trnN, rrnS, trnV, rrnL, trnL(yaa), trnL(nag), nad1, -trnS(nga), -cob, -nad6, trnP, -trnT, nad4L, nad4, trnH, nad5, trnF, -trnE, -trnS(nct), trnI, -trnQ, trnM, nad2, trnW, -trnC, -trnY; where '-' indicates genes transcribed on the opposite strand. The gene order, although unique amongst Metazoa, shared the greatest number of gene boundaries and the longest contiguous fragments with the chelicerate Limulus polyphemus. The mt genomes of these taxa differed only by a single inversion of 18 contiguous genes bounded by rrnS and trnS(nct). Other arthropods and nematodes shared fewer gene boundaries but considerably more than the most similar non-ecdysozoan.     

Complete mitogenome of Nycteribia allotopa Speiser, 1901 (Diptera,Hippoboscoidea, Nycteribiidae) and comparative analysis of mitochondrial genomes of Nycteribiidae

In recent years, the global pandemic of bat-associated pathogens has led to increasing attention on bat ectoparasites. Numerous studies have identified human-associated pathogens in Nycteribiidae, indicating their potential as vectors. In this study, the first complete sequencing of the mitochondrial genome of Nycteribia allotopa Speiser, 1901 was sequenced and analyzed. We also compared the mitochondrial sequences of N. allotopa with those available in the database for other Nycteribiidae species. The complete mitochondrial genome of N. allotopa was found to be 15,161 bp in size with an A + T content of 82.49%. Nucleotide polymorphism analysis of 13 protein-coding genes from five species of Nycteribiidae showed that nad6 exhibited the most significant variation, while cox1 was the most conserved. Furthermore, selection pressure analysis revealed cox1 to exhibit the strongest purifying selection, while atp8, nad2, nad4L, and nad5 showed slightly looser purifying selection. Pairwise genetic distances indicated that cox1 and cox2 were evolving comparatively slowly, whereas atp8, nad2, and nad6 were evolving comparatively quickly. Phylogenetic trees constructed using Bayesian inference and maximum likelihood methods demonstrated that all four families within the superfamily Hippoboscoidea clustered into one branch each, indicating their monophyly. N. allotopa was found to be most closely related to the same genus N. parvula. This study significantly enriches the molecular database for Nycteribiidae and provides invaluable reference data for future species identification, phylogenetic analysis, and exploration of their potential as vectors for human-associated pathogens.     

The rearranged mitochondrial genome of Podagrion sp. (Hymenoptera: Torymidae), a parasitoid wasp of mantis

The complete sequence of the mitochondrial genome of Podagrion sp. (Hymenoptera: Torymidae) is described. The mitogenome was 15,845 bp in size, and contained typical sets of mitochondrial genes. The base composition of the Podagrion sp. mitogenome was also biased toward A + T bases (81.8%). The mitochondrial genome of Podagrion sp. has a weak AT skew (0.07) and a strong GC skew (?0.26). Podagrion sp. exhibits a novel rearrangement compared with the ancestral order, including six protein-coding genes (nad3, cox3, atp6, atp8, cox2 and cox1), which have inverted to the minor strand from the major strand. The A + T-rich region of Podagrion sp., which is located between trnN and trnI, have five tandem repeats. The apomorphic rearrangements, including the conserved block “cox3-atp6-atp8-cox2-cox1-nad5-nad4-nad4l-nad6-cob” and the special locations of trnV and trnA, were mapped onto the phylogeny of Proctotrupomorpha.     

Phylogeny of Rhigonematomorpha based on the complete mitochondrial genome of Rhigonema thysanophora (Nematoda: Chromadorea)

We determined the complete mitochondrial genome sequence of Rhigonema thysanophora, the first representative of Rhigonematomorpha, and used this sequence along with 57 other nematode species for phylogenetic analyses. The R. thysanophora mtDNA is 15 015 bp and identical to all other chromadorean nematode mtDNAs published to date in that it contains 36 genes (lacking atp8) encoded in the same direction. Phylogenetic analyses of nucleotide and amino acid sequence data for the 12 protein‐coding genes recovered Rhigonematomorpha as the sister group to the heterakoid species, Ascaridia columbae (Ascaridomorpha). The organization of R. thysanophora mtDNA resembles the most common pattern for the Rhabditomorpha+Ascaridomorpha+Diplogasteromorpha clade in gene order, but with some substantial gene rearrangements. This similarity in gene order is in agreement with the sequence‐based analyses that indicate a close relationship between Rhigonematomorpha and Rhabditomorpha+Ascaridomorpha+Diplogasteromorpha. These results are consistent with certain analyses of nuclear SSU rDNA for R. thysanophora and some earlier classification systems that asserted phylogenetic affinity between Rhigonematomorpha and Ascaridomorpha, but inconsistent with morphology‐based phylogenetic hypotheses that suggested a close (taxonomic) relationship between rhigonematomorphs and oxyuridomorphs (pinworms). These observations must be tempered by noting that few rhigonematomorph species have been sequenced and included in phylogenetic analyses, and preliminary studies based on SSU rDNA suggest the group is not monophyletic. Additional mitochondrial genome sequences of rhigonematids are needed to characterize their phylogenetic relationships within Chromadorea, and to increase understanding of mitochondrial genome evolution.     

Complete mitochondrial genome of Coelomactra antiquata (Mollusca: Bivalvia): The first representative from the family Mactridae with novel gene order and unusual tandem repeats

The complete mitochondrial genome plays an important role in the accurate inference of phylogenetic relationships among metazoans. Mactridae, also known as trough shells or duck clams, is an important family of marine bivalve clams in the order Veneroida. Here we present the complete mitochondrial genome sequence of the Xishishe Coelomactra antiquata (Mollusca: Bivalvia), which is the first representative from the family Mactridae. The mitochondrial genome of C. antiquata is of 17,384 bp in length, and encodes 35 genes, including 12 protein-coding, 21 transfer RNA, and 2 ribosomal RNA genes. Compared with the typical gene content of animal mitochondrial genomes, atp8 and tRNAS₂ are missing. Gene order of the mitochondrial genome of C. antiquata is unique compared with others from Veneroida. In the mitochondrial genome of the C. antiquata, a total of 2189 bp of non-coding nucleotides are scattered among 26 non-coding regions. The largest non-coding region contains one section of tandem repeats (99 bp × 11), which is the second largest tandem repeats found in the mitochondrial genomes from Veneroida. The phylogenetic trees based on mitochondrial genomes support the monophyly of Veneridae and Lucinidae, and the relationship at the family level: ((Veneridae + Mactridae) + (Cardiidae + Solecurtidae)) + Lucinidae. The phylogenetic result is consistent with the morphological classification. Meanwhile, bootstrap values are very high (BP = 94–100), suggesting that the evolutionary relationship based on mitochondrial genomes is very reliable.     

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.     

Characterization of mitochondrial genome of sea cucumber <Emphasis Type="Italic">Stichopus horrens</Emphasis>: A novel gene arrangement in Holothuroidea

The complete mitochondrial DNA sequence contains useful information for phylogenetic analyses of metazoa. In this study, the complete mitochondrial DNA sequence of sea cucumber Stichopus horrens (Holothuroidea: Stichopodidae: Stichopus) is presented. The complete sequence was determined using normal and long PCRs. The mitochondrial genome of Stichopus horrens is a circular molecule 16257 bps long, composed of 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes. Most of these genes are coded on the heavy strand except for one protein-coding gene (nad6) and five tRNA genes (tRNA ^Ser(UCN) , tRNA ^Gln , tRNA ^Ala , tRNA ^Val , tRNA ^Asp ) which are coded on the light strand. The composition of the heavy strand is 30.8% A, 23.7% C, 16.2% G, and 29.3% T bases (AT skew=0.025; GC skew=−0.188). A non-coding region of 675 bp was identified as a putative control region because of its location and AT richness. The intergenic spacers range from 1 to 50 bp in size, totaling 227 bp. A total of 25 overlapping nucleotides, ranging from 1 to 10 bp in size, exist among 11 genes. All 13 protein-coding genes are initiated with an ATG. The TAA codon is used as the stop codon in all the protein coding genes except nad3 and nad4 that use TAG as their termination codon. The most frequently used amino acids are Leu (16.29%), Ser (10.34%) and Phe (8.37%). All of the tRNA genes have the potential to fold into typical cloverleaf secondary structures. We also compared the order of the genes in the mitochondrial DNA from the five holothurians that are now available and found a novel gene arrangement in the mitochondrial DNA of Stichopus horrens.     

Characterization of the complete mitochondrial genome of Japanagallia spinosa and Durgades nigropicta (Hemiptera: Cicadellidae: Megophthalminae)

In this study, we determined and analyzed the complete mitochondrial genomes (mitogenomes) of Japanagallia spinosa and Durgades nigropicta (Hemiptera: Megophthalminae). The circular genome were 15,655 bp long in J. spinosa (GenBank: KY123686) and 15,974 bp long in D. nigropicta (GenBank: KY123687). The J. spinosa and D. nigropicta mitogenomes both contained 37 genes and the gene order was similar to that in other leafhoppers. All of the protein-coding genes started with ATN. In the J. spinosa mitogenome, the nad3, nad4L, and cytb genes used TAG as a stop codon, the atp8 and nad1 genes used TGA, and the cox2 gene used a single T. However, in the D. nigropicta mitogenome, three genes used a single T as the stop codon, whereas the nad3 gene used TAG. We predicted the secondary structures of the rRNAs in J. spinosa and D. nigropicta. The secondary structure of rrnL comprised six domains (domain III is absent in arthropods) with 42 helices and that of rrnS comprised three structural domains with 26 helices. Comparisons of J. spinosa and D. nigropicta detected some differences in H577 and H673. We determined the structural organization of the control regions in the mitogenomes of leafhoppers, where three types of repeat regions were found in most. The phylogenetic relationships between J. spinosa and D. nigropicta with related lineages were reconstructed using Bayesian inference and maximum likelihood analyses. The monophyly of each superfamily considered in this study was confirmed by the clades in the phylogenetic tree. And in this study, Cicadellidae was resolved as monophyletic by the phylogenetic analysis. This mitogenome information for J. spinosa and D. nigropicta could facilitate future studies of mitogenomic diversity and the evolution of related insect lineages.