Genomic diversity in Onchocerca volvulus and its Wolbachia endosymbiont

Whole genome amplification and sequencing of individual Dirofilaria immitis microfilariae

Dirofilaria immitis is a filarial parasitic nematode of veterinary significance. With the emergence of drug-resistant isolates in the USA, it is imperative to determine the likelihood of resistance occurring in other regions of the world. One approach is to conduct population genetic studies across an extensive geographical range, and to sequence the genomes of individual worms to understand genome-wide genetic variation associated with resistance. The immature life stages of D. immitis found in the host blood are more accessible and less invasive to sample compared to extracting adult stages from the host heart. To assess the use of immature life stages for population genetic analyses, we have performed whole genome amplification and whole-genome sequencing on nine (n = 9) individual D. immitis microfilaria samples isolated from dog blood. On average, less than 1% of mapped reads aligned to each D. immitis genome (nuclear, mitochondrial, and Wolbachia endosymbiont). For the dog genome, an average of over 99% of mapped reads aligned to the nuclear genome and less than 1% aligned to the mitochondrial genome. The average coverage for all D. immitis genomes and the dog nuclear genome was less than 1, while the dog mitochondrial genome had an average coverage of 2.87. The overwhelming proportion of sequencing reads mapping to the dog host genome can be attributed to residual dog blood cells in the microfilariae samples. These results demonstrate the challenges of conducting genome-wide studies on individual immature parasite life stages, particularly in the presence of extraneous host DNA.

Distinguishing recrudescence from reinfection in lymphatic filariasis

BACKGROUND

The Global Program to Eliminate Lymphatic Filariasis (GPELF) is the largest public health program based on mass drug administration (MDA). Despite decades of MDA, ongoing transmission in some countries remains a challenge. To optimise interventions, it is critical to differentiate between recrudescence and new infections. Since adult filariae are inaccessible in humans, deriving a method that relies on the offspring microfilariae (mf) is necessary.

METHODS

We developed a genome amplification and kinship analysis-based approach using Brugia malayi samples from gerbils, and applied it to analyse Wuchereria bancrofti mf from humans in Côte d'Ivoire. We examined the pre-treatment genetic diversity in 269 mf collected from 18 participants, and further analysed 1-year post-treatment samples of 74 mf from 4 participants. Hemizygosity of the male X-chromosome allowed for direct inference of haplotypes, facilitating robust maternal parentage inference. To enrich parasite DNA from samples contaminated with host DNA, a whole-exome capture panel was created for W. bancrofti.

FINDINGS

By reconstructing and temporally tracking sibling relationships across pre- and post-treatment samples, we differentiated between new and established maternal families, suggesting reinfection in one participant and recrudescence in three participants. The estimated number of reproductively active adult females ranged between 3 and 11 in the studied participants. Population structure analysis revealed genetically distinct parasites in Côte d'Ivoire compared to samples from other countries. Exome capture identified protein-coding variants with ∼95% genotype concordance rate.

INTERPRETATION

We have generated resources to facilitate the development of molecular genetic tools that can estimate adult worm burdens and monitor parasite populations, thus providing essential information for the successful implementation of GPELF.

FUNDING

