Complex ecological dynamics and eradicability of the vector borne macroparasitic disease, lymphatic filariasis

Lymphatic filariasis endgame strategies: Using GEOFIL to model mass drug administration and targeted surveillance and treatment strategies in American Samoa

American Samoa underwent seven rounds of mass drug administration (MDA) for lymphatic filariasis (LF) from 2000-2006, but subsequent surveys found evidence of ongoing transmission. American Samoa has since undergone further rounds of MDA in 2018, 2019, and 2021; however, recent surveys indicate that transmission is still ongoing. GEOFIL, a spatially-explicit agent-based LF model, was used to compare the effectiveness of territory-wide triple-drug MDA (3D-MDA) with targeted surveillance and treatment strategies. Both approaches relied on treatment with ivermectin, diethylcarbamazine, and albendazole. We simulated three levels of whole population coverage for 3D-MDA: 65%, 73%, and 85%, while the targeted strategies relied on surveillance in schools, workplaces, and households, followed by targeted treatment. In the household-based strategies, we simulated 1-5 teams travelling village-to-village and offering antigen (Ag) testing to randomly selected households in each village. If an Ag-positive person was identified, treatment was offered to members of all households within 100m-1km of the positive case. All simulated interventions were finished by 2027 and their effectiveness was judged by their 'control probability'-the proportion of simulations in which microfilariae prevalence decreased between 2030 and 2035. Without future intervention, we predict Ag prevalence will rebound. With 3D-MDA, a 90% control probability required an estimated ≥ 4 further rounds with 65% coverage, ≥ 3 rounds with 73% coverage, or ≥ 2 rounds with 85% coverage. While household-based strategies were substantially more testing-intensive than 3D-MDA, they could offer comparable control probabilities with substantially fewer treatments; e.g. three teams aiming to test 50% of households and offering treatment to a 500m radius had approximately the same control probability as three rounds of 73% 3D-MDA, but used < 40% the number of treatments. School- and workplace-based interventions proved ineffective. Regardless of strategy, reducing Ag prevalence below the 1% target threshold recommended by the World Health Organization was a poor indicator of the interruption of LF transmission, highlighting the need to review blanket elimination targets.

Evaluating elimination thresholds and stopping criteria for interventions against the vector-borne macroparasitic disease, lymphatic filariasis, using mathematical modelling

We leveraged the ability of EPIFIL transmission models fit to field data to evaluate the use of the WHO Transmission Assessment Survey (TAS) for supporting Lymphatic Filariasis (LF) intervention stopping decisions. Our results indicate that understanding the underlying parasite extinction dynamics, particularly the protracted transient dynamics involved in shifts to the extinct state, is crucial for understanding the impacts of using TAS for determining the achievement of LF elimination. These findings warn that employing stopping criteria set for operational purposes, as employed in the TAS strategy, without a full consideration of the dynamics of extinction could seriously undermine the goal of achieving global LF elimination.

Neglected tropical diseases and systemic racism especially in Brazil: from socio-economic aspects to the development of new drugs

Neglected tropical diseases (NTDs) are highly prevalent communicable diseases in tropical and subtropical countries, generally not economically attractive for drug development and related to poverty. In Brazil, more specifically, socioeconomic inequalities and health indicators are strongly influenced by skin color, race, and ethnicity, due to the historical process of slavery. In this context, it is important to understand the concept of systemic racism: a form of indirect racial discrimination present in many institutions, which determines the process of illness and death of the black population, the ethnic group most affected by these diseases. The main objective of this paper was to carry out a literature review on the socioeconomic aspects of these diseases, relating them to institutional racism, and to encourage reflection on the influence of this type of racism in the NTDs context. Therefore, we present a paper that brings a evident correlation between racism versus neglected populations, which are affected by equally neglected diseases. A more humane and comprehensive view is needed to realize that these illnesses affect neglected and vulnerable populations, who require decent living conditions, health, and social justice. We hope to provide, with this paper, enough, but not exhaust, knowledge to initiate the discussion about neglected diseases, their socioeconomic aspects and institutional racism.

Lymphatic filariasis in 2016 in American Samoa: Identifying clustering and hotspots using non-spatial and three spatial analytical methods

Background

American Samoa completed seven rounds of mass drug administration from 2000–2006 as part of the Global Programme to Eliminate Lymphatic Filariasis (LF). However, resurgence was confirmed in 2016 through WHO-recommended school-based transmission assessment survey and a community-based survey. This paper uses data from the 2016 community survey to compare different spatial and non-spatial methods to characterise clustering and hotspots of LF.

Method

