Improved phylogenomic sampling of free-living nematodes enhances resolution of higher-level nematode phylogeny

Pristionchus - Beetle associations: Towards a new natural history

The free-living nematode Pristionchus pacificus has been established as a model system in integrative evolutionary biology by combining laboratory studies with field work and evolutionary biology. Multiple genetic, molecular and experimental tools and a collection of more than 2,500 P. pacificus strains and more than 50 Pristionchus species, which are available as living cultures or frozen stock collections, support research on various life history traits. Species of Pristionchus exhibit a number of complex traits unknown from Caenorhabditis elegans and most other free-living nematodes. First, P. pacificus can form two alternative mouth forms, an example of developmental plasticity that is increasingly studied to investigate the role of plasticity as a facilitator of evolutionary novelty. More than a decade of work has identified associated genetic and epigenetic mechanisms and revealed the evolutionary and ecological significance of feeding structure plasticity. Second, one of the two mouth morphs results in predatory behavior against other nematodes and is currently used to investigate the neurobiology of predation. Third, potential predation results in the risk of cannibalism among conspecifics. Strikingly, Pristionchus nematodes have developed a self-recognition system that allows the distinction of self (kin) and non-self. Given all these organismal features, this nematode has recently been considered a key example for research towards a new natural history (West-Eberhard, 2024). Here, I summarize recent work on Pristionchus with a focus on a 'new natural history'. In addition, I review some recent studies that indicate an interaction of Pristionchus with EPNs that was suggested based on various surveys in different ecological habitats.

SID-2 is a conserved extracellular vesicle protein that is not associated with environmental RNAi in parasitic nematodes

In the free-living nematode Caenorhabditis elegans, the transmembrane protein SID-2 imports double-stranded RNA into intestinal cells to trigger systemic RNA interference (RNAi), allowing organisms to sense and respond to environmental cues such as the presence of pathogens. This process, known as environmental RNAi, has not been observed in the most closely related parasites that are also within clade V. Previous sequence-based searches failed to identify sid-2 orthologues in available clade V parasite genomes. In this study, we identified sid-2 orthologues in these parasites using genome synteny and protein structure-based comparison, following identification of a SID-2 orthologue in extracellular vesicles from the murine intestinal parasitic nematode Heligmosomoides bakeri. Expression of GFP-tagged H. bakeri SID-2 in C. elegans showed similar localization to the intestinal apical membrane as seen for GFP-tagged C. elegans SID-2, and further showed mobility in intestinal cells in vesicle-like structures. We tested the capacity of H. bakeri SID-2 to functionally complement environmental RNAi in a C. elegans SID-2 null mutant and show that H. bakeri SID-2 does not rescue the phenotype in this context. Our work identifies SID-2 as a highly abundant EV protein whose ancestral function may be unrelated to environmental RNAi, and rather highlights an association with extracellular vesicles in nematodes.

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

In this study, we report a biological 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 complete mitochondrial genome sequence of this organism and show through dN/dS analysis evidence of positive selection in H. mephisto cytochrome c oxidase subunits. Seventeen of these positively selected amino acid substitutions were located in proximity to the H- and K-pathway proton channels of the complex. Surprisingly, the H. mephisto cytochrome c oxidase completely shuts down at low temperatures (20 °C), leading to a 4.8-fold reduction in the transmembrane proton gradient (ΔΨm) compared to optimal temperature (37 °C). Direct measurement of oxygen consumption found a corresponding 4.6-fold drop at 20 °C compared to 37 °C. Correspondingly, the lifecycle of H. mephisto takes four times longer at low temperature than at higher. This elegant evolutionary adaptation creates a finely-tuned mitochondrial temperature sensor, allowing this ectothermic organism to maximize its reproductive success across varying environmental temperatures.

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.

A set of nematode rRNA cistron databases and a primer assessment tool to enable more flexible and comprehensive metabarcoding

The ITS-2-rRNA has been particularly useful for nematode metabarcoding but does not resolve all phylogenetic relationships, and reference sequences are not available for many nematode species. This is a particular issue when metabarcoding complex communities such as wildlife parasites or terrestrial and aquatic free-living nematode communities. We have used markerDB to produce four databases of distinct regions of the rRNA cistron: the 18S rRNA gene, the 28S rRNA gene, the ITS-1 intergenic spacer and the region spanning ITS-1_5.8S_ITS-2. These databases comprise 2645, 254, 13,461 and 10,107 unique full-length sequences representing 1391, 204, 1837 and 1322 nematode species, respectively. The comparative analysis illustrates the complementary value but also reveals a better representation of Clade III, IV and V than Clade I and Clade II nematodes in each case. Although the ITS-1 database includes the largest number of unique full-length sequences, the 18S rRNA database provides the widest taxonomic coverage. We also developed PrimerTC, a tool to assess primer sequence conservation across any reference sequence database, and have applied it to evaluate a large number of previously published rRNA cistron primers. We identified sets of primers that currently provide the broadest taxonomic coverage for each rRNA marker across the nematode phylum. These new resources will facilitate more comprehensive metabarcoding of nematode communities using either short-read or long-read sequencing platforms. Further, PrimerTC is available as a simple WebApp to guide or assess PCR primer design for any genetic marker and/or taxonomic group beyond the nematode phylum.

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.

Genome assembly and annotation of the mermithid nematode Mermis nigrescens

Genetic studies of nematodes have been dominated by Caenorhabditis elegans as a model species. A lack of genomic resources has limited the expansion of genetic research to other groups of nematodes. Here, we report a draft genome assembly of a mermithid nematode, Mermis nigrescens. Mermithidae are insect parasitic nematodes with hosts including a wide range of terrestrial arthropods. We sequenced, assembled, and annotated the whole genome of M. nigrescens using nanopore long reads and 10X Chromium link reads. The assembly is 524 Mb in size consisting of 867 scaffolds. The N50 value is 2.42 Mb, and half of the assembly is in the 30 longest scaffolds. The assembly BUSCO score from the eukaryotic database (eukaryota_odb10) indicates that the genome is 86.7% complete and 5.1% partial. The genome has a high level of heterozygosity (6.6%) with a repeat content of 83.98%. mRNA-seq reads from different sized nematodes (≤2 cm, 3.5-7 cm, and >7 cm body length) representing different developmental stages were also generated and used for the genome annotation. Using ab initio and evidence-based gene model predictions, 12,313 protein-coding genes and 24,186 mRNAs were annotated. These genomic resources will help researchers investigate the various aspects of the biology and host-parasite interactions of mermithid nematodes.

Newly identified nematodes from the Great Salt Lake are associated with microbialites and specially adapted to hypersaline conditions

