Mapping resistance-associated anthelmintic interactions in the model nematode Caenorhabditis elegans.

Parasitic nematodes infect billions of people and are mainly controlled by anthelmintic mass drug administration (MDA). While there are growing efforts to better understand mechanisms of anthelmintic resistance in human and animal populations, it is unclear how resistance mechanisms that alter susceptibility to one drug affect the interactions and efficacy of drugs used in combination. Mutations that alter drug permeability across primary nematode barriers have been identified as potential resistance mechanisms using the model nematode Caenorhabditis elegans. We leveraged high-throughput assays in this model system to measure altered anthelmintic susceptibility in response to genetic perturbations of potential cuticular, amphidial, and alimentary routes of drug entry. Mutations in genes associated with these tissue barriers differentially altered susceptibility to the major anthelmintic classes (macrocyclic lactones, benzimidazoles, and nicotinic acetylcholine receptor agonists) as measured by animal development. We investigated two-way anthelmintic interactions across C. elegans genetic backgrounds that confer resistance or hypersensitivity to one or more drugs. We observe that genetic perturbations that alter susceptibility to a single drug can shift the drug interaction landscape and lead to the appearance of novel synergistic and antagonistic interactions. This work establishes a framework for investigating combinatorial therapies in model nematodes that can potentially be translated to amenable parasite species.

A high-throughput sensory assay for parasitic and free-living nematodes.

Sensory pathways first elucidated in Caenorhabditis elegans are conserved across free-living and parasitic nematodes, even though each species responds to a diverse array of compounds. Most nematode sensory assays are performed by tallying observations of worm behavior on two-dimensional planes using agarose plates. These assays have been successful in the study of volatile sensation but are poorly suited for investigation of water-soluble gustation or parasitic nematodes without a free-living stage. In contrast, gustatory assays tend to be tedious, often limited to the manipulation of a single individual at a time. We have designed a nematode sensory assay using a microfluidics device that allows for the study of gustation in a 96-well, three-dimensional environment. This device is suited for free-living worms and parasitic worms that spend their lives in an aqueous environment, and we have used it to show that ivermectin inhibits the gustatory ability of vector-borne parasitic nematodes. Insight box Nematodes are powerful model organisms for understanding the sensory biology of multicellular eukaryotes, and many parasitic species cause disease in humans. Simple sensory assays performed on agarose plates have been the bedrock for establishing the neuronal, genetic, and developmental foundations for many sensory modalities in nematodes. However, these classical assays are poorly suited for translational movement of many parasitic nematodes and the sensation of water-soluble molecules (gustation). We have designed a device for high-throughput nematode sensory assays in a gel matrix. This 'gustatory microplate' is amenable to several species and reveals novel responses by free-living and parasitic nematodes to cues and drugs.

Extracellular vesicles secreted by Brugia malayi microfilariae modulate the melanization pathway in the mosquito host.

Vector-borne, filarial nematode diseases cause significant disease burdens in humans and domestic animals worldwide. Although there is strong direct evidence of parasite-driven immunomodulation of mammalian host responses, there is less evidence of parasite immunomodulation of the vector host. We have previously reported that all life stages of Brugia malayi, a filarial nematode and causative agent of Lymphatic filariasis, secrete extracellular vesicles (EVs). Here we investigate the immunomodulatory effects of microfilariae-derived EVs on the vector host Aedes aegypti. RNA-seq analysis of an Ae. aegypti cell line treated with B. malayi microfilariae EVs showed differential expression of both mRNAs and miRNAs. AAEL002590, an Ae. aegypti gene encoding a serine protease, was shown to be downregulated when cells were treated with biologically relevant EV concentrations in vitro. Injection of adult female mosquitoes with biologically relevant concentrations of EVs validated these results in vivo, recapitulating the downregulation of AAEL002590 transcript. This gene was predicted to be involved in the mosquito phenoloxidase (PO) cascade leading to the canonical melanization response and correspondingly, both suppression of this gene using RNAi and parasite EV treatment reduced PO activity in vivo. Our data indicate that parasite-derived EVs interfere with critical immune responses in the vector host, including melanization.

