The genome of Onchocerca volvulus, agent of river blindness

Fexinidazole and Corallopyronin A target Wolbachia -infected sheath cells present in filarial nematodes

The discovery of the endosymbiotic bacteria Wolbachia as an obligate symbiont of filarial nematodes has led to antibiotic-based treatments for filarial diseases. While lab and clinical studies have yielded promising results, recent animal studies reveal that Wolbachia levels may rebound following treatment with suboptimal doses of the antibiotic rifampicin. Previous work showed that a likely source of the bacterial rebound in females were dense clusters of Wolbachia in ovarian tissue. The number, size, and density of these Wolbachia clusters were not diminished despite antibiotic treatment. Here we define the cellular characteristics of the Wolbachia clusters in Brugia pahangi (wBp) and identify drugs that also target them. We have evidence that the Wolbachia clusters originate from newly formed sheath cells adjacent to the ovarian Distal Tip Cells. The dramatically enlarged volume of an infected sheath cell is strikingly similar to endosymbiont-induced bacteriocytes found in many insect species. Ultrastructural analysis reveals that the clustered Wolbachia present within the sheath cells exhibit a distinct morphology and form direct connections with the oocyte membrane and possibly the cytoplasm. This includes membrane-based channels providing a connection between Wolbachia -infected sheath cells and oocytes. We also determined that the Wolbachia within the sheath cells are either quiescent or replicating at a very low rate. Screens of known antibiotics and other drugs revealed that two drugs, Fexinidazole and Corallopyronin A, significantly reduced the number of clustered Wolbachia located within the sheath cells.

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.

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.

Diverse RNA viruses of parasitic nematodes can elicit antibody responses in vertebrate hosts

Parasitic nematodes have an intimate, chronic and lifelong exposure to vertebrate tissues. Here we mined 41 published parasitic nematode transcriptomes from vertebrate hosts and identified 91 RNA viruses across 13 virus orders from 24 families in ~70% (28 out of 41) of parasitic nematode species, which include only 5 previously reported viruses. We observe widespread distribution of virus-nematode associations across multiple continents, suggesting an ancestral acquisition event and host-virus co-evolution. Characterization of viruses of Brugia malayi (BMRV1) and Onchocerca volvulus (OVRV1) shows that these viruses are abundant in reproductive tissues of adult parasites. Importantly, the presence of BMRV1 RNA in B. malayi parasites mounts an RNA interference response against BMRV1 suggesting active viral replication. Finally, BMRV1 and OVRV1 were found to elicit antibody responses in serum samples from infected jirds and infected or exposed humans, indicating direct exposure to the immune system.

Spirocerca lupi draft genome, vaccine and anthelmintic targets

Spirocerca lupi is a parasitic nematode affecting predominantly domestic dogs. It causes spirocercosis, a disease that is often fatal. The assembled draft genome of S. lupi consists of 13,627 predicted protein-coding genes and is approximately 150 Mb in length. Several known anthelmintic gene targets such as for β-Tubulin, glutamate, and GABA receptors as well as known vaccine gene targets such as cysteine protease inhibitor and cytokines were identified in S. lupi by comparing orthologs of C. elegans anthelmintic gene targets as well as orthologs to known vaccine candidates. New anthelmintic targets were predicted through an inclusion-exclusion strategy and new vaccine targets were predicted through an immunoinformatics approach. New anthelminthic targets include DNA-directed RNA polymerases, chitin synthase, polymerases, and other enzymes. New vaccine targets include cuticle collagens. These gene targets provide a starting platform for new drug identification and vaccine design.

Comparative development of human filariae Loa loa, Onchocerca volvulus and Mansonella perstans in immunocompromised mouse strains

Introduction

Mouse models of human filarial infections are not only urgently needed to investigate the biology of the nematodes and their modulation of the host's immunity, but will also provide a platform to screen and test novel anti-filarial drugs. Recently, murine Loa loa infection models have been stablished using immunocompromised mouse strains, whereas murine Mansonella perstans infections have not been implemented until now.

Methods

Therefore, we aim to establish experimental M. perstans infections using the immunocompromised mouse strains RAG2IL-2Rγ-/- (lack B, T and natural killer cells), IL-4Rα/IL-5-/- (impaired IL-4/5 signalling and eosinophil activation) and NOD.Cg-PrkdcscidIl2rgtm1Wj l/SzJ (NOD scid gamma, NSG) BALB/c mice (lack mature lymphocytes) through subcutaneous (s.c.) or intraperitoneal (i.p.) inoculation of infective stage 3 larvae (L3) isolated from engorged vectors.

Results

In total, 145 immunocompromised mice have been inoculated with 3,250 M. perstans, 3,337 O. volvulus, and 2,720 Loa loa L3 to comparatively analyse which immunocompromised mouse strain is susceptible to human filarial infections. Whereas, no M. perstans and O. volvulus L3 could be recovered upon 2-63 days post-inoculation, a 62-66% Loa loa L3 recovery rate could be achieved in the different mouse strains. Gender of mice, type of inoculation (s.c. or i.p.) or time point of analysis (2-63 days post inoculation) did not interfere with the success of L3 recovery. In addition, administration of the immune suppressants hydrocortisone, prednisolone and cyclophosphamide did not restore M. perstans L3 recovery rates.

Discussion

These findings show that RAG2IL-2Rg-/-BALB/c and C57BL/6, IL-4Rα/IL-5-/- BALB/c and NSG mice were not susceptible to M. perstans and O. volvulus L3 inoculation using the applied methods, whereas Loa loa infection could be maintained. Further studies should investigate if humanized immunocompromised mice might be susceptible to M. perstans. and O. volvulus.

Genome structure and population genomics of the canine heartworm Dirofilaria immitis

