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

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.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.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 transmission in Ghana: Persistence under different control strategies and the role of the simuliid vectors

Background: The World Health Organization (WHO) aims at eliminating onchocerciasis by 2020 in selected African countries. Current control focuses on community-directed treatment with ivermectin (CDTI). In Ghana, persistent transmission has been reported despite long-term control. We present spatial and temporal patterns of onchocerciasis transmission in relation to ivermectin treatment history. Methodology/Principal Findings: Host-seeking and ovipositing blackflies were collected from seven villages in four regions of Ghana with 3–24 years of CDTI at the time of sampling. A total of 16,443 flies was analysed for infection; 5,812 (35.3%) were dissected for parity (26.9% parous). Heads and thoraces of 12,196 flies were dissected for Onchocerca spp. and DNA from 11,122 abdomens was amplified using Onchocerca primers. A total of 463 larvae (0.03 larvae/fly) from 97 (0.6%) infected and 62 (0.4%) infective flies was recorded; 258 abdomens (2.3%) were positive for Onchocerca DNA. Infections (all were O. volvulus) were more likely to be detected in ovipositing flies. Transmission occurred, mostly in the wet season, at Gyankobaa and Bosomase, with transmission potentials of, respectively, 86 and 422 L3/person/month after 3 and 6 years of CDTI. The numbers of L3/1,000 parous flies at these villages were over 100times the WHO threshold of one L3/1,000 for transmission control. Vector species influenced transmission parameters. At Asubende, the number of L3/1,000 ovipositing flies (1.4, 95% CI = 0–4) also just exceeded the threshold despite extensive vector control and 24 years of ivermectin distribution, but there were no infective larvae in host-seeking flies. Conclusions/Significance: Despite repeated ivermectin treatment, evidence of O. volvulus transmission was documented in all seven villages and above the WHO threshold in two. Vector species influences transmission through biting and parous rates and vector competence, and should be included in transmission models. Oviposition traps could augment vector collector methods for monitoring and surveillance

Model-based geostatistical mapping of the prevalence of onchocerca volvulus in West Africa.

Background: The initial endemicity (pre-control prevalence) of onchocerciasis has been shown to be an important determinant of the feasibility of elimination by mass ivermectin distribution. We present the first geostatistical map of microfilarial prevalence in the former Onchocerciasis Control Programme in West Africa (OCP) before commencement of antivectorial and antiparasitic interventions. Methods and Findings: Pre-control microfilarial prevalence data from 737 villages across the 11 constituent countries in the OCP epidemiological database were used as ground-truth data. These 737 data points, plus a set of statistically selected environmental covariates, were used in a Bayesian model-based geostatistical (B-MBG) approach to generate a continuous surface (at pixel resolution of 5 km x 5km) of microfilarial prevalence in West Africa prior to the commencement of the OCP. Uncertainty in model predictions was measured using a suite of validation statistics, performed on bootstrap samples of held-out validation data. The mean Pearson’s correlation between observed and estimated prevalence at validation locations was 0.693; the mean prediction error (average difference between observed and estimated values) was 0.77%, and the mean absolute prediction error (average magnitude of difference between observed and estimated values) was 12.2%. Within OCP boundaries, 17.8 million people were deemed to have been at risk, 7.55 million to have been infected, and mean microfilarial prevalence to have been 45% (range: 2–90%) in 1975. Conclusions and Significance: This is the first map of initial onchocerciasis prevalence in West Africa using B-MBG. Important environmental predictors of infection prevalence were identified and used in a model out-performing those without spatial random effects or environmental covariates. Results may be compared with recent epidemiological mapping efforts to find areas of persisting transmission. These methods may be extended to areas where data are sparse, and may be used to help inform the feasibility of elimination with current and novel tools

Designing antifilarial drug trials using clinical trial simulators

Lymphatic filariasis and onchocerciasis are neglected tropical diseases (NTDs) targeted for elimination by mass (antifilarial) drug administration. These drugs are predominantly active against the microfilarial progeny of adult worms. New drugs or combinations are needed to improve patient therapy and to enhance the effectiveness of interventions in persistent hotspots of transmission. Several therapies and regimens are currently in (pre-)clinical testing. Clinical trial simulators (CTSs) project patient outcomes to inform the design of clinical trials but have not been widely applied to NTDs, where their resource-saving payoffs could be highly beneficial. We demonstrate the utility of CTSs using our individual-based onchocerciasis transmission model (EPIONCHO-IBM) that projects trial outcomes of a hypothetical macrofilaricidal drug. We identify key design decisions that influence the power of clinical trials, including participant eligibility criteria and post-treatment follow-up times for measuring infection indicators. We discuss how CTSs help to inform target product profiles

Oxamniquine resistance alleles are widespread in Old World Schistosoma mansoni and predate drug deployment

Do mutations required for adaptation occur de novo, or are they segregating within populations as standing genetic variation? This question is key to understanding adaptive change in nature, and has important practical consequences for the evolution of drug resistance. We provide evidence that alleles conferring resistance to oxamniquine (OXA), an antischistosomal drug, are widespread in natural parasite populations under minimal drug pressure and predate OXA deployment. OXA has been used since the 1970s to treat Schistosoma mansoni infections in the New World where S. mansoni established during the slave trade. Recessive loss-of-function mutations within a parasite sulfotransferase (SmSULT-OR) underlie resistance, and several verified resistance mutations, including a deletion (p.E142del), have been identified in the New World. Here we investigate sequence variation in SmSULT-OR in S. mansoni from the Old World, where OXA has seen minimal usage. We sequenced exomes of 204 S. mansoni parasites from West Africa, East Africa and the Middle East, and scored variants in SmSULT-OR and flanking regions. We identified 39 non-synonymous SNPs, 4 deletions, 1 duplication and 1 premature stop codon in the SmSULT-OR coding sequence, including one confirmed resistance deletion (p.E142del). We expressed recombinant proteins and used an in vitro OXA activation assay to functionally validate the OXA-resistance phenotype for four predicted OXA-resistance mutations. Three aspects of the data are of particular interest: (i) segregating OXA-resistance alleles are widespread in Old World populations (4.29–14.91% frequency), despite minimal OXA usage, (ii) two OXA-resistance mutations (p.W120R, p.N171IfsX28) are particularly common (>5%) in East African and Middle-Eastern populations, (iii) the p.E142del allele has identical flanking SNPs in both West Africa and Puerto Rico, suggesting that parasites bearing this allele colonized the New World during the slave trade and therefore predate OXA deployment. We conclude that standing variation for OXA resistance is widespread in S. mansoni