The genome of Onchocerca volvulus, agent of river blindness

Identification of small molecule enzyme inhibitors as broad-spectrum anthelmintics

Targeting chokepoint enzymes in metabolic pathways has led to new drugs for cancers, autoimmune disorders and infectious diseases. This is also a cornerstone approach for discovery and development of anthelmintics against nematode and flatworm parasites. Here, we performed omics-driven knowledge-based identification of chokepoint enzymes as anthelmintic targets. We prioritized 10 of 186 phylogenetically conserved chokepoint enzymes and undertook a target class repurposing approach to test and identify new small molecules with broad spectrum anthelmintic activity. First, we identified and tested 94 commercially available compounds using an in vitro phenotypic assay, and discovered 11 hits that inhibited nematode motility. Based on these findings, we performed chemogenomic screening and tested 32 additional compounds, identifying 6 more active hits. Overall, 6 intestinal (single-species), 5 potential pan-intestinal (whipworm and hookworm) and 6 pan-Phylum Nematoda (intestinal and filarial species) small molecule inhibitors were identified, including multiple azoles, Tadalafil and Torin-1. The active hit compounds targeted three different target classes in humans, which are involved in various pathways, including carbohydrate, amino acid and nucleotide metabolism. Last, using representative inhibitors from each target class, we demonstrated in vivo efficacy characterized by negative effects on parasite fecundity in hamsters infected with hookworms.

Genome expansion of an obligate parthenogenesis-associated Wolbachia poses an exception to the symbiont reduction model

Background

Theory predicts that dependency within host-endosymbiont interactions results in endosymbiont genome size reduction. Unexpectedly, the largest Wolbachia genome was found in the obligate, parthenogenesis-associated wFol. In this study, we investigate possible processes underlying this genome expansion by comparing a re-annotated wFol genome to other Wolbachia genomes. In addition, we also search for candidate genes related to parthenogenesis induction (PI).

Results

Within wFol, we found five phage WO regions representing 25.4% of the complete genome, few pseudogenized genes, and an expansion of DNA-repair genes in comparison to other Wolbachia. These signs of genome conservation were mirrored in the wFol host, the springtail F. candida, which also had an expanded DNA-repair gene family and many horizontally transferred genes. Across all Wolbachia genomes, there was a strong correlation between gene numbers of Wolbachia strains and their hosts. In order to identify genes with a potential link to PI, we assembled the genome of an additional PI strain, wLcla. Comparisons between four PI Wolbachia, including wFol and wLcla, and fourteen non-PI Wolbachia yielded a small set of potential candidate genes for further investigation.

Conclusions

The strong similarities in genome content of wFol and its host, as well as the correlation between host and Wolbachia gene numbers suggest that there may be some form of convergent evolution between endosymbiont and host genomes. If such convergent evolution would be strong enough to overcome the evolutionary forces causing genome reduction, it would enable expanded genomes within long-term obligate endosymbionts.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5492-9) contains supplementary material, which is available to authorized users.

Dirofilaria and Wolbachia in mosquitoes (Diptera: Culicidae) in central European Russia and on the Black Sea coast

Dirofilariasis is endemic in Russia, as well as in many other European countries. The aim of this study was to assess the ability of mosquitoes to transfer Dirofilaria immitis and Dirofilaria repens in regions with temperate and subtropical climates. The possible impact of the symbiotic bacterium Wolbachia on Dirofilaria transmission was also investigated. 5333 female mosquitoes were collected at 11 points in central European Russia and on the Black Sea coast during the period 2013–2017. Out of 20 mosquito species examined, 14 were infected with D. repens and 13 with D. immitis. Both species of Dirofilaria were found in different climatic regions. The total Dirofilaria spp. estimated infection rate (EIR) in the central part of Russia varied from 3.1% to 3.7% and, in the southern region, from 1.1% to 3.0%. The highest estimated infection rate was found in Anopheles messeae, the lowest in Culex pipiens. The greatest epidemiological danger was represented by Aedes aegypti, Ae. geniculatus, An. messeae and Ae. communis. Six out of 20 mosquito species were infected with Wolbachia. Pools of Aedes albopictus, Cx. pipiens and Coquillettidia richiardii were simultaneously infected with Dirofilaria and Wolbachia. After checking mosquitoes individually, it was found that there was no development of Dirofilaria to the infective larval stage in specimens infected with Wolbachia. Twenty-two Dirofilaria-infective pools were Wolbachia-free and only two mosquito pools were Wolbachia-infected. The potential for transmission of Dirofilaria in mosquito species naturally uninfected with the symbiotic bacterium Wolbachia is higher than in species infected with the bacterium.

Comparative genomics of the major parasitic worms