The heartworm, Dirofilaria immitis, is a filarial parasitic nematode responsible for significant morbidity and mortality in wild and domesticated canids. Resistance to macrocyclic lactone drug prevention represents a significant threat to parasite control and has prompted investigations to understand the genetic determinants of resistance. This study aimed to improve the genomic resources of D. immitis to enable a more precise understanding of how genetic variation is distributed within and between parasite populations worldwide, which will inform the likelihood and rate by which parasites, and in turn, resistant alleles, might spread. We have guided the scaffolding of a recently published genome assembly for D. immitis (ICBAS_JMDir_1.0) using the chromosomal-scale reference genomes of Brugia malayi and Onchocerca volvulus, resulting in an 89.5 Mb assembly composed of four autosomal- and one sex-linked chromosomal-scale scaffolds representing 99.7% of the genome. Publicly available and new whole-genome sequencing data from 32 D. immitis samples from Australia, Italy and the USA were assessed using principal component analysis, nucleotide diversity (Pi) and absolute genetic divergence (Dxy) to characterise the global genetic structure and measure within- and between-population diversity. These population genetic analyses revealed broad-scale genetic structure among globally diverse samples and differences in genetic diversity between populations; however, fine-scale subpopulation analysis was limited and biased by differences between sample types. Finally, we mapped single nucleotide polymorphisms previously associated with macrocyclic lactone resistance in the new genome assembly, revealing the physical linkage of high-priority variants on chromosome 3, and determined their frequency in the studied populations. This new chromosomal assembly for D. immitis now allows for a more precise investigation of selection on genome-wide genetic variation and will enhance our understanding of parasite transmission and the spread of genetic variants responsible for resistance to treatment.

Advancing Strongyloides omics data: bridging the gap with Caenorhabditis elegans

Among nematodes, the free-living model organism Caenorhabditis elegans boasts the most advanced portfolio of high-quality omics data. The resources available for parasitic nematodes, including Strongyloides spp., however, are lagging behind. While C. elegans remains the most tractable nematode and has significantly advanced our understanding of many facets of nematode biology, C. elegans is not suitable as a surrogate system for the study of parasitism and it is important that we improve the omics resources available for parasitic nematode species. Here, we review the omics data available for Strongyloides spp. and compare the available resources to those for C. elegans and other parasitic nematodes. The advancements in C. elegans omics offer a blueprint for improving omics-led research in Strongyloides. We suggest areas of priority for future research that will pave the way for expansions in omics resources and technologies. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

The genome of Litomosoides sigmodontis illuminates the origins of Y chromosomes in filarial nematodes

Heteromorphic sex chromosomes are usually thought to have originated from a pair of autosomes that acquired a sex-determining locus and subsequently stopped recombining, leading to degeneration of the sex-limited chromosome. The majority of nematode species lack heteromorphic sex chromosomes and determine sex using an X-chromosome counting mechanism, with males being hemizygous for one or more X chromosomes (XX/X0). Some filarial nematode species, including important parasites of humans, have heteromorphic XX/XY karyotypes. It has been assumed that sex is determined by a Y-linked locus in these species. However, karyotypic analyses suggested that filarial Y chromosomes are derived from the unfused homologue of an autosome involved in an X-autosome fusion event. Here, we generated a chromosome-level reference genome for Litomosoides sigmodontis, a filarial nematode with the ancestral filarial karyotype and sex determination mechanism (XX/X0). By mapping the assembled chromosomes to the rhabditid nematode ancestral linkage (or Nigon) elements, we infer that the ancestral filarial X chromosome was the product of a fusion between NigonX (the ancestrally X-linked element) and NigonD (ancestrally autosomal). In the two filarial lineages with XY systems, there have been two independent X-autosome chromosome fusion events involving different autosomal Nigon elements. In both lineages, the region shared by the neo-X and neo-Y chromosomes is within the ancestrally autosomal portion of the X, confirming that the filarial Y chromosomes are derived from the unfused homologue of the autosome. Sex determination in XY filarial nematodes therefore likely continues to operate via the ancestral X-chromosome counting mechanism, rather than via a Y-linked sex-determining locus.

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

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

Optimized strategy for real-time qPCR detection of Onchocerca volvulus DNA in pooled Simulium sp. blackfly vectors

BACKGROUND

Onchocerca volvulus is a filarial parasite that is a major cause of dermatitis and blindness in endemic regions primarily in sub-Saharan Africa. Widespread efforts to control the disease caused by O. volvulus infection (onchocerciasis) began in 1974 and in recent years, following successful elimination of transmission in much of the Americas, the focus of efforts in Africa has moved from control to the more challenging goal of elimination of transmission in all endemic countries. Mass drug administration (MDA) with ivermectin has reached more than 150 million people and elimination of transmission has been confirmed in four South American countries, with at least two African countries having now stopped MDA as they approach verification of elimination. It is essential that accurate data for active transmission are used to assist in making the critical decision to stop MDA, since missing low levels of transmission and infection can lead to continued spread or recrudescence of the disease.

METHODOLOGY/PRINCIPAL FINDINGS

Current World Health Organization guidelines for MDA stopping decisions and post-treatment surveillance include screening pools of the Simulium blackfly vector for the presence of O. volvulus larvae using a PCR-ELISA-based molecular technique. In this study, we address the potential of an updated, practical, standardized molecular diagnostic tool with increased sensitivity and species-specificity by comparing several candidate qPCR assays. When paired with heat-stable reagents, a qPCR assay with a mitochondrial DNA target (OvND5) was found to be more sensitive and species-specific than an O150 qPCR, which targets a non-protein coding repetitive DNA sequence. The OvND5 assay detected 19/20 pools of 100 blackfly heads spiked with a single L3, compared to 16/20 for the O150 qPCR assay.

CONCLUSIONS/SIGNIFICANCE

Given the improved sensitivity, species-specificity and resistance to PCR inhibitors, we identified OvND5 as the optimal target for field sample detection. All reagents for this assay can be shipped at room temperature with no loss of activity. The qPCR protocol we propose is also simpler, faster, and more cost-effective than the current end-point molecular assays.

Ancient diversity in host-parasite interaction genes in a model parasitic nematode

Host-parasite interactions exert strong selection pressures on the genomes of both host and parasite. These interactions can lead to negative frequency-dependent selection, a form of balancing selection that is hypothesised to explain the high levels of polymorphism seen in many host immune and parasite antigen loci. Here, we sequence the genomes of several individuals of Heligmosomoides bakeri, a model parasite of house mice, and Heligmosomoides polygyrus, a closely related parasite of wood mice. Although H. bakeri is commonly referred to as H. polygyrus in the literature, their genomes show levels of divergence that are consistent with at least a million years of independent evolution. The genomes of both species contain hyper-divergent haplotypes that are enriched for proteins that interact with the host immune response. Many of these haplotypes originated prior to the divergence between H. bakeri and H. polygyrus, suggesting that they have been maintained by long-term balancing selection. Together, our results suggest that the selection pressures exerted by the host immune response have played a key role in shaping patterns of genetic diversity in the genomes of parasitic nematodes.

