How important is the spatial movement of people in attempts to eliminate the transmission of human helminth infections by mass drug administration?

Human mobility contributes to the spatial dynamics of many infectious diseases, and understanding these dynamics helps us to determine the most effective ways to intervene and plan surveillance. In this paper, we describe a novel transmission model for the spatial dynamics of hookworm, a parasitic worm which is a common infection across sub-Saharan Africa, East Asia and the Pacific islands. We fit our model, with and without mobility, to data obtained from a sub-county in Kenya, and validate the model's predictions against the decline in prevalence observed over the course of a clustered randomized control trial evaluating methods of administering mass chemotherapy. We find that our model which incorporates human mobility is able to reproduce the observed patterns in decline of prevalence during the TUMIKIA trial, and additionally, that the widespread bounce-back of infection may be possible over many years, depending on the rates of people movement between villages. The results have important implications for the design of mass chemotherapy programmes for the elimination of human helminth transmission. This article is part of the theme issue 'Challenges and opportunities in the fight against neglected tropical diseases: a decade from the London Declaration on NTDs'.

A prevalence-based transmission model for the study of the epidemiology and control of soil-transmitted helminthiasis.

Much effort has been devoted by the World Health Organization (WHO) to eliminate soil-transmitted helminth (STH) infections by 2030 using mass drug administration targeted at particular risk groups alongside the availability to access water, sanitation and hygiene services. The targets set by the WHO for the control of helminth infections are typically defined in terms of the prevalence of infection, whereas the standard formulation of STH transmission models typically describe dynamic changes in the mean-worm burden. We develop a prevalence-based deterministic model to investigate the transmission dynamics of soil-transmitted helminthiasis in humans, subject to continuous exposure to infection over time. We analytically determine local stability criteria for all equilibria and find bifurcation points. Our model predicts that STH infection will either be eliminated (if the initial prevalence value, y(0), is sufficiently small) or remain endemic (if y(0) is sufficiently large), with the two stable points of endemic infection and parasite eradication separated by a transmission breakpoint. Two special cases of the model are analysed: (1) the distribution of the STH parasites in the host population is highly aggregated following a negative binomial distribution, and (2) no density-dependent effects act on the parasite population. We find that disease extinction is always possible for Case (1), but it is not so for Case (2) if y(0) is sufficiently large. However, by introducing stochastic perturbation into the deterministic model, we discover that chance effects can lead to outcomes not predicted by the deterministic model alone, with outcomes highly dependent on the degree of worm clumping, k. Specifically, we show that if the reproduction number and clumping are sufficiently bounded, then stochasticity will cause the parasite to die out. It follows that control of soil-transmitted helminths will be more difficult if the worm distribution tends towards clumping.

The observed relationship between the degree of parasite aggregation and the prevalence of infection within human host populations for soil-transmitted helminth and schistosome infections.

BACKGROUND: Soil-transmitted helminths (STH) and schistosome parasites are highly aggregated within the human population. The probability distribution of worms per person is described well by the negative binomial probability distribution with aggregation parameter, k, which varies inversely with parasite clustering. The relationship between k and prevalence in defined populations subject to mass drug administration is not well understood. METHODS AND RESULTS: We use statistical methods to estimate k using two large independent datasets for STH and schistosome infections from India and Niger, respectively, both of which demonstrate increased aggregation of parasites in a few hosts, as the prevalence of infections declines across the dataset. CONCLUSIONS: A greater attention needs to be given in monitoring and evaluation programmes to find and treat the remaining aggregates of parasites.

Analysing pneumococcal invasiveness using Bayesian models of pathogen progression rates.

