WormBase ParaSite - a comprehensive resource for helminth genomics

Transcriptome of Taenia solium during in vitro cyst activation and initial growth into the tapeworm stage

The cestode Taenia solium develops as a tapeworm solely in the human intestine, starting from a larva (cyst). Upon maturing, it produces hundreds of thousands of infectious eggs. When ingested by pigs or humans, the eggs develop as cysts that lodge in various tissues, including the brain, leading to neurocysticercosis. Despite advances in understanding cestode biology through genomic and transcriptomic studies, particularly in model organisms, much remains unknown about the activation of T. solium cysts in the human digestive tract and the events that drive the development into adult worms-the stage responsible for dispersing the parasite. We present a transcriptome generated by Next Generation Sequencing from T. solium cysts activated in culture and collected at three different in vitro growth phases, defined by their morphology. Differentially expressed genes and biological processes relevant to activation and growth can be explored with the dataset. The information is valuable for identifying genes that regulate the molecular, metabolic, and cellular events leading to parasite maturation or elements driving its transmission.

Vaccination against helminth IL-33 modulators permits immune-mediated parasite ejection

The murine intestinal nematode Heligmosomoides polygyrus bakeri (Hpb) modulates the host immune response via the Hpb alarmin release inhibitor (HpARI) family (HpARI1/2/3), which acts on interleukin (IL)-33, and the Hpb binds alarmin receptor and inhibits (HpBARI) family (HpBARI and HpBARI_Hom2), which acts on the IL-33 receptor ST2. Here, we find that this immunomodulation is evident only in the first week of infection and affects local and distal tissues. Vaccination with HpARI or HpBARI proteins raises antibody responses that block their immunomodulatory activities: HpARI2 vaccination results in significantly increased type 2 innate lymphoid cells (ILC2s), T helper (Th)2, and serum IL-4 and IL-5 responses, while HpBARI + HpBARI_Hom2 vaccination reverses infection-mediated ST2 suppression and increases Th2 immunity. A cocktail of HpARI2 + HpBARI + HpBARI_Hom2 gives robust protection against infection, associated with stunting of adult parasites, reduced egg burden, increased type 2 immune responses, and intestinal goblet cell expansion. Therefore, vaccination with immunomodulatory proteins can protect the host against infection and can be used as a tool for blocking the effects of specific parasite-derived proteins.

Gastrointestinal nematode parasites of grazing ruminants: a comprehensive literature review of diagnostic methods for quantifying parasitism, larval differentiation and measuring anthelmintic resistance

This review summarises up-to-date research on the diagnosis of gastrointestinal nematode (GIN) infection in livestock and anthelmintic resistance in GIN. It was commissioned to assist funding bodies to prioritise and guide research and extension efforts to improve the health, welfare and productivity of grazing ruminants in the face of challenge with GIN. A comprehensive review of published articles from journals, books and websites was undertaken, with a focus on peer-reviewed articles published between 2000 and 2024 involving genera of GIN in grazing sheep and cattle with economic importance to New Zealand. Suggestions for articles to include were received from 14 experts in GIN diagnostics. This review is a summary of a longer report submitted to the sponsoring organisation. Clinical signs of GIN infection in grazing ruminants in temperate grazing systems are inadequate as triggers for management interventions including anthelmintic treatment as they are visible only after economically significant pathological changes have occurred. Livestock producers benefit from monitoring GIN burdens using faecal egg counts (FEC) or associated signals such as weight gain. In future, they may use remote monitoring devices for activity in animals, as well as estimating pasture larval contamination. Methods of diagnosing GIN infections using automated FEC devices have improved the convenience of monitoring parasite burdens compared with traditional laboratory methods. However, a lack of quality control measures and a gap in training of skilled technicians for larval differentiation may lead to a shortage of diagnostic capability. Current methods of diagnosing anthelmintic resistance, particularly FEC reduction tests, are not likely to be replaced by laboratory assays in the near future and attention should be focused on facilitating application of new FEC technologies for both animal monitoring and resistance diagnosis. Extension and application of currently available methods and technology will improve animal health and productivity in ruminant grazing systems in the short term. Adoption of novel technologies for remote animal monitoring, practical tools for estimating pasture larval contamination and promoting genetic selection for immunity and resilience to GIN in both sheep and cattle will further enhance productivity in the long term.

Identification and characterization of a novel ApeC-containing transmembrane protein family in parasitic flatworms

The Apextrin C-terminal (ApeC) domain is widely distributed across various animal clades. Several ApeC domains have been identified as carbohydrate-binding domains, contributing to extracellular adhesion and mucosal immunity. However, the functions and evolutionary history of most ApeC-containing protein (ACP) families remain largely unexplored. In this study, we identified 213 ACPs in flatworms (Platyhelminthes), with each species containing one to two such proteins. All flatworm ACPs belong to a unique transmembrane protein family, characterized by a length of 1700-2200 amino acids and a distinctive domain architecture (SP-[Kringle]1-2-[ApeC]7-[EGF]1-4-TM) unlike any found in other phyla. This ACP is conserved across all major parasitic flatworm lineages, including flukes (Trematoda), tapeworms (Cestoda) and monogeneans (Monogenea), despite their diverse morphologies and habitats. Notably, it is also present in one group of free-living flatworms, the planarians (Tricladida), suggesting that this ACP originated in the free-living ancestor of parasitic flatworms. This ACP contains seven consecutive ApeC domains, an unparalleled number among known animal proteins. Five of these ApeC domains are highly divergent, necessitating the definition of a new domain model (ApeC_Pla; Pfam: PF24148) for accurate classification. Structural predictions indicate that these ApeC domains adopt a conserved three β-sheet structure. Furthermore, transcriptomic analysis revealed that flatworm ACPs, along with several important serpins and proteases, are predominantly expressed in parenchymal cells and feeding organs, suggesting that ACPs serve as novel marker genes for parenchymal tissue and may be involved in cell adhesion, oral immunity, and parasite-host interactions. Taken together, our findings indicate that this flatworm ACP represents a promising target for vaccine development and provides key insights into the physiology of the parenchyma, a unique flatworm tissue that functions as a substitute for a body cavity.

Transposable Elements Drive Regulatory and Functional Innovation of F-box Genes

Protein domains of transposable elements (TEs) and viruses increase the protein diversity of host genomes by recombining with other protein domains. By screening 10 million eukaryotic proteins, we identified several domains that define multicopy gene families and frequently co-occur with TE/viral domains. Among these, a Tc1/Mariner transposase helix-turn-helix (HTH) domain was captured by F-box genes in the Caenorhabditis genus, creating a new class of F-box genes. For specific members of this class, like fbxa-215, we found that the HTH domain is required for diverse processes including germ granule localization, fertility, and thermotolerance. Furthermore, we provide evidence that Heat Shock Factor 1 (HSF-1) mediates the transcriptional integration of fbxa-215 into the heat shock response by binding to Helitron TEs directly upstream of the fbxa-215 locus. The interactome of HTH-bearing F-box factors suggests roles in post-translational regulation and proteostasis, consistent with established functions of F-box proteins. Based on AlphaFold2 multimer proteome-wide screens, we propose that the HTH domain may diversify the repertoire of protein substrates that F-box factors regulate post-translationally. We also describe an independent capture of a TE domain by F-box genes in zebrafish. In conclusion, we identify two independent TE domain captures by F-box genes in eukaryotes and provide insights into how these novel proteins are integrated within host gene regulatory networks.

Robust Binding Capability and Occasional Gene Loss of Telomere-Binding Proteins Underlying Telomere Evolution in Nematoda

Telomeres, the nucleoprotein complexes that protect the ends of linear chromosomes, are essential for maintaining the stability of eukaryotic genomes. As telomeres generally consist of repetitive DNA associated with specifically bound proteins, telomeric repeat motifs are thought to be difficult to evolve. However, a recent study identified nematodes with telomeric repeats distinct from the canonical TTAGGC motif. Here, we investigated how telomere repeats could have evolved despite the challenge posed by the specificity of telomere-binding proteins (TBPs) to the telomeric DNA in Nematoda. We performed a phylogenetic analysis and electrophoresis mobility shift assays to assess the binding affinities of two TBPs, which displayed different conservation patterns. Our results revealed that the well-conserved protein CEH-37 exhibits limited specificity, unable to distinguish telomeric repeats found in nematodes except for the TTAGGG motif, while the less conserved POT proteins displayed rigid specificity. These findings suggest that the emergence of novel telomeric repeat motifs correlated with the characteristics and evolutionary outcomes of TBPs in Nematoda. Our study not only revealed the dynamics of telomere evolution but also enhanced the understanding of the evolutionary relationship between proteins and DNAs.

Diet-regulated transcriptional plasticity of plant parasites in plant-mutualist environments