A high-throughput nematode sensory assay reveals an inhibitory effect of ivermectin on parasite gustation.

Sensory pathways first elucidated in Caenorhabditis elegans are conserved across free-living and parasitic nematodes, even though each species responds to a diverse array of compounds. Most nematode sensory assays are performed by tallying observations of worm behavior on two-dimensional planes using agarose plates. These assays have been successful in the study of volatile sensation but are poorly suited for investigation of water-soluble gustation or parasitic nematodes without a free-living stage. In contrast, gustatory assays tend to be tedious, often limited to the manipulation of a single individual at a time. We have designed a nematode sensory assay using a microfluidics device that allows for the study of gustation in a 96-well, three-dimensional environment. This device is suited for free-living worms and parasitic worms that spend their lives in an aqueous environment, and we have used it to show that ivermectin inhibits the gustatory ability of vector-borne parasitic nematodes.

Mapping resistance-associated anthelmintic interactions in the model nematode Caenorhabditis elegans.

Parasitic nematodes infect billions of people and are mainly controlled by anthelmintic mass drug administration (MDA). While there are growing efforts to better understand mechanisms of anthelmintic resistance in human and animal populations, it is unclear how resistance mechanisms that alter susceptibility to one drug affect the interactions and efficacy of drugs used in combination. Mutations that alter drug permeability across primary nematode barriers have been identified as potential resistance mechanisms using the model nematode Caenorhabditis elegans. We leveraged high-throughput assays in this model system to measure altered anthelmintic susceptibility in response to genetic perturbations of potential cuticular, amphidial, and alimentary routes of drug entry. Mutations in genes associated with these tissue barriers differentially altered susceptibility to the major anthelmintic classes (macrocyclic lactones, benzimidazoles, and nicotinic acetylcholine receptor agonists) as measured by animal development. We investigated two-way anthelmintic interactions across C. elegans genetic backgrounds that confer resistance or hypersensitivity to one or more drugs. We observe that genetic perturbations that alter susceptibility to a single drug can shift the drug interaction landscape and lead to the appearance of novel synergistic and antagonistic interactions. This work establishes a framework for investigating combinatorial therapies in model nematodes that can potentially be translated to amenable parasite species.

Phenotypic Profiling of Macrocyclic Lactones on Parasitic Schistosoma Flatworms.

Macrocyclic lactones are front-line therapies for parasitic roundworm infections; however, there are no comprehensive studies on the activity of this drug class against parasitic flatworms. Ivermectin is well known to be inactive against flatworms. However, the structure-activity relationship of macrocyclic lactones may vary across phyla, and it is entirely possible other members of this drug class do in fact show antiparasitic activity on flatworms. For example, there are several reports hinting at the anti-schistosomal activity of doramectin and moxidectin. To explore this class further, we developed an automated imaging assay combined with measurement of lactate levels from worm media. This assay was applied to the screening of 21 macrocyclic lactones (avermectins, milbemycins, and others such as spinosyns) against adult schistosomes. These in vitro assays identified several macrocyclic lactones (emamectin, milbemycin oxime, and the moxidectin metabolite 23-ketonemadectin) that caused contractile paralysis and lack of lactate production. Several of these were also active against miracidia, which infect the snail intermediate host. Hits prioritized from these in vitro assays were administered to mice harboring patent schistosome infections. However, no reduction in worm burden was observed. Nevertheless, these data show the utility of a multiplexed in vitro screening platform to quantitatively assess drug action and exclude inactive compounds from a chemical series before proceeding to in vivo studies. While the prototypical macrocyclic lactone ivermectin displays minimal activity against adult Schistosoma mansoni, this family of compounds does contain schistocidal compounds which may serve as a starting point for development of new anti-flatworm chemotherapies.

