A research agenda for helminth diseases of humans: modelling for control and elimination

Spatial transferability of an agent-based model to simulate Taenia solium control interventions

BACKGROUND

Models can be used to study and predict the impact of interventions aimed at controlling the spread of infectious agents, such as Taenia solium, a zoonotic parasite whose larval stage causes epilepsy and economic loss in many rural areas of the developing nations. To enhance the credibility of model estimates, calibration against observed data is necessary. However, this process may lead to a paradoxical dependence of model parameters on location-specific data, thus limiting the model's geographic transferability.

METHODS

In this study, we adopted a non-local model calibration approach to assess whether it can improve the spatial transferability of CystiAgent, our agent-based model of local-scale T. solium transmission. The calibration dataset for CystiAgent consisted of cross-sectional data on human taeniasis, pig cysticercosis and pig serology collected in eight villages in Northwest Peru. After calibration, the model was transferred to a second group of 21 destination villages in the same area without recalibrating its parameters. Model outputs were compared to pig serology data collected over a period of 2 years in the destination villages during a trial of T. solium control interventions, based on mass and spatially targeted human and pig treatments.

RESULTS

Considering the uncertainties associated with empirical data, the model produced simulated pre-intervention pig seroprevalences that were successfully validated against data collected in 81% of destination villages. Furthermore, the model outputs were able to reproduce validated pig seroincidence values in 76% of destination villages when compared to the data obtained after the interventions. The results demonstrate that the CystiAgent model, when calibrated using a non-local approach, can be successfully transferred without requiring additional calibration.

CONCLUSIONS

This feature allows the model to simulate both baseline pre-intervention transmission conditions and the outcomes of control interventions across villages that form geographically homogeneous regions, providing a basis for developing large-scale models representing T. solium transmission at a regional level.

Diagnosis of helminths depends on worm fecundity and the distribution of parasites within hosts

Helminth transmission and morbidity are dependent on the number of mature parasites within a host; however, observing adult worms is impossible for many natural infections. An outstanding challenge is therefore relating routine diagnostics, such as faecal egg counts, to the underlying worm burden. This relationship is complicated by density-dependent fecundity (egg output per worm reduces due to crowding at high burdens) and the skewed distribution of parasites (majority of helminths aggregated in a small fraction of hosts). We address these questions for the carcinogenic liver fluke Opisthorchis viverrini, which infects approximately 10 million people across Southeast Asia, by analysing five epidemiological surveys (n = 641) where adult flukes were recovered. Using a mechanistic model, we show that parasite fecundity varies between populations, with surveys from Thailand and Laos demonstrating distinct patterns of egg output and density-dependence. As the probability of observing faecal eggs increases with the number of mature parasites within a host, we quantify diagnostic sensitivity as a function of the worm burden and find that greater than 50% of cases are misdiagnosed as false negative in communities close to elimination. Finally, we demonstrate that the relationship between observed prevalence from routine diagnostics and true prevalence is nonlinear and strongly influenced by parasite aggregation.

On Spatiotemporal Overdispersion and Macroparasite Accumulation in Hosts Leading to Aggregation: A Quantitative Framework

In many host-parasite systems, overdispersion in the distribution of macroparasites leads to parasite aggregation in the host population. This overdispersed distribution is often characterized by the negative binomial or by the power law. The aggregation is shown by a clustering of parasites in certain hosts, while other hosts have few or none. Here, I present a theory behind the overdispersion in complex spatiotemporal systems as well as a computational analysis for tracking the behavior of transmissible diseases with this kind of dynamics. I present a framework where heterogeneity and complexity in host-parasite systems are related to aggregation. I discuss the problem of focusing only on the average parasite burden without observing the risk posed by the associated variance; the consequences of under- or overestimation of disease transmission in a heterogenous system and environment; the advantage of including the network of social interaction in epidemiological modeling; and the implication of overdispersion in the management of health systems during outbreaks.

Sensitivity Analysis of a Transmission Interruption Model for the Soil-Transmitted Helminth Infections in Kenya

As the world rallies toward the endgame of soil-transmitted helminths (STH) elimination by the year 2030, there is a need for efficient and robust mathematical models that would enable STH programme managers to target the scarce resources and interventions, increase treatment coverage among specific sub-groups of the population, and develop reliable surveillance systems that meet sensitivity and specificity requirements for the endgame of STH elimination. However, the considerable complexities often associated with STH-transmission models underpin the need for specifying a large number of parameters and inputs, which are often available with considerable degree of uncertainty. Additionally, the model may behave counter-intuitive especially when there are non-linearities in multiple input-output relationships. In this study, we performed a global sensitivity analysis (GSA), based on a variance decomposition method: extended Fourier Amplitude Sensitivity Test (eFAST), to a recently developed STH-transmission model in Kenya (an STH endemic country) to; (1) robustly compute sensitivity index (SI) for each parameter, (2) rank the parameters in order of their importance (from most to least influential), and (3) quantify the influence of each parameter, singly and cumulatively, on the model output. The sensitivity analysis (SA) results demonstrated that the model outcome (STH worm burden elimination in the human host) was significantly sensitive to some key parameter groupings: combined effect of improved water source and sanitation (ϕ), rounds of treatment offered (τ), efficacy of the drug used during treatment (h), proportion of the adult population treated (g_a: akin to community-wide treatment), mortality rate of the mature worms in the human host (μ), and the strength of the -dependence of worm egg production (γ). For STH control programmes to effectively reach the endgame (STH elimination in the entire community), these key parameter groupings need to be targeted since together they contribute to a strategic public health intervention.

Predicting future community-level ocular Chlamydia trachomatis infection prevalence using serological, clinical, molecular, and geospatial data

Trachoma is an infectious disease characterized by repeated exposures to Chlamydia trachomatis (Ct) that may ultimately lead to blindness. Efficient identification of communities with high infection burden could help target more intensive control efforts. We hypothesized that IgG seroprevalence in combination with geospatial layers, machine learning, and model-based geostatistics would be able to accurately predict future community-level ocular Ct infections detected by PCR. We used measurements from 40 communities in the hyperendemic Amhara region of Ethiopia to assess this hypothesis. Median Ct infection prevalence among children 0–5 years old increased from 6% at enrollment, in the context of recent mass drug administration (MDA), to 29% by month 36, following three years without MDA. At baseline, correlation between seroprevalence and Ct infection was stronger among children 0–5 years old (ρ = 0.77) than children 6–9 years old (ρ = 0.48), and stronger than the correlation between active trachoma and Ct infection (0-5y ρ = 0.56; 6-9y ρ = 0.40). Seroprevalence was the strongest concurrent predictor of infection prevalence at month 36 among children 0–5 years old (cross-validated R2 = 0.75, 95% CI: 0.58–0.85), though predictive performance declined substantially with increasing temporal lag between predictor and outcome measurements. Geospatial variables, a spatial Gaussian process, and stacked ensemble machine learning did not meaningfully improve predictions. Serological markers among children 0–5 years old may be an objective tool for identifying communities with high levels of ocular Ct infections, but accurate, future prediction in the context of changing transmission remains an open challenge.

Modeling the Interruption of the Transmission of Soil-Transmitted Helminths Infections in Kenya: Modeling Deworming, Water, and Sanitation Impacts

Kenya, just like other countries with endemic soil-transmitted helminths (STH), has conducted regular mass drug administration (MDA) program for the last 5 years among school aged children as a way to reduce STH infections burden in the country. However, the point of interruption of transmission of these infections still remains unclear. In this study, we developed and analyzed an age structured mathematical model to predict the elimination period (i.e., time taken to interrupt STH transmission) of these infections in Kenya. The study utilized a deterministic age structured model of the STH population dynamics under a regular treatment program. The model was applied to three main age groups: pre-school age children (2-4 years), school age children (5-14 years), and adult populations (≥15 years) and compared the impact of two interventions on worm burden and elimination period. The model-simulated results were compared with the 5 year field data from the Kenyan deworming program for all the three types of STH (Ascaris lumbricoides, Trichuris trichiura, and hookworm). The model demonstrated that the reduction of worm burden and elimination period depended heavily on four parameter groups; drug efficacy, number of treatment rounds, MDA and water, sanitation and hygiene (WASH) coverage. The analysis showed that for STH infections to be eliminated using MDA alone in a short time period, 3-monthly MDA plan is desired. However, complementation of MDA with WASH at an optimal (95%) coverage level was most effective. These results are important to the Kenyan STH control program as it will guide the recently launched Breaking Transmission Strategy.