Parasitic nematodes (roundworms) and platyhelminths (flatworms) cause debilitating chronic infections of humans and animals, decimate crop production and are a major impediment to socioeconomic development. Here we report a broad comparative study of 81 genomes of parasitic and non-parasitic worms. We have identified gene family births and hundreds of expanded gene families at key nodes in the phylogeny that are relevant to parasitism. Examples include gene families that modulate host immune responses, enable parasite migration though host tissues or allow the parasite to feed. We reveal extensive lineage-specific differences in core metabolism and protein families historically targeted for drug development. From an in silico screen, we have identified and prioritized new potential drug targets and compounds for testing. This comparative genomics resource provides a much-needed boost for the research community to understand and combat parasitic worms.

Highlights on the Application of Genomics and Bioinformatics in the Fight Against Infectious Diseases: Challenges and Opportunities in Africa

Genomics and bioinformatics are increasingly contributing to our understanding of infectious diseases caused by bacterial pathogens such as Mycobacterium tuberculosis and parasites such as Plasmodium falciparum. This ranges from investigations of disease outbreaks and pathogenesis, host and pathogen genomic variation, and host immune evasion mechanisms to identification of potential diagnostic markers and vaccine targets. High throughput genomics data generated from pathogens and animal models can be combined with host genomics and patients’ health records to give advice on treatment options as well as potential drug and vaccine interactions. However, despite accounting for the highest burden of infectious diseases, Africa has the lowest research output on infectious disease genomics. Here we review the contributions of genomics and bioinformatics to the management of infectious diseases of serious public health concern in Africa including tuberculosis (TB), dengue fever, malaria and filariasis. Furthermore, we discuss how genomics and bioinformatics can be applied to identify drug and vaccine targets. We conclude by identifying challenges to genomics research in Africa and highlighting how these can be overcome where possible.

Novel Findings of Anti-Filarial Drug Target and Structure-Based Virtual Screening for Drug Discovery

Lymphatic filariasis and onchocerciasis caused by filarial nematodes are important diseases leading to considerable morbidity throughout tropical countries. Diethylcarbamazine (DEC), albendazole (ALB), and ivermectin (IVM) used in massive drug administration are not highly effective in killing the long-lived adult worms, and there is demand for the development of novel macrofilaricidal drugs affecting new molecular targets. A Ca2+ binding protein, calumenin, was identified as a novel and nematode-specific drug target for filariasis, due to its involvement in fertility and cuticle development in nematodes. As sterilizing and killing effects of the adult worms are considered to be ideal profiles of new drugs, calumenin could be an eligible drug target. Indeed, the Caenorhabditis elegans mutant model of calumenin exhibited enhanced drug acceptability to both microfilaricidal drugs (ALB and IVM) even at the adult stage, proving the roles of the nematode cuticle in efficient drug entry. Molecular modeling revealed that structural features of calumenin were only conserved among nematodes (C. elegans, Brugia malayi, and Onchocerca volvulus). Structural conservation and the specificity of nematode calumenins enabled the development of drugs with good target selectivity between parasites and human hosts. Structure-based virtual screening resulted in the discovery of itraconazole (ITC), an inhibitor of sterol biosynthesis, as a nematode calumenin-targeting ligand. The inhibitory potential of ITC was tested using a nematode mutant model of calumenin.

Prediction and validation of the structural features of Ov58GPCR, an immunogenic determinant of Onchocerca volvulus

Onchocerciasis is a severely debilitating yet neglected tropical disease (NTD) that creates social stigma, generates and perpetuates poverty, and leads ultimately in some cases to irreversible unilateral or bilateral blindness if untreated. Consequently, the disease is a major impediment to socioeconomic development. Many control programs have been launched for the disease with moderate successes achieved. This mitigated hit is partially due to the lingering need for reliable, non-invasive and easily applicable tools for mapping endemic regions and post-elimination surveillance. In this work, bioinformatics analyses combined with immunological assays were applied in a bid to develop potential tools for diagnosis and assessing the success of drug treatment programs. We report that (i) the O. volvulus antigen, Ov58GPCR is a G-protein coupled receptor (GPCR) conserved in related nematodes, (ii) synthetic peptides predicted to be in the extracellular domain (ECD) of Ov58GPCR are indeed immunogenic epitopes in actively-infected individuals, (iii) synthetic peptide cocktails discriminate between actively-infected individuals, treated individuals and healthy African controls, (iv) polyclonal antibodies against one of the peptides or against the bacterially-expressed ECD reacted specifically with the native antigen of O. volvulus total and surface extracts, (v) Ov58GPCR is transcribed in both larvae and adult parasite stages, (vi) IgG and IgE responses to the recombinant ECD decline with ivermectin treatment. All these findings suggest that the extracellular domain and synthetic peptides of Ov58GPCR, as well as the specific immune response generated could be harnessed in the context of disease diagnosis and surveillance.