Extreme environments enable the study of simplified food-webs and serve as models for evolutionary bottlenecks and early Earth ecology. We investigated the biodiversity of invertebrate meiofauna in the benthic zone of the Great Salt Lake (GSL), Utah, USA, one of the most hypersaline lake systems in the world. The hypersaline bays within the GSL are currently thought to support only two multicellular animals: brine fly larvae and brine shrimp. Here, we report the presence, habitat, and microbial interactions of novel free-living nematodes. Nematode diversity drops dramatically along a salinity gradient from a freshwater river into the south arm of the lake. In Gilbert Bay, nematodes primarily inhabit reef-like organosedimentary structures built by bacteria called microbialites. These structures likely provide a protective barrier to UV and aridity, and bacterial associations within them may support life in hypersaline environments. Notably, sampling from Owens Lake, another terminal lake in the Great Basin that lacks microbialites, did not recover nematodes from similar salinities. Phylogenetic divergence suggests that GSL nematodes represent previously undescribed members of the family Monhysteridae-one of the dominant fauna of the abyssal zone and deep-sea hydrothermal vents. These findings update our understanding of halophile ecosystems and the habitable limit of animals.

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.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and the sex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

Acquired stress resilience through bacteria-to-nematode interdomain horizontal gene transfer

Natural selection drives the acquisition of organismal resilience traits to protect against adverse environments. Horizontal gene transfer (HGT) is an important evolutionary mechanism for the acquisition of novel traits, including metazoan acquisitions in immunity, metabolic, and reproduction function via interdomain HGT (iHGT) from bacteria. Here, we report that the nematode gene rml-3 has been acquired by iHGT from bacteria and that it enables exoskeleton resilience and protection against environmental toxins in Caenorhabditis elegans. Phylogenetic analysis reveals that diverse nematode RML-3 proteins form a single monophyletic clade most similar to bacterial enzymes that biosynthesize L-rhamnose, a cell-wall polysaccharide component. C. elegans rml-3 is highly expressed during larval development and upregulated in developing seam cells upon heat stress and during the stress-resistant dauer stage. rml-3 deficiency impairs cuticle integrity, barrier functions, and nematode stress resilience, phenotypes that can be rescued by exogenous L-rhamnose. We propose that interdomain HGT of an ancient bacterial rml-3 homolog has enabled L-rhamnose biosynthesis in nematodes, facilitating cuticle integrity and organismal resilience to environmental stressors during evolution. These findings highlight a remarkable contribution of iHGT on metazoan evolution conferred by the domestication of a bacterial gene.

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.

Telomeric repeat evolution in the phylum Nematoda revealed by high-quality genome assemblies and subtelomere structures

Telomeres are composed of tandem arrays of telomeric-repeat motifs (TRMs) and telomere-binding proteins (TBPs), which are responsible for ensuring end-protection and end-replication of chromosomes. TRMs are highly conserved owing to the sequence specificity of TBPs, although significant alterations in TRM have been observed in several taxa, except Nematoda. We used public whole-genome sequencing data sets to analyze putative TRMs of 100 nematode species and determined that three distinct branches included specific novel TRMs, suggesting that evolutionary alterations in TRMs occurred in Nematoda. We focused on one of the three branches, the Panagrolaimidae family, and performed a de novo assembly of four high-quality draft genomes of the canonical (TTAGGC) and novel TRM (TTAGAC) isolates; the latter genomes revealed densely clustered arrays of the novel TRM. We then comprehensively analyzed the subtelomeric regions of the genomes to infer how the novel TRM evolved. We identified DNA damage-repair signatures in subtelomeric sequences that were representative of consequences of telomere maintenance mechanisms by alternative lengthening of telomeres. We propose a hypothetical scenario in which TTAGAC-containing units are clustered in subtelomeric regions and pre-existing TBPs capable of binding both canonical and novel TRMs aided the evolution of the novel TRM in the Panagrolaimidae family.

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.

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.

Acquired stress resilience through bacteria-to-nematode horizontal gene transfer

Natural selection drives acquisition of organismal resilience traits to protect against adverse environments. Horizontal gene transfer (HGT) is an important evolutionary mechanism for the acquisition of novel traits, including metazoan acquisition of functions in immunity, metabolism, and reproduction via interdomain HGT (iHGT) from bacteria. We report that the nematode gene rml-3, which was acquired by iHGT from bacteria, enables exoskeleton resilience and protection against environmental toxins in C. elegans. Phylogenetic analysis reveals that diverse nematode RML-3 proteins form a single monophyletic clade most highly similar to bacterial enzymes that biosynthesize L-rhamnose to build cell wall polysaccharides. C. elegans rml-3 is regulated in developing seam cells by heat stress and stress-resistant dauer stage. Importantly, rml-3 deficiency impairs cuticle integrity, barrier functions and organismal stress resilience, phenotypes that are rescued by exogenous L-rhamnose. We propose that iHGT of an ancient bacterial rml-3 homolog enables L-rhamnose biosynthesis in nematodes that facilitates cuticle integrity and organismal resilience in adaptation to environmental stresses during evolution. These findings highlight the remarkable contribution of iHGT on metazoan evolution that is conferred by the domestication of bacterial genes.

Systematic and parasite-host relationship by Baruscapillaria appendiculata in Phalacrocorax brasilianus collected from Marajó Island, State of Pará, Brazil

The genus Baruscapillaria Moravec, 1982 has six valid species recorded in birds Phalacrocoracidae, namely Baruscapillaria appendiculata Freitas, 1933, B. spiculata Freitas, 1933, B. carbonis (Dubinin & Dubinina, 1940), B. jaenschi (Johnston & Mawson, 1945), B. phalacrocoraxi (Borgarenko, 1975) and B. rudolphii Moravec, Scholz and Našincová, 1994. Helminthological tests carried out on cormorants of the species Phalacrocorax brasilianus (Gmelin), a migratory bird that occurs in the northeast of the State of Pará, Brazil, demonstrate B. appendiculata parasitizing the cloaca of these birds, through light microscopy, scanning electron microscopy and molecular biology. These studies allowed a redescription of males and females of this nematode in these hosts and in this geographical area through integrative taxonomy. The occurrence of lesions in the cloaca caused by this nematode parasite was registered using histological analysis. This is a new geographic report for this nematode.

Transcriptomics of Cruznema velatum (Nematoda: Rhabditidae) with a redescription of the species

