Trichinella spiralis mtDNA: a nematode mitochondrial genome that encodes a putative ATP8 and normally structured tRNAS and has a gene arrangement relatable to those of coelomate metazoans

High-Throughput Sequencing Enables Rapid Analyses of Nematode Mitochondrial Genomes from an Environmental Sample

Mitochondrial genomes serve as essential tools in evolutionary biology, phylogenetics, and population genetics due to their maternal inheritance, lack of recombination, and conserved structure. Traditional morphological methods for identifying nematodes are often insufficient for distinguishing cryptic species complexes. This study highlights recent advancements in nematode mitochondrial genome research, particularly the impact of long-read sequencing technologies such as Oxford Nanopore. These technologies have facilitated the assembly of mitochondrial genomes from mixed soil samples, overcoming challenges associated with designing specific primers for long PCR amplification across different groups of parasitic nematodes. In this study, we successfully recovered and assembled eleven nematode mitochondrial genomes using long-read sequencing, including those of two plant-parasitic nematode species. Notably, we detected Heterodera cruciferae in Victoria, expanding its known geographic range within Australia. Additionally, short-read sequencing data from a previous draft genome study revealed the presence of the mitochondrial genome of Heterodera filipjevi. Comparative analyses of Heterodera mitogenomes revealed conserved protein-coding genes essential for oxidative phosphorylation, as well as gene rearrangements and variations in transfer RNA placement, which may reflect adaptations to parasitic lifestyles. The consistently high A+T content and strand asymmetry observed across species align with trends reported in related genera. This study demonstrates the utility of long-read sequencing for identifying coexisting nematode species in agricultural fields, providing a rapid, accurate, and comprehensive alternative to traditional diagnostic methods. By incorporating non-target endemic species into public databases, this approach enhances biodiversity records and informs biosecurity strategies. These findings reinforce the potential of mitochondrial genomics to strengthen Australia's as well as the global biosecurity framework against plant-parasitic nematode threats.

Comparative mitogenomics of Cheiracanthium species (Araneae: Cheiracanthiidae) with phylogenetic implication and evolutionary insights

The genus Cheiracanthium C. L. Koch, 1839 is the most species-rich genus of the family Cheiracanthiidae. Given the unavailability of information on the evolutionary biology and molecular taxonomy of this genus, here we sequenced nine mitochondrial genomes (mitogenomes) of Cheiracanthium species, four of which were fully annotated, and conducted comparative analyses with other well-characterized Araneae mitogenomes. We also provide phylogenetic insights on the genus Cheiracanthium. The circular mitogenomes of the Cheiracanthium contain 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and one putative control region (CR). All genes show a high A+T bias, characterized by a negative AT skew and positive GC skew, along with numerous overlapped regions and intergenic spacers. Approximately half of the tRNAs lack TΨC and/or dihydrouracil (DHU) arm and are characterized with unpaired amino acid acceptor arms. Most PCGs used the standard ATN start codons and TAR termination codons. The mitochondrial gene order of Cheiracanthium differs significantly from the putative ancestral gene order (Limulus polyphemus). Our novel phylogenetic analyses infer Cheiracanthiidae to be the sister group of Salticidae in BI analysis, but as sister to the node with Miturgidae, Viridasiidae, Corinnidae, Selenopidae, Salticidae, and Philodromidae in ML analysis. We confirm that Cheiracanthium is paraphyletic, for the first time using molecular phylogenetic approaches, with the earliest divergence estimated at 67 Ma. Our findings enhance our understanding of Cheiracanthium taxonomy and evolution.

Characterization and phylogenetic analysis of the mitochondrial genome of Cylicostephanus longibursatus

Cylicostephanus longibursatus is a common parasite in equine animals. Hosts infected by these nematodes might face disease or death. This study utilized next-generation sequencing technology to sequence the complete mitochondrial genome (mt genome) of C. longibursatus. Through bioinformatics techniques, the genomic base composition, codon usage, tRNA secondary structures, evolutionary relationships, and taxonomic status were analyzed. The results revealed that the mitochondrial genome of C. longibursatus is a double-stranded, 13,807-bp closed circular molecule with an AT content of 76.0%, indicating a clear preference for AT bases. The mitochondrial genome consisted of a total of 12 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 2 non-coding regions. Among the 12 protein-coding genes, TTG and ATT were the common start codons. TAA was the predominant termination codon, except for the ND3 and ND6 coding genes, and the COШ genes used TAG and "T" as termination codons, respectively. All tRNAs exhibited atypical clover-leaf secondary structures, except for tRNALys and tRNALeu2, where two tRNASer genes lacked DHU arms and DHU loops, tRNAmet lacked the TΨC-arm, tRNAIle lacked the TΨC-loop, and the remaining 16 tRNAs lacked the TΨC-arm and TΨC loop, which were substituted by the "TV-replacement loop". Phylogenetic analyses, based on the 12 protein-coding genes and utilizing maximum likelihood (ML) and Bayesian inference (BI) analyses, indicated that C. longibursatus did not form a monophyletic group with other Cylicostephanus but was instead more closely related to Cyathostomum. These research findings provide fundamental data for exploring the population classification and phylogeny of strongylid nematodes.

Mitogenome selection shaped the terrestrial adaptation of Grapsidae (Decapoda: Brachyura)

The colonization of aquatic to terrestrial habitats by brachyuran crabs requires genetic innovations as well as morphological adaptations to adapt to terrestrial environments. The genetic basis of such adaptive evolution, however, is largely unknown. This study focuses on terrestrialization in Geograpsus (Grapsidae) the only highly terrestrial genus in this family, which represents a notable example of terrestrial adaptive radiation. Here, we sequenced the mitogenomes of two Geograpsus species and used the mitogenomes of 215 representative crabs to construct phylogenetic and time frameworks that we used to infer terrestrial origins and evolution. Using mitochondrial genomic data, we demonstrated that marine crab ancestors began to settle on land during the early Eocene. Ocean acidification, the Paleocene-Eocene Thermal Maximum (PETM), and mangrove expansion at that time may have driven the diversification and ecological expansion of these terrestrial crabs. Evolutionary analyses reveal strong positive selection signals on monophyletic lineages of Grapsidae, especially the terrestrial species of Geograpsus. Positively selected sites in functionally important regions of ND5 and ND4 may imply enhanced energy metabolism in Grapsidae compared to other crabs, and may have played an important role in their terrestrial adaptation. Overall, our work provides valuable resources and opportunities to reveal the adaptation of crabs to complex terrestrial environments.

Both-strand gene coding in a plastome-like mitogenome of an enoplid nematode

The phylum Nematoda remains very poorly sampled for mtDNA, with a strong bias toward parasitic, economically important or model species of the Chromadoria lineage. Most chromadorian mitogenomes share a specific order of genes encoded on one mtDNA strand. However, the few sequenced representatives of the Dorylaimia lineage exhibit a variable order of mtDNA genes encoded on both strands. While the ancestral arrangement of nematode mitogenome remains undefined, no evidence has been reported for Enoplia, the phylum's third early divergent major lineage. We describe the first mitogenome of an enoplian nematode, Campydora demonstrans, and contend that the complete 37-gene repertoire and both-strand gene encoding are ancestral states preserved in Enoplia and Dorylaimia versus the derived mitogenome arrangement in some Chromadoria. The C. demonstrans mitogenome is 17,018 bp in size and contains a noncoding perfect inverted repeat with 2013 bp-long arms, subdividing the mitogenome into two coding regions. This mtDNA arrangement is very rare among animals and instead resembles that of chloroplast genomes in land plants. Our report broadens mtDNA taxonomic sampling of the phylum Nematoda and adds support to the applicability of cox1 gene as a phylogenetic marker for establishing nematode relationships within higher taxa.

Molecular phylogeny of the family Rhabdiasidae (Nematoda: Rhabditida), with morphology, genetic characterization and mitochondrial genomes of Rhabdias kafunata and R. bufonis

BACKGROUND

The family Rhabdiasidae (Nematoda: Rhabditida) is a globally distributed group of nematode parasites, with over 110 species parasitic mainly in amphibians and reptiles. However, the systematic position of the family Rhabdiasidae in the order Rhabditida remains unsolved, and the evolutionary relationships among its genera are still unclear. Moreover, the present knowledge of the mitochondrial genomes of rhabdiasids remains limited.