Assessing Onchocerca volvulus Intensity of Infection and Genetic Diversity Using Mitochondrial Genome Sequencing of Single Microfilariae Obtained before and after Ivermectin Treatment

Onchocerciasis is a neglected tropical disease targeted for elimination using ivermectin mass administration. Ivermectin kills the microfilariae and temporarily arrests microfilariae production by the macrofilariae. We genotyped 436 microfilariae from 10 people each in Ituri, Democratic Republic of the Congo (DRC), and Maridi County, South Sudan, collected before and 4-5 months after ivermectin treatment. Population genetic analyses identified 52 and 103 mitochondrial DNA haplotypes among the microfilariae from DRC and South Sudan, respectively, with few haplotypes shared between people. The percentage of genotype-based correct assignment to person within DRC was ~88% and within South Sudan ~64%. Rarefaction and extrapolation analysis showed that the genetic diversity in DRC, and even more so in South Sudan, was captured incompletely. The results indicate that the per-person adult worm burden is likely higher in South Sudan than DRC. Analyses of haplotype data from a subsample (n = 4) did not discriminate genetically between pre- and post-treatment microfilariae, confirming that post-treatment microfilariae are not the result of new infections. With appropriate sampling, mitochondrial haplotype analysis could help monitor changes in the number of macrofilariae in a population as a result of treatment, identify cases of potential treatment failure, and detect new infections as an indicator of continuing transmission.

MetaPro: a scalable and reproducible data processing and analysis pipeline for metatranscriptomic investigation of microbial communities

BACKGROUND

Whole microbiome RNASeq (metatranscriptomics) has emerged as a powerful technology to functionally interrogate microbial communities. A key challenge is how best to process, analyze, and interpret these complex datasets. In a typical application, a single metatranscriptomic dataset may comprise from tens to hundreds of millions of sequence reads. These reads must first be processed and filtered for low quality and potential contaminants, before being annotated with taxonomic and functional labels and subsequently collated to generate global bacterial gene expression profiles.

RESULTS

Here, we present MetaPro, a flexible, massively scalable metatranscriptomic data analysis pipeline that is cross-platform compatible through its implementation within a Docker framework. MetaPro starts with raw sequence read input (single-end or paired-end reads) and processes them through a tiered series of filtering, assembly, and annotation steps. In addition to yielding a final list of bacterial genes and their relative expression, MetaPro delivers a taxonomic breakdown based on the consensus of complementary prediction algorithms, together with a focused breakdown of enzymes, readily visualized through the Cytoscape network visualization tool. We benchmark the performance of MetaPro against two current state-of-the-art pipelines and demonstrate improved performance and functionality.

CONCLUSIONS

MetaPro represents an effective integrated solution for the processing and analysis of metatranscriptomic datasets. Its modular architecture allows new algorithms to be deployed as they are developed, ensuring its longevity. To aid user uptake of the pipeline, MetaPro, together with an established tutorial that has been developed for educational purposes, is made freely available at https://github.com/ParkinsonLab/MetaPro . The software is freely available under the GNU general public license v3. Video Abstract.

Morphological and molecular characterization of Haemonchus contortus isolated from the small ruminants of south Gujarat, India

The successful design of strategic control measures against the blood-sucking gastrointestinal nematode, Haemonchus contortus in small ruminants can be facilitated by revealing its general features from morphology to the molecular level. In the south Gujarat region of India, a total of 2408 H. contortus were collected from 84 slaughtered sheep's abomasum, consisting of 347 males and 2061 females (1:6 ratio) (p<0.05). Furthermore, 726 H. contortus were collected from 61 goats, comprising 145 males and 581 females (1:4 ratio) (p<0.05). The male worms were approximately 12±0.06 mm long, while female worms were about 20±0.09 mm long. The vulvar morphotypes of the female worms were found to be 17.7% linguiform, 76.6 % knobbed/button (p<0.05), and 5.7 % smooth type, demonstrating common features of H. contortus. The nucleotide sequences of the Internal Transcribed Spacer 1 (ITS-1) of 165 bp or ITS-2 plus of 256 bp were aligned, and it was found that the genotypes of male and female specimens of either sheep or goat origin were identical, with a 100 % match. The present isolates shared >95 % and >94 % homology with published sequences of ITS-1 and ITS-2 plus of H. contortus, respectively, with more nucleotide transitions than transversions in the aligned sequences. The reconstructed phylogram of either ITS-1 or ITS-2 plus revealed two major clades, one for H. contortus and another for other nematodes, with Haemonchus placei showing its proximity with the clade of H. contortus. The study established the role of morphological and molecular features in identifying and differentiating H. contortus parasite at the local level.

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.

Chromosome fusions repatterned recombination rate and facilitated reproductive isolation during Pristionchus nematode speciation

Large-scale genome-structural evolution is common in various organisms. Recent developments in speciation genomics revealed the importance of inversions, whereas the role of other genome-structural rearrangements, including chromosome fusions, have not been well characterized. We study genomic divergence and reproductive isolation of closely related nematodes: the androdioecious (hermaphroditic) model Pristionchus pacificus and its dioecious sister species Pristionchus exspectatus. A chromosome-level genome assembly of P. exspectatus using single-molecule and Hi-C sequencing revealed a chromosome-wide rearrangement relative to P. pacificus. Strikingly, genomic characterization and cytogenetic studies including outgroup species Pristionchus occultus indicated two independent fusions involving the same chromosome, ChrIR, between these related species. Genetic linkage analysis indicated that these fusions altered the chromosome-wide pattern of recombination, resulting in large low-recombination regions that probably facilitated the coevolution between some of the ~14.8% of genes across the entire genomes. Quantitative trait locus analyses for hybrid sterility in all three sexes revealed that major quantitative trait loci mapped to the fused chromosome ChrIR. While abnormal chromosome segregations of the fused chromosome partially explain hybrid female sterility, hybrid-specific recombination that breaks linkage of genes in the low-recombination region was associated with hybrid male sterility. Thus, recent chromosome fusions repatterned recombination rate and drove reproductive isolation during Pristionchus speciation.

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.

Rational design of a novel multi-epitope peptide-based vaccine against Onchocerca volvulus using transmembrane proteins