Cruznema velatum isolated from soil in a chestnut orchard located at Guangdong province, China, is redescribed with morphology, molecular barcoding sequences, and transcriptome data. The morphological comparison for C. velatum and six other valid species is provided. Phylogeny analysis suggests genus Cruznema is monophyletic. The species is amphimix, can be cultured with Escherichia coli in 7-9 days from egg to egg-laying adult, and has a lifespan of 11 to 14 days at 20°C. The transcription data generated 45,366 unigenes; 29.9%, 31.3%, 24.8%, and 18.6% of unigenes were annotated in KOG, SwissProt, GO, and KEGG, respectively. Further gene function analysis demonstrated that C. velatum share the same riboflavin, lipoic acid, and vitamin B6 metabolic pathways with Caenorhabditis elegans and Pristionchus pacificus.

A novel nematode species from the Siberian permafrost shares adaptive mechanisms for cryptobiotic survival with C. elegans dauer larva

Some organisms in nature have developed the ability to enter a state of suspended metabolism called cryptobiosis when environmental conditions are unfavorable. This state-transition requires execution of a combination of genetic and biochemical pathways that enable the organism to survive for prolonged periods. Recently, nematode individuals have been reanimated from Siberian permafrost after remaining in cryptobiosis. Preliminary analysis indicates that these nematodes belong to the genera Panagrolaimus and Plectus. Here, we present precise radiocarbon dating indicating that the Panagrolaimus individuals have remained in cryptobiosis since the late Pleistocene (~46,000 years). Phylogenetic inference based on our genome assembly and a detailed morphological analysis demonstrate that they belong to an undescribed species, which we named Panagrolaimus kolymaensis. Comparative genome analysis revealed that the molecular toolkit for cryptobiosis in P. kolymaensis and in C. elegans is partly orthologous. We show that biochemical mechanisms employed by these two species to survive desiccation and freezing under laboratory conditions are similar. Our experimental evidence also reveals that C. elegans dauer larvae can remain viable for longer periods in suspended animation than previously reported. Altogether, our findings demonstrate that nematodes evolved mechanisms potentially allowing them to suspend life over geological time scales.

Description, molecular identification and pathological lesions of Huffmanela persica sp. nov. (Nematoda: Trichosomoididae: Huffmanelinae) from the daggertooth pike conger Muraenesox cinereus

BACKGROUND

The genus Huffmanela Moravec, 1987 (Nematoda, Trichosomoididae, Huffmanelinae), represents a group of nematodes that infect both marine and freshwater fish, and the main gross feature of infection with different species of the genus is the presence of noticeable dark spots or tracks within the parasitized tissues. The purpose of this study was to describe morphologically and morphometrically the eggs of a new marine species of Huffmanela (Huffmanela persica sp. nov.), which was found in the form of black spots in the ovary and the tunica serosa of the stomach of the daggertooth pike conger (Muraenesox cinereus). The new species differs from Huffmanela hamo, another species reported from musculature of this host in Japan, in egg metrics, eggshell features and targeted organ. Molecular identification and pathological examination of the lesions caused by the new species are also reported.

METHODS

Nematode eggs with varying degrees of development were separated from the infected tissues (ovary and tunica serosa of stomach) and investigated using light and scanning electron microscopy. Different species-specific markers (small subunit ribosomal DNA, 18S; large subunit ribosomal DNA, 28S; internal transcribed spacer, ITS) were used for molecular identification and phylogenetic study of the new species. Infected tissues were fixed in buffered formalin for pathological investigations.

RESULTS

The fully developed eggs of H. persica sp. nov. are distinguished from those previously described from this host on the basis of their measurements (size, 54-68 × 31-43 µm; polar plugs, 6.4-9.7 × 8.4-12 µm; shell thickness, 3.5-6.1 µm) and a delicate but ornate uterine layer (UL) covering the entire eggshell including the polar plugs. Histopathological examination revealed a fibro-granulomatous inflammation in the ovary and the serosal layer of the stomach of infected fish. Maximum-likelihood (ML) phylogenetic analysis recovered a sister relationship between the new species of marine origin and Huffmanela species previously collected from freshwater hosts.

CONCLUSIONS

The present study is the first to report the molecular characterization and phylogenetic position of a teleost-associated marine species of the genus Huffmanela. A comprehensive list of nominal and innominate populations of Huffmanela is also provided.

Transcriptional patterns of sexual dimorphism and in host developmental programs in the model parasitic nematode Heligmosomoides bakeri

BACKGROUND

Heligmosomoides bakeri (often mistaken for Heligmosomoides polygyrus) is a promising model for parasitic nematodes with the key advantage of being amenable to study and manipulation within a controlled laboratory environment. While draft genome sequences are available for this worm, which allow for comparative genomic analyses between nematodes, there is a notable lack of information on its gene expression.

METHODS

We generated biologically replicated RNA-seq datasets from samples taken throughout the parasitic life of H. bakeri. RNA from tissue-dwelling and lumen-dwelling worms, collected under a dissection microscope, was sequenced on an Illumina platform.

RESULTS

We find extensive transcriptional sexual dimorphism throughout the fourth larval and adult stages of this parasite and identify alternative splicing, glycosylation, and ubiquitination as particularly important processes for establishing and/or maintaining sex-specific gene expression in this species. We find sex-linked differences in transcription related to aging and oxidative and osmotic stress responses. We observe a starvation-like signature among transcripts whose expression is consistently upregulated in males, which may reflect a higher energy expenditure by male worms. We detect evidence of increased importance for anaerobic respiration among the adult worms, which coincides with the parasite's migration into the physiologically hypoxic environment of the intestinal lumen. Furthermore, we hypothesize that oxygen concentration may be an important driver of the worms encysting in the intestinal mucosa as larvae, which not only fully exposes the worms to their host's immune system but also shapes many of the interactions between the host and parasite. We find stage- and sex-specific variation in the expression of immunomodulatory genes and in anthelmintic targets.

CONCLUSIONS

We examine how different the male and female worms are at the molecular level and describe major developmental events that occur in the worm, which extend our understanding of the interactions between this parasite and its host. In addition to generating new hypotheses for follow-up experiments into the worm's behavior, physiology, and metabolism, our datasets enable future more in-depth comparisons between nematodes to better define the utility of H. bakeri as a model for parasitic nematodes in general.

On the origin and evolution of RNA editing in metazoans

Extensive adenosine-to-inosine (A-to-I) editing of nuclear-transcribed mRNAs is the hallmark of metazoan transcriptional regulation. Here, by profiling the RNA editomes of 22 species that cover major groups of Holozoa, we provide substantial evidence supporting A-to-I mRNA editing as a regulatory innovation originating in the last common ancestor of extant metazoans. This ancient biochemistry process is preserved in most extant metazoan phyla and primarily targets endogenous double-stranded RNA (dsRNA) formed by evolutionarily young repeats. We also find intermolecular pairing of sense-antisense transcripts as an important mechanism for forming dsRNA substrates for A-to-I editing in some but not all lineages. Likewise, recoding editing is rarely shared across lineages but preferentially targets genes involved in neural and cytoskeleton systems in bilaterians. We conclude that metazoan A-to-I editing might first emerge as a safeguard mechanism against repeat-derived dsRNA and was later co-opted into diverse biological processes due to its mutagenic nature.