Annotation of G-Protein Coupled Receptors in the Genomes of Parasitic Blood Flukes.

Infection with Schistosoma parasitic flatworms ( Schistosoma haematobium, Schistosoma mansoni and Schistosoma japonicum ) causes the neglected tropical disease schistosomiasis. There is a need to identify new chemotherapies to treat these parasites, and G-protein coupled receptors (GPCRs) are a logical druggable targets to explore given they control key aspects of schistosome biology such as neuromuscular function and reproduction. Updated chromosome level genome assemblies for each of the three major species have recently been released. However, studies on these GPCRs require accurate, updated genome annotations. Here, we have re-annotated the GPCRs present in each of the three major schistosome species.

Multivariate chemogenomic screening prioritizes new macrofilaricidal leads.

Development of direct acting macrofilaricides for the treatment of human filariases is hampered by limitations in screening throughput imposed by the parasite life cycle. In vitro adult screens typically assess single phenotypes without prior enrichment for chemicals with antifilarial potential. We developed a multivariate screen that identified dozens of compounds with submicromolar macrofilaricidal activity, achieving a hit rate of >50% by leveraging abundantly accessible microfilariae. Adult assays were multiplexed to thoroughly characterize compound activity across relevant parasite fitness traits, including neuromuscular control, fecundity, metabolism, and viability. Seventeen compounds from a diverse chemogenomic library elicited strong effects on at least one adult trait, with differential potency against microfilariae and adults. Our screen identified five compounds with high potency against adults but low potency or slow-acting microfilaricidal effects, at least one of which acts through a novel mechanism. We show that the use of microfilariae in a primary screen outperforms model nematode developmental assays and virtual screening of protein structures inferred with deep learning. These data provide new leads for drug development, and the high-content and multiplex assays set a new foundation for antifilarial discovery.

Pharmacological Profiling of a Brugia malayi Muscarinic Acetylcholine Receptor as a Putative Antiparasitic Target.

The diversification of anthelmintic targets and mechanisms of action will help ensure the sustainable control of nematode infections in response to the growing threat of drug resistance. G protein-coupled receptors (GPCRs) are established drug targets in human medicine but remain unexploited as anthelmintic substrates despite their important roles in nematode neuromuscular and physiological processes. Bottlenecks in exploring the druggability of parasitic nematode GPCRs include a limited helminth genetic toolkit and difficulties establishing functional heterologous expression. In an effort to address some of these challenges, we profile the function and pharmacology of muscarinic acetylcholine receptors in the human parasite Brugia malayi, an etiological agent of human lymphatic filariasis. While acetylcholine-gated ion channels are intensely studied as targets of existing anthelmintics, comparatively little is known about metabotropic receptor contributions to parasite cholinergic signaling. Using multivariate phenotypic assays in microfilariae and adults, we show that nicotinic and muscarinic compounds disparately affect parasite fitness traits. We identify a putative G protein-linked acetylcholine receptor of B. malayi (Bma-GAR-3) that is highly expressed across intramammalian life stages and adapt spatial RNA in situ hybridization to map receptor transcripts to critical parasite tissues. Tissue-specific expression of Bma-gar-3 in Caenorhabditis elegans (body wall muscle, sensory neurons, and pharynx) enabled receptor deorphanization and pharmacological profiling in a nematode physiological context. Finally, we developed an image-based feeding assay as a reporter of pharyngeal activity to facilitate GPCR screening in parasitized strains. We expect that these receptor characterization approaches and improved knowledge of GARs as putative drug targets will further advance the study of GPCR biology across medically important nematodes.

wrmXpress: A modular package for high-throughput image analysis of parasitic and free-living worms.