Almost a decade ago, it was recognized that the global elimination of onchocerciasis by 2030 will not be feasible without, at least, an effective prophylactic and/or therapeutic vaccine to complement chemotherapy and vector control strategies. Recent advances in computational immunology (immunoinformatics) have seen the design of novel multi-epitope onchocerciasis vaccine candidates which are however yet to be evaluated in clinical settings. Still, continued research to increase the pool of vaccine candidates, and therefore the chance of success in a clinical trial remains imperative. Here, we designed a multi-epitope vaccine candidate by assembling peptides from 14 O. volvulus (Ov) proteins using an immunoinformatics approach. An initial 126 Ov proteins, retrieved from the Wormbase database, and at least 90% similar to orthologs in related nematode species of economic importance, were screened for localization, presence of transmembrane domain, and antigenicity using different web servers. From the 14 proteins retained after the screening, 26 MHC-1 and MHC-II (T-cell) epitopes, and linear B-lymphocytes epitopes were predicted and merged using suitable linkers. The Mycobacterium tuberculosis Resuscitation-promoting factor E (RPFE_MYCTU), which is an agonist of TLR4, was then added to the N-terminal of the vaccine candidate as a built-in adjuvant. Immune simulation analyses predicted strong B-cell and IFN-γ based immune responses which are necessary for protection against O. volvulus infection. Protein-protein docking and molecular dynamic simulation predicted stable interactions between the 3D structure of the vaccine candidate and human TLR4. These results show that the designed vaccine candidate has the potential to stimulate both humoral and cellular immune responses and should therefore be subject to further laboratory investigation.

Identification of Novel Genes and Proteoforms in Angiostrongylus costaricensis through a Proteogenomic Approach

RNA sequencing (RNA-Seq) and mass-spectrometry-based proteomics data are often integrated in proteogenomic studies to assist in the prediction of eukaryote genome features, such as genes, splicing, single-nucleotide (SNVs), and single-amino-acid variants (SAAVs). Most genomes of parasite nematodes are draft versions that lack transcript- and protein-level information and whose gene annotations rely only on computational predictions. Angiostrongylus costaricensis is a roundworm species that causes an intestinal inflammatory disease, known as abdominal angiostrongyliasis (AA). Currently, there is no drug available that acts directly on this parasite, mostly due to the sparse understanding of its molecular characteristics. The available genome of A. costaricensis, specific to the Costa Rica strain, is a draft version that is not supported by transcript- or protein-level evidence. This study used RNA-Seq and MS/MS data to perform an in-depth annotation of the A. costaricensis genome. Our prediction improved the reference annotation with (a) novel coding and non-coding genes; (b) pieces of evidence of alternative splicing generating new proteoforms; and (c) a list of SNVs between the Brazilian (Crissiumal) and the Costa Rica strain. To the best of our knowledge, this is the first time that a multi-omics approach has been used to improve the genome annotation of A. costaricensis. We hope this improved genome annotation can assist in the future development of drugs, kits, and vaccines to treat, diagnose, and prevent AA caused by either the Brazil strain (Crissiumal) or the Costa Rica strain.

Development and validation of a quantitative PCR for the detection of Guinea worm (Dracunculus medinensis)

Dracunculus medinensis (Guinea worm) is a parasitic nematode that can cause the debilitating disease dracunculiasis (Guinea worm disease) in humans. The global Guinea Worm Eradication Program has led intervention and eradication efforts since the 1980s, and Guinea worm infections in people have decreased >99.99%. With the final goal of eradication drawing nearer, reports of animal infections from some remaining endemic countries pose unique challenges. Currently, confirmation of suspected Guinea worm infection relies on conventional molecular techniques such as polymerase chain reaction (PCR), which is not specific to Guinea worm and, therefore, requires sequencing of the PCR products to confirm the identity of suspect samples, a process that often takes a few weeks. To decrease the time required for species confirmation, we developed a quantitative PCR assay targeting the mitochondrial cytochrome b (cytb) gene of Guinea worm. Our assay has a limit of detection of 10 copies per reaction. The mean analytical parameters (± SE) were as follows: efficiency = 93.4 ± 7.7%, y-intercept = 40.93 ± 1.11, slope = -3.4896 ± 0.12, and the R² = 0.999 ± 0.004. The assay did not amplify other nematodes found in Guinea worm-endemic regions and demonstrated 100% diagnostic sensitivity and specificity. Implementation of this quantitative PCR assay for Guinea worm identification could eliminate the need for DNA sequencing to confirm species. Thus, this approach can be implemented to provide more rapid confirmation of Guinea worm infections, leading to faster execution of Guinea worm interventions while increasing our understanding of infection patterns.

Guinea worm (Dracunculus medinensis) is a parasitic nematode that causes debilitating disease in humans. The Guinea Worm Eradication Program would benefit from having a rapid molecular test that can confirm species identification without time-consuming DNA sequencing. We developed a qPCR protocol targeting the mitochondrial cytochrome b (cytb) gene of Guinea worm. The assay was validated analytically over 12 experiments using a standard serial dilution as well as diagnostically on DNA samples from non-target host species and other parasites (n = 180) and Guinea worm samples (n = 200) from a diversity of hosts and geographic regions. This assay could reliably detect 10 copies of the target DNA sequence and had a mean efficiency of 93.4% with 100% diagnostic sensitivity and specificity.

Improving helminth genome resources in the post-genomic era

Rapid advancement in high-throughput sequencing and analytical approaches has seen a steady increase in the generation of genomic resources for helminth parasites. Now, helminth genomes and their annotations are a cornerstone of numerous efforts to compare genetic and transcriptomic variation, from single cells to populations of globally distributed parasites, to genome modifications to understand gene function. Our understanding of helminths is increasingly reliant on these genomic resources, which are primarily static once published and vary widely in quality and completeness between species. This article seeks to highlight the cause and effect of this variation and argues for the continued improvement of these genomic resources - even after their publication - which is necessary to provide a more accurate and complete understanding of the biology of these important pathogens.

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

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

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

The genome sequence of the Australian filarial nematode, Cercopithifilaria johnstoni

We present a genome assembly and annotation of an individual female Cercopithifilaria johnstoni, a parasitic filarial nematode that is transmitted by hard ticks (Ixodidae) to infect a broad range of native Australian murid and marsupial hosts. The genome sequence is 76.9 Mbp in length, and although in draft form (N50 = 99 kbp, N50[n] = 232), is largely complete based on universally conserved orthologs (BUSCOs; genome = 94.9%, protein = 96.5%) and relative to other related filarial species. These data represent the first genomic resources for the genus Cercopithifilaria, a group of parasites with a broad host range, and form the basis for comparative analysis with the human-infective parasite, Onchocerca volvulus, both of which are responsible for similar eye and skin pathologies in their respective hosts.

