Twenty Years after De Ley and Blaxter-How Far Did We Progress in Understanding the Phylogeny of the Phylum Nematoda?

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.

Evolution of neuropeptide Y/RFamide-like receptors in nematodes

The Neuropeptide Y/RFamide-like receptors belong to the Rhodopsin-like G protein-coupled receptors G protein-coupled receptors (GPCRs) and are involved in functions such as locomotion, feeding and reproduction. With 41 described receptors they form the best-studied group of neuropeptide GPCRs in Caenorhabditis elegans. In order to understand the expansion of the Neuropeptide Y/RFamide-like receptor family in nematodes, we started from the sequences of selected receptor paralogs in C. elegans as query and surveyed the corresponding orthologous sequences in another 159 representative nematode target genomes. To this end we employed a automated pipeline based on ExonMatchSolver, a tool that solves the paralog-to-contig assignment problem. Utilizing subclass-specific HMMs we were able to detect a total of 1557 Neuropeptide Y/RFamide-like receptor sequences (1100 NPRs, 375 FRPRs and 82 C09F12.3) in the 159 target nematode genomes investigated here. These sequences demonstrate a good conservation of the Neuropeptide Y/RFamide-like receptors across the Nematoda and highlight the diversification of the family in nematode evolution. No other genus shares all Neuropeptide Y/RFamide-like receptors with the genus Caenorhabditis. At the same time, we observe large numbers of clade specific duplications and losses of family members across the phylum Nematoda.

Systematics and phylogeny of the entomopathogenic nematobacterial complexes Steinernema-Xenorhabdus and Heterorhabditis-Photorhabdus

Entomopathogenic nematodes of the genera Steinernema and Heterorhabditis, along with their bacterial symbionts from the genera Xenorhabdus and Photorhabdus, respectively, are important biological control agents against agricultural pests. Rapid progress in the development of genomic tools has catalyzed a transformation of the systematics of these organisms, reshaping our understanding of their phylogenetic and cophlylogenetic relationships. In this review, we discuss the major historical events in the taxonomy and systematics of this group of organisms, highlighting the latest advancements in these fields. Additionally, we synthesize information on nematode-bacteria associations and assess the existing evidence regarding their cophylogenetic relationships.

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.

The dilemma of underestimating freshwater biodiversity: morphological and molecular approaches

BACKGROUND

Anthropogenic impacts on freshwater habitats are causing a recent biodiversity decline far greater than that documented for most terrestrial ecosystems. However, knowledge and description of freshwater biodiversity is still limited, especially targeting all size classes to uncover the distribution of biodiversity between different trophic levels. We assessed the biodiversity of the Lower Rhine and associated water bodies in the river's flood plain including the river's main channel, oxbows and gravel-pit lakes, spanning from the level of protists up to the level of larger invertebrate predators and herbivores organized in size classes (nano-, micro, meio- and macrofauna). Morphological diversity was determined by morphotypes, while the molecular diversity (amplicon sequencing variants, ASVs) was assessed through eDNA samples with metabarcoding targeting the V9 region of the 18S rDNA.

RESULTS

Considering all four investigated size classes, the percentage of shared taxa between both approaches eDNA (ASVs with 80-100% sequence similarity to reference sequences) and morphology (morphotypes), was always below 15% (5.4 ± 3.9%). Even with a more stringent filtering of ASVs (98-100% similarity), the overlap of taxa could only reach up to 43% (18.3 ± 12%). We observed low taxonomic resolution of reference sequences from freshwater organisms in public databases for all size classes, especially for nano-, micro-, and meiofauna, furthermore lacking metainformation if species occur in freshwater, marine or terrestrial ecosystems.

CONCLUSIONS

In our study, we provide a combination of morphotype detection and metabarcoding that particularly reveals the diversity in the smaller size classes and furthermore highlights the lack of genetic resources in reference databases for this diversity. Especially for protists (nano- and microfauna), a combination of molecular and morphological approaches is needed to gain the highest possible community resolution. The assessment of freshwater biodiversity needs to account for its sub-structuring in different ecological size classes and across compartments in order to reveal the ecological dimension of diversity and its distribution.

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.

The complete mitochondrial genome of Wellcomia compar (Spirurina: Oxyuridae) and its genome characterization and phylogenetic analysis