Onchodermatitis: Where Are We Now?

Onchocerciasis causes debilitating pruritus and rashes as well as visual impairment and blindness. Prior to control measures, eye disease was particularly prominent in savanna areas of sub-Saharan Africa whilst skin disease was more common across rainforest regions of tropical Africa. Mass drug distribution with ivermectin is changing the global scene of onchocerciasis. There has been successful progressive elimination in Central and Southern American countries and the World Health Organization has set a target for elimination in Africa of 2025. This literature review was conducted to examine progress regarding onchocercal skin disease. PubMed searches were performed using keywords 'onchocerciasis', 'onchodermatitis' and 'onchocercal skin disease' over the past eight years. Articles in English, or with an English abstract, were assessed for relevance, including any pertinent references within the articles. Recent progress in awareness of, understanding and treatment of onchocercal skin disease is reviewed with particular emphasis on publications within the past five years. The global burden of onchodermatitis is progressively reducing and is no longer seen in children in many formerly endemic foci.

Cysteine proteases during larval migration and development of helminths in their final host

Neglected tropical diseases caused by metazoan parasites are major public health concerns, and therefore, new methods for their control and elimination are needed. Research over the last 25 years has revealed the vital contribution of cysteine proteases to invasion of and migration by (larval) helminth parasites through host tissues, in addition to their roles in embryogenesis, molting, egg hatching, and yolk degradation. Their central function to maintaining parasite survival in the host has made them prime intervention targets for novel drugs and vaccines. This review focuses on those helminth cysteine proteases that have been functionally characterized during the varied early stages of development in the human host and embryogenesis.

A Role for Epitope Networking in Immunomodulation by Helminths

Helminth infections, by nematodes, trematodes, or cestodes, can lead to the modulation of host immune responses. This allows long-duration parasite infections and also impacts responses to co-infections. Surface, secreted, excreted, and shed proteins are thought to play a major role in modulation. A commonly reported feature of such immune modulation is the role of T regulatory (Treg) cells and IL-10. Efforts to identify helminth proteins, which cause immunomodulation, have identified candidates but not provided clarity as to a uniform mechanism driving modulation. In this study, we applied a bioinformatics systems approach, allowing us to analyze predicted T-cell epitopes of 17 helminth species and the responses to their surface proteins. In addition to major histocompatibility complex (MHC) binding, we analyzed amino acid motifs that would be recognized by T-cell receptors [T-cell-exposed motifs (TCEMs)]. All the helminth species examined have, within their surface proteins, peptides, which combine very common TCEMs with predicted high affinity binding to many human MHC alleles. This combination of features would result in large cognate T cell and a high probability of eliciting Treg responses. The TCEMs, which determine recognition by responding T-cell clones, are shared to a high degree between helminth species and with Plasmodium falciparum and Mycobacterium tuberculosis, both common co-infecting organisms. The implication of our observations is not only that Treg cells play a significant role in helminth-induced immune modulation but also that the epitope specificities of Treg responses are shared across species and genera of helminth. Hence, the immune response to a given helminth cannot be considered in isolation but rather forms part of an epitope ecosystem, or microenvironment, in which potentially immunosuppressive peptides in the helminth network via their common T-cell receptor recognition signals with T-cell epitopes in self proteins, microbiome, other helminths, and taxonomically unrelated pathogens. Such a systems approach provides a high-level view of the antigen-immune system signaling dynamics that may bias a host's immune response to helminth infections toward immune modulation. It may indicate how helminths have evolved to select for peptides that favor long-term parasite host coexistence.

Draft genome sequence of the Wolbachia endosymbiont of Wuchereria bancrofti wWb

The draft genome assembly of the Wolbachia endosymbiont of Wuchereria bancrofti (wWb) consists of 1060 850 bp in 100 contigs and contains 961 ORFs, with a single copy of the 5S rRNA, 16S rRNA and 23S rRNA and each of the 34 tRNA genes. Phylogenetic core genome analyses show wWb to cluster with other strains in supergroup D of the Wolbachia phylogeny, while being most closely related to the Wolbachia endosymbiont of Brugia malayi strain TRS (wBm). The wWb and wBm genomes share 779 orthologous clusters with wWb having 101 unclustered genes and wBm having 23 unclustered genes. The higher number of unclustered genes in the wWb genome likely reflects the fragmentation of the draft genome.

Single-Molecule Sequencing Reveals the Chromosome-Scale Genomic Architecture of the Nematode Model Organism Pristionchus pacificus

