Modelling variability in lymphatic filariasis: macrofilarial dynamics in the Brugia pahangi--cat model

Sensitivity analysis and parameters estimation for the transmission of lymphatic filariasis

Lymphatic filariasis is a neglected tropical disease which poses public health concern and socio-economic challenges in developing and low-income countries. In this paper, we formulate a deterministic mathematical model for transmission dynamics of lymphatic filariasis to generate data by white noise and use least square method to estimate parameter values. The validity of estimated parameter values is tested by Gaussian distribution method. The residuals of model outputs are normally distributed and hence can be used to study the dynamics of Lymphatic filariasis. After deriving the basic reproduction number, R 0 by the next generation matrix approach, the Partial Rank Correlation Coefficient is employed to explore which parameters significantly affect and most influential to the model outputs. The analysis for equilibrium states shows that the Lymphatic free equilibrium is globally asymptotically stable when the basic reproduction number is less a unity and endemic equilibrium is globally asymptotically stable when R 0 ≥ 1 . The findings reveal that rate of human infection, recruitment rate of mosquitoes increase the average new infections for Lymphatic filariasis. Moreover, asymptomatic individuals contribute significantly in the transmission of Lymphatic filariasis.

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.

Pathogen growth in insect hosts: inferring the importance of different mechanisms using stochastic models and response-time data

Pathogen population dynamics within individual hosts can alter disease epidemics and pathogen evolution, but our understanding of the mechanisms driving within-host dynamics is weak. Mathematical models have provided useful insights, but existing models have only rarely been subjected to rigorous tests, and their reliability is therefore open to question. Most models assume that initial pathogen population sizes are so large that stochastic effects due to small population sizes, so-called demographic stochasticity, are negligible, but whether this assumption is reasonable is unknown. Most models also assume that the dynamic effects of a host's immune system strongly affect pathogen incubation times or "response times," but whether such effects are important in real host-pathogen interactions is likewise unknown. Here we use data for a baculovirus of the gypsy moth to test models of within-host pathogen growth. By using Bayesian statistical techniques and formal model-selection procedures, we are able to show that the response time of the gypsy moth virus is strongly affected by both demographic stochasticity and a dynamic response of the host immune system. Our results imply that not all response-time variability can be explained by host and pathogen variability, and that immune system responses to infection may have important effects on population-level disease dynamics.

Bayesian experimental design for models with intractable likelihoods

In this paper we present a methodology for designing experiments for efficiently estimating the parameters of models with computationally intractable likelihoods. The approach combines a commonly used methodology for robust experimental design, based on Markov chain Monte Carlo sampling, with approximate Bayesian computation (ABC) to ensure that no likelihood evaluations are required. The utility function considered for precise parameter estimation is based upon the precision of the ABC posterior distribution, which we form efficiently via the ABC rejection algorithm based on pre-computed model simulations. Our focus is on stochastic models and, in particular, we investigate the methodology for Markov process models of epidemics and macroparasite population evolution. The macroparasite example involves a multivariate process and we assess the loss of information from not observing all variables.

Population genetics of concurrent selection with albendazole and ivermectin or diethylcarbamazine on the possible spread of albendazole resistance in Wuchereria bancrofti

The Global Program for the Elimination of Lymphatic Filariasis (GPELF) intends to achieve its aims through yearly mass treatments with albendazole (ABZ) combined with ivermectin (IVM) or diethylcarbamazine (DEC). The use of ABZ and IVM separately to combat parasites of veterinary importance has, on many occasions, resulted in widespread drug resistance. In order to help predict the spread of potential ABZ resistance alleles through a population of Wuchereria bancrofti, we have developed a mathematical model that incorporates population genetics into EPIFIL, a model which examines the transmission dynamics of the parasite. Our model considers the effect of the combined treatments on the frequency of a recessive allele, which confers ABZ resistance. The model predicts that after 10 yearly treatments with ALB and DEC, 85% coverage and an initial resistance allele frequency of 5%, the frequency of the resistance genotype will increase from 0.25 to 12.7%. If non-random mating is assumed, the initial genotype frequency will be 2.34% and will increase to 62.7%. ABZ and IVM combination treatment may lead to weaker selection for this genotype. Treatment coverage, initial allele frequencies and number of treatments also affect the rate of selection.

Prevalence, intensity and risk factors for soil-transmitted helminth infection in a South Indian fishing village