The genome sequence of the Australian filarial nematode, Cercopithifilaria johnstoni

We present a genome assembly and annotation of an individual female Cercopithifilaria johnstoni, a parasitic filarial nematode that is transmitted by hard ticks (Ixodidae) to infect a broad range of native Australian murid and marsupial hosts. The genome sequence is 76.9 Mbp in length, and although in draft form (N50 = 99 kbp, N50[n] = 232), is largely complete based on universally conserved orthologs (BUSCOs; genome = 94.9%, protein = 96.5%) and relative to other related filarial species. These data represent the first genomic resources for the genus Cercopithifilaria, a group of parasites with a broad host range, and form the basis for comparative analysis with the human-infective parasite, Onchocerca volvulus, both of which are responsible for similar eye and skin pathologies in their respective hosts.

The evolution of the metazoan Toll receptor family and its expression during protostome development

BACKGROUND

Toll-like receptors (TLRs) play a crucial role in immunity and development. They contain leucine-rich repeat domains, one transmembrane domain, and one Toll/IL-1 receptor domain. TLRs have been classified into V-type/scc and P-type/mcc TLRs, based on differences in the leucine-rich repeat domain region. Although TLRs are widespread in animals, detailed phylogenetic studies of this gene family are lacking. Here we aim to uncover TLR evolution by conducting a survey and a phylogenetic analysis in species across Bilateria. To discriminate between their role in development and immunity we furthermore analyzed stage-specific transcriptomes of the ecdysozoans Priapulus caudatus and Hypsibius exemplaris, and the spiralians Crassostrea gigas and Terebratalia transversa.

RESULTS

We detected a low number of TLRs in ecdysozoan species, and multiple independent radiations within the Spiralia. V-type/scc and P-type/mcc type-receptors are present in cnidarians, protostomes and deuterostomes, and therefore they emerged early in TLR evolution, followed by a loss in xenacoelomorphs. Our phylogenetic analysis shows that TLRs cluster into three major clades: clade α is present in cnidarians, ecdysozoans, and spiralians; clade β in deuterostomes, ecdysozoans, and spiralians; and clade γ is only found in spiralians. Our stage-specific transcriptome and in situ hybridization analyses show that TLRs are expressed during development in all species analyzed, which indicates a broad role of TLRs during animal development.

CONCLUSIONS

Our findings suggest that a clade α TLR gene (TLR-Ca) and a clade β/γ TLR gene (TLR-Cβ/γ) were already present in the cnidarian-bilaterian common ancestor. However, although TLR-Ca was conserved in cnidarians, TLR-Cβ/γ was lost during the early evolution of these taxa. Moreover, TLR-Cβ/γ duplicated to generate TLR-Cβ and TLR-Cγ in the lineage to the last common protostome-deuterostome ancestor. TLR-Ca, TLR-Cβ and TLR-Cγ further expanded generating the three major TLR clades. While all three clades radiated in several spiralian lineages, specific TLRs clades have been presumably lost in other lineages. Furthermore, the expression of the majority of these genes during protostome ontogeny suggests a likely role in development.

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.

Development of emodepside as a possible adulticidal treatment for human onchocerciasis-The fruit of a successful industrial-academic collaboration

Current mass drug administration (MDA) programs for the treatment of human river blindness (onchocerciasis) caused by the filarial worm Onchocerca volvulus rely on ivermectin, an anthelmintic originally developed for animal health. These treatments are primarily directed against migrating microfilariae and also suppress fecundity for several months, but fail to eliminate adult O. volvulus. Therefore, elimination programs need time frames of decades, well exceeding the life span of adult worms. The situation is worsened by decreased ivermectin efficacy after long-term therapy. To improve treatment options against onchocerciasis, a drug development candidate should ideally kill or irreversibly sterilize adult worms. Emodepside is a broad-spectrum anthelmintic used for the treatment of parasitic nematodes in cats and dogs (Profender and Procox). Our current knowledge of the pharmacology of emodepside is the result of more than 2 decades of intensive collaborative research between academia and the pharmaceutical industry. Emodepside has a novel mode of action with a broad spectrum of activity, including against extraintestinal nematode stages such as migrating larvae or macrofilariae. Therefore, emodepside is considered to be among the most promising candidates for evaluation as an adulticide treatment against onchocerciasis. Consequently, in 2014, Bayer and the Drugs for Neglected Diseases initiative (DNDi) started a collaboration to develop emodepside for the treatment of patients suffering from the disease. Macrofilaricidal activity has been demonstrated in various models, including Onchocerca ochengi in cattle, the parasite most closely related to O. volvulus. Emodepside has now successfully passed Phase I clinical trials, and a Phase II study is planned. This Bayer–DNDi partnership is an outstanding example of “One World Health,” in which experience gained in veterinary science and drug development is translated to human health and leads to improved tools to combat neglected tropical diseases (NTDs) and shorten development pathways and timelines in an otherwise neglected area.

Onchocerca volvulus is the causative agent of human river blindness, and current elimination programs rely on the use of ivermectin to kill microfilariae. Since no adulticidal drug is available and adult worms have a life span of up to 15 years, elimination programs need to be sustained over several decades. Emodepside is an anthelmintic that is licensed as a dewormer for cats and dogs. Due to its ability to eliminate nematodes located in various extraintestinal host tissues, including migrating larvae and adult filarial worms, it is considered to be an excellent candidate for the treatment of onchocerciasis. Intense collaboration between academia and the pharmaceutical industry has led to a deep understanding of the novel mode of action of the drug and of its parasite target spectrum. Phase I clinical trials with emodepside have demonstrated its safety and adulticide activity against the closely related cattle parasite Onchocerca ochengi. Currently, Phase II clinical trials are planned to confirm that emodepside, developed initially to improve animal health, has also the potential to improve human health by tackling a very important neglected tropical disease (NTD).

Abdominal angiostrongyliasis in the Americas: fifty years since the discovery of a new metastrongylid species, Angiostrongylus costaricensis