The improved genome of the nematode Parapristionchus giblindavisi provides insights into lineage-specific gene family evolution

Nematodes such as Caenorhabditis elegans and Pristionchus pacificus are extremely successful model organisms for comparative biology. Several studies have shown that phenotypic novelty but also conserved processes are controlled by taxon-restricted genes. To trace back the evolution of such new or rapidly evolving genes, a robust phylogenomic framework is indispensable. Here, we present an improved version of the genome of Parapristionchus giblindavisi which is the only known member of the sister group of Pristionchus. Relative to the previous short-read assembly, the new genome is based on long reads and displays higher levels of contiguity, completeness, and correctness. Specifically, the number of contigs dropped from over 7,303 to 735 resulting in an N50 increase from 112 to 791 kb. We made use of the new genome to revisit the evolution of multiple gene families. This revealed Pristionchus-specific expansions of several environmentally responsive gene families and a Pristionchus-specific loss of the de novo purine biosynthesis pathway. Focusing on the evolution of sulfatases and sulfotransferases, which control the mouth form plasticity in P. pacificus, reveals differences in copy number and genomic configurations between the genera Pristionchus and Parapristionchus. Altogether, this demonstrates the utility of the P. giblindavisi genome to date and polarizes lineage-specific patterns.

Paraphyly of Marimermithida refines primary routes of transition to parasitism in roundworms

Parasitic life-strategies in the phylum Nematoda (roundworms) are remarkably diverse and intricate in terms of evolution and taxonomy. By analysing novel rDNA data obtained on rare host-associated groups with unusual biology, we reveal paraphyly of the last major taxon with uncertain higher-rank classification that united solely parasitic nematodes (Marimermithida) to show that primarily marine parasitism only emerged independently and repeatedly in a few free-living lineages. We report secondary seaward ingression of land-based parasites (Mermithida) via invading hosts in the subtidal zone to illustrate the host-borne scenario of oceanic fish and mammal colonization by primarily terrestrial parasites (Spiruria). We also present the first molecular data on marine nematodes from unicellular hosts (foraminiferan protozoans) to demonstrate the independent origins of exploitative nematode associations at a microscopic scale. We argue that, in contrast with primarily intestinal associations arising from saprotrophy and commensalism, non-intestinal host capture (colonization of host body cavity or internal organs) is likely to be a primary route of transition to truly exploitative parasitism in roundworms. Predispositions to host capture in nematode morphology, ecology and life cycles imply its evolution as part of innate pre-adaptations to crossing environmental boundaries to enable multiple successful transitions to parasitism in the phylum history.

Ecological lipidology

Dietary lipids (DLs), particularly sterols and fatty acids, are precursors for endogenous lipids that, unusually for macronutrients, shape cellular and organismal function long after ingestion. These functions – cell membrane structure, intracellular signalling, and hormonal activity – vary with the identity of DLs, and scale up to influence health, survival, and reproductive fitness, thereby affecting evolutionary change. Our Ecological Lipidology approach integrates biochemical mechanisms and molecular cell biology into evolution and nutritional ecology. It exposes our need to understand environmental impacts on lipidomes, the lipid specificity of cell functions, and predicts the evolution of lipid-based diet choices. Broad interdisciplinary implications of Ecological Lipidology include food web alterations, species responses to environmental change, as well as sex differences and lifestyle impacts on human nutrition, and opportunities for DL-based therapies.

Genetic diversity and identity of Ascaris worms from human and pig hosts in Thailand

Ascaris roundworms are of public health and socio-economic importance worldwide. They are conventionally attributed to two taxa - A. lumbricoides infecting principally human and A. suum infecting principally pig. Phylogenomic analysis has revealed that Ascaris worms from both human and pig are represented in Clades A and B. A recent study indicates that the Ascaris worms from human and pig in Thailand belong to Clade A. We examined adult Ascaris worms from human and pig in Thailand by means of the partial sequences of three mitochondrial genes (cox1, cox2 and nad1) and concatenation of these genes. Phylogenomic analysis indicates that two isolates (H1,H2) of A. lumbricoides from human belonged to Clade B; one isolate (H3) belonged to Clade A (based on cox1, cox2 and concatenated sequences) or as an outlier to Clades A and B (based on nad1 sequences). All the eight isolates of A. suum from pig clustered in Clade A. The partial nad1 and the concatenated sequences revealed two lineages of A. suum isolates which were distinct from the two A. lumbricoides isolates of Clade B. It is evident that greater genetic diversity, and a more robust phylogeny, could be uncovered by the application of multiple genes. In sum, the present study reveals the presence in Thailand of A. lumbricoides from human in Clades A and B which necessitates appropriate treatment and control measures; Clades A and B have been reported to contain haplotypes of Ascaris worms from both human and pig in other parts of the world. A country wide study is needed to elucidate the identity, distribution, prevalence, cross transmission, genetic diversity and phylogeny of the Ascaris worms in Thailand.

Multiple Pristionchus pacificus genomes reveal distinct evolutionary dynamics between de novo candidates and duplicated genes

The birth of new genes is a major molecular innovation driving phenotypic diversity across all domains of life. Although repurposing of existing protein-coding material by duplication is considered the main process of new gene formation, recent studies have discovered thousands of transcriptionally active sequences as a rich source of new genes. However, differential loss rates have to be assumed to reconcile the high birth rates of these incipient de novo genes with the dominance of ancient gene families in individual genomes. Here, we test this rapid turnover hypothesis in the context of the nematode model organism Pristionchus pacificus We extended the existing species-level phylogenomic framework by sequencing the genomes of six divergent P. pacificus strains. We used these data to study the evolutionary dynamics of different age classes and categories of origin at a population level. Contrasting de novo candidates with new families that arose by duplication and divergence from known genes, we find that de novo candidates are typically shorter, show less expression, and are overrepresented on the sex chromosome. Although the contribution of de novo candidates increases toward young age classes, multiple comparisons within the same age class showed significantly higher attrition in de novo candidates than in known genes. Similarly, young genes remain under weak evolutionary constraints with de novo candidates representing the fastest evolving subcategory. Altogether, this study provides empirical evidence for the rapid turnover hypothesis and highlights the importance of the evolutionary timescale when quantifying the contribution of different mechanisms toward new gene formation.

Huffmanela cf. huffmani (Nematoda: Trichosomoididae) infecting swim bladder, peritoneum, and gonad of variable platyfish, Xiphophorus variatus (Cyprinodontiformes: Poeciliidae) and eastern mosquitofish, Gambusia holbrooki (Poeciliidae) in Florida; taxonomy, phylogenetic analysis, and pathological changes