Epidemiological Impact of SARS-CoV-2 Vaccination: Mathematical Modeling Analyses

This study aims to inform SARS-CoV-2 vaccine development/licensure/decision-making/implementation, using mathematical modeling, by determining key preferred vaccine product characteristics and associated population-level impacts of a vaccine eliciting long-term protection. A prophylactic vaccine with efficacy against acquisition (VE_S) ≥70% can eliminate the infection. A vaccine with VE_S <70% may still control the infection if it reduces infectiousness or infection duration among those vaccinated who acquire the infection, if it is supplemented with <20% reduction in contact rate, or if it is complemented with herd-immunity. At VE_S of 50%, the number of vaccinated persons needed to avert one infection is 2.4, and the number is 25.5 to avert one severe disease case, 33.2 to avert one critical disease case, and 65.1 to avert one death. The probability of a major outbreak is zero at VE_S ≥70% regardless of the number of virus introductions. However, an increase in social contact rate among those vaccinated (behavior compensation) can undermine vaccine impact. In addition to the reduction in infection acquisition, developers should assess the natural history and disease progression outcomes when evaluating vaccine impact.

Epidemiological Impact of SARS-CoV-2 Vaccination: Mathematical Modeling Analyses

This study aims to inform SARS-CoV-2 vaccine development/licensure/decision-making/implementation, using mathematical modeling, by determining key preferred vaccine product characteristics and associated population-level impacts of a vaccine eliciting long-term protection. A prophylactic vaccine with efficacy against acquisition (VES) ≥70% can eliminate the infection. A vaccine with VES <70% may still control the infection if it reduces infectiousness or infection duration among those vaccinated who acquire the infection, if it is supplemented with <20% reduction in contact rate, or if it is complemented with herd-immunity. At VES of 50%, the number of vaccinated persons needed to avert one infection is 2.4, and the number is 25.5 to avert one severe disease case, 33.2 to avert one critical disease case, and 65.1 to avert one death. The probability of a major outbreak is zero at VES ≥70% regardless of the number of virus introductions. However, an increase in social contact rate among those vaccinated (behavior compensation) can undermine vaccine impact. In addition to the reduction in infection acquisition, developers should assess the natural history and disease progression outcomes when evaluating vaccine impact.

Biting rates and relative abundance of Simulium flies under different climatic conditions in an onchocerciasis endemic community in Ghana

BACKGROUND

Knowledge of the relative abundance and biting rates of riverine blackflies (vectors of onchocerciasis) is essential as these entomological indices affect transmission of the disease. However, transmission patterns vary from one ecological zone to another and this may be due to differences in species of blackfly vectors and the climatic conditions in the area. This study investigated the effects of climate variability on the relative abundance and biting rates of blackflies in the Tanfiano community (Nkoranza North District, Bono East Region, Ghana). Such information will help to direct policy on effective timing of the annual mass drug administration of ivermectin in the area.

METHODS

The study employed human landing collections and locally built Esperanza window traps to collect blackflies from March 2018 to February 2019. The relative abundance and biting rates of the Simulium vectors as well as the monthly climatic conditions of the study area were monitored. Correlation analysis and Poisson regression were used to establish the relationships between the variables.

RESULTS

The relative abundance and biting rates of the Simulium vectors were highest in the drier months of March, April and August, characterized by high temperatures, low humidity, longer hours of sunshine and stronger winds. The rainy months of May, June and July, characterized by low temperatures, high humidity, few hours of sunshine and weaker winds, had relatively low blackfly abundance and biting activity. Correlation analysis showed that only temperature was significantly, positively correlated with the relative abundance of blackflies (r = 0.617, n = 12, P = 0.033) and monthly biting rates (r = 0.612, n = 12, P = 0.034). A model to predict relative abundance and monthly biting rates using climatological variables was developed.

CONCLUSIONS

This study demonstrated that Simulium species in the study area preferred higher temperature, lower humidity and rainfall, more hours of sunshine and relatively stronger winds for survival. It is thus recommended that for the study district and others with similar climatological characteristics, mass drug administration of ivermectin should take place in April and September when the abundance of vectors has begun to decline after peaking.

Assessing the impact of aggregating disease stage data in model predictions of human African trypanosomiasis transmission and control activities in Bandundu province (DRC)

Since the turn of the century, the global community has made great progress towards the elimination of gambiense human African trypanosomiasis (HAT). Elimination programs, primarily relying on screening and treatment campaigns, have also created a rich database of HAT epidemiology. Mathematical models calibrated with these data can help to fill remaining gaps in our understanding of HAT transmission dynamics, including key operational research questions such as whether integrating vector control with current intervention strategies is needed to achieve HAT elimination. Here we explore, via an ensemble of models and simulation studies, how including or not disease stage data, or using more updated data sets affect model predictions of future control strategies.

Currently Available Monitoring and Surveillance Systems for Taenia spp., Echinococcus spp., Schistosoma spp., and Soil-Transmitted Helminths at the Control/Elimination Stage: A Systematic Review

An increasing global focus on neglected tropical diseases (NTDs) has resulted in the set up of numerous control and elimination activities worldwide. This is partly true for Taenia solium taeniasis/cysticercosis, the most important foodborne parasitic infection. Despite substantial progress, adequate monitoring and surveillance (M&S) are required to sustain a status of control/elimination. This is often lacking, especially for T. solium. Therefore, the objective was to conduct a systematic literature review of the currently available M&S systems at the control/elimination stage of the four top-ranked helminth NTDs. Specifically, Taenia spp., Echinococcus spp., Schistosoma spp., and soil-transmitted helminths (STHs) were considered to determine if there are any similarities between their M&S systems and whether certain approaches can be adopted from each other. The systematic review demonstrated that rigorous M&S systems have been designed for the control/elimination stage of both STHs and schistosomiasis, particularly in China. On the other hand, a concept of M&S for Taenia spp. and Echinococcus spp. has not been fully developed yet, due to a lack of epidemiological data and the fact that many endemic countries are far away from reaching control/elimination. Moreover, accurate diagnostic tools for all four diseases are still imperfect, which complicates proper M&S. Finally, there is an urgent need to develop and harmonize/standardize M&S activities in order to reliably determine and compare the epidemiological situation worldwide.

Zoonotic transmission of intestinal helminths in southeast Asia: Implications for control and elimination

Intestinal helminths are extremely widespread and highly prevalent infections of humans, particularly in rural and poor urban areas of low and middle-income countries. These parasites have chronic and often insidious effects on human health and child development including abdominal problems, anaemia, stunting and wasting. Certain animals play a fundamental role in the transmission of many intestinal helminths to humans. However, the contribution of zoonotic transmission to the overall burden of human intestinal helminth infection and the relative importance of different animal reservoirs remains incomplete. Moreover, control programmes and transmission models for intestinal helminths often do not consider the role of zoonotic reservoirs of infection. Such reservoirs will become increasingly important as control is scaled up and there is a move towards interruption and even elimination of parasite transmission. With a focus on southeast Asia, and the Philippines in particular, this review summarises the major zoonotic intestinal helminths, risk factors for infection and highlights knowledge gaps related to their epidemiology and transmission. Various methodologies are discussed, including parasite genomics, mathematical modelling and socio-economic analysis, that could be employed to improve understanding of intestinal helminth spread, reservoir attribution and the burden associated with infection, as well as assess effectiveness of interventions. For sustainable control and ultimately elimination of intestinal helminths, there is a need to move beyond scheduled mass deworming and to consider animal and environmental reservoirs. A One Health approach to control of intestinal helminths is proposed, integrating interventions targeting humans, animals and the environment, including improved access to water, hygiene and sanitation. This will require coordination and collaboration across different sectors to achieve best health outcomes for all.

Mapping Schistosoma mansoni endemicity in Rwanda: a critical assessment of geographical disparities arising from circulating cathodic antigen versus Kato-Katz diagnostics