Angiostrongylus costaricensis is a zoonotic parasitic nematode described for the first time in 1971 by Pedro Morera and Rodolfo Céspedes in Costa Rica. This parasite causes an infection known as abdominal angiostrongyliasis, affecting mainly school-aged children and young adults. Infection with A. costaricensis has been associated with a myriad of rodent and mollusk species in the Americas and the Caribbean, as its natural hosts and reservoirs. In this commemorative review, we highlight the extensive research collected through a 50-year journey, which includes ecological, pathological, and molecular studies on A. costaricensis and its implicated disease. We also identify major knowledge gaps in its evolutionary history, the ecological role of imported and invasive mollusk species, and immune response. We propose that the advent of -omics analyses will allow us to gather novel information regarding A. costaricensis biology and infection dynamics, as well as to promote the design of much-needed sensitive and specific diagnostic tools.

A telomere-to-telomere assembly of Oscheius tipulae and the evolution of rhabditid nematode chromosomes

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

Triploid pregnancy-Clinical implications

Triploidy is a life-limiting genetic aberration resulting from an extra haploid set of chromosomes of paternal (diandric triploidy) or maternal origin (digynic triploidy). Triploidy affects around 1%-2% of all conceptions. The majority of cases is miscarried at early developmental stages. In consequence of genomic imprinting, parental origin affects the phenotype of triploid pregnancies as well as the prevalence and spectrum of related maternal complications. Distinctive ultrasound features of both triploid phenotypes as well as characteristic patterns of biochemical markers may be useful in diagnosis. Molecular confirmation of the parental origin allows to predict the risk of complications, such as gestational trophoblastic neoplasia, hyperthyroidism, hypertension, or preeclampsia associated with the paternal origin of triploidy. Diagnosis of partial hydatidiform mole associated with diandric triploidy is challenging especially in the first trimester pregnancy loss due to the limitations of both histopathology and ultrasound. We present important clinical aspects of triploid pregnancies and indicate unresolved issues demanding further studies.

Onchocerciasis (river blindness) - more than a century of research and control

This review summarises more than a century of research on onchocerciasis, also known as river blindness, and its control. River blindness is an infection caused by the tissue filaria Onchocerca volvulus affecting the skin, subcutaneous tissue and eyes and leading to blindness in a minority of infected persons. The parasite is transmitted by its intermediate hosts Simulium spp. which breed in rivers. Featured are history and milestones in onchocerciasis research and control, state-of-the-art data on the parasite, its endobacteria Wolbachia, on the vectors, previous and current prevalence of the infection, its diagnostics, the interaction between the parasite and its host, immune responses and the pathology of onchocerciasis. Detailed information is documented on the time course of control programmes in the afflicted countries in Africa and the Americas, a long road from previous programmes to current successes in control of the transmission of this infectious disease. By development, adjustment and optimization of the control measures, transmission by the vector has been interrupted in foci of countries in the Americas, in Uganda, in Sudan and elsewhere, followed by onchocerciasis eliminations. The current state and future perspectives for control, elimination and eradication within the next 20-30 years are described and discussed. This review contributes to a deeper comprehension of this disease by a tissue-dwelling filaria and it will be helpful in efforts to control and eliminate other filarial infections.

Wolbachia: endosymbiont of onchocercid nematodes and their vectors

Background

Wolbachia is an obligate intracellular maternally transmitted, gram-negative bacterium which forms a spectrum of endosymbiotic relationships from parasitism to obligatory mutualism in a wide range of arthropods and onchocercid nematodes, respectively. In arthropods Wolbachia produces reproductive manipulations such as male killing, feminization, parthenogenesis and cytoplasmic incompatibility for its propagation and provides an additional fitness benefit for the host to protect against pathogens, whilst in onchocercid nematodes, apart from the mutual metabolic dependence, this bacterium is involved in moulting, embryogenesis, growth and survival of the host.

Methods

This review details the molecular data of Wolbachia and its effect on host biology, immunity, ecology and evolution, reproduction, endosymbiont-based treatment and control strategies exploited for filariasis. Relevant peer-reviewed scientic papers available in various authenticated scientific data bases were considered while writing the review.

Conclusions

The information presented provides an overview on Wolbachia biology and its use in the control and/or treatment of vectors, onchocercid nematodes and viral diseases of medical and veterinary importance. This offers the development of new approaches for the control of a variety of vector-borne diseases.

Graphic Abstract

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Genome editing as control tool for filarial infections

Human filarial infections are vector-borne nematode infections, which include lymphatic filariasis, onchocerciasis, loiasis, and mansonella filariasis. With a high prevalence in developing countries, filarial infections are responsible for some of the most debilitating morbidities and a vicious cycle of poverty and disease. Global initiatives set to eradicate these infections include community mass treatments, vector control, provision of care for morbidity, and search for vaccines. However, there are growing challenges associated with mass treatments, vector control, and antifilarial vaccine development. With the emergence of genome editing tools and successful applications in other infectious diseases, the integration of genetic editing techniques in future control strategies for filarial infections would offer the best option for eliminating filarial infections. In this review, we briefly discuss the mechanisms of the three main genetic editing techniques and explore the potential applications of these powerful tools to control filarial infections.

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.

An Integrated Approach to Identify New Anti-Filarial Leads to Treat River Blindness, a Neglected Tropical Disease

Filarial worms cause multiple debilitating diseases in millions of people worldwide, including river blindness. Currently available drugs reduce transmission by killing larvae (microfilariae), but there are no effective cures targeting the adult parasites (macrofilaricides) which survive and reproduce in the host for very long periods. To identify effective macrofilaricides, we carried out phenotypic screening of a library of 2121 approved drugs for clinical use against adult Brugia pahangi and prioritized the hits for further studies by integrating those results with a computational prioritization of drugs and associated targets. This resulted in the identification of 18 hits with anti-macrofilaricidal activity, of which two classes, azoles and aspartic protease inhibitors, were further expanded upon. Follow up screening against Onchocerca spp. (adult Onchocerca ochengi and pre-adult O. volvulus) confirmed activity for 13 drugs (the majority having IC50 < 10 μM), and a counter screen of a subset against L. loa microfilariae showed the potential to identify selective drugs that prevent adverse events when co-infected individuals are treated. Stage specific activity was also observed. Many of these drugs are amenable to structural optimization, and also have known canonical targets, making them promising candidates for further optimization that can lead to identifying and characterizing novel anti-macrofilarial drugs.

Diminutive, degraded but dissimilar: Wolbachia genomes from filarial nematodes do not conform to a single paradigm

