Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection

BACKGROUND

Nematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures.

RESULTS

Generally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection.

CONCLUSIONS

Nematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.

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.

Evolution of a biological thermocouple by adaptation of cytochrome c oxidase in a subterrestrial metazoan

In this study we report a naturally evolved temperature-sensing electrical regulator in the cytochrome c oxidase of the Devil Worm, Halicephalobus mephisto. This extremophile metazoan was isolated 1.3 km underground in a South African goldmine, where it adapted to heat and potentially to hypoxia, making its mitochondrial sequence a likely target of adaptational change. We obtained the full mitochondrial genome sequence of this organism, and show through dN/dS analysis statistically robust evidence of positive selection in H. mephisto cytochrome c oxidase subunits. Seventeen of these positively-selected amino acid substitutions were localized in proximity to the H- and K-pathway proton channels of the complex. Surprisingly, the H. mephisto cytochrome c oxidase proton pump completely shuts down at low temperatures (20°C) leading to approximately a 4.8-fold reduction in the transmembrane proton gradient voltage (ΔΨm) compared to optimal temperature (37°C). Direct measurement of oxygen consumption found a corresponding 4.7-fold drop at 20°C compared to 37°C. Correspondingly, the lifecycle of H. mephisto takes four-fold longer at the low temperature compared to higher. This elegant evolutionary adaptation creates a finely-tuned mitochondrial temperature sensor, allowing this ectothermic organism to maximize its reproductive success in varying environmental temperatures. Our study shows that evolutionary innovation may remodel core metabolism to make it more accurately map onto environmental variation.

Complete mitochondrial genome of Metathelazia capsulata (Pneumospiruridae) and comparison with other Spiruromorpha species

Metathelazia capsulata (family Pneumospiruridae) is a lungworm parasitizing the bronchi and bronchioles, described in four species of wild carnivores. Very little molecular data are available on this nematode and none on other species of the Pneumospiruridae family. In this work, we describe for the first time the complete mitogenome (mitochondrial genome) of M. capsulata, being the first described of the family Pneumospiruridae. The mitogenome of M. capsulata has 13,659 bp in length, an A + T content of 79.2%. The mitogenome included 12 protein-coding genes (PCGs) (lacking the atp8 gene), 22 tRNA genes, 2 rRNA genes (all the genes are coded by the heavy strand), and an AT-rich region. The PCGs varied in size (232 bp-1645 bp). Only the tRNA-Trp has the standard cloverleaf secondary structure, while the other 21 do not. The AT-rich region, with a 90.5% A + T content and a length of 389 bp, is located between the cox3 and tRNA-Ala genes. Comparison with the mitogenomes of 29 species of Spiruromorpha infraorder, belonging to different families, demonstrates that M. capsulata mitogenome shared the common characteristics of most of them. The phylogeny constructions yielded phylogenies that were in agreement with the obtained previously by using sequences and gene order data of mitogenomes.

Single-worm long-read sequencing reveals genome diversity in free-living nematodes

Obtaining sufficient genetic material from a limited biological source is currently the primary operational bottleneck in studies investigating biodiversity and genome evolution. In this study, we employed multiple displacement amplification (MDA) and Smartseq2 to amplify nanograms of genomic DNA and mRNA, respectively, from individual Caenorhabditis elegans. Although reduced genome coverage was observed in repetitive regions, we produced assemblies covering 98% of the reference genome using long-read sequences generated with Oxford Nanopore Technologies (ONT). Annotation with the sequenced transcriptome coupled with the available assembly revealed that gene predictions were more accurate, complete and contained far fewer false positives than de novo transcriptome assembly approaches. We sampled and sequenced the genomes and transcriptomes of 13 nematodes from early-branching species in Chromadoria, Dorylaimia and Enoplia. The basal Chromadoria and Enoplia species had larger genome sizes, ranging from 136.6 to 738.8 Mb, compared with those in the other clades. Nine mitogenomes were fully assembled, and displayed a complete lack of synteny to other species. Phylogenomic analyses based on the new annotations revealed strong support for Enoplia as sister to the rest of Nematoda. Our result demonstrates the robustness of MDA in combination with ONT, paving the way for the study of genome diversity in the phylum Nematoda and beyond.