Crop pathogens often lack exclusive access to their host and must interact with plants concurrently engaged with numerous other symbionts. Here, we demonstrate that the colonization of hosts by plant-mutualistic mycorrhizal fungi can indirectly induce transcriptional responses of a major plant parasite, the nematode Globodera pallida, via a modified host resource profile. A shift in the resource profile of the root, where the parasite feeds, is perceived and responded to by the parasite through transcriptional changes, potentially to optimize resource intake. Specifically, G. pallida react to reduced host-photosynthate influx due to concurrent mycorrhizal-host symbiosis by upregulating the expression of a sugar transporter (SWEET3) in the nematode intestine. We identify this gene's role in parasite growth and development, regulated by the putative diet-responsive transcription factor Gp-HBL1. Overall, our data unveil a mechanism by which a parasitic animal responds to fluctuations in host plant quality that is induced by a plant-mutualistic fungus, to enhance parasitism and reproduction.

Neuron-specific repression of alternative splicing by the conserved CELF protein UNC-75 in Caenorhabditis elegans

Tissue-regulated alternative exons are dictated by the interplay between cis-elements and trans-regulatory factors such as RNA-binding proteins (RBPs). Despite extensive research on splicing regulation, the full repertoire of these cis and trans features and their evolutionary dynamics across species are yet to be fully characterized. Members of the CUG-binding protein and ETR-like family (CELF) of RBPs are known to play a key role in the regulation of tissue-biased splicing patterns, and when mutated, these proteins have been implicated in a number of neurological and muscular disorders. In this study, we sought to characterize specific mechanisms that drive tissue-specific splicing in vivo of a model switch-like exon regulated by the neuronal-enriched CELF ortholog in Caenorhabditis elegans, UNC-75. Using sequence alignments, we identified deeply conserved intronic UNC-75 binding motifs overlapping the 5' splice site and upstream of the 3' splice site, flanking a strongly neural-repressed alternative exon in the Zonula Occludens gene zoo-1. We confirmed that loss of UNC-75 or mutations in either of these cis-elements lead to substantial de-repression of the alternative exon in neurons. Moreover, mis-expression of UNC-75 in muscle cells is sufficient to induce the neuron-like robust skipping of this alternative exon. Lastly, we demonstrate that overlapping an UNC-75 motif within a heterologous 5' splice site leads to increased skipping of the adjacent alternative exon in an unrelated splicing event. Together, we have demonstrated that a specific configuration and combination of cis elements bound by this important family of RBPs can achieve robust splicing outcomes in vivo.

Bioinformatics Analysis of the Glutamate-Gated Chloride Channel Family in Bursaphelenchus xylophilus

Glutamate-gated chloride channels (GluCls), a class of ion channels found in the nerve and muscle cells of invertebrates, are involved in vital life processes. Bursaphelenchus xylophilus, the pathogen of pine wilt disease, has induced major economic and ecological losses in invaded areas of Asia and Europe. We identified 33 GluCls family members by sequence alignment analysis. A subsequent bioinformatic analysis revealed the physicochemical properties, protein structure, and gene expression patterns in different developmental stages. The results showed that GluCls genes are distributed across all six chromosomes of B. xylophilus. These proteins indicated a relatively conserved structure by NCBI-conserved domains and InterPro analysis. A gene structure analysis revealed that GluCls genes consist of 5 to 14 exons. Expression pattern analysis revealed BxGluCls were extensively involved in the development of second instar larvae of B. xylophilus. Furthermore, BxGluCls15, BxGluCls25, and BxGluCls28 were mainly associated with the development of eggs of B. xylophilus. BxGluCls12, BxGluCls18, and BxGluCls32 were predominantly linked to nematode resistance and adaptation. Investigation the structure and expression patterns of BxGluCls is crucial to understand the developmental trends of B. xylophilus. It also helps identify molecular targets for the development of biopesticides or drugs designed to control this nematode.

Transformer-based deep learning enables improved B-cell epitope prediction in parasitic pathogens: A proof-of-concept study on Fasciola hepatica

BACKGROUND

The identification of B-cell epitopes (BCEs) is fundamental to advancing epitope-based vaccine design, therapeutic antibody development, and diagnostics, such as in neglected tropical diseases caused by parasitic pathogens. However, the structural complexity of parasite antigens and the high cost of experimental validation present certain challenges. Advances in Artificial Intelligence (AI)-driven protein engineering, particularly through machine learning and deep learning, offer efficient solutions to enhance prediction accuracy and reduce experimental costs.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we present deepBCE-Parasite, a Transformer-based deep learning model designed to predict linear BCEs from peptide sequences. By leveraging a state-of-the-art self-attention mechanism, the model achieved remarkable predictive performance, achieving an accuracy of approximately 81% and an AUC of 0.90 in both 10-fold cross-validation and independent testing. Comparative analyses against 12 handcrafted features and four conventional machine learning algorithms (GNB, SVM, RF, and LGBM) highlighted the superior predictive power of the model. As a case study, deepBCE-Parasite predicted eight BCEs from the leucine aminopeptidase (LAP) protein in Fasciola hepatica proteomic data. Dot-blot immunoassays confirmed the specific binding of seven synthetic peptides to positive sera, validating their IgG reactivity and demonstrating the model's efficacy in BCE prediction.

CONCLUSIONS/SIGNIFICANCE

deepBCE-Parasite demonstrates excellent performance in predicting BCEs across diverse parasitic pathogens, offering a valuable tool for advancing the design of epitope-based vaccines, antibodies, and diagnostic applications in parasitology.

Global analysis of ligand-gated ion channel conservation across Platyhelminthes

Ligand-gated ion channels (LGICs) are critical for neurotransmission, mediating responses to neurotransmitters and hormones, and influencing diverse physiological processes. This study identifies and classifies LGICs across Platyhelminthes, with a particular focus on parasitic neodermatans, which impact human and animal health. Using bioinformatics tools, we analyzed LGICs from 41 neodermatan species and expanded our investigation to encompass vertebrates, other invertebrates, and non-bilaterians to trace LGIC evolutionary pathways across Metazoa. We identified 2,269 putative LGICs within neodermatan species, which we classified into the cys-loop, ASIC/Deg/ENaC, iGluR, and P2X families. Our phylogenetic and clustering analyses reveal lineage-specific patterns with distinct evolutionary trajectories for each LGIC family in neodermatans compared to free-living platyhelminths and other taxa. Notably, the ASIC/Deg/ENaC family displayed the greatest degree of neodermatan-specific divergence, while cys-loop and P2X families were more conserved across taxa. To provide insight into their potential physiological roles, we analyzed LGIC expression patterns in Schistosoma mansoni, revealing widespread expression across neuronal and muscle cell types. The distribution of acid-sensing ion channels (ASICs) in both neurons and muscles suggests a role in neuromuscular signalling, while the P2X receptor (Smp_333600) exhibited sex-specific expression, potentially indicating distinct functional roles in males and females. Additionally, several cys-loop acetylcholine and GABA receptors showed differential neuronal and muscle expression, highlighting their likely contributions to cholinergic and inhibitory neurotransmission. These findings underscore the relevance of LGICs in parasite physiology, particularly in neuromuscular and sensory processes, and suggest potential targets for antiparasitic interventions.

Transcriptional phenotype of the anti-parasitic benzodiazepine meclonazepam on the blood fluke Schistosoma mansoni

There are limited control measures for the disease schistosomiasis, despite the fact that infection with parasitic blood flukes affects hundreds of millions of people worldwide. The current treatment, praziquantel, has been in use since the 1980's and there is a concern that drug resistance may emerge with continued monotherapy. Given the need for additional antischistosomal drugs, we have re-visited an old lead, meclonazepam. In comparison to praziquantel, there has been relatively little work on its antiparasitic mechanism. Recent findings indicate that praziquantel and meclonazepam act through distinct receptors, making benzodiazepines a promising chemical series for further exploration. Previous work has profiled the transcriptional changes evoked by praziquantel treatment. Here, we examine in detail schistosome phenotypes evoked by in vitro and in vivo meclonazepam treatment. These data confirm that meclonazepam causes extensive tegument damage and directly kills parasites, as measured by pro-apoptotic caspase activation. In vivo meclonazepam exposure results in differential expression of many genes that are divergent in parasitic flatworms, as well as several gene products implicated in blood feeding and regulation of hemostasis in other parasites. Many of these transcripts are also differentially expressed with praziquantel exposure, which may reflect a common schistosome response to the two drugs. However, despite these similarities in drug response, praziquantel-resistant parasites retain susceptibility to meclonazepam's schistocidal effects. These data provide new insight into the mechanism of antischistosomal benzodiazepines, resolving similarities and differences with the current frontline therapy, praziquantel.

Genome assembly variation and its implications for gene discovery in nematodes

Genome assemblers are a critical component of genome science, but the choice of assembly software and protocols can be daunting. Here, we investigate genome assembly variation and its implications for gene discovery across three nematode species-Caenorhabditis bovis, Haemonchus contortus, and Heligmosomoides bakeri-highlighting the critical interplay between assembly choice and downstream genomic analysis. Selecting commonly used genome assemblers, we generated multiple assemblies for each species, analyzing their structure, completeness, and effect on gene family analysis. Our findings demonstrate that assembly variations can significantly affect gene family composition, with notable differences in gene families important in anthelmintic discovery and immunomodulation. Despite broadly similar performance using various assembly metrics, comparisons of assemblies with a single species revealed underlying structural rearrangements and inconsistencies in gene content, which would affect downstream analyses. This emphasizes the need for continuous refinement of genome assemblies and their annotations.