Variable platyfish, Xiphophorus variatus (Meek, 1904) (Cyprinodontiformes: Poeciliidae) and eastern mosquitofish, Gambusia holbrooki Girard, 1859 (Poeciliidae) from earthen ponds in west central Florida were examined for parasitic infections. At necropsy, we observed myriad nematodes (adults and eggs), which we identified as Huffmanela cf. huffmani, infecting the swim bladder, gonad, and visceral peritoneum. Nucleotide sequences (small subunit ribosomal DNA, 18S) of H. cf. huffmani from variable platyfish and eastern mosquitofish were identical; likewise for newly obtained 18S sequences of Huffmanela huffmani Moravec, 1987 from the swim bladder of red breast sunfish, Lepomis auritus (Linnaeus, 1758) (Centrarchiformes: Centrarchidae) and warmouth, Lepomis gulosus (Cuvier, 1829) from the San Marcos River (type locality for Huffmanela huffmani Moravec, 1987), Texas. The sequences of H. huffmani and H. cf. huffmani differed by 7 (1%) nucleotides. Pathological changes comprised proliferation of the tunica externa of the swim bladder in low-intensity infections in addition to inflammation, proliferation, and tissue necrosis of swim bladder, peritoneum, and gonad in high-intensity infections. The lesion was severe, affecting the cellular constituents of the swim bladder wall and reducing the size of the swim bladder lumen; potentially reducing swim bladder physiological efficiency. The present study is the first record of a freshwater species of Huffmanela Moravec, 1987 from beyond the San Marcos River, first record of a species of Huffmanela from a livebearer, first nucleotide sequences and phylogenetic analysis for Huffmanela, and first evidence that an infection by a species of Huffmanela causes pathological changes that could impact organ function.

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.

Evolution of sexual systems, sex chromosomes and sex-linked gene transcription in flatworms and roundworms

Many species with separate male and female individuals (termed ‘gonochorism’ in animals) have sex-linked genome regions. Here, we investigate evolutionary changes when genome regions become completely sex-linked, by analyses of multiple species of flatworms (Platyhelminthes; among which schistosomes recently evolved gonochorism from ancestral hermaphroditism), and roundworms (Nematoda) which have undergone independent translocations of different autosomes. Although neither the evolution of gonochorism nor translocations fusing ancestrally autosomal regions to sex chromosomes causes inevitable loss of recombination, we document that formerly recombining regions show genomic signatures of recombination suppression in both taxa, and become strongly genetically degenerated, with a loss of most genes. Comparisons with hermaphroditic flatworm transcriptomes show masculinisation and some defeminisation in schistosome gonad gene expression. We also find evidence that evolution of sex-linkage in nematodes is accompanied by transcriptional changes and dosage compensation. Our analyses also identify sex-linked genes that could assist future research aimed at controlling some of these important parasites.

Transitions between hermaphroditic and separate sexes are relatively understudied in animals compared to pants. Here, Wang et al. reconstruct the evolution of separate sexes in the flatworms and complex changes of sex chromosomes in the roundworms.

Myrica cerifera, a Medicinal Plant of the Lumbee Tribe, has Antibacterial and Nematicidal Properties

Currently threatening the world of medicine is a growing number of antibiotic-resistant diseases. More specifically, bacteria and nematodes have gained resistance to many of the world's leading antibiotics and nematicides, respectively, making infections more difficult to treat. Subsequently, these parasitic organisms are able to continue damaging crops and other living organisms like humans without strong interference. To help people and the environment, the development of new and novel antibiotics is vital. Previous research suggests that phytochemicals are a potential solution that will not only help inhibit bacterial growth but also reduce nematode survival. We hypothesized that Myrica cerifera, a plant often used by the Lumbee tribe to treat illness, possesses antibacterial and nematicidal properties. To answer our hypothesis, we began by collecting plant specimens to extract material for biological assays and to subsequently isolate and elucidate the structures of active components. The extract was evaluated for antibacterial properties with an agar diffusion assay and then nematicidal properties using Caenorhabditis elegans. M. cerifera extract was added onto an agar lawn at various doses, and the nematodes' lifespans were scored. The findings of this study show that extracts of this plant, more commonly referred to as 'wax myrtle', do significantly decrease the lifespan of C. elegans and increase the zone of inhibition for Staphylococcus epidermidis and Staphylococcus aureus. In addition, two compounds were isolated and characterized through chemical extraction, chromatographic separation, and spectroscopic analysis. These compounds could potentially be used to treat bacterial and nematode infections.

Non-canonical nematode endogenous retroviruses resulting from RNA virus glycoprotein gene capture by a metavirus

Reverse-transcribing retroviruses exist as horizontally transmitted infectious agents or vertically transmitted endogenous retroviruses (ERVs) resident in eukaryotic genomes, and they are phylogenetically related to the long terminal repeat (LTR) class of retrotransposons. ERVs and retrotransposons are often distinguished only by the presence or absence of a gene encoding the envelope glycoprotein (env). Endogenous elements of the virus family Metaviridae include the insect-restricted Errantivirus genus of ERVs, for which some members possess env, and the pan-eukaryotic Metavirus genus that lacks an envelope glycoprotein gene. Here we report a novel Nematoda endogenous retrovirus (NERV) clade with core retroviral genes arranged uniquely as a continuous gag-env-pro-pol ORF. Reverse transcriptase sequences were phylogenetically related to metaviruses, but envelope glycoprotein sequences resembled those of the Nyamiviridae and Chrysoviridae RNA virus families, suggesting env gene capture during host cell infection by an RNA virus. NERVs were monophyletic, restricted to the nematode subclass Chromadoria, and included additional ORFs for a small hypothetical protein or a large Upf1-like RNA-dependent AAA-ATPase/helicase indicative of viral transduction of a host gene. Provirus LTR identity, low copy number, ORF integrity and segregation of three loci in Meloidogyne incognita, taken together with detection of NERV transcriptional activity, support potential infectivity of NERVs, along with their recent emergence and integration. Altogether, NERVs constitute a new and distinct Metaviridae lineage demonstrating retroviral evolution through sequential heterologous gene capture events.

Spatial transcriptomics reveals antiparasitic targets associated with essential behaviors in the human parasite Brugia malayi