The disease burden attributable to opportunistic pathogens depends on their prevalence in asymptomatic colonisation and the rate at which they progress to cause symptomatic disease. Increases in infections caused by commensals can result from the emergence of "hyperinvasive" strains. Such pathogens can be identified through quantifying progression rates using matched samples of typed microbes from disease cases and healthy carriers. This study describes Bayesian models for analysing such datasets, implemented in an RStan package (https://github.com/nickjcroucher/progressionEstimation). The models converged on stable fits that accurately reproduced observations from meta-analyses of Streptococcus pneumoniae datasets. The estimates of invasiveness, the progression rate from carriage to invasive disease, in cases per carrier per year correlated strongly with the dimensionless values from meta-analysis of odds ratios when sample sizes were large. At smaller sample sizes, the Bayesian models produced more informative estimates. This identified historically rare but high-risk S. pneumoniae serotypes that could be problematic following vaccine-associated disruption of the bacterial population. The package allows for hypothesis testing through model comparisons with Bayes factors. Application to datasets in which strain and serotype information were available for S. pneumoniae found significant evidence for within-strain and within-serotype variation in invasiveness. The heterogeneous geographical distribution of these genotypes is therefore likely to contribute to differences in the impact of vaccination in between locations. Hence genomic surveillance of opportunistic pathogens is crucial for quantifying the effectiveness of public health interventions, and enabling ongoing meta-analyses that can identify new, highly invasive variants.

What is the impact of acquired immunity on the transmission of schistosomiasis and the efficacy of current and planned mass drug administration programmes?

Schistosomiasis causes severe morbidity in many countries with endemic infection with the schistosome digenean parasites in Africa and Asia. To control and eliminate the disease resulting from infection, regular mass drug administration (MDA) is used, with a focus on school-aged children (SAC; 5-14 years of age). In some high transmission settings, the World Health Organization (WHO) also recommends the inclusion of at-risk adults in MDA treatment programmes. The question of whether ecology (age-dependant exposure) or immunity (resistance to reinfection), or some combination of both, determines the form of observed convex age-intensity profile is still unresolved, but there is a growing body of evidence that the human hosts acquire some partial level of immunity after a long period of repeated exposure to infection. In the majority of past research modelling schistosome transmission and the impact of MDA programmes, the effect of acquired immunity has not been taken into account. Past work has been based on the assumption that age-related contact rates generate convex horizontal age-intensity profiles. In this paper, we use an individual based stochastic model of transmission and MDA impact to explore the effect of acquired immunity in defined MDA programmes. Compared with scenarios with no immunity, we find that acquired immunity makes the MDA programme less effective with a slower decrease in the prevalence of infection. Therefore, the time to achieve morbidity control and elimination as a public health problem is longer than predicted by models with just age-related exposure and no build-up of immunity. The level of impact depends on the baseline prevalence prior to treatment (the magnitude of the basic reproductive number R0) and the treatment frequency, among other factors. We find that immunity has a larger impact within moderate to high transmission settings such that it is very unlikely to achieve morbidity and transmission control employing current MDA programmes.

Spatial scales in human movement between reservoirs of infection.

The life cycle of parasitic organisms that are the cause of much morbidity in humans often depend on reservoirs of infection for transmission into their hosts. Understanding the daily, monthly and yearly movement patterns of individuals between reservoirs is therefore of great importance to implementers of control policies seeking to eliminate various parasitic diseases as a public health problem. This is due to the fact that the underlying spatial extent of the reservoir of infection, which drives transmission, can be strongly affected by inputs from external sources, i.e., individuals who are not spatially attributed to the region defined by the reservoir itself can still migrate and contribute to it. In order to study the importance of these effects, we build and examine a novel theoretical model of human movement between spatially-distributed focal points for infection clustered into regions defined as 'reservoirs of infection'. Using our model, we vary the spatial scale of human moment defined around focal points and explicitly calculate how varying this definition can influence the temporal stability of the effective transmission dynamics - an effect which should strongly influence how control measures, e.g., mass drug administration (MDA), define evaluation units (EUs). Considering the helminth parasites as our main example, by varying the spatial scale of human movement, we demonstrate that a critical scale exists around infectious focal points at which the migration rate into their associated reservoir can be neglected for practical purposes. This scale varies by species and geographic region, but is generalisable as a concept to infectious reservoirs of varying spatial extents and shapes. Our model is designed to be applicable to a very general pattern of infectious disease transmission modified by the migration of infected individuals between clustered communities. In particular, it may be readily used to study the spatial structure of hosts for macroparasites with temporally stationary distributions of infectious focal point locations over the timescales of interest, which is viable for the soil-transmitted helminths and schistosomes. Additional developments will be necessary to consider diseases with moving reservoirs, such as vector-born filarial worm diseases.

Deworming women of reproductive age during adolescence and pregnancy: what is the impact on morbidity from soil-transmitted helminths infection?

BACKGROUND: Soil-transmitted helminths (STHs) are a major cause of poor health in low- and middle-income countries. In particular, hookworm is known to cause anaemia in children and women of reproductive age (WRA). One goal of the World Health Organization's (WHO) 2030 roadmap for neglected tropical diseases is to reduce STH-related morbidity in WRA. As a minimal intervention, the WHO recommends deworming adolescent girls annually during human papilloma virus vaccination programmes and WRA during pregnancy and lactation. These routine interventions are low cost and can be implemented even by the most basic health services in endemic countries. In this study we use a cohort model to investigate the potential impact on STH-related morbidity in WRA. RESULTS: Annual deworming treatment of adolescent girls reduces the prevalence of moderate- and heavy-intensity infections in this age group by up to 60% in moderate transmission settings and by 12-27% in high transmission settings. Treatment of WRA during pregnancy and lactation on its own has a small (< 20%) but significant effect on morbidity although it does not lead to the achievement of the morbidity target (< 2% moderate- to high-intensity infections) in this age group. However, depending on the age-intensity profile of infection, which may vary geographically, and assumptions on the density-dependence of egg production by fertilised female worms, continued school-based treatment may be able to reduce the force of infection acting on WRA, both through an indirect effect on the overall population-based force of infection and via reducing the burden of infection as children age and move into the WRA age classes. As a result, morbidity in WRA may be eliminated. CONCLUSION: While deworming during pregnancy and lactation does not lead to the achievement of the morbidity target in WRA and its efficacy may vary by setting, it is still expected to be beneficial for maternity and child health. Monitoring of any WRA-based intervention is recommended to evaluate its effectiveness.

Stochastic challenges to interrupting helminth transmission.

Predicting the effect of different programmes designed to control both the morbidity induced by helminth infections and parasite transmission is greatly facilitated by the use of mathematical models of transmission and control impact. In such models, it is essential to account for the many sources of uncertainty - natural, or otherwise - to ensure robustness in prediction and to accurately depict variation around an expected outcome. In this paper, we investigate how well the standard deterministic models match the predictions made using individual-based stochastic simulations. We also explore how well concepts which derive from deterministic models, such as 'breakpoints' in transmission, apply in the stochastic world. Employing an individual-based stochastic model framework we also investigate how transmission and control are affected by the migration of infected people into a defined community. To give our study focus we consider the control of soil-transmitted helminths (STH) by mass drug administration (MDA), though our methodology is readily applicable to the other helminth species such as the schistosome parasites and the filarial worms. We show it is possible to theoretically define a 'stochastic breakpoint' where much noise surrounds the expected deterministic breakpoint. We also discuss the concept of the 'interruption of transmission' independent of the 'breakpoint' concept where analyses of model behaviour illustrate the current limitations of deterministic models to account for the 'fade-out' or transmission extinction behaviour in simulations. Our analysis of migration confirms a relationship between the critical infected human migration rate scale (i.e., order of magnitude) per unit of time and the death rate of infective stages that are released into the free-living environment. This relationship is shown to determine the likelihood that control activities aim at chemotherapeutic treatment of the human host will eliminate transmission. The development of a new stochastic simulation code for STH in the form of a publicly-available open-source python package which includes features to incorporate many population stratifications, different control interventions including mass drug administration (with defined frequency, coverage levels and compliance patterns) and inter-village human migration is also described.

Individual adherence to mass drug administration in neglected tropical disease control: A probability model conditional on past behaviour.

We present a general framework which describes the systematic (binary) scenario of individuals either taking treatment or not for any reason, over the course of mass drug administration (MDA)-which we refer to as 'adherence' and 'non-adherence'. The probability models developed can be informed by observed adherence behaviour as well as employed to explore how different patterns influence the impact of MDA programmes, by the use of mathematical models of transmission and control. We demonstrate the interpretative value of the developed probability model employing a dataset collected in the TUMIKIA project, a randomised trial of deworming strategies to control soil-transmitted helminths (STH) by MDA conducted in coastal Kenya. We stratify our analysis by age and sex, although the framework which we introduce here may be readily adapted to accommodate other stratifications. Our findings include the detection of specific patterns of non-adherence in all age groups to varying extents. This is particularly apparent in men of ages 30+. We then demonstrate the use of the probability model in stochastic individual-based simulations by running two example forecasts for the elimination of STH transmission employing MDA within the TUMIKIA trial setting with different adherence patterns. This suggested a substantial reduction in the probability of elimination (between 23-43%) when comparing observed adherence patterns with an assumption of independence, with important implications for programmes. The results here demonstrate the considerable impact and utility of considering non-adherence on the success of MDA programmes to control neglected tropical diseases (NTDs).

Forecasting the effectiveness of the DeWorm3 trial in interrupting the transmission of soil-transmitted helminths in three study sites in Benin, India and Malawi.

BACKGROUND: The DeWorm3 project is an ongoing cluster-randomised trial assessing the feasibility of interrupting the transmission of soil-transmitted helminths (STH) through mass drug administration (MDA) using study sites in India, Malawi and Benin. In this article, we describe an approach which uses a combination of statistical and mathematical methods to forecast the outcome of the trial with respect to its stated goal of reducing the prevalence of infection to below 2%. METHODS: Our approach is first to define the local patterns of transmission within each study site, which is achieved by statistical inference of key epidemiological parameters using the baseline epidemiological measures of age-related prevalence and intensity of STH infection which have been collected by the DeWorm3 trials team. We use these inferred parameters to calibrate an individual-based stochastic simulation of the trial at the cluster and study site level, which is subsequently run to forecast the future prevalence of STH infections. The simulator takes into account both the uncertainties in parameter estimation and the variability inherent in epidemiological and demographic processes in the simulator. We interpret the forecast results from our simulation with reference to the stated goal of the DeWorm3 trial, to achieve a target of [Formula: see text] prevalence at a point 24 months post-cessation of MDA. RESULTS: Simulated output predicts that the two arms will be distinguishable from each other in all three country sites at the study end point. In India and Malawi, measured prevalence in the intervention arm is below the threshold with a high probability (90% and 95%, respectively), but in Benin the heterogeneity between clusters prevents the arm prevalence from being reduced below the threshold value. At the level of individual study arms within each site, heterogeneity among clusters leads to a very low probability of achieving complete elimination in an intervention arm, yielding a post-study scenario with widespread elimination but a few 'hot spot' areas of persisting STH transmission. CONCLUSIONS: Our results suggest that geographical heterogeneities in transmission intensity and worm aggregation have a large impact on the effect of MDA. It is important to accurately assess cluster-level, or even smaller scale, heterogeneities in factors which influence transmission and aggregation for a clearer perspective on projecting the outcomes of MDA control of STH and other neglected tropical diseases.