Biotransformation of anthelmintics in nematodes in relation to drug resistance

In all organisms, the biotransformation of xenobiotics to less toxic and more hydrophilic compounds represents an effective defense strategy. In pathogens, the biotransformation of drugs (used for their elimination from the host) may provide undesirable protective effects that could potentially compromise the drug's efficacy. Accordingly, increased drug deactivation via accelerated biotransformation is now considered as one of the mechanisms of drug resistance. The present study summarizes the current knowledge regarding the biotransformation of anthelmintics, specifically drugs used to treat mainly nematodes, a group of parasites that are a significant health concern for humans and animals. The main biotransformation enzymes are introduced and their roles in anthelmintics metabolism in nematodes are discussed with a particular focus on their potential participation in drug resistance. Similarly, the inducibility of biotransformation enzymes with sublethal doses of anthelmintics is presented in view of its potential contribution to drug resistance development. In the conclusion, the main tasks awaiting scientists in this area are outlined.

Heterochromatin protein 1 (HP1) of Schistosoma mansoni: non-canonical chromatin landscape and oviposition effects

BACKGROUND

Heterochromatin protein 1 (HP1) is widespread in several organisms playing a role in control of gene expression by heterochromatin formation and maintenance of silent chromatin. Schistosoma mansoni is a human parasite that is responsible for Schistosomiasis, a tropical neglected disease in the tropical and subtropical areas in the world, where the intermediate host Biomphalaria glabrata is present.

OBJECTIVES

In this study we attempted to investigate if the SmHP1 is enriched in S. mansoni chromatin in cercariae larvae stage, compared with another larvae stage sporocysts and its importance for S. mansoni life cycle progression and parasite oviposition.

METHODS

We used ChIPmentation with commercial antibody ab109028 that passed in-house quality control. We also used RNA interference, mice infection and histology.

FINDINGS

Our data show that S. mansoni HP1 enrichment is non-canonical with a peak at the transcription end sites of protein coding genes. We did not find strong differences in SmHP1 chromatin landscapes between sporocysts and cercariae. Knock- down of SmHP1 in schistosomula and in vivo experiments in mice unexpectedly increased parasite oviposition.

MAIN CONCLUSIONS

Our results suggest that SmHP1 may influence chromatin structure in a non-canonical way in S. mansoni stages and may play a role in regulation of parasite oviposition.

Protocol to identify transcription factor target genes using TargetOrtho2

TargetOrtho2 uses transcription factor binding site information to predict transcription factor targets in C. elegans, based on an in silico phylogenetic footprinting approach. Here, we present a protocol to identify transcription factor target genes using a new version of TargetOrtho2. We provide instructions for installing TargetOrtho2 and its required suite of programs, for predicting transcription factor target genes, and for updating and adding new genomes to TargetOrtho2.

The SUbventral-Gland Regulator (SUGR-1) of nematode virulence

Pathogens must precisely tailor their gene expression to cause infection. However, a signaling cascade from host signal to effector production has remained elusive for metazoan pathogens. Here, we show that plants contain molecular signals, termed effectostimulins, that activate the first identified regulator of plant-parasitic nematode effectors. SUGR-1 directly binds effector promoters, and is central to a transcriptional network that activates 58 effector genes. Importantly, we demonstrate that downregulation of sugr-1 inhibits parasitism, underlining SUGR-1 signaling as a valuable target for crop protection and food security. This, in the wider context of nematodes as parasites of humans and other animals, has scope for potentially broader impact: Disrupting effector production could, in principle, be applied to any pathogen that secrets effectors.

CelEst: a unified gene regulatory network for estimating transcription factor activities in C. elegans

Transcription factors (TFs) play a pivotal role in orchestrating critical intricate patterns of gene regulation. Although gene expression is complex, differential expression of hundreds of genes is often due to regulation by just a handful of TFs. Despite extensive efforts to elucidate TF-target regulatory relationships in Caenorhabditis elegans, existing experimental datasets cover distinct subsets of TFs and leave data integration challenging. Here, I introduce CelEst, a unified gene regulatory network designed to estimate the activity of 487 distinct C. elegans TFs-∼58% of the total-from gene expression data. To integrate data from ChIP-seq, DNA-binding motifs, and eY1H screens, optimal processing of each data type was benchmarked against a set of TF perturbation RNA-seq experiments. Moreover, I showcase how leveraging TF motif conservation in target promoters across genomes of related species can distinguish highly informative interactions, a strategy which can be applied to many model organisms. Integrated analyses of data from commonly studied conditions including heat shock, bacterial infection, and sex differences validates CelEst's performance and highlights overlooked TFs that likely play major roles in coordinating the transcriptional response to these conditions. CelEst can infer TF activity on a standard laptop computer within minutes. Furthermore, an R Shiny app with a step-by-step guide is provided for the community to perform rapid analysis with minimal coding required. I anticipate that widespread adoption of CelEsT will significantly enhance the interpretive power of transcriptomic experiments, both present and retrospective, thereby advancing our understanding of gene regulation in C. elegans and beyond.

A genome-scale drug discovery pipeline uncovers new therapeutic targets and a unique p97 allosteric binding site in Schistosoma mansoni

Schistosomes are parasitic flatworms that infect more than 200 million people globally. However, there is a shortage of molecular tools that enable the discovery of potential drug targets within schistosomes. Thus, praziquantel has remained the frontline treatment for schistosomiasis despite known liabilities. Here, we have conducted a genome-wide study in S. mansoni using the human druggable genome as a bioinformatic template to identify essential genes within schistosomes bearing similarity to catalogued drug targets. Then, we assessed these candidate targets in silico using a set of unbiased criteria to determine which possess ideal characteristics for a ready-made drug discovery campaign. Following this prioritization, we pursued a parasite p97 ortholog as a bona-fide drug target for the development of therapeutics to treat schistosomiasis. From this effort, we identified a covalent inhibitor series that kills schistosomes through an on-target killing mechanism by disrupting the ubiquitin proteasome system. Fascinatingly, these inhibitors induce a conformational change in the conserved D2 domain P-loop of schistosome p97 upon modification of Cys519. This conformational change reveals an allosteric binding site adjacent to the D2 domain active site reminiscent of the 'DFG' flip in protein kinases. This allosteric binding site can potentially be utilized to generate new classes of species-selective p97 inhibitors. Furthermore, these studies provide a resource for the development of alternative therapeutics for schistosomiasis and a workflow to identify potential drug targets in similar systems with few available molecular tools.

Chromatin remodeling protein BPTF mediates chromatin accessibility at gene promoters in planarian stem cells

BACKGROUND

The regulation of chromatin accessibility is essential in eukaryotic cells as one of several mechanisms that ensure gene activation occurs at appropriate times and in appropriate cell types. Accordingly, mutations in chromatin remodeling proteins are linked to many different developmental disorders and cancers. One example of a chromatin protein that has been linked to both developmental abnormalities and cancer is BPTF/NURF301, the largest subunit of the Nucleosome Remodeling Factor (NuRF) complex. The BPTF subunit is not only important for the formation of NuRF but also helps direct its activity to particular regions of chromatin by preferentially binding histone H3 lysine four trimethylation (H3K4me3). Notably, defects caused by knockdown of bptf in Xenopus embryos mimic those caused by knockdown of wdr5, a core subunit of all H3K4me3 methyltransferase complexes. However, the mechanistic details of how and where BPTF/NuRF is recruited to regulate gene expression vary between studies and have been largely tested in vitro and/or in cultured cells. Improving our understanding of how this chromatin remodeling complex targets specific gene loci and regulates their expression in an organismal context will provide important insight into how pathogenic mutations disrupt its normal, in vivo, cellular functions.

RESULTS

Here, we report our findings on the role of BPTF in maintaining chromatin accessibility and essential function in planarian (Schmidtea mediterranea) stem cells. We find that depletion of planarian BPTF primarily affects accessibility at gene promoters near transcription start sites (TSSs). BPTF-dependent loss of accessibility did not correlate with decreased gene expression when we considered all affected loci. However, we found that genes marked by Set1-dependent H3K4me3, but not MLL1/2-dependent H3K4me3, showed increased sensitivity to the loss of BPTF-dependent accessibility. In addition, knockdown of bptf (Smed-bptf) produces loss-of-function phenotypes similar to those caused by knockdown of Smed-set1.

CONCLUSIONS

The S.mediterranea homolog of NuRF protein BPTF (SMED-BPTF) is essential for normal homeostasis in planarian tissues, potentially through its role in maintaining chromatin accessibility at a specific subset of gene promoters in planarian stem cells. By identifying loci that lose both chromatin accessibility and gene expression after depletion of BPTF, we have identified a cohort of genes that may have important functions in stem cell biology.

Spatial proteomics of Onchocerca volvulus with pleomorphic neoplasms shows local and systemic dysregulation of protein expression