Advances in high-throughput and high-content imaging technologies require concomitant development of analytical software capable of handling large datasets and generating relevant phenotypic measurements. Several tools have been developed to analyze drug response phenotypes in parasitic and free-living worms, but these are siloed and often limited to specific instrumentation, worm species, and single phenotypes. No unified tool exists to analyze diverse high-content phenotypic imaging data of worms and provide a platform for future extensibility. We have developed wrmXpress, a unified framework for analyzing a variety of phenotypes matched to high-content experimental assays of free-living and parasitic nematodes and flatworms. We demonstrate its utility for analyzing a suite of phenotypes, including motility, development/size, fecundity, and feeding, and establish the package as a platform upon which to build future custom phenotypic modules. We show that wrmXpress can serve as an analytical workhorse for anthelmintic screening efforts across schistosomes, filarial nematodes, and free-living model nematodes and holds promise for enabling collaboration among investigators with diverse interests.

Making sense of sensory behaviors in vector-borne helminths.

Migrations performed by helminths are impressive and diverse, and accumulating evidence shows that many are controlled by sophisticated sensory programs. The migrations of vector-borne helminths are particularly complex, requiring precise, stage-specific regulation. We review the contrasting states of knowledge on snail-borne schistosomes and mosquito-borne filarial nematodes. Rich observational data exist for the chemosensory behaviors of schistosomes, while the molecular sensory pathways in nematodes are well described. Recent investigations on the molecular mechanisms of sensation in schistosomes and filarial nematodes have revealed some features conserved within their respective phyla, but adaptations correlated with parasitism are pronounced. Technological developments are likely to extend these advances, and we forecast how these technologies may be applied.

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.

Ivermectin inhibits extracellular vesicle secretion from parasitic nematodes.

Lymphatic filariasis (LF) is a disease caused by parasitic filarial nematodes that is endemic in 49 countries of the world and affects or threatens over 890 million people. Strategies to control LF rely heavily on mass administration of anthelmintic drugs including ivermectin (IVM), a macrocyclic lactone drug considered an Essential Medicine by the WHO. However, despite its widespread use the therapeutic mode of action of IVM against filarial nematodes is not clear. We have previously reported that filarial nematodes secrete extracellular vesicles (EVs) and that their cargo has immunomodulatory properties. Here we investigate the effects of IVM and other anti-filarial drugs on parasitic nematode EV secretion, motility, and protein secretion. We show that inhibition of EV secretion was a specific property of IVM, which had consistent and significant inhibitory effects across nematode life stages and species, with the exception of male parasites. IVM inhibited EV secretion, but not parasite motility, at therapeutically relevant concentrations. Protein secretion was inhibited by IVM in the microfilariae stage, but not in any other stage tested. Our data provides evidence that inhibiting the secretion of immunomodulatory EVs by parasitic nematodes could explain, at least in part, IVM mode of action and provides a phenotype for novel drug discovery.

Long-read RNA sequencing of human and animal filarial parasites improves gene models and discovers operons.

Filarial parasitic nematodes (Filarioidea) cause substantial disease burden to humans and animals around the world. Recently there has been a coordinated global effort to generate, annotate, and curate genomic data from nematode species of medical and veterinary importance. This has resulted in two chromosome-level assemblies (Brugia malayi and Onchocerca volvulus) and 11 additional draft genomes from Filarioidea. These reference assemblies facilitate comparative genomics to explore basic helminth biology and prioritize new drug and vaccine targets. While the continual improvement of genome contiguity and completeness advances these goals, experimental functional annotation of genes is often hindered by poor gene models. Short-read RNA sequencing data and expressed sequence tags, in cooperation with ab initio prediction algorithms, are employed for gene prediction, but these can result in missing clade-specific genes, fragmented models, imperfect mapping of gene ends, and lack of isoform resolution. Long-read RNA sequencing can overcome these drawbacks and greatly improve gene model quality. Here, we present Iso-Seq data for B. malayi and Dirofilaria immitis, etiological agents of lymphatic filariasis and canine heartworm disease, respectively. These data cover approximately half of the known coding genomes and substantially improve gene models by extending untranslated regions, cataloging novel splice junctions from novel isoforms, and correcting mispredicted junctions. Furthermore, we validated computationally predicted operons, manually curated new operons, and merged fragmented gene models. We carried out analyses of poly(A) tails in both species, leading to the identification of non-canonical poly(A) signals. Finally, we prioritized and assessed known and putative anthelmintic targets, correcting or validating gene models for molecular cloning and target-based anthelmintic screening efforts. Overall, these data significantly improve the catalog of gene models for two important parasites, and they demonstrate how long-read RNA sequencing should be prioritized for ongoing improvement of parasitic nematode genome assemblies.