De novo assembled mitogenome analysis of Trichuris trichiura from Korean individuals using nanopore-based long-read sequencing technology

Knowledge about mitogenomes has been proven to be essential in human parasite diagnostics and understanding of their diversity. However, the lack of substantial data for comparative analysis is still a challenge in Trichuris trichiura research. To provide high quality mitogenomes, we utilized long-read sequencing technology of Oxford Nanopore Technologies (ONT) to better resolve repetitive regions and to construct de novo mitogenome assembly minimizing reference biases. In this study, we got three de novo assembled mitogenomes of T. trichiura isolated from Korean individuals. These circular complete mitogenomes of T. trichiura are 14,508 bp, 14,441 bp, and 14,440 bp in length. A total of 37 predicted genes were identified consisting of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNAs) genes, two ribosomal RNA (rRNA) genes (rrnS and rrnL), and two non-coding regions. Interestingly, the assembled mitogenome has up to six times longer AT-rich regions than previous reference sequences, thus proving the advantage of long-read sequencing in resolving unreported non-coding regions. Furthermore, variant detection and phylogenetic analysis using concatenated protein coding genes, cox1, rrnL, and nd1 genes confirmed the distinct molecular identity of this newly assembled mitogenome while at the same time showing high genetic relationship with sequences from China or Tanzania. Our study provided a new set of reference mitogenome with better contiguity and resolved repetitive regions that could be used for meaningful phylogenetic analysis to further understand disease transmission and parasite biology.

Evolutionary Rates, Divergence Rates, and Performance of Individual Mitochondrial Genes Based on Phylogenetic Analysis of Copepoda

Copepoda is a large and diverse group of crustaceans, which is widely distributed worldwide. It encompasses roughly 9 orders, whose phylogeny remains unresolved. We sequenced the complete mitochondrial genome (mitogenome) of Sinergasilus major (Markevich, 1940) and used it to explore the phylogeny and mitogenomic evolution of Copepoda. The mitogenome of S. major (14,588 bp) encodes the standard 37 genes as well as a putative control region, and molecular features are highly conserved compared to other Copepoda mitogenomes. Comparative analyses indicated that the nad2 gene has relatively high nucleotide diversity and evolutionary rate, as well as the largest amount of phylogenetic information. These results indicate that nad2 may be a better marker to investigate phylogenetic relationships among closely related species in Copepoda than the commonly used cox1 gene. The sister-group relationship of Siphonostomatoida and Cyclopoida was recovered with strong support in our study. The only topological ambiguity was found within Cyclopoida, which might be caused by the rapid evolution and sparse taxon sampling of this lineage. More taxa and genes should be used to reconstruct the Copepoda phylogeny in the future.

Context-dependent responses of terrestrial invertebrates to anthropogenic nitrogen enrichment: A meta-analysis

Anthropogenic increases in nitrogen (N) concentrations in the environment are affecting plant diversity and ecosystems worldwide, but relatively little is known about N impacts on terrestrial invertebrate communities. Here, we performed an exploratory meta-analysis of 4365 observations from 126 publications reporting on the richness (number of taxa) or abundance (number of individuals per taxon) of terrestrial arthropods or nematodes in relation to N addition. We found that the response of invertebrates to N enrichment is highly dependent on both species' traits and local climate. The abundance of arthropods with incomplete metamorphosis, including agricultural pest species, increased in response to N enrichment. In contrast, arthropods exhibiting complete or no metamorphosis, including pollinators and detritivores, showed a declining abundance trend with increasing N enrichment, particularly in warmer climates. These contrasting and context-dependent responses may explain why we detected no overall response of arthropod richness. For nematodes, the abundance response to N enrichment was dependent on mean annual precipitation and varied between feeding guilds. We found a declining trend in abundance with N enrichment in dry areas and an increasing trend in wet areas, with slopes differing between feeding guilds. For example, at mean levels of precipitation, bacterivore abundance showed a positive trend in response to N addition while fungivore abundance declined. We further observed an overall decline in nematode richness with N addition. These N-induced changes in invertebrate communities could have negative consequences for various ecosystem functions and services, including those contributing to human food production.

Comparative Mitogenomics of Jumping Spiders with First Complete Mitochondrial Genomes of Euophryini (Araneae: Salticidae)