A study of the prevalence, intensity and risk factors for soil-transmitted helminth infection was undertaken among school children aged 5-9 years attending a primary school in the fishing village in Peda Jalaripet, Visakhapatnam, South India. One hundred and eighty nine (92.6%) of 204 children were infected with one or more soil transmitted helminth parasites. The predominant parasite was Ascaris lumbricoides (prevalence of 91%), followed by Trichuris trichiura (72%) and hookworm (54%). Study of age-specific prevalence and intensity of infection revealed that the prevalence and intensity of A. lumbricoides infection was higher among younger children than older children. While aggregation of parasite infection was observed, hookworm infection was more highly aggregated than either A. lumbricoides or T. trichiura. Multivariate analysis identified parental occupation, child's age and mother's education as the potential risk factors contributing to the high intensity of A. lumbricoides infection. Children from fishing families with low levels of education of the mother had the highest intensity of A. lumbricoides infection. As the outcome of chemotherapy programs to control soil transmitted helminth infection is dependant on the dynamics of their transmission, there is a need for further studies to better define the role of specific factors that determine their prevalence, intensity and aggregation in different epidemiological settings.

The dynamics ofWuchereria bancroftiinfection: a model-based analysis of longitudinal data from Pondicherry, India

This paper presents a model-based analysis of longitudinal data describing the impact of integrated vector management on the intensity ofWuchereria bancroftiinfection in Pondicherry, India. The aims of this analysis were (1) to gain insight into the dynamics of infection, with emphasis on the possible role of immunity, and (2) to develop a model that can be used to predict the effects of control. Using the LYMFASIM computer simulation program, two models with different types of immunity (anti-L3 larvae or anti-adult worm fecundity) were compared with a model without immunity. Parameters were estimated by fitting the models to data from 5071 individuals with microfilaria-density measurement before and after cessation of a 5-year vector management programme. A good fit, in particular of the convex shape of the age-prevalence curve, required inclusion of anti-L3 or anti-fecundity immunity in the model. An individual's immune-responsiveness was found to halve in ~10 years after cessation of boosting. Explanation of the large variation in Mf-density required considerable variation between individuals in exposure and immune responsiveness. The mean life-span of the parasite was estimated at about 10 years. For the post-control period, the models predict a further decline in Mf prevalence, which agrees well with observations made 3 and 6 years after cessation of the integrated vector management programme.

Robust parameter estimation techniques for stochastic within-host macroparasite models

We present a stochastic model of the within-host population dynamics of lymphatic filariasis, and use a simulated goodness-of-fit (GOF) method to estimate immunological parameters and their confidence intervals from experimental data. A variety of deterministic moment closure approximations to the stochastic system are explored and compared with simulation results. For the maximum GOF parameter estimates, none of the methods of closure accurately reproduce the behaviour of the stochastic model. However, direct analysis of the stochastic model demonstrates that the high levels of variation observed in the data can be reproduced without requiring parameters to vary between hosts. This indicates that the observed aggregation of parasite load may be dynamically generated by random variation in the development of an effective immune response against parasite larvae.

Spatial parasite transmission, drug resistance, and the spread of rare genes

The transmission of many parasitic worms involves aggregated movement between hosts of "packets" of infectious larvae. We use a generic metapopulation model to show that this aggregation naturally promotes the preferential spread of rare recessive genes, compared with the expectations of traditional nonspatial models. A more biologically realistic model also demonstrates that this effect could explain the rapid observed spread of recessive or weakly dominant drug-resistant genotypes in nematode parasites of sheep. This promotion of a recessive trait arises from a novel mechanism of inbreeding arising from the metapopulation dynamics of transmission.

Immunity, antigenic heterogeneity, and aggregation of helminth parasites

Empirical studies of helminth parasites reveal that the distribution of parasite burdens in their host populations is highly aggregated. This aggregation is fundamental to the ecology and epidemiology of helminth parasites. Results from a stochastic model predict that aggregation of helminth parasites is inversely related to the intensity of host immunity. Aggregation also decreases with antigenic heterogeneity and increases with heterogeneity in transmissibility among parasite strains. It is also found that the degree of aggregation is greater when immunity affects parasite fecundity than when immunity acts on host susceptibility. Potential relevance of this result for assessing the influence of vaccines that target either host susceptibility or parasite fecundity on the level of aggregation and consequent effects on drug resistance and disease prevalence are discussed.

Transmission dynamics of lymphatic filariasis: density-dependence in the uptake of Wuchereria bancrofti microfilariae by vector mosquitoes

Gaining a better understanding of parasite infection dynamics in the vector mosquito (Diptera: Culicidae) population is central to improving knowledge regarding the transmission, persistence and hence control of lymphatic filariasis. Here, we use data on mosquito feeding experiments collated from the published literature to examine the available evidence regarding the functional form of the first component of this parasite-vector relationship for Wuchereria bancrofti (Filarioidea: Onchocercidae) causing Bancroftian filariasis, i.e. the rate of microfilariae (mf) uptake from the blood of infected humans by the feeding mosquito vector. Using a simple logarithmic regression model for describing the observed relationships between the mean numbers of mf ingested per mosquito and parasite load in humans in each study, and a linear mixed-effects meta-analytical framework for synthesizing the observed regressions across studies, we show here for the first time clear evidence for the existence of density-dependence in this process for all the three major filariasis transmitting mosquito vectors. An important finding of this study is that this regulation of mf uptake also varies significantly between the vector genera, being weakest in Culex, comparatively stronger in Aedes and most severe and occurring at significantly lower human mf loads in Anopheles mosquitoes. The analysis of the corresponding mf uptake prevalence data has further highlighted how density-dependence in mf uptake may influence the observed distributions of mf in vector populations. These results show that whereas strong regulation of mf uptake, especially when it leads to saturation in uptake at low human parasite intensities, can lead to static distributions of mf per mosquito with host parasite intensity, a weaker regulation of mf ingestion can give rise to changes in both mean mf loads and in the frequency distribution of parasites/mosquito with increasing human parasite intensity. These findings highlight the importance of considering local vector infection dynamics when attempting to predict the impacts of community-based filariasis control. They also emphasize the value of developing and applying robust meta-analytic methods for estimating functional relationships regarding parasitic infection from population ecological data.