Lymphatic filariasis (LF) is a chronic debilitating neglected tropical disease (NTD) caused by mosquito-transmitted nematodes that afflicts over 60 million people. Control of LF relies on routine mass drug administration with antiparasitics that clear circulating larval parasites but are ineffective against adults. The development of effective adulticides is hampered by a poor understanding of the processes and tissues driving parasite survival in the host. The adult filariae head region contains essential tissues that control parasite feeding, sensory, secretory, and reproductive behaviors, which express promising molecular substrates for the development of antifilarial drugs, vaccines, and diagnostics. We have adapted spatial transcriptomic approaches to map gene expression patterns across these prioritized but historically intractable head tissues. Spatial and tissue-resolved data reveal distinct biases in the origins of known drug targets and secreted antigens. These data were used to identify potential new drug and vaccine targets, including putative hidden antigens expressed in the alimentary canal, and to spatially associate receptor subunits belonging to druggable families. Spatial transcriptomic approaches provide a powerful resource to aid gene function inference and seed antiparasitic discovery pipelines across helminths of relevance to human and animal health.

Nematode chromosomes

The nematode Caenorhabditis elegans has shed light on many aspects of eukaryotic biology, including genetics, development, cell biology, and genomics. A major factor in the success of C. elegans as a model organism has been the availability, since the late 1990s, of an essentially gap-free and well-annotated nuclear genome sequence, divided among 6 chromosomes. In this review, we discuss the structure, function, and biology of C. elegans chromosomes and then provide a general perspective on chromosome biology in other diverse nematode species. We highlight malleable chromosome features including centromeres, telomeres, and repetitive elements, as well as the remarkable process of programmed DNA elimination (historically described as chromatin diminution) that induces loss of portions of the genome in somatic cells of a handful of nematode species. An exciting future prospect is that nematode species may enable experimental approaches to study chromosome features and to test models of chromosome evolution. In the long term, fundamental insights regarding how speciation is integrated with chromosome biology may be revealed.

A High-Throughput Phenotypic Screen of the 'Pandemic Response Box' Identifies a Quinoline Derivative with Significant Anthelmintic Activity

Parasitic nematodes cause diseases in livestock animals and major economic losses to the agricultural industry worldwide. Nematodes of the order Strongylida, including Haemonchus contortus, are particularly important. The excessive use of anthelmintic compounds to treat infections and disease has led to widespread resistance to these compounds in nematodes, such that there is a need for new anthelmintics with distinctive mechanisms of action. With a focus on discovering new anthelmintic entities, we screened 400 chemically diverse compounds within the 'Pandemic Response Box' (from Medicines for Malaria Venture, MMV) for activity against H. contortus and its free-living relative, Caenorhabditis elegans-a model organism. Using established phenotypic assays, test compounds were evaluated in vitro for their ability to inhibit the motility and/or development of H. contortus and C. elegans. Dose-response evaluations identified a compound, MMV1581032, that significantly the motility of H. contortus larvae (IC₅₀ = 3.4 ± 1.1 μM) and young adults of C. elegans (IC₅₀ = 7.1 ± 4.6 μM), and the development of H. contortus larvae (IC₅₀ = 2.2 ± 0.7 μM). The favourable characteristics of MMV1581032, such as suitable physicochemical properties and an efficient, cost-effective pathway to analogue synthesis, indicates a promising candidate for further evaluation as a nematocide. Future work will focus on a structure-activity relationship investigation of this chemical scaffold, a toxicity assessment of potent analogues and a mechanism/mode of action investigation.

Small RNA pathways in the nematode Ascaris in the absence of piRNAs

Small RNA pathways play key and diverse regulatory roles in C. elegans, but our understanding of their conservation and contributions in other nematodes is limited. We analyzed small RNA pathways in the divergent parasitic nematode Ascaris. Ascaris has ten Argonautes with five worm-specific Argonautes (WAGOs) that associate with secondary 5’-triphosphate 22-24G-RNAs. These small RNAs target repetitive sequences or mature mRNAs and are similar to the C. elegans mutator, nuclear, and CSR-1 small RNA pathways. Even in the absence of a piRNA pathway, Ascaris CSR-1 may still function to “license” as well as fine-tune or repress gene expression. Ascaris ALG-4 and its associated 26G-RNAs target and likely repress specific mRNAs during testis meiosis. Ascaris WAGO small RNAs demonstrate target plasticity changing their targets between repeats and mRNAs during development. We provide a unique and comprehensive view of mRNA and small RNA expression throughout spermatogenesis. Overall, our study illustrates the conservation, divergence, dynamics, and flexibility of small RNA pathways in nematodes.

The parasitic nematode Ascaris lacks piRNAs. Here the authors compare Argonaute proteins and small RNAs from C. elegans and Ascaris, expanding our understanding of the conservation, divergence, and flexibility of Argonautes and small RNA pathways in nematodes.

Comprehensive species sampling and sophisticated algorithmic approaches refute the monophyly of Arachnida

Deciphering the evolutionary relationships of Chelicerata (arachnids, horseshoe crabs, and allied taxa) has proven notoriously difficult, due to their ancient rapid radiation and the incidence of elevated evolutionary rates in several lineages. Although conflicting hypotheses prevail in morphological and molecular data sets alike, the monophyly of Arachnida is nearly universally accepted, despite historical lack of support in molecular data sets. Some phylotranscriptomic analyses have recovered arachnid monophyly, but these did not sample all living orders, whereas analyses including all orders have failed to recover Arachnida. To understand this conflict, we assembled a data set of 506 high-quality genomes and transcriptomes, sampling all living orders of Chelicerata with high occupancy and rigorous approaches to orthology inference. Our analyses consistently recovered the nested placement of horseshoe crabs within a paraphyletic Arachnida. This result was insensitive to variation in evolutionary rates of genes, complexity of the substitution models, and alternative algorithmic approaches to species tree inference. Investigation of sources of systematic bias showed that genes and sites that recover arachnid monophyly are enriched in noise and exhibit low information content. To test the impact of morphological data, we generated a 514-taxon morphological data matrix of extant and fossil Chelicerata, analyzed in tandem with the molecular matrix. Combined analyses recovered the clade Merostomata (the marine orders Xiphosura, Eurypterida, and Chasmataspidida), but merostomates appeared nested within Arachnida. Our results suggest that morphological convergence resulting from adaptations to life in terrestrial habitats has driven the historical perception of arachnid monophyly, paralleling the history of numerous other invertebrate terrestrial groups.

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.

Phylogenomic analyses of mud dragons (Kinorhyncha)