BACKGROUND

Schistosomiasis is a neglected tropical disease caused by Schistosoma parasites. Intervention relies on identifying high-risk regions, yet rapid Schistosoma diagnostics (Kato-Katz stool assays (KK) and circulating cathodic antigen urine assays (CCA)) yield different prevalence estimates. We mapped S. mansoni prevalence and delineated at-risk regions using a survey of schoolchildren in Rwanda, where S. mansoni is an endemic parasite. We asked if different diagnostics resulted in disparities in projected infection risk.

METHODS

Infection data was obtained from a 2014 Rwandan school-based survey that used KK and CCA diagnostics. Across 386 schools screened by CCA (N = 19,217). To allow for uncertainty when interpreting ambiguous CCA trace readings, which accounted for 28.8% of total test results, we generated two presence-absence datasets: CCA trace as positive and CCA trace as negative. Samples (N = 9,175) from 185 schools were also screened by KK. We included land surface temperature (LST) and the Normalized Difference Vegetation and Normalized Difference Water Indices (NDVI, NDWI) as predictors in geostatistical regressions.

FINDINGS

Across 8,647 children tested by both methods, prevalence was 35.93% for CCA trace as positive, 7.21% for CCA trace as negative and 1.95% for KK. LST was identified as a risk factor using KK, whereas NDVI was a risk factor for CCA models. Models predicted high endemicity in Northern and Western regions of Rwanda, though the CCA trace as positive model identified additional high-risk areas that were overlooked by the other methods. Estimates of current burden for children at highest risk (boys aged 5-9 years) varied by an order of magnitude, with 671,856 boys projected to be infected by CCA trace as positive and only 60,453 projected by CCA trace as negative results.

CONCLUSIONS

Our findings show that people in Rwanda's Northern, Western and capital regions are at high risk of S. mansoni infection. However, variation in identification of environmental risk factors and delineation of at-risk regions using different diagnostics likely provides confusing messages to disease intervention managers. Further research and statistical analyses, such as latent class analysis, can be used to improve CCA result classification and assess its use in guiding treatment regimes.

Mapping Soil-Transmitted Helminth Parasite Infection in Rwanda: Estimating Endemicity and Identifying At-Risk Populations

Soil-transmitted helminth (STH) infections are globally distributed intestinal parasite infections caused by Ascaris lumbricoides, Trichuris trichiura, and hookworms (Ancylostoma duodenale and Necator americanus). STH infection constitutes a major public health threat, with heavy burdens observed in many of the world’s tropical and subtropical regions. Mass drug administration and sanitation improvements can drastically reduce STH prevalence and associated morbidity. However, identifying targeted areas in need of treatment is hampered by a lack of knowledge on geographical and population-level risk factors. In this study, we applied Bayesian geostatistical modelling to data from a national school-based STH infection survey in Rwanda to (1) identify ecological and population-level risk factors and (2) provide comprehensive precision maps of infection burdens. Our results indicated that STH infections were heterogeneously distributed across the country and showed signatures of spatial clustering, though the magnitude of clustering varied among parasites. The highest rates of endemic clustering were attributed to A. lumbricoides infection. Concordant infection patterns among the three parasite groups highlighted populations currently most at-risk of morbidity. Population-dense areas in the Western and North-Western regions of Rwanda represent areas that have continued to exhibit high STH burden across two surveys and are likely in need of targeted interventions. Our maps support the need for an updated evaluation of STH endemicity in western Rwanda to evaluate progress in MDA efforts and identify communities that need further local interventions to further reduce morbidity caused by STH infections.

Strategies for tackling Taenia solium taeniosis/cysticercosis: A systematic review and comparison of transmission models, including an assessment of the wider Taeniidae family transmission models

Background

The cestode Taenia solium causes the neglected (zoonotic) tropical disease cysticercosis, a leading cause of preventable epilepsy in endemic low and middle-income countries. Transmission models can inform current scaling-up of control efforts by helping to identify, validate and optimise control and elimination strategies as proposed by the World Health Organization (WHO).

Methodology/Principal findings

A systematic literature search was conducted using the PRISMA approach to identify and compare existing T. solium transmission models, and related Taeniidae infection transmission models. In total, 28 modelling papers were identified, of which four modelled T. solium exclusively. Different modelling approaches for T. solium included deterministic, Reed-Frost, individual-based, decision-tree, and conceptual frameworks. Simulated interventions across models agreed on the importance of coverage for impactful effectiveness to be achieved.

Other Taeniidae infection transmission models comprised force-of-infection (FoI), population-based (mainly Echinococcus granulosus) and individual-based (mainly E. multilocularis) modelling approaches. Spatial structure has also been incorporated (E. multilocularis and Taenia ovis) in recognition of spatial aggregation of parasite eggs in the environment and movement of wild animal host populations.

Conclusions/Significance

Gaps identified from examining the wider Taeniidae family models highlighted the potential role of FoI modelling to inform model parameterisation, as well as the need for spatial modelling and suitable structuring of interventions as key areas for future T. solium model development. We conclude that working with field partners to address data gaps and conducting cross-model validation with baseline and longitudinal data will be critical to building consensus-led and epidemiological setting-appropriate intervention strategies to help fulfil the WHO targets.

Taenia solium infection in humans (taeniosis and neurocysticercosis) and pigs (cysticercosis) presents a significant global public health and economic challenge. The World Health Organization has called for validated strategies and wider consensus on which strategies are suitable for different epidemiological settings to support successful T. solium control and elimination efforts. Transmission models can be used to inform these strategies. Therefore, a modelling review was undertaken to assess the current state and gaps relating to T. solium epidemiological modelling. The literature surrounding models for other Taeniidae family infections was also considered, identifying approaches to aid further development of existing T. solium models. A variety of different modelling approaches have been used for T. solium including differences in structural and parametric assumptions associated with T. solium transmission biology. Despite these differences, all models agreed on the importance of coverage on intervention effectiveness. Other Taeniidae family models highlighted the need for incorporating spatial structure when necessary to capture aggregation of transmission stages in the environment and movement of animal hosts.

Modelling the elimination of river blindness using long-term epidemiological and programmatic data from Mali and Senegal

•

Onchocerciasis is earmarked for elimination in some African countries by 2020/2025.

•

15+ years of ivermectin treatment drove infection prevalence to zero in areas of Mali & Senegal.

•

Data-driven model projections are used to evaluate the risk of infection resurgence.

•

Latent infections can initiate slow resurgence in communities with high transmission propensity.

•

Highly sensitive and long-term surveillance will be necessary to verify elimination.

The onchocerciasis transmission models EPIONCHO and ONCHOSIM have been independently developed and used to explore the feasibility of eliminating onchocerciasis from Africa with mass (annual or biannual) distribution of ivermectin within the timeframes proposed by the World Health Organization (WHO) and endorsed by the 2012 London Declaration on Neglected Tropical Diseases (i.e. by 2020/2025). Based on the findings of our previous model comparison, we implemented technical refinements and tested the projections of EPIONCHO and ONCHOSIM against long-term epidemiological data from two West African transmission foci in Mali and Senegal where the observed prevalence of infection was brought to zero circa 2007–2009 after 15–17 years of mass ivermectin treatment. We simulated these interventions using programmatic information on the frequency and coverage of mass treatments and trained the model projections using longitudinal parasitological data from 27 communities, evaluating the projected outcome of elimination (local parasite extinction) or resurgence. We found that EPIONCHO and ONCHOSIM captured adequately the epidemiological trends during mass treatment but that resurgence, while never predicted by ONCHOSIM, was predicted by EPIONCHO in some communities with the highest (inferred) vector biting rates and associated pre-intervention endemicities. Resurgence can be extremely protracted such that low (microfilarial) prevalence between 1% and 5% can be maintained for 3–5 years before manifesting more prominently. We highlight that post-treatment and post-elimination surveillance protocols must be implemented for long enough and with high enough sensitivity to detect possible residual latent infections potentially indicative of resurgence. We also discuss uncertainty and differences between EPIONCHO and ONCHOSIM projections, the potential importance of vector control in high-transmission settings as a complementary intervention strategy, and the short remaining timeline for African countries to be ready to stop treatment safely and begin surveillance in order to meet the impending 2020/2025 elimination targets.