Salticidae is the most species-rich family of spiders with diverse morphology, ecology and behavior. However, the characteristics of the mitogenomes within this group are poorly understood with relatively few well-characterized complete mitochondrial genomes. In this study, we provide completely annotated mitogenomes for Corythalia opima and Parabathippus shelfordi, which represent the first complete mitogenomes of the tribe Euophryini of Salticidae. The features and characteristics of the mitochondrial genomes are elucidated for Salticidae by thoroughly comparing the known well-characterized mitogenomes. The gene rearrangement between trnL2 and trnN was found in two jumping spider species, Corythalia opima and Heliophanus lineiventris Simon, 1868. Additionally, the rearrangement of nad1 to between trnE and trnF found in Asemonea sichuanensis Song & Chai, 1992 is the first protein-coding gene rearrangement in Salticidae, which may have an important phylogenetic implication for the family. Tandem repeats of various copy numbers and lengths were discovered in three jumping spider species. The codon usage analyses showed that the evolution of codon usage bias in salticid mitogenomes was affected by both selection and mutational pressure, but selection may have played a more important role. The phylogenetic analyses provided insight into the taxonomy of Colopsus longipalpis (Żabka, 1985). The data presented in this study will improve our understanding of the evolution of mitochondrial genomes within Salticidae.

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.

The Complete Mitogenome of Toxocara vitulorum: Novel In-Sights into the Phylogenetics in Toxocaridae

Toxocara vitulorum (Ascaridida: Nematoda) is one of the most common intestinal nematodes of cattle and buffalos and, therefore, represents a serious threat to their populations worldwide. Despite its significance in veterinary health the epidemiology, population genetics, and molecular ecology of this nematode remain poorly understood. The mitogenome can yield a foundation for studying these areas and assist in the surveillance and control of T. vitulorum. Herein, the first whole mitogenome of T. vitulorum was sequenced utilizing Illumina technology and characterized with bioinformatic pipeline analyses. The entire genome of T. vitulorum was 15,045 bp in length and contained 12 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), and two ribosomal RNAs (rRNAs). The gene arrangement (GA) of T. vitulorum was similar to those of other Toxocara species under GA3. The whole genome showed significant levels of AT and GC skew. Comparative mitogenomics including sequence identities, Ka/Ks, and sliding window analysis, indicated a purifying selection of 12 PCGs with cox1 and nad6 having the lowest and highest evolutionary rate, respectively. Whole amino acid sequence-based phylogenetic analysis supported a novel sister-species relationship of T. vitulorum with the congeneric species Toxocara canis, Toxocara cati, and Toxocara malaysiensis in the family Toxocaridae. Further, 12 (PCGs) single gene-based phylogenies suggested that nad4 and nad6 genes shared same topological trees with that of the whole genome, suggesting that these genes were suitable as novel genetic markers for phylogenetic and evolutionary studies of Ascaridida species. This complete mitogenome of T. vitulorum refined phylogenetic relationships in Toxocaridae and provided the resource of markers for population genetics, systematics, and epidemiology of this bovine nematode.

The Complete Mitochondrial Genome of Entemnotrochus rumphii, a Living Fossil for Vetigastropoda (Mollusca: Gastropoda)

Pleurotomarioidea represents a truly isolated and basally diverging lineage in Vetigastropoda (Mollusca: Gastropoda) whose fossil record can date back to the late Cambrian, thus providing rare insights into the evolutionary history of molluscs. Here, we sequenced and assembled the complete mitochondrial genome of one representative species from Pleurotomarioidea-Entemnotrochus rumphii (Schepman, 1879)-of which the mitogenome is 15,795 bp in length, including 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes. The nucleotide composition was biased toward AT, and A + T content reached 65.2%. E. rumphii was recovered as sister to all other living vetigastropods according to mitogenome-based phylogenetic analysis. The mitochondrial gene order was consistent with major vetigastropods and the hypothetical ancestral gastropoda, suggesting the deep conservation of mitogenome arrangement in Vetigastropoda.

Comprehensive Molecular Characterization of the Mitochondrial Genome of the Takin Lungworm Varestrongylus eleguneniensis (Strongylida: Protostrongylidae)