METHODS

Two rhabdiasid species: Rhabdias kafunata Sata, Takeuchi & Nakano, 2020 and R. bufonis (Schrank, 1788) collected from the Asiatic toad Bufo gargarizans Cantor (Amphibia: Anura) in China, were identified based on morphology (light and scanning electron microscopy) and molecular characterization (sequencing of the nuclear 28S and ITS regions and mitochondrial cox1 and 12S genes). The complete mitochondrial genomes of R. kafunata and R. bufonis were also sequenced and annotated for the first time. Moreover, phylogenetic analyses based on the amino acid sequences of 12 protein-coding genes (PCGs) of the mitochondrial genomes were performed to clarify the systematic position of the family Rhabdiasidae in the order Rhabditida using maximum likelihood (ML) and Bayesian inference (BI). The phylogenetic analyses based on the 28S + ITS sequences, were also inferred to assess the evolutionary relationships among the genera within Rhabdiasidae.

RESULTS

The detailed morphology of the cephalic structures, vulva and eggs in R. kafunata and R. bufonis was revealed using scanning electron microscopy (SEM) for the first time. The characterization of 28S and ITS regions of R. kafunata was reported for the first time. The mitogenomes of R. kafunata and R. bufonis are 15,437 bp and 15,128 bp long, respectively, and both contain 36 genes, including 12 PCGs (missing atp8). Comparative mitogenomics revealed that the gene arrangement of R. kafunata and R. bufonis is different from all of the currently available mitogenomes of nematodes. Phylogenetic analyses based on the ITS + 28S data showed Neoentomelas and Kurilonema as sister lineages, and supported the monophyly of Entomelas, Pneumonema, Serpentirhabdias and Rhabdias. Mitochondrial phylogenomic results supported Rhabdiasidae as a member of the superfamily Rhabditoidea in the suborder Rhabditina, and its occurrance as sister to the family Rhabditidae.

CONCLUSIONS

The complete mitochondrial genome of R. kafunata and R. bufonis were reported for the first time, and two new gene arrangements of mitogenomes in Nematoda were revealed. Mitogenomic phylogenetic results indicated that the family Rhabdiasidae is a member of Rhabditoidea in Rhabditina, and is closely related to Rhabditidae. Molecular phylogenies based on the ITS + 28S sequence data supported the validity of Kurilonema, and showed that Kurilonema is sister to Neoentomelas. The present phylogenetic results also indicated that the ancestors of rhabdiasids seem to have initially infected reptiles, then spreading to amphibians.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Rhabdias kafunata (Rhabditida: Rhabdiasidae)

Rhabdias kafunata (Rhabditida: Rhabdiasidae) is a parasitic nematode that significantly affects bufonids. To better understand the genome-level characteristics of related species, Illumina sequencing was used to identify mitochondrial genes and analyze their basic characteristics and gene arrangements. The mitogenome of R. kafunata is 14,068 bp in length and contains 36 genes, including 12 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one noncoding region (NCR). The nucleotide composition is highly biased toward A + T, accounting for 75.5% of the entire mitochondrial genome. The cox1 sequence is relatively conserved in Ka/Ks analyses and can be used as a gene fragment for species identification. While 8 of the 12 PCGs use the typical ATN initiation codon, nad1-2, nad4, and cox3 utilize a TTG initiation codon. Most stop codons end with the standard TAA or TAG, except for cytb, which ends with an incomplete TA. Additionally, trnM, trnK, and trnI have the typical clover-leaf secondary structure, while the remaining tRNAs lack the DHU arm or TΨC arm. Phylogenetic analysis indicates that R. kafunata belongs to the Rhabditidae family and is closely related to Litoditis marina and Caenorhabditis angaria among sequenced lepidopteran mitochondrial genomes.

Two-Fold ND5 Genes, Three-Fold Control Regions, lncRNA, and the "Missing" ATP8 Found in the Mitogenomes of Polypedates megacephalus (Rhacophridae: Polypedates)

In prior research on the mitochondrial genome (mitogenome) of Polypedates megacephalus, the one copy of ND5 gene was translocated to the control region (CR) and the ATP8 gene was not found. Gene loss is uncommon among vertebrates. However, in this study, we resequenced the mitogenomes of P. megacephalus from different regions using a "primer bridging" approach with Sanger sequencing technologies, which revealed the "missing" ATP8 gene in P. megacephalus as well as three other previously published Polypedates. The mitogenome of this species was found to contain two copies of the ND5 genes and three copies of the control regions. Furthermore, multiple tandem repeats were identified in the control regions. Notably, we observed that there was no correlation between genetic divergence and geographic distance. However, using the mitogenome, gene expression analysis was performed via RT-qPCR of liver samples and it was thus determined that COIII, ND2, ND4, and ND6 were reduced to 0.64 ± 0.24, 0.55 ± 0.34, 0.44 ± 0.21 and 0.65 ± 0.17, respectively, under low-temperature stress (8 °C) as compared with controls (p < 0.05). Remarkably, the transcript of long non-coding RNA (lncRNA) between positions 8029 and 8612 decreased significantly with exposure to low-temperature stress (8 °C). Antisense ND6 gene expression showed a downward trend, but this was not significant. These results reveal that modulations of protein-coding mitochondrial genes and lncRNAs of P. megacephalus play a crucial role in the molecular response to cold stress.

Mitogenomic phylogenies suggest the resurrection of the subfamily Porrocaecinae and provide insights into the systematics of the superfamily Ascaridoidea (Nematoda: Ascaridomorpha), with the description of a new species of Porrocaecum

BACKGROUND

The family Toxocaridae is a group of zooparasitic nematodes of veterinary, medical and economic significance. However, the evolutionary relationship of Porrocaecum and Toxocara, both genera currently classified in Toxocaridae, and the monophyly of the Toxocaridae remain under debate. Moreover, the validity of the subgenus Laymanicaecum in the genus Porrocaecum is open to question. Due to the scarcity of an available genetic database, molecular identification of Porrocaecum nematodes is still in its infancy.

METHODS

A number of Porrocaecum nematodes collected from the Eurasian marsh harrier Circus aeruginosus (Linnaeus) (Falconiformes: Accipitridae) in the Czech Republic were identified using integrated morphological methods (light and scanning electron microscopy) and molecular techniques (sequencing and analyzing the nuclear 18S, 28S and ITS regions). The complete mitochondrial genomes of the collected nematode specimens and of Porrocaecum (Laymanicaecum) reticulatum (Linstow, 1899) were sequenced and annotated for the first time. Phylogenetic analyses of ascaridoid nematodes based on the amino acid sequences of 12 protein-coding genes of mitochondrial genomes were performed using maximum likelihood and Bayesian inference.

RESULTS

A new species of Porrocaecum, named P. moraveci n. sp., is described based on the morphological and genetic evidence. The mitogenomes of P. moraveci n. sp. and P. reticulatum both contain 36 genes and are 14,517 and 14,210 bp in length, respectively. Comparative mitogenomics revealed that P. moraveci n. sp. represents the first known species with three non-coding regions and that P. reticulatum has the lowest overall A + T content in the mitogenomes of ascaridoid nematodes tested to date. Phylogenetic analyses showed the representatives of Toxocara clustered together with species of the family Ascarididae rather than with Porrocaecum and that P. moraveci n. sp. is a sister to P. reticulatum.

CONCLUSIONS

The characterization of the complete mitochondrial genomes of P. moraveci n. sp. and P. reticulatum is reported for the first time. Mitogenomic phylogeny analyses indicated that the family Toxocaridae is non-monophyletic and that the genera Porrocaecum and Toxocara do not have an affinity. The validity of the subgenus Laymanicaecum in Porrocaecum was also rejected. Our results suggest that: (i) Toxocaridae should be degraded to a subfamily of the Ascarididae that includes only the genus Toxocara; and (ii) the subfamily Porrocaecinae should be resurrected to include only the genus Porrocaecum. The present study enriches the database of ascaridoid mitogenomes and provides a new insight into the systematics of the superfamily Ascaridoidea.

The complete mitochondrial genomes of the Leuciscus baicalensis and Rutilus rutilus: a detailed genomic comparison among closely related species of the Leuciscinae subfamily