Review of mathematical models of HSV-2 vaccination: Implications for vaccine development

Development of a vaccine against herpes simplex virus type 2 (HSV-2), a life-long sexually-transmitted infection (STI), would be a major step forward in improving global sexual and reproductive health. In this review, we identified published literature of dynamic mathematical models assessing the impact of either prophylactic or therapeutic HSV-2 vaccination at the population level. We compared each study's model structure and assumptions as well as predicted vaccination impact. We examined possible causes of heterogeneity across model predictions, key gaps, and the implications of these findings for future modelling efforts. Only eight modelling studies have assessed the potential public health impact of HSV-2 vaccination, with the majority focusing on impact of prophylactic vaccines. The studies showed that even an imperfect prophylactic HSV-2 vaccine could have an important public health impact on HSV-2 incidence, and could also impact HIV indirectly in high HIV prevalence settings. Therapeutic vaccines also may provide public health benefits, though they have been explored less extensively. However, there was substantial variation in predicted population-level impact for both types of vaccine, reflecting differences in assumptions between model scenarios. Importantly, many models did not account for heterogeneity in infection rates such as by age, sex and sexual activity. Future modelling work to inform decisions on HSV vaccine development and implementation should consider cost-effectiveness, account for additional HSV-2 sequelae such as neonatal transmission, and model greater heterogeneity in infection rates between individuals, more realistic vaccine deployment, and more thorough sensitivity and uncertainty analyses.

Opportunities and challenges for modelling epidemiological and evolutionary dynamics in a multihost, multiparasite system: Zoonotic hybrid schistosomiasis in West Africa

Multihost multiparasite systems are evolutionarily and ecologically dynamic, which presents substantial trans-disciplinary challenges for elucidating their epidemiology and designing appropriate control. Evidence for hybridizations and introgressions between parasite species is gathering, in part in line with improvements in molecular diagnostics and genome sequencing. One major system where this is becoming apparent is within the Genus Schistosoma, where schistosomiasis represents a disease of considerable medical and veterinary importance, the greatest burden of which occurs in sub-Saharan Africa. Interspecific hybridizations and introgressions bring an increased level of complexity over and above that already inherent within multihost, multiparasite systems, also representing an additional source of genetic variation that can drive evolution. This has the potential for profound implications for the control of parasitic diseases, including, but not exclusive to, widening host range, increased transmission potential and altered responses to drug therapy. Here, we present the challenging case example of haematobium group Schistosoma spp. hybrids in West Africa, a system involving multiple interacting parasites and multiple definitive hosts, in a region where zoonotic reservoirs of schistosomiasis were not previously considered to be of importance. We consider how existing mathematical model frameworks for schistosome transmission could be expanded and adapted to zoonotic hybrid systems, exploring how such model frameworks can utilize molecular and epidemiological data, as well as the complexities and challenges this presents. We also highlight the opportunities and value such mathematical models could bring to this and a range of similar multihost, multi and cross-hybridizing parasites systems in our changing world.

Geographic variation in sexual behavior can explain geospatial heterogeneity in the severity of the HIV epidemic in Malawi

BACKGROUND

In sub-Saharan Africa, where ~ 25 million individuals are infected with HIV and transmission is predominantly heterosexual, there is substantial geographic variation in the severity of epidemics. This variation has yet to be explained. Here, we propose that it is due to geographic variation in the size of the high-risk group (HRG): the group with a high number of sex partners. We test our hypothesis by conducting a geospatial analysis of data from Malawi, where ~ 13% of women and ~ 8% of men are infected with HIV.

METHODS

We used georeferenced HIV testing and behavioral data from ~ 14,000 participants of a nationally representative population-level survey: the 2010 Malawi Demographic and Health Survey (MDHS). We constructed gender-stratified epidemic surface prevalence (ESP) maps by spatially smoothing and interpolating the HIV testing data. We used the behavioral data to construct gender-stratified risk maps that reveal geographic variation in the size of the HRG. We tested our hypothesis by fitting gender-stratified spatial error regression (SER) models to the MDHS data.

RESULTS

The ESP maps show considerable geographic variation in prevalence: 1-29% (women), 1-20% (men). Risk maps reveal substantial geographic variation in the size of the HRG: 0-40% (women), 16-58% (men). Prevalence and the size of the HRG are highest in urban centers. However, the majority of HIV-infected individuals (~75% of women, ~ 80% of men) live in rural areas, as does most of the HRG (~ 80% of women, ~ 85% of men). We identify a significant (P < 0.001) geospatial relationship linking the size of the HRG with prevalence: the greater the size, the higher the prevalence. SER models show HIV prevalence in women is expected to exceed the national average in districts where > 20% of women are in the HRG. Most importantly, the SER models show that geographic variation in the size of the HRG can explain a substantial proportion (73% for women, 67% for men) of the geographic variation in epidemic severity.

CONCLUSIONS

Taken together, our results provide substantial support for our hypothesis. They provide a potential mechanistic explanation for the geographic variation in the severity of the HIV epidemic in Malawi and, potentially, in other countries in sub-Saharan Africa.

Collaborate or Collapse: Capacity Building in Zoonotic and Neglected Tropical Disease Modelling

We share the insights from a successful collaboration in organizing and implementing an international scientific capacity-building workshop in Malaysia titled Mathematical Modelling of Neglected Infectious Diseases: Capacity Building in Southeast Asia. This workshop focused on the delivery of technical know-how and on essential soft skills related to effective grant proposal writing and networking.

Spatiotemporal distribution and population at risk of soil-transmitted helminth infections following an eight-year school-based deworming programme in Burundi, 2007-2014

Background

Investigating the effect of successive annual deworming rounds on the spatiotemporal distribution of infection prevalence and numbers at risk for soil-transmitted helminths (STHs) can help identify communities nearing elimination and those needing further interventions. In this study, we aim to quantify the impact of an 8-year mass drug administration (MDA) programme (from 2007 to 2014) on the spatiotemporal distribution of prevalence of STH infections and to estimate the number of school-aged children infected with STHs in Burundi.

Methods

During annual longitudinal school-based surveys in Burundi between 2007 and 2011, STH infection and anthropometric data for a total of 40,656 children were collected; these data were supplemented with data from a national survey conducted in 2014. Bayesian model based geostatistics (MBG) were used to generate predictive prevalence maps for each STH species and year. The numbers of children at-risk of infection per district between 2008 and 2014 were estimated as the product of the predictive prevalence maps and population density maps.

Results

Overall, the degree of spatial clustering of STH infections decreased between 2008 and 2011; in 2014 the geographical clusters of all STH infections reappeared. The reduction in prevalence was small for Ascaris lumbricoides and Trichuris trichiura in the centre and central north of the country. Our predictive prevalence maps for hookworm indicate a reduction in prevalence along the periphery of the country. The predicted number of children infected with any STH species decreased substantially between 2007 and 2011, but in 2014 there was an increase in the predicted number of children infected with A. lumbricoides and T. trichiura. In 2014, the districts with the highest predicted number of children infected with A. lumbricoides, T. trichiura and hookworms were Kibuye district (n = 128,903), Mabayi district (n = 35,302) and Kiremba (n = 87,511), respectively.

Conclusions

While the MDA programme in Burundi resulted in a reduction in STH prevalence, this reduction was spatiotemporally heterogeneous, with some pockets of high prevalence remaining, suggesting that treatment coverage and complementary interventions should be evaluated to improve impact.

Electronic supplementary material

The online version of this article (10.1186/s13071-017-2505-x) contains supplementary material, which is available to authorized users.

Systematic review of studies generating individual participant data on the efficacy of drugs for treating soil-transmitted helminthiases and the case for data-sharing

Background