The nematode Pristionchus pacificus is an established model for integrative evolutionary biology and comparative studies with Caenorhabditis elegans. While an existing genome draft facilitated the identification of several genes controlling various developmental processes, its high degree of fragmentation complicated virtually all genomic analyses. Here, we present a de novo genome assembly from single-molecule, long-read sequencing data consisting of 135 P. pacificus contigs. When combined with a genetic linkage map, 99% of the assembly could be ordered and oriented into six chromosomes. This allowed us to robustly characterize chromosomal patterns of gene density, repeat content, nucleotide diversity, linkage disequilibrium, and macrosynteny in P. pacificus. Despite widespread conservation of synteny between P. pacificus and C. elegans, we identified one major translocation from an autosome to the sex chromosome in the lineage leading to C. elegans. This highlights the potential of the chromosome-scale assembly for future genomic studies of P. pacificus.

The zinc transporter ZIPT-7.1 regulates sperm activation in nematodes

Sperm activation is a fascinating example of cell differentiation, in which immotile spermatids undergo a rapid and dramatic transition to become mature, motile sperm. Because the sperm nucleus is transcriptionally silent, this transition does not involve transcriptional changes. Although Caenorhabditis elegans is a leading model for studies of sperm activation, the mechanisms by which signaling pathways induce this transformation remain poorly characterized. Here we show that a conserved transmembrane zinc transporter, ZIPT-7.1, regulates the induction of sperm activation in Caenorhabditis nematodes. The zipt-7.1 mutant hermaphrodites cannot self-fertilize, and males reproduce poorly, because mutant spermatids are defective in responding to activating signals. The zipt-7.1 gene is expressed in the germ line and functions in germ cells to promote sperm activation. When expressed in mammalian cells, ZIPT-7.1 mediates zinc transport with high specificity and is predominantly located on internal membranes. Finally, genetic epistasis places zipt-7.1 at the end of the spe-8 sperm activation pathway, and ZIPT-7.1 binds SPE-4, a presenilin that regulates sperm activation. Based on these results, we propose a new model for sperm activation. In spermatids, inactive ZIPT-7.1 is localized to the membranous organelles, which contain higher levels of zinc than the cytoplasm. When sperm activation is triggered, ZIPT-7.1 activity increases, releasing zinc from internal stores. The resulting increase in cytoplasmic zinc promotes the phenotypic changes characteristic of activation. Thus, zinc signaling is a key step in the signal transduction process that mediates sperm activation, and we have identified a zinc transporter that is central to this activation process.

Sperm are specialized cells with transcriptionally silent DNA that has been packaged for delivery into the egg. In their final step of development, immature sperm undergo a rapid transition from nonmotile cells to mature, motile sperm capable of fertilization. The signals that trigger this change are not clearly understood. By identifying mutants in the roundworm Caenorhabditis elegans that are defective in sperm activation, we discovered a conserved transmembrane protein, ZIPT-7.1, that transports zinc and promotes sperm activation in both sexes. ZIPT-7.1 is expressed in the germ line and functions there to control sperm activation. When expressed ectopically in mammalian cells, the protein specifically transports zinc across membranes and localizes primarily to membranes within the cell. Previous genetic studies had identified two pathways that mediate sperm activation in C. elegans, and our results suggest that zipt-7.1 acts at the end of one of these two, the spe-8 pathway. We propose that when this pathway triggers sperm activation, it acts through ZIPT-7.1, which mediates the release of zinc from internal stores in the immature sperm. This released zinc functions as a second messenger to promote the differentiation of mature, motile sperm.

G-protein-coupled receptor genes of Dirofilaria immitis

The diversity and uniqueness of nematode heterotrimeric G-protein-coupled receptors (GPCRs) provides impetus for identifying ligands that can be used as therapeutics for treating diseases caused by parasitic nematode infections. In human medicine, GPCRs have represented the largest group of 'drugable' targets exploited in the market today. In the filarial nematode Dirofilaria immitis, which causes heartworm disease, the macrocyclic lactones (ML) have been used as the sole preventatives for more than 25 years and now there is confirmed ML resistance in this parasite. A novel anthelmintic emodepside, with antifilarial activity, can act on a GPCR. In view of the ML resistance, there is an urgent need to identify new drug targets and GPCRs of D. immitis may be promising receptors. Knowledge of polymorphism within the GPCR superfamily is of interest. A total of 127 GPCR genes have been identified, so far, in the genome of D. immitis. Whole genome sequencing data from four ML susceptible and four ML loss of efficacy populations was used to identify 393 polymorphic loci in 35 D. immitis GPCR genes. Out of 57 SNPs in exonic regions, 36 of them caused a change in an amino acid, out of which 2 changed the predicted secondary structure of the protein. Knowledge about GPCR genes and their polymorphism is valuable information for drug design processes. Further studies need to be carried out to more fully understand the implications of each of the SNPs identified by this study.

Inferring synteny between genome assemblies: a systematic evaluation

BACKGROUND