Cyprinidae is the largest family in the order of freshwater fish Cypriniformes. Increased subfamily members of Cyprinidae have been suggested to be re-classified for decades. In this study, we sequenced the mitochondrial genomes (mitogenomes) of Leuciscus baicalensis and Rutilus rutilus collected from northwest China and compared with other closely related species to determine their associated family or subfamily. We used Illumina NovaSeq to sequence the entire mitochondrial genomes of Leuciscus baicalensis and Rutilus rutilus and characterized the mitogenomes by the gene structure, gene order, and the secondary structures of the 22 tRNA genes. We compared mitogenome features of Leuciscinae with other subfamilies in Cyprinidae. We used the analytic Bayesian Information and Maximum Likelihood methods to determine phylogenetic trees of 13 PCGs. The mitogenomes of Leuciscus baicalensis and Rutilus rutilus were 16,607 bp and 16,606 bp, respectively. Organization and location of these genes were consistent with already studied Leuciscinae fishes. Synonymous codon usage was conservative in Leuciscinae as compared with other subfamilies in Cyprinidae. Phylogenetic analysis indicated that Leuciscinae was a monophyletic group, and genus Leuciscus was a paraphyletic group. Our approach, for the first time, of studying comparative mitochondrial genomics and phylogenetics together provided a supportive platform to the analysis of population genetics and phylogeny for Leuciscinae. Our results indicated a promising potential of comparative mitochondrial genomics in the manifestation of phylogenetic relationships between fishes, leading us to a suggestion that mitogenomes should be routinely considered in clarifying phylogenetics of family and subfamily members of fish.

The metazoan landscape of mitochondrial DNA gene order and content is shaped by selection and affects mitochondrial transcription

Mitochondrial DNA (mtDNA) harbors essential genes in most metazoans, yet the regulatory impact of the multiple evolutionary mtDNA rearrangements has been overlooked. Here, by analyzing mtDNAs from ~8000 metazoans we found high gene content conservation (especially of protein and rRNA genes), and codon preferences for mtDNA-encoded tRNAs across most metazoans. In contrast, mtDNA gene order (MGO) was selectively constrained within but not between phyla, yet certain gene stretches (ATP8-ATP6, ND4-ND4L) were highly conserved across metazoans. Since certain metazoans with different MGOs diverge in mtDNA transcription, we hypothesized that evolutionary mtDNA rearrangements affected mtDNA transcriptional patterns. As a first step to test this hypothesis, we analyzed available RNA-seq data from 53 metazoans. Since polycistron mtDNA transcripts constitute a small fraction of the steady-state RNA, we enriched for polycistronic boundaries by calculating RNA-seq read densities across junctions between gene couples encoded either by the same strand (SSJ) or by different strands (DSJ). We found that organisms whose mtDNA is organized in alternating reverse-strand/forward-strand gene blocks (mostly arthropods), displayed significantly reduced DSJ read counts, in contrast to organisms whose mtDNA genes are preferentially encoded by one strand (all chordates). Our findings suggest that mtDNA rearrangements are selectively constrained and likely impact mtDNA regulation.

Aonchotheca (Nematoda: Capillariidae) is validated as a separated genus from Capillaria by both mitochondrial and nuclear ribosomal DNA

BACKGROUND

The family Capillariidae is a group of thread-like nematodes of 27 genera and over 300 species that infect a great variety of hosts including humans. Among these, some taxa such as the genus Aonchotheca have remained controversial regarding their systematic status for decades. The aim of the current study was to verify Aonchotheca's systemic status and to further determine whether it is a distinct genus from Capillaria using molecular and phylogenetic analyses.

RESULTS

We sequenced the mitochondrial (mt) genome and nuclear small subunit (18S) rRNA gene of Aonchotheca putorii, a representative species of the genus, and investigated its systematic status in Trichinellida using maximum likelihood and Bayesian inference. The differences in amino acid sequences of 13 protein-coding genes were 12.69-67.35% among Aonchotheca, Capillaria, Eucoleus, and Pseudocapillaria with cox1 (12.69%) and atp8 (67.35%) as the most and the least conserved gene, respectively, and the difference of two mt rRNAs was 18.61-34.15%. Phylogenetic analyses of the complete mt genome and 18S rRNAs unequivocally showed that Aonchotheca was a distinct genus from Capillaria.

CONCLUSIONS

Large difference exists among Aonchotheca, Capillaria, Eucoleus, and Pseudocapillarias. Aonchotheca putorii is the first species in the genus Aonchotheca for which a complete mitogenome has been sequenced. These data are useful for phylogenetics, systematics and the evolution of Capillariidae.

Confirmation of the first intronic sequence in the bivalvian mitochondrial genome of Macoma balthica (Linnaeus, 1758)

In 2020, the first male-type mitochondrial genome from the clam Macoma balthica was published. Apart from the unusual doubly uniparental inheritance of mtDNA, scientists observed a unique (over 4k bp long) extension in the middle of the cox2 gene. We have attempted to replicate these data by NGS DNA sequencing and explore further the expression of the long cox2 gene. In our study, we report an even longer cox2 gene (over 5.5 kbp) with no stop codon separating conserved cox2 domains, as well as, based on the rtPCR, a lower relative gene expression pattern of the middle part of the gene (5' = 1; mid = 0.46; 3' = 0.89). Lastly, we sequenced the cox2 gene transcript proving the excision of the intronic sequence.

Characterization of the complete mitochondrial genome of the fluke of turdus, Plagiorchis elegans, and phylogenetic implications

Plagiorchis elegans (Trematoda: Digenea) is mainly parasitic in the intestines of vertebrate animals, including humans, causing irreversible pathological damage and herd-spherical influences. However, little information is available about its molecular epidemiology, population genetics, and phylogeny. In the present study, we sequenced, assembled, and annotated the complete mitochondrial (mt) genome of P. elegans. Combining with the available mitochondrial data of subclass Digenea, phylogenetic analysis was performed based on Bayesian inference (BI). The results showed that the complete length of P. elegans is 13,862 bp, including 12 PCGs, 2 ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), and one non-coding gene (NCR). There was an obvious A + T content from 61.0% to 71.3% and the values of the Ka/Ks ratio ranged from 0.119 (cox1) to 1.053 (nad6). In the BI analysis, different from previous studies, phylogenetic analysis showed genus Glypthelmins was paraphyletic rather than monophyletic and had a closer relationship with Plagiorchis and Orientocreadium. Additionally, the BI tree also presented that the genus Echinostoma was monophyletic. Our results provided molecular data in the family Plagiorchiidae proposing new insight within Xiphidiata and Echinostomata.

Complete Mitogenomes of Polypedates Tree Frogs Unveil Gene Rearrangement and Concerted Evolution within Rhacophoridae

New developments in sequencing technology and nucleotide analysis have allowed us to make great advances in reconstructing anuran phylogeny. As a clade of representative amphibians that have radiated from aquatic to arboreal habitats, our understanding of the systematic status and molecular biology of rhacophorid tree frogs is still limited. We determined two new mitogenomes for the genus Polypedates (Rhacophoridae): P. impresus and P. mutus. We conducted comparative and phylogenetic analyses using our data and seven other rhacophorid mitogenomes. The mitogenomes of the genera Polypedates, Buergeria, and Zhangixalus were almost identical, except that the ATP8 gene in Polypedates had become a non-coding region; Buergeria maintained the legacy "LTPF" tRNA gene cluster compared to the novel "TLPF" order in the other two genera; and B. buergeri and Z. dennysi had no control region (CR) duplication. The resulting phylogenetic relationship supporting the above gene rearrangement pathway suggested parallel evolution of ATP8 gene loss of function (LoF) in Polypedates and CR duplication with concerted evolution of paralogous CRs in rhacophorids. Finally, conflicting topologies in the phylograms of 185 species reflected the advantages of phylogenetic analyses using multiple loci.

Inverted base composition skews and discontinuous mitochondrial genome architecture evolution in the Enoplea (Nematoda)

Background

Within the class Enoplea, the earliest-branching lineages in the phylum Nematoda, the relatively highly conserved ancestral mitochondrial architecture of Trichinellida is in stark contrast to the rapidly evolving architecture of Dorylaimida and Mermithida. To better understand the evolution of mitogenomic architecture in this lineage, we sequenced the mitogenome of a fish parasite Pseudocapillaria tomentosa (Trichinellida: Capillariidae) and compared it to all available enoplean mitogenomes.

Results