Preventive chemotherapy and transmission control (PCT) by mass drug administration is the cornerstone of the World Health Organization (WHO)’s policy to control soil-transmitted helminthiases (STHs) caused by Ascaris lumbricoides (roundworm), Trichuris trichiura (whipworm) and hookworm species (Necator americanus and Ancylostama duodenale) which affect over 1 billion people globally. Despite consensus that drug efficacies should be monitored for signs of decline that could jeopardise the effectiveness of PCT, systematic monitoring and evaluation is seldom implemented. Drug trials mostly report aggregate efficacies in groups of participants, but heterogeneities in design complicate classical meta-analyses of these data. Individual participant data (IPD) permit more detailed analysis of drug efficacies, offering increased sensitivity to identify atypical responses potentially caused by emerging drug resistance.

Methodology

We performed a systematic literature review to identify studies concluding after 2000 that collected IPD suitable for estimating drug efficacy against STH. We included studies that administered a variety of anthelmintics with follow ups less than 60 days after treatment. We estimated the number of IPD and extracted cohort- and study-level meta-data.

Principal findings

We estimate that there exist individual data on approximately 35,000 participants from 129 studies conducted in 39 countries, including 34 out of 103 countries where PCT is recommended. We find significant heterogeneity in diagnostic methods, times of outcome assessment, and the reported measure of efficacy. We also quantify cohorts comprising pre-school age children, pregnant women, and co-infected participants, including with HIV.

Conclusions

We argue that establishing a global IPD repository would improve the capacity to monitor and evaluate the efficacy of anthelmintic drugs, respond to changes and safeguard the ongoing effectiveness of PCT. Establishing a fair, transparent data governance policy will be key for the engagement of the global STH community.

Soil-transmitted helminthiases (STHs) caused by roundworm, whipworm or hookworm affect over one billion of the world’s poorest people mostly living in low and middle income countries, exerting a major health and economic toll. These infections are controlled by regular mass drug distribution to affected populations. But with very few alternative medicines, the effectiveness of treatment programmes is vulnerable to the potential emergence of drug resistance. Despite a recent scale-up of mass drug distribution, systematic monitoring and evaluation of the efficacy of treatment is too rarely undertaken and our knowledge of how the drugs are performing is largely based on information from clinical trials. However, the design and reporting of information from these trials is very variable which makes it difficult to form a comprehensive picture of the status and trends in drug efficacy. Here, we present a systematic review of published studies completed since 2000, characterise variation in their design, implementation and reporting and estimate the abundance of individual participant data. We argue that the co-ordinated sharing of these individual data would greatly increase the capacity of the global health community to monitor effectively drug efficacy, to respond accordingly to changes, and thereby to safeguard the effectiveness of STH control.

Experimental and modelling investigations of Opisthorchis viverrini miracidia transmission over time and across temperatures: implications for control

Transmissibility is a significant factor in parasite fitness. The rate and magnitude of parasite transmission affect prevalence and infection intensity in individual hosts and are influenced by environmental factors. In this context, the objectives of this study were: (i) to experimentally assess Opisthorchis viverrini miracidia survival and infectivity over time and across temperatures; and (ii) to combine these experimental results with environmental data to build a key component of a transmission model, identifying seasonal windows of transmission risk in hyper-endemic northeastern Thailand. Five replicates of 50 O. viverrini eggs were randomly distributed and maintained under four temperature conditions (25°C, 30°C, 35°C, 40°C). Microscopic observations were performed on all experimental units over a period of 3months to record miracidia motility and mortality trends. Six infection trials were also conducted to assess infectivity of miracidia over time and across temperatures, using observations of egg hatching success and infection rates. Upon completion of experiments, data were integrated into a transmission model to create a transmission risk index and to simulate seasonal transmission risk. Miracidia survival rate and motility decreased steadily with 50% mortality observed after 2weeks. Hatching and infection success also decreased significantly after 3weeks. Temperatures over 30°C were associated with increased mortality and decreased infectivity. When incorporating local environmental parameters into our model, we observed low transmission risk during the dry season and increasing transmission risk at the onset of the rainy season, culminating with the highest risk in September. We believe that our results provide the first estimates of O. viverrini miracidia survival and transmission potential under variable temperature conditions and suggest that high temperature treatment (>40°C) of fecal waste could be an efficient control strategy.

Assessing the impact of intervention strategies against Taenia solium cysticercosis using the EPICYST transmission model

BACKGROUND

The pork tapeworm, Taenia solium, and associated human infections, taeniasis, cysticercosis and neurocysticercosis, are serious public health problems, especially in developing countries. The World Health Organization (WHO) has set goals for having a validated strategy for control and elimination of T. solium taeniasis/cysticercosis by 2015 and interventions scaled-up in selected countries by 2020. Timely achievement of these internationally-endorsed targets requires that the relative benefits and effectiveness of potential interventions be explored rigorously within a quantitative framework.

METHODS

A deterministic, compartmental transmission model (EPICYST) was developed to capture the dynamics of the taeniasis/cysticercosis disease system in the human and pig hosts. Cysticercosis prevalence in humans, an outcome of high epidemiological and clinical importance, was explicitly modelled. A next generation matrix approach was used to derive an expression for the basic reproduction number, R 0. A full sensitivity analysis was performed using a methodology based on Latin-hypercube sampling partial rank correlation coefficient index.

RESULTS

EPICYST outputs indicate that chemotherapeutic intervention targeted at humans or pigs would be highly effective at reducing taeniasis and cysticercosis prevalence when applied singly, with annual chemotherapy of humans and pigs resulting, respectively, in 94 and 74% of human cysticercosis cases averted. Improved sanitation, meat inspection and animal husbandry are less effective but are still able to reduce prevalence singly or in combination. The value of R 0 for taeniasis was estimated at 1.4 (95% Credible Interval: 0.5-3.6).

CONCLUSIONS

Human- and pig-targeted drug-focussed interventions appear to be the most efficacious approach from the options currently available. The model presented is a forward step towards developing an informed control and elimination strategy for cysticercosis. Together with its validation against field data, EPICYST will be a valuable tool to help reach the WHO goals and to conduct economic evaluations of interventions in varying epidemiological settings.

Balancing biomedical and ecological perspectives in research framing of liver fluke and cholangiocarcinoma in NE Thailand

This review examines the association of Asian liver flukes and cholangiocarcinoma (CCA) from the standpoint of two contrasting research perspectives: that aligned with the biomedical model predominantly employed to date; and, that aligned with ecological (and evolutionary) thinking increasingly being used to frame research questions that address this association in Northeast Thailand. An examination of the assumptions that underlie most of this research, requisite of evidence-based health research, shows how a broadened research frame that incorporates 'ecologic' perspectives provides alternatives to the prevailing scientific interpretations and public narrative. A more balanced and integrative research approach that combines elements of the biomedical model and ecologic models of health is suggested to overcome the limited progress toward the reduction of liver fluke infection prevalence and CCA incidence in this region. Similarly, this approach presents an opportunity to further enhance collaborative research programs involving Parasitology and the complementary fields in the health sciences.

Mathematical modelling of lymphatic filariasis elimination programmes in India: required duration of mass drug administration and post-treatment level of infection indicators

Background

India has made great progress towards the elimination of lymphatic filariasis. By 2015, most endemic districts had completed at least five annual rounds of mass drug administration (MDA). The next challenge is to determine when MDA can be stopped. We performed a simulation study with the individual-based model LYMFASIM to help clarify this.

Methods

We used a model-variant for Indian settings. We considered different hypotheses on detectability of antigenaemia (Ag) in relation to underlying adult worm burden, choosing the most likely hypothesis by comparing the model predicted association between community-level microfilaraemia (Mf) and antigenaemia (Ag) prevalence levels to observed data (collated from literature). Next, we estimated how long MDA must be continued in order to achieve elimination in different transmission settings and what Mf and Ag prevalence may still remain 1 year after the last required MDA round. The robustness of key-outcomes was assessed in a sensitivity analysis.

Results

Our model matched observed data qualitatively well when we assumed an Ag detection rate of 50 % for single worm infections, which increases with the number of adult worms (modelled by relating detection to the presence of female worms). The required duration of annual MDA increased with higher baseline endemicity and lower coverage (varying between 2 and 12 rounds), while the remaining residual infection 1 year after the last required treatment declined with transmission intensity. For low and high transmission settings, the median residual infection levels were 1.0 % and 0.4 % (Mf prevalence in the 5+ population), and 3.5 % and 2.0 % (Ag prevalence in 6–7 year-old children).