Wellcomia compar (Spirurina: Oxyuridae) is a pinworm that infects wild and captive porcupines. Despite clear records of its morphological structure, its genetics, systematics, and biology are poorly understood. This study aimed to determine the complete mitochondrial (mt) genome of W. compar and reconstruct its phylogenetic relationship with other nematodes. We sequenced the complete mt genome of W. comparand conducted phylogenetic analyses using concatenated coding sequences of 12 protein-coding genes (PCGs) by maximum likelihood and Bayesian inference. The complete mt genome is 14,373 bp in size and comprises 36 genes, including 12 protein-coding, two rRNA and 22 tRNA genes. Apart from 28 intergenic regions, one non-coding region and one overlapping region also occur. A comparison of the gene arrangements of Oxyuridomorpha revealed relatively similar features in W. compar and Wellcomia siamensis. Phylogenetic analysis also showed that W. compar and W. siamensis formed a sister group. In Oxyuridomorpha the genetic distance between W. compar and W. siamensis was 0.0805. This study reports, for the first time, the complete W. compar mt genome sequence obtained from Chinese porcupines. It provides genetic markers for investigating the taxonomy, population genetics, and phylogenetics of pinworms from different hosts and has implications for the diagnosis, prevention, and control of parasitic diseases in porcupines and other animals.

Microbiomes of microscopic marine invertebrates do not reveal signatures of phylosymbiosis

Animals and microorganisms often establish close ecological relationships. However, much of our knowledge about animal microbiomes comes from two deeply studied groups: vertebrates and arthropods. To understand interactions on a broader scale of diversity, we characterized the bacterial microbiomes of close to 1,000 microscopic marine invertebrates from 21 phyla, spanning most of the remaining tree of metazoans. Samples were collected from five temperate and tropical locations covering three marine habitats (sediment, water column and intertidal macroalgae) and bacterial microbiomes were characterized using 16S ribosomal RNA gene sequencing. Our data show that, despite their size, these animals harbour bacterial communities that differ from those in the surrounding environment. Distantly related but coexisting invertebrates tend to share many of the same bacteria, suggesting that guilds of microorganisms preferentially associated with animals, but not tied to any specific host lineage, are the main drivers of the ecological relationship. Host identity is a minor factor shaping these microbiomes, which do not show the same correlation with host phylogeny, or 'phylosymbiosis', observed in many large animals. Hence, the current debate on the varying strength of phylosymbiosis within selected lineages should be reframed to account for the possibility that such a pattern might be the exception rather than the rule.

Six First Reports of Pin Nematodes from Portugal, with an Update of the Systematics, Genetic Diversity, and Phylogeny of the Genus Paratylenchus (Nematoda: Tylenchulidae)

Pin nematodes (Paratylenchus spp.) currently comprise 132 species of polyphagous plant ectoparasites with at least seven species that are plant pathogenic emphasizing the need for correct identification to establish an appropriate management strategy. Sequences of highly conserved regions of ribosomal and mitochondrial RNA genes are a powerful species-level diagnostic tool within Tylenchulidae. A nematological survey was conducted from 2019 to 2021 in commercial vineyards distributed in four major wine-producing regions in the Central and South Portugal. Nine populations of Paratylenchus spp. were isolated from the rhizosphere of grapevines that were characterized from morphological data and molecular phylogenetic analysis using two rRNA genes (D2–D3 expansion segments of the 28S, and partial 18S) and a region partial of the COI mtRNA gene. Contrasting morphological hypotheses with molecular data provided rapid detection of six species, specifically P. goodeyi, P. hamatus, P. pedrami, P. tenicaudatus, P. variabilis, and P. veruculatus. Bayesian inference (BI) phylogenetic trees for these molecular markers established phylogenetic relationships underscore the importance of using genomic and molecular phylogenetic data for accurate pin nematode identification at the species level. To our knowledge, this is the first report of Paratylenchus spp. in Portugal, extending the geographical distribution of these species in the Mediterranean Basin, and the first record of P. goodeyi, P. pedrami, P.variabilis, P. veruculatus and P. tenicaudatus parasitizing grapevine.

Comparative reconstruction of the predatory feeding structures of the polyphenic nematode Pristionchus pacificus

Pristionchus pacificus is a nematode model for the developmental genetics of morphological polyphenism, especially at the level of individual cells. Morphological polyphenism in this species includes an evolutionary novelty, moveable teeth, which have enabled predatory feeding in this species and others in its family (Diplogastridae). From transmission electron micrographs of serial thin sections through an adult hermaphrodite of P. pacificus, we three‐dimensionally reconstructed all epithelial and myoepithelial cells and syncytia, corresponding to 74 nuclei, of its face, mouth, and pharynx. We found that the epithelia that produce the predatory morphology of P. pacificus are identical to Caenorhabditis elegans in the number of cell classes and nuclei. However, differences in cell form, spatial relationships, and nucleus position correlate with gross morphological differences from C. elegans and outgroups. Moreover, we identified fine‐structural features, especially in the anteriormost pharyngeal muscles, that underlie the conspicuous, left‐right asymmetry that characterizes the P. pacificus feeding apparatus. Our reconstruction provides an anatomical map for studying the genetics of polyphenism, feeding behavior, and the development of novel form in a satellite model to C. elegans.

All cells making the dimorphic, novel form of an animal with cell constancy were identified. Although the number of cells is fully conserved, divergence in form and connectivity—including fixed asymmetries—sheds light on the origins of this trait.

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.