Wolbachia are alpha-proteobacteria symbionts infecting a large range of arthropod species and two different families of nematodes. Interestingly, these endosymbionts are able to induce diverse phenotypes in their hosts: they are reproductive parasites within many arthropods, nutritional mutualists within some insects and obligate mutualists within their filarial nematode hosts. Defining Wolbachia 'species' is controversial and so they are commonly classified into 17 different phylogenetic lineages, termed supergroups, named A-F, H-Q and S. However, available genomic data remain limited and not representative of the full Wolbachia diversity; indeed, of the 24 complete genomes and 55 draft genomes of Wolbachia available to date, 84 % belong to supergroups A and B, exclusively composed of Wolbachia from arthropods. For the current study, we took advantage of a recently developed DNA-enrichment method to produce four complete genomes and two draft genomes of Wolbachia from filarial nematodes. Two complete genomes, wCtub and wDcau, are the smallest Wolbachia genomes sequenced to date (863 988 bp and 863 427 bp, respectively), as well as the first genomes representing supergroup J. These genomes confirm the validity of this supergroup, a controversial clade due to weaknesses of the multilocus sequence typing approach. We also produced the first draft Wolbachia genome from a supergroup F filarial nematode representative (wMhie), two genomes from supergroup D (wLsig and wLbra) and the complete genome of wDimm from supergroup C. Our new data confirm the paradigm of smaller Wolbachia genomes from filarial nematodes containing low levels of transposable elements and the absence of intact bacteriophage sequences, unlike many Wolbachia from arthropods, where both are more abundant. However, we observe differences among the Wolbachia genomes from filarial nematodes: no global co-evolutionary pattern, strong synteny between supergroup C and supergroup J Wolbachia, and more transposable elements observed in supergroup D Wolbachia compared to the other supergroups. Metabolic pathway analysis indicates several highly conserved pathways (haem and nucleotide biosynthesis, for example) as opposed to more variable pathways, such as vitamin B biosynthesis, which might be specific to certain host-symbiont associations. Overall, there appears to be no single Wolbachia-filarial nematode pattern of co-evolution or symbiotic relationship.

Population genomic evidence that human and animal infections in Africa come from the same populations of Dracunculus medinensis

BACKGROUND

Guinea worm-Dracunculus medinensis-was historically one of the major parasites of humans and has been known since antiquity. Now, Guinea worm is on the brink of eradication, as efforts to interrupt transmission have reduced the annual burden of disease from millions of infections per year in the 1980s to only 54 human cases reported globally in 2019. Despite the enormous success of eradication efforts to date, one complication has arisen. Over the last few years, hundreds of dogs have been found infected with this previously apparently anthroponotic parasite, almost all in Chad. Moreover, the relative numbers of infections in humans and dogs suggests that dogs are currently the principal reservoir on infection and key to maintaining transmission in that country.

PRINCIPAL FINDINGS

In an effort to shed light on this peculiar epidemiology of Guinea worm in Chad, we have sequenced and compared the genomes of worms from dog, human and other animal infections. Confirming previous work with other molecular markers, we show that all of these worms are D. medinensis, and that the same population of worms are causing both infections, can confirm the suspected transmission between host species and detect signs of a population bottleneck due to the eradication efforts. The diversity of worms in Chad appears to exclude the possibility that there were no, or very few, worms present in the country during a 10-year absence of reported cases.

CONCLUSIONS

This work reinforces the importance of adequate surveillance of both human and dog populations in the Guinea worm eradication campaign and suggests that control programs aiming to interrupt disease transmission should stay aware of the possible emergence of unusual epidemiology as pathogens approach elimination.

Speciation and adaptive evolution reshape antioxidant enzymatic system diversity across the phylum Nematoda

Background

Nematodes have evolved to survive in diverse ecological niches and can be a serious burden on agricultural economy, veterinary medicine, and public health. Antioxidant enzymes in parasitic nematodes play a critical role in defending against host oxidative stress. However, the features of the evolution of antioxidant enzymes in the phylum Nematoda remain elusive.

Results

Here, we systematically investigated the evolution and gene expression of antioxidant enzymes in the genomes of 59 nematodes and transcriptomes of 20 nematodes. Catalase has been independently lost in several orders, suggesting that it is unnecessary for some nematodes. Unlike in mammals, phospholipid hydroperoxide glutathione peroxidase is widely distributed in nematodes, among which it has evolved independently. We found that superoxide dismutase (SOD) has been present throughout nematode evolutionary process, and the extracellular isoform (SOD3) is diverged from the corresponding enzyme in mammals and has undergone duplication and differentiation in several nematodes. Moreover, the evolution of intracellular and extracellular SOD isoforms in filaria strongly indicates that extracellular SOD3 originated from intracellular SOD1 and underwent rapid evolution to form the diversity of extracellular SOD3. We identify a novel putative metal-independent extracellular SOD presenting independently in Steinernema and Strongyloididae lineage that featured a high expression level in Strongyloides larvae. Sequence divergence of SOD3 between parasitic nematodes and their closest free-living nematode, the specifically high expression in the parasitic female stage, and presence in excretory-secretory proteome of Strongyloides suggest that SOD3 may be related with parasitism.

Conclusions

This study advances our understanding of the complex evolution of antioxidant enzymes across Nematoda and provides targets for controlling parasitic nematode diseases.

The Functional Parasitic Worm Secretome: Mapping the Place of Onchocerca volvulus Excretory Secretory Products

Nematodes constitute a very successful phylum, especially in terms of parasitism. Inside their mammalian hosts, parasitic nematodes mainly dwell in the digestive tract (geohelminths) or in the vascular system (filariae). One of their main characteristics is their long sojourn inside the body where they are accessible to the immune system. Several strategies are used by parasites in order to counteract the immune attacks. One of them is the expression of molecules interfering with the function of the immune system. Excretory-secretory products (ESPs) pertain to this category. This is, however, not their only biological function, as they seem also involved in other mechanisms such as pathogenicity or parasitic cycle (molting, for example). We will mainly focus on filariae ESPs with an emphasis on data available regarding Onchocerca volvulus, but we will also refer to a few relevant/illustrative examples related to other worm categories when necessary (geohelminth nematodes, trematodes or cestodes). We first present Onchocerca volvulus, mainly focusing on the aspects of this organism that seem relevant when it comes to ESPs: life cycle, manifestations of the sickness, immunosuppression, diagnosis and treatment. We then elaborate on the function and use of ESPs in these aspects.

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.