Conclusion

To achieve elimination in high transmission settings, MDA must be continued longer and infection levels must be reduced to lower levels than in low-endemic communities. Although our simulations were for Indian settings, qualitatively similar patterns are also expected in other areas. This should be taken into account in decision algorithms to define whether MDA can be interrupted. Transmission assessment surveys should ideally be targeted to communities with the highest pre-control transmission levels, to minimize the risk of programme failure.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1768-y) contains supplementary material, which is available to authorized users.

Multi-host model and threshold of intermediate host Oncomelania snail density for eliminating schistosomiasis transmission in China

Schistosomiasis remains a serious public health issue in many tropical countries, with more than 700 million people at risk of infection. In China, a national integrated control strategy, aiming at blocking its transmission, has been carried out throughout endemic areas since 2005. A longitudinal study was conducted to determine the effects of different intervention measures on the transmission dynamics of S. japonicum in three study areas and the data were analyzed using a multi-host model. The multi-host model was also used to estimate the threshold of Oncomelania snail density for interrupting schistosomiasis transmission based on the longitudinal data as well as data from the national surveillance system for schistosomiasis. The data showed a continuous decline in the risk of human infection and the multi-host model fit the data well. The 25th, 50th and 75th percentiles, and the mean of estimated thresholds of Oncomelania snail density below which the schistosomiasis transmission cannot be sustained were 0.006, 0.009, 0.028 and 0.020 snails/0.11 m², respectively. The study results could help develop specific strategies of schistosomiasis control and elimination tailored to the local situation for each endemic area.

Estimating the Future Impact of a Multi-Pronged Intervention Strategy on Ocular Disease Sequelae Caused by Trachoma: A Modeling Study

Purpose: Trachoma control programs are underway in endemic regions worldwide. They are based on the SAFE strategy (Surgery for trichiasis, Antibiotic distribution, Facial cleanliness, and Environmental improvement). Although much is known about the effect of community-wide treatment with antibiotics on the prevalence of Chlamydia trachomatis, the impact of the SAFE strategy on severe ocular disease sequelae (the main focus of the Global Elimination of blinding Trachoma by 2020 program) remains largely unknown.

Methods: We use a mathematical model to explore the impact of each of the components of the SAFE strategy, individually and together, on disease sequelae, arising from repeat infection and subsequent conjunctival scarring. We ask whether two elimination goals, to reduce the prevalence of trachomatous trichiasis to 1 per 1000 persons, and the incidence of corneal opacity to 1 per 10,000 persons per annum, are achievable, and which combinations of interventions have the greatest impact on these indicators.

Results: In high prevalence communities (here, >20% infection of children aged 1–9 years), a combination of efforts is needed to bring down sustainably the prevalence and incidence of ocular disease sequelae.

Conclusion: The mass delivery of antibiotics is highly beneficial for the clearance of infection, inflammation and prevention of subsequent scarring, but needs to be supplemented with sustained reductions in transmission and surgery to consider realistically the elimination of blindness by the year 2020.

MODELING THE SPREAD OF HOOKWORM DISEASE AND ASSESSING CHEMOTHERAPY PROGRAMS: MATHEMATICAL ANALYSIS AND COMPARISON WITH SURVEILLANCE DATA

Soil-transmitted helminthes including hookworm infections result in a major disease burden worldwide. Periodic chemotherapy treatments with anthelminthic drugs remain major intervention strategy in highly endemic areas and currently there is no approved human vaccine. We developed a model to study the spread of the hookworm infection in the population level having a variety of parasite development stages. We investigated long-term effectiveness of selective and mass chemotherapy, scenarios in which the therapy becomes less effective, and is interrupted. We analyzed the mathematical model, which includes determining the basic reproductive number and steady states, defining model initial conditions, proving positiveness and boundedness of solutions, local and global stability of the disease-free and endemic steady states, and uniform persistence of the system. We fitted our model using field data from hookworm control programs in endemic areas in China and Zimbabwe. The model provided accurate predictions of disease prevalence for all considered locations for multiple years of utilizing deworming programs. We demonstrated that administration of the deworming therapy is an effective method in terms of controlling the hookworm infection. Nevertheless, the chemotherapy programs do not eradicate the disease in the villages and the infection occurrence promptly increases once the program is interrupted, thus ultimately leading to the pre-treatment levels. Our modeling results suggest that the deworming programs should be maintained at least once per year to control the infection.

New approaches to measuring anthelminthic drug efficacy: parasitological responses of childhood schistosome infections to treatment with praziquantel

Background

By 2020, the global health community aims to control and eliminate human helminthiases, including schistosomiasis in selected African countries, principally by preventive chemotherapy (PCT) through mass drug administration (MDA) of anthelminthics. Quantitative monitoring of anthelminthic responses is crucial for promptly detecting changes in efficacy, potentially indicative of emerging drug resistance. Statistical models offer a powerful means to delineate and compare efficacy among individuals, among groups of individuals and among populations.

Methods

We illustrate a variety of statistical frameworks that offer different levels of inference by analysing data from nine previous studies on egg counts collected from African children before and after administration of praziquantel.

Results

We quantify responses to praziquantel as egg reduction rates (ERRs), using different frameworks to estimate ERRs among population strata, as average responses, and within strata, as individual responses. We compare our model-based average ERRs to corresponding model-free estimates, using as reference the World Health Organization (WHO) 90 % threshold of optimal efficacy. We estimate distributions of individual responses and summarize the variation among these responses as the fraction of ERRs falling below the WHO threshold.

Conclusions

Generic models for evaluating responses to anthelminthics deepen our understanding of variation among populations, sub-populations and individuals. We discuss the future application of statistical modelling approaches for monitoring and evaluation of PCT programmes targeting human helminthiases in the context of the WHO 2020 control and elimination goals.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1312-0) contains supplementary material, which is available to authorized users.

River Blindness: Mathematical Models for Control and Elimination

Human onchocerciasis (river blindness) is one of the few neglected tropical diseases (NTDs) whose control strategies have been informed by mathematical modelling. With the change in focus from elimination of the disease burden to elimination of Onchocerca volvulus, much remains to be done to refine, calibrate and validate existing models. Under the impetus of the NTD Modelling Consortium, the teams that developed EPIONCHO and ONCHOSIM have joined forces to compare and improve these frameworks to better assist ongoing elimination efforts. We review their current versions and describe how they are being used to address two key questions: (1) where can onchocerciasis be eliminated with current intervention strategies by 2020/2025? and (2) what alternative/complementary strategies could help to accelerate elimination where (1) cannot be achieved? The control and elimination of onchocerciasis from the African continent is at a crucial crossroad. The African Programme for Onchocerciasis Control closed at the end of 2015, and although a new platform for support and integration of NTD control has been launched, the disease will have to compete with a myriad of other national health priorities at a pivotal time in the road to elimination. However, never before had onchocerciasis control a better arsenal of intervention strategies as well as diagnostics. It is, therefore, timely to present two models of different geneses and modelling traditions as they come together to produce robust decision-support tools. We start by describing the structural and parametric assumptions of EPIONCHO and ONCHOSIM; we continue by summarizing the modelling of current treatment strategies with annual (or biannual) mass ivermectin distribution and introduce a number of alternative strategies, including other microfilaricidal therapies (such as moxidectin), macrofilaricidal (anti-wolbachial) treatments, focal vector control and the possibility of an onchocerciasis vaccine. We conclude by discussing challenges, opportunities and future directions.

The Role of More Sensitive Helminth Diagnostics in Mass Drug Administration Campaigns: Elimination and Health Impacts

