Draft genome of the filarial nematode parasite Brugia malayi

Transcriptional responses to in vitro macrocyclic lactone exposure in Toxocara canis larvae using RNA-seq

Toxocara canis, the causative agent of zoonotic toxocariasis in humans, is a parasitic roundworm of canids with a complex lifecycle. While macrocyclic lactones (MLs) are successful at treating adult T. canis infections when used at FDA-approved doses in dogs, they fail to kill somatic third-stage larvae. In this study, we profiled the transcriptome of third-stage larvae derived from larvated eggs and treated in vitro with 10 μM of the MLs - ivermectin and moxidectin with Illumina sequencing. We analyzed transcriptional changes in comparison with untreated control larvae. In ivermectin-treated larvae, we identified 608 differentially expressed genes (DEGs), of which 453 were upregulated and 155 were downregulated. In moxidectin-treated larvae, we identified 1,413 DEGs, of which 902 were upregulated and 511 were downregulated. Notably, many DEGs were involved in critical biological processes and pathways including transcriptional regulation, energy metabolism, neuronal structure and function, physiological processes such as reproduction, excretory/secretory molecule production, host-parasite response mechanisms, and parasite elimination. We also assessed the expression of known ML targets and transporters, including glutamate-gated chloride channels (GluCls), and ATP-binding cassette (ABC) transporters, subfamily B, with a particular focus on P-glycoproteins (P-gps). We present gene names for previously uncharacterized T. canis GluCl genes using phylogenetic analysis of nematode orthologs to provide uniform gene nomenclature. Our study revealed that the expression of Tca-glc-3 and six ABCB genes, particularly four P-gps, were significantly altered in response to ML treatment. Compared to controls, Tca-glc-3, Tca-Pgp-11.2, and Tca-Pgp-13.2 were downregulated in ivermectin-treated larvae, while Tca-abcb1, Tca-abcb7, Tca-Pgp-11.2, and Tca-Pgp-13.2 were downregulated in moxidectin-treated larvae. Conversely, Tca-abcb9.1 and Tca-Pgp-11.3 were upregulated in moxidectin-treated larvae. These findings suggest that MLs broadly impact transcriptional regulation in T. canis larvae.

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.

Visualizing genomic evolution in Caenorhabditis through WormSynteny

Understanding the syntenic relationships among genomes is crucial to elucidate the genomic mechanisms that drive the evolution of species. The nematode Caenorhabditis is a good model for studying genomic evolution due to the well-established biology of Caenorhabditis elegans and the availability of > 50 genomes in the genus. However, effective alignment of more than ten species in Caenorhabditis has not been conducted before, and there is currently no tool to visualize the synteny of more than two species. In this study, we used Progressive Cactus, a recently developed multigenome aligner, to align the genomes of eleven Caenorhabditis species. Through the progressive alignment, we reconstructed nine ancestral genomes, analyzed the mutational types that cause genomic rearrangement during speciation, and found that insertion and duplication are the major driving forces for genome expansion. Dioecious species appear to expand their genomes more than androdioecious species. We then built an online interactive app called WormSynteny to visualize the syntenic relationship among the eleven species. Users can search the alignment dataset using C. elegans query sequences, construct synteny plots at different genomic scales, and use a set of options to control alignment output and plot presentation. We showcased the use of WormSynteny to visualize the syntenic conservation of one-to-one orthologues among species, tandem and dispersed gene duplication in C. elegans, and the evolution of exon and intron structures. Importantly, the integration of orthogroup information with synteny linkage in WormSynteny allows the easy visualization of conserved genomic blocks and disruptive rearrangement. In conclusion, WormSynteny provides immediate access to the syntenic relationships among the most widely used Caenorhabditis species and can facilitate numerous comparative genomics studies. This pilot study with eleven species also serves as a proof-of-concept to a more comprehensive larger-scale analysis using hundreds of nematode genomes, which is expected to reveal mechanisms that drive genomic evolution in the Nematoda phylum. Finally, the WormSynteny software provides a generalizable solution for visualizing the output of Progressive Cactus with interactive graphics, which would be useful for a broad community of genome researchers.

Global analysis of neuropeptide receptor conservation across phylum Nematoda

BACKGROUND

The phylum Nematoda is incredibly diverse and includes many parasites of humans, livestock, and plants. Peptide-activated G protein-coupled receptors (GPCRs) are central to the regulation of physiology and numerous behaviors, and they represent appealing pharmacological targets for parasite control. Efforts are ongoing to characterize the functions and define the ligands of nematode GPCRs, with already most peptide GPCRs known or predicted in Caenorhabditis elegans. However, comparative analyses of peptide GPCR conservation between C. elegans and other nematode species are limited, and many nematode GPCRs remain orphan. A phylum-wide perspective on peptide GPCR profiles will benefit functional and applied studies of nematode peptide GPCRs.

RESULTS

We constructed a pan-phylum resource of C. elegans peptide GPCR orthologs in 125 nematode species using a semi-automated pipeline for analysis of predicted proteome datasets. The peptide GPCR profile varies between nematode species of different phylogenetic clades and multiple C. elegans peptide GPCRs have orthologs across the phylum Nematoda. We identified peptide ligands for two highly conserved orphan receptors, NPR-9 and NPR-16, that belong to the bilaterian galanin/allatostatin A (Gal/AstA) and somatostatin/allatostatin C (SST/AstC) receptor families. The AstA-like NLP-1 peptides activate NPR-9 in cultured cells and are cognate ligands of this receptor in vivo. In addition, we discovered an AstC-type peptide, NLP-99, that activates the AstC-type receptor NPR-16. In our pan-phylum resource, the phylum-wide representation of NPR-9 and NPR-16 resembles that of their cognate ligands more than those of allatostatin-like peptides that do not activate these receptors.

CONCLUSIONS

The repertoire of C. elegans peptide GPCR orthologs varies across phylogenetic clades and several peptide GPCRs show broad conservation in the phylum Nematoda. Our work functionally characterizes the conserved receptors NPR-9 and NPR-16 as the respective GPCRs for the AstA-like NLP-1 peptides and the AstC-related peptide NLP-99. NLP-1 and NLP-99 are widely conserved in nematodes and their representation matches that of their receptor in most species. These findings demonstrate the conservation of a functional Gal/AstA and SST/AstC signaling system in nematodes. Our dataset of C. elegans peptide GPCR orthologs also lays a foundation for further functional studies of peptide GPCRs in the widely diverse nematode phylum.

Genome-wide RNA interference of the nhr gene family in barber's pole worm identified members crucial for larval viability in vitro