P. tomentosa exhibited highly reduced noncoding regions (the largest was 98 bp), and a unique base composition among the Enoplea. We attributed the latter to the inverted GC skew (0.08) in comparison to the ancestral skew in Trichinellidae (-0.43 to -0.37). Capillariidae, Trichuridae and Longidoridae (Dorylaimida) generally exhibited low negative or low positive skews (-0.1 to 0.1), whereas Mermithidae exhibited fully inverted low skews (0 to 0.05). This is indicative of inversions in the strand replication order or otherwise disrupted replication mechanism in the lineages with reduced/inverted skews. Among the Trichinellida, Trichinellidae and Trichuridae have almost perfectly conserved architecture, whereas Capillariidae exhibit multiple rearrangements of tRNA genes. In contrast, Mermithidae (Mermithida) and Longidoridae (Dorylaimida) exhibit almost no similarity to the ancestral architecture.

Conclusions

Longidoridae exhibited more rearranged mitogenomic architecture than the hypervariable Mermithidae. Similar to the Chromadorea, the evolution of mitochondrial architecture in enoplean nematodes exhibits a strong discontinuity: lineages possessing a mostly conserved architecture over tens of millions of years are interspersed with lineages exhibiting architectural hypervariability. As Longidoridae also have some of the smallest metazoan mitochondrial genomes, they contradict the prediction that compact mitogenomes should be structurally stable. Lineages exhibiting inverted skews appear to represent the intermediate phase between the Trichinellidae (ancestral) and fully derived skews in Chromadorean mitogenomes (GC skews = 0.18 to 0.64). Multiple lines of evidence (CAT-GTR analysis in our study, a majority of previous mitogenomic results, and skew disruption scenarios) support the Dorylaimia split into two sister-clades: Dorylaimida + Mermithida and Trichinellida. However, skew inversions produce strong base composition biases, which can hamper phylogenetic and other evolutionary studies, so enoplean mitogenomes have to be used with utmost care in evolutionary studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08607-4.

Molecular detection of a novel Ancylostoma sp. by whole mtDNA sequence from pangolin Manis javanica

Background

Ancylostoma species are hematophagous parasites that cause chronic hemorrhage in various animals and humans. Pangolins, also known as scaly anteaters, are mammals that live in soil environments where they are readily exposed to soil-borne parasitic nematodes. However, only a limited number of helminth species have been identified in this animal host so far.

Methods

Ancylostoma sp. was isolated from a wild pangolin, and the complete mitochondrial (mt) genome of Ancylostoma sp. was obtained by Illumina sequencing of total genomic DNA.

Results

The circular complete mt genome that was assembled had a total length of 13,757 bp and comprised 12 protein-coding genes (PCGs), 22 transfer ribosomal RNAs, two ribosomal RNAs (rRNAs), two non-coding regions and one AT-rich region, but lacked the gene coding for ATPase subunit 8 (atp8). The overall AT content of the mt genome of Ancylostoma sp. was 76%, which is similar to that of other nematodes. The PCGs used two start codons (ATT and TTG) and three stop codons (TAA, TAG, and T). The nucleotide identity of the 12 PCGs ranged from 83.1% to 89.7% and had the highest sequence identity with Ancylostoma caninum among species in the Ancylostomatidae family. Also, the pangolin-derived Ancylostoma sp. lacked repeat sequences in the non-coding regions and in the unique sequence of the short non-coding regions, which differentiated it from other Ancylostoma species. In addition, phylogenetic analyses of 18S rRNA and mtDNA sequences revealed that the Ancylostoma sp. was positioned in a separate branch in the subfamily Ancylostomatinae along with other Ancylostoma species.

Conclusions

The Ancylostoma sp. isolated from a pangolin in this study was identified as a possible new Ancylostoma species. The identification of this Ancylostoma sp. from pangolin enriches our knowledge of the species in the Ancylostomatidae family and provides information that will lead to a better understanding of the taxonomy, diagnostics, and biology of hookworms.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05191-0.

Complete Mitogenome of Cruznema Tripartitum Confirms Highly Conserved Gene Arrangement within Family Rhabditidae

Mitochondrial genomes have widely been used as molecular markers in understanding the patterns and processes of nematode evolution. The species in genus Cruznema are free-living bacterivores as well as parasites of crickets and mollusks. The complete mitochondrial genome of C. tripartitum was determined through high-throughput sequencing as the first sequenced representative of the genus Cruznema. The genome is comprised of 14,067 bp nucleotides, and includes 12 protein-coding, two rRNA, and 22 tRNA genes. Phylogenetic analyses based on amino acid data support C. tripartitum as a sister to the clade containing Caenorhabditis elegans and Oscheius chongmingensis. The analysis of gene arrangement suggested that C. tripartitum shares the same gene order with O. chongmingensis, Litoditis marina, Diplocapter coronatus, genus Caenorhabditis, and Pristionchus pacificus. Thus, the mitochondrial gene arrangement is highly conserved in the family Rhabditidae as well as some species in Diplogasteridae.

Morphological and Molecular Characterization of Trichuris sp. (Nematoda: Trichuridae) in Crested Porcupines (Hystrix cristata; Rodentia: Hystricidae) from Italy

Adult specimens of Trichuris sp. collected from crested porcupines (Hystrix cristata) from Italy were characterized using an integrative taxonomic approach involving morphological and molecular tools. The morphological features of this Trichuris sp. were compared to data already available for Trichuris spp. from Hystrix sp., revealing diagnostic traits, such as spicule length in males or vulva shape in females, which distinguish this Trichuris sp. from the other species. Evidence from sequences analysis of the partial mitochondrial COX1 region indicated that the taxon under study is a distinct lineage. Biometrical and genetic data suggested this Trichuris sp. to be a valid and separated taxon. However, since molecular data from other Trichuris spp. infecting Hystrix, such as T. infundibulus, T. hystricis, T. javanica, T.landak and T. lenkorani, are missing in public repositories, the number and identity of distinct lineages able to infect porcupines remain only partially defined.

Sequencing of the Complete Mitochondrial Genome of Pingus sinensis (Spirurina: Quimperiidae): Gene Arrangements and Phylogenetic Implications

Despite several decades of intensive research on spirurine nematodes, molecular data on some of the main lineages are still absent, which makes taxonomic classification insufficiently resolved. In the present study, we sequenced the first complete mitogenome for the family Quimperiidae, belonging to P. sinensis (Spirurina: Quimperiidae), a parasite living in the intestines of snakehead (Ophiocephalus argus). The circular mitogenome is 13,874 bp long, and it contains the standard nematode gene set: 22 transfer RNAs, 2 ribosomal RNAs and 12 protein-coding genes. There are also two long non-coding regions (NCR), in addition to only 8 other intergenic regions, ranging in size from 1 to 58 bp. To investigate its phylogenetic position and study the relationships among other available Spirurina, we performed the phylogenetic analysis using Bayesian inference and maximum likelihood approaches by concatenating the nucleotide sequences of all 36 genes on a dataset containing all available mitogenomes of the suborder Spirurina from NCBI and compared with gene order phylogenies using the MLGO program. Both supported the closer relationship of Ascaridoidea to Seuratoidea than to Spiruroidea. Pingus formed a sister-group with the Cucullanus genus. The results provide a new insights into the relationships within Spirurina.

Infection, genetics, and evolution of Trichinella: Historical insights and applications to molecular epidemiology

Genetic variation in pathogen populations provides the means to answer questions in disease ecology and transmission, illuminating interactions between genetic traits, environmental exposures, and disease. Such studies elucidate the phylogeny, evolution, transmission and pathogenesis of viruses, bacteria and parasites. Here, we review how such studies have fostered understanding of the biology and epidemiology of zoonotic nematode parasites in the genus Trichinella spp., which impose considerable economic and health burdens by infecting wildlife, livestock, and people. To use such data to define ongoing chains of local transmission and source traceback, researchers first must understand the extent and distribution of genetic variation resident in regional parasite populations. Thus, genetic variability illuminates a population's past as well as its present. Here we review how such data have helped define population dynamics of Trichinella spp. in wild and domesticated hosts, creating opportunities to harness genetic variation in the quest to prevent, track, and contain future outbreaks.

The mitochondrial genome sequence analysis of Ophidascaris baylisi from the Burmese python (Python molurus bivittatus)