Onchocerca volvulus is the agent of onchocerciasis (river blindness) and targeted by WHO for elimination though mass drug administration with ivermectin. A small percentage of adult female worms develop pleomorphic neoplasms (PN) which occur more frequently after ivermectin treatment. Worms with PN have a lower life expectancy and improved understanding of proteins expressed in PN and their impact on different tissues could help elucidate the mechanisms of macrofilaricidal activity of ivermectin. Within paraffin embedded nodules removed after ivermectin treatment, we detected 24 (5.6%) O. volvulus females with PN. To assess the protein inventory of the PN and identify proteins potentially linked with tumor development, we used laser capture microdissection and highly sensitive mass spectrometry analysis. Three female worms were used to compare the protein profiles of three tissue types (body wall, uterus, and intestine) to the PN, and then to healthy female worms without PN. The healthy females showed all normal embryogenesis. In PN worms, 151 proteins were detected in the body wall, 215 proteins in the intestine, 47 proteins in the uterus and 1,577 proteins in the PN. Only the uterus of one PN female with some stretched intrauterine microfilariae had an elevated number of proteins (601) detectable, while in the uteri of the healthy females 1,710 proteins were detected. Even in tissues that were not directly affected by PN (intestine, body wall), fewer proteins were detected compared to the corresponding tissue of the healthy controls. Immunolocalization of calcium binding protein OvDig-1 (OVOC8391), which was identified through mass spectrometry as one of the proteins with the highest spectral counts in the PN tissue triplicates, allowed us to confirm the results using an independent method. In conclusion we identified proteins that are potentially linked to the development of PN, and systemic dysregulation of protein expression may contribute to worm mortality.

A Krüppel-like factor establishes cellular heterogeneity during schistosome tegumental maintenance

Schistosomes are blood dwelling parasitic flatworms that can survive in the circulation of their human hosts for decades. These parasites possess a unique syncytial skin-like surface tissue known as the tegument that is thought to be uniquely adapted for survival in the blood by mediating evasion of host defenses. Previous studies have shown that cell bodies within the tegumental syncytium are turned over and perpetually replaced by new tegumental cells derived from a pool of somatic stem cells called neoblasts. Thus, neoblast-driven tegumental homeostasis has been suggested to be a key part of the parasite's strategy for long-term survival in the blood. However, the comprehensive set of molecular programs that control the specification of tegumental cells are not defined. To better understand these programs, we characterized a homolog of a Krüppel-like factor 4 (klf4) transcription factor that was identified in previous single-cell RNA sequencing (scRNAseq) studies to be expressed in a putative tegument related lineage (TRL) of Schistosoma mansoni. Here, using a combination of RNAi, coupled with scRNAseq and bulk RNAseq approaches, we show that klf4 is essential for the maintenance of an entire TRL. Loss of this klf4+ TRL resulted in loss of a subpopulation of molecularly unique tegument cells, without altering the total number of mature tegumental cells. Thus, klf4 is critical for regulating the balance between different cell populations within the tegumental progenitor pool and thereby influences tegumental production dynamics and the fine-tuning of the molecular identity of the mature tegument. Understanding the functions of distinct populations of cells within the tegumental syncytium is expected to provide insights into parasite defense mechanisms and new avenues for combatting the disease these worms cause.

Unraveling the mechanisms of benzimidazole resistance in hookworms: A molecular docking and dynamics study

BACKGROUND

Benzimidazole resistance is an emerging challenge among parasitic helminths. It is caused by single nucleotide polymorphisms (SNPs) in specific loci in helminths' β-tubulin genes. Field studies and laboratory investigations reported resistance-associated SNPs in 4 codon locations with 7 allelic variations among hookworms. This study aimed to determine the effects of these mutations on the binding efficiency and behavior of the β-tubulin protein in four hookworm species against four benzimidazole drugs.

METHODS

β-tubulin gene coding sequences of Ancylostoma caninum, A. duodenale, A. ceylanicum, and Necator americanus were retrieved, assessed phylogenetically, and used to construct the 3D structure models of the proteins. The modeled protein structures were verified and edited to contain the reported SNPs: Q134H, F167Y, E198A, E198K, E198V, F200L, and F200Y. Benzimidazole drugs such as albendazole (ABZ), fenbendazole (FBZ), mebendazole (MBZ) and oxfendazole (OBZ) were used as ligands. Molecular docking experiments were performed with the wild-type and mutated proteins. Molecular dynamics simulation assessed the dynamic behavior of the β-tubulin-benzimidazole complex.

RESULTS

In silico docking assessments showed that various amino acid substitutions due to resistance-associated SNPs cause alterations in binding affinities and positions. E198K and Q134H in hookworm β-tubulins substantially weakened the binding affinities and altered the binding positions of benzimidazole drugs. Molecular dynamics analysis revealed that these mutations also caused marked reductions in the binding free energies owing to diminished hydrogen bond contacts with the benzimidazole ligands.

CONCLUSION

The evidence shown herein indicates that mutations at positions 198 and 134 are detrimental to conferring benzimidazole resistance among hookworms. The presence of these mutations may alter the efficacy of pharmacological interventions. Hence, further studies should be conducted to assess their emergence among hookworms in endemic areas with histories of chemotherapy.

Widespread occurrence of benzimidazole resistance single nucleotide polymorphisms in the canine hookworm, Ancylostoma caninum, in Australia

Canine hookworm (Ancylostoma caninum), a gastrointestinal nematode of domestic dogs, principally infects the small intestine of dogs and has the potential to cause zoonotic disease. In greyhounds and pet dogs in the USA, A. caninum has been shown to be resistant to multiple anthelmintics. We conducted a molecular survey of benzimidazole resistance in A. caninum from dogs at veterinary diagnostic centers in Australia and New Zealand. First, we implemented an internal transcribed spacer (ITS)-2 rDNA deep amplicon metabarcoding sequencing approach to ascertain the species of hookworms infecting dogs in the region. Then, we evaluated the frequency of the canonical F167Y and Q134H isotype-1 β-tubulin mutations, which confer benzimidazole resistance, using the same sequencing approach. The most detected hookworm species in diagnostic samples was A. caninum (90%; 83/92); the related Northern hookworm (Uncinaria stenocephala) was identified in 11% (10/92) of the diagnostic samples. There was a single sample with coinfection by A. caninum and U. stenocephala. Both isotype-1 β-tubulin mutations were present in A. caninum, 49% and 67% for Q134H and F167Y, respectively. Mutation F167Y in the isotype-1 β-tubulin mutation was recorded in U. stenocephala for the first known time. Canonical benzimidazole resistance codons 198 and 200 mutations were absent. Egg hatch assays performed on a subset of the A. caninum samples showed significant correlation between 50% inhibitory concentration (IC50) to thiabendazole and F167Y, with an increased IC50 for samples with > 75% F167Y mutation. We detected 14% of dogs with > 75% F167Y mutation in A. caninum. Given that these samples were collected from dogs across various regions of Australia, the present study suggests that benzimidazole resistance in A. caninum is widespread. Therefore, to mitigate the risk of resistance selection and further spread, adoption of a risk assessment-based approach to limit unnecessary anthelmintic use should be a key consideration for future parasite control.

Typing of Echinococcus multilocularis by Region-Specific Extraction and Next-Generation Sequencing of the mitogenome

Background

Infection by the fox tapeworm Echinococcus multilocularis may lead to a severe zoonosis in humans, alveolar echinococcosis, which may be fatal if left untreated. Typing is important to understand the epidemiology of this parasite, yet there is limited knowledge on the microdiversity of E. multilocularis on the local scale, since the typing resolution of established methods is restricted.

Methods

The mitogenome of E. multilocularis was used as the target regions to modify, apply and validate the Region-Specific Extraction (RSE) method in combination with Next-Generation Sequencing (NGS). Single Nucleotide Polymorphisms (SNPs) were detected in the mitochondrial DNA (mtDNA) and analysed bioinformatically. To validate the success and the accuracy of the RSE protocol, the mitogenomes of some E. multilocularis isolates were also analysed by the Whole-Genome Sequencing (WGS).

Results

With the chosen combination of methods, the entire mitogenome (~13 kb) of E. multilocularis could be captured and amplified. The read depth (median ≥ 156X) was sufficient to detect existing SNPs. The comparison of mitogenome sequences extracted by RSE with mitogenome sequences obtained by WGS showed that the accuracy of the RSE method was consistently comparable to direct Whole-Genome Sequencing.

Conclusion

The results demonstrate that the RSE method in combination with NGS is suitable to analyse the microdiversity of E. multilocularis at the whole mitogenome level. For the capture and sequencing of large (several kb) genomic regions of E. multilocularis and other applications, this method can be very helpful.

Identification through a transcriptomic approach of candidate genes involved in the adaptation of the cyst nematode Globodera pallida to the potato resistance factor GpaVvrn

BACKGROUND

Since the banning of chemical products used to control plant-parasitic nematode populations, the use of resistant plants has become the most effective management approach against the potato cyst nematode Globodera pallida. However, some populations, from experimental evolution setups and field samplings, are able to overcome these resistances. Herein, a transcriptomics approach was used to disentangle the mechanisms by which G. pallida adapts to the plant resistant factor GpaVvrn, and to elucidate the functions involved in this adaptation.