The takin lungworm Varestrongylus eleguneniensis (Strongylida: Protostrongylidae) causes lethal bronchopneumonia and represents severe threats to captive and wild populations. However, until now there has been very limited information available concerning the molecular epidemiology and evolutionary biology of V. eleguneniensis. Mitochondrial genomes (mtDNAs) can provide resources for investigations in these areas and, therefore, can assist with the surveillance and control of this lungworm. Herein, the complete mtDNA of V. eleguneniensis was sequenced and characterized with Illumina pipeline analyses. This circular genome (13,625 bp) encoded twelve protein-coding genes (PCGs), two rRNAs, and twenty-two tRNAs, with notable levels of AT and GC skews. Comparative genomics revealed a purifying selection among PCGs, with cox1 and nad6 having the lowest and the highest evolutionary rate, respectively. Genome-wide phylogenies showed a close relationship between V. eleguneniensis and Protostrongylus rufescens in Strongylida. Single gene (PCGs or rRNAs)-based phylogenies indicated that cox1 and nad5 genes shared the same family-level topology with that inferred from genomic datasets, suggesting that both genes could be suitable genetic markers for evolutionary and phylogenetic studies of Strongylida species. This was the first mtDNA of any member of the genus Varestrongylus, and its comprehensive molecular characterization represents a new resource for systematic, population genetic and evolutionary biological studies of Varestrongylus lungworms in wildlife.

pSATdb: a database of mitochondrial common, polymorphic, and unique microsatellites

The polymorphic microSATellites database (pSATdb) provides information on common, polymorphic, and unique mitochondrial microsatellites.