Non-spatial clustering of infection markers (antigen [Ag], microfilaraemia [Mf], and antibodies (Ab [Wb123, Bm14, Bm33]) was assessed using intra-cluster correlation coefficients (ICC) at household and village levels. Spatial dependence, clustering and hotspots were examined using semivariograms, Kulldorf’s scan statistic and Getis-Ord Gi* statistics based on locations of surveyed households.

Results

The survey included 2671 persons (750 households, 730 unique locations in 30 villages). ICCs were higher at household (0.20–0.69) than village levels (0.10–0.30) for all infection markers. Semivariograms identified significant spatial dependency for all markers (range 207–562 metres). Using Kulldorff’s scan statistic, significant spatial clustering was observed in two previously known locations of ongoing transmission: for all markers in Fagali’i and all Abs in Vaitogi. Getis-Ord Gi* statistic identified hotspots of all markers in Fagali’i, Vaitogi, and Pago Pago-Anua areas. A hotspot of Ag and Wb123 Ab was identified around the villages of Nua-Seetaga-Asili. Bm14 and Bm33 Ab hotspots were seen in Maleimi and Vaitogi-Ili’ili-Tafuna.

Conclusion

Our study demonstrated the utility of different non-spatial and spatial methods for investigating clustering and hotspots, the benefits of using multiple infection markers, and the value of triangulating results between methods.

Complexity and critical thresholds in the dynamics of visceral leishmaniasis

We study a general multi-host model of visceral leishmaniasis including both humans and animals, and where host and vector characteristics are captured via host competence along with vector biting preference. Additionally, the model accounts for spatial heterogeneity in human population and heterogeneity in biting behaviour of sandflies. We then use parameters for visceral leishmaniasis in the Indian subcontinent as an example and demonstrate that the model exhibits backward bifurcation, i.e. it has a human infection and a sandfly population threshold, characterized by a bi-stable region. These thresholds shift as a function of host competence, host population size, vector feeding preference, spatial heterogeneity, biting heterogeneity and control efforts. In particular, if control is applied through human treatment a new and lower human infection threshold is created, making elimination difficult to achieve, before eventually the human infection threshold no longer exists, making it impossible to control the disease by only reducing the infection levels below a certain threshold. A better strategy would be to reduce the human infection below a certain threshold potentially by early diagnosis, control animal population levels and keep the vector population under check. Spatial heterogeneity in human populations lowers the overall thresholds as a result of weak migration between patches.

Statistical methods for linking geostatistical maps and transmission models: Application to lymphatic filariasis in East Africa

•

Novel methodology for combining geostatistical mapping and transmission modelling.

•

Guide the planning of spatial control programmes by identifying affected areas.

•

Current intervention strategy will not be sufficient to eliminate LF in most areas.

•

Alternative strategies will be required to accelerate LF elimination in East Africa.

Infectious diseases remain one of the major causes of human mortality and suffering. Mathematical models have been established as an important tool for capturing the features that drive the spread of the disease, predicting the progression of an epidemic and hence guiding the development of strategies to control it. Another important area of epidemiological interest is the development of geostatistical methods for the analysis of data from spatially referenced prevalence surveys. Maps of prevalence are useful, not only for enabling a more precise disease risk stratification, but also for guiding the planning of more reliable spatial control programmes by identifying affected areas. Despite the methodological advances that have been made in each area independently, efforts to link transmission models and geostatistical maps have been limited. Motivated by this fact, we developed a Bayesian approach that combines fine-scale geostatistical maps of disease prevalence with transmission models to provide quantitative, spatially-explicit projections of the current and future impact of control programs against a disease. These estimates can then be used at a local level to identify the effectiveness of suggested intervention schemes and allow investigation of alternative strategies. The methodology has been applied to lymphatic filariasis in East Africa to provide estimates of the impact of different intervention strategies against the disease.

Simulating a Transmission Assessment Survey: An evaluation of current methods used in determining the elimination of the neglected tropical disease, Lymphatic Filariasis

INTRODUCTION

The World Health Organization (WHO) recommends Transmission Assessment Surveys (TAS) to determine when an evaluation unit (EU) (a designated population survey area) has achieved elimination of transmission of the vector-borne macroparasitic disease Lymphatic Filariasis (LF). These determinations are based on combining data from multiple survey units within an EU; it is unclear how underlying cluster-level variation influences the outcome of the TAS at EU level. We simulate LF infection distribution in an EU and compare three methods for assessing whether LF elimination has occurred based on currently recommended decision thresholds and sampling methods.

METHODS

We simulate an EU divided into clusters of varying size and disease prevalence. We produce 1000 samples according to LF TAS examples and WHO guidelines and compare three decision-making approaches: lot quality assurance sampling (LQAS) (recommended by WHO), one-sided interval estimate (CI), and n^th order statistic (MAX). Summary statistics demonstrating the "pass" rate for the EU under different disease transmission conditions are generated using a versatile SAS® macro.

RESULTS

As the prevalence of LF decreases, the LQAS and CI approaches produce increased likelihood of a pass outcome for an EU while some cluster units may still have a high likelihood of transmission. The MAX provides an alternative that increases the likelihood of determining a pass only once the whole area has a low likelihood of transmission. LQAS and CI approaches designed to estimate the LF prevalence in the EU miss hotspots that will continue to transmit infection while the MAX approach focuses on identifying clusters with high risk of transmission.

CONCLUSIONS

The current TAS methodology has a flaw that may result in false predictions of LF transmission interruption throughout an EU. Modifying the TAS methodology to address results from extreme clusters rather than being based on mean prevalence over an EU will result in greater success for global elimination of LF.

Defining a prevalence level to describe the elimination of Lymphatic Filariasis (LF) transmission and designing monitoring & evaluating (M&E) programmes post the cessation of mass drug administration (MDA)

The global decline in prevalence of lymphatic filariasis has been one of the major successes of the WHO's NTD programme. The recommended strategy of intensive, community-wide mass drug administration, aims to break localised transmission by either reducing the prevalence of microfilaria positive infections to below 1%, or antigen positive infections to below 2%. After the threshold is reached, and mass drug administration is stopped, geographically defined evaluation units must pass Transmission Assessment Surveys to demonstrate that transmission has been interrupted. In this study, we use an empirically parameterised stochastic transmission model to investigate the appropriateness of 1% microfilaria-positive prevalence as a stopping threshold, and statistically evaluate how well various monitoring prevalence-thresholds predict elimination or disease resurgence in the future by calculating their predictive value. Our results support the 1% filaremia prevalence target as appropriate stopping criteria. However, because at low prevalence-levels random events dominate the transmission dynamics, we find single prevalence measurements have poor predictive power for predicting resurgence, which suggests alternative criteria for restarting MDA may be beneficial.

Elimination or Resurgence: Modelling Lymphatic Filariasis After Reaching the 1% Microfilaremia Prevalence Threshold

The low prevalence levels associated with lymphatic filariasis elimination pose a challenge for effective disease surveillance. As more countries achieve the World Health Organization criteria for halting mass treatment and move on to surveillance, there is increasing reliance on the utility of transmission assessment surveys (TAS) to measure success. However, the long-term disease outcomes after passing TAS are largely untested. Using 3 well-established mathematical models, we show that low-level prevalence can be maintained for a long period after halting mass treatment and that true elimination (0% prevalence) is usually slow to achieve. The risk of resurgence after achieving current targets is low and is hard to predict using just current prevalence. Although resurgence is often quick (<5 years), it can still occur outside of the currently recommended postintervention surveillance period of 4-6 years. Our results highlight the need for ongoing and enhanced postintervention monitoring, beyond the scope of TAS, to ensure sustained success.

Data-driven modelling and spatial complexity supports heterogeneity-based integrative management for eliminating Simulium neavei-transmitted river blindness

Concern is emerging regarding the challenges posed by spatial complexity for modelling and managing the area-wide elimination of parasitic infections. While this has led to calls for applying heterogeneity-based approaches for addressing this complexity, questions related to spatial scale, the discovery of locally-relevant models, and its interaction with options for interrupting parasite transmission remain to be resolved. We used a data-driven modelling framework applied to infection data gathered from different monitoring sites to investigate these questions in the context of understanding the transmission dynamics and efforts to eliminate Simulium neavei- transmitted onchocerciasis, a macroparasitic disease that causes river blindness in Western Uganda and other regions of Africa. We demonstrate that our Bayesian-based data-model assimilation technique is able to discover onchocerciasis models that reflect local transmission conditions reliably. Key management variables such as infection breakpoints and required durations of drug interventions for achieving elimination varied spatially due to site-specific parameter constraining; however, this spatial effect was found to operate at the larger focus level, although intriguingly including vector control overcame this variability. These results show that data-driven modelling based on spatial datasets and model-data fusing methodologies will be critical to identifying both the scale-dependent models and heterogeneity-based options required for supporting the successful elimination of S. neavei-borne onchocerciasis.

Accelerating river blindness elimination by supplementing MDA with a vegetation "slash and clear" vector control strategy: a data-driven modeling analysis

Attention is increasingly focusing on how best to accelerate progress toward meeting the WHO's 2030 goals for neglected tropical diseases (NTDs). For river blindness, a major NTD targeted for elimination, there is a long history of using vector control to suppress transmission, but traditional larvicide-based approaches are limited in their utility. One innovative and sustainable approach, "slash and clear", involves clearing vegetation from breeding areas, and recent field trials indicate that this technique very effectively reduces the biting density of Simulium damnosum s.s. In this study, we use a Bayesian data-driven mathematical modeling approach to investigate the potential impact of this intervention on human onchocerciasis infection. We developed a novel "slash and clear" model describing the effect of the intervention on seasonal black fly biting rates and coupled this with our population dynamics model of Onchocerca volvulus transmission. Our results indicate that supplementing annual drug treatments with "slash and clear" can significantly accelerate the achievement of onchocerciasis elimination. The efficacy of the intervention is not very sensitive to the timing of implementation, and the impact is meaningful even if vegetation is cleared only once per year. As such, this community-driven technique will represent an important option for achieving and sustaining O. volvulus elimination.

Evaluating the Evidence for Lymphatic Filariasis Elimination

In the global drive for elimination of lymphatic filariasis (LF), 15 countries have achieved validation of elimination as a public health problem (EPHP). Recent empirical evidence has demonstrated that EPHP does not always lead to elimination of transmission (EOT). Here we show how the probability of elimination explicitly depends on key biological parameters, many of which have been poorly characterized, leading to a poor evidence base for the elimination threshold. As more countries progress towards EPHP it is essential that this process is well-informed, as prematurely halting treatment and surveillance programs could pose a serious threat to global progress. We highlight that refinement of the weak empirical evidence base is vital to understand drivers of elimination and inform long-term policy.

Understanding heterogeneities in mosquito-bite exposure and infection distributions for the elimination of lymphatic filariasis

It is well known that individuals in the same community can be exposed to a highly variable number of mosquito bites. This heterogeneity in bite exposure has consequences for the control of vector-borne diseases because a few people may be contributing significantly to transmission. However, very few studies measure sources of heterogeneity in a way which is relevant to decision-making. We investigate the relationship between two classic measures of heterogeneity, spatial and individual, within the context of lymphatic filariasis, a parasitic mosquito-borne disease. Using infection and mosquito-bite data for five villages in Papua New Guinea, we measure biting characteristics to model what impact bed-nets have had on control of the disease. We combine this analysis with geospatial modelling to understand the spatial relationship between disease indicators and nightly mosquito bites. We found a weak association between biting and infection heterogeneity within villages. The introduction of bed-nets increased biting heterogeneity, but the reduction in mean biting more than compensated for this, by reducing prevalence closer to elimination thresholds. Nightly biting was explained by a spatial heterogeneity model, while parasite load was better explained by an individual heterogeneity model. Spatial and individual heterogeneity are qualitatively different with profoundly different policy implications.

Mathematical model of the life cycle of taenia-cysticercosis: transmission dynamics and chemotherapy (Part 1)

BACKGROUND

Taenia solium is the aetiological agent of human taeniasis, pig cysticercosis and human neurocysticercosis, which are serious public health problems, especially in developing countries.

METHODS

A mathematical model of the transmission dynamics of taeniasis-cysticercosis is formulated. The model consists of a coupled system of differential equations, which are density-dependent equations for describing the flow of the parasite through the life cycle. The model is hybrid since it comprises deterministic equations with stochastic elements which describe changes in the mean parasite burden and incorporates the overall pattern of the parasites' distribution.

RESULTS

Sensitivity and bifurcation analyses were carried out to determine the range of values of the model. The model can reproduce the observed epidemiological patterns of human taeniasis, pig and human cysticercosis. For example, for a wide range of parameter values, the mean intensity of adult worms tends to rapidly stabilize in one parasite per individual host. From this model, we also derived a Susceptible-Infected model to describe the prevalence of infection in humans and pigs. Chemotherapeutic interventions against pig cysticercosis or human taeniasis may reduce rapidly and effectively the mean intensity of human taeniasis, pig cysticercosis and human cysticercosis. This effect can be achieved even if the protective efficacy of the drug is of the order of 90% and the coverage rate is 90%. This means that health in humans infected either with adult worms or cysticerci may be achieved by the application of anthelmintic drugs against pig cysticercosis. However, treatment against human cysticercosis alone, does not influence neither human teniasis nor pig cysticercosis. This is because human cysticercosis infection does not influence the value of the basic reproductive number (Ro).

CONCLUSIONS

Even coverage of 100% in the administration of anthelmintics did not eliminate the infection. Then elimination of the infection in all hosts does not seem a feasible goal to achieve by administering only chemotherapeutic interventions. Throughout the manuscript a discussion of our model in the context of other models of taeniasis-cysticercosis is presented.

Estimating the elimination feasibility in the 'end game' of control efforts for parasites subjected to regular mass drug administration: Methods and their application to schistosomiasis

Progress towards controlling and eliminating parasitic worms, including schistosomiasis, onchocerciasis, and lymphatic filariasis, is advancing rapidly as national governments, multinational NGOs, and pharmaceutical companies launch collaborative chemotherapeutic control campaigns. Critical questions remain regarding the potential for achieving elimination of these infections, and analytical methods can help to quickly estimate progress towards-and the probability of achieving-elimination over specific timeframes. Here, we propose the effective reproduction number, Reff, as a proxy of elimination potential for sexually reproducing worms that are subject to poor mating success at very low abundance (positive density dependence, or Allee effects). Reff is the number of parasites produced by a single reproductive parasite at a given stage in the transmission cycle, over the parasite's lifetime-it is the generalized form of the more familiar basic reproduction number, R0, which only applies at the beginning of an epidemic-and it can be estimated in a 'model-free' manner by an estimator ('ε'). We introduce ε, demonstrate its estimation using simulated data, and discuss how it may be used in planning and evaluation of ongoing elimination efforts for a range of parasitic diseases.

Substantiating freedom from parasitic infection by combining transmission model predictions with disease surveys

Stopping interventions is a critical decision for parasite elimination programmes. Quantifying the probability that elimination has occurred due to interventions can be facilitated by combining infection status information from parasitological surveys with extinction thresholds predicted by parasite transmission models. Here we demonstrate how the integrated use of these two pieces of information derived from infection monitoring data can be used to develop an analytic framework for guiding the making of defensible decisions to stop interventions. We present a computational tool to perform these probability calculations and demonstrate its practical utility for supporting intervention cessation decisions by applying the framework to infection data from programmes aiming to eliminate onchocerciasis and lymphatic filariasis in Uganda and Nigeria, respectively. We highlight a possible method for validating the results in the field, and discuss further refinements and extensions required to deploy this predictive tool for guiding decision making by programme managers.

The decision when to stop an intervention is a critical component of parasite elimination programmes, but reliance on surveillance data alone can be inaccurate. Here, Michael et al. combine parasite transmission model predictions with disease survey data to more reliably determine when interventions can be stopped.

Continental-scale, data-driven predictive assessment of eliminating the vector-borne disease, lymphatic filariasis, in sub-Saharan Africa by 2020

BACKGROUND

There are growing demands for predicting the prospects of achieving the global elimination of neglected tropical diseases as a result of the institution of large-scale nation-wide intervention programs by the WHO-set target year of 2020. Such predictions will be uncertain due to the impacts that spatial heterogeneity and scaling effects will have on parasite transmission processes, which will introduce significant aggregation errors into any attempt aiming to predict the outcomes of interventions at the broader spatial levels relevant to policy making. We describe a modeling platform that addresses this problem of upscaling from local settings to facilitate predictions at regional levels by the discovery and use of locality-specific transmission models, and we illustrate the utility of using this approach to evaluate the prospects for eliminating the vector-borne disease, lymphatic filariasis (LF), in sub-Saharan Africa by the WHO target year of 2020 using currently applied or newly proposed intervention strategies. METHODS AND RESULTS: We show how a computational platform that couples site-specific data discovery with model fitting and calibration can allow both learning of local LF transmission models and simulations of the impact of interventions that take a fuller account of the fine-scale heterogeneous transmission of this parasitic disease within endemic countries. We highlight how such a spatially hierarchical modeling tool that incorporates actual data regarding the roll-out of national drug treatment programs and spatial variability in infection patterns into the modeling process can produce more realistic predictions of timelines to LF elimination at coarse spatial scales, ranging from district to country to continental levels. Our results show that when locally applicable extinction thresholds are used, only three countries are likely to meet the goal of LF elimination by 2020 using currently applied mass drug treatments, and that switching to more intensive drug regimens, increasing the frequency of treatments, or switching to new triple drug regimens will be required if LF elimination is to be accelerated in Africa. The proportion of countries that would meet the goal of eliminating LF by 2020 may, however, reach up to 24/36 if the WHO 1% microfilaremia prevalence threshold is used and sequential mass drug deliveries are applied in countries.

CONCLUSIONS

We have developed and applied a data-driven spatially hierarchical computational platform that uses the discovery of locally applicable transmission models in order to predict the prospects for eliminating the macroparasitic disease, LF, at the coarser country level in sub-Saharan Africa. We show that fine-scale spatial heterogeneity in local parasite transmission and extinction dynamics, as well as the exact nature of intervention roll-outs in countries, will impact the timelines to achieving national LF elimination on this continent.

Assessing endgame strategies for the elimination of lymphatic filariasis: A model-based evaluation of the impact of DEC-medicated salt

Concern is growing regarding the prospects of achieving the global elimination of lymphatic filariasis (LF) by 2020. Apart from operational difficulties, evidence is emerging which points to unique challenges that could confound achieving LF elimination as extinction targets draw near. Diethylcarbamazine (DEC)-medicated salt may overcome these complex challenges posed by the endgame phase of parasite elimination. We calibrated LF transmission models using Bayesian data-model assimilation techniques to baseline and follow-up infection data from 11 communities that underwent DEC salt medication. The fitted models were used to assess the utility of DEC salt treatment for achieving LF elimination, in comparison with other current and proposed drug regimens, during the endgame phase. DEC-medicated salt consistently reduced microfilaria (mf) prevalence from 1% mf to site-specific elimination thresholds more quickly than the other investigated treatments. The application of DEC salt generally required less than one year to achieve site-specific LF elimination, while annual and biannual MDA options required significantly longer durations to achieve the same task. The use of DEC-medicated salt also lowered between-site variance in extinction timelines, especially when combined with vector control. These results indicate that the implementation of DEC-medicated salt, where feasible, can overcome endgame challenges facing LF elimination programs.

Low Transmission to Elimination: Rural Development as a Key Determinant of the End-Game Dynamics of Schistosoma japonicum in China

Rural development has been a critical component of China’s economic miracle since the start of economic reform in the early 1980s, both benefiting from and contributing to the nation’s rapid economic growth. This development has yielded substantial improvements of public health relevance, including contributing to major reductions in schistosomiasis prevalence. The history of schistosomiasis elimination in Japan suggests that development played a dominant causal role in that nation. We argue that it is highly probable that a similar story is playing out in at least some large regions of China. In particular, we summarize evidence from Sichuan Province which supports the case that economic development has led to improvements in rural irrigation and water supply which, together with changes in crop selection and agricultural mechanization, have all contributed to sustainable reductions in the prevalence of Schistosoma japonicum. The two major factors that have experienced major reductions are the area of snail habitat and the degree of human exposure, both through a variety of mechanisms which differ by region and economic circumstance. However, hotspots of transmission remain. Overall, however, economic development in traditionally endemic areas has provided the resources to carry out projects that have had major beneficial impacts on disease transmission that are likely to be sustainable.

Modeling the Parasitic Filariasis Spread by Mosquito in Periodic Environment

In this paper a mosquito-borne parasitic infection model in periodic environment is considered. Threshold parameter R0 is given by linear next infection operator, which determined the dynamic behaviors of system. We obtain that when R0 < 1, the disease-free periodic solution is globally asymptotically stable and when R0 > 1 by Poincaré map we obtain that disease is uniformly persistent. Numerical simulations support the results and sensitivity analysis shows effects of parameters on R0, which provided references to seek optimal measures to control the transmission of lymphatic filariasis.

Socio-ecological dynamics and challenges to the governance of Neglected Tropical Disease control

The current global attempts to control the so-called “Neglected Tropical Diseases (NTDs)” have the potential to significantly reduce the morbidity suffered by some of the world’s poorest communities. However, the governance of these control programmes is driven by a managerial rationality that assumes predictability of proposed interventions, and which thus primarily seeks to improve the cost-effectiveness of implementation by measuring performance in terms of pre-determined outputs. Here, we argue that this approach has reinforced the narrow normal-science model for controlling parasitic diseases, and in doing so fails to address the complex dynamics, uncertainty and socio-ecological context-specificity that invariably underlie parasite transmission. We suggest that a new governance approach is required that draws on a combination of non-equilibrium thinking about the operation of complex, adaptive, systems from the natural sciences and constructivist social science perspectives that view the accumulation of scientific knowledge as contingent on historical interests and norms, if more effective control approaches sufficiently sensitive to local disease contexts are to be devised, applied and managed. At the core of this approach is an emphasis on the need for a process that assists with the inclusion of diverse perspectives, social learning and deliberation, and a reflexive approach to addressing system complexity and incertitude, while balancing this flexibility with stability-focused structures. We derive and discuss a possible governance framework and outline an organizational structure that could be used to effectively deal with the complexity of accomplishing global NTD control. We also point to examples of complexity-based management structures that have been used in parasite control previously, which could serve as practical templates for developing similar governance structures to better manage global NTD control. Our results hold important wider implications for global health policy aiming to effectively control and eradicate parasitic diseases across the world.

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.

Modelling the health impact and cost-effectiveness of lymphatic filariasis eradication under varying levels of mass drug administration scale-up and geographic coverage

BACKGROUND

A global programme to eliminate lymphatic filariasis (GPELF) is underway, yet two key programmatic features are currently still lacking: (1) the extension of efforts to all lymphatic filariasis (LF) endemic countries, and (2) the expansion of geographic coverage of mass drug administration (MDA) within countries. For varying levels of scale-up of MDA, we assessed the health benefits and the incremental cost-effectiveness ratios (ICERs) associated with LF eradication, projected the potential savings due to decreased morbidity management needs, and estimated potential household productivity gains as a result of reduced LF-related morbidity.

METHODS

We extended an LF transmission model to track hydrocele and lymphoedema incidence in order to obtain estimates of the disability adjusted life years (DALYs) averted due to scaling up MDA over a period of 50 years. We then estimated the ICERs and the cost-effectiveness acceptability curves associated with different rates of MDA scale-up. Health systems savings were estimated by considering the averted morbidity, treatment-seeking behaviour and morbidity management costs. Gains in worker productivity were estimated by multiplying estimated working days lost as a result of morbidity with country-specific per-worker agricultural wages.

RESULTS

Our projections indicate that a massive scaling-up of MDA could lead to 4.38 million incremental DALYs averted over a 50-year time horizon compared to a scenario which mirrors current efforts against LF. In comparison to maintaining the current rate of progress against LF, massive scaling-up of MDA-pursuing LF eradication as soon as possible-was most likely to be cost-effective above a willingness to pay threshold of US$71.5/DALY averted. Intensified MDA scale-up was also associated with lower ICERs. Furthermore, this could result in health systems savings up to US$483 million. Extending coverage to all endemic areas could generate additional economic benefits through gains in worker productivity between US$3.4 and US$14.4 billion.

CONCLUSIONS

In addition to ethical and political motivations for scaling-up MDA rapidly, this analysis provides economic support for increasing the intensity of MDA programmes.

Heterogeneous dynamics, robustness/fragility trade-offs, and the eradication of the macroparasitic disease, lymphatic filariasis

BACKGROUND

The current WHO-led initiative to eradicate the macroparasitic disease, lymphatic filariasis (LF), based on single-dose annual mass drug administration (MDA) represents one of the largest health programs devised to reduce the burden of tropical diseases. However, despite the advances made in instituting large-scale MDA programs in affected countries, a challenge to meeting the goal of global eradication is the heterogeneous transmission of LF across endemic regions, and the impact that such complexity may have on the effort required to interrupt transmission in all socioecological settings.

METHODS

Here, we apply a Bayesian computer simulation procedure to fit transmission models of LF to field data assembled from 18 sites across the major LF endemic regions of Africa, Asia and Papua New Guinea, reflecting different ecological and vector characteristics, to investigate the impacts and implications of transmission heterogeneity and complexity on filarial infection dynamics, system robustness and control.

RESULTS

We find firstly that LF elimination thresholds varied significantly between the 18 study communities owing to site variations in transmission and initial ecological parameters. We highlight how this variation in thresholds lead to the need for applying variable durations of interventions across endemic communities for achieving LF elimination; however, a major new result is the finding that filarial population responses to interventions ultimately reflect outcomes of interplays between dynamics and the biological architectures and processes that generate robustness/fragility trade-offs in parasite transmission. Intervention simulations carried out in this study further show how understanding these factors is also key to the design of options that would effectively eliminate LF from all settings. In this regard, we find how including vector control into MDA programs may not only offer a countermeasure that will reliably increase system fragility globally across all settings and hence provide a control option robust to differential locality-specific transmission dynamics, but by simultaneously reducing transmission regime variability also permit more reliable macroscopic predictions of intervention effects.

CONCLUSIONS

Our results imply that a new approach, combining adaptive modelling of parasite transmission with the use of biological robustness as a design principle, is required if we are to both enhance understanding of complex parasitic infections and delineate options to facilitate their elimination effectively.

Filarial hydropneumothorax: a strange journey

Filarial infection can have varied manifestations, but hydropneumothorax at presentation has not yet been reported. A 28-year-old man presented to our hospital with heaviness of the left chest for the past 10 days, which was preceded by a sudden, short stabbing pain in the left chest after straining. Chest X-ray revealed left-sided hydropneumothorax. A peripheral blood picture revealed significant eosinophilia. A pleural fluid report also showed eosinophilia and a few motile microfilaria of Wuchereria bancrofti. Microfilaria was also documented in peripheral blood. There was no evidence of other organ system involvement. The patient was diagnosed with 'Filarial Hydropneumothorax'. After treatment with a temporary chest drain and oral diethylcarbamazine citrate, there was dramatic relief of symptoms and radiological improvement. The patient has been symptom free with no features of recurrence through 8 months of follow-up.

Modelling the distribution and transmission intensity of lymphatic filariasis in sub-Saharan Africa prior to scaling up interventions: integrated use of geostatistical and mathematical modelling

BACKGROUND

Lymphatic filariasis (LF) is one of the neglected tropical diseases targeted for global elimination. The ability to interrupt transmission is, partly, influenced by the underlying intensity of transmission and its geographical variation. This information can also help guide the design of targeted surveillance activities. The present study uses a combination of geostatistical and mathematical modelling to predict the prevalence and transmission intensity of LF prior to the implementation of large-scale control in sub-Saharan Africa.

METHODS

A systematic search of the literature was undertaken to identify surveys on the prevalence of Wuchereria bancrofti microfilaraemia (mf), based on blood smears, and on the prevalence of antigenaemia, based on the use of an immuno-chromatographic card test (ICT). Using a suite of environmental and demographic data, spatiotemporal multivariate models were fitted separately for mf prevalence and ICT-based prevalence within a Bayesian framework and used to make predictions for non-sampled areas. Maps of the dominant vector species of LF were also developed. The maps of predicted prevalence and vector distribution were linked to mathematical models of the transmission dynamics of LF to infer the intensity of transmission, quantified by the basic reproductive number (R0).

RESULTS

The literature search identified 1267 surveys that provide suitable data on the prevalence of mf and 2817 surveys that report the prevalence of antigenaemia. Distinct spatial predictions arose from the models for mf prevalence and ICT-based prevalence, with a wider geographical distribution when using ICT-based data. The vector distribution maps demonstrated the spatial variation of LF vector species. Mathematical modelling showed that the reproduction number (R0) estimates vary from 2.7 to 30, with large variations between and within regions.

CONCLUSIONS

LF transmission is highly heterogeneous, and the developed maps can help guide intervention, monitoring and surveillance strategies as countries progress towards LF elimination.

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.

Modelling strategies to break transmission of lymphatic filariasis--aggregation, adherence and vector competence greatly alter elimination

BACKGROUND

With ambitious targets to eliminate lymphatic filariasis over the coming years, there is a need to identify optimal strategies to achieve them in areas with different baseline prevalence and stages of control. Modelling can assist in identifying what data should be collected and what strategies are best for which scenarios.

METHODS

We develop a new individual-based, stochastic mathematical model of the transmission of lymphatic filariasis. We validate the model by fitting to a first time point and predicting future timepoints from surveillance data in Kenya and Sri Lanka, which have different vectors and different stages of the control programme. We then simulate different treatment scenarios in low, medium and high transmission settings, comparing once yearly mass drug administration (MDA) with more frequent MDA and higher coverage. We investigate the potential impact that vector control, systematic non-compliance and different levels of aggregation have on the dynamics of transmission and control.

RESULTS

In all settings, increasing coverage from 65 to 80 % has a similar impact on control to treating twice a year at 65 % coverage, for fewer drug treatments being distributed. Vector control has a large impact, even at moderate levels. The extent of aggregation of parasite loads amongst a small portion of the population, which has been estimated to be highly variable in different settings, can undermine the success of a programme, particularly if high risk sub-communities are not accessing interventions.

CONCLUSION

Even moderate levels of vector control have a large impact both on the reduction in prevalence and the maintenance of gains made during MDA, even when parasite loads are highly aggregated, and use of vector control is at moderate levels. For the same prevalence, differences in aggregation and adherence can result in very different dynamics. The novel analysis of a small amount of surveillance data and resulting simulations highlight the need for more individual level data to be analysed to effectively tailor programmes in the drive for elimination.

What Is Needed to Eradicate Lymphatic Filariasis? A Model-Based Assessment on the Impact of Scaling Up Mass Drug Administration Programs

Background

Lymphatic filariasis (LF) is a neglected tropical disease for which more than a billion people in 73 countries are thought to be at-risk. At a global level, the efforts against LF are designed as an elimination program. However, current efforts appear to aim for elimination in some but not all endemic areas. With the 2020 goal of elimination looming, we set out to develop plausible scale-up scenarios to reach global elimination and eradication. We predict the duration of mass drug administration (MDA) necessary to reach local elimination for a variety of transmission archetypes using an existing model of LF transmission, estimate the number of treatments required for each scenario, and consider implications of rapid scale-up.

Methodology

We have defined four scenarios that differ in their geographic coverage and rate of scale-up. For each scenario, country-specific simulations and calculations were performed that took into account the pre-intervention transmission intensity, the different vector genera, drug regimen, achieved level of population coverage, previous progress toward elimination, and potential programmatic delays due to mapping, operations, and administration.

Principal Findings

Our results indicate that eliminating LF by 2020 is unlikely. If MDA programs are drastically scaled up and expanded, the final round of MDA for LF eradication could be delivered in 2028 after 4,159 million treatments. However, if the current rate of scale-up is maintained, the final round of MDA to eradicate LF may not occur until 2050.

Conclusions/Significance

Rapid scale-up of MDA will decrease the amount of time and treatments required to reach LF eradication. It may also propel the program towards success, as the risk of failure is likely to increase with extended program duration.

Lymphatic filariasis (LF) is a disease caused by filarial worms transmitted by different types of mosquitos that can lead to massive disability, including elephantiasis and hydrocele. LF has no significant zoonotic reservoir and is thought to be a potentially eradicable disease through once yearly treatment distributed by mass drug administration (MDA). In this study, we set out to determine how many treatments and over how much time it might take to globally eliminate and eradicate LF under different levels of treatment intensities. We created a model that took into account country-specific and disease-specific variables, and found that if the current intensity of MDA is maintained, 3,409 million treatments distributed over the next 37 years will be required. However, if treatment is rapidly expanded to the entire at-risk population in all endemic countries, eradication could be achieved with 4,159 million treatments and in less than half the time. While our estimates suggest more time may be needed to reach LF elimination than what is currently projected, with continued commitment, eradicating LF is within reach.

Venom allergen 5 is Associated With Deltamethrin Resistance in Culex pipiens pallens (Diptera: Culicidae)

The mosquito, Culex pipiens pallens (L.), is an important vector of encephalitis and filariasis in northern China. The control of these mosquitoes occurs primarily via the use of pyrethroid insecticides, such as deltamethrin. The widespread and improper application of pyrethroid has resulted in the evolution of pyrethroid resistance amongst many mosquito populations, including Cx. pipiens pallens. Previous studies using high-throughput transcriptome sequencing have identified that the venom allergen 5 gene is differentially expressed between deltamethrin-susceptible and deltamethrin-resistant Cx. pipiens pallens. In this study, quantitative real-time polymerase chain reaction analyses revealed that venom allergen 5 was significantly overexpressed in adult females of both deltamethrin-resistant laboratory populations and two field populations. The transcriptional level of venom allergen 5 in the laboratory populations was elevated as the levels of deltamethrin resistance increased. Full-length cDNAs of the venom allergen 5 gene were cloned from Cx. pipiens pallens, and contained an open reading frame of 765 bp, encoding a protein with 254 amino acids. The deduced amino acid sequence shared 100% identity with the ortholog in Culex quinquefasciatus Say. The overexpression of venom allergen 5 decreased the susceptibility of mosquito cells to deltamethrin, while knockdown of this gene by RNAi increased the susceptibility of mosquitoes to deltamethrin. This study provides the first evidence of the association between the venom allergen 5 gene and deltamethrin resistance in mosquitoes.

Exploring the contribution of host susceptibility to epidemiological patterns of Schistosoma japonicum infection using an individual-based model

We recently reported the analysis of epidemiological data suggesting variability in individual susceptibility to infection by Schistosoma japonicum among rural villagers who reside in Sichuan Province of southwestern China. By supplementing the data used in the earlier analysis from other studies we have reported from this region, we presented improved estimates of cercarial exposure, which in turn, result in stronger evidence of susceptibility. This analysis was conducted using an individual-based mathematical model (IBM) whose use was motivated by the nature and extent of field data from the low-transmission environments exemplified by one of our datasets and typical of the current situation in most endemic areas of China. In addition to individual susceptibility and water contact, the model includes stochastic aspects of cercarial exposure as well as of diagnostic procedures, the latter being particularly relevant to the low-transmission environment. The simulation studies show that, to produce key aspects of the epidemiological findings, the distribution of susceptibility ranges over several orders of magnitude and is highly right skewed. We found no compelling evidence that the distribution of susceptibility differed between the two populations that underlie both the epidemiological and simulation results.

The Allee effect and elimination of neglected tropical diseases: a mathematical modelling study

Elimination and control programmes for neglected tropical diseases (NTDs) are underway around the world, yet they are generally informed by epidemiological modelling only to a rudimentary degree. Chief among the modelling-derived predictors of disease emergence or controllability is the basic reproduction number R0. The ecological systems of several of the NTDs include density-dependent processes--which alter the rate of e.g. parasite establishment or fecundity--that complicate the calculation of R0. Here we show how the forms of the density-dependent functions for a model of the NTD lymphatic filariasis affect the effective reproduction number Reff. We construct infection transmission models containing various density-dependent functions and show how they alter the shape of the Reff profile, affecting two important epidemiological outcome variables that relate to elimination and control programmes: the parasite transmission breakpoint (or extinction threshold) and the reproduction fitness, as measured by Reff. The current drive to control, eliminate or eradicate several parasitic infections would be substantially aided by the existence of ecological Allee effects. For these control programmes, the findings of this paper are encouraging, since a single positive density dependency (DD) can introduce a reasonable chance of achieving elimination; however, there are diminishing returns to additional positive DDs.

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.

How effective is integrated vector management against malaria and lymphatic filariasis where the diseases are transmitted by the same vector?

BACKGROUND

The opportunity to integrate vector management across multiple vector-borne diseases is particularly plausible for malaria and lymphatic filariasis (LF) control where both diseases are transmitted by the same vector. To date most examples of integrated control targeting these diseases have been unanticipated consequences of malaria vector control, rather than planned strategies that aim to maximize the efficacy and take the complex ecological and biological interactions between the two diseases into account.

METHODOLOGY/PRINCIPAL FINDINGS

We developed a general model of malaria and LF transmission and derived expressions for the basic reproductive number (R0) for each disease. Transmission of both diseases was most sensitive to vector mortality and biting rate. Simulating different levels of coverage of long lasting-insecticidal nets (LLINs) and larval control confirms the effectiveness of these interventions for the control of both diseases. When LF was maintained near the critical density of mosquitoes, minor levels of vector control (8% coverage of LLINs or treatment of 20% of larval sites) were sufficient to eliminate the disease. Malaria had a far greater R0 and required a 90% population coverage of LLINs in order to eliminate it. When the mosquito density was doubled, 36% and 58% coverage of LLINs and larval control, respectively, were required for LF elimination; and malaria elimination was possible with a combined coverage of 78% of LLINs and larval control.

CONCLUSIONS/SIGNIFICANCE

Despite the low level of vector control required to eliminate LF, simulations suggest that prevalence of LF will decrease at a slower rate than malaria, even at high levels of coverage. If representative of field situations, integrated management should take into account not only how malaria control can facilitate filariasis elimination, but strike a balance between the high levels of coverage of (multiple) interventions required for malaria with the long duration predicted to be required for filariasis elimination.

The global distribution and transmission limits of lymphatic filariasis: past and present

BACKGROUND

Lymphatic filariasis (LF) is one of the neglected tropical diseases targeted for global elimination by 2020 and to guide elimination efforts countries have, in recent years, conducted extensive mapping surveys. Documenting the past and present distribution of LF and its environmental limits is important for a number of reasons. Here, we present an initiative to develop a global atlas of LF and present a new global map of the limits of LF transmission.

METHODS

We undertook a systematic search and assembly of prevalence data worldwide and used a suite of environmental and climatic data and boosted regression trees (BRT) modelling to map the transmission limits of LF.

RESULTS

Data were identified for 66 of the 72 countries currently endemic and for a further 17 countries where LF is no longer endemic. Our map highlights a restricted and highly heterogeneous distribution in sub-Saharan Africa, with transmission more widespread in West Africa compared to east, central and southern Africa where pockets of transmission occur. Contemporary transmission occurs across much of south and South-east Asia and the Pacific. Interestingly, the risk map reflects environmental conditions suitable for LF transmission across Central and South America, including the southern States of America, although active transmission is only known in a few isolated foci. In countries that have eliminated LF, our predictions of environmental suitability are consistent with historical distribution.

CONCLUSIONS

The global distribution of LF is highly heterogeneous and geographically targeted and sustained control will be required to achieve elimination. This first global map can help evaluate the progress of interventions and guide surveillance activities.

Exploring the impact of infection-induced immunity on the transmission of Schistosoma japonicum in hilly and mountainous environments in China

Schistosomiasis has long been a threat to villagers in hilly and mountainous areas of southwestern China where the intermediate snail host is abundant. In recent years our group has focused on the development and parameterization of a community-level mathematical model of S. japonicum transmission that accounts for the role of environmental determinants of transmission intensity in Sichuan Province. To date the model has not incorporated acquired immunity. A review of previous epidemiologic data from our study area in Sichuan suggested modeling of acquired immunity as a function of history of infection. To explore the potential impacts on the dynamics of transmission, a mathematical representation of acquired immunity was incorporated, and parameterized based on this epidemiological evidence. It is shown through simulation that the effect of immunity is to reduce the rate of worm development and thereby lower the endemic level significantly. The effect was more striking at increasing levels of a village'tm)s basic reproductive number. Further, residual immunity modestly alters the threshold of external parasite input necessary to trigger re-emergence of transmission and its subsequent rate of development. Despite limitations in our quantitative knowledge of the immunity function, these findings, along with the uncertainties in transmission dynamics at low infection levels, underscore the need for improved diagnostic methods for disease control, especially in potentially re-emergent settings.

Recasting the theory of mosquito-borne pathogen transmission dynamics and control

Mosquito-borne diseases pose some of the greatest challenges in public health, especially in tropical and sub-tropical regions of the world. Efforts to control these diseases have been underpinned by a theoretical framework developed for malaria by Ross and Macdonald, including models, metrics for measuring transmission, and theory of control that identifies key vulnerabilities in the transmission cycle. That framework, especially Macdonald's formula for R₀ and its entomological derivative, vectorial capacity, are now used to study dynamics and design interventions for many mosquito-borne diseases. A systematic review of 388 models published between 1970 and 2010 found that the vast majority adopted the Ross–Macdonald assumption of homogeneous transmission in a well-mixed population. Studies comparing models and data question these assumptions and point to the capacity to model heterogeneous, focal transmission as the most important but relatively unexplored component in current theory. Fine-scale heterogeneity causes transmission dynamics to be nonlinear, and poses problems for modeling, epidemiology and measurement. Novel mathematical approaches show how heterogeneity arises from the biology and the landscape on which the processes of mosquito biting and pathogen transmission unfold. Emerging theory focuses attention on the ecological and social context for mosquito blood feeding, the movement of both hosts and mosquitoes, and the relevant spatial scales for measuring transmission and for modeling dynamics and control.

Insecticidal bed nets and filariasis transmission in Papua New Guinea

BACKGROUND

Global efforts to eliminate lymphatic filariasis are based on the annual mass administration of antifilarial drugs to reduce the microfilaria reservoir available to the mosquito vector. Insecticide-treated bed nets are being widely used in areas in which filariasis and malaria are coendemic.

METHODS

We studied five villages in which five annual mass administrations of antifilarial drugs, which were completed in 1998, reduced the transmission of Wuchereria bancrofti, one of the nematodes that cause lymphatic filariasis. A total of 21,899 anopheles mosquitoes were collected for 26 months before and 11 to 36 months after bed nets treated with long-lasting insecticide were distributed in 2009. We evaluated the status of filarial infection and the presence of W. bancrofti DNA in anopheline mosquitoes before and after the introduction of insecticide-treated bed nets. We then used a model of population dynamics to estimate the probabilities of transmission cessation.

RESULTS

Village-specific rates of bites from anopheline mosquitoes ranged from 6.4 to 61.3 bites per person per day before the bed-net distribution and from 1.1 to 9.4 bites for 11 months after distribution (P<0.001). During the same period, the rate of detection of W. bancrofti in anopheline mosquitoes decreased from 1.8% to 0.4% (P=0.005), and the rate of detection of filarial DNA decreased from 19.4% to 14.9% (P=0.13). The annual transmission potential was 5 to 325 infective larvae inoculated per person per year before the bed-net distribution and 0 after the distribution. Among all five villages with a prevalence of microfilariae of 2 to 38%, the probability of transmission cessation increased from less than 1.0% before the bed-net distribution to a range of 4.9 to 95% in the 11 months after distribution.

CONCLUSIONS

Vector control with insecticide-treated bed nets is a valuable tool for W. bancrofti elimination in areas in which anopheline mosquitoes transmit the parasite. (Funded by the U.S. Public Health Service and the National Institutes of Health.).

Mapping, bayesian geostatistical analysis and spatial prediction of lymphatic filariasis prevalence in Africa

There is increasing interest to control or eradicate the major neglected tropical diseases. Accurate modelling of the geographic distributions of parasitic infections will be crucial to this endeavour. We used 664 community level infection prevalence data collated from the published literature in conjunction with eight environmental variables, altitude and population density, and a multivariate Bayesian generalized linear spatial model that allows explicit accounting for spatial autocorrelation and incorporation of uncertainty in input data and model parameters, to construct the first spatially-explicit map describing LF prevalence distribution in Africa. We also ran the best-fit model against predictions made by the HADCM3 and CCCMA climate models for 2050 to predict the likely distributions of LF under future climate and population changes. We show that LF prevalence is strongly influenced by spatial autocorrelation between locations but is only weakly associated with environmental covariates. Infection prevalence, however, is found to be related to variations in population density. All associations with key environmental/demographic variables appear to be complex and non-linear. LF prevalence is predicted to be highly heterogenous across Africa, with high prevalences (>20%) estimated to occur primarily along coastal West and East Africa, and lowest prevalences predicted for the central part of the continent. Error maps, however, indicate a need for further surveys to overcome problems with data scarcity in the latter and other regions. Analysis of future changes in prevalence indicates that population growth rather than climate change per se will represent the dominant factor in the predicted increase/decrease and spread of LF on the continent. We indicate that these results could play an important role in aiding the development of strategies that are best able to achieve the goals of parasite elimination locally and globally in a manner that may also account for the effects of future climate change on parasitic infection.

Sequential modelling of the effects of mass drug treatments on anopheline-mediated lymphatic filariasis infection in Papua New Guinea

BACKGROUND

Lymphatic filariasis (LF) has been targeted by the WHO for global eradication leading to the implementation of large scale intervention programs based on annual mass drug administrations (MDA) worldwide. Recent work has indicated that locality-specific bio-ecological complexities affecting parasite transmission may complicate the prediction of LF extinction endpoints, casting uncertainty on the achievement of this initiative. One source of difficulty is the limited quantity and quality of data used to parameterize models of parasite transmission, implying the important need to update initially-derived parameter values. Sequential analysis of longitudinal data following annual MDAs will also be important to gaining new understanding of the persistence dynamics of LF. Here, we apply a Bayesian statistical-dynamical modelling framework that enables assimilation of information in human infection data recorded from communities in Papua New Guinea that underwent annual MDAs, into our previously developed model of parasite transmission, in order to examine these questions in LF ecology and control.

RESULTS

Biological parameters underlying transmission obtained by fitting the model to longitudinal data remained stable throughout the study period. This enabled us to reliably reconstruct the observed baseline data in each community. Endpoint estimates also showed little variation. However, the updating procedure showed a shift towards higher and less variable values for worm kill but not for any other drug-related parameters. An intriguing finding is that the stability in key biological parameters could be disrupted by a significant reduction in the vector biting rate prevailing in a locality.

CONCLUSIONS

Temporal invariance of biological parameters in the face of intervention perturbations indicates a robust adaptation of LF transmission to local ecological conditions. The results imply that understanding the mechanisms that underlie locally adapted transmission dynamics will be integral to identifying points of system fragility, and thus countermeasures to reliably facilitate LF extinction both locally and globally.

Modelling co-infection with malaria and lymphatic filariasis

Malaria and lymphatic filariasis (LF) continue to cause a considerable public health burden globally and are co-endemic in many regions of sub-Saharan Africa. These infections are transmitted by the same mosquito species which raises important questions about optimal vector control strategies in co-endemic regions, as well as the effect of the presence of each infection on endemicity of the other; there is currently little consensus on the latter. The need for comprehensive modelling studies to address such questions is therefore significant, yet very few have been undertaken to date despite the recognised explanatory power of reliable dynamic mathematical models. Here, we develop a malaria-LF co-infection modelling framework that accounts for two key interactions between these infections, namely the increase in vector mortality as LF mosquito prevalence increases and the antagonistic Th1/Th2 immune response that occurs in co-infected hosts. We consider the crucial interplay between these interactions on the resulting endemic prevalence when introducing each infection in regions where the other is already endemic (e.g. due to regional environmental change), and the associated timescale for such changes, as well as effects on the basic reproduction number R₀ of each disease. We also highlight potential perverse effects of vector controls on human infection prevalence in co-endemic regions, noting that understanding such effects is critical in designing optimal integrated control programmes. Hence, as well as highlighting where better data are required to more reliably address such questions, we provide an important framework that will form the basis of future scenario analysis tools used to plan and inform policy decisions on intervention measures in different transmission settings.

Field transmission intensity of Schistosoma japonicum measured by basic reproduction ratio from modified Barbour's model

BACKGROUND

Schistosomiasis japonica, caused by infection with Schistosoma japonicum, is still recognized as a major public health problem in the Peoples' Republic of China. Mathematical modelling of schistosomiasis transmission has been undertaken in order to assess and project the effects of various control strategies for elimination of the disease. Seasonal fluctuations in transmission may have the potential to impact on the population dynamics of schistosomiasis, yet no model of S. japonicum has considered such effects. In this paper, we characterize the transmission dynamics of S. japonicum using a modified version of Barbour's model to account for seasonal variation (SV), and investigate the effectiveness of the control strategy adopted in Liaonan village of Xingzi county, Jiangxi Province.

METHODS

We use mathematical tools for stability analysis of periodic systems and derive expressions for the basic reproduction ratio of S. japonicum in humans; we parameterise such expressions with surveillance data to investigate the conditions for persistence or elimination of the disease in the study village. We perform numerical simulations and parametric sensitivity analysis to understand local transmission conditions and compare values of the basic reproductive ratio with and without seasonal fluctuations.

RESULTS

The explicit formula of the basic reproduction ratio for the SV-modified Barbour's model is derived. Results show that the value of the basic reproduction ratio, R0, of Liaonan village, Xingzi county is located between 1.064 and 1.066 (very close to 1), for schistosomiasis transmission during 2006 to 2010, after intensification of control efforts.

CONCLUSIONS

Our modified version of the Barbour model to account for seasonal fluctuations in transmission has the potential to provide better estimations of infection risk than previous models. Ignoring seasonality tends to underestimate R0 values albeit only marginally. In the absence of simultaneous R0 estimations for villages not under control interventions (such villages do not currently exist in China), it is difficult to assess whether control strategies have had a substantial impact on levels of transmission, as the parasite population would still be able to maintain itself at an endemic level, highlighting the difficulties faced by elimination efforts.

Filariasis in an infant with B-cell acute lymphoblastic leukemia: a rare enigma

Filariasis, a tropical parasitic infection, is a common public health problem in the Indian sub-continent. Occurrence of filariasis with acute lymphoblastic leukemia (ALL) is unusual, particularly in infants, since filariasis has a long incubation period of about 1 year. Though, there are case reports of leishmaniasis, malaria and other vector borne diseases seen in association with leukemias, filariasis co-existing with ALL has not been documented to the best of our knowledge. We report an incidental finding of Bancroftian filariasis in an 11 month old female already diagnosed as B-cell ALL.

Functional characterization of an arrestin gene on insecticide resistance of Culex pipiens pallens

Background

Continuous and excessive application of insecticides has resulted in the rapid development of insecticide resistance in several mosquito species, including Culex pipiens pallens. Previous studies in our laboratory found that arrestin gene expression was higher in the deltamethrin-resistant (DR) strain than in the deltamethrin-susceptible (DS) strain of Cx. pipiens pallens. Similarly, other studies reported that arrestin was highly expressed in permethrin-resistant Cx. quinquefasciatus and in dichlorodiphenyltrichloroethane (DDT)-resistant Drosophila melanogaster.

Methods

Full-length cDNAs of an arrestin gene were cloned from Cx. pipiens pallens via polymerase chain reaction (PCR) and rapid amplification of cDNA end (RACE). The mRNA levels of the arrestin gene in the whole life cycle of DR and DS strains of Cx. pipiens pallens were investigated via quantitative real-time PCR. In addition, the relationship between arrestin and deltamethrin (DM) resistance were identified using genetic overexpression strategies and arrestin RNAi in mosquito cells. Cell viability was analyzed with cholecystokinin octapeptide after DM treatment. Moreover, the mRNA levels of cytochrome P450 6A1 (CYP6A1) and opsin in the transfected cells and controls were analyzed.

Results

Complete arrestin gene sequence was cloned and expressed throughout the life cycle of Cx. pipiens pallens. Moreover, arrestin was significantly upregulated in the DR strain, compared with that in the DS strain at the egg, pupae, and adult stages. Arrestin overexpression comparably increased the mosquito cell viability, whereas arrestin knockdown by siRNA decreased mosquito cell viability with deltamethrin (DM) treatment. Meanwhile, the mRNA levels of CYP6A1 and opsin were upregulated in mosquito cells transfected with arrestin and downregulated in mosquito cells with arrestin knockdown.

Conclusion

This study presented the first evidence that arrestin might be associated with insecticide resistance in Cx. pipiens pallens.

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

Mathematical modelling of helminth infections has the potential to inform policy and guide research for the control and elimination of human helminthiases. However, this potential, unlike in other parasitic and infectious diseases, has yet to be realised. To place contemporary efforts in a historical context, a summary of the development of mathematical models for helminthiases is presented. These efforts are discussed according to the role that models can play in furthering our understanding of parasite population biology and transmission dynamics, and the effect on such dynamics of control interventions, as well as in enabling estimation of directly unobservable parameters, exploration of transmission breakpoints, and investigation of evolutionary outcomes of control. The Disease Reference Group on Helminth Infections (DRG4), established in 2009 by the Special Programme for Research and Training in Tropical Diseases (TDR), was given the mandate to review helminthiases research and identify research priorities and gaps. A research and development agenda for helminthiasis modelling is proposed based on identified gaps that need to be addressed for models to become useful decision tools that can support research and control operations effectively. This agenda includes the use of models to estimate the impact of large-scale interventions on infection incidence; the design of sampling protocols for the monitoring and evaluation of integrated control programmes; the modelling of co-infections; the investigation of the dynamical relationship between infection and morbidity indicators; the improvement of analytical methods for the quantification of anthelmintic efficacy and resistance; the determination of programme endpoints; the linking of dynamical helminth models with helminth geostatistical mapping; and the investigation of the impact of climate change on human helminthiases. It is concluded that modelling should be embedded in helminth research, and in the planning, evaluation, and surveillance of interventions from the outset. Modellers should be essential members of interdisciplinary teams, propitiating a continuous dialogue with end users and stakeholders to reflect public health needs in the terrain, discuss the scope and limitations of models, and update biological assumptions and model outputs regularly. It is highlighted that to reach these goals, a collaborative framework must be developed for the collation, annotation, and sharing of databases from large-scale anthelmintic control programmes, and that helminth modellers should join efforts to tackle key questions in helminth epidemiology and control through the sharing of such databases, and by using diverse, yet complementary, modelling approaches.

The challenge of effective surveillance in moving from low transmission to elimination of schistosomiasis in China

Until recently, intensified efforts in China to suppress the transmission of Schistosoma japonicum relied principally on routine praziquantel treatment, extensive use of molluscicides and health education programs. These efforts, now supplemented by a broader range of control measures, have been quite successful in reducing the prevalence and intensity of human infection to very low levels. However, re-emergent transmission has occurred in formerly endemic areas of several provinces, signalling the need for more locally effective, integrated control strategies. We argue that these low but persistent levels of transmission also require important changes in both the tactics and strategy of disease surveillance to move forward towards elimination. Here we present recent data exemplifying the low transmission environment which suggests that we are reaching limits of detection of current diagnostic techniques used for human infection surveillance in these communities. However, both epidemiological data and theoretical results indicate that (i) transmission in the human population can persist at very low infection intensities even in the presence of routine control activities; (ii) the parasite can be reintroduced into parasite-free environments by very modest external inputs; and (iii) transmission at these low infection intensities exhibits very slow inter-year dynamics. These observations motivate the need for new, sensitive tools to identify low-level infections in mammalian or snail hosts, or the presence of S. japonicum in environmental media. Environmental monitoring offers an alternative, and perhaps more efficient, approach to large-scale surveillance of human infections in low transmission regions.

Neglected tropical diseases: moving beyond mass drug treatment to understanding the science

Neglected tropical diseases (NTDs) represent a major health burden in many developing countries. To date, global efforts to control thirteen parasitic and bacterial infections that affect more than 1.4 billion people have relied on mass drug administration. This singular approach should now be expanded to a more comprehensive suite of tools including coordinated community-based programs, vector control, local training, education, and environmental change. In addition, an intensive basic research agenda is urgently needed to develop effective diagnostic, preventive, and therapeutic interventions to stay one step ahead of the evolutionary adaptation tactics of disease-causing microbes and parasites.