Diagnostics play a crucial role in determining treatment protocols and evaluating success of mass drug administration (MDA) programmes used to control soil-transmitted helminths (STHs). The current diagnostic, Kato-Katz, relies on inexpensive, reusable materials and can be used in the field, but only trained microscopists can read slides. This diagnostic always underestimates the true prevalence of infection, and the accuracy worsens as the true prevalence falls. We investigate how more sensitive diagnostics would impact on the management and life cycle of MDA programmes, including number of mass treatment rounds, health impact, number of unnecessary treatments and probability of elimination. We use an individual-based model of STH transmission within the current World Health Organization (WHO) treatment guidelines which records individual disability-adjusted life years (DALY) lost. We focus on Ascaris lumbricoides due to the availability of high-quality data on existing diagnostics. We show that the effect of improving the sensitivity of diagnostics is principally determined by the precontrol prevalence in the community. Communities at low true prevalence (<30%) and high true prevalence (>70%) do not benefit greatly from improved diagnostics. Communities with intermediate prevalence benefit greatly from increased chemotherapy application, both in terms of reduced DALY loss and increased probability of elimination. Our results suggest that programmes should be extended beyond school-age children, especially in high prevalence communities. Finally, we argue against using apparent or measured prevalence as an uncorrected proxy for true prevalence.

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.

Potential effects of warmer worms and vectors on onchocerciasis transmission in West Africa

Development times of eggs, larvae and pupae of vectors of onchocerciasis (Simulium spp.) and of Onchocerca volvulus larvae within the adult females of the vectors decrease with increasing temperature. At and above 25°C, the parasite could reach its infective stage in less than 7 days when vectors could transmit after only two gonotrophic cycles. After incorporating exponential functions for vector development into a novel blackfly population model, it was predicted that fly numbers in Liberia and Ghana would peak at air temperatures of 29°C and 34°C, about 3°C and 7°C above current monthly averages, respectively; parous rates of forest flies (Liberia) would peak at 29°C and of savannah flies (Ghana) at 30°C. Small temperature increases (less than 2°C) might lead to changes in geographical distributions of different vector taxa. When the new model was linked to an existing framework for the population dynamics of onchocerciasis in humans and vectors, transmission rates and worm loads were projected to increase with temperature to at least 33°C. By contrast, analyses of field data on forest flies in Liberia and savannah flies in Ghana, in relation to regional climate change predictions, suggested, on the basis of simple regressions, that 13–41% decreases in fly numbers would be expected between the present and before 2040. Further research is needed to reconcile these conflicting conclusions.

A Theoretical Analysis of the Geography of Schistosomiasis in Burkina Faso Highlights the Roles of Human Mobility and Water Resources Development in Disease Transmission

We study the geography of schistosomiasis across Burkina Faso by means of a spatially explicit model of water-based disease dynamics. The model quantitatively addresses the geographic stratification of disease burden in a novel framework by explicitly accounting for drivers and controls of the disease, including spatial information on the distributions of population and infrastructure, jointly with a general description of human mobility and climatic/ecological drivers. Spatial patterns of disease are analysed by the extraction and the mapping of suitable eigenvectors of the Jacobian matrix subsuming the stability of the disease-free equilibrium. The relevance of the work lies in the novel mapping of disease burden, a byproduct of the parametrization induced by regional upscaling, by model-guided field validations and in the predictive scenarios allowed by exploiting the range of possible parameters and processes. Human mobility is found to be a primary control at regional scales both for pathogen invasion success and the overall distribution of disease burden. The effects of water resources development highlighted by systematic reviews are accounted for by the average distances of human settlements from water bodies that are habitats for the parasite's intermediate host. Our results confirm the empirical findings about the role of water resources development on disease spread into regions previously nearly disease-free also by inspection of empirical prevalence patterns. We conclude that while the model still needs refinements based on field and epidemiological evidence, the proposed framework provides a powerful tool for large-scale public health planning and schistosomiasis management.

Bayesian calibration of simulation models for supporting management of the elimination of the macroparasitic disease, Lymphatic Filariasis

Background

Mathematical models of parasite transmission can help integrate a large body of information into a consistent framework, which can then be used for gaining mechanistic insights and making predictions. However, uncertainty, spatial variability and complexity, can hamper the use of such models for decision making in parasite management programs.

Methods

We have adapted a Bayesian melding framework for calibrating simulation models to address the need for robust modelling tools that can effectively support management of lymphatic filariasis (LF) elimination in diverse endemic settings. We applied this methodology to LF infection and vector biting data from sites across the major LF endemic regions in order to quantify model parameters, and generate reliable predictions of infection dynamics along with credible intervals for modelled output variables. We used the locally calibrated models to estimate breakpoint values for various indicators of parasite transmission, and simulate timelines to parasite extinction as a function of local variations in infection dynamics and breakpoints, and effects of various currently applied and proposed LF intervention strategies.

Results

We demonstrate that as a result of parameter constraining by local data, breakpoint values for all the major indicators of LF transmission varied significantly between the sites investigated. Intervention simulations using the fitted models showed that as a result of heterogeneity in local transmission and extinction dynamics, timelines to parasite elimination in response to the current Mass Drug Administration (MDA) and various proposed MDA with vector control strategies also varied significantly between the study sites. Including vector control, however, markedly reduced the duration of interventions required to achieve elimination as well as decreased the risk of recrudescence following stopping of MDA.

Conclusions

We have demonstrated how a Bayesian data-model assimilation framework can enhance the use of transmission models for supporting reliable decision making in the management of LF elimination. Extending this framework for delivering predictions in settings either lacking or with only sparse data to inform the modelling process, however, will require development of procedures to estimate and use spatio-temporal variations in model parameters and inputs directly, and forms the next stage of the work reported here.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1132-7) contains supplementary material, which is available to authorized users.

Cost and cost-effectiveness of soil-transmitted helminth treatment programmes: systematic review and research needs

BACKGROUND

In this time of rapidly expanding mass drug administration (MDA) coverage and the new commitments for soil-transmitted helminth (STH) control, it is essential that resources are allocated in an efficient manner to have the greatest impact. However, many questions remain regarding how best to deliver STH treatment programmes; these include which age-groups should be targeted and how often. To perform further analyses to investigate what the most cost-effective control strategies are in different settings, accurate cost data for targeting different age groups at different treatment frequencies (in a range of settings) are essential.

METHODS

Using the electronic databases PubMed, MEDLINE, and ISI Web of Knowledge, we perform a systematic review of costing studies and cost-effectiveness evaluations for potential STH treatment strategies. We use this review to highlight research gaps and outline the key future research needs.

RESULTS

We identified 29 studies reporting costs of STH treatment and 17 studies that investigated its cost-effectiveness. The majority of these pertained to programmes only targeting school-aged children (SAC), with relatively few studies investigating alternative preventive chemotherapy (PCT) treatment strategies. The methods of cost data collection, analysis and reporting were highly variable among the different studies. Only four of the costing studies were found to have high applicability for use in forthcoming economic evaluations. There are also very few studies quantifying the costs of increasing the treatment frequency.

CONCLUSIONS

The absence of cost data and inconsistencies in the collection and analysis methods constitutes a major research gap for STH control. Detailed and accurate costs of targeting different age groups or increasing treatment frequency will be essential to formulate cost-effective public health policy. Defining the most cost-effective control strategies in different settings is of high significance during this period of expanding MDA coverage and new resource commitments for STH control.

Modelling the effects of mass drug administration on the molecular epidemiology of schistosomes

As national governments scale up mass drug administration (MDA) programs aimed to combat neglected tropical diseases (NTDs), novel selection pressures on these parasites increase. To understand how parasite populations are affected by MDA and how to maximize the success of control programmes, it is imperative for epidemiological, molecular and mathematical modelling approaches to be combined. Modelling of parasite population genetic and genomic structure, particularly of the NTDs, has been limited through the availability of only a few molecular markers to date. The landscape of infectious disease research is being dramatically reshaped by next-generation sequencing technologies and our understanding of how repeated selective pressures are shaping parasite populations is radically altering. Genomics can provide high-resolution data on parasite population structure, and identify how loci may contribute to key phenotypes such as virulence and/or drug resistance. We discuss the incorporation of genetic and genomic data, focussing on the recently sequenced Schistosoma spp., into novel mathematical transmission models to inform our understanding of the impact of MDA and other control methods. We summarize what is known to date, the models that exist and how population genetics has given us an understanding of the effects of MDA on the parasites. We consider how genetic and genomic data have the potential to shape future research, highlighting key areas where data are lacking, and how future molecular epidemiology knowledge can aid understanding of transmission dynamics and the effects of MDA, ultimately informing public health policy makers of the best interventions for NTDs.