Models for the population biology and control of human onchocerciasis

The absence of animal models in which to reproduce successfully the complete life cycle of Onchocerca volvulus has hindered progress towards unravelling the processes involved in the regulation of parasite abundance in the vertebrate host. Mathematical frameworks have been developed to explore the consequences of such processes in determining parasite population dynamics and the effect on these of control interventions. Post-control predictions are strongly influenced by the assumptions concerning the reproductive life span of the adult female worm (the longest-lived parasite stage) and the distribution of its survival times, and this notion is important to all frameworks. Here, we review the development of models concerning onchocerciasis and discuss the various approaches that have been used, presenting a deterministic framework with parameter values estimated from the Mexican onchocerciasis control programme. This model is used to evaluate interventions combining the removal of adult worms (nodulectomy) and the microfilaricidal and possibly sterilizing effect of ivermectin.

Transmission intensity and the immunoepidemiology of bancroftian filariasis in East Africa

Previous attempts to determine the interactions between filariasis transmission intensity, infection and chronic disease have been limited by a lack of a theoretical framework that allows the explicit examination of mechanisms that may link these variables at the community level. Here, we show how deterministic mathematical models, in conjunction with analyses of standardized field data from communities with varying parasite transmission intensities, can provide a particularly powerful framework for investigating this topic. These models were based on adult worm population dynamics, worm initiated chronic disease and two major forms of acquired immunity (larval- versus adult-worm generated) explicitly linked to community transmission intensity as measured by the Annual Transmission Potential (ATP). They were then fitted to data from low, moderate and moderately high transmission communities from East Africa to determine the mechanistic relationships between transmission, infection and observed filarial morbidity. The results indicate a profound effect of transmission intensity on patent infection and chronic disease, and on the generation and impact of immunity on these variables. For infection, the analysis indicates that in areas of higher parasite transmission, community-specific microfilarial rates may increase proportionately with transmission intensity until moderated by the generation of herd immunity. This supports recent suggestions that acquired immunity in filariasis is transmission driven and may be significant only in areas of high transmission. In East Africa, this transmission threshold is likely to be higher than an ATP of at least 100. A new finding from the analysis of the disease data is that per capita worm pathogenicity could increase with transmission intensity such that the prevalences of both hydrocele and lymphoedema, even without immunopathological involvement, may increase disproportionately with transmission intensity. For lymphoedema, this rise may be further accelerated with the onset of immunopathology. An intriguing finding is that there may be at least two types of immunity operating in filariasis: one implicated in anti-infection immunity and generated by past experience of adult worms, the other involved in immune-mediated pathology and based on cumulative experience of infective larvae. If confirmed, these findings have important implications for the new global initiative to achieve control of this disease.

Epidemiology of soil-transmitted nematode infections in Ha Nam Province, Vietnam

Baseline epidemiological data are presented from a parasitological survey conducted in Thuy Loi commune, Ha Nam province, Vietnam; a farming community where night soil is routinely used as fertilizer for crops. 177 households were visited and 543 individuals (aged 1-88 years) recruited to the study. Helminth infection intensity was assessed by Kato-Katz to determine the density of parasite eggs per gram of stool (epg). Ascaris lumbricoides, Trichuris trichiura and hookworm infections were the only species detected. 83% of individuals were infected with A. lumbricoides (mean epg = 11,971), 94% with T. trichiura (mean epg = 793) and 59% with hookworm (mean epg = 302). Age-dependent patterns of infection prevalence and intensity were similar for A. lumbricoides and T. trichiura, but markedly different for hookworm infection. Similarly, age-dependency in the k-values for the three infections was due to covariance with the respective mean intensities with age rather than to independent age effects, with similar patterns for A. lumbricoides and T. trichiura, and a different pattern for hookworm. Three major conclusions can be drawn from the multiple-species analyses: There is positive interaction between A. lumbricoides and T. trichiura infections; high-intensity A. lumbricoides infections are significantly associated with high-intensity T. trichiura infections; and there is positive interaction between these two species such that infection intensity of A. lumbricoides is higher in individuals with concurrent T. trichiura infection than in individuals without and vice versa.