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

Abstract

BackgroundTreatment 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 findingsPooled 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/SignificanceThis 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.

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

“…With a few notable exceptions (Choi et al, 2016;Small et al, 2016;Doyle et al, 2017;Small et al, 2019), most genomic research for filarial nematodes has focused on predictions of gene function. The aim of many of these "functional" predictions is to identify targets for new drug discovery, usually in the absence of methods for functional genomic validation of targets.…”

A Whole Genome Re-Sequencing Based GWA Analysis Reveals Candidate Genes Associated with Ivermectin Resistance in Haemonchus contortus

“…However, the extensive and uncontrolled use of IVM in livestock treatment has led to widespread resistance in gastrointestinal nematodes, such as H. contortus [7,8]. Alongside of veterinary sector where the rapid emergence of resistance against both single and multiple drug classes has been widely reported [7,9], there are also growing concerns regarding the emergence of resistance against the drugs used in helminths control in human [10][11][12]. The problem of IVM resistance has also been reported for pests [10,[13][14][15].…”

Economic evaluations of onchocerciasis interventions: a systematic review and research needs

“…) have undergone extensive refinement in recent years to better capture parasite population dynamics during interventions and to produce more robust projections on the likelihood of elimination. These modelling efforts will be particularly useful in identifying epidemiological and programmatic circumstances in which alternative strategies will be required to reach elimination (for example in highly endemic settings, settings with suboptimal responses to ivermectin , or where the initiation of programmes has previously been delayed ). Moreover, transmission models (as opposed to so‐called ‘static’ models) capture explicitly our current understanding of the changing parasite dynamics during interventions .…”