Ophidascaris species are parasitic roundworms that inhabit the python gut, resulting in severe granulomatous lesions or even death. However, the classification and nomenclature of these roundworms are still controversial. Our study aims to identify a snake roundworm from the Burmese python (Python molurus bivittatus) and analyze the mitochondrial genome. We identified this roundworm as Ophidascaris baylisi based on the morphology and cytochrome c oxidase subunit I (cox1) sequence. Ophidascaris baylisi complete mitochondrial genome was 14,784 bp in length, consisting of two non-coding regions and 36 mitochondrial genes (12 protein-coding genes, 22 tRNA genes, and two rRNA genes). The protein-coding genes used TTG, ATG, ATT, or TTA as start codons and TAG, TAA, or T as stop codons. All tRNA genes showed a TV-loop structure, except trnS1^AGN and trnS2^UCN revealed a D-loop structure. The mitochondrial large ribosomal subunit 16S (rrnL) and small ribosomal subunit 12S (rrnS) were 956 bp and 700 bp long, respectively. Phylogenetic analysis based on O. baylisi mitochondrial protein-coding genes demonstrated that O. baylisi clustered with the family Ascarididae members and was most closely related to Ophidascaris wangi. These results may enhance the nematode mitochondrial genome database and provide valuable molecular markers for further research on the taxonomy, phylogeny, and genetic relationships of Ophidascaris nematodes.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

The complete mitochondrial genome of capillariid nematodes (Eucoleus annulatus): A novel gene arrangement and phylogenetic implications

Capillariid nematode is a group of endoparasites of vertebrates with a complex taxonomy, causing significant economic losses to poultry industry. The taxonomic status of the genus Eucoleus remained controversial for several decades. Mitochondrial (mt) DNA provides useful genetic markers for accurate identification of species, but complete mt genome sequences have been lacking for any Capillariid nematodes. In the present study, we decoded the complete mt genome of E. annulatus and examined its phylogenetic relationship with selected members of the class Enoplea nematodes. The circular mt genome of E. annulatus was 14,118 bp, encoded 37 genes with a single non-coding region and showed substantial gene rearrangements (especially tRNA genes) compared to other nematodes studied to date. The complete mt genome of E. annulatus showed a clear A + T bias in nucleotide composition. The number of A (5404) was approximately equal to T (5405) and the GC-skew was negative on average (-0.073). Phylogenetic analyses based on 18S rDNA placed Eucoleus spp. well apart from each other and supported the proposal that Eucoleus and Capillaria are two distinct genera. Similarly, Bayesian inference (BI) and Maximum likelihood (ML) phylogenies based on mtDNA sequences revealed that the family Capillariidae is more closely related to the family Trichuridae than to the family Trichinellidae. This is the first report of the complete mt genome of capillariid nematodes, and it will provide additional genetic markers for studying the molecular epidemiology, population genetics and systematics of capillariid nematodes and should have implications for the molecular diagnosis, prevention, and control of capillariosis in animals.

How are the mitochondrial genomes reorganized in Hexapoda? Differential evolution and the first report of convergences within Hexapoda

The evolution of the Hexapoda mitochondrial genome has been the focus of several genetic and evolutionary studies over the last decades. However, they have concentrated on certain taxonomic orders of economic or health importance. The recent increase of mitochondrial genomes sequencing of diverse taxonomic orders generates an important opportunity to clarify the evolution of this group of organisms. However, there is no comparative study that investigates the evolution of the Hexapoda mitochondrial genome. In order to verify the level of rearrangement and the mitochondrial genome evolution, we performed a comparative genomic analysis of the Hexapoda mitochondrial genome available in the NCBI database. Using a combination of bioinformatics methods to carefully examine the mitochondrial gene rearrangements in 1198 Hexapoda species belonging to 32 taxonomic orders, we determined that there is a great variation in the rate of rearrangement by gene and by taxonomic order. A higher rate of genetic reassortment is observed in Phthiraptera, Thysanoptera, Protura, and Hymenoptera; compared to other taxonomic orders. Twenty-four events of convergence in the genetic order between different taxonomic orders were determined, most of them not previously reported; which proves the great evolutionary dynamics within Hexapoda.

Complete Mitochondrial Genome of Trichuristrichiura from Macaca sylvanus and Papio papio

Trichuriasis is among the most prevalent worldwide parasitism caused by helminths. For many years, Trichuris spp. have been described with a relatively narrow range of both morphological and biometrical features. The use of the complete mitochondrial genome (mitogenome) is an alternative and powerful molecular method for inferring phylogenies. Here, we present an overview of the contributions of mitogenome for Trichuris spp. from human and non-human primates. In addition, we carry out structural and phylogenetic comparative analyses with genomes of Trichuris species available in public datasets. The complete mt genomes of Trichuris trichiura and Trichuris sp. from Macaca sylvanus and T. trichiura from Papio papio are 14,091 bp, 14,047 bp and 14,089 bp in length, respectively. The three mt genomes are circular and consist of 37 genes—13 PCGs (cox1–3, nad1–6, nad4L, atp6, atp8 and cob), 22 transfer RNA genes (tRNAs), and two rRNAs (rrnL and rrnS). The molecular evidence presented here supports the hypothesis that T. trichiura de M. sylvanus (TMF31) and T. trichiura de P. papio (TPM1) were similar but genetically different with respect to Trichuris sp. from macaques (TMM5). The phylogenetic study also supported the evolution of the different Trichuris species. In conclusion, we suggest the existence of two cryptic species parasitizing M. sylvanus.

Complete mitochondrial genomes of Epeorus carinatus and E. dayongensis (Ephemeroptera: Heptageniidae): Genomic comparison and phylogenetic inference

The current research on Ephemeroptera is mainly based on its morphology, since only small numbers of mitogenomes have been reported. In this study, the mitogenomes of Epeorus carinatus (15,338 bp) and E. dayongensis (15,609 bp) were sequenced, annotated and compared to genome data from congeners. Both mitogenomes had 23 tRNA genes including standard 22 and one extra tRNA^Met. The duplicated tRNA^Met gene had been found in other heptageniid species except Paegniodes cupulatus, suggesting it could be used as a molecular synapomorphy for partial Heptageniidae. The phylogenetic analyses based on Bayesian Inference (BI) and Maximum Likelihood (ML) showed that Heptageniidae was monophyletic and the relationships among known Epeorus species were ((E. carinatus + E. herklotsi) + (E. dayongensis + E. sp. 1)), which implied the focal species E. carinatus and E. dayongensis should be grouped into different subgenera.

Comparative mitogenome analyses uncover mitogenome features and phylogenetic implications of the subfamily Cobitinae

BACKGROUND

Loaches of Cobitinae, widely distributed in Eurasian continent, have high economic, ornamental and scientific value. However, the phylogeny of Cobitinae fishes within genera or family level remains complex and controversial. Up to now, about 60 Cobitinae mitogenomes had been deposited in GenBank, but their integrated characteristics were not elaborated.

RESULTS

In this study, we sequenced and analyzed the complete mitogenomes of a female Cobits macrostigma. Then we conducted a comparative mitogenome analysis and revealed the conserved and unique characteristics of 58 Cobitinae mitogenomes, including C. macrostigma. Cobitinae mitogenomes display highly conserved tRNA secondary structure, overlaps and non-coding intergenic spacers. In addition, distinct base compositions were observed among different genus and significantly negative linear correlation between AT% and AT-skew were found among Cobitinae, genus Cobitis and Pangio mitogenomes, respectively. A specific 3 bp insertion (GCA) in the atp8-atp6 overlap was identified as a unique feature of loaches, compared to other Cypriniformes fish. Additionally, all protein coding genes underwent a strong purifying selection. Phylogenetic analysis strongly supported the paraphyly of Cobitis and polyphyly of Misgurnus. The strict molecular clock predicted that Cobitinae might have split into northern and southern lineages in the late Eocene (42.11 Ma), furthermore, mtDNA introgression might occur (14.40 Ma) between ancestral species of Cobitis and ancestral species of Misgurnus.

CONCLUSIONS

The current study represents the first comparative mitogenomic and phylogenetic analyses within Cobitinae and provides new insights into the mitogenome features and evolution of fishes belonging to the cobitinae family.

Nuclear and Mitochondrial Data on Trichuris from Macaca fuscata Support Evidence of Host Specificity