Filarial parasitic nematodes are vector-borne parasites that infect humans and animals. Brugia malayi and Dirofilaria immitis are transmitted by mosquitoes and cause human lymphatic filariasis and canine heartworm disease, respectively. Recent years have seen a dramatic increase in genomic and transcriptomic data sets and the concomitant increase in innovative strategies for drug target identification, validation, and screening. However, while the completeness of genome assemblies of filarial parasitic nematodes has seen steady improvements, the reliability of gene models has not kept pace, hindering cloning efforts. Long-read RNA sequencing technologies are uniquely able to improve gene models, but have not been widely used for the causative agents of neglected tropical diseases. Here, we report the improvement of gene models in both B. malayi and D. immitis by long-read RNA sequencing. We identified novel operons, deprecated false positive operons, identified dozens of novel genes, and described the parameters of polyadenylation. We also focused on putative anthelmintic targets, identifying novel isoforms and correcting gene models. These data substantially increase the trustworthiness of gene models in these two species and demonstrate how long-read sequencing approaches should be prioritized in the continued improvement of genome assemblies and their gene annotations.

Genomic and transcriptomic variation defines the chromosome-scale assembly of Haemonchus contortus, a model gastrointestinal worm

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

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.

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.

Prediction pipeline for discovery of regulatory motifs associated with Brugia malayi molting

Filarial nematodes can cause debilitating diseases in humans. They have complicated life cycles involving an insect vector and mammalian hosts, and they go through a number of developmental molts. While whole genome sequences of parasitic worms are now available, very little is known about transcription factor (TF) binding sites and their cognate transcription factors that play a role in regulating development. To address this gap, we developed a novel motif prediction pipeline, Emotif Alpha, that integrates ten different motif discovery algorithms, multiple statistical tests, and a comparative analysis of conserved elements between the filarial worms Brugia malayi and Onchocerca volvulus, and the free-living nematode Caenorhabditis elegans. We identified stage-specific TF binding motifs in B. malayi, with a particular focus on those potentially involved in the L3-L4 molt, a stage important for the establishment of infection in the mammalian host. Using an in vitro molting system, we tested and validated three of these motifs demonstrating the accuracy of the motif prediction pipeline.

Diseases caused by parasitic worms such as the filariae are among the leading causes of morbidity in the developing world. Very little is known about how development is regulated in these vector-transmitted parasites. We have developed a computational method to identify motifs that correspond to transcription factor binding sites in the genome of the parasitic worm, Brugia malayi, one of the causative agents of lymphatic filariasis. Using this approach, we were able to predict stage-specific transcription factor binding sites involved in a stage of the molting process important for the establishment of the infection. We validated the role of these motifs using an in vitro molting system.

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.

Sex chromosome evolution in parasitic nematodes of humans

Sex determination mechanisms often differ even between related species yet the evolution of sex chromosomes remains poorly understood in all but a few model organisms. Some nematodes such as Caenorhabditis elegans have an XO sex determination system while others, such as the filarial parasite Brugia malayi, have an XY mechanism. We present a complete B. malayi genome assembly and define Nigon elements shared with C. elegans, which we then map to the genomes of other filarial species and more distantly related nematodes. We find a remarkable plasticity in sex chromosome evolution with several distinct cases of neo-X and neo-Y formation, X-added regions, and conversion of autosomes to sex chromosomes from which we propose a model of chromosome evolution across different nematode clades. The phylum Nematoda offers a new and innovative system for gaining a deeper understanding of sex chromosome evolution.

The insufficiency of circulating miRNA and DNA as diagnostic tools or as biomarkers of treatment efficacy for Onchocerca volvulus

Skin snip evaluation for onchocerciasis has insufficient sensitivity when skin microfilarial (mf) densities are low, such as following ivermectin treatment. Mf density is suitable for assessing microfilaricidal efficacy but only serves as an indirect indicator of macrofilaricidal activity. We assessed circulating nucleic acids from Onchocerca volvulus as an alternative to skin snips. We screened a plasma sample set of infected individuals followed up at four, 12 and 21 months after microfilaricidal (ivermectin, n = four), macrofilaricidal (doxycycline, n = nine), or combination treatment (n = five). Two parasite-derived miRNAs, cel-miR-71-5p and bma-lin-4, and O-150 repeat DNA were assessed. Highly abundant DNA repeat families identified in the O. volvulus genome were also evaluated. miRNAs were detected in two of 72 plasma samples (2.8%) and two of 47 samples (4.3%) with microfilaridermia using RT-qPCR. O-150 DNA was detected in eight (44.4%) baseline samples by qPCR and the number of positives declined post-treatment. One doxycycline-treated individual remained O-150 positive. However, only 11 (23.4%) samples with microfilaridermia were qPCR-positive. Analysis by qPCR showed novel DNA repeat families were comparatively less abundant than the O-150 repeat. Circulating parasite-derived nucleic acids are therefore insufficient as diagnostic tools or as biomarkers of treatment efficacy for O. volvulus.

Genomic Epidemiology in Filarial Nematodes: Transforming the Basis for Elimination Program Decisions

Onchocerciasis and lymphatic filariasis are targeted for elimination, primarily using mass drug administration at the country and community levels. Elimination of transmission is the onchocerciasis target and global elimination as a public health problem is the end point for lymphatic filariasis. Where program duration, treatment coverage, and compliance are sufficiently high, elimination is achievable for both parasites within defined geographic areas. However, transmission has re-emerged after apparent elimination in some areas, and in others has continued despite years of mass drug treatment. A critical question is whether this re-emergence and/or persistence of transmission is due to persistence of local parasites-i.e., the result of insufficient duration or drug coverage, poor parasite response to the drugs, or inadequate methods of assessment and/or criteria for determining when to stop treatment-or due to re-introduction of parasites via human or vector movement from another endemic area. We review recent genetics-based research exploring these questions in Onchocerca volvulus, the filarial nematode that causes onchocerciasis, and Wuchereria bancrofti, the major pathogen for lymphatic filariasis. We focus in particular on the combination of genomic epidemiology and genome-wide associations to delineate transmission zones and distinguish between local and introduced parasites as the source of resurgence or continuing transmission, and to identify genetic markers associated with parasite response to chemotherapy. Our ultimate goal is to assist elimination efforts by developing easy-to-use tools that incorporate genetic information about transmission and drug response for more effective mass drug distribution, surveillance strategies, and decisions on when to stop interventions to improve sustainability of elimination.