Modelling exposure heterogeneity and density dependence in onchocerciasis using a novel individual-based transmission model, EPIONCHO-IBM: Implications for elimination and data needs.

Background Density dependence in helminth establishment and heterogeneity in exposure to infection are known to drive resilience to interventions based on mass drug administration (MDA). However, the interaction between these processes is poorly understood. We developed a novel individual-based model for onchocerciasis transmission, EPIONCHO-IBM, which accounts for both processes. We fit the model to pre-intervention epidemiological data and explore parasite dynamics during MDA with ivermectin. Methodology/Principal findings Density dependence and heterogeneity in exposure to blackfly (vector) bites were estimated by fitting the model to matched pre-intervention microfilarial prevalence, microfilarial intensity and vector biting rate data from savannah areas of Cameroon and Cote d'Ivoire/Burkina Faso using Latin hypercube sampling. Transmission dynamics during 25 years of annual and biannual ivermectin MDA were investigated. Density dependence in parasite establishment within humans was estimated for different levels of (fixed) exposure heterogeneity to understand how parametric uncertainty may influence treatment dynamics. Stronger overdispersion in exposure to blackfly bites results in the estimation of stronger density-dependent parasite establishment within humans, consequently increasing resilience to MDA. For all levels of exposure heterogeneity tested, the model predicts a departure from the functional forms for density dependence assumed in the deterministic version of the model. Conclusions/Significance This is the first, stochastic model of onchocerciasis, that accounts for and estimates densitydependent parasite establishment in humans alongside exposure heterogeneity. Capturing the interaction between these processes is fundamental to our understanding of resilience to MDA interventions. Given that uncertainty in these processes results in very different treatment dynamics, collecting data on exposure heterogeneity would be essential for improving model predictions during MDA. We discuss possible ways in which such data may be collected as well as the importance of better understanding the effects of immunological responses on establishing parasites prior to and during ivermectin treatment.