Genome assemblies across all domains of life are being produced routinely. Initial analysis of a new genome usually includes annotation and comparative genomics. Synteny provides a framework in which conservation of homologous genes and gene order is identified between genomes of different species. The availability of human and mouse genomes paved the way for algorithm development in large-scale synteny mapping, which eventually became an integral part of comparative genomics. Synteny analysis is regularly performed on assembled sequences that are fragmented, neglecting the fact that most methods were developed using complete genomes. It is unknown to what extent draft assemblies lead to errors in such analysis.

RESULTS

We fragmented genome assemblies of model nematodes to various extents and conducted synteny identification and downstream analysis. We first show that synteny between species can be underestimated up to 40% and find disagreements between popular tools that infer synteny blocks. This inconsistency and further demonstration of erroneous gene ontology enrichment tests raise questions about the robustness of previous synteny analysis when gold standard genome sequences remain limited. In addition, assembly scaffolding using a reference guided approach with a closely related species may result in chimeric scaffolds with inflated assembly metrics if a true evolutionary relationship was overlooked. Annotation quality, however, has minimal effect on synteny if the assembled genome is highly contiguous.

CONCLUSIONS

Our results show that a minimum N50 of 1 Mb is required for robust downstream synteny analysis, which emphasizes the importance of gold standard genomes to the science community, and should be achieved given the current progress in sequencing technology.

Mining Filarial Genomes for Diagnostic and Therapeutic Targets

Filarial infections of humans cause some of the most important neglected tropical diseases. The global efforts for eliminating filarial infections by mass drug administration programs may require additional tools (safe macrofilaricidal drugs, vaccines, and diagnostic biomarkers). The accurate and sensitive detection of viable parasites is essential for diagnosis and for surveillance programs. Current community-wide treatment modalities do not kill the adult filarial worms effectively; hence, there is a need to identify and develop safe macrofilaricidal drugs. High-throughput sequencing, mass spectroscopy methods and advances in computational biology have greatly accelerated the discovery process. Here, we describe post-genomic developments toward the identification of diagnostic biomarkers and drug targets for the filarial infection of humans.

Genome Architecture and Evolution of a Unichromosomal Asexual Nematode

Asexual reproduction in animals, though rare, is the main or exclusive mode of reproduction in some long-lived lineages. The longevity of asexual clades may be correlated with the maintenance of heterozygosity by mechanisms that rearrange genomes and reduce recombination. Asexual species thus provide an opportunity to gain insight into the relationship between molecular changes, genome architecture, and cellular processes. Here we report the genome sequence of the parthenogenetic nematode Diploscapter pachys with only one chromosome pair. We show that this unichromosomal architecture is shared by a long-lived clade of asexual nematodes closely related to the genetic model organism Caenorhabditis elegans. Analysis of the genome assembly reveals that the unitary chromosome arose through fusion of six ancestral chromosomes, with extensive rearrangement among neighboring regions. Typical nematode telomeres and telomeric protection-encoding genes are lacking. Most regions show significant heterozygosity; homozygosity is largely concentrated to one region and attributed to gene conversion. Cell-biological and molecular evidence is consistent with the absence of key features of meiosis I, including synapsis and recombination. We propose that D. pachys preserves heterozygosity and produces diploid embryos without fertilization through a truncated meiosis. As a prelude to functional studies, we demonstrate that D. pachys is amenable to experimental manipulation by RNA interference.

Insights Into Onchocerca volvulus Population Biology Through Multilocus Immunophenotyping

We have developed a serologically based immunophenotyping approach to study Onchocerca volvulus (Ov) population diversity. Using genomic sequence data and polymerase chain reaction-based genotyping, we identified nonsynonymous single-nucleotide polymorphisms (SNPs) in the genes of 16 major immunogenic Ov proteins: Ov-CHI-1/Ov-CHI-2, Ov16, Ov-FAR-1, Ov-CPI-1, Ov-B20, Ov-ASP-1, Ov-TMY-1, OvSOD1, OvGST1, Ov-CAL-1, M3/M4, Ov-RAL-1, Ov-RAL-2, Ov-ALT-1, Ov-FBA-1, and Ov-B8. We assessed the immunoreactivity of onchocerciasis patient sera (n = 152) from the Americas, West Africa, Central Africa, and East Africa against peptides derived from 10 of these proteins containing SNPs. Statistically significant variation in immunoreactivity among the regions was seen in SNP-containing peptides derived from 8 of 10 proteins tested: OVOC1192(1-15), OVOC9988(28-42), OVOC9225(320-334), OVOC7453(22-36), OVOC11517(14-28), OVOC3177(283-297), OVOC7911(594-608), and OVOC12628(174-188). Our data show that differences in immunoreactivity to variant antigenic peptides may be used to characterize Ov populations, thereby elucidating features of Ov population biology previously inaccessible because of the limited availability of parasite material.

Onchocerca volvulus: The Road from Basic Biology to a Vaccine

