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 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.

Parasitism is an important life history strategy in many metazoan taxa. This is particularly true of the Phylum Nematoda, in which parasitism has evolved independently at least nine times. The apparent ease with which parasitism has evolved amongst nematodes may, in part, be due to a feature of nematode development acting as a pre-adaptation for the transition from a free-living to a parasitic life history. One candidate pre-adaptive feature for evolution in terrestrial nematodes is the dauer larva, a developmentally arrested morph formed in response to environmental signals.

Wolbachia are intracellular bacteria found in arthropods and several filarial nematode species. The filarial Wolbachia have been proposed to be involved in the immunopathology associated with onchocerciasis. Higher Wolbachia-to-nematode ratios have been reported in the savannah-ecotype compared to the forest-ecotype, and have been interpreted as consistent with a correlation between Wolbachia density and disease severity. However, factors such as geographic stratification and ivermectin drug exposure can lead to significant genetic heterogeneity in the nematode host populations, so we investigated whether Wolbachia copy number variation is also associated with these underlying factors. Genomic DNA was prepared from single adult nematodes representing forest and savannah ecotypes sampled from Togo, Ghana, Côte d’Ivoire and Mali. A qPCR assay was developed to measure the number of Wolbachia genome(s) per nematode genome. Next-generation sequencing (NGS) was also used to measure relative Wolbachia copy number, and independently verify the qPCR assay. Significant variation was observed within the forest (range: 0.02 to 452.99; median: 10.58) and savannah (range: 0.01 to 1106.25; median: 9.10) ecotypes, however, no significant difference between ecotypes (P = 0.645) was observed; rather, strongly significant Wolbachia variation was observed within and between the nine study communities analysed (P = 0.021), independent of ecotype. Analysis of ivermectin-treated and untreated nematodes by qPCR showed no correlation (P = 0.869); however, an additional analysis of a subset of the nematodes by qPCR and NGS revealed a correlation between response to ivermectin treatment and Wolbachia copy number (P = 0.020). This study demonstrates that extensive within and between population variation exists in the Wolbachia content of individual adult O. volvulus. The origin and functional significance of such variation (up to ~ 100,000-fold between worms; ~10 to 100-fold between communities) in the context of the proposed mutualistic relationship between the worms and the bacteria, and between the presence of Wolbachia and clinical outcome of infection, remains unclear. These data do not support a correlation between Wolbachia copy number and forest or savannah ecotype, and may have implications for the development of anti-Wolbachia drugs as a macrofilaricidal treatment of onchocerciasis. The biological significance of a correlation between variation in Wolbachia copy number and ivermectin response remains unexplained.

Genetic surveillance of the human filarial parasite, Onchocerca volvulus, from onchocerciasis endemic regions will ideally focus on genotyping individual infective larval stages collected from their intermediate host, Simuliid blackflies. However, blackflies also transmit other Onchocerca species, including the cattle parasite O. ochengi, which are difficult to distinguish from the human parasite based on morphological characteristics alone. This study describes a versatile approach to discriminate between O. volvulus and O. ochengi that is demonstrated using parasite infective larvae dissected from blackflies. A speciation assay was designed based on genetic differentiation between O. volvulus and O. ochengi mitochondrial genome sequences that can be performed in high-throughput high-resolution melt (HRM)- or lower throughput conventional restriction fragment length polymorphism (RFLP) analyses. This assay was validated on 185 Onchocerca larvae dissected from blackflies captured from 14 communities in Ghana throughout 2011−2013. The frequency of O. ochengi was approximately 67 % of all larvae analysed, which is significantly higher than previously reported in this region. Furthermore, the species distribution was not uniform throughout the study region, with 25 %, 47 % and 93 % of O. volvulus being found in the western-most (Black Volta, Tain and Tombe), the central (Pru) and eastern-most (Daka) river basins, respectively. This tool provides a simple and cost-effective approach to determine the identity and distribution of two Onchocerca species, and will be valuable for future genetic studies that focus on parasites collected from blackflies. The results presented highlight the need to discriminate Onchocerca species in transmission studies, as the frequency of each species varied significantly between the communities studied.

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.