Whipworms are parasitic intestinal nematodes infecting mammals, and traditionally humans and other primates that have so far been considered infected by Trichuris trichiura. Recent molecular studies report a more complex scenario suggesting the presence of a species complex with several Trichuris taxa specifically infecting only one primate species as well as taxa able to infect a range of primate species. The systematics of the group is important for taxonomic inference, to estimate the relative zoonotic potential, and for conservation purposes. In fact, captive animals living in zoological gardens are usually infected by persistent monoxenous intestinal parasites. Here, two Japanese macaques living in the Bioparco Zoological Garden of Rome were found infected by Trichuris sp. Nematodes were characterized at the molecular level using nuclear (btub and 18S) and mitochondrial (16S and cytb) markers and then compared to Trichuris collected previously in the same location, and to other Trichuris infecting primates. Evidences from mitochondrial and nuclear markers allowed for the identification of Trichuris sp. specific to Macaca fuscata. Results obtained here also described a uniform taxonomic unit of Trichuris, separated but closely related to Trichuris trichiura, thus, emphasizing its zoonotic potential for workers and visitors.

The revised complete mitogenome sequence of the tree frog Polypedatesmegacephalus (Anura, Rhacophoridae) by next-generation sequencing and phylogenetic analysis

The mitochondrial genome (mitogenome) sequence of the tree frog Polypedates megacephalus (16,473 bp) was previously reported as having the unusual characteristic of lacking the ND5 gene. In this study, a new mitogenome of P. megacephalus (19,952 bp) was resequenced using the next-generation sequencing (NGS) and standard Sanger sequencing technologies. It was discovered that the ND5 gene was not lost but translocated to the control region (CR) from its canonical location between the ND4 and ND6 genes. In addition, a duplicated control region was found in the new mitogenome of this species. Conservative region identification of the ND5 gene and phylogenetic analysis confirmed that the ND5 gene was located between two control regions. The phylogenetic relationship among 20 related species of anura revealed a rearrangement of the ND5 gene during the evolutionary process. These results also highlighted the advantages of next-generation sequencing. It will not only decrease the time and cost of sequencing, but also will eliminate the errors in published mitogenome databases.

The complete mitochondrial genome of Echinostoma miyagawai: Comparisons with closely related species and phylogenetic implications

Echinostoma miyagawai (Trematoda: Echinostomatidae) is a common parasite of poultry that also infects humans. Es. miyagawai belongs to the "37 collar-spined" or "revolutum" group, which is very difficult to identify and classify based only on morphological characters. Molecular techniques can resolve this problem. The present study, for the first time, determined, and presented the complete Es. miyagawai mitochondrial genome. A comparative analysis of closely related species, and a reconstruction of Echinostomatidae phylogeny among the trematodes, is also presented. The Es. miyagawai mitochondrial genome is 14,416 bp in size, and contains 12 protein-coding genes (cox1-3, nad1-6, nad4L, cytb, and atp6), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and one non-coding region (NCR). All Es. miyagawai genes are transcribed in the same direction, and gene arrangement in Es. miyagawai is identical to six other Echinostomatidae and Echinochasmidae species. The complete Es. miyagawai mitochondrial genome A + T content is 65.3%, and full-length, pair-wise nucleotide sequence identity between the six species within the two families range from 64.2-84.6%. The Es. miyagawai sequences is most similar to Echinostoma caproni. Sequence difference are 15.0-33.5% at the nucleotide level, and 8.6-44.2% at the amino acid level, among the six species, for the 12 protein-coding genes. ATG and TAG are the most common initiation and termination codons, respectively. Twenty of the Es. miyagawai transfer RNA genes transcribe products of the conventional cloverleaf structure, while two of the transfer RNA genes, namely trnS1^(AGC) and trnS2^(UGA), have unpaired D-arms. Phylogenetic analyses using our mitochondrial data indicate that Es. miyagawai is closely related to other Echinostomatidae species, except for Echinostoma hortense, which forms a distinct paraphyletic branch, and Echinochasmus japonicus, which is outside the clade containing all other Echinostomatidae species. These phylogenetic results support the elevation of subfamily Echinostomatidae. Our dataset also provides a significant resource of molecular markers to study the taxonomy, population genetics, and systematics of the echinostomatids.

Characterization of the Complete Mitochondrial Genome of a Whipworm Trichuris skrjabini (Nematoda: Trichuridae)

The complete mitochondrial (mt) genome of Trichuris skrjabini has been determined in the current study and subsequently compared with closely related species by phylogenetic analysis based on concatenated datasets of mt amino acid sequences. The whole mt genome of T. skrjabini is circular and 14,011 bp in length. It consists of a total of 37 genes including 13 protein coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNAs) genes, and two non-coding regions. The gene arrangement and contents were consistent with other members of the Trichuridae family including Trichuris suis, Trichuris trichiura, Trichuris ovis, and Trichuris discolor. Phylogenetic analysis based on concatenated datasets of amino acids of the 12 PCGs predicted the distinctiveness of Trichuris skrjabini as compared to other members of the Trichuridae family. Overall, our study supports the hypothesis that T. skrjabini is a distinct species. The provision of molecular data of whole mt genome of T. skrjabini delivers novel genetic markers for future studies of diagnostics, systematics, population genetics, and molecular epidemiology of T. skrjabini.

Unusual AT-skew of Sinorhodeus microlepis mitogenome provides new insights into mitogenome features and phylogenetic implications of bitterling fishes

Sinorhodeus microlepis (S. microlepis) is recently described as a new species and represents a new genus Sinorhodeu of the subfamily Acheilognathinae. In this study, we first sequenced the complete mitogenome of S. microlepis and compared with the other 29 bitterling mitogenomes. The S. microlepis mitogenome is 16,591 bp in length and contains 37 genes. Gene distribution pattern is identical among 30 bitterling mitogenomes. A significant linear correlation between A+T% and AT-skew were found among 29 bitterling mitogenomes, except S. microlepis shows unusual AT-skew with slightly negative in tRNAs and PCGs. Bitterling mitogenomes exhibit highly conserved usage bias of start codon, relative synonymous codons and amino acids, overlaps and non-coding intergenic spacers. Phylogenetic trees constructed by 13 PCGs strongly support the polyphyly of the genus Acheilognathus and the paraphyly of Rhodeus and Tanakia. Together with the unusual characters of S. microlepis mitogenomes and phylogenetic trees, S. microlepis should be a sister species to the genus Rhodeu that might diverge about 13.69 Ma (95% HPD: 12.96-14.48 Ma).

High genetic diversity in the Dirofilaria repens species complex revealed by mitochondrial genomes of feline microfilaria samples from Narathiwat, Thailand

Dirofilaria repens is a zoonotic, mosquito-borne filaria infecting carnivores, particularly dogs. It is expanding its range in Europe but epidemiological information is sparse for other Eurasian regions. In Hong Kong and India, the closely related species Candidatus Dirofilaria hongkongensis was proposed. Previous analysis of 2.5 kb partial mitochondrial genome sequences containing the particularly variable non-coding control region revealed low diversity in European D. repens while Asian nematodes showed high diversity. Sequences derived from feline blood samples from Narathiwat (Thailand) led to the proposal of a third potential species, Dirofilaria sp. "Thailand II". To avoid bias from rapidly evolving non-coding regions, this study aimed to compare Dirofilaria sp. "Thailand II" with D. repens and C. D. hongkongensis based on complete mitochondrial genomes. Using PCRs and Sanger sequencing, three complete mitochondrial genomes (13,651 bp) were assembled from DNA obtained from different feline blood samples. Mitochondrial genome organization was identical to other onchocercids with eleven protein-coding, two rRNA and 22 tRNA genes and no atp-8 gene. All genes were on the same strand showing an extremely high thymidine content (56.7%). Maximum-likelihood phylogenetic analysis using protein and rRNA sequences confirmed closer relationship of Dirofilaria sp. "Thailand II" to C. D. hongkongensis than to D. repens. All distances between these three putative species were considerably larger than the distance between the valid sibling species Onchocerca volvulus and Onchocerca ochengi. Sequencing of a 2.5 kb fragment containing the control region from microfilarial DNA from additional feline blood samples from Narathiwat 3-4 years later revealed that these also fell into the C. D. hongkongensis clade but were remarkably different from C. D. hongkongensis and Dirofilaria sp. "Thailand II". Since D. repens-like filaria are absent from dogs in Narathiwat, further field studies are required to confirm if these genotypes represent locally circulating cat-specific Dirofilaria genotypes or species.