RESULTS

Differential gene expression analysis between virulent and avirulent lineages originating from experimental evolution experiments identified candidate genes involved in the adaptation to GpaVvrn. GO enrichment analyses showed that virulent lineages up-regulated genes involved in cell wall destruction and stress response compared to avirulent lineages. In virulent lineages, a set of genes was up-regulated later in the parasitism stages and are thus potentially involved in adaptation. These genes encode effectors of the VAP and SPRYSEC families contributing to the suppression of plant immunity.

CONCLUSION

These results will have a major impact on our understanding of the mechanisms by which nematodes adapt to resistant plants, and will contribute to identify effective and sustainable management strategies.

Eggs-posed: revision of Schistosoma mansoni venom allergen-like proteins unveils new genes and offers new insights into egg-host interactions

BACKGROUND

Venom allergen-like proteins (VALs) are abundant in the excretory-secretory products (ESPs) of numerous parasitic helminths and have been extensively studied for over 30 years because of their potential to interact with host systems. Despite substantial research, however, the precise functions of these proteins remain largely unresolved. Schistosomes, parasites of the circulatory system, are no exception, with 29 SmVAL genes identified in the genome of Schistosoma mansoni to date. The eggs of these parasites, as primary pathogenic agents, interact directly with host tissues and release excretory-secretory products that aid their egress from the host. Although SmVALs have been detected in the egg secretome in the past, direct evidence of their secretion and functional interaction with host molecules has never been demonstrated. These findings fuel the ongoing debate as to whether egg-expressed SmVALs interact with the mammalian host or are rather miracidial proteins synthesized within the egg during larval development.

RESULTS

Based on complete revision of the SmVAL family and an associated robust transcriptomic meta-analysis of gene expression across the life cycle, we show that many of SmVAL genes, including 6 newly identified genes, are expressed in the infective larvae-producing stages (eggs and sporocysts). Following localization of two "egg-specific" SmVAL9 and SmVAL29 did not prove active secretion of these molecules into surrounding tissues but were aligned with miracidial structures interfacing with the molluscan host, specifically the larval surface and penetration glands. Finally, we show the complete lack of interactions between candidate SmVAL proteins and an array of 755 human cell receptors via a state-of-the-art SAVEXIS screen.

CONCLUSIONS

Overall, we conclude that these "egg" SmVALs are not involved in direct host‒parasite interactions in the mammalian host and are rather proteins employed during intermediate host invasion. Our study revisits and updates the SmVAL gene family, highlighting the limitations of in silico protein function predictions while emphasizing the need for up-to-date datasets and tools together with experimental validation in host-parasite interactions. By uncovering the diversity, expression patterns, and interaction dynamics of SmVALs, we open new avenues for understanding host manipulation and reevaluating orthologous proteins in other helminths.

The retinoic acid family-like nuclear receptor SmRAR identified by single-cell transcriptomics of ovarian cells controls oocyte differentiation in Schistosoma mansoni

Studies on transcription regulation in platyhelminth development are scarce, especially for parasitic flatworms. Here, we employed single-cell transcriptomics to identify genes involved in reproductive development in the trematode model Schistosoma mansoni. This parasite causes schistosomiasis, a major neglected infectious disease affecting >240 million people worldwide. The pathology of schistosomiasis is closely associated with schistosome eggs deposited in host organs including the liver. Unlike other trematodes, schistosomes exhibit distinct sexes, with egg production reliant on the pairing-dependent maturation of female reproductive organs. Despite this significance, the molecular mechanisms underlying ovary development and oocyte differentiation remain largely unexplored. Utilizing an organ isolation approach for S. mansoni, we extracted ovaries of paired females followed by single-cell RNA sequencing (RNA-seq) with disassociated oocytes. A total of 1967 oocytes expressing 7872 genes passed quality control (QC) filtering. Unsupervised clustering revealed four distinct cell clusters: somatic, germ cells and progeny, intermediate and late germ cells. Among distinct marker genes for each cluster, we identified a hitherto uncharacterized transcription factor of the retinoic acid receptor family, SmRAR. Functional analyses of SmRAR and associated genes like Smmeiob (meiosis-specific, oligonucleotide/oligosaccharide binding motif (OB) domain-containing) demonstrated their pairing-dependent and ovary-preferential expression and their decisive roles in oocyte differentiation of S. mansoni.

Hookworm genes encoding intestinal excreted-secreted proteins are transcriptionally upregulated in response to the host's immune system

Hookworms are intestinal parasitic nematodes that chronically infect ~500 million people, with reinfection common even after clearance by drugs. How infecting hookworms successfully overcome host protective mechanisms is unclear, but it may involve hookworm proteins that digest host tissues, or counteract the host's immune system, or both. To find such proteins in the zoonotic hookworm Ancylostoma ceylanicum, we identified hookworm genes encoding excreted-secreted (ES) proteins, hookworm genes preferentially expressed in the hookworm intestine, and hookworm genes whose transcription is stimulated by the host immune system. We collected ES proteins from adult hookworms harvested from hamsters; mass spectrometry identified 565 A. ceylanicum genes encoding ES proteins. We also used RNA-seq to identify A. ceylanicum genes expressed both in young adults (12 days post-infection) and in intestinal and non-intestinal tissues dissected from mature adults (19 days post-infection), with hamster hosts that either had normal immune systems or were immunosuppressed by dexamethasone. In adult A. ceylanicum, we observed 1,670 and 1,196 genes with intestine- and non-intestine-biased expression, respectively. Comparing hookworm gene activity in normal versus immunosuppressed hosts, we observed almost no changes of gene activity in 12-day young adults or non-intestinal 19-day adult tissues. However, in intestinal 19-day adult tissues, we observed 1,951 positively immunoregulated genes (upregulated at least two-fold in normal hosts versus immunosuppressed hosts), and 137 genes that were negatively immunoregulated. Thus, immunoregulation was observed primarily in mature adult hookworm intestine directly exposed to host blood; it may include hookworm genes activated in response to the host immune system in order to neutralize the host immune system. We observed 153 ES genes showing positive immunoregulation in 19-day adult intestine; of these genes, 69 had ES gene homologs in the closely related hookworm Ancylostoma caninum, 24 in the human hookworm Necator americanus, and 24 in the more distantly related strongylid parasite Haemonchus contortus. Such a mixture of rapidly evolving and conserved genes could comprise virulence factors enabling infection, provide new targets for drugs or vaccines against hookworm, and aid in developing therapies for autoimmune diseases.

Prognostic markers and molecular pathways in primary colorectal cancer with a high potential of liver metastases: a systems biology approach

Background and purpose

Colorectal cancer (CRC) holds the position of being the third most prevalent cancer and the second primary cause of cancer-related fatalities on a global scale. Approximately 65% of CRC patients survive for 5 years following diagnosis. Metastasis and recurrence frequently occur in half of CRC patients diagnosed at the late stage. This study used bioinformatics analysis to identify key signaling pathways, hub genes, transcription factors, and protein kinases involved in transforming primary CRC with liver metastasis potential. Prognostic markers in CRC were also identified.

Experimental approach

The GSE81582 dataset was re-analyzed to identify differentially expressed genes (DEGs) in early CRC compared to non-tumoral tissues. A protein interaction network (PIN) was constructed, revealing significant modules and hub genes. Prognostic markers, transcription factors, and protein kinases were determined. Boxplot and gene set enrichment analyses were performed.

Findings/Results

This study identified 1113 DEGs in primary CRC compared to healthy controls. PIN analysis revealed 75 hub genes and 8 significant clusters associated with early CRC. The down-regulation of SUCLG2 and KPNA2 correlated with poor prognosis. SIN3A and CDK6 played crucial roles in early CRC transformation, affecting rRNA processing pathways.

Conclusion and implications

This study demonstrated several pathways, biological processes, and genes mediating the malignant transformation of healthy colorectal tissues to primary CRC and may help the prognosis and treatment of patients with early CRC.

High-fidelity annotated triploid genome of the quarantine root-knot nematode, Meloidogyne enterolobii

Root-knot nematodes (RKN) of the genus Meloidogyne are obligatory plant endoparasites that cause substantial economic losses to agricultural production and impact the global food supply. These plant parasitic nematodes belong to the most widespread and devastating genus worldwide, yet few measures of control are available. The most efficient way to control RKN is deployment of resistance genes in plants. However, current resistance genes that control other Meloidogyne species are mostly inefficient on Meloidogyne enterolobii. Consequently, M. enterolobii was listed as a European Union quarantine pest requiring regulation. To gain insight into the molecular characteristics underlying its parasitic success, exploring the genome of M. enterolobii is essential. Here, we report a high-quality genome assembly of M. enterolobii using the high-fidelity long-read sequencing technology developed by Pacific Biosciences, combined with a gap-aware sequence transformer, DeepConsensus. The resulting triploid genome assembly spans 285.4 Mb with 556 contigs, a GC% of 30 ± 0.042 and an N50 value of 2.11 Mb, constituting a useful platform for comparative, population and functional genomics.

Carbon dioxide shapes parasite-host interactions in a human-infective nematode