Human onchocerciasis - commonly known as river blindness - is one of the most devastating yet neglected tropical diseases, leaving many millions in sub-Saharan Africa blind and/or with chronic disabilities. Attempts to eliminate onchocerciasis, primarily through the mass drug administration of ivermectin, remains challenging and has been heightened by the recent news that drug-resistant parasites are developing in some populations after years of drug treatment. Needed, and needed now, in the fight to eliminate onchocerciasis are new tools, such as preventive and therapeutic vaccines. This review summarizes the progress made to advance the onchocerciasis vaccine from the research laboratory into the clinic.

Validation of onchocerciasis biomarker N-acetyltyramine-O-glucuronide (NATOG)

The Neglected Tropical Disease onchocerciasis is a parasitic disease. Despite many control programmes by the World Health Organization (WHO), large communities in West and Central Africa are still affected. Besides logistic challenges during biannual mass drug administration, the lack of a robust, point-of-care diagnostic is limiting successful eradication of onchocerciasis. Towards the implementation of a non-invasive and point-of-care diagnostic, we have recently reported the discovery of the biomarker N-acetyltyramine-O-glucuronide (NATOG) in human urine samples using a metabolomics-mining approach. NATOG's biomarker value was enhanced during an investigation in a rodent model. Herein, we further detail the specificity of NATOG in active onchocerciasis infections as well as the co-infecting parasites Loa loa and Mansonella perstans. Our results measured by liquid chromatography coupled with mass spectrometry (LC-MS) reveal elevated NATOG values in mono- and co-infection samples only in the presence of the nematode Onchocerca volvulus. Metabolic pathway investigation of l-tyrosine/tyramine in all investigated nematodes uncovered an important link between the endosymbiotic bacterium Wolbachia and O. volvulus for the biosynthesis of NATOG. Based on these extended studies, we suggest NATOG as a biomarker for tracking active onchocerciasis infections and provide a threshold concentration value of NATOG for future diagnostic tool development.

Toward Universal Forward Genetics: Using a Draft Genome Sequence of the Nematode Oscheius tipulae To Identify Mutations Affecting Vulva Development

Mapping-by-sequencing has become a standard method to map and identify phenotype-causing mutations in model species. Here, we show that a fragmented draft assembly is sufficient to perform mapping-by-sequencing in nonmodel species. We generated a draft assembly and annotation of the genome of the free-living nematode Oscheius tipulae, a distant relative of the model Caenorhabditis elegans We used this draft to identify the likely causative mutations at the O. tipulae cov-3 locus, which affect vulval development. The cov-3 locus encodes the O. tipulae ortholog of C. elegans mig-13, and we further show that Cel-mig-13 mutants also have an unsuspected vulval-development phenotype. In a virtuous circle, we were able to use the linkage information collected during mutant mapping to improve the genome assembly. These results showcase the promise of genome-enabled forward genetics in nonmodel species.

Genome-wide analysis of ivermectin response by Onchocerca volvulus reveals that genetic drift and soft selective sweeps contribute to loss of drug sensitivity

Background

Treatment of onchocerciasis using mass ivermectin administration has reduced morbidity and transmission throughout Africa and Central/South America. Mass drug administration is likely to exert selection pressure on parasites, and phenotypic and genetic changes in several Onchocerca volvulus populations from Cameroon and Ghana—exposed to more than a decade of regular ivermectin treatment—have raised concern that sub-optimal responses to ivermectin's anti-fecundity effect are becoming more frequent and may spread.

Methodology/Principal findings

Pooled next generation sequencing (Pool-seq) was used to characterise genetic diversity within and between 108 adult female worms differing in ivermectin treatment history and response. Genome-wide analyses revealed genetic variation that significantly differentiated good responder (GR) and sub-optimal responder (SOR) parasites. These variants were not randomly distributed but clustered in ~31 quantitative trait loci (QTLs), with little overlap in putative QTL position and gene content between the two countries. Published candidate ivermectin SOR genes were largely absent in these regions; QTLs differentiating GR and SOR worms were enriched for genes in molecular pathways associated with neurotransmission, development, and stress responses. Finally, single worm genotyping demonstrated that geographic isolation and genetic change over time (in the presence of drug exposure) had a significantly greater role in shaping genetic diversity than the evolution of SOR.

Conclusions/Significance

This study is one of the first genome-wide association analyses in a parasitic nematode, and provides insight into the genomics of ivermectin response and population structure of O. volvulus. We argue that ivermectin response is a polygenically-determined quantitative trait (QT) whereby identical or related molecular pathways but not necessarily individual genes are likely to determine the extent of ivermectin response in different parasite populations. Furthermore, we propose that genetic drift rather than genetic selection of SOR is the underlying driver of population differentiation, which has significant implications for the emergence and potential spread of SOR within and between these parasite populations.