Genetic characterisation and phylogenetic status of whipworms (Trichuris spp.) from captive non-human primates in China, determined by nuclear and mitochondrial sequencing

Background

Whipworms (Nematoda: Trichuridae), among the most common soil-transmitted helminths (STHs), can cause the socioeconomically important disease trichuriasis in various mammalian hosts including humans and non-human primates. For many years, Trichuris from non-human primates has been assigned to the same species as the one infecting humans Trichuris trichiura. More recently, several molecular reports challenged this assumption following recognition of a Trichuris species complex observed in humans and non-human primates. A refined concept for species limits within Trichuris contributes to an understanding of diversity and the potential (zoonotic) transmission among humans and non-human primates. In this study, we expanded previous investigations by exploring the diversity of Trichuris among eight primates including three Asian autochthonous species (i.e. Rhinopithecus roxellana, Rhinopithecus bieti and Nomascus leucogenys). Species-level identification, whether novel or assignable to known lineages of Trichuris, was based on analyses of nuclear internal transcribed spacers (ITS) and mitochondrial cytochrome c oxidase subunit 1 (cox1) genes.

Results

In total, seven genetically distinct subgroups of whipworms were determined to be present among the primates sampled. Most Trichuris lineages, including Subgroups 1, 1’, 3, 5 and 6, showed a broad host range and were not restricted to particular primate species; in addition to T. trichiura, a complex of Trichuris species was shown infecting primates. Furthermore, it was assumed that Trichuris spp. from either N. leucogenys and P. hamadryas or R. roxellana and R. bieti, respectively, were conspecific. Each pair was indicated to be a discrete lineage of Trichuris, designated, respectively, as Subgroups 1 or 1’ and 2, based on integrated genetic and phylogenetic evidence.

Conclusion

These results emphasise that the taxonomy and genetic variations of Trichuris are more complicated than previously acknowledged. These cumulative molecular and phylogenetic data provide a better understanding of the taxonomy, genetics and evolutionary biology of the whipworms.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3100-5) contains supplementary material, which is available to authorized users.

Complete mitochodrial genome of the crab spider Ebrechtella tricuspidata (Araneae: Thomisidae): A novel tRNA rearrangement and phylogenetic implications for Araneae

Mitochondrial genomes are widely used for phylogenetic and phylogeographic analyses among arthropods, but there is a lack of sufficient mitochondrial genome sequence data for spiders. Herein, we sequenced and characterized the complete mitochondrial genome of a crab spider Ebrechtella tricuspidata (Araneae: Thomisidae). The circular mitochondrial genome is 14,352 bp long, including a standard set of 37 genes and an A + T-rich region. Nucleotide composition is highly biased toward A + T nucleotides (77.3%). A novel gene order rearrangement is detected by a tRNA (trnL1) translocation. Tandem repeats are not identified in the A + T-rich region. Most of the tRNAs are greatly reduced in size and cannot be folded into typical cloverleaf-shaped secondary structures. The phylogenetic analysis confirms that the mitochondrial genome sequences are useful in resolving higher-level relationship of Araneae. Overall, our data present in this study will elevate our knowledge on the architecture and evolution of spider mitochondrial genome.

Comparative analyses of the complete mitochondrial genomes of Cyathostomum pateratum and Cyathostomum catinatum provide new molecular data for the evolution of Cyathostominae nematodes

The parasite Cyathostomum pateratum, which occurs in the large intestine of equines, is a common species of the subfamily Cyathostominae. Cyathostominae nematodes are a complex nematode group for which only limited genetic information has been reported. To re-examine the phylogenetic relationships among Cyathostominae nematodes, we sequenced the complete mitochondrial (mt) genome of Cy. pateratum and compared it with the mt genome of the congeneric species Cyathostomum catinatum. The complete mtDNA sequence of Cy. pateratum was 13,822 bp in length, 16 bp shorter than that of Cy. catinatum. The mtDNA sequences of both species contained 12 protein-coding genes, two rRNA genes and 22 tRNA genes, and all 36 genes were transcribed in the same direction and in the same strand. Pairwise comparisons of the 12 predicted amino acid sequences between Cy. catinatum and Cy. pateratum revealed differences of 0.4-3.1%; the least conserved sequence was that of cytochrome c oxidase subunit 3 (cox3). Phylogenetic analysis of the concatenated amino acid sequences using Bayesian inference and maximum likelihood methods showed that Cy. catinatum and Cy. pateratum clustered together with very high nodal support, and Cylicostephanus goldi was closer to the Cyathostomum nematodes than to other Cyathostominae nematodes. The mtDNA sequence of Cy. pateratum is reported here for the first time. The study will shed some light on the genetic evolution among parasitic nematodes in Cyathostomum.

Actively transcribed and expressed atp8 gene in Mytilus edulis mussels

Background

Animal mitochondrial genomes typically encode 37 genes: 13 proteins, 22 tRNAs and two rRNAs. However, many species represent exceptions to that rule. Bivalvia along with Nematoda and Platyhelminthes are often suspected to fully or partially lack the ATP synthase subunit 8 (atp8) gene. This raises the question as to whether they are really lacking this gene or is this maybe an annotation problem? Among bivalves, Mytilus edulis has been inferred to lack an ATP8 gene since the characterization of its mitochondrial genome in 1992. Even though recent bioinformatic analyses suggested that atp8 is present in Mytilus spp., due to high divergence in predicted amino acid sequences, the existence of a functional atp8 gene in this group remains controversial.

Results

Here we demonstrate that M. edulis mitochondrial open reading frames suggested to be atp8 (in male and female mtDNAs) are actively translated proteins. We also provide evidence that both proteins are an integral part of the ATP synthase complex based on in-gel detection of ATP synthase activity and two-dimensional Blue-Native and SDS polyacrylamide electrophoresis.

Conclusion

Many organisms (e.g., Bivalvia along with Nematoda and Platyhelminthes) are considered to be lacking certain mitochondrial genes often only based on poor similarity between protein coding gene sequences in genetically closed species. In some situations, this may lead to the inference that the ATP8 gene is absent, when it is in fact present, but highly divergent. This shows how important complementary role protein-based approaches, such as those in the present study, can provide to bioinformatic, genomic studies (i.e., ability to confirm the presence of a gene).

The bipartite mitochondrial genome of Ruizia karukerae (Rhigonematomorpha, Nematoda)

Mitochondrial genes and whole mitochondrial genome sequences are widely used as molecular markers in studying population genetics and resolving both deep and shallow nodes in phylogenetics. In animals the mitochondrial genome is generally composed of a single chromosome, but mystifying exceptions sometimes occur. We determined the complete mitochondrial genome of the millipede-parasitic nematode Ruizia karukerae and found its mitochondrial genome consists of two circular chromosomes, which is highly unusual in bilateral animals. Chromosome I is 7,659 bp and includes six protein-coding genes, two rRNA genes and nine tRNA genes. Chromosome II comprises 7,647 bp, with seven protein-coding genes and 16 tRNA genes. Interestingly, both chromosomes share a 1,010 bp sequence containing duplicate copies of cox2 and three tRNA genes (trnD, trnG and trnH), and the nucleotide sequences between the duplicated homologous gene copies are nearly identical, suggesting a possible recent genesis for this bipartite mitochondrial genome. Given that little is known about the formation, maintenance or evolution of abnormal mitochondrial genome structures, R. karukerae mtDNA may provide an important early glimpse into this process.

Full mitochondrial and nuclear genome comparison confirms that Onchocerca sp. "Siisa" is Onchocerca ochengi