This work was financially supported by the Bill and Melinda Gates Foundation (https://www.gatesfoundation.org) under grant OPP1201530 (Co-PIs PUF & Gary J. Weil). B. malayi parasite material was generated with support of the Foundation for Barnes Jewish Hospital (PUF). In addition, the development of computational methods was supported by the National Institutes of Health under grants AI144161 (MM) and AI146353 (MM). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Genome-based tools for onchocerciasis elimination: utility of the mitochondrial genome for delineating Onchocerca volvulus transmission zones

National programs in Africa have expanded their objectives from control of onchocerciasis (river blindness) as a public health problem to elimination of parasite transmission, motivated by the reduction of Onchocerca volvulus infection prevalence in many African meso- and hyperendemic areas due to mass drug administration of ivermectin (MDAi). Given the large, contiguous hypo-, meso-, and hyperendemic areas, sustainable elimination of onchocerciasis in sub-Saharan Africa requires delineation of geographic boundaries for parasite transmission zones, so that programs can consider the risk of parasite re-introduction through vector or human migration from areas with ongoing transmission when making decisions to stop MDAi. We propose that transmission zone boundaries can be delineated by characterising the parasite genetic population structure within and between potential zones. We analysed whole mitochondrial genome sequences of 189 O. volvulus adults to determine the pattern of genetic similarity across three West African countries: Ghana, Mali, and Côte d'Ivoire. Population genetic structure indicates that parasites from villages near the Pru, Daka, and Black Volta rivers in central Ghana belong to one parasite population, indicating that the assumption that river basins constitute individual transmission zones is not supported by the data. Parasites from Mali and Côte d'Ivoire are genetically distinct from those from Ghana. This research provides the basis for developing tools for elimination programs to delineate transmission zones, to estimate the risk of parasite re-introduction via vector or human movement when intervention is stopped in one area while transmission is ongoing in others, to identify the origin of infections detected post-treatment cessation, and to investigate whether persisting prevalence despite ongoing interventions in one area is due to parasites imported from others.

Chorioretinitis among Immigrant and Travellers. Comment on Mansour et al. Presumed Onchocerciasis Chorioretinitis Spilling over into North America, Europe and Middle East. Diagnostics 2023, 13, 3626

We read, with interest, the paper by Mansour et al [...].

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.

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.

Distinguishing recrudescence from reinfection in lymphatic filariasis: a genomics-based approach for monitoring worm burden

Background

The Global Program to Eliminate Lymphatic Filariasis is the largest public health program based on mass drug administration (MDA). Despite decades of MDA, ongoing transmission in some countries remains a challenge. To optimize interventions, it is essential to differentiate between recrudescence (poor drug response and persistent infection) and new infections (ongoing transmission). Since adult filariae are inaccessible in humans, an approach that relies on genotyping the offspring microfilariae (mf) is required.

Methods

We utilized Brugia malayi adults and mf obtained from gerbils with a known pedigree to develop and validate our whole-genome amplification and kinship analysis approach. We then sequenced the genomes of Wuchereria bancrofti mf from infected humans from Côte d'Ivoire (CDI), characterized the population genetic diversity, and made inferences about the adult breeders. We developed a whole-exome capture panel for W. bancrofti to enrich parasite nuclear DNA from lower-quality samples contaminated with host DNA.

Results

We established a robust analysis pipeline using B. malayi adult and mf. We estimated the pre-treatment genetic diversity in W. bancrofti from 269 mf collected from 18 individuals, and further analyzed 1-year post-treatment samples of 74 mf from 4 individuals. By reconstructing and temporally tracking sibling relationships across pre- and post-treatment samples, we differentiated between new and established maternal families, suggesting reinfection in one subject and recrudescence in three subjects. Estimated reproductively active adult females ranged between 3 and 9 in the studied subjects. Hemizygosity of the male X-chromosome allowed for direct inference of haplotypes, facilitating robust maternal parentage inference, even when the genetic diversity was low. Population structure analysis revealed genetically distinct parasites among our CDI samples. Sequence composition and variant analysis of whole-exome libraries showed that the hybridization capture approach can effectively enrich parasite nuclear DNA and identify protein-coding variants with ∼95% genotype concordance rate.

Conclusions

We have generated resources to facilitate development of field-deployable genotyping tools that can estimate worm burdens and monitor parasite populations. These tools are essential for the success of lymphatic filariasis MDA programs. With further expansion of the databases to include geographically diverse samples, we will be able to spatially track parasite movement associated with host/vector migration.

Procedure for Handling and Storage of Onchocerca volvulus Microfilariae Obtained from Skin Snips for Downstream Genetic Work

WHO and endemic countries target elimination of transmission of Onchocerca volvulus, the parasite causing onchocerciasis. Population genetic analysis of O. volvulus may provide data to improve the evidence base for decisions on when, where, and for how long to deploy which interventions and post-intervention surveillance to achieve elimination. Development of necessary methods and tools requires parasites suitable for genetic analysis. Based on our experience with microfilariae obtained from different collaborators, we developed a microfilariae transfer procedure for large-scale studies in the Democratic Republic of Congo (DRC) comparing safety and efficacy of ivermectin, the mainstay of current onchocerciasis elimination strategies, and moxidectin, a new drug. This procedure is designed to increase the percentage of microfilariae in skin snips suitable for genetic analysis, improve assignment to metadata, and minimize time and materials needed by the researchers collecting the microfilariae. Among 664 microfilariae from South Sudan, 35.7% and 39.5% failed the mitochondrial and nuclear qPCR assay. Among the 576 microfilariae from DRC, 16.0% and 16.7% failed these assays, respectively. This difference may not only be related to the microfilariae transfer procedure but also to other factors, notably the ethanol concentration in the tubes in which microfilariae were stored (64% vs. ≥75%).

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.

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

High Prevalence of Glaucoma among Patients in an Onchocerciasis Endemic Area (Mahenge, Tanzania)

Onchocerciasis is known to cause skin lesions and blindness, but there is also epidemiological evidence that onchocerciasis is associated with epilepsy, including nodding syndrome. We carried out ocular exams in persons with epilepsy in Mahenge, an onchocerciasis endemic area with a high prevalence of epilepsy in Tanzania. We recruited 278 consecutive persons with epilepsy attending the epilepsy clinic at Mahenge hospital and satellite clinics in rural villages. They underwent a general physical and a detailed ocular examination and were tested for onchocerciasis Ov16 IgG4 antibodies. Glaucoma was defined by a raised intraocular pressure above 21 mmHg with evidence of typical glaucomatous disc changes in one or both eyes. Among the 278 participants, median age 27 (IQR 21–38) years, 55.4% were female; 151/210 (71.9%) (95% CI: 65.3–77.9) were Ov16 positive. The most frequent ophthalmic lesions were glaucoma (33.1%), vitreous opacities (6.5%) and cataracts (2.9%). In multivariate analysis, glaucoma (adjusted IRR = 1.46; 95% CI: 1.24–1.70) and age (adjusted IRR = 1.01; 95% CI: 1.01–1.02) were significantly associated with onchocerciasis. In conclusion, a high prevalence of glaucoma was observed among Ov16 positive persons with epilepsy. Persons with epilepsy with O. volvulus infection should undergo screening for glaucoma to prevent one of the causes of preventable blindness.

Population genomics of ancient and modern Trichuris trichiura

The neglected tropical disease trichuriasis is caused by the whipworm Trichuris trichiura, a soil-transmitted helminth that has infected humans for millennia. Today, T. trichiura infects as many as 500 million people, predominantly in communities with poor sanitary infrastructure enabling sustained faecal-oral transmission. Using whole-genome sequencing of geographically distributed worms collected from human and other primate hosts, together with ancient samples preserved in archaeologically-defined latrines and deposits dated up to one thousand years old, we present the first population genomics study of T. trichiura. We describe the continent-scale genetic structure between whipworms infecting humans and baboons relative to those infecting other primates. Admixture and population demographic analyses support a stepwise distribution of genetic variation that is highest in Uganda, consistent with an African origin and subsequent translocation with human migration. Finally, genome-wide analyses between human samples and between human and non-human primate samples reveal local regions of genetic differentiation between geographically distinct populations. These data provide insight into zoonotic reservoirs of human-infective T. trichiura and will support future efforts toward the implementation of genomic epidemiology of this globally important helminth.

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.

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.

Predicting the environmental suitability for onchocerciasis in Africa as an aid to elimination planning

Recent evidence suggests that, in some foci, elimination of onchocerciasis from Africa may be feasible with mass drug administration (MDA) of ivermectin. To achieve continental elimination of transmission, mapping surveys will need to be conducted across all implementation units (IUs) for which endemicity status is currently unknown. Using boosted regression tree models with optimised hyperparameter selection, we estimated environmental suitability for onchocerciasis at the 5 × 5-km resolution across Africa. In order to classify IUs that include locations that are environmentally suitable, we used receiver operating characteristic (ROC) analysis to identify an optimal threshold for suitability concordant with locations where onchocerciasis has been previously detected. This threshold value was then used to classify IUs (more suitable or less suitable) based on the location within the IU with the largest mean prediction. Mean estimates of environmental suitability suggest large areas across West and Central Africa, as well as focal areas of East Africa, are suitable for onchocerciasis transmission, consistent with the presence of current control and elimination of transmission efforts. The ROC analysis identified a mean environmental suitability index of 0·71 as a threshold to classify based on the location with the largest mean prediction within the IU. Of the IUs considered for mapping surveys, 50·2% exceed this threshold for suitability in at least one 5 × 5-km location. The formidable scale of data collection required to map onchocerciasis endemicity across the African continent presents an opportunity to use spatial data to identify areas likely to be suitable for onchocerciasis transmission. National onchocerciasis elimination programmes may wish to consider prioritising these IUs for mapping surveys as human resources, laboratory capacity, and programmatic schedules may constrain survey implementation, and possibly delaying MDA initiation in areas that would ultimately qualify.

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.

Integrating Multiple Biomarkers to Increase Sensitivity for the Detection of Onchocerca volvulus Infection

BACKGROUND

Serological assessments for human onchocerciasis are based on IgG4 reactivity against the OV-16 antigen, with sensitivities of 60-80%. We have previously identified 7 novel proteins that could improve serodiagnosis.

METHODS

IgG4 responses to these 7 proteins were assessed by luciferase immunoprecipitation (LIPS) and enzyme-linked immunosorbent (ELISA) immunoassays.

RESULTS

OVOC10469 and OVOC3261 were identified as the most promising candidates by IgG4-based immunoassays with sensitivities of 53% for rOVOC10469 and 78% for rOVOC3261 while specificity for each was >99%. These 2 antigens in combination with OV-16 increased the sensitivity for patent infections to 94%. The kinetics of appearance of these IgG4 responses based on experimentally infected non-human primates indicated that they were microfilarial- driven. Further, the IgG4 responses to both OVOC10469 and OVOC3261 (as well as to OV-16) drop significantly (p<0.05) following successful treatment for onchocerciasis. A prototype lateral flow rapid diagnostic test to detect IgG4 to both Ov-16 and OVOC3261 was developed and tested demonstrating an overall 94% sensitivity.

CONCLUSION

The combined use of rOVOC3261 with OV-16 improved serologic assessment of O. volvulus infection, a current unmet need toward the goal of elimination of transmission of O. volvulus.

Taking the strain out of onchocerciasis? A reanalysis of blindness and transmission data does not support the existence of a savannah blinding strain of onchocerciasis in West Africa

Onchocerciasis (also known as 'river blindness'), is a neglected tropical disease (NTD) caused by the (Simulium-transmitted) filarial nematode Onchocerca volvulus. The occurrence of 'blinding' (savannah) and non-blinding (forest) parasite strains and the existence of corresponding, locally adapted Onchocerca-Simulium complexes were postulated to explain greater blindness prevalence in savannah than in forest foci. As a result, the World Health Organization (WHO) Onchocerciasis Control Programme in West Africa (OCP) focused anti-vectorial and anti-parasitic interventions in savannah endemic areas. In this paper, village-level data on blindness prevalence, microfilarial prevalence, and transmission intensity (measured by the annual transmission potential, the number of infective, L3, larvae per person per year) were extracted from 16 West-Central Africa-based publications, and analysed according to habitat (forest, forest-savannah mosaic, savannah) to test the dichotomous strain hypothesis in relation to blindness. When adjusting for sample size, there were no statistically significant differences in blindness prevalence between the habitats (one-way ANOVA, P=0.68, mean prevalence for forest=1.76±0.37 (SE); mosaic=1.49±0.38; savannah=1.89±0.26). The well-known relationship between blindness prevalence and annual transmission potential for savannah habitats was confirmed and shown to hold for (but not to be statistically different from) forest foci (excluding data from southern Côte d'Ivoire, in which blindness prevalence was significantly lower than in other West African forest communities, but which had been the focus of studies leading to the strain-blindness hypothesis that was accepted by OCP planners). We conclude that the evidence for a savannah blinding onchocerciasis strain in simple contrast with a non-blinding forest strain is equivocal. A re-appraisal of the strain hypothesis to explain patterns of ocular disease is needed to improve understanding of onchocerciasis epidemiology and disease burden estimates in the light of the WHO 2030 goals for onchocerciasis.

Essential Amino Acid Enrichment and Positive Selection Highlight Endosymbiont's Role in a Global Virus-Vectoring Pest

Host-associated microbes display remarkable convergence in genome repertoire resulting from selection to supplement missing host functions. Nutritional supplementation has been proposed in the verrucomicrobial endosymbiont Xiphinematobacter sp., which lives within a globally widespread group of plant-parasitic nematodes that vector damaging nepoviruses to plants. Only one genome sequence has been published from this symbiont, leaving unanswered questions about its diversity, host range, role, and selective pressures within its hosts. Because its hosts are exceptionally resistant to culturing, this symbiont is best studied through advanced genomic approaches. To analyze the role of Xiphinematobacter sp. in its host, sequencing was performed on nematode communities, and then genomes were extracted for comparative genomics, gene ontology enrichment tests, polymorphism analysis, de Bruijn-based genome-wide association studies, and tests of pathway- and site-specific selection on genes predicted play a role in the symbiosis. Results showed a closely clustered set of Xiphinematobacter isolates with reduced genomes of ∼917 kbp, for which a new species was proposed. Symbionts shared only 2.3% of genes with outgroup Verrucomicrobia, but comparative analyses showed high conservation of all 10 essential amino acid (EAA) biosynthesis pathways plus several vitamin pathways. These findings were supported by gene ontology enrichment tests and high polymorphisms in these pathways compared with background. Genome-wide association analysis confirmed high between-species fixation of alleles with significant functional enrichment for EAA and thiamine synthesis. Strong positive selection was detected on sites within these pathways, despite several being under increased purifying selection. Together, these results suggest that supplementation of EAAs missing in the host diet may drive this widespread symbiosis.IMPORTANCE Xiphinematobacter spp. are distinctly evolved intracellular symbionts in the phylum Verrucomicrobia, which includes the important human gut-associated microbe Akkermansia muciniphila and many highly abundant free-living soil microbes. Like Akkermansia sp., Xiphinematobacter sp. is obligately associated with the gut of its hosts, which in this case consists of a group of plant-parasitic nematodes that are among the top 10 most destructive species to global agriculture, by vectoring plant viruses. This study examined the hypothesis that the key to this symbiont's stable evolutionary association with its host is through provisioning nutrients that its host cannot make that may be lacking in the nematode's plant phloem diet, such as essential amino acids and several vitamins. The significance of our research is in demonstrating, using population genomics, the signatures of selective pressure on these hypothesized roles to ultimately learn how this independently evolved symbiont functionally mirrors symbionts of phloem-feeding insects.

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.

Genetic epidemiology of lymphatic filariasis in American Samoa after mass drug administration

Over 892 million people in 48 countries are at risk of infection by nematodes that cause lymphatic filariasis. As part of the Global Programme to Eliminate Lymphatic Filariasis, mass drug administration is distributed to communities until surveillance indicates infection rates are below target prevalence thresholds. In some countries, including American Samoa, lymphatic filariasis transmission persists despite years of mass drug administration and/or has resurged after cessation. Nothing is known about the population genetics of Wuchereria bancrofti worms in Polynesia, or whether local transmission is persisting and/or increasing due to inadequate mass drug administration coverage, expansion from residual hotspots, reintroduction from elsewhere, or a combination. We extracted DNA from microfilariae on blood slides collected during prevalence surveys in 2014 and 2016, comprising 31 pools of five microfilariae from 22 persons living in eight villages. We sequenced 1104 bp across three mitochondrial markers (ND4, COI, CYTB). We quantified parasite genetic differentiation using variant calls and estimated haplotypes using principal components analysis, F-statistics, and haplotype networks. Of the variants called, all but eight were shared across the main island of Tutuila, and three of those were from a previously described hotspot village, Fagali'i. Genotypic data did not support population genetic structure among regions or villages in 2016, although differences were observed between worms collected in Fagali'i in 2014 and those from 2016. Because estimated haplotype frequency varied between villages, these statistics suggested genetic differentiation, but were not consistent among villages. Finally, haplotype networks demonstrated American Samoan sequence clusters were related to previously published sequences from Papua New Guinea. These are, to our knowledge, the first reports of W. bancrofti genetic variation in Polynesia. The resurgent parasites circulating on the main island of American Samoa represent a single population. This study is the first step towards investigating how parasite population structure might inform strategies to manage resurgence and elimination of lymphatic filariasis.

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.

Large scale genome reconstructions illuminate Wolbachia evolution

Wolbachia is an iconic example of a successful intracellular bacterium. Despite its importance as a manipulator of invertebrate biology, its evolutionary dynamics have been poorly studied from a genomic viewpoint. To expand the number of Wolbachia genomes, we screen over 30,000 publicly available shotgun DNA sequencing samples from 500 hosts. By assembling over 1000 Wolbachia genomes, we provide a substantial increase in host representation. Our phylogenies based on both core-genome and gene content provide a robust reference for future studies, support new strains in model organisms, and reveal recent horizontal transfers amongst distantly related hosts. We find various instances of gene function gains and losses in different super-groups and in cytoplasmic incompatibility inducing strains. Our Wolbachia-host co-phylogenies indicate that horizontal transmission is widespread at the host intraspecific level and that there is no support for a general Wolbachia-mitochondrial synchronous divergence.

By greatly expanding the number of assembled genomes for Wolbachia (a group of intracellular bacteria) and constructing robust phylogenies, this study finds strong rate heterogeneity among Wolbachiapopulations and no support for synchronous divergence between Wolbachia and host mitochondria.

Gene duplications, divergence and recombination shape adaptive evolution of the fish ectoparasite Gyrodactylus bullatarudis

Determining the molecular basis of parasite adaptation to its host is an important component in understanding host-parasite coevolution and the epidemiology of parasitic infections. Here, we investigate short- and long-term adaptive evolution in the eukaryotic parasite Gyrodactylus bullatarudis infecting Caribbean guppies (Poecilia reticulata), by comparing the reference genome of Tobagonian G. bullatarudis with other Platyhelminthes, and by analysing resequenced samples from local Trinidadian populations. At the macroevolutionary timescale, we observed duplication of G-protein and serine proteases genes, which are probably important in host-parasite arms races. Serine protease also showed strong evidence of ongoing, diversifying selection at the microevolutionary timescale. Furthermore, our analyses revealed that a hybridization event, involving two divergent genomes, followed by recombination has dramatically affected the genetic composition of Trinidadian populations. The recombinant genotypes invaded Trinidad and replaced local parasites in all populations. We localized more than 300 genes in regions fixed in local populations for variants of different origin, possibly due to diversifying selection pressure from local host populations. In addition, around 70 genes were localized in regions identified as heterozygous in some, but not all, individuals. This pattern is consistent with a very recent spread of recombinant parasites. Overall, our results are consistent with the idea that recombination between divergent genomes can result in particularly successful parasites.

Conserved novel ORFs in the mitochondrial genome of the ctenophore Beroe forskalii

To date, five ctenophore species' mitochondrial genomes have been sequenced, and each contains open reading frames (ORFs) that if translated have no identifiable orthologs. ORFs with no identifiable orthologs are called unidentified reading frames (URFs). If truly protein-coding, ctenophore mitochondrial URFs represent a little understood path in early-diverging metazoan mitochondrial evolution and metabolism. We sequenced and annotated the mitochondrial genomes of three individuals of the beroid ctenophore Beroe forskalii and found that in addition to sharing the same canonical mitochondrial genes as other ctenophores, the B. forskalii mitochondrial genome contains two URFs. These URFs are conserved among the three individuals but not found in other sequenced species. We developed computational tools called pauvre and cuttlery to determine the likelihood that URFs are protein coding. There is evidence that the two URFs are under negative selection, and a novel Bayesian hypothesis test of trinucleotide frequency shows that the URFs are more similar to known coding genes than noncoding intergenic sequence. Protein structure and function prediction of all ctenophore URFs suggests that they all code for transmembrane transport proteins. These findings, along with the presence of URFs in other sequenced ctenophore mitochondrial genomes, suggest that ctenophores may have uncharacterized transmembrane proteins present in their mitochondria.

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.

De novo Assembly of the Brugia malayi Genome Using Long Reads from a Single MinION Flowcell

Filarial nematode infections cause a substantial global disease burden. Genomic studies of filarial worms can improve our understanding of their biology and epidemiology. However, genomic information from field isolates is limited and available reference genomes are often discontinuous. Single molecule sequencing technologies can reduce the cost of genome sequencing and long reads produced from these devices can improve the contiguity and completeness of genome assemblies. In addition, these new technologies can make generation and analysis of large numbers of field isolates feasible. In this study, we assessed the performance of the Oxford Nanopore Technologies MinION for sequencing and assembling the genome of Brugia malayi, a human parasite widely used in filariasis research. Using data from a single MinION flowcell, a 90.3 Mb nuclear genome was assembled into 202 contigs with an N50 of 2.4 Mb. This assembly covered 96.9% of the well-defined B. malayi reference genome with 99.2% identity. The complete mitochondrial genome was obtained with individual reads and the nearly complete genome of the endosymbiotic bacteria Wolbachia was assembled alongside the nuclear genome. Long-read data from the MinION produced an assembly that approached the quality of a well-established reference genome using comparably fewer resources.

The global diversity of Haemonchus contortus is shaped by human intervention and climate

Haemonchus contortus is a haematophagous parasitic nematode of veterinary interest. We have performed a survey of its genome-wide diversity using single-worm whole genome sequencing of 223 individuals sampled from 19 isolates spanning five continents. We find an African origin for the species, together with evidence for parasites spreading during the transatlantic slave trade and colonisation of Australia. Strong selective sweeps surrounding the β-tubulin locus, a target of benzimidazole anthelmintic drug, are identified in independent populations. These sweeps are further supported by signals of diversifying selection enriched in genes involved in response to drugs and other anthelmintic-associated biological functions. We also identify some candidate genes that may play a role in ivermectin resistance. Finally, genetic signatures of climate-driven adaptation are described, revealing a gene acting as an epigenetic regulator and components of the dauer pathway. These results begin to define genetic adaptation to climate in a parasitic nematode.

Based on single worm whole genome sequencing, the authors here characterise the global evolution of the gastrointestinal parasite Haemonchus contortus and identify genes that play a role in drug resistance as well as climate-driven adaptations involving an epigenetic regulator.

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.

Imported onchocerciasis in migrants and travelers

PURPOSE OF REVIEW

With increasing international travel and mass global population migration, clinicians in nonendemic countries must be familiar with imported neglected tropical diseases including onchocerciasis, which is commonly known as 'river blindness'.

RECENT FINDINGS

Imported onchocerciasis manifests differently in travelers compared with migrants from endemic areas and is likely underdiagnosed in both groups. Recent clinical studies confirm that eosinophilia is not a sensitive marker for Onchocerca volvulus, with one-third of patients having a normal eosinophil count. Novel diagnostics measuring antibodies to multiple recombinant O. volvulus antigens maintain a high sensitivity while improving specificity compared with conventional pan-filarial serologic testing. A 6-week course of doxycycline has macrofilaricidal activity through Wolbachia depletion and may be useful in nonendemic areas in addition to standard serial ivermectin.

SUMMARY

Recent studies characterizing distinct clinical presentations in travelers and migrants may enable clinicians to better recognize imported onchocerciasis. Although novel diagnostics have improved specificity, most remain restricted to tropical disease reference laboratories and to date there is no marker of cure. Prolonged doxycycline treatment may reduce the need for serial ivermectin, though more potent short-course macrofilaricidal drugs are being developed.

Multi-locus sequence typing of Ixodes ricinus and its symbiont Candidatus Midichloria mitochondrii across Europe reveals evidence of local co-cladogenesis in Scotland

Ticks have relatively complex microbiomes, but only a small proportion of the bacterial symbionts recorded from ticks are vertically transmitted. Moreover, co-cladogenesis between ticks and their symbionts, indicating an intimate relationship over evolutionary history driven by a mutualistic association, is the exception rather than the rule. One of the most widespread tick symbionts is Candidatus Midichloria, which has been detected in all of the major tick genera of medical and veterinary importance. In some species of Ixodes, such as the sheep tick Ixodes ricinus (infected with Candidatus Midichloria mitochondrii), the symbiont is fixed in wild adult female ticks, suggesting an obligate mutualism. However, almost no information is available on genetic variation in Candidatus M. mitochondrii or possible co-cladogenesis with its host across its geographic range. Here, we report the first survey of Candidatus M. mitochondrii in I. ricinus in Great Britain and a multi-locus sequence typing (MLST) analysis of tick and symbiont between British ticks and those collected in continental Europe. We show that while the prevalence of the symbiont in nymphs collected in England is similar to that reported from the continent, a higher prevalence in nymphs and adult males is apparent in Wales. In general, Candidatus M. mitochondrii exhibits very low levels of sequence diversity, although a consistent signal of host-symbiont coevolution was apparent in Scotland. Moreover, the tick MLST scheme revealed that Scottish specimens form a clade that is partially separated from other British ticks, with almost no contribution of continental sequence types in this north-westerly border of the tick's natural range. The low diversity of Candidatus M. mitochondrii, in contrast with previously reported high rates of polymorphism in I. ricinus mitogenomes, suggests that the symbiont may have swept across Europe recently via a horizontal, rather than vertical, transmission route.

Multi-locus phylogenetics of the Midichloria endosymbionts reveals variable specificity of association with ticks

Candidatus Midichloria mitochondrii is a maternally inherited bacterium of ticks with a unique intra-mitochondrial lifestyle. Here, we investigate on the evolutionary history of these associations and the degree of Midichloria-tick specificity. While previous surveys used the 16S rRNA gene as an exclusive molecular marker, we rather developed a multi-locus typing method based on four more variable housekeeping genes (groEL, rpoB, dnaK and ftsZ) and on one flagellum gene (fliC) present in Midichloria genomes. Using this method, multi-locus phylogenetic analyses revealed the structuring of a wide Midichloria genetic diversity into three distinct lineages associated with ticks. Overall, two distinct evolutionary strategies are obvious depending on lineage: two Midichloria lineages are generalists with infections acquired through horizontal transfers between distantly related tick species but one other Midichloria lineage rather show a high specificity degree to the Ixodes tick genus. This pattern suggests a capacity of certain Midichloria strains to maintain infections in only limited range of related tick species. These different infection strategies of Midichloria highlight an unexpected variability in their dependency to their tick hosts. We further conjecture that this pattern is also likely to indicate variability in their effects on ticks.

A Multi-Omics Database for Parasitic Nematodes and Trematodes

Helminth.net ( www.helminth.net ) is a web-based resource that was launched in 2000 as simply " Nematode.net " to host and investigate gene sequences from nematode genomes. Over the years it has evolved to become the moniker for a collection of databases: Nematode.net and Trematode.net . These databases host information for 73 nematode (roundworms) and 17 trematode (flatworms) species and serve as backbone for a number of tools that allow users to query slices of the data for multifactorial combinations of species-omics properties. Recent focus has been on inclusion of gene and protein expression data, population genomics and cross-kingdom interactions (metagenomics datasets). This chapter describes the website, the available tools and some of the new features.

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.

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.

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.

Pre-control relationship of onchocercal skin disease with onchocercal infection in Guinea Savanna, Northern Nigeria

Background

Onchocerca volvulus infection can result in blindness, itching and skin lesions. Previous research concentrated on blindness.

Methods

A clinical classification system of the cutaneous changes in onchocerciasis was used for the first time in this study within the context of an early ivermectin drug trial in the savanna region of Kaduna State, northern Nigeria. Skin examinations were performed in 6,790 individuals aged 5+ years in endemic communities and 1,343 individuals in nonendemic communities.

Results / Discussion

There was increased risk for all forms of onchocercal skin disease in endemic communities with the most common finding being the presence of nodules (1,438 individuals, 21.2%), followed by atrophy (367, 6.1% of those < 50 years), acute papular onchodermatitis, APOD (233, 3.4%), depigmentation (216, 3.2%) and chronic papular onchodermatitis, CPOD (155, 2.3%). A further 645 individuals (9.5%) complained of pruritus but had completely normal skin. APOD was more common in males whereas atrophy, hanging groin and nodules were more common in females. After controlling for age and sex, microfilarial positivity was a risk factor for CPOD, depigmentation, hanging groin and nodules (OR 1.54, p = 0.046; OR 2.29, p = 0.002; OR 2.18, p = 0.002 and OR 3.80, p <0.001 respectively). Comparable results were found using presence of nodules as the marker for infection. Microfilarial load showed similar, though weaker, results. A total of 2621(38.6%) endemic residents had itching with normal skin, or had one or more types of onchocercal skin disease including nodules, which may be considered as a composite index of the overall prevalence of onchocercal skin disease.

Conclusion

Significant levels of onchocercal skin disease were documented in this savanna area, which subsequently resulted in a reassessment of the true burden of skin disease in onchocerciasis. This paper represents the first detailed report of the association of onchocercal skin disease with markers for onchocercal infection.

Onchocerciasis is a tropical parasitic infection caused by the nematode worm Onchocerca volvulus. The disease mainly occurs across tropical Africa and infection can result in blindness, debilitating itching and a variety of skin changes. Initial research concentrated mainly on the problem of blindness. A number of studies on onchocercal skin disease were performed but were difficult to interpret and compare because of the use of inconsistent terminology. Within the setting of one of the early trials of ivermectin in a savanna area of northern Nigeria, where there were known high rates of onchocercal blindness, we used a novel clinical classification of the skin changes in onchocerciasis. We identified significant levels of itching and various forms of onchocercal skin disease within these endemic communities. A positive skin-snip result proved to be a significant risk factor for the presence of chronic papular onchodermatitis (CPOD), depigmentation, hanging groin and onchocercal nodules. Comparable results were found when the presence of nodules was used as the marker for infection and similar, though weaker odds ratios were found with microfilarial load per se. The findings triggered a reassessment of the true burden of skin disease in onchocerciasis. It is the first detailed report of the association between onchocercal skin disease and markers of infection.

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.