Mud dragons (Kinorhyncha) are microscopic invertebrates, inhabiting marine sediments across the globe from intertidal to hadal depths. They are segmented, moulting animals like arthropods, but grouping with the unsegmented priapulans and loriciferans within Ecdysozoa. There are more than 300 species of kinorhynchs described within 31 genera and 11 families, however, their evolutionary relationships have so far only been investigated using morphology and a few molecular markers. Here we aim to resolve the relationships and classification of major clades within Kinorhyncha using transcriptomic data. In addition, we wish to revisit the position of three indistinctly segmented, aberrant genera in order to reconstruct the evolution of distinct segmentation within the group. We conducted a phylogenomic analysis of Kinorhyncha including 21 kinorhynch transcriptomes (of which 18 are new) representing 15 genera, and seven outgroups including priapulan, loriciferan, nematode and nematomorph transcriptomes. Results show a congruent and robust tree that supports the division of Kinorhyncha into two major clades: Cyclorhagida and Allomalorhagida. Cyclorhagida is composed of three subclades: Xenosomata, Kentrorhagata comb. nov. (including the aberrant Zelinkaderes) and Echinorhagata. Allomalorhagida is composed of two subclades: Pycnophyidae and Anomoirhaga nom. nov. Anomoirhaga nom. nov. accommodates the aberrant genera Cateria (previously nested within Cyclorhagida) and Franciscideres together with five additional genera. The distant and derived positions of the aberrant Zelinkaderes, Cateria and Franciscideres species suggest that their less distinct trunk segmentation evolved convergently, and that segmentation evolved among kinorhynch stem groups.

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

Molecular phylogenetics brought radical changes to our understanding of nematode evolution, resulting in substantial modifications to nematode classification implemented by De Ley and Blaxter and widely accepted now. Numerous phylogenetic studies were subsequently published that both improved and challenged this classification. Here we present a summary of these changes. We created cladograms that summarise phylogenetic relationships within Nematoda using phylum-wide to superfamily-wide molecular phylogenies published in since 2005, and supplemented with the phylogenetic analyses for Enoplia and Chromadoria with the aim of clarifying the position of several taxa. The results show which parts of the Nematode tree are well resolved and understood, and which parts require more research, either by adding taxa that have not been included yet (increasing taxon coverage), or by changing the phylogenetic approach (improving data quality, using different types of data or different methods of analysis). The currently used classification of the phylum Nematoda in many cases does not reflect the phylogeny and in itself requires numerous improvements and rearrangements.

What Is an “Arachnid”? Consensus, Consilience, and Confirmation Bias in the Phylogenetics of Chelicerata

The basal phylogeny of Chelicerata is one of the opaquest parts of the animal Tree of Life, defying resolution despite application of thousands of loci and millions of sites. At the forefront of the debate over chelicerate relationships is the monophyly of Arachnida, which has been refuted by most analyses of molecular sequence data. A number of phylogenomic datasets have suggested that Xiphosura (horseshoe crabs) are derived arachnids, refuting the traditional understanding of arachnid monophyly. This result is regarded as controversial, not least by paleontologists and morphologists, due to the widespread perception that arachnid monophyly is unambiguously supported by morphological data. Moreover, some molecular datasets have been able to recover arachnid monophyly, galvanizing the belief that any result that challenges arachnid monophyly is artefactual. Here, we explore the problems of distinguishing phylogenetic signal from noise through a series of in silico experiments, focusing on datasets that have recently supported arachnid monophyly. We assess the claim that filtering by saturation rate is a valid criterion for recovering Arachnida. We demonstrate that neither saturation rate, nor the ability to assemble a molecular phylogenetic dataset supporting a given outcome with maximal nodal support, is a guarantor of phylogenetic accuracy. Separately, we review empirical morphological phylogenetic datasets to examine characters supporting Arachnida and the downstream implication of a single colonization of terrestrial habitats. We show that morphological support of arachnid monophyly is contingent upon a small number of ambiguous or incorrectly coded characters, most of these tautologically linked to adaptation to terrestrial habitats.

Assessing the suitability of mitochondrial and nuclear DNA genetic markers for molecular systematics and species identification of helminths

BACKGROUND

Genetic markers are employed widely in molecular studies, and their utility depends on the degree of sequence variation, which dictates the type of application for which they are suited. Consequently, the suitability of a genetic marker for any specific application is complicated by its properties and usage across studies. To provide a yardstick for future users, in this study we assess the suitability of genetic markers for molecular systematics and species identification in helminths and provide an estimate of the cut-off genetic distances per taxonomic level.

METHODS

We assessed four classes of genetic markers, namely nuclear ribosomal internal transcribed spacers, nuclear rRNA, mitochondrial rRNA and mitochondrial protein-coding genes, based on certain properties that are important for species identification and molecular systematics. For molecular identification, these properties are inter-species sequence variation; length of reference sequences; easy alignment of sequences; and easy to design universal primers. For molecular systematics, the properties are: average genetic distance from order/suborder to species level; the number of monophyletic clades at the order/suborder level; length of reference sequences; easy alignment of sequences; easy to design universal primers; and absence of nucleotide substitution saturation. Estimation of the cut-off genetic distances was performed using the 'K-means' clustering algorithm.

RESULTS

The nuclear rRNA genes exhibited the lowest sequence variation, whereas the mitochondrial genes exhibited relatively higher variation across the three groups of helminths. Also, the nuclear and mitochondrial rRNA genes were the best possible genetic markers for helminth molecular systematics, whereas the mitochondrial protein-coding and rRNA genes were suitable for molecular identification. We also revealed that a general gauge of genetic distances might not be adequate, using evidence from the wide range of genetic distances among nematodes.

CONCLUSION

This study assessed the suitability of DNA genetic markers for application in molecular systematics and molecular identification of helminths. We provide a novel way of analyzing genetic distances to generate suitable cut-off values for each taxonomic level using the 'K-means' clustering algorithm. The estimated cut-off genetic distance values, together with the summary of the utility and limitations of each class of genetic markers, are useful information that can benefit researchers conducting molecular studies on helminths.

A draft transcriptome of a parasite Neocamacolaimus parasiticus (Camacolaimidae, Plectida)

Camacolaimidae is a clade of nematodes that include both free-living epistrate feeding forms and parasites of marine protozoans and invertebrates. Neocamacolaimus parasiticus is a parasite of marine polychaete worms. Given its phylogenetic affinities to free-living species, Neocamacolaimus can be a reference for research of the origin of parasitism in an aquatic environment. Here, we present a draft transcriptome obtained from a single post-parasitic juvenile individual of this species. The final assembly consists of 19,180 protein coding sequences (including isoforms) with the following BUSCO scores for Nematoda: 65.38% complete, 9.06% partial, and 25.56% missing, and for Metazoa: 79.45% complete, 3.17% partial, and 17.38% missing.

The community-curated Pristionchus pacificus genome facilitates automated gene annotation improvement in related nematodes

Background

The nematode Pristionchus pacificus is an established model organism for comparative studies with Caenorhabditis elegans. Over the past years, it developed into an independent animal model organism for elucidating the genetic basis of phenotypic plasticity. Community-based curations were employed recently to improve the quality of gene annotations of P. pacificus and to more easily facilitate reverse genetic studies using candidate genes from C. elegans.