Onchocerca ochengi is a nodule-forming filarial nematode parasite of cattle. It is the closest known relative of the human parasite Onchocerca volvulus, with which it shares the black fly vector Simulium damnosum. Onchocerca sp. "Siisa" was described in black flies and in cattle and, based on limited mitochondrial sequence information, appeared to be about equally phylogenetically distant from O. ochengi and O. volvulus. Based on molecular genetic markers and apparent interbreeding, we later proposed that O. sp. "Siisa" belongs to the species O. ochengi. However, we did not demonstrate directly that the hybrids were fertile, and we were still unable to resolve the phylogenetic relationship of O. ochengi, O. sp. "Siisa," and O. volvulus, leaving some concerns with the conclusion mentioned above. Here, we present fully assembled, manually curated mitochondrial genomes of O. ochengi and O. sp. "Siisa," and we compare multiple individuals of these two taxa with respect to their whole mitochondrial and nuclear genomes. Based on the mitochondrial genomes, O. ochengi and O. sp. "Siisa" are phylogenetically much closer to each other than to O. volvulus. The differences between them are well within the range of what is expected for within-species variation. The nuclear genome comparison provided no indication of genetic separation of O. ochengi and O. sp. "Siisa." From this, in combination with the earlier literature, we conclude that O. ochengi and O. sp. "Siisa" should be considered one species.

Morphological-molecular characterization and phylogenetic relationships of a new Trichuris species (Nematoda: Trichuridae) parasitic on Holochilus chacarius (Cricetidae: Sigmodontinae) from the Chaco ecoregion (Argentina)

A new Trichuris species isolated from Holochilus chacarius (Cricetidae: Sigmodontinae) from the Chaco ecoregion in Argentina is described based on morphological characteristics and mitochondrial (cox1, cob) and nuclear (ITS2) markers. The new species is distinguished from 27 species of Trichuris from North and South American rodents based on morphological and biometrical features, such as the absence of a spicular tube, presence of a cylindrical spicular sheath, non-protrusive vulva, length of spicule, and proximal and distal cloacal tube. In addition, the results based on three molecular markers of the new species and three Trichuris species previously analysed from sigmodontine rodents: Trichuris pardinasi, Trichuris navonae and Trichuris bainae confirmed that the specimens here studied belong to a different species. Molecular data are further used to discuss the phylogenetic relationships among the Trichuris species of rodents from Argentina. The combined analysis of mitochondrial genes (cox1 and cob) revealed four clades corresponding with four different species of Trichuris. T. navonae (Akodontini rodents) and T. massoiai n. sp. as a sister group related to T. bainae (Oryzomyini rodents) and separated of T. pardinasi (Phyllotini rodents).

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class

BACKGROUND

Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes.

RESULTS

The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder.

CONCLUSIONS

Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.

Mitochondrial DNA Evidence Supports the Hypothesis that Triodontophorus Species Belong to Cyathostominae

Equine strongyles, the significant nematode pathogens of horses, are characterized by high quantities and species abundance, but classification of this group of parasitic nematodes is debated. Mitochondrial (mt) genome DNA data are often used to address classification controversies. Thus, the objectives of this study were to determine the complete mt genomes of three Cyathostominae nematode species (Cyathostomum catinatum, Cylicostephanus minutus, and Poteriostomum imparidentatum) of horses and reconstruct the phylogenetic relationship of Strongylidae with other nematodes in Strongyloidea to test the hypothesis that Triodontophorus spp. belong to Cyathostominae using the mt genomes. The mt genomes of Cy. catinatum, Cs. minutus, and P. imparidentatum were 13,838, 13,826, and 13,817 bp in length, respectively. Complete mt nucleotide sequence comparison of all Strongylidae nematodes revealed that sequence identity ranged from 77.8 to 91.6%. The mt genome sequences of Triodontophorus species had relatively high identity with Cyathostominae nematodes, rather than Strongylus species of the same subfamily (Strongylinae). Comparative analyses of mt genome organization for Strongyloidea nematodes sequenced to date revealed that members of this superfamily possess identical gene arrangements. Phylogenetic analyses using mtDNA data indicated that the Triodontophorus species clustered with Cyathostominae species instead of Strongylus species. The present study first determined the complete mt genome sequences of Cy. catinatum, Cs. minutus, and P. imparidentatum, which will provide novel genetic markers for further studies of Strongylidae taxonomy, population genetics, and systematics. Importantly, sequence comparison and phylogenetic analyses based on mtDNA sequences supported the hypothesis that Triodontophorus belongs to Cyathostominae.

Characterization of the complete mitochondrial genome of Ortleppascaris sinensis (Nematoda: Heterocheilidae) and comparative mitogenomic analysis of eighteen Ascaridida nematodes

Ortleppascaris sinensis (Nematoda: Ascaridida) is a dominant intestinal nematode of the captive Chinese alligator. However, the epidemiology, molecular ecology and population genetics of this parasite remain largely unexplored. In this study, the complete mitochondrial (mt) genome sequence of O. sinensis was first determined using a polymerase chain reaction (PCR)-based primer-walking strategy, and this is also the first sequencing of the complete mitochondrial genome of a member of the genus Ortleppascaris. The circular mitochondrial genome (13,828 bp) of O. sinensis contained 12 protein-coding, 22 transfer RNA and 2 ribosomal RNA genes, but lacked the ATP synthetase subunit 8 gene. Finally, phylogenetic analysis of mtDNAs indicated that the genus Ortleppascaris should be attributed to the family Heterocheilidae. It is necessary to sequence more mtNDAs of Ortleppascaris nematodes in the future to test and confirm our conclusion. The complete mitochondrial genome sequence of O. sinensis reported here should contribute to molecular diagnosis, epidemiological investigations and ecological studies of O. sinensis and other related Ascaridida nematodes.

Duplication of Drosophila melanogaster mitochondrial EF-Tu: pre-adaptation to T-arm truncation and exclusion of bulky aminoacyl residues

Translation elongation factor Tu (EF-Tu) delivers aminoacyl-tRNA (aa-tRNA) to ribosomes in protein synthesis. EF-Tu generally recognizes aminoacyl moieties and acceptor- and T-stems of aa-tRNAs. However, nematode mitochondrial (mt) tRNAs frequently lack all or part of the T-arm that is recognized by canonical EF-Tu. We previously reported that two distinct EF-Tu species, EF-Tu1 and EF-Tu2, respectively, recognize mt tRNAs lacking T-arms and D-arms in the mitochondria of the chromadorean nematode Caenorhabditis elegansC. elegans EF-Tu2 specifically recognizes the seryl moiety of serylated D-armless tRNAs. Mitochondria of the enoplean nematode Trichinella possess three structural types of tRNAs: T-armless tRNAs, D-armless tRNAs, and cloverleaf tRNAs with a short T-arm. Trichinella mt EF-Tu1 binds to all three types and EF-Tu2 binds only to D-armless Ser-tRNAs, showing an evolutionary intermediate state from canonical EF-Tu to chromadorean nematode (e.g. C. elegans) EF-Tu species. We report here that two EF-Tu species also participate in Drosophila melanogaster mitochondria. Both D. melanogaster EF-Tu1 and EF-Tu2 bound to cloverleaf and D-armless tRNAs. D. melanogaster EF-Tu1 has the ability to recognize T-armless tRNAs that do not evidently exist in D. melanogaster mitochondria, but do exist in related arthropod species. In addition, D. melanogaster EF-Tu2 preferentially bound to aa-tRNAs carrying small amino acids, but not to aa-tRNAs carrying bulky amino acids. These results suggest that the Drosophila mt translation system could be another intermediate state between the canonical and nematode mitochondria-type translation systems.

Mitochondrial genome diversity in dagger and needle nematodes (Nematoda: Longidoridae)

Dagger and needle nematodes included in the family Longidoridae (viz. Longidorus, Paralongidorus, and Xiphinema) are highly polyphagous plant-parasitic nematodes in wild and cultivated plants and some of them are plant-virus vectors (nepovirus). The mitochondrial (mt) genomes of the dagger and needle nematodes, Xiphinema rivesi, Xiphinema pachtaicum, Longidorus vineacola and Paralongidorus litoralis were sequenced in this study. The four circular mt genomes have an estimated size of 12.6, 12.5, 13.5 and 12.7 kb, respectively. Up to date, the mt genome of X. pachtaicum is the smallest genome found in Nematoda. The four mt genomes contain 12 protein-coding genes (viz. cox1-3, nad1-6, nad4L, atp6 and cob) and two ribosomal RNA genes (rrnL and rrnS), but the atp8 gene was not detected. These mt genomes showed a gene arrangement very different within the Longidoridae species sequenced, with the exception of very closely related species (X. americanum and X. rivesi). The sizes of non-coding regions in the Longidoridae nematodes were very small and were present in a few places in the mt genome. Phylogenetic analysis of all coding genes showed a closer relationship between Longidorus and Paralongidorus and different phylogenetic possibilities for the three Xiphinema species.