Onchocerciasis is a human parasitic disease endemic across large areas of Sub-Saharan Africa, where more than 99% of the estimated 100 million people globally at-risk live. The microfilarial stage of Onchocerca volvulus causes pathologies ranging from mild itching to visual impairment and ultimately, irreversible blindness. Mass administration of ivermectin kills microfilariae and has an anti-fecundity effect on adult worms by temporarily inhibiting the development in utero and/or release into the skin of new microfilariae, thereby reducing morbidity and transmission. Phenotypic and genetic changes in some parasite populations that have undergone multiple ivermectin treatments in Cameroon and Ghana have raised concern that sub-optimal response to ivermectin's anti-fecundity effect may increase in frequency, reducing the impact of ivermectin-based control measures. We used next generation sequencing of small pools of parasites to define genome-wide genetic differences between phenotypically characterised good and sub-optimal responder parasites from Cameroon and Ghana, and identified multiple regions of the genome that differentiated the response types. These regions were largely different between parasites from these two countries but revealed common molecular pathways that might be involved in determining the extent of response to ivermectin's anti-fecundity effect. These data reveal a more complex than previously described pattern of genetic diversity among O. volvulus populations that differ in their geography and response to ivermectin treatment.

Metabolite profiling of infection-associated metabolic markers of onchocerciasis

The global efforts for onchocerciasis elimination may require additional tools (safe micro and macrofilaricidal drugs, vaccines and biomarkers) as elimination efforts move toward the "end game". Efforts toward the identification of suitable biomarkers have focused on specific protein(s) and/or nucleic acids, but metabolites present an alternative option as they have limited half-lives and are the result of combinatorial effects. In comparison to previously used methodology of LC-MS for metabolomic approaches, we used a non-targeted capillary electrophoresis time-of-flight mass spectrometry (CE-TOFMS) to analyze the serum metabolic profiles of Ov-infected and -uninfected individuals (n=20). We identified 286 known metabolites (167 in the cation mode and 119 in the anion mode). In addition, putative metabolites were identified based on KEGG (51), HMDB (37) and HMT (6) databases. One hundred ten of these putative metabolites were quantified based on peak areas of internal standards and their ability to be mapped to known pathways (primary-, carbon-, lipid-, amino acid-, nucleotide and coenzyme-metabolism). Multivariate analysis demonstrated clustering and segregation of some of these metabolites to either the infected or control groups. The levels of serotonin, hypoxanthine, pipecolic acid and inosine were significantly elevated in those with onchocerciasis, whereas the levels of glycerophosphocholine, choline and adenine were significantly lower. This non-targeted metabolomic approach provides a global view of the metabolic variations that occur during Ov infection and thus allow the discovery of key metabolites (and associated pathways) that may serve as useful biomarkers in human onchocerciasis.

Discovery of Specific Antigens That Can Predict Microfilarial Intensity in Loa loa Infection

Antigen-based immunoassays are currently needed for point-of-care quantification of Loa loa microfilariae (mf). Coupling transcriptomic approaches with bioinformatic analysis, we have identified 11 specific putative proteins (coding mRNAs) with potential utility as biomarkers of patent (mf ⁺ ) L. loa infection. We successfully developed antigen capture immunoassays to quantify 2 (LOAG_14221 and LOAG_15846) of these proteins in individual plasma/serum samples. Of the 2 quantifiable circulating biomarkers, LOAG_14221 showed the highest degree of specificity, particularly with a monoclonal antibody-based immunoassay. Moreover, the levels of LOAG_14221 in L. loa mf ⁺ patients were positively correlated to the mf densities in the corresponding blood samples (r = 0.53 and P = 0.008 for polyclonal assay; r = 0.54 and P = 0.004 for monoclonal assay). Thus, LOAG_14221 is a very promising biomarker that will be exploited in a quantitative point-of-care immunoassay for determination of L. loa mf densities.

Factors affecting onchocerciasis transmission: lessons for infection control

INTRODUCTION

Onchocerca volvulus infects in excess of 15 million people. The vectors are Simulium blackflies, varieties of which differ in their ecologies, behavior and vectorial abilities. Control of the vectors and mass administrations of ivermectin have succeeded in reducing prevalences with elimination achieved in some foci, particularly in Central and southern America. In Africa, progress towards elimination has been less successful. Areas covered: Even with community directed treatment with ivermectin (CDTI), control has been difficult in African areas with initial prevalences in excess of 55%, especially if only annual treatments are dispensed. This is partly attributable to insufficient coverage, but the appearance of incipiently resistant non-responding parasites and lack of attention to vector biology in modeling and planning outcomes of intervention programmes have also played their parts, with recrudescence now appearing in some treated areas. Expert commentary: The biology of onchocerciasis is complex involving different vectors with differing abilities to transmit parasites, diverse pathologies related to geographical and parasite variations and endosymbionts in both parasite and vector. Modeling to predict epidemiological and control outcomes is addressing this complexity but more attention needs to be given to the vectors' roles to further understanding of where and when control measures will succeed.