Microsatellites, also termed as simple sequence repeats, are repetitive tracts in a DNA sequence, typically consisting of one to six nucleotides. These repeats are found in all genomes and play key roles in phylogeny and species identification. Microsatellites are highly polymorphic, and their length may differ from species to species. There are several online resources dedicated to mitochondria; however, comprehensive information is not available about the length variation of mitochondrial microsatellites. Therefore, to explore it between species among a genus, we have developed a database named pSATdb (polymorphic microSATellites database; https://lms.snu.edu.in/pSATdb/). pSATdb contains 28,710 perfect microsatellites identified across 5,976 mitochondrial genome (mt-genome) sequences from 1,576 genera which includes 1,535 (5,846 mt-genome) and 41 (130 mt-genome) genera of Metazoa and Viridiplantae, respectively. pSATdb is the only database which provides genus-wise information about the length variation of mitochondrial microsatellites. Because of the emerging role of microsatellites in genomics studies, the identified common, polymorphic, and unique microsatellites stored in pSATdb will be effectively useful in various studies including genetic diversity, mapping, marker-assisted selection, and comparative population studies.

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.

Phylogenomic Analysis of the Phylum Nematoda: Conflicts and Congruences With Morphology, 18S rRNA, and Mitogenomes

Phylogenetic relationships within many lineages of the phylum Nematoda remain unresolved, despite numerous morphology-based and molecular analyses. We performed several phylogenomic analyses using 286 published genomes and transcriptomes and 19 new transcriptomes by focusing on Trichinellida, Spirurina, Rhabditina, and Tylenchina separately, and by analyzing a selection of species from the whole phylum Nematoda. The phylogeny of Trichinellida supported the division of Trichinella into encapsulated and non-encapsulated species and placed them as sister to Trichuris. The Spirurina subtree supported the clades formed by species from Ascaridomorpha and Spiruromorpha respectively, but did not support Dracunculoidea. The analysis of Tylenchina supported a clade that included all sampled species from Tylenchomorpha and placed it as sister to clades that included sampled species from Cephalobomorpha and Panagrolaimomorpha, supporting the hypothesis that postulates the single origin of the stomatostylet. The Rhabditina subtree placed a clade composed of all sampled species from Diplogastridae as sister to a lineage consisting of paraphyletic Rhabditidae, a single representative of Heterorhabditidae and a clade composed of sampled species belonging to Strongylida. It also strongly supported all suborders within Strongylida. In the phylum-wide analysis, a clade composed of all sampled species belonging to Enoplia were consistently placed as sister to Dorylaimia + Chromadoria. The topology of the Nematoda backbone was consistent with previous studies, including polyphyletic placement of sampled representatives of Monhysterida and Araeolaimida.

The re-description of Synoecnemahirsutum Timm, 1959 (Synoecneminae, Ungellidae, Drilonematoidea) from a pheretimoid earthworm in Vietnam with the analysis of its phylogenetic relationships

Synoecnemahirsutum Timm, 1959 (Ungellidae, Drilonematoidea), found in the body cavity of the pheretimoid earthworm at the border of Laos and Vietnam, was re-described and illustrated. The mitochondrial genome of S.hirsutum obtained with Illumina HiSeq sequencing is the first annotated mitochondrial genome as a representative of the superfamily Drilonematoidea. The phylogeny inferred from the analysis of 12 mitochondrial genes has shown some similarity of S.hirsutum with a cephalobid Acrobeloidesvarius.

What lies behind the curtain: Cryptic diversity in helminth parasites of human and veterinary importance

Parasite cryptic species are morphologically indistinguishable but genetically distinct organisms, leading to taxa with unclear species boundaries. Speciation mechanisms such as cospeciation, host colonization, taxon pulse, and oscillation may lead to the emergence of cryptic species, influencing host-parasite interactions, parasite ecology, distribution, and biodiversity. The study of cryptic species diversity in helminth parasites of human and veterinary importance has gained relevance, since their distribution may affect clinical and epidemiological features such as pathogenicity, virulence, drug resistance and susceptibility, mortality, and morbidity, ultimately affecting patient management, course, and outcome of treatment. At the same time, the need for recognition of cryptic species diversity has implied a transition from morphological to molecular diagnostic methods, which are becoming more available and accessible in parasitology. Here, we discuss the general approaches for cryptic species delineation and summarize some examples found in nematodes, trematodes and cestodes of medical and veterinary importance, along with the clinical implications of their taxonomic status. Lastly, we highlight the need for the correct interpretation of molecular information, and the correct use of definitions when reporting or describing new cryptic species in parasitology, since molecular and morphological data should be integrated whenever possible.

•

Cryptic diversity has been described in helminths of human and animal importance.

•

Cryptic species are morphologically indistinguishable but genetically distinct organisms.

•

These entities emerge by different evolutionary and speciation mechanisms.

•

Analysis of molecular and morphological data is needed for cryptic species delimitation.

•

Cryptic diversity may affect pathogenicity, virulence and drug resistance of helminths.

Molecular typing of Strongyloides stercoralis in Latin America, the clinical connection

This study analysed Strongyloides stercoralis genetic variability based on a 404 bp region of the cox1 gene from Latin-American samples in a clinical context including epidemiological, diagnosis and follow-up variables. A prospective, descriptive, observational study was conducted to evaluate clinical and parasitological evolution after ivermectin treatment of 41 patients infected with S. stercoralis. Reactivation of the disease was defined both by clinical symptoms appearance and/or direct larvae detection 30 days after treatment or later. We described 10 haplotypes organized in two clusters. Most frequent variants were also described in the Asian continent in human (HP24 and HP93) and canine (HP24) samples. Clinical presentation (intestinal, severe, cutaneous and asymptomatic), immunological status and eosinophil count were not associated with specific haplotypes or clusters. Nevertheless, presence of cluster 1 haplotypes during diagnosis increased the risk of reactivation with an odds ratio (OR) of 7.51 [confidence interval (CI) 95% 1.38–44.29, P = 0.026]. In contrast, reactivation probability was 83 times lower if cluster 2 (I152V mutation) was detected (OR = 0.17, CI 95% 0.02–0.80, P = 0.02). This is the first analysis of S. stercoralis cox1 diversity in the clinical context. Determination of clusters during the diagnosis could facilitate and improve the design of follow-up strategies to prevent severe reactivations of this chronic disease.

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.

Deep in the systematics of Camallanidae (Nematoda): using integrative taxonomy to better understand the phylogeny and consistency of diagnostic traits

Due to conflicts between classic and molecular systematics of Camallanidae, different data types were used for the first time, to better understand the evolutionary history and taxa consistency within this family. Genetic [18S and 28S rDNA; cytochrome c oxidase subunit I (COI) mtDNA], morphological and life history traits were used to infer phylogenies using Bayesian inference, reconstructed from separated and concatenated datasets. The consistency of tree and morphological traits was evaluated using the consistency index. Characters were mapped on the trees and the phylogenetic informativeness of genetic markers was estimated. Phylogenetic informativeness of 18S provided better resolution for outer nodes, COI for inners and 28S had an intermediate profile. New sequences for two camallanid species were obtained. Phylogenies of genetic and concatenated data largely agreed, showing more divergence in the COI dataset, due to its higher mutation rate vs stable morphology for diagnosing higher taxa. No genus sustained monophyly. The lack of autapomorphy and phylogenetic proximity supported the partition of Batrachocamallanus as synonym of Procamallanus and Spirocamallanus, which should not be considered as subgenera. Although traits of buccal capsule, male tail, habitat, host and biogeographic were highly consistent, intrinsic patterns varied according to different taxa assemblages. Morphological systematics of Camallanidae, based on buccal capsule, is artificial for certain taxa.

Mitogenomics and Evolutionary History of Rodent Whipworms (Trichuris spp.) Originating from Three Biogeographic Regions

Trichuris spp. is a widespread nematode which parasitizes a wide range of mammalian hosts including rodents, the most diverse mammalian order. However, genetic data on rodent whipworms are still scarce, with only one published whole genome (Trichuris muris) despite an increasing demand for whole genome data. We sequenced the whipworm mitogenomes from seven rodent hosts belonging to three biogeographic regions (Palearctic, Afrotropical, and Indomalayan), including three previously described species: Trichuris cossoni, Trichurisarvicolae, and Trichurismastomysi. We assembled and annotated two complete and five almost complete mitogenomes (lacking only the long non-coding region) and performed comparative genomic and phylogenetic analyses. All the mitogenomes are circular, have the same organisation, and consist of 13 protein-coding, 2 rRNA, and 22 tRNA genes. The phylogenetic analysis supports geographical clustering of whipworm species and indicates that T. mastomysi found in Eastern Africa is able to infect multiple closely related rodent hosts. Our results are informative for species delimitation based on mitochondrial markers and could be further used in studies on phylogeny, phylogeography, and population genetics of rodent whipworms

Molluscan mitochondrial genomes break the rules

The first animal mitochondrial genomes to be sequenced were of several vertebrates and model organisms, and the consistency of genomic features found has led to a 'textbook description'. However, a more broad phylogenetic sampling of complete animal mitochondrial genomes has found many cases where these features do not exist, and the phylum Mollusca is especially replete with these exceptions. The characterization of full mollusc mitogenomes required considerable effort involving challenging molecular biology, but has created an enormous catalogue of surprising deviations from that textbook description, including wide variation in size, radical genome rearrangements, gene duplications and losses, the introduction of novel genes, and a complex system of inheritance dubbed 'doubly uniparental inheritance'. Here, we review the extraordinary variation in architecture, molecular functioning and intergenerational transmission of molluscan mitochondrial genomes. Such features represent a great potential for the discovery of biological history, processes and functions that are novel for animal mitochondrial genomes. This provides a model system for studying the evolution and the manifold roles that mitochondria play in organismal physiology, and many ways that the study of mitochondrial genomes are useful for phylogeny and population biology. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

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.

Phylogenetic Analysis of Mitogenomic Data Sets Resolves the Relationship of Seven Macropostrongyloides Species from Australian Macropodid and Vombatid Marsupials

Nematodes of the genus Macropostrongyloides inhabit the large intestines or stomachs of macropodid (kangaroos and wallabies) and vombatid (wombats) marsupials. This study established the relationships of seven species of Macropostrongyloides using mitochondrial (mt) protein amino acid sequence data sets. Phylogenetic analyses revealed that species of Macropostrongyloides (M. lasiorhini, M. baylisi, M. yamagutii, M. spearei, M. mawsonae and M. woodi) from the large intestines of their hosts formed a monophyletic assemblage with strong nodal support to the exclusion of M. dissimilis from the stomach of the swamp wallaby. Furthermore, the mitochondrial protein-coding genes provided greater insights into the diversity and phylogeny of the genus Macropostrongyloides; such data sets could potentially be used to elucidate the relationships among other parasitic nematodes of Australian marsupials.