Results

Here, I demonstrate that the reannotation of phylogenomic data from nine related nematode species using the community-curated P. pacificus gene set as homology data substantially improves the quality of gene annotations. Benchmarking of universal single copy orthologs (BUSCO) estimates a median completeness of 84% which corresponds to a 9% increase over previous annotations. Nevertheless, the ability to infer gene models based on homology already drops beyond the genus level reflecting the rapid evolution of nematode lineages. This also indicates that the highly curated C. elegans genome is not optimally suited for annotating non-Caenorhabditis genomes based on homology. Furthermore, comparative genomic analysis of apparently missing BUSCO genes indicates a failure of ortholog detection by the BUSCO pipeline due to the insufficient sample size and phylogenetic breadth of the underlying OrthoDB data set. As a consequence, the quality of multiple divergent nematode genomes might be underestimated.

Conclusions

This study highlights the need for optimizing gene annotation protocols and it demonstrates the benefit of a high quality genome for phylogenomic data of related species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07529-x.

SLIDR and SLOPPR: flexible identification of spliced leader trans-splicing and prediction of eukaryotic operons from RNA-Seq data

BACKGROUND

Spliced leader (SL) trans-splicing replaces the 5' end of pre-mRNAs with the spliced leader, an exon derived from a specialised non-coding RNA originating from elsewhere in the genome. This process is essential for resolving polycistronic pre-mRNAs produced by eukaryotic operons into monocistronic transcripts. SL trans-splicing and operons may have independently evolved multiple times throughout Eukarya, yet our understanding of these phenomena is limited to only a few well-characterised organisms, most notably C. elegans and trypanosomes. The primary barrier to systematic discovery and characterisation of SL trans-splicing and operons is the lack of computational tools for exploiting the surge of transcriptomic and genomic resources for a wide range of eukaryotes.

RESULTS

Here we present two novel pipelines that automate the discovery of SLs and the prediction of operons in eukaryotic genomes from RNA-Seq data. SLIDR assembles putative SLs from 5' read tails present after read alignment to a reference genome or transcriptome, which are then verified by interrogating corresponding SL RNA genes for sequence motifs expected in bona fide SL RNA molecules. SLOPPR identifies RNA-Seq reads that contain a given 5' SL sequence, quantifies genome-wide SL trans-splicing events and predicts operons via distinct patterns of SL trans-splicing events across adjacent genes. We tested both pipelines with organisms known to carry out SL trans-splicing and organise their genes into operons, and demonstrate that (1) SLIDR correctly detects expected SLs and often discovers novel SL variants; (2) SLOPPR correctly identifies functionally specialised SLs, correctly predicts known operons and detects plausible novel operons.

CONCLUSIONS

SLIDR and SLOPPR are flexible tools that will accelerate research into the evolutionary dynamics of SL trans-splicing and operons throughout Eukarya and improve gene discovery and annotation for a wide range of eukaryotic genomes. Both pipelines are implemented in Bash and R and are built upon readily available software commonly installed on most bioinformatics servers. Biological insight can be gleaned even from sparse, low-coverage datasets, implying that an untapped wealth of information can be retrieved from existing RNA-Seq datasets as well as from novel full-isoform sequencing protocols as they become more widely available.

Rhabditophanes diutinus a parthenogenetic clade IV nematode with dauer larvae

Comparative studies using non-parasitic model species such as Caenorhabditis elegans, have been very helpful in investigating the basic biology and evolution of parasitic nematodes. However, as phylogenetic distance increases, these comparisons become more difficult, particularly when outside of the nematode clade to which C. elegans belongs (V). One of the reasons C. elegans has nevertheless been used for these comparisons, is that closely related well characterized free-living species that can serve as models for parasites of interest are frequently not available. The Clade IV parasitic nematodes Strongyloides are of great research interest due to their life cycle and other unique biological features, as well as their medical and veterinary importance. Rhabditophanes, a closely related free-living genus, forms part of the Strongyloidoidea nematode superfamily. Rhabditophanes diutinus (= R. sp. KR3021) was included in the recent comparative genomic analysis of the Strongyloididae, providing some insight into the genomic nature of parasitism. However, very little is known about this species, limiting its usefulness as a research model. Here we provide a species description, name the species as R. diutinus and investigate its life cycle and subsequently gene expression in multiple life stages. We identified two previously unreported starvation induced life stages: dauer larvae and arrested J2 (J2A) larvae. The dauer larvae are morphologically similar to and are the same developmental stage as dauers in C. elegans and infective larvae in Strongyloides. As in C. elegans and Strongyloides, dauer formation is inhibited by treatment with dafachronic acid, indicating some genetic control mechanisms are conserved. Similarly, the expression patterns of putative dauer/infective larva control genes resemble each other, in particular between R. diutinus and Strongyloides spp. These findings illustrate and increase the usefulness of R. diutinus as a non-parasitic, easy to work with model species for the Strongyloididae for studying the evolution of parasitism as well as many aspects of the biology of Strongyloides spp, in particular the formation of infective larvae.

Parasitic worms are an issue of great medical, veterinary, agricultural and economic importance, yet little is known about how worms become parasites. Comparative studies with non-parasitic model species like C. elegans have been useful, however, this usefulness decreases as the evolutionary distance between the species increases. One way to combat this is by having more well-studied closely related species to parasites of interest. To address this, we provide information about Rhabditophanes diutinus, a free-living nematode that is part of the same superfamily as the medically and veterinary important Strongyloides parasites. We provide analysis on its life cycle, in particular on two starvation induced life stages, along with gene expression data. Overall, this important information illustrates and improves the use of R. diutinus, as a non-parasitic model species for studying parasite evolution and basic biology within Strongyloides.

Evolution of Developmental GATA Factors in Nematodes

GATA transcription factors are found in animals, plants, and fungi. In animals, they have important developmental roles in controlling specification of cell identities and executing tissue-specific differentiation. The Phylum Nematoda is a diverse group of vermiform animals that inhabit ecological niches all over the world. Both free-living and parasitic species are known, including those that cause human infectious disease. To date, GATA factors in nematodes have been studied almost exclusively in the model system C. elegans and its close relatives. In this study, we use newly available sequences to identify GATA factors across the nematode phylum. We find that most species have fewer than six GATA factors, but some species have 10 or more. Comparisons of gene and protein structure suggest that there were at most two GATA factors at the base of the phylum, which expanded by duplication and modification to result in a core set of four factors. The high degree of structural similarity with the corresponding orthologues in C. elegans suggests that the nematode GATA factors share similar functions in development.