Skin-penetrating nematodes infect nearly one billion people worldwide. The developmentally arrested infective larvae (iL3s) seek out hosts, invade hosts via skin penetration, and resume development inside the host in a process called activation. Activated infective larvae (iL3as) traverse the host body, ending up as parasitic adults in the small intestine. Skin-penetrating nematodes respond to many chemosensory cues, but how chemosensation contributes to host seeking and intra-host navigation-two crucial steps of the parasite-host interaction-remains poorly understood. Here, we investigate the role of carbon dioxide (CO2) in promoting host seeking and intra-host navigation in the human-infective threadworm Strongyloides stercoralis. We show that S. stercoralis exhibits life-stage-specific behavioral preferences for CO2: iL3s are repelled, non-infective larvae and adults are neutral, and iL3as are attracted. CO2 repulsion in iL3s may prime them for host seeking by stimulating dispersal from host feces, while CO2 attraction in iL3as may direct worms toward high-CO2 areas of the body, such as the lungs and intestine. We also identify sensory neurons that detect CO2; these neurons display CO2-evoked calcium activity, promote behavioral responses to CO2, and express the receptor guanylate cyclase Ss-GCY-9. Finally, we develop an approach for generating stable knockout lines in S. stercoralis and use this approach to show that Ss-gcy-9 is required for CO2-evoked behavioral responses in both iL3s and iL3as. Our results highlight chemosensory mechanisms that shape the interaction between parasitic nematodes and their human hosts and may aid in the design of novel anthelmintics that target the CO2-sensing pathway.

Antibody Responses and the Vaccine Efficacy of Recombinant Glycosyltransferase and Nicastrin Against Schistosoma japonicum

Schistosomiasis is a neglected tropical disease and the second most common parasitic disease after malaria. While praziquantel remains the primary treatment, concerns about drug resistance highlight the urgent need for new drugs and effective vaccines to achieve sustainable control. Previous proteomic studies from our group revealed that the expression of Schistosoma japonicum glycosyltransferase and nicastrin as proteins was higher in single-sex males than mated males, suggesting their critical roles in parasite reproduction and their potential as vaccine candidates. In this study, bioinformatic tools were employed to analyze the structural and functional properties of these proteins, including their signal peptide regions, transmembrane domains, tertiary structures, and protein interaction networks. Recombinant forms of glycosyltransferase and nicastrin were expressed and purified, followed by immunization experiments in BALB/c mice. Immunized mice exhibited significantly elevated specific IgG antibody levels after three immunizations compared to adjuvant and PBS controls. Furthermore, immunization with recombinant glycosyltransferase and nicastrin significantly reduced the reproductive capacity of female worms and liver egg burden, though egg hatchability and adult worm survival were unaffected. These findings demonstrate that recombinant glycosyltransferase and nicastrin are immunogenic and reduce female worm fecundity, supporting their potential as vaccine candidates against schistosomiasis.

Adaptive evolution of stress response genes in parasites aligns with host niche diversity

BACKGROUND

Stress responses are key the survival of parasites and, consequently, also the evolutionary success of these organisms. Despite this importance, our understanding of the evolution of molecular pathways dealing with environmental stressors in parasitic animals remains limited. Here, we tested the link between adaptive evolution of parasite stress response genes and their ecological diversity and species richness. We comparatively investigated antioxidant, heat shock, osmoregulatory, and behaviour-related genes (foraging) in two model parasitic flatworm lineages with contrasting ecological diversity, Cichlidogyrus and Kapentagyrus (Platyhelminthes: Monopisthocotyla), through whole-genome sequencing of 11 species followed by in silico exon bait capture as well as phylogenetic and codon analyses.

RESULTS

We assembled the sequences of 48 stress-related genes and report the first foraging (For) gene orthologs in flatworms. We found duplications of heat shock (Hsp) and oxidative stress genes in Cichlidogyrus compared to Kapentagyrus. We also observed positive selection patterns in genes related to mitochondrial protein import (Hsp) and behaviour (For) in species of Cichlidogyrus infecting East African cichlids-a host lineage under adaptive radiation. These patterns are consistent with a potential adaptation linked to a co-radiation of these parasites and their hosts. Additionally, the absence of cytochrome P450 and kappa and sigma-class glutathione S-transferases in monogenean flatworms is reported, genes considered essential for metazoan life.

CONCLUSIONS

This study potentially identifies the first molecular function linked to a flatworm radiation. Furthermore, the observed gene duplications and positive selection indicate the potentially important role of stress responses for the ecological adaptation of parasite species.

Evaluating the Immunoprotective and Diagnostic Potential of Schistosoma mansoni Epitopes from Sm050890 and Sm141290 Proteins Identified Through Reverse Vaccinology

PURPOSE

Schistosomiasis remains a parasitic disease affecting millions of people worldwide, requiring interventions like vaccination. In previous work, our group used reverse vaccinology to identify two epitopes from the Schistosoma mansoni proteins, Sm050890 (44-58) and Sm141290 (225-239). This study evaluated the immune response profile and protection induced by peptides, as a mixture of immunogens, in murine vaccination trials. Additionally, the diagnostic potential of these peptides was assessed on immunoassays.

METHODS

Mice were immunized with a formulation containing the mixture of the peptides, subsequently infected, and perfused for worm burden recovery and quantification. Liver and blood samples from animals were used to evaluate the effect of immunization on the formation of granulomas and specific anti-peptide antibodies (IgG). Additionally, cytokine measurement was performed in splenocyte cultures from immunized mice, and peripheral blood serum from individuals infected with S. mansoni was used to assess the recognition of the peptides by IgG antibodies.

RESULTS

The vaccine stimulated an increase in the production of IgG and IgG2c antibodies, associated with a significant reduction of 44 - 29% in worm burden. Although the vaccine did not reduce liver pathology, it enhanced the production of IFN-γ while decreasing IL-10 production by splenocytes. Furthermore, the peptides Sm050890 (44-58) and Sm141290 (225-239) were not recognized by IgG antibodies in the serum from infected individuals.

CONCLUSION

Overall, our data suggest that the peptides Sm050890 (44-58) and Sm141290 (225-239) are promising vaccine candidates against schistosomiasis and can be used to compose a multiepitope/chimeric vaccine in future studies.

UniProt: the Universal Protein Knowledgebase in 2025

The aim of the UniProt Knowledgebase (UniProtKB; https://www.uniprot.org/) is to provide users with a comprehensive, high-quality and freely accessible set of protein sequences annotated with functional information. In this publication, we describe ongoing changes to our production pipeline to limit the sequences available in UniProtKB to high-quality, non-redundant reference proteomes. We continue to manually curate the scientific literature to add the latest functional data and use machine learning techniques. We also encourage community curation to ensure key publications are not missed. We provide an update on the automatic annotation methods used by UniProtKB to predict information for unreviewed entries describing unstudied proteins. Finally, updates to the UniProt website are described, including a new tab linking protein to genomic information. In recognition of its value to the scientific community, the UniProt database has been awarded Global Core Biodata Resource status.

Ivermectin induces oxidative stress and mitochondrial damage in Haemonchus contortus

Haemonchus contortus poses a severe threat to livestock, leading to substantial economic losses. The extensive use of ivermectin (IVM), an effective anthelmintic drug, has resulted in the development of resistant strains, reducing the drug's effectiveness and making disease management more challenging. Hence, understanding the underlying mechanisms of IVM resistance is crucial. IVM has been shown to induce apoptosis and oxidative stress imbalance in cancer cells, thereby inhibiting tumor growth; however, whether there is a similar impact on H. contortus has not been well-documented. In this study, the biological functions of previously identified resistance-associated genes were investigated, revealing their involvement in apoptosis and oxidative stress. Experiments were conducted to compare cell apoptosis, mitochondrial function, and oxidative stress markers in IVM-treated and untreated sensitive and resistant strains isolated from the field. The findings demonstrated that mitochondrial dysfunction and increased apoptosis were induced by IVM treatment in sensitive strains, as indicated by elevated mitochondrial membrane potential and apoptosis rates. Conversely, stable mitochondrial function and apoptosis levels under IVM exposure were observed in resistant strains. These results provide novel insights into the resistance mechanisms of H. contortus and offer a theoretical basis for future research on overcoming IVM resistance.

Proteomic characterization and comparison of the infective and adult life stage secretomes from Necator americanus and Ancylostoma ceylanicum

More than 470 million people globally are infected with the hookworms Ancylostoma ceylanicum and Necator americanus, resulting in an annual loss of 2.1 to 4 million disability-adjusted-life-years. Current infection management approaches are limited by modest drug efficacy, the costs associated with frequent mass drug administration campaigns, and the risk of reinfection and burgeoning drug resistance. Subunit vaccines based on proteins excreted and secreted (ES) by hookworms that reduce worm numbers and associated disease burden are a promising management strategy to overcome these limitations. However, studies on the ES proteomes of hookworms have mainly described proteins from the adult life stage which may preclude the opportunity to target the infective larva. Here, we employed high resolution mass spectrometry to identify 103 and 57 ES proteins from the infective third larvae stage (L3) as well as 106 and 512 ES proteins from the adult N. americanus and A. ceylanicum respectively. Comparisons between these developmental stages identified 91 and 41 proteins uniquely expressed in the L3 ES products of N. americanus and A. ceylanicum, respectively. We characterized these proteins based on functional annotation, KEGG pathway analysis, InterProScan signature and gene ontology. We also performed reciprocal BLAST analysis to identify orthologs across species for both the L3 and adult stages and identified five orthologous proteins in both life stages and 15 proteins that could be detected only in the L3 stage of both species. Last, we performed a three-way reciprocal BLAST on the L3 proteomes from both hookworm species together with a previously reported L3 proteome from the rodent hookworm Nippostrongylus brasiliensis, and identified eight L3 proteins that could be readily deployed for testing using well established rodent models. This novel characterization of L3 proteins and taxonomic conservation across hookworm species provides a raft of potential candidates for vaccine discovery for prevention of hookworm infection and disease.

Exploitation of phylum-spanning omics resources reveals complexity in the nematode FLP signalling system and provides insights into flp-gene evolution

BACKGROUND

Parasitic nematodes significantly undermine global human and animal health and productivity. Parasite control is reliant on anthelmintic administration however over-use of a limited number of drugs has resulted in escalating parasitic nematode resistance, threatening the sustainability of parasite control and underscoring an urgent need for the development of novel therapeutics. FMRFamide-like peptides (FLPs), the largest family of nematode neuropeptides, modulate nematode behaviours including those important for parasite survival, highlighting FLP receptors (FLP-GPCRs) as appealing putative novel anthelmintic targets. Advances in omics resources have enabled the identification of FLPs and neuropeptide-GPCRs in some parasitic nematodes, but remaining gaps in FLP-ligand libraries hinder the characterisation of receptor-ligand interactions, which are required to drive the development of novel control approaches.

RESULTS

In this study we exploited recent expansions in nematode genome data to identify 2143 flp-genes in > 100 nematode species across free-living, entomopathogenic, plant, and animal parasitic lifestyles and representing 7 of the 12 major nematode clades. Our data reveal that: (i) the phylum-spanning flps, flp-1, -8, -14, and - 18, may be representative of the flp profile of the last common ancestor of nematodes; (ii) the majority of parasitic nematodes have a reduced flp complement relative to free-living species; (iii) FLP prepropeptide architecture is variable within and between flp-genes and across nematode species; (iv) FLP prepropeptide signatures facilitate flp-gene discrimination; (v) FLP motifs display variable length, amino acid sequence, and conservation; (vi) CLANS analysis provides insight into the evolutionary history of flp-gene sequelogues and reveals putative flp-gene paralogues and, (vii) flp expression is upregulated in the infective larval stage of several nematode parasites.

CONCLUSIONS

These data provide the foundation required for phylum-spanning FLP-GPCR deorphanisation screens in nematodes to seed the discovery and development of novel parasite control approaches.

Glycolipidomics of Liver Flukes and Host Tissues during Fascioliasis: Insights from Mass Spectrometry Imaging

Fascioliasis, a zoonotic disease caused by liver flukes of the genus Fasciola, poses significant health threats to both humans and livestock. While some infections remain asymptomatic, others can lead to fatal outcomes, particularly during the acute phase characterized by the migration of immature parasites causing severe liver damage. Through the combination of data acquired via high-spatial-resolution atmospheric-pressure scanning microprobe matrix-assisted laser desorption/ionization mass spectrometry imaging (AP-SMALDI MSI) and nanohydrophilic interaction chromatography tandem mass spectrometry, we investigated glycosphingolipids (GSLs) in both adult and immature parasite stages as well as the host liver and bile duct to unravel the intricacies of the host-pathogen interplay and associated pathology. Several GSLs showed characteristic distribution patterns within the parasite depending on the fatty acid composition of their ceramides, notably including GSLs carrying very long-chain fatty acids. Additionally, GSL compositions within the tegument of immature versus adult parasites varied, suggestive of tissue remodeling upon maturation. AP-SMALDI MSI further enabled the identification of GSLs potentially involved in in vivo interactions between the host and immature parasites. Moreover, our experiments unveiled alterations in other lipid classes during Fasciola infection, providing a broader understanding of lipidomic changes associated with the disease. Collectively, our findings contribute to a deeper comprehension of the molecular intricacies underlying fascioliasis, with a specific focus on GSLs.

Inference of essential genes in Brugia malayi and Onchocerca volvulus by machine learning and the implications for discovering new interventions

Detailed explorations of the model organisms Caenorhabditis elegans (elegant worm) and Drosophila melanogaster (vinegar fly) have substantially improved our knowledge and understanding of biological processes and pathways in metazoan organisms. Extensive functional genomic and multi-omic data sets have enabled the discovery and characterisation of 'essential' genes that are critical for the survival of these organisms. Recently, we showed that a machine learning (ML)-based pipeline could be utilised to predict essential genes in both C. elegans and D. melanogaster using features from DNA, RNA, protein and/or cellular data or associated information. As these distantly-related species are within the Ecdysozoa, we hypothesised that this approach could be suited for non-model organisms within the same group (phylum) of protostome animals. In the present investigation, we cross-predicted essential genes within the phylum Nematoda - between C. elegans and the parasitic filarial nematodes Brugia malayi and Onchocerca volvulus, and then ranked and prioritised these genes. Highly ranked genes were linked to key biological pathways or processes, such as ribosome biogenesis, translation and RNA processing, and were expressed at relatively high levels in the germline, gonad, hypodermis and/or nerves. The present in silico workflow is hoped to expedite the identification of drug targets in parasitic organisms for subsequent experimental validation in the laboratory.

Analysis of Haemonchus embryos at single cell resolution identifies two eukaryotic elongation factors as intervention target candidates

Advances in single cell technologies are allowing investigations of a wide range of biological processes and pathways in animals, such as the multicellular model organism Caenorhabditis elegans - a free-living nematode. However, there has been limited application of such technology to related parasitic nematodes which cause major diseases of humans and animals worldwide. With no vaccines against the vast majority of parasitic nematodes and treatment failures due to drug resistance or inefficacy, new intervention targets are urgently needed, preferably informed by a deep understanding of these nematodes' cellular and molecular biology - which is presently lacking for most worms. Here, we created the first single cell atlas for an early developmental stage of Haemonchus contortus - a highly pathogenic, C. elegans-related parasitic nematode. We obtained and curated RNA sequence (snRNA-seq) data from single nuclei from embryonating eggs of H. contortus (150,000 droplets), and selected high-quality transcriptomic data for > 14,000 single nuclei for analysis, and identified 19 distinct clusters of cells. Guided by comparative analyses with C. elegans, we were able to reproducibly assign seven cell clusters to body wall muscle, hypodermis, neuronal, intestinal or seam cells, and identified eight genes that were transcribed in all cell clusters/types, three of which were inferred to be essential in H. contortus. Two of these genes (i.e. Hc-eef-1A and Hc-eef1G), coding for eukaryotic elongation factors (called Hc-eEF1A and Hc-eEF1G), were also demonstrated to be transcribed and expressed in all key developmental stages of H. contortus. Together with these findings, sequence- and structure-based comparative analyses indicated the potential of Hc-eEF1A and/or Hc-eEF1G as intervention targets within the protein biosynthesis machinery of H. contortus. Future work will focus on single cell studies of all key developmental stages and tissues of H. contortus, and on evaluating the suitability of the two elongation factor proteins as drug targets in H. contortus and related nematodes, with a view to finding new nematocidal drug candidates.

Dissection of the synthesis of polyunsaturated fatty acids in nematodes and Collembola of the soil fauna

It is becoming increasingly clear that not only unicellular, photoautotrophic eukaryotes, plants, and fungi, but also invertebrates are capable of synthesizing ω3 long-chain polyunsaturated fatty acids (LC-PUFA) de novo. However, the distribution of this anabolic capacity among different invertebrate groups and its implementation at the gene and protein level are often still unknown. This study investigated the PUFA pathways in common soil fauna, i.e. two nematode and two Collembola species. Of these, one species each (Panagrellus redivivus, Folsomia candida) was assumed to produce ω3 LC-PUFA de novo, while the others (Acrobeloides bodenheimeri, Isotoma caerulea) were supposed to be unable to do so. A highly labeled oleic acid (99 % 13C) was supplemented and the isotopic signal was used to trace its metabolic path. All species followed the main pathway of lipid biosynthesis. However, in A. bodenheimeri this terminated at arachidonic acid (ω6 PUFA), whereas the other three species continued the pathway to eicosapentaenoic acid (ω3 PUFA), including I. caerulea. For the nematode P. redivivus the identification and functional characterization of four new fatty acid desaturase (FAD) genes was performed. These genes encode the FAD activities Δ9, Δ6, and Δ5, respectively. Additionally, the Δ12 desaturase was analyzed, yet the observed activity of an ω3 FAD could not be attributed to a coding gene. In the Collembola F. candida, 11 potential first desaturases (Δ9) and 13 front-end desaturases (Δ6 or Δ5 FADs) have been found. Further sequence analysis indicates the presence of omega FADs, specifically Δ12, which are likely derived from Δ9 FADs.

The N terminus-only (trans) function of the adhesion G protein-coupled receptor latrophilin-1 controls multiple processes in reproduction of Caenorhabditis elegans

Adhesion G protein-coupled receptors are unique molecules. They are able to transmit classical signals via G protein activation as well as mediate functions solely through their extracellular N termini, completely independently of the seven transmembrane helices domain and the C terminus. This dual mode of action is highly unusual for G protein-coupled receptors and allows for a plethora of possible cellular consequences. However, the physiological implications and molecular details of this N terminus-mediated signaling are poorly understood. Here, we show that several distinct seven transmembrane helices domain-independent/trans functions of the adhesion G protein-coupled receptor latrophilin homolog latrophilin-1 in the nematode Caenorhabditis elegans together regulate reproduction: sperm guidance, ovulation, and germ cell apoptosis. In these contexts, the receptor elicits its functions in a noncell autonomous manner. The functions might be realized through alternative splicing of the receptor specifically generating N terminus-only variants. Thus, our findings shed light on the versatility of seven transmembrane helices domain-independent/N terminus-only/trans functions of adhesion G protein-coupled receptor and discuss possible molecular details.

The Raf/LIN-45 C-terminal distal tail segment negatively regulates signaling in Caenorhabditis elegans

Raf protein kinases act as Ras-GTP sensing components of the ERK signal transduction pathway in animal cells, influencing cell proliferation, differentiation, and survival. In humans, somatic and germline mutations in the genes BRAF and RAF1 are associated with malignancies and developmental disorders. Recent studies shed light on the structure of activated Raf, a heterotetramer consisting of Raf and 14-3-3 dimers, and raised the possibility that a Raf C-terminal distal tail segment (DTS) regulates activation. We investigated the role of the DTS using the Caenorhabditis elegans Raf ortholog lin-45. Truncations removing the DTS strongly enhanced lin-45(S312A), a weak gain-of-function allele equivalent to RAF1 mutations found in patients with Noonan Syndrome. We genetically defined three elements of the LIN-45 DTS, which we termed the active site binding sequence (ASBS), the KTP motif, and the aromatic cluster. In the context of lin-45(S312A), the mutation of each of these elements enhanced activity. We used AlphaFold to predict DTS protein interactions for LIN-45, fly Raf, and human BRAF within the activated heterotetramer complex. We propose the following distinct functions for the LIN-45 DTS elements: (1) the ASBS binds the kinase active site as an inhibitor; (2) phosphorylation of the KTP motif modulates the DTS-kinase domain interaction; and (3) the aromatic cluster anchors the DTS in an inhibitory conformation. Human RASopathy-associated variants in BRAF affect residues of the DTS, consistent with these predictions. This work establishes that the Raf/LIN-45 DTS negatively regulates signaling in C. elegans and provides a model for its function in other Raf proteins.

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.

Muscular remodeling and anteroposterior patterning during tapeworm segmentation

BACKGROUND

Tapeworms are parasitic flatworms that independently evolved a segmented body plan, historically confounding comparisons with other animals. Anteroposterior (AP) patterning in free-living flatworms and in tapeworm larvae is associated with canonical Wnt signaling and positional control genes (PCGs) are expressed by their musculature in regionalized domains along the AP axis. Here, we extend investigations of PCG expression to the adult of the mouse bile-duct tapeworm Hymenolepis microstoma, focusing on the growth zone of the neck region and the initial establishment of segmental patterning.

RESULTS

We show that the adult musculature includes new, segmental elements that first appear in the neck and that the spatial patterns of Wnt factors are consistent with expression by muscle cells. Wnt factor expression is highly regionalized and becomes AP-polarized in segments, marking them with axes in agreement with the polarity of the main body axis, while the transition between the neck and strobila is specifically demarcated by the expression domain of a Wnt11 paralog.

CONCLUSION

We suggest that segmentation could originate in the muscular system and participate in patterning the AP axis through regional and polarized expression of PCGs, akin to the gene regulatory networks employed by free-living flatworms and other animals.

Chromosome-contiguous genome for the Haecon-5 strain of Haemonchus contortus reveals marked genetic variability and enables the discovery of essential gene candidates

Millions of livestock animals worldwide are infected with the haematophagous barber's pole worm, Haemonchus contortus, the aetiological agent of haemonchosis. Despite the major significance of this parasite worldwide and its widespread resistance to current treatments, the lack of a high-quality genome for the well-defined strain of this parasite from Australia, called Haecon-5, has constrained research in a number of areas including host-parasite interactions, drug discovery and population genetics. To enable research in these areas, we report here a chromosome-contiguous genome (∼280 Mb) for Haecon-5 with high-quality models for 19,234 protein-coding genes. Comparative genomic analyses show significant genomic similarity (synteny) with a UK strain of H. contortus, called MHco3(ISE).N1 (abbreviated as "ISE"), but we also discover marked differences in genomic structure/gene arrangements, distribution of nucleotide variability (single nucleotide polymorphisms (SNPs) and indels) and orthology between Haecon-5 and ISE. We used the genome and extensive transcriptomic resources for Haecon-5 to predict a subset of essential single-copy genes employing a "cross-species" machine learning (ML) approach using a range of features from nucleotide/protein sequences, protein orthology, subcellular localisation, single-cell RNA-seq and/or histone methylation data available for the model organisms Caenorhabditis elegans and Drosophila melanogaster. From a set of 1,464 conserved single copy genes, transcribed in key life-cycle stages of H. contortus, we identified 232 genes whose homologs have critical functions in C. elegans and/or D. melanogaster, and prioritised 10 of them for further characterisation; nine of the 10 genes likely play roles in neurophysiological processes, germline, hypodermis and/or respiration, and one is an unknown (orphan) gene for which no detailed functional information exists. Future studies of these genes/gene products are warranted to elucidate their roles in parasite biology, host-parasite interplay and/or disease. Clearly, the present Haecon-5 reference genome and associated resources now underpin a broad range of fundamental investigations of H. contortus and could assist in accelerating the discovery of novel intervention targets and drug candidates to combat haemonchosis.

cGMP-dependent pathway and a GPCR kinase are required for photoresponse in the nematode Pristionchus pacificus

Light sensing is a critical function in most organisms and is mediated by photoreceptor proteins and phototransduction. Although most nematodes lack eyes, some species exhibit phototaxis. In the nematode Caenorhabditis elegans, the unique photoreceptor protein Cel-LITE-1, its downstream G proteins, and cyclic GMP (cGMP)-dependent pathways are required for phototransduction. However, the mechanism of light-sensing in other nematodes remains unknown. To address this question, we used the nematode Pristionchus pacificus, which was established as a satellite model organism for comparison with C. elegans. Similar to C. elegans, illumination with short-wavelength light induces avoidance behavior in P. pacificus. Opsin, cryptochrome/photolyase, and lite-1 were not detected in the P. pacificus genome using orthology and domain prediction-based analyses. To identify the genes related to phototransduction in P. pacificus, we conducted forward genetic screening for light-avoidance behavior and isolated five light-unresponsive mutants. Whole-genome sequencing and genetic mapping revealed that the cGMP-dependent pathway and Ppa-grk-2, which encodes a G protein-coupled receptor kinase (GRK) are required for light avoidance. Although the cGMP-dependent pathway is conserved in C. elegans phototransduction, GRK is not necessary for light avoidance in C. elegans. This suggests similarities and differences in light-sensing mechanisms between the two species. Using a reverse genetic approach, we showed that gamma-aminobutyric acid (GABA) and glutamate were involved in light avoidance. Through reporter analysis and suppression of synapse transmission, we identified candidate photosensory neurons. These findings advance our understanding of the diversity of phototransduction in nematodes even in the absence of eyes.

Genome Assembly of the Nematode Rhabditoides Inermis From a Complex Microbial Community

Free-living nematodes such as Caenorhabditis elegans and Pristionchus pacificus are powerful model systems for linking specific traits to their underlying genetic basis. To trace the evolutionary history of specific traits or genes, a robust phylogenomic framework is indispensable. In the context of the nematode family Diplogastridae to which P. pacificus belongs, the identity of a sister group has long been debated. In this work, we generated a pseudochromosome level genome assembly of the nematode Rhabditoides inermis, which has previously been proposed as the sister taxon. The genome was assembled from a complex microbial community that is stably associated with R. inermis isolates and that consists of multiple bacteria and a fungus, which we identified as a strain of Vanrija albida. The R. inermis genome spans 173.5Mb that are largely assembled into five pseudochromosomes. This chromosomal configuration likely arose from two recent fusions of different Nigon elements. Phylogenomic analysis did not support a sister group relationship between R. inermis and diplogastrids, but rather supports a sister group relationship between the monophyletic Diplogastridae and a group of genera of Rhabditidae including C. elegans and R. inermis. Thus, our work addresses for the first time the long lasting question about the sister group to diplogastrids at the phylogenomic level and provides with the genomes of R. inermis and the associated fungus V. albida valuable resources for future genomic comparisons.