Modelling lymphatic filariasis transmission and control: modelling frameworks, lessons learned and future directions

Mathematical modelling provides a useful tool for policy making and planning in lymphatic filariasis control programmes, by providing trend forecasts based on sound scientific knowledge and principles. This is now especially true, in view of the ambitious target to eliminate lymphatic filariasis as a public health problem globally by the year 2020 and the short remaining timeline to achieve this. To meet this target, elimination programmes need to be accelerated, requiring further optimization of strategies and tailoring to local circumstances. Insights from epidemiological transmission models provide a useful basis. Two general models of lymphatic filariasis transmission and control are nowadays in use to support decision-making, namely a population-based deterministic model (EPIFIL) and an individual-based stochastic model (LYMFASIM). Model predictions confirm that lymphatic filariasis transmission can be interrupted by annual mass drug administration (MDA), but this may need to be continued much longer than the initially suggested 4-6 years in areas with high transmission intensity or poor treatment coverage. However, the models have not been validated against longitudinal data describing the impact of MDA programmes. Some critical issues remain to be incorporated in one or both of the models to make predictions on elimination more realistic, including the possible occurrence of systematic noncompliance, the risk of emerging parasite resistance to anthelmintic drugs, and spatial heterogeneities. Rapid advances are needed to maximize the utility of models in decision-making for the ongoing ambitious lymphatic filariasis elimination programmes.

Toward integrated opisthorchiasis control in northeast Thailand: the Lawa project

Human liver fluke, Opisthorchis viverrini, a food-borne trematode is a significant public health problem in Southeast Asia, particularly in Thailand. Despite a long history of control programmes in Thailand and a nationwide reduction, O. viverrini infection prevalence remains high in the northeastern provinces. Therefore, a new strategy for controlling the liver fluke infection using the EcoHealth/One Health approach was introduced into the Lawa Lake area in Khon Kaen province where the liver fluke is endemic. A programme has been carried using anthelminthic treatment, novel intensive health education methods both in the communities and in schools, ecosystem monitoring and active community participation. As a result, the infection rate in the more than 10 villages surrounding the lake has declined to approximate one third of the average of 50% as estimated by a baseline survey. Strikingly, the Cyprinoid fish species in the lake, which are the intermediate host, now showed less than 1% prevalence compared to a maximum of 70% at baseline. This liver fluke control programme, named "Lawa model," is now recognised nationally and internationally, and being expanding to other parts of Thailand and neighbouring Mekong countries. Challenges to O. viverrini disease control, and lessons learned in developing an integrative control programme using a community-based, ecosystem approach, and scaling-up regionally based on Lawa as a model are described.

Integrated control programmes for schistosomiasis and other helminth infections in P.R. China

The prevalence of human schistosomiasis and soil-transmitted helminthiasis (STH) has decreased significantly in the People's Republic of China (P.R. China), particularly after 2005 when the national control programmes were reinforced by forming of integrated control strategies. Furthermore, social-economic development also contributed to the decrease of schistosome and soil-transmitted helminth infections. The prevalence of the zoonotic helminthiasis, including clonorchiasis and echinococcosis, on the other hand, is either underestimated or has in fact increased due to changes in social and environmental factors. In comparison with the control strategies in force and their effects on those four kinds of helminthiasis, the challenges and control priorities for the potential transfer from control to elimination of each disease is reviewed, to provide evidence for policy-makers to act upon.

Neglected tools for neglected diseases: mathematical models in economic evaluations

Despite many current interventions against neglected tropical diseases (NTDs) being highly cost-effective, new strategies are needed to reach the WHO's control and elimination goals. Here we argue for the importance of incorporating economic evaluations of new strategies in decisions regarding resource allocation. Such evaluation should ideally be conducted using dynamic transmission models that capture inherent nonlinearities in transmission and the indirect benefits ('herd effects') of interventions. A systematic review of mathematical models that have been used for economic analysis of interventions against the ten NTDs covered by the London Declaration reveals that only 16 out of 49 studies used dynamic transmission models, highlighting a fundamental--but addressable--gap in the evaluation of interventions against NTDs.

Seven challenges for modelling indirect transmission: vector-borne diseases, macroparasites and neglected tropical diseases

Many of the challenges which face modellers of directly transmitted pathogens also arise when modelling the epidemiology of pathogens with indirect transmission – whether through environmental stages, vectors, intermediate hosts or multiple hosts. In particular, understanding the roles of different hosts, how to measure contact and infection patterns, heterogeneities in contact rates, and the dynamics close to elimination are all relevant challenges, regardless of the mode of transmission. However, there remain a number of challenges that are specific and unique to modelling vector-borne diseases and macroparasites. Moreover, many of the neglected tropical diseases which are currently targeted for control and elimination are vector-borne, macroparasitic, or both, and so this article includes challenges which will assist in accelerating the control of these high-burden diseases. Here, we discuss the challenges of indirect measures of infection in humans, whether through vectors or transmission life stages and in estimating the contribution of different host groups to transmission. We also discuss the issues of “evolution-proof” interventions against vector-borne disease.

Neglected tropical diseases: now more than just 'other diseases'--the post-2015 agenda

Neglected tropical diseases (NTDs) have become recognised as important health problems facing at least a billion people in the low-income countries and the poorest communities in middle-income countries. WHO plays a leading role in developing strategies to address these diseases, pharmaceutical companies provide drug donations to treat or control the NTDs and many partners from different constituencies have become increasingly committed to their control or elimination. This review looks to the post-2015 agenda and emphasises that despite the progress made over recent years, if the targets established are to be achieved, then not only will additional financial resources be required to up-scale treatments and increase access, but increased applied and operational research will be necessary to address problems and human capacity in NTD skills will need to be strengthened. Continuing advocacy for the relevance of control or elimination of NTDs must be placed in the context of universal health coverage and access to donated essential medicines for the poor as a right. The evidence that investment in NTD interventions are cost-effective and impact not only on health, but also to enhance socio-economic development, must be refined and promulgated. The global burden of disease attributable to NTDs requires reassessment to appropriately define the true burden, while the potential for unexpected events, political, climatic, environmental as well as biological, have the potential to reduce future progress towards the agreed post-2015 targets. NTD progress towards the WHO Roadmap targets and the fulfilment of the World Health Assembly Resolution 66.12 of 2013 demand continued commitment from all partner constituencies when challenges emerge.

Modelling the impact of ivermectin on River Blindness and its burden of morbidity and mortality in African Savannah: EpiOncho projections

Background

The African Programme for Onchocerciasis Control (APOC) has refocused its goals on the elimination of infection where possible, seemingly achievable by 15–17 years of annual mass distribution of ivermectin in some African foci. Previously, APOC had focused on the elimination of onchocerciasis as a public health problem. Timeframes have been set by the World Health Organization, the London Declaration on Neglected Tropical Diseases and the World Bank to achieve these goals by 2020–2025.

Methods

A novel mathematical model of the dynamics of onchocercal disease is presented which links documented associations between Onchocerca volvulus infection and the prevalence and incidence of morbidity and mortality to model outputs from our host age- and sex-structured onchocerciasis transmission framework (EpiOncho). The model is calibrated for African savannah settings, and used to assess the impact of long-term annual mass administration of ivermectin on infection and ocular and skin disease and to explore how this depends on epidemiological and programmatic variables.

Results

Current onchocerciasis disease projections, which do not account for excess mortality of sighted individuals with heavy microfilarial loads, underestimate disease burden. Long-term annual ivermectin treatment is highly effective at reducing both the morbidity and mortality associated with onchocerciasis, and this result is not greatly influenced by treatment coverage and compliance. By contrast, impact on microfilarial prevalence and intensity is highly dependent on baseline endemicity, treatment coverage and systematic non-compliance.

Conclusions

The goals of eliminating morbidity and infection with ivermectin alone are distinctly influenced by epidemiological and programmatic factors. Whilst the former goal is most certainly achievable, reaching the latter will strongly depend on initial endemicity (the higher the endemicity, the greater the magnitude of inter-treatment transmission), advising caution when generalising the applicability of successful elimination outcomes to other areas. The proportion of systematic non-compliers will become far more influential in terms of overall success in achieving elimination goals.