Nuclear hormone receptors (NHRs) are emerging target candidates against nematode infection and resistance. However, there is a lack of comprehensive information on NHR-coding genes in parasitic nematodes. In this study, we curated the nhr gene family for 60 major parasitic nematodes from humans and animals. Compared with the free-living model organism Caenorhabditis elegans, a remarkable contraction of the nhr family was revealed in parasitic species, with genetic diversification and conservation unveiled among nematode Clades I (10-13), III (16-42), IV (33-35) and V (25-64). Using an in vitro biosystem, we demonstrated that 40 nhr genes in a blood-feeding nematode Haemonchus contortus (clade V; barber's pole worm) were responsive to host serum and one nhr gene (i.e., nhr-64) was consistently stimulated by anthelmintics (i.e., ivermectin, thiabendazole and levamisole); Using a high-throughput RNA interference platform, we knocked down 43 nhr genes of H. contortus and identified at least two genes that are required for the viability (i.e., nhr-105) and development (i.e., nhr-17) of the infective larvae of this parasitic nematode in vitro. Harnessing this preliminary functional atlas of nhr genes for H. contortus will prime the biological studies of this gene family in nematode genetics, infection, and anthelmintic metabolism within host animals, as well as the promising discovery of novel intervention targets.

Distinguishing recrudescence from reinfection in lymphatic filariasis

BACKGROUND

The Global Program to Eliminate Lymphatic Filariasis (GPELF) is the largest public health program based on mass drug administration (MDA). Despite decades of MDA, ongoing transmission in some countries remains a challenge. To optimise interventions, it is critical to differentiate between recrudescence and new infections. Since adult filariae are inaccessible in humans, deriving a method that relies on the offspring microfilariae (mf) is necessary.

METHODS

We developed a genome amplification and kinship analysis-based approach using Brugia malayi samples from gerbils, and applied it to analyse Wuchereria bancrofti mf from humans in Côte d'Ivoire. We examined the pre-treatment genetic diversity in 269 mf collected from 18 participants, and further analysed 1-year post-treatment samples of 74 mf from 4 participants. Hemizygosity of the male X-chromosome allowed for direct inference of haplotypes, facilitating robust maternal parentage inference. To enrich parasite DNA from samples contaminated with host DNA, a whole-exome capture panel was created for W. bancrofti.

FINDINGS

By reconstructing and temporally tracking sibling relationships across pre- and post-treatment samples, we differentiated between new and established maternal families, suggesting reinfection in one participant and recrudescence in three participants. The estimated number of reproductively active adult females ranged between 3 and 11 in the studied participants. Population structure analysis revealed genetically distinct parasites in Côte d'Ivoire compared to samples from other countries. Exome capture identified protein-coding variants with ∼95% genotype concordance rate.

INTERPRETATION

We have generated resources to facilitate the development of molecular genetic tools that can estimate adult worm burdens and monitor parasite populations, thus providing essential information for the successful implementation of GPELF.

FUNDING

This work was financially supported by the Bill and Melinda Gates Foundation (https://www.gatesfoundation.org) under grant OPP1201530 (Co-PIs PUF & Gary J. Weil). B. malayi parasite material was generated with support of the Foundation for Barnes Jewish Hospital (PUF). In addition, the development of computational methods was supported by the National Institutes of Health under grants AI144161 (MM) and AI146353 (MM). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Common analysis of direct RNA sequencinG CUrrently leads to misidentification of m5C at GCU motifs

RNA modifications, such as methylation, can be detected with Oxford Nanopore Technologies direct RNA sequencing. One commonly used tool for detecting 5-methylcytosine (m5C) modifications is Tombo, which uses an "Alternative Model" to detect putative modifications from a single sample. We examined direct RNA sequencing data from diverse taxa including viruses, bacteria, fungi, and animals. The algorithm consistently identified a m5C at the central position of a GCU motif. However, it also identified a m5C in the same motif in fully unmodified in vitro transcribed RNA, suggesting that this is a frequent false prediction. In the absence of further validation, several published predictions of m5C in a GCU context should be reconsidered, including those from human coronavirus and human cerebral organoid samples.

Understanding anthelmintic resistance in livestock using "omics" approaches

Widespread and improper use of various anthelmintics, genetic, and epidemiological factors has resulted in anthelmintic-resistant (AR) helminth populations in livestock. This is currently quite common globally in different livestock animals including sheep, goats, and cattle to gastrointestinal nematode (GIN) infections. Therefore, the mechanisms underlying AR in parasitic worm species have been the subject of ample research to tackle this challenge. Current and emerging technologies in the disciplines of genomics, transcriptomics, metabolomics, and proteomics in livestock species have advanced the understanding of the intricate molecular AR mechanisms in many major parasites. The technologies have improved the identification of possible biomarkers of resistant parasites, the ability to find actual causative genes, regulatory networks, and pathways of parasites governing the AR development including the dynamics of helminth infection and host-parasite infections. In this review, various "omics"-driven technologies including genome scan, candidate gene, quantitative trait loci, transcriptomic, proteomic, and metabolomic approaches have been described to understand AR of parasites of veterinary importance. Also, challenges and future prospects of these "omics" approaches are also discussed.

Present status with impacts and roles of miRNA on Soil Transmitted Helminthiosis control: A review

Soil-Transmitted Helminthiasis (STH) is one of the most widespread Neglected Tropical Diseases (NTDs), and almost 1.5 billion of the global population is affected, mostly in the indigent, countryside sectors of tropics/subtropics. STH, commonly caused by various nematodes, adversely affects the hosts' growth, cognatic development, and immunity. Albendazole is most commonly used against STH (Soil-Transmitted Helminths) but resistance has already been reported in different countries. To date, no effective vaccine is present against STH. miRNAs are a unique class of small non-coding RNA, regulating various biological activities indulging host immune responses in host-pathogen interaction of STH. Dysregulation of miRNAs are being considered as one of the most important aspect of host-parasite interactions. Thus, it is the prime importance to identify and characterize parasite-specific as well as host-derived miRNAs to understand the STH infection at the molecular level. Systematic bibliometric analysis reveals a huge knowledge gap in understanding the disease by using both host and parasitic miRNAs as a potential biomarker. In this study, we addressed the present status of the STH prevalence, and therapy under the light of miRNAs. This would further help in designing new inhibitors and therapeutic strategies to control STH.

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.

NHR-23 activity is necessary for C. elegans developmental progression and apical extracellular matrix structure and function

Nematode molting is a remarkable process where animals must repeatedly build a new apical extracellular matrix (aECM) beneath a previously built aECM that is subsequently shed. The nuclear hormone receptor NHR-23 (also known as NR1F1) is an important regulator of C. elegans molting. NHR-23 expression oscillates in the epidermal epithelium, and soma-specific NHR-23 depletion causes severe developmental delay and death. Tissue-specific RNAi suggests that nhr-23 acts primarily in seam and hypodermal cells. NHR-23 coordinates the expression of factors involved in molting, lipid transport/metabolism and remodeling of the aECM. NHR-23 depletion causes dampened expression of a nas-37 promoter reporter and a loss of reporter oscillation. The cuticle collagen ROL-6 and zona pellucida protein NOAH-1 display aberrant annular localization and severe disorganization over the seam cells after NHR-23 depletion, while the expression of the adult-specific cuticle collagen BLI-1 is diminished and frequently found in patches. Consistent with these localization defects, the cuticle barrier is severely compromised when NHR-23 is depleted. Together, this work provides insight into how NHR-23 acts in the seam and hypodermal cells to coordinate aECM regeneration during development.

Common Analysis of Direct RNA SequencinG CUrrently Leads to Misidentification of 5-Methylcytosine Modifications at GCU Motifs

RNA modifications, such as méthylation, can be detected with Oxford Nanopore Technologies direct RNA sequencing. One commonly used tool for detecting 5-methylcytosine (m5C) modifications is Tombo, which uses an "Alternative Model" to detect putative modifications from a single sample. We examined direct RNA sequencing data from diverse taxa including virus, bacteria, fungi, and animals. The algorithm consistently identified a 5-methylcytosine at the central position of a GCU motif. However, it also identified a 5-methylcytosine in the same motif in fully unmodified in vitro transcribed RNA, suggesting that this a frequent false prediction. In the absence of further validation, several published predictions of 5-methylcytosine in human coronavirus and human cerebral organoid RNA in a GCU context should be reconsidered.

Wolbachia Ferrochelatase as a potential drug target against filarial infections

Filarial infections are among the world's most disturbing diseases caused by 3 major parasitic worms; Onchocerca volvulus, Wuchereria bancrofti, and Brugia malayi, affecting more than 500 million people worldwide. Currently used drugs for mass drug administration (MDA) have been met with several challenges including the development of complications in individuals with filaria co-infections and parasitic drug resistance. The filarial endosymbiont, Wolbachia, has emerged as an attractive therapeutic target for filariasis elimination, due to the dependence of the filaria on this endosymbiont for survival. Here, we target an important enzyme in the Wolbachia heme biosynthetic pathway (ferrochelatase), using high-throughput virtual screening and molecular dynamics with MM-PBSA calculations. We identified four drug candidates; Nilotinib, Ledipasvir, 3-benzhydryloxy-8-methyl-8-azabicyclo[3.2.1]octane, and 2-(4-Amino-piperidin-1-yl)-ethanol as potential small molecules inhibitors as they could compete with the enzyme's natural substrate (Protoporphyrin IX) for active pocket binding. This prevents the worm from receiving the heme molecule from Wolbachia for their growth and survival, resulting in their death. This study which involved targeting enzymes in biosynthetic pathways of the parasitic worms' endosymbiont (Wolbachia), has proven to be an alternative therapeutic option leading to the discovery of new drugs, which will help facilitate the elimination of parasitic infections.

Genomes of the human filarial parasites Mansonella perstans and Mansonella ozzardi

The filarial parasites Mansonella ozzardi and Mansonella perstans, causative agents of mansonellosis, infect hundreds of millions of people worldwide, yet remain among the most understudied of the human filarial pathogens. M. ozzardi is highly prevalent in Latin American countries and Caribbean Islands, while M. perstans is predominantly found in sub-Saharan Africa as well as in a few areas in South America. In addition to the differences in their geographical distribution, the two parasites are transmitted by different insect vectors, as well as exhibit differences in their responses to commonly used anthelminthic drugs. The lack of genome information has hindered investigations into the biology and evolution of Mansonella parasites and understanding the molecular basis of the clinical differences between species. In the current study, high quality genomes of two independent clinical isolates of M. perstans from Cameroon and two M. ozzardi isolates one from Brazil and one from Venezuela are reported. The genomes are approximately 76 Mb in size, encode about 10,000 genes each, and are largely complete based on BUSCO scores of about 90%, similar to other completed filarial genomes. These sequences represent the first genomes from Mansonella parasites and enabled a comparative genomic analysis of the similarities and differences between Mansonella and other filarial parasites. Horizontal DNA transfers (HDT) from mitochondria (nuMTs) as well as transfers from genomes of endosymbiotic Wolbachia bacteria (nuWTs) to the host nuclear genome were identified and analyzed. Sequence comparisons and phylogenetic analysis of known targets of anti-filarial drugs diethylcarbamazine (DEC), ivermectin and mebendazole revealed that all known target genes were present in both species, except for the DEC target encoded by gon-2 gene, which is fragmented in genome assemblies from both M. ozzardi isolates. These new reference genome sequences will provide a valuable resource for further studies on biology, symbiosis, evolution and drug discovery.

Wolbachia Effect on Drosophila melanogaster Lipid and Carbohydrate Metabolism

The effect of maternally inherited endosymbiotic bacteria Wolbachia on triglyceride and carbohydrate metabolism, starvation resistance and feeding behavior of Drosophila melanogaster females was studied. Eight D. melanogaster lines of the same nuclear background were investigated; one had no infection and served as the control, and seven others were infected with different Wolbachia strains pertaining to wMel and wMelCS groups of genotypes. Most of the infected lines had a higher overall lipid content and triglyceride level than the control line and their expression of the bmm gene regulating triglyceride catabolism was reduced. The glucose content was higher in the infected lines compared to that in the control, while their trehalose levels were similar. It was also found that the Wolbachia infection reduced the level of tps1 gene expression (coding for enzyme for trehalose synthesis from glucose) and had no effect on treh gene expression (coding for trehalose degradation enzyme). The infected lines exhibited lower appetite but higher survival under starvation compared to the control. The data obtained may indicate that Wolbachia foster their hosts' energy exchange through increasing its lipid storage and glucose content to ensure the host's competitive advantage over uninfected individuals. The scheme of carbohydrate and lipid metabolism regulation under Wolbachia's influence was suggested.

A proteasomal β5 subunit of Haemonchus contortus with a role in the growth, development and life span

BACKGROUND

The proteasome in eukaryotic cells can degrade a variety of proteins and plays an important role in regulating the cell cycle, cell survival and apoptosis. The proteasome receives much attention as a potential chemotherapeutic target for treatment of a variety of infectious parasitic diseases, but few studies of proteasomes have been done on parasitic nematodes.

METHODS

A proteasomal β5 subunit encoding gene (named Hc-pbs-5) and its inferred product (Hc-PBS-5) in Haemonchus contortus were identified and characterized in this study. Then, the transcriptional profiles and anatomical expression were studied using an integrated molecular approach. Finally, a specific proteasome inhibitor bortezomib (BTZ), together with RNA interference (RNAi), was employed to assess the function of Hc-PBS-5.

RESULTS

Bioinformatic analysis revealed that the coding sequence of Hc-pbs-5 was 855 bp long and encoded 284 amino acids (aa). The predicted protein (Hc-PBS-5) had core conservative sequences (65-250 aa) belonging to N-terminal nucleophile (Ntn) family of hydrolases. Real-time PCR results revealed that Hc-pbs-5 was continuously transcribed in eight developmental stages with higher levels at the infective third-stage larvae (L3s) and adult males of H. contortus. Immunohistochemical results revealed that Hc-PBS-5 was expressed in intestine, outer cuticle, muscle cells under the outer cuticle, cervical glands and seminal vesicles of male adults and also in intestine, outer cuticle, cervical glands, uterine wall, eggs and ovaries of female adults of H. contortus. BTZ could reduce proportions of egg hatching, and the fourth-stage larvae (L4s) developed from the exsheathed L3s (xL3s) of H. contortus. In addition, silencing Hc-pbs-5 by soaking the specific double-stranded RNA (dsRNA) could decrease the transcription of Hc-pbs-5 and result in fewer xL3s developing to L4s in vitro.

CONCLUSIONS

These results indicate that proteasomal β5 subunit plays an important role in the growth, development and life span of H. contortus.

Beyond bacteria: Reconstructing microorganism connections and deciphering the predicted mutualisms in mammalian gut metagenomes

Numerous gut microbial studies have focused on bacteria. However, archaea, viruses, fungi, protists, and nematodes are also regular residents of the gut ecosystem. Little is known about the composition and potential interactions among these six kingdoms in the same samples. Here, we unraveled the complex connection among them using approximately 123 gut metagenomes from 42 mammalian species (including carnivores, omnivores, and herbivores). We observed high variation in bacterial and fungal families and relatively low variation in archaea, viruses, protists, and nematodes. We found that some fungi in the mammalian intestine might come from environmental sources (e.g., soil and dietary plants), and some might be native to the intestine (e.g., the occurrence of Neocallimastigomycetes). The Methanobacteriaceae and Plasmodiidae families (archaea and protozoa, respectively) were predominant in these metagenomes, whereas Onchocercidae and Trichuridae were the two most common nematodes, and Siphoviridae and Myoviridae the two most common virus families in these mammalian gut metagenomes. Interestingly, most of the pairwise co-occurrence patterns were significantly positive among these six kingdoms, and significantly negative networks mainly occurred between fungi and prokaryotes (both bacteria and archaea). Our study revealed some inconvenient characteristics in the mammalian gut microorganism ecosystem: (1) the community formed by members of the analyzed kingdoms reflects the life history of the host and the potential threat posed by pathogenic protists and nematodes in mammals; and (2) the networks suggest the existence of predicted mutualism among members of these six kingdoms and of the predicted competition, mainly among fungi and other kingdoms.

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.

Transcriptomics of ivermectin response in Caenorhabditis elegans: Integrating abamectin quantitative trait loci and comparison to the Ivermectin-exposed DA1316 strain

Parasitic nematodes pose a significant threat to human and animal health, as well as cause economic losses in the agricultural sector. The use of anthelmintic drugs, such as Ivermectin (IVM), to control these parasites has led to widespread drug resistance. Identifying genetic markers of resistance in parasitic nematodes can be challenging, but the free-living nematode Caenorhabditis elegans provides a suitable model. In this study, we aimed to analyze the transcriptomes of adult C. elegans worms of the N2 strain exposed to the anthelmintic drug Ivermectin (IVM), and compare them to those of the resistant strain DA1316 and the recently identified Abamectin Quantitative Trait Loci (QTL) on chromosome V. We exposed pools of 300 adult N2 worms to IVM (10-7 and 10-8 M) for 4 hours at 20°C, extracted total RNA and sequenced it on the Illumina NovaSeq6000 platform. Differentially expressed genes (DEGs) were determined using an in-house pipeline. The DEGs were compared to genes from a previous microarray study on IVM-resistant C. elegans and Abamectin-QTL. Our results revealed 615 DEGs (183 up-regulated and 432 down-regulated genes) from diverse gene families in the N2 C. elegans strain. Of these DEGs, 31 overlapped with genes from IVM-exposed adult worms of the DA1316 strain. We identified 19 genes, including the folate transporter (folt-2) and the transmembrane transporter (T22F3.11), which exhibited an opposite expression in N2 and the DA1316 strain and were deemed potential candidates. Additionally, we compiled a list of potential candidates for further research including T-type calcium channel (cca-1), potassium chloride cotransporter (kcc-2), as well as other genes such as glutamate-gated channel (glc-1) that mapped to the Abamectin-QTL.

Reconstitution of an N-AChR from Brugia malayi an evolved change in acetylcholine receptor accessory protein requirements in filarial parasites

Neurotransmission is an important target for anthelmintic drugs, where receptor characteristics and response can be examined through reconstitution ex vivo in Xenopus laevis oocytes. The homomeric ACR-16 nicotine sensitive acetylcholine receptors (N-AChRs) of several helminth species have been characterized in this way. Our efforts to reconstitute the N-AChR from the clade III filarial parasite, Brugia malayi using similar conditions, initially produced no detectable response. A robust response to acetylcholine is obtained from the closely related clade III parasite Ascaris suum, suggesting that specific changes have occurred between Ascaris and Brugia. N-AChRs from three species intermediate between A. suum and B. malayi were characterized to provide information on the cause. Maximal response to acetylcholine did not change abruptly, consistent with a discrete event, but rather decreased progressively from A. suum through Dracunculus medinensis, Gonglylonema pulchrum and Thelazia callipaeda. Receptor responses to the characteristic nicotine, and other agonists were generally similar. The decrease in maximal current did correlate with a delayed time to reach larger response. Together, this suggested that the failure to reconstitute the B. malayi N-AChR was one extreme of a progressive decrease and that an issue with synthesis of the receptor in oocytes was responsible. Addition of accessory proteins EMC-6, NRA-2 and NRA-4, in addition to RIC-3, produced a small, but measurable B. malayi N-AChR response. Pharmacological properties of a chimeric B. malayi N-AChR were equivalent to the other species, confirming the receptor response remains unchanged while its production is increasingly dependent on accessory proteins. One possibility is that loss of many subunits for acetylcholine receptors from the filarial nematode genome is linked to new subunit combinations that lead to such a dependence. This novel phylogenetic approach allowed the first characterization of a B. malayi AChR ex vivo and in doing so, provides a framework for the successful characterization of other receptors that have yet to be reconstituted.

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

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

Sirt5 Inhibits BmNPV Replication by Promoting a Relish-Mediated Antiviral Pathway in Bombyx mori

Silent information regulators (Sirtuins) belong to the family of nicotinamide adenine dinucleotide (NAD⁺)-dependent histone deacetylases (HDACs) that have diverse functions in cells. Mammalian Sirtuins have seven isoforms (Sirt1–7) which have been found to play a role in viral replication. However, Sirtuin members of insects are very different from mammals, and the function of insect Sirtuins in regulating virus replication is unclear. The silkworm, Bombyx mori, as a model species of Lepidoptera, is also an important economical insect. B. mori nucleopolyhedrovirus (BmNPV) is a major pathogen that specifically infects silkworms and causes serious losses in the sericulture industry. Here, we used the infection of the silkworm by BmNPV as a model to explore the effect of Sirtuins on virus replication. We initially knocked down all silkworm Sirtuins, and then infected with BmNPV to analyze its replication. Sirt2 and Sirt5 were found to have potential antiviral functions in the silkworm. We further confirmed the antiviral function of silkworm Sirt5 through its effects on viral titers during both knockdown and overexpression experiments. Additionally, Suramin, a Sirt5 inhibitor, was found to promote BmNPV replication. In terms of molecular mechanism, it was found that silkworm Sirt5 might promote the immune pathway mediated by Relish, thereby enhancing the host antiviral response. This study is the first to explore the role of Sirtuins in insect-virus interactions, providing new insights into the functional role of members of the insect Sirtuin family.

Advances in our understanding of nematode ion channels as potential anthelmintic targets

Ion channels are specialized multimeric proteins that underlie cell excitability. These channels integrate with a variety of neuromuscular and biological functions. In nematodes, the physiological behaviors including locomotion, navigation, feeding and reproduction, are regulated by these protein entities. Majority of the antinematodal chemotherapeutics target the ion channels to disrupt essential biological functions. Here, we have summarized current advances in our understanding of nematode ion channel pharmacology. We review cys-loop ligand gated ion channels (LGICs), including nicotinic acetylcholine receptors (nAChRs), acetylcholine-chloride gated ion channels (ACCs), glutamate-gated chloride channels (GluCls), and GABA (γ-aminobutyric acid) receptors, and other ionotropic receptors (transient receptor potential (TRP) channels and potassium ion channels). We have provided an update on the pharmacological properties of these channels from various nematodes. This article catalogs the differences in ion channel composition and resulting pharmacology in the phylum Nematoda. This diversity in ion channel subunit repertoire and pharmacology emphasizes the importance of pursuing species-specific drug target research. In this review, we have provided an overview of recent advances in techniques and functional assays available for screening ion channel properties and their application.

Current perspective of new anti-Wolbachial and direct-acting macrofilaricidal drugs as treatment strategies for human filariasis

Filarial diseases like lymphatic filariasis and onchocerciasis belong to the Neglected Tropical Diseases and remain a public health problem in endemic countries. Lymphatic filariasis and onchocerciasis can lead to stigmatizing pathologies and present a socio-economic burden for affected people and their endemic countries. Current treatment recommendations by the WHO include mass drug administration with ivermectin for the treatment of onchocerciasis and a combination of ivermectin, albendazole and diethylcarbamazine (DEC) for the treatment of lymphatic filariasis in areas that are not co-endemic for onchocerciasis or loiasis. Limitations of these treatment strategies are due to potential severe adverse events in onchocerciasis and loiasis patients following DEC or ivermectin treatment, respectively, the lack of a macrofilaricidal efficacy of those drugs and the risk of drug resistance development. Thus, to achieve the elimination of transmission of onchocerciasis and the elimination of lymphatic filariasis as a public health problem by 2030, the WHO defined in its roadmap that new alternative treatment strategies with macrofilaricidal compounds are required. Within a collaboration of the non-profit organizations Drugs for Neglected Diseases initiative (DNDi), the Bill & Melinda Gates Foundation, and partners from academia and industry, several new promising macrofilaricidal drug candidates were identified, which will be discussed in this review.

Transgenesis in parasitic helminths: a brief history and prospects for the future

Helminth infections impact the health of hundreds of millions of persons globally and also cause important economic losses in livestock farming. Methodological limitations as well as the low attention given to the study of helminths have impacted biological research and, thus, the procurement of accurate diagnosis and effective treatments. Understanding the biology of helminths using genomic and proteomic approaches could contribute to advances in understanding host-helminth interactions and lead to new vaccines, drugs and diagnostics. Despite the significant advances in genomics in the last decade, the lack of methodological adaptation of current transgenesis techniques has hampered the progression of post-genomic research in helminthology. However, the application of new techniques, such as CRISPR, to the study of trematodes and nematodes has opened new avenues for genome editing-powered functional genomics for these pathogens. This review summarises the historical advances in functional genomics in parasitic helminths and highlights pending limitations that will need to be overcome to deploy transgenesis tools.

From Extrapolation to Precision Chemical Hazard Assessment: The Ecdysone Receptor Case Study

Hazard assessment strategies are often supported by extrapolation of damage probabilities, regarding chemical action and species susceptibilities. Yet, growing evidence suggests that an adequate sampling of physiological responses across a representative taxonomic scope is of paramount importance. This is particularly relevant for Nuclear Receptors (NR), a family of transcription factors, often triggered by ligands and thus, commonly exploited by environmental chemicals. Within NRs, the ligand-induced Ecdysone Receptor (EcR) provides a remarkable example. Long regarded as arthropod specific, this receptor has been extensively targeted by pesticides, seemingly innocuous to non-target organisms. Yet, current evidence clearly suggests a wider presence of EcR orthologues across metazoan lineages, with unknown physiological consequences. Here, we address the state-of-the-art regarding the phylogenetic distribution and functional characterization of metazoan EcRs and provide a critical analysis of the potential disruption of such EcRs by environmental chemical exposure. Using EcR as a case study, hazard assessment strategies are also discussed in view of the development of a novel "precision hazard assessment paradigm.

Targeted transcriptomics reveals signatures of large-scale independent origins and concerted regulation of effector genes in Radopholus similis

The burrowing nematode, Radopholus similis, is an economically important plant-parasitic nematode that inflicts damage and yield loss to a wide range of crops. This migratory endoparasite is widely distributed in warmer regions and causes extensive destruction to the root systems of important food crops (e.g., citrus, banana). Despite the economic importance of this nematode, little is known about the repertoire of effectors owned by this species. Here we combined spatially and temporally resolved next-generation sequencing datasets of R. similis to select a list of candidates for the identification of effector genes for this species. We confirmed spatial expression of transcripts of 30 new candidate effectors within the esophageal glands of R. similis by in situ hybridization, revealing a large number of pioneer genes specific to this nematode. We identify a gland promoter motif specifically associated with the subventral glands (named Rs-SUG box), a putative hallmark of spatial and concerted regulation of these effectors. Nematode transcriptome analyses confirmed the expression of these effectors during the interaction with the host, with a large number of pioneer genes being especially abundant. Our data revealed that R. similis holds a diverse and emergent repertoire of effectors, which has been shaped by various evolutionary events, including neofunctionalization, horizontal gene transfer, and possibly by de novo gene birth. In addition, we also report the first GH62 gene so far discovered for any metazoan and putatively acquired by lateral gene transfer from a bacterial donor. Considering the economic damage caused by R. similis, this information provides valuable data to elucidate the mode of parasitism of this nematode.

Advances in the treatment, diagnosis, control and scientific understanding of taeniid cestode parasite infections over the past 50 years

In the past 50 years, enormous progress has been made in the diagnosis, treatment and control of taeniid cestode infections/diseases and in the scientific understanding thereof. Most interest in this group of parasites stems from the serious diseases that they cause in humans. It is through this lens that we summarize here the most important breakthroughs that have made a difference to the treatment of human diseases caused by these parasites, reduction in transmission of the taeniid species associated with human disease, or understanding of the parasites' biology likely to impact diagnosis or treatment in the foreseeable future. Key topics discussed are the introduction of anti-cestode drugs, including benzimidazoles and praziquantel, and the development of new imaging modalities that have transformed the diagnosis and post-treatment monitoring of human echinococcoses and neurocysticercosis. The availability of new anti-cestode drugs for use in dogs and a detailed understanding of the transmission dynamics of Echinococcus granulosus sensu lato have underpinned successful programs that have eliminated cystic echinococcosis in some areas of the world and greatly reduced the incidence of infection in others. Despite these successes, cystic and alveolar echinococcosis and neurocysticercosis continue to be prevalent in many parts of the world, requiring new or renewed efforts to prevent the associated taeniid infections. Major advances made in the development of practical vaccines against E. granulosus and Taenia solium will hopefully assist in this endeavour, as might the understanding of the parasites' biology that have come from an elucidation of the nuclear genomes of each of the most important taeniid species causing human diseases.

Comparative mitogenomics of Spirocerca lupi from South Africa and China: Variation and possible heteroplasmy

The complete mitochondrial genome of Spirocerca lupi isolated from a dog in South Africa was sequenced using next generation sequencing (NGS) technology and the 12 protein coding genes along with the two rRNA genes were compared to 18 other nematode species as well as S. lupi from China. The mitochondrial genome of S. lupi South Africa had a mean genetic diversity of 6.1 % compared to S. lupi China with some variation in nucleotide composition, gene positioning and size. Pairwise distance results indicated slightly higher variation when compared to the pairwise distances of other closely related species, however, this variation was not high enough for it to be considered a cryptic species. Phylogenetic analysis indicated that S. lupi from the two continents are very similar. In addition, single nucleotide polymorphisms were detected in the nad2 gene with ten sequence variants identified from 10 clones from a single nematode, suggesting possible heteroplasmy. The origin of the heteroplasmy is currently unknown but it is speculated to have arisen from accumulated mutations in the mitochondria during somatic replication.

From Magic Bullet to Magic Bomb: Reductive Bioactivation of Antiparasitic Agents

Paul Ehrlich coined the term "magic bullet" to describe how a drug kills the parasite inside its human host without harming the host itself. Ehrlich concluded that the drug must have a greater affinity to the parasite than to human cells. Today, the specificity of drug action is understood in terms of the drug target. An ideal target is a protein that is essential for the proliferation of the pathogen but absent in human cells. Examples are the enzymes of folate synthesis or of the nonmevalonate pathway in the malaria parasites. However, there are other ways how a drug can kill selectively. Of particular relevance is the specific activation of a prodrug inside the pathogen but not in the host, as this is how the current frontrunners of parasite chemotherapy work. Artemisinins for malaria, fexinidazole for human African trypanosomiasis, benznidazole for Chagas' disease, metronidazole for intestinal protozoa: these molecules are "magic bombs" that are triggered selectively. They are prodrugs that need to be activated by chemical reduction, i.e., the acquisition of an electron, which occurs in the parasite. Such a mode of action is shared by the novel antimalarial peroxides arterolane and artefenomel, which are activated by reduction of the endoperoxide bond with ferrous heme as the likely electron donor, a metabolic end-product of Plasmodium falciparum. Here we provide an overview on the molecular basis of selectivity of antiparasitic drug action with particular reference to the ozonides, the new generation of antimalarial peroxides designed by Jonathan Vennerstrom.

X-treme loss of sequence diversity linked to neo-X chromosomes in filarial nematodes

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication.

Chromosome-scale assembly and whole-genome sequencing of 266 giant panda roundworms provide insights into their evolution, adaptation and potential drug targets

Helminth diseases have long been a threat to the health of humans and animals. Roundworms are important organisms for studying parasitic mechanisms, disease transmission and prevention. The study of parasites in the giant panda is of importance for understanding how roundworms adapt to the host. Here, we report a high‐quality chromosome‐scale genome of Baylisascaris schroederi with a genome size of 253.60 Mb and 19,262 predicted protein‐coding genes. We found that gene families related to epidermal chitin synthesis and environmental information processes in the roundworm genome have expanded significantly. Furthermore, we demonstrated unique genes involved in essential amino acid metabolism in the B. schroederi genome, inferred to be essential for the adaptation to the giant panda‐specific diet. In addition, under different deworming pressures, we found that four resistance‐related genes (glc‐1, nrf‐6, bre‐4 and ced‐7) were under strong positive selection in a captive population. Finally, 23 known drug targets and 47 potential drug target proteins were identified. The genome provides a unique reference for inferring the early evolution of roundworms and their adaptation to the host. Population genetic analysis and drug sensitivity prediction provide insights revealing the impact of deworming history on population genetic structure of importance for disease prevention.

Brugia malayi Glycoproteins Detected by the Filariasis Test Strip Antibody AD12.1

Background

Rapid and accurate prevalence mapping of lymphatic filariasis (LF) is necessary to eliminate this disfiguring and disabling neglected tropical disease. Unfortunately, rapid tests such as the filariasis test strip (FTS) for Wuchereria bancrofti, the causative agent of LF in Africa, can cross-react with antigens circulating in some persons infected by the African eye worm, Loa loa, rendering the test unreliable in eleven co-endemic nations. The intended target of the FTS is a heavily glycosylated W. bancrofti circulating filarial antigen (Wb-CFA). Previously, we determined that the FTS monoclonal antibody, AD12.1, which detects a carbohydrate epitope on Wb-CFA, also detects multiple L. loa proteins in cross-reactive sera from persons with loiasis. Since the carbohydrate epitope recognized by AD12.1 is present on glycoproteins of other parasitic nematodes, including Brugia species, it is unclear why reactive glycoproteins are not detected in infections with other filarial parasites.

Methods

To gain a better understanding of the proteins recognized by the FTS diagnostic antibody, we used proteomics and lectin array technology to characterize filarial glycoproteins that are bound by the AD12.1 antibody using Brugia malayi as a model.

Results

Distinct but overlapping sets of AD12 glycoproteins were identified from somatic and excretory/secretory worm products. One of the identified proteins, Bm18019 was confirmed as a secreted AD12-reactive glycoprotein by in-gel proteomics and immunoassays. Based on lectin binding patterns, Brugia AD12-reactive glycoproteins express glycans including core fucose, galactose, N-acetylglucosamine and galactose (β1-3)N-acetylgalactosamine in addition to the epitope recognized by AD12.1. None of the lectins that bound B. malayi AD12 glycoproteins had affinity for the Wb-CFA, highlighting a key difference between it and other AD12 glycoproteins.

Conclusions

B. malayi somatic and excretory/secretory proteins are similar to L. loa antigens found in FTS-positive human sera, bolstering the hypothesis that circulating L. loa AD12 antigens result from worm tissue damage or death. The difference in glycan and protein composition between the Wb-CFA and other AD12 glycoproteins can be used to differentiate LF from cross-reactive loiasis.

Heartworm disease - Overview, intervention, and industry perspective

Dirofilaria immitis, also known as heartworm, is a major parasitic threat for dogs and cats around the world. Because of its impact on the health and welfare of companion animals, heartworm disease is of huge veterinary and economic importance especially in North America, Europe, Asia and Australia. Within the animal health market many different heartworm preventive products are available, all of which contain active components of the same drug class, the macrocyclic lactones. In addition to compliance issues, such as under-dosing or irregular treatment intervals, the occurrence of drug-resistant heartworms within the populations in the Mississippi River areas adds to the failure of preventive treatments. The objective of this review is to provide an overview of the disease, summarize the current disease control measures and highlight potential new avenues and best practices for treatment and prevention.

Transcriptome Analysis and miRNA Target Profiling at Various Stages of Root-Knot Nematode Meloidogyne incognita Development for Identification of Potential Regulatory Networks

Root-knot nematodes (RKNs) are a group of plant-parasitic nematodes that cause damage to various plant species and extensive economical losses. In this study, we performed integrated analysis of miRNA and mRNA expression data to explore the regulation of miRNA and mRNA in RKNs. In particular, we aimed to elucidate the mRNA targets of Meloidogyne incognita miRNAs and variations of the RKN transcriptome during five stages of its life cycle. Stage-wise RNA sequencing of M. incognita resulted in clean read numbers of 56,902,902, 50,762,456, 40,968,532, 47,309,223, and 51,730,234 for the egg, J2, J3, J4, and female stages, respectively. Overall, stage-dependent mRNA sequencing revealed that 17,423 genes were expressed in the transcriptome of M. incognita. The egg stage showed the maximum number of transcripts, and 12,803 gene transcripts were expressed in all stages. Functional Gene Ontology (GO) analysis resulted in three main GO classes: biological process, cellular components, and molecular function; the detected sequences were longer than sequences in the reference genome. Stage-wise selected fragments per kilobase of transcript per million mapped reads (FPKM) values of the top 10 stage-specific and common mRNAs were used to construct expression profiles, and 20 mRNAs were validated through quantitative real-time PCR (qRT-PCR). Next, we used three target prediction programs (miRanda, RNAhybrid, and PITA) to obtain 2431 potential target miRNA genes in RKNs, which regulate 8331 mRNAs. The predicted potential targets of miRNA were generally involved in cellular and metabolic processes, binding of molecules in the cell, membranes, and organelles. Stage-wise miRNA target analysis revealed that the egg stage contains heat shock proteins, transcriptional factors, and DNA repair proteins, whereas J2 includes DNA replication, heat shock, and ubiquitin-conjugating pathway-related proteins; the J3 and J4 stages are represented by the major sperm protein domain and translation-related proteins, respectively. In the female stage, we found proteins related to the maintenance of molybdopterin-binding domain-containing proteins and ubiquitin-mediated protein degradation. In total, 29 highly expressed stage-specific mRNA-targeting miRNAs were analyzed using qRT-PCR to validate the sequence analysis data. Overall, our findings provide new insights into the molecular mechanisms occurring at various developmental stages of the RKN life cycle, thus aiding in the identification of potential control strategies.

Parasitic Helminth-Derived microRNAs and Extracellular Vesicle Cargos as Biomarkers for Helminthic Infections

As an adaption to their complex lifecycles, helminth parasites garner a unique repertoire of genes at different developmental stages with subtle regulatory mechanisms. These parasitic worms release differential components such as microRNAs (miRNAs) and extracellular vesicles (EVs) as mediators which participate in the host-parasite interaction, immune regulation/evasion, and in governing processes associated with host infection. MiRNAs are small (~ 22-nucleotides) non-coding RNAs that regulate gene expression at the post-transcriptional level, and can exist in stable form in bodily fluids such as serum/plasma, urine, saliva and bile. In addition to reports focusing on the identification of miRNAs or in the probing of differentially expressed miRNA profiles in different development stages/sexes or in specific tissues, a number of studies have focused on the detection of helminth-derived miRNAs in the mammalian host circulatory system as diagnostic biomarkers. Extracellular vesicles (EVs), small membrane-surrounded structures secreted by a wide variety of cell types, contain rich cargos that are important in cell-cell communication. EVs have attracted wide attention due to their unique functional relevance in host-parasite interactions and for their potential value in translational applications such as biomarker discovery. In the current review, we discuss the status and potential of helminth parasite-derived circulating miRNAs and EV cargos as novel diagnostic tools.

Setaria cervi (Filarioidea, Onchocercidae) undressing in ungulates: altered morphology of developmental stages, their molecular detection and complete sequence cox1 gene

This work introduces new morphological and molecular information on the filaroid nematode Setaria cervi (Rudolphi, 1819) obtained from 13 infected game ungulates out of 96 dissected. The hosts comprised the following: a single moose (Alces alces), ten red deer (Cervus elaphus) and two sika deer (Cervus nippon) originating from the western and northern regions of the Czech Republic. Based on the complete sequences of the gene encoding mitochondrial cytochrome c oxidase subunit 1 (cox1), all 20 females and four males belonged to the species S. cervi. We detected three developmental female stages (adult fertile females, juvenile L5 females and L4 female larvae) differing in size and some morphological traits as the subtle structure of peribuccal crown and shape and features of tail knob. Such differences were described in detail for the first time. The phylogenetic relationships within the family Onchocercidae have been evaluated using new information on the cox1 sequence of S. cervi (maximum likelihood method, GTR + I + G model). In accordance with the latest phylogenetic studies, the present analysis confirmed the ancient separation of the subclass Setariinae from the remaining two onchocercid lineages Dirofilariinae and Onchocerinae.

Genomics of the Parasitic Nematode Ascaris and Its Relatives

Nematodes of the genus Ascaris are important parasites of humans and swine, and the phylogenetically related genera (Parascaris, Toxocara, and Baylisascaris) infect mammals of veterinary interest. Over the last decade, considerable genomic resources have been established for Ascaris, including complete germline and somatic genomes, comprehensive mRNA and small RNA transcriptomes, as well as genome-wide histone and chromatin data. These datasets provide a major resource for studies on the basic biology of these parasites and the host-parasite relationship. Ascaris and its relatives undergo programmed DNA elimination, a highly regulated process where chromosomes are fragmented and portions of the genome are lost in embryonic cells destined to adopt a somatic fate, whereas the genome remains intact in germ cells. Unlike many model organisms, Ascaris transcription drives early development beginning prior to pronuclear fusion. Studies on Ascaris demonstrated a complex small RNA network even in the absence of a piRNA pathway. Comparative genomics of these ascarids has provided perspectives on nematode sex chromosome evolution, programmed DNA elimination, and host-parasite coevolution. The genomic resources enable comparison of proteins across diverse species, revealing many new potential drug targets that could be used to control these parasitic nematodes.

Extensive non-redundancy in a recently duplicated developmental gene family

BACKGROUND

It has been proposed that recently duplicated genes are more likely to be redundant with one another compared to ancient paralogues. The evolutionary logic underpinning this idea is simple, as the assumption is that recently derived paralogous genes are more similar in sequence compared to members of ancient gene families. We set out to test this idea by using molecular phylogenetics and exploiting the genetic tractability of the model nematode, Caenorhabditis elegans, in studying the nematode-specific family of Hedgehog-related genes, the Warthogs. Hedgehog is one of a handful of signal transduction pathways that underpins the development of bilaterian animals. While having lost a bona fide Hedgehog gene, most nematodes have evolved an expanded repertoire of Hedgehog-related genes, ten of which reside within the Warthog family.

RESULTS

We have characterised their evolutionary origin and their roles in C. elegans and found that these genes have adopted new functions in aspects of post-embryonic development, including left-right asymmetry and cell fate determination, akin to the functions of their vertebrate counterparts. Analysis of various double and triple mutants of the Warthog family reveals that more recently derived paralogues are not redundant with one another, while a pair of divergent Warthogs do display redundancy with respect to their function in cuticle biosynthesis.

CONCLUSIONS

We have shown that newer members of taxon-restricted gene families are not always functionally redundant despite their recent inception, whereas much older paralogues can be, which is considered paradoxical according to the current framework in gene evolution.

Mitochondrial Genome of Brugia malayi Microfilariae Isolated From a Clinical Sample

Lymphatic filariasis is a neglected parasitic disease that is a leading cause of long-term disability. Information obtained from genome sequencing of filarial worm can help us identify systems in the worm that are likely to be useful for novel drug design. Brugia (B.) malayi is still the only lymphatic-dwelling filarial parasite with a nearly complete, fully annotated, and published genome. However, most previous studies were based on the FR3 strain of B. malayi, which originally was isolated from a human patient, and was adapted to the rodent model, then maintained in laboratories for more than 60 years. It is uncertain whether genetic variation exists, thus, sequencing of clinical isolates of lymphatic dwelling filarial parasites is a high priority. Here, we report for the first time the complete mitochondrial genome of B. malayi microfilariae from clinical isolate. Complete mitochondrial (mt) genome of the microfilariae isolated from a blood sample taken from a Thai subject living in Narathiwat Province, which is an endemic area of brugian filariasis, was assembled with sequencing reads obtained by Illumina sequencing. Gene annotation, phylogenetic analysis and single nucleotide polymorphism (SNP) were deployed. A complete 13,658-bp mt genome of B. malayi microfilaria was obtained, and it shows 68x coverage. Based on gene annotation, the mt genome consists of 12 protein-coding, two rRNA, and 23 tRNA genes. Phylogenetic analysis using all protein sequences of DNA sequences of mt genome or cytochrome c oxidase subunit I (COX1) revealed a close relationship among three lymphatic filariae (i.e., B. timori, zoonotic B. pahangi, and Wuchereria spp.). The SNPs in the COX1 gene can differentiate microfilariae of B. malayi in human from those found in canine. Furthermore, the number, order and transcription, and direction of B. malayi microfilariae mitochondrial genes were the same as those found in the FR3 strain of B. malayi. The comparison on mitochondrial genome of B. malayi could have important implications on the development of a new intervention or vaccine to treat or prevent this disease in endemic areas/regions around the world.

Dual RNAseq analyses at soma and germline levels reveal evolutionary innovations in the elephantiasis-agent Brugia malayi, and adaptation of its Wolbachia endosymbionts

Brugia malayi is a human filarial nematode responsible for elephantiasis, a debilitating condition that is part of a broader spectrum of diseases called filariasis, including lymphatic filariasis and river blindness. Almost all filarial nematode species infecting humans live in mutualism with Wolbachia endosymbionts, present in somatic hypodermal tissues but also in the female germline which ensures their vertical transmission to the nematode progeny. These α-proteobacteria potentially provision their host with essential metabolites and protect the parasite against the vertebrate immune response. In the absence of Wolbachia wBm, B. malayi females become sterile, and the filarial nematode lifespan is greatly reduced. In order to better comprehend this symbiosis, we investigated the adaptation of wBm to the host nematode soma and germline, and we characterized these cellular environments to highlight their specificities. Dual RNAseq experiments were performed at the tissue-specific and ovarian developmental stage levels, reaching the resolution of the germline mitotic proliferation and meiotic differentiation stages. We found that most wBm genes, including putative effectors, are not differentially regulated between infected tissues. However, two wBm genes involved in stress responses are upregulated in the hypodermal chords compared to the germline, indicating that this somatic tissue represents a harsh environment to which wBm have adapted. A comparison of the B. malayi and C. elegans germline transcriptomes reveals a poor conservation of genes involved in the production of oocytes, with the filarial germline proliferative zone relying on a majority of genes absent from C. elegans. The first orthology map of the B. malayi genome presented here, together with tissue-specific expression enrichment analyses, indicate that the early steps of oogenesis are a developmental process involving genes specific to filarial nematodes, that likely result from evolutionary innovations supporting the filarial parasitic lifestyle.

Resolution of polycistronic RNA by SL2 trans-splicing is a widely conserved nematode trait

Spliced leader trans-splicing is essential for the processing and translation of polycistronic RNAs generated by eukaryotic operons. In C. elegans, a specialized spliced leader, SL2, provides the 5' end for uncapped pre-mRNAs derived from polycistronic RNAs. Studies of other nematodes suggested that SL2-type trans-splicing is a relatively recent innovation, confined to Rhabditina, the clade containing C. elegans and its close relatives. Here we conduct a survey of transcriptome-wide spliced leader trans-splicing in Trichinella spiralis, a distant relative of C. elegans with a particularly diverse repertoire of 15 spliced leaders. By systematically comparing the genomic context of trans-splicing events for each spliced leader, we identified a subset of T. spiralis spliced leaders that are specifically used to process polycistronic RNAs-the first examples of SL2-type spliced leaders outside of Rhabditina. These T. spiralis spliced leader RNAs possess a perfectly conserved stem-loop motif previously shown to be essential for SL2-type trans-splicing in C. elegans We show that genes trans-spliced to these SL2-type spliced leaders are organized in operonic fashion, with short intercistronic distances. A subset of T. spiralis operons show conservation of synteny with C. elegans operons. Our work substantially revises our understanding of nematode spliced leader trans-splicing, showing that SL2 trans-splicing is a major mechanism for nematode polycistronic RNA processing, which may have evolved prior to the radiation of the Nematoda. This work has important implications for the improvement of genome annotation pipelines in nematodes and other eukaryotes with operonic gene organization.

Molecular evidence of hybridization between pig and human Ascaris indicates an interbred species complex infecting humans

Human ascariasis is a major neglected tropical disease caused by the nematode Ascaris lumbricoides. We report a 296 megabase (Mb) reference-quality genome comprised of 17,902 protein-coding genes derived from a single, representative Ascaris worm. An additional 68 worms were collected from 60 human hosts in Kenyan villages where pig husbandry is rare. Notably, the majority of these worms (63/68) possessed mitochondrial genomes that clustered closer to the pig parasite Ascaris suum than to A. lumbricoides. Comparative phylogenomic analyses identified over 11 million nuclear-encoded SNPs but just two distinct genetic types that had recombined across the genomes analyzed. The nuclear genomes had extensive heterozygosity, and all samples existed as genetic mosaics with either A. suum-like or A. lumbricoides-like inheritance patterns supporting a highly interbred Ascaris species genetic complex. As no barriers appear to exist for anthroponotic transmission of these 'hybrid' worms, a one-health approach to control the spread of human ascariasis will be necessary.

CRISPR/Cas9: A new tool for the study and control of helminth parasites

Recent reports of CRISPR/Cas9 genome editing in parasitic helminths open up new avenues for research on these dangerous pathogens. However, the complex morphology and life cycles inherent to these parasites present obstacles for the efficient application of CRISPR/Cas9-targeted mutagenesis. This is especially true with the trematode flukes where only modest levels of gene mutation efficiency have been achieved. Current major challenges in the application of CRISPR/Cas9 for study of parasitic worms thus lie in enhancing gene mutation efficiency and overcoming issues involved in host passage so that mutated parasites survive. Strategies developed for CRISPR/Cas9 studies on Caenorhabditis elegans, protozoa and mammalian cells, including novel delivery methods, the choice of selectable markers, and refining mutation precision represent novel tactics whereby these impediments can be overcome. Furthermore, employing CRISPR/Cas9-mediated gene drive to interfere with vector transmission represents a novel approach for the control of parasitic worms that is worthy of further exploration.

Expression and functionality of allergenic genes regulated by simulated gastric juice in Anisakis pegreffii

To better understand the molecular mechanisms underlying allergens and parasite immunity and discover the stage-enriched gene expression of fish-borne zoonotic nematodes in the stomach, we used RNA-seq to study the transcriptome profiles of Anisakis pegreffii (Nematoda: Anisakidae, AP) in simulated gastric juice. Mobile L3 larvae were incubated in simulated medium at 37 °C in 5% CO₂ (AP-GJ) and the control group larvae were collected in PBS under the same conditions (AP-PBS). We found that the sequences of A. pegreffii were highly similar to Toxocara canis sequences. Among the transcripts, there would be 138 up-regulated putative genes and 251 down-regulated putative genes in AP-GJ group. Several lipid binging-related genes were more highly expressed in AP-GJ larvae. Moreover, 17 allergen genes were up-regulated and 29 were down-regulated in AP-GJ larvae. Eleven allergen genes belonged to one or more of the following three categories: biological process, cellular component, and molecular function. According to KEGG analysis, the main pathways that were represented included protein processing in transcription, immune system, cancer, and infectious disease. In particular, the most significant changes in the expression of parasite-derived allergen products occurred in AP-GJ larvae. This study helps us to extend our understanding of the biology of the fish-borne zoonotic parasite A. pegreffii and could be helpful for more precise risk assessment and providing guidelines for allergic consumers.

Comprehensive Chromosome End Remodeling during Programmed DNA Elimination

Germline and somatic genomes are in general the same in a multicellular organism. However, programmed DNA elimination leads to a reduced somatic genome compared to germline cells. Previous work on the parasitic nematode Ascaris demonstrated that programmed DNA elimination encompasses high-fidelity chromosomal breaks and loss of specific genome sequences including a major tandem repeat of 120 bp and ~1,000 germline-expressed genes. However, the precise chromosomal locations of these repeats, breaks regions, and eliminated genes remained unknown. We used PacBio long-read sequencing and chromosome conformation capture (Hi-C) to obtain fully assembled chromosomes of Ascaris germline and somatic genomes, enabling a complete chromosomal view of DNA elimination. We found that all 24 germline chromosomes undergo comprehensive chromosome end remodeling with DNA breaks in their subtelomeric regions and loss of distal sequences including the telomeres at both chromosome ends. All new Ascaris somatic chromosome ends are recapped by de novo telomere healing. We provide an ultrastructural analysis of Ascaris DNA elimination and show that eliminated DNA is incorporated into double membrane-bound structures, similar to micronuclei, during telophase of a DNA elimination mitosis. These micronuclei undergo dynamic changes including loss of active histone marks and localize to the cytoplasm following daughter nuclei formation and cytokinesis where they form autophagosomes. Comparative analysis of nematode chromosomes suggests that chromosome fusions occurred, forming Ascaris sex chromosomes that become independent chromosomes following DNA elimination breaks in somatic cells. These studies provide the first chromosomal view and define novel features and functions of metazoan programmed DNA elimination.

Structure, development, and evolutive patterns of spermatozoa in rhabditid nematodes (Nematoda: Rhabditida)

Spermatogenesis of five rhabditid nematodes was studied using transmission electron microscopy and is described herein. Structure and development of nematode sperm in all available representatives of the extensive order Rhabditida have been analysed and the main characteristics of each infraorder are discussed. The ancestral sperm of the order Rhabditida was reconstructed using maximum likelihood and Bayesian methods based on 44 ultrastructural sperm characters. The hypothetical ancestral spermatogenesis of the order Rhabditida agrees with the previously suggested "rhabditid" pattern and appears to be conserved throughout the order Rhabditida. Despite the enormous variation of rhabditid nematodes, few groups deviate from the ancestral pattern. This conserved pattern can be informative within the phylum Nematoda at order level, but poses limitations when used in taxonomic and phylogenetic analysis within Rhabditida.

Modeling the metabolic interplay between a parasitic worm and its bacterial endosymbiont allows the identification of novel drug targets

The filarial nematode Brugia malayi represents a leading cause of disability in the developing world, causing lymphatic filariasis in nearly 40 million people. Currently available drugs are not well-suited to mass drug administration efforts, so new treatments are urgently required. One potential vulnerability is the endosymbiotic bacteria Wolbachia-present in many filariae-which is vital to the worm. Genome scale metabolic networks have been used to study prokaryotes and protists and have proven valuable in identifying therapeutic targets, but have only been applied to multicellular eukaryotic organisms more recently. Here, we present iDC625, the first compartmentalized metabolic model of a parasitic worm. We used this model to show how metabolic pathway usage allows the worm to adapt to different environments, and predict a set of 102 reactions essential to the survival of B. malayi. We validated three of those reactions with drug tests and demonstrated novel antifilarial properties for all three compounds.

Advances in Omic Studies Drive Discoveries in the Biology of Anisakid Nematodes

Advancements in technologies employed in high-throughput next-generation sequencing (NGS) methods are supporting the spread of studies that, combined with advances in computational biology and bioinformatics, have greatly accelerated discoveries within basic and biomedical research for many parasitic diseases. Here, we review the most updated "omic" studies performed on anisakid nematodes, a family of marine parasites that are causative agents of the fish-borne zoonosis known as anisakiasis or anisakidosis. Few deposited data on Anisakis genomes are so far available, and this still hinders the deep and highly accurate characterization of biological aspects of interest, even as several transcriptomic and proteomic studies are becoming available. These have been aimed at discovering and characterizing molecules specific to peculiar developmental parasitic stages or tissues, as well as transcripts with pathogenic potential as toxins and allergens, with a broad relevance for a better understanding of host-pathogen relationships and for the development of reliable diagnostic tools.