Lessons from the genomes and transcriptomes of filarial nematodes

Human filarial infections are a leading cause of morbidity in the developing world. While a small arsenal of drugs exists to treat these infections, there remains a tremendous need for the development of additional interventions. Recent genome sequences and transcriptome analyses of filarial nematodes have provided novel biological insight and allowed for the prediction of novel drug targets as well as potential vaccine candidates. In this review, we discuss the currently available data, insights gained into the metabolism of these organisms, and how the filaria field can move forward by leveraging these data.

WormBase ParaSite - a comprehensive resource for helminth genomics

•

WormBase ParaSite is a new resource for helminth genomics.

•

The resource provides access to over 100 nematode and platyhelminth genomes.

•

The genomes are consistently annotated, organised and presented.

•

A variety of views and tools for exploring and querying the data are provided.

The number of publicly available parasitic worm genome sequences has increased dramatically in the past three years, and research interest in helminth functional genomics is now quickly gathering pace in response to the foundation that has been laid by these collective efforts. A systematic approach to the organisation, curation, analysis and presentation of these data is clearly vital for maximising the utility of these data to researchers. We have developed a portal called WormBase ParaSite (http://parasite.wormbase.org) for interrogating helminth genomes on a large scale. Data from over 100 nematode and platyhelminth species are integrated, adding value by way of systematic and consistent functional annotation (e.g. protein domains and Gene Ontology terms), gene expression analysis (e.g. alignment of life-stage specific transcriptome data sets), and comparative analysis (e.g. orthologues and paralogues). We provide several ways of exploring the data, including genome browsers, genome and gene summary pages, text search, sequence search, a query wizard, bulk downloads, and programmatic interfaces. In this review, we provide an overview of the back-end infrastructure and analysis behind WormBase ParaSite, and the displays and tools available to users for interrogating helminth genomic data.

Stage-Specific Transcriptome and Proteome Analyses of the Filarial Parasite Onchocerca volvulus and Its Wolbachia Endosymbiont

Onchocerciasis (river blindness) is a neglected tropical disease that has been successfully targeted by mass drug treatment programs in the Americas and small parts of Africa. Achieving the long-term goal of elimination of onchocerciasis, however, requires additional tools, including drugs, vaccines, and biomarkers of infection. Here, we describe the transcriptome and proteome profiles of the major vector and the human host stages (L1, L2, L3, molting L3, L4, adult male, and adult female) of Onchocerca volvulus along with the proteome of each parasitic stage and of its Wolbachia endosymbiont (wOv). In so doing, we have identified stage-specific pathways important to the parasite’s adaptation to its human host during its early development. Further, we generated a protein array that, when screened with well-characterized human samples, identified novel diagnostic biomarkers of O. volvulus infection and new potential vaccine candidates. This immunomic approach not only demonstrates the power of this postgenomic discovery platform but also provides additional tools for onchocerciasis control programs.

The global onchocerciasis (river blindness) elimination program will have to rely on the development of new tools (drugs, vaccines, biomarkers) to achieve its goals by 2025. As an adjunct to the completed genomic sequencing of O. volvulus, we used a comprehensive proteomic and transcriptomic profiling strategy to gain a comprehensive understanding of both the vector-derived and human host-derived parasite stages. In so doing, we have identified proteins and pathways that enable novel drug targeting studies and the discovery of novel vaccine candidates, as well as useful biomarkers of active infection.

Genomic diversity in Onchocerca volvulus and its Wolbachia endosymbiont

Ongoing elimination efforts have altered the global distribution of Onchocerca volvulus, the agent of river blindness, and further population restructuring is expected as efforts continue. Therefore, a better understanding of population genetic processes and their effect on biogeography is needed to support elimination goals. We describe O. volvulus genome variation in 27 isolates from the early 1990s (before widespread mass treatment) from four distinct locales: Ecuador, Uganda, the West African forest and the West African savanna. We observed genetic substructuring between Ecuador and West Africa and between the West African forest and savanna bioclimes, with evidence of unidirectional gene flow from savanna to forest strains. We identified forest:savanna-discriminatory genomic regions and report a set of ancestry informative loci that can be used to differentiate between forest, savanna and admixed isolates, which has not previously been possible. We observed mito-nuclear discordance possibly stemming from incomplete lineage sorting. The catalogue of the nuclear, mitochondrial and endosymbiont DNA variants generated in this study will support future basic and translational onchocerciasis research, with particular relevance for ongoing control programmes, and boost efforts to characterize drug, vaccine and diagnostic targets.