Genetic and functional diversification of chemosensory pathway receptors in mosquito-borne filarial nematodes.

Lymphatic filariasis (LF) afflicts over 60 million people worldwide and leads to severe pathological outcomes in chronic cases. The nematode parasites (Nematoda: Filarioidea) that cause LF require both arthropod (mosquito) intermediate hosts and mammalian definitive hosts for their propagation. The invasion and migration of filarial worms through host tissues are complex and critical to survival, yet little is known about the receptors and signaling pathways that mediate directed migration in these medically important species. In order to better understand the role of chemosensory signaling in filarial worm taxis, we employ comparative genomics, transcriptomics, reverse genetics, and chemical approaches to identify putative chemosensory receptor proteins and perturb chemotaxis phenotypes in filarial worms. We find that chemoreceptor family size is correlated with the presence of environmental (extrahost) stages in nematode life cycles, and that filarial worms contain compact and highly diverged chemoreceptor complements and lineage-specific ion channels that are predicted to operate downstream of chemoreceptor activation. In Brugia malayi, an etiological agent of LF, chemoreceptor expression patterns correspond to distinct parasite migration events across the life cycle. To interrogate the role of chemosensation in the migration of larval worms, arthropod and mammalian infectious stage Brugia parasites were incubated in nicotinamide, an agonist of the nematode transient receptor potential (TRP) channel OSM-9. Exposure of microfilariae to nicotinamide alters intramosquito migration, and exposure of L3s reduces chemotaxis toward host-associated cues in vitro. Nicotinamide also potently modulates thermosensory responses in L3s, suggesting a polymodal sensory role for Brugia osm-9. Reverse genetic studies implicate both Brugia osm-9 and the cyclic nucleotide-gated (CNG) channel subunit tax-4 in larval chemotaxis toward host serum, and these ion channel subunits partially rescue sensory defects in Caenorhabditis elegans osm-9 and tax-4 knock-out strains. Together, these data reveal genetic and functional diversification of chemosensory signaling proteins in filarial worms and encourage a more thorough investigation of clade- and parasite-specific facets of nematode sensory receptor biology.

Correction to: Sequence and structural variation in the genome of the Biomphalaria glabrata embryonic (Bge) cell line.

Following publication of the original article [], the authors reported an error in figure 1.

Sequence and structural variation in the genome of the Biomphalaria glabrata embryonic (Bge) cell line.

BACKGROUND: The aquatic pulmonate snail Biomphalaria glabrata is a significant vector and laboratory host for the parasitic flatworm Schistosoma mansoni, an etiological agent for the neglected tropical disease schistosomiasis. Much is known regarding the host-parasite interactions of these two organisms, and the B. glabrata embryonic (Bge) cell line has been an invaluable resource in these studies. The B. glabrata BB02 genome sequence was recently released, but nothing is known of the sequence variation between this reference and the Bge cell genome, which has likely accumulated substantial genetic variation in the ~50 years since its isolation. RESULTS: Here, we report the genome sequence of our laboratory subculture of the Bge cell line (designated Bge3), which we mapped to the B. glabrata BB02 reference genome. Single nucleotide variants (SNVs) were predicted and focus was given to those SNVs that are most likely to affect the structure or expression of protein-coding genes. Furthermore, we have highlighted and validated high-impact SNVs in genes that have often been studied using Bge cells as an in vitro model, and other genes that may have contributed to the immortalization of this cell line. We also resolved representative karyotypes for the Bge3 subculture, which revealed a mixed population exhibiting substantial aneuploidy, in line with previous reports from other Bge subcultures. CONCLUSIONS: The Bge3 genome differs from the B. glabrata BB02 reference genome in both sequence and structure, and these are likely to have significant biological effects. The availability of the Bge3 genome sequence, and an awareness of genomic differences with B. glabrata, will inform the design of experiments to understand gene function in this unique in vitro snail cell model. Additionally, this resource will aid in the development of new technologies and molecular approaches that promise to reveal more about this schistosomiasis-transmitting snail vector.

Functional analysis of Girardia tigrina transcriptome seeds pipeline for anthelmintic target discovery.

BACKGROUND: Neglected diseases caused by helminth infections impose a massive hindrance to progress in the developing world. While basic research on parasitic flatworms (platyhelminths) continues to expand, researchers have yet to broadly adopt a free-living model to complement the study of these important parasites. METHODS: We report the high-coverage sequencing (RNA-Seq) and assembly of the transcriptome of the planarian Girardia tigrina across a set of dynamic conditions. The assembly was annotated and extensive orthology analysis was used to seed a pipeline for the rational prioritization and validation of putative anthelmintic targets. A small number of targets conserved between parasitic and free-living flatworms were comparatively interrogated. RESULTS: 240 million paired-end reads were assembled de novo to produce a strictly filtered predicted proteome consisting of over 22,000 proteins. Gene Ontology annotations were extended to 16,467 proteins. 2,693 sequences were identified in orthology groups spanning flukes, tapeworms and planaria, with 441 highlighted as belonging to druggable protein families. Chemical inhibitors were used on three targets in pharmacological screens using both planaria and schistosomula, revealing distinct motility phenotypes that were shown to correlate with planarian RNAi phenotypes. CONCLUSIONS: This work provides the first comprehensive and annotated sequence resource for the model planarian G. tigrina, alongside a prioritized list of candidate drug targets conserved among parasitic and free-living flatworms. As proof of principle, we show that a simple RNAi and pharmacology pipeline in the more convenient planarian model system can inform parasite biology and serve as an efficient screening tool for the identification of lucrative anthelmintic targets.

Novel RNAi-mediated approach to G protein-coupled receptor deorphanization: proof of principle and characterization of a planarian 5-HT receptor.

G protein-coupled receptors (GPCRs) represent the largest known superfamily of membrane proteins extending throughout the Metazoa. There exists ample motivation to elucidate the functional properties of GPCRs given their role in signal transduction and their prominence as drug targets. In many target organisms, these efforts are hampered by the unreliable nature of heterologous receptor expression platforms. We validate and describe an alternative loss-of-function approach for ascertaining the ligand and G protein coupling properties of GPCRs in their native cell membrane environment. Our efforts are focused on the phylum Platyhelminthes, given the heavy health burden exacted by pathogenic flatworms, as well as the role of free-living flatworms as model organisms for the study of developmental biology. RNA interference (RNAi) was used in conjunction with a biochemical endpoint assay to monitor cAMP modulation in response to the translational suppression of individual receptors. As proof of principle, this approach was used to confirm the neuropeptide GYIRFamide as the cognate ligand for the planarian neuropeptide receptor GtNPR-1, while revealing its endogenous coupling to Gα(i/o). The method was then extended to deorphanize a novel Gα(s)-coupled planarian serotonin receptor, DtSER-1. A bioinformatics protocol guided the selection of receptor candidates mediating 5-HT-evoked responses. These results provide functional data on a neurotransmitter central to flatworm biology, while establishing the great potential of an RNAi-based deorphanization protocol. Future work can help optimize and adapt this protocol for higher-throughput platforms as well as other phyla.