Email updates

Keep up to date with the latest news and content from Parasites & Vectors and BioMed Central.

This article is part of the series Mathematical models for parasites and vectors.

Open Access Research

Slaving and release in co-infection control

Laith Yakob1*, Gail M Williams1, Darren J Gray12, Kate Halton3, Juan Antonio Solon4 and Archie CA Clements1

Author Affiliations

1 Infectious Disease Epidemiology Unit, School of Population Health, University of Queensland, Herston, Brisbane, Qld 4006, Australia

2 Molecular Parasitology Laboratory, Queensland Institute of Medial Research, Herston, Brisbane, Qld 4029, Australia

3 Institute of Health and Biomedical Innovation, Queensland University of Technology, Kelvingrove, Qld 4059, Australia

4 College of Public Health, University of the Philippines Manila, 625 Pedro Gil Street, Ermita, Manila, Philippines

For all author emails, please log on.

Parasites & Vectors 2013, 6:157  doi:10.1186/1756-3305-6-157


The electronic version of this article is the complete one and can be found online at: http://www.parasitesandvectors.com/content/6/1/157


Received:1 March 2013
Accepted:23 May 2013
Published:31 May 2013

© 2013 Yakob et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Animal and human infection with multiple parasite species is the norm rather than the exception, and empirical studies and animal models have provided evidence for a diverse range of interactions among parasites. We demonstrate how an optimal control strategy should be tailored to the pathogen community and tempered by species-level knowledge of drug sensitivity with use of a simple epidemiological model of gastro-intestinal nematodes.

Methods

We construct a fully mechanistic model of macroparasite co-infection and use it to explore a range of control scenarios involving chemotherapy as well as improvements to sanitation.

Results

Scenarios are presented whereby control not only releases a more resistant parasite from antagonistic interactions, but risks increasing co-infection rates, exacerbating the burden of disease. In contrast, synergisms between species result in their becoming epidemiologically slaved within hosts, presenting a novel opportunity for controlling drug resistant parasites by targeting co-circulating species.

Conclusions

Understanding the effects on control of multi-parasite species interactions, and vice versa, is of increasing urgency in the advent of integrated mass intervention programmes.

Keywords:
Nematode; Mathematical model; Infectious disease; Epidemiology; Next generation matrix; Drug resistance

Background

The epidemiological and economic impact of multi-species co-infection is largely unknown and seldom considered in strategising control. However, many pathogens are typically found to co-circulate because hosts are a rare and patchy resource, and because of shared geographical distribution of parasites, and similar transmission routes [1]. The one-host one-pathogen paradigm that constitutes the backbone of infectious disease understanding requires expansion to account for this ubiquitous epidemiological setting. A framework is needed for optimising control when an infectious agent constitutes one of many interacting pathogens. Using gastro-intestinal nematode parasites as an illustrative example, we simulate different control strategies and their effects on infection prevalence when there are multiple, interacting parasites.

Gastro-intestinal nematodes are direct lifecycle parasites of veterinary and human health significance. Economically important nematode parasites of livestock include Haemonchus spp., Teladorsagia spp. and Nematodirus spp., which infect ruminants, Ascaris suum, a pig parasite, and Strongylus vulgaris and Strongyloides westeri in horses, among others. Human worm burden is mainly associated with Ascaris lumbricoides (large roundworm), Trichuris trichiura (whipworm) and the hookworms Ancylostoma duodenale and Necator americanus. These soil-transmitted helminths (STHs) infect over 1 billion people, incurring 40 million disability-adjusted life years, mostly in sub-Saharan Africa, East Asia, China, India and South America [2-4]. Gastro-intestinal nematode parasite transmission is through contact with contaminated soil or consumption of contaminated food or water, and co-infection with multiple species is widely-documented in humans and animals [5]. Co-infections have important consequences to parasite ecology [6] and associated host morbidity [1,7], and they warrant careful consideration when strategising control.

Control of gastro-intestinal nematodes is largely dependent on drugs of three main chemical classes: benzimidazoles, imidazothiazoles and macrocyclic lactones. Given the prevalence of co-infection and common approach to control for these parasites, there is increasing momentum in public health towards a model of integration, whereby multiple species are targeted simultaneously [8-13]. However, there are two key issues related to drug dependence that hold implications for this approach. First, chemotherapeutic efficacy is highly species specific. For example, cure rates of STHs following treatment with Albendazole are 88% for A. lumbricoides, 78% for the hookworms and only 28% for T. trichiura[14]. Furthermore, drug and multidrug resistance has become rife in livestock infections [15,16] and threatens the success of future mass drug application programmes in humans [17]. Second, the implications for control when host populations are co-infected with multiple, interacting parasite species are largely unknown. We developed a simple epidemiological model of co-infection to explore both these issues.

Given difficulties in obtaining empirical evidence, mathematical models have been developed to fill in data gaps and inform deworming programme strategies [18,19]. Most studies have followed in the footsteps of Anderson and May (1985) by analysing the population biology of a single helminth species [20]. Within-species interactions can be direct (e.g., via resource competition) or indirect (e.g., via immunomodulation) and generally are assumed to act upon two distinct life stages of the parasite: establishment of new infection or adult parasite fecundity [21]. The crux of our investigation, however, lies with the combinations of interactions that can arise from multi-species infections. By investigating the diverse range of interspecific interactions documented in laboratory animals [22-25] and humans [26,27], we discover epidemiological settings that can strongly enhance, attenuate, or even reverse, the intended impact of antihelminthic control programmes. In the advent of integrated parasite control we highlight the importance of developing a more holistic epidemiological framework for understanding the effects of perturbing parasite assemblages.

Methods

Rationale for a microparasite approach to modelling macroparasite control

Since Anderson and May (1978), there has been a divergence in the way by which ‘microparasites’ (e.g., viruses, bacteria, protozoa) and ‘macroparasites’ (e.g., nematodes, trematodes) have been modelled mathematically [28]. Typically, microparasite infection transmission is modelled using derivations of the classic SIR (Susceptible Infected Recovered) framework in which rates of change in proportions of the population between the epidemiological compartments are tracked using ordinary differential equations. This provides a constantly updated record of the prevalence of infection. Tracking prevalence of infection has been deemed insufficient for understanding macroparasite epidemiology because these types of parasites are highly aggregated among hosts, with hosts harbouring higher worm burdens (intensity) suffering more severe symptoms. Consequently, studies of macroparasite disease transmission, including co-infection models, have used a negative binomial forcing function to describe the aggregation of parasites within the host population [29-33].

While it would have been convenient to follow these examples, it would not have been appropriate for the intended purposes of our analysis. Key to our study is the interplay between parasite interactions and control. Although the negative binomial parasite aggregation method has proven very useful in describing endemic parasite distribution among hosts in the absence of control, it consists of a post hoc method of forcing the model output to better resemble the endemic equilibrium data. This method is not mechanistic - it does not attempt to describe the mechanism by which parasites become aggregated within host populations. Indeed, the precise mechanisms by which parasites aggregate within their hosts are unknown and likely to vary with species. Of paramount importance is the fact that the different mechanisms by which aggregation comes about will affect, and be affected by, control differently. Descriptive, phenomenological approaches (such as the negative binomial aggregation method) therefore become inappropriate for simulating non-equilibrial conditions, such as control scenarios.

One alternative modelling approach is to use individual-based simulations, such as those described by Bottomley et al. (2005) and Fenton et al. (2010) [34,35]. In the model of Bottomley et al. (2005), exposure to the parasite(s) is indistinguishable from host susceptibility to infection. The authors’ justification for this conflation was that, for their purposes, it did not matter whether a host remained uninfected because they were not exposed to a parasite or because the host was not susceptible to infection. While it has been useful in exploring interspecific effects of parasites at endemic equilibrium, this framework cannot be used for the simulation of control scenarios because the dynamics of free-living parasites and hosts’ infection status/susceptibility require explicit decoupling.

Fenton et al. (2010) also adopted an individual based modelling approach. These authors used a fixed lifelong probability of being infected with a particular parasite (with the probability assigned from a negative binomial distribution), in order to explore interspecific interactions during co-infection. This assumption was necessary because “allowing hosts to vary in their infection rates throughout their lifetime could introduce considerable noise to the data, greatly hampering the ability to detect any clear signal of interspecific interaction”. Unfortunately, this methodology cannot be applied to situations in which infection rates are anticipated to vary temporally, such as would be expected during, and following, control programmes. Therefore, to date, there is no modelling framework that can simulate macroparasite co-infection transmission dynamics in order to assess disease mitigation strategies. This was the motivation for our study.

The mathematical model of parasite co-infection

Two co-circulating gastro-intestinal nematode species were simulated using a simple adaptation of classic epidemiological model design whereby the host population was compartmentalised into Susceptible (S), Infected with nematode species 1 (I1) or 2 (I2) or both (I12). Infected hosts contaminate (β) the environment, and the subsequent free-living stages then either infect secondary hosts (λ) or perish (ν). The proportion of the total environment (E) that is contaminated with either (E1 or E2), or both (E12), nematode species was modelled explicitly. Traditionally, contamination of the environment is not modelled explicitly and, for mathematical convenience, exposure of hosts to free-living stages is assumed to be directly proportional to the adult worm population harboured by infected hosts [20]. We explicitly model the environment for three main reasons. First, it incorporates a more realistic time delay between the shedding of eggs and exposure of secondary hosts to infective stages. Second, it allows us to track separately the patterns of exposure and the subsequent infection status of the hosts (the importance of measuring both phenomena is exemplified in Figure 1). Third, explicitly modelling the environment bridges the gap between previous studies and future attempts to simulate more realistic scenarios that include spatial heterogeneities.

thumbnailFigure 1. Observations of co-infection prevalence can misrepresent between-parasite species interactions. In the above example, exposure to antagonistic parasites has resulted in no infection, or mono-specific infection, in most hosts. 5/10 people are infected with Parasite 1 (white) and 5/10 people are infected with Parasite 2 (black). More hosts are co-infected (4/10) than might be expected from random chance (5/10 x 5/10 = 2.5/10). This would typically lead to the incorrect assertion that the parasites exhibit synergistic interactions. Only by knowing both the patterns of infection status and exposure to contaminated environment can inferences be made of within-host between-parasite species interactions.

Infected individuals can return to the Susceptible class through recovery from infection (γ) under the assumption of no lasting immunity. Table 1 shows the symbology and definitions of the model parameters. Because many nematode infections can be long-lasting, host demographics were incorporated and a stable population was assumed (births = deaths). The set of ordinary differential equations used to describe two nematode species co-circulating within a host population and the environment is as follows:

<a onClick="popup('http://www.parasitesandvectors.com/content/6/1/157/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/6/1/157/mathml/M1">View MathML</a>

(1.1)

<a onClick="popup('http://www.parasitesandvectors.com/content/6/1/157/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/6/1/157/mathml/M2">View MathML</a>

(1.2)

<a onClick="popup('http://www.parasitesandvectors.com/content/6/1/157/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/6/1/157/mathml/M3">View MathML</a>

(1.3)

<a onClick="popup('http://www.parasitesandvectors.com/content/6/1/157/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/6/1/157/mathml/M4">View MathML</a>

(1.4)

<a onClick="popup('http://www.parasitesandvectors.com/content/6/1/157/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/6/1/157/mathml/M5">View MathML</a>

(1.5)

<a onClick="popup('http://www.parasitesandvectors.com/content/6/1/157/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/6/1/157/mathml/M6">View MathML</a>

(1.6)

Table 1. The parameters used in the model along with their definitions and simulated values

Where, S = 1 – (I1 + I2 + I12) under the assumption that hosts are born infection-free, and E = 1 – (E1 + E2 + E12). Our analysis explored the effect of synergistic and antagonistic interspecific interactions by adjusting the susceptibility to heterologous infection (respectively, α1 or α2 >1; α1 or α2 <1) and by adjusting the fecundity of heterologous parasites during co-infection (respectively, β1|12 or β2|12 >1; β1|12 or β2|12 <1).

A key characteristic of gastro-intestinal nematode epidemiology is the overdispersed distribution of parasites within a population: a minority of individual hosts harbour the majority of adult worms [36,37]. Reasons for this clumping are incompletely understood but are normally attributed to genetic predisposition, between-individual immunological differences and/or heterogeneities in behavioural risk factors [40]. Experimental work carried out with numerous parasite species has demonstrated strongly negative density-dependent effects, whereby worm survival and/or fecundity is compromised within hosts suffering higher burdens [41-43]. Other mathematical studies have described in great detail these single-species intra-specific effects on parasite population biology [20,21,44]. Higher per capita proliferation experienced at lower worm burdens has the inevitable effect of reducing the extent of heterogeneity in individual-host level parasite transmission potential that would otherwise be expected [45,46]. This generalisation facilitates our analysis by allowing us to focus on the prevalence (and not the heterogeneous intensities) of multiple parasite species that are co-circulating at the host-population level. Moreover, it provides a framework that is more analytically tractable and, therefore, more adaptable to a wider range of specific scenarios. An in-depth, individual-based, stochastic mathematical exploration of intensity-prevalence relationships in multi-helminth infections has been described elsewhere [34]. Further, in the context of an elimination programme, it will be the trickling effect of drug-attenuated low-intensity infections that perpetuates these parasites, as high-intensity infections will become vanishingly rare and substantially less influential on the population-level epidemiology.

Starting with the simplest case in which there is only one parasite species (Parasite ‘i’), the basic reproduction number is the product of the rate of environmental contamination and rate of host infection (accounting for host recovery as well as host and parasite mortality):

<a onClick="popup('http://www.parasitesandvectors.com/content/6/1/157/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/6/1/157/mathml/M7">View MathML</a>

(2)

Because we will be using next generation matrix methods for outbreak threshold calculations in the more complex stages of this analysis (Supporting Material), the above expression actually corresponds to the squared dominant eigenvalue of the next generation matrix. In order to calculate the invasion threshold of P2 following its introduction into a host population that already harbours P1 at endemic equilibrium, the rates of change for E1 and I1 are set to equal zero and the resulting simultaneous equations are then solved, leading to:

<a onClick="popup('http://www.parasitesandvectors.com/content/6/1/157/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/6/1/157/mathml/M8">View MathML</a>

(3.1)

<a onClick="popup('http://www.parasitesandvectors.com/content/6/1/157/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/6/1/157/mathml/M9">View MathML</a>

(3.2)

The hats over Î and Ê indicate that they are at their endemic equilibrium level. These can then be incorporated into the outbreak threshold calculation:

<a onClick="popup('http://www.parasitesandvectors.com/content/6/1/157/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/6/1/157/mathml/M10">View MathML</a>

(4)

In words, the introduced P2 either manifested in a susceptible host (1st term within square brackets), in a host already infected with P1 (2nd term) or a host that has since cleared P1 following P2 establishment (3rd term). This introduced case can then go on to infect a susceptible host (1st term in second set of square brackets) or a host already carrying P1 (2nd term). Initially, only the uni-directional effect of P1 on P2 was explored (α2 = β1|12 = 1); later scenarios allowed for bi-directionality in parasite interactions (Supporting Material).

The design of the model means that it is also applicable to the exploration of within-parasite species i.e. different strains of same species, instead of different species, interactions. Asymmetric competitive interactions between strains can be used to explore evolutionary selectivity with differences in susceptibility to chemotherapy representing the development of drug resistance. In this circumstance, calculations of the successful spread of a secondary parasite strain into a host population already endemic for a primary strain (R2|1) demonstrate the parasitological life history trade-offs that need to be negotiated in order for drug resistance to spread.

Results

Synergisms and antagonisms in animal and human co-infection are typically inferred through (co-)infection prevalence data. Figure 1 illustrates how this can misrepresent the true nature of between-pathogen interactions, which can only be achieved with knowledge of exposure levels, as well as host infection status. We demonstrate the profound public health implications posed by this popular misconception with a simple mathematical model that simulates the attempted intervention of two co-circulating parasite species.

Unidirectional effects of Parasite species 1 on Parasite species 2

First, we consider a co-infection transmission model that only accounts for one direction of interactions (the effects of parasite species 1 ‘P1’ on parasite species 2 ‘P2’). Through host immunological modulation and/or direct resource competition, an extant parasite can enhance or suppress both the host susceptibility to, and the fecundity of, heterologous parasites. Figure 2 shows the four combinations of antagonism/synergism for P2 establishment/fecundity resulting from interaction with P1 (the effect on outbreak thresholds is shown in Additional file 1: Figure S1). The prevalence of P1, P2 and co-infections are plotted for a range of P1-specific control levels. This simulates a species-specific infection treatment or, alternatively, the administering of a drug to which P2 is resistant. When P1 facilitates the establishment of P2 infection as well as its subsequent fecundity, the pathogens become epidemiologically slaved whereby the disproportionate majority of infections are tethered in co-infections (Figure 2D). ‘Allee effects’ are then experienced by P2 whereby its prevalence diminishes with increasing levels of P1 control. In other words, inadvertent suppression of a secondary pathogen occurs by the administration of a drug that would be assumed to have no effect in a one-host one-pathogen mindset.

thumbnailFigure 2. Parasite prevalence following chemotherapeutic control of P1. P2 and co-infection are the solid and broken lines on the left y-axis, P1 is the thick black line on the right y-axis. The between-parasite species interactions for these four scenarios are: A) P1 reduces both host susceptibility to, and fecundity of, P2 (−/−); B) P1 reduces host susceptibility to, but increases fecundity of, P2 (−/+); C) P1 increases host susceptibility to, but reduces fecundity of, P2 (+/−); D) P1 increases host susceptibility to, and increases fecundity of, P2 (+/+).

Additional file 1. Yakob et al. Slaving and Release in Co-infection Control. Allowing for bi-directionality in the inter-specific interactions of P1 and P2.

Format: DOC Size: 725KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

The effect of intervention on (co-)infection prevalence is very different when mixed synergistic/antagonistic pressures are exerted by P1 (Figure 2B,C). Here, moderate levels of P1 control exacerbate the prevalence of P2 as well as co-infections. This is of major concern given the enhanced morbidity often associated with co-infection [7]. There is a trade-off between P2 being released from suppression by P1 at one stage of its lifecycle and its transmission being enhanced at a different stage. Therefore, when elimination is unachievable, infectious disease suppression is not necessarily the next-best option; the pros of suppressing the more prevalent parasite must be weighed against the cons of potentially exacerbating co-infections and associated co-morbidities.

Bidirectional effects between Parasite species 1 and Parasite species 2

Until now, we have only considered uni-directional pathogen interactions (the effect of P1 on P2). However, incorporating bi-directionality is important because, typically, nematode co-infections are neither uni-directional, mutually synergistic nor mutually antagonistic, but rather, a more complex blend [47-49]. As a result of this additional flexibility, the number of qualitatively distinct combinations increases, resulting in a richer diversity in epidemiological outcomes. Figure 3 shows the prevalence of these pathogens following treatment of the host population for P1 infection. When the parasites are mutually beneficial (Figure 3D iv), drug control of P1 is attenuated, necessitating aggressive levels of P1 control before any drop in overall prevalence is actually achieved. In accordance with the simpler, uni-directional simulations, there are non-monotonic effects between control effort and pathogen transmission for mixed synergistic/antagonistic combinations (e.g., Figure 3B iii), and moderate levels of control risk increasing the overall prevalence of co-infections. The qualitatively different responses of co-circulating parasite epidemiology following moderate P1 control are summarised in Figure 4. In general, three different adverse outcomes can result from attempting P1 control in a co-infection system: 1) The presence of co-circulating P2 can buffer the effect of chemotherapeutics, providing an indirect resistance of P1 to control (when the parasites are mutually synergistic e.g., top scenario of Figure 4); 2) Controlling P1 can result in the competitive release of P2 whose transmission would otherwise be suppressed by P1 (when P1 reduces host susceptibility to heterologous infection but enhances heterologous parasite fecundity e.g., top scenario of synergistic/antagonistic mix of Figure 4); and, 3) Co-infection exacerbation can occur (when the previous situation also involves a positive effect of P2 infection on the host susceptibility to P1). This last combination is the worst case scenario as it yields all three adverse effects.

thumbnailFigure 3. Parasite prevalence under control of Parasite 1 and allowing for bi-directionality in their effect. P2 and co-infection are the solid and broken lines on the left y-axis, P1 is the thick black line on the right y-axis; the x-axis is the fold reduction in Parasite 1 infection duration. The effect of Parasite 1 on Parasite 2 is depicted by the four larger quadrants: A) P1 reduces both host susceptibility to, and fecundity of, P2 (−/−); B) P1 reduces host susceptibility to, but increases fecundity of, P2 (−/+); C) P1 increases host susceptibility to, but reduces fecundity of, P2 (+/−); D) P1 increases host susceptibility to, and increases fecundity of, P2 (+/+). The equivalent effect of Parasite 2 on Parasite 1 is represented in the smaller sub-plots (respectively i, ii, iii and iv).

thumbnailFigure 4. The qualitative effects of moderate P1 control are dependent on interspecific interactions between P1 and P2. “Moderate control” is 50% of the level of chemotherapeutic control required for elimination. Interactions can be synergistic (+) or antagonistic (˗) and can occur at different stages of the transmission cycle (the left interaction symbol denotes host susceptibility to heterologous infection, and, the right symbol is fecundity of heterologous adult parasites following infection establishment). Where the polarity of these interactions does not affect qualitative behaviour, the interaction is left blank. The prevalence of P1 and P2 before and after moderate P1 control is represented by the circle size, and the prevalence of co-infection by the amount of intersection. Where P2 fails to persist, it is represented by a point. To exemplify the link between this conceptual figure and the corresponding simulation results, the purely synergistic scenario (top interaction) corresponds to half way across the x-axis of Figure 3Div, and the purely antagonistic scenario (second from top interaction) corresponds to half way across the x-axis of Figure 3Ai.

Chemotherapeutic control using drugs that are equivalently effective on both pathogens, and reducing environmental contamination with pre-adult parasite stages by improving sanitation, yield very similar results to each other (Additional file 1: Figure S2 & S3). The range of epidemiological patterns is constrained when controls have symmetrical effects on each of the pathogens. In general, presence of a secondary parasite only affects control of the dominant parasite when it enhances its transmission. Drugs and sanitation then more effectively control P2 (the pathogen with a lower intrinsic rate of transmission) and co-infections, while suppression of P1 is not achieved until the heterologous pathogen is eliminated.

Discussion

The results we present encompass a broad continuum of interspecific interactions. Determining the extent to which co-infecting parasite species interact in nature is not trivial. Using log-linear regression, Howard et al. (2001) analysed metadata collected on 215000 individuals in Africa and Asia [50]. Although most significant associations between the four main STH species were positive, significant negative associations were also found for each of the different combinations of co-infection. Poulin (2005) and Fenton et al. (2010) recently described the inherent difficulties involved in tweezing apart macroparasite interactions from ecological data [35,51]. These authors suggest that current methods that are routinely implemented for detecting interactions are unreliable and likely underestimate the extent to which they occur in reality. Moreover, ecological data of association cannot distinguish between the different mechanisms of interaction occurring at different stages of the transmission cycle – a feature that we and a previous, single species modelling study demonstrate to be critical in determining the impact of control [21].

Much of our understanding, therefore, comes from experimental animal models. Cox (2001) has compiled a comprehensive list of different types of mixed-species parasite interactions mediated by the immune systems of animal models [5]. Examples include Trichuris muris and Heligmosomoides polygyrus in which the trichurid benefits during mouse co-infection, and Schistosoma mansoni (a trematode) and Strongyloides venezuelensis (a nematode) resulting in a detrimental effect for the nematode [52,53]. However, it is difficult to know precisely how representative these animal models are of real-world epidemiology when qualitatively inconsistent findings are obtained from highly related mammalian host species [54,55].

Efforts to ameliorate morbidity and mortality from infectious disease require a shift in the way in which pathogenic transmission is conceptualised. The one-host one-parasite epidemiological paradigm does not apply in a world where infections are largely multi-host and concomitant [1,56]. The mathematical extensions to epidemiological models used in this study are very simplistic; it was our intention to demonstrate transparently the unintuitive, and sometimes counter-intuitive, repercussions of perturbing a pathogen when it constitutes merely a single component of a more complex community. The complex spectrum of qualitative behaviours presented here will only become increasingly multifarious when host heterogeneities and age-structure are included, or a third competing parasite species is considered [34,57]. As we continue to build upon our model, however, we do not anticipate revising the conclusions of our study.

Conclusions

Data on co-infection prevalence is insufficient to infer upon between-species parasite interactions. Moreover, population-level intervention efficacy can be seriously compromised if between-parasite species interactions are not taken into account. We have shown how knowledge of parasite interactions can be used to facilitate infection control, even in the presence of drug resistance. Although this demonstration of indirect control through targeting a co-circulating synergist offers a novel tool for combating drug insensitivity, of greater immediacy is elucidating when and where public health might actually be compromised by injudicious control. In following the construction of classic epidemiological models, our framework offers an easy transition between the current paradigm and future studies that consider explicitly the pathogen species assemblage.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors developed the concept of the study. LY produced and analysed the mathematical model and wrote the first draft of the manuscript. All authors read and approved the final version of the manuscript.

Acknowledgements

Funding is provided by the National Health and Medical Research Council project grant 1006254. Darren Gray is an Australian Research Council (DECRA) Fellow. Archie Clements is supported by a NHMRC Career Development Fellowship.

References

  1. Petney TN, Andrews RH: Multiparasite communities in animals and humans: frequency, structure and pathogenic significance.

    Int J Parasitol 1998, 28:377-393. OpenURL

  2. Chan MS: The global burden of intestinal nematodes infections – fifty years on.

    Parasitol Today 1997, 13:438-443. OpenURL

  3. Brooker S: Estimating the global distribution and disease burden of intestinal nematode infections: Adding up the numbers – A review.

    Int J Parasitol 2010, 40:1137-1144. OpenURL

  4. Pullan RBS: The global limits and population at risk of soil-transmitted helminth infections in 2010.

    Parasit Vectors 2012, 5(1):81. OpenURL

  5. Cox FEG: Concomitant infections, parasites and immune responses.

    Parasitology 2001, 122(Suppl):S23-S38. OpenURL

  6. Graham AL: Ecological rules governing helminth–micro parasite co infection.

    Proc Natl Acad Sci 2008, 105:566-570. OpenURL

  7. Pullan R, Brooker S: The health impact of polyparasitism in humans: are we under-estimating the burden of parasitic diseases?

    Parasitology 2008, 135:783-794. OpenURL

  8. Brooker S, Kabatereine NB, Fleming F, Devlin N: Cost and cost–effectiveness of nationwide school-based helminth control in Uganda: intra-country variation and effects of scaling-up.

    Health Policy Plan 2008, 23:24-35. OpenURL

  9. Mathers CD, Ezzatti M, Lopez AD: Measuring the burden of neglected tropical diseases: the global burden of disease framework.

    PLoS Negl Trop Dis 2007, 1:E114. OpenURL

  10. Brooker S, Clements ACA: Spatial heterogeneity of parasite co-infection: determinants and geostatistical prediction at regional scales.

    Int J Parasitol 2009, 39:591-597. OpenURL

  11. Clements ACA, Deville MA, Ndayishimiye O, Fenwick A: Spatial co-distribution of neglected tropical diseases in the East African Great Lakes Region: Revisiting the justification for integrated control.

    Trop Med Int Health 2010, 15:198-207. OpenURL

  12. Utzinger J, De Savigny D: Control of neglected tropical diseases: integrated chemotherapy and beyond.

    PLoS Med 2006, 3:e112. OpenURL

  13. Mwinzi PMS, Owaga C, Mwanje M, Muok E, Ayisi J, Laserson K, Muchiri E, Secor W, Karanja D: Integrated community-directed intervention for schistosomiasis and soil transmitted helminths in western Kenya - a pilot study.

    Parasit Vectors 2012, 5(1):182. OpenURL

  14. Keiser J, Utzinger J: Efficacy of current drugs against soil-transmitted helminth infections: systematic review and meta-analysis.

    J Am Med Assoc 2008, 299:1937-1948. OpenURL

  15. Besier B: New anthelmintics for livestock: the time is right.

    Trends Parasitol 2007, 23:21-24. OpenURL

  16. Waller PJ: Anthelmintic resistance.

    Vet Parasitol 1997, 72:391-412. OpenURL

  17. Albonico M, Engels D, Savioli L: Monitoring drug efficacy and early detection of drug resistance in human soil-transmitted nematodes: a pressing public health agenda for helminth control.

    Int J Parasitol 2004, 34:1205-1210. OpenURL

  18. Medley GF, Guyatt HL, Bundy DAP: A quantitative framework for evaluating the effect of community treatment on the morbidity due to ascariasis.

    Parasitology 1993, 106:201-221. OpenURL

  19. Smith G: A mathematical model for the evolution of anthelmintic resistance in a direct life cycle nematode parasite.

    Int J Parasitol 1990, 20:913-921. OpenURL

  20. Anderson RM, May RM: Helminth Infections of Humans: Mathematical Models, Population Dynamics, and Control.

    Adv Parasitol 1985, 24:1-101. OpenURL

  21. Churcher TS, Filipe JAN, Basáñez M-G: Density dependence and the control of helminth parasites.

    J Anim Ecol 2006, 75:1313-1320. OpenURL

  22. Behnke JM, Wahid FN, Grencis RK, Else KJ, Ben-Smith A, Goyal PK: Immunological relationships during primary infection with Heligmosomoides polygyrus (Nematospiroides dubius): down regulation of specific cytokine secretion (IL-9 and IL-10) correlates with poor mastocytosis and chronic survival of adult worms.

    Parasite Immunol 1993, 15:415-421. OpenURL

  23. Behnke JM: Structure in parasite component communities in wild rodents: predictability, stability, associations and interactions or pure randomness?

    Parasitology 2008, 135:751-766. OpenURL

  24. Hermanek J, Goyal PK, Wakelin D: Lymphocyte, antibody and cytokine responses during concurrent infections between helminths that selectively promote T-helper-1 or T-helper-2 activity.

    Parasite Immunol 1994, 16:111-117. OpenURL

  25. Christensen NO, Nansen P, Fagbemi BO, Monrad J: Heterologous antagonistic and synergistic interactions between helminths and between helminths and protozoans in concurrent experimental-infection of mammalian hosts.

    Parasitol Res 1987, 73:387-410. OpenURL

  26. Fleming FM, Brooker S, Geiger SM, Caldas IR, Correa-Oliveira R, Hotez PJ, Bethony JM: Synergistic associations between hookworm and other helminth species in a rural community in Brazil.

    Trop Med Int Health 2006, 11:56-64. OpenURL

  27. Webster M, Correa-Oliveira R, Gazzinelli G, Viana IRC, Fraga LADO, Silveira AMS, Dunne DW: Factors affecting high and low human IgE responses to schistosome worm antigens in an area of Brazil endemic for Schistosoma mansoni and hookworm.

    AmJTrop Med Hyg 1997, 57(4):487-494. OpenURL

  28. Anderson RM, May RM: Regulation and stability of host-parasite population interactions: I. Regulatory processes.

    J Anim Ecol 1978, 47:219-247. OpenURL

  29. Dobson AP: The population-dynamics of competition between parasites.

    Parasitology 1985, 91:317-347. OpenURL

  30. Roberts MG, Dobson AP: The population-dynamics of communities of parasitic helminths.

    Math Biosci 1995, 126:191-214. OpenURL

  31. Gatto M, DeLeo GA: Interspecific competition among macro parasites in a density-dependent host population.

    J Math Biol 1998, 37:467-490. OpenURL

  32. Pugliese A: Coexistence of macro parasites without direct interactions.

    Theor Popul Biol 2000, 57:145-165. OpenURL

  33. Fenton A: Worms and germs: the population dynamic consequences of microparasite-macroparasite co-infection.

    Parasitology 2008, 135:1545-1560. OpenURL

  34. Bottomley C, Isham VS, Basáñez M-G: Population biology of multispecies helminth infection: interspecific interactions and parasite distribution.

    Parasitology 2005, 131:417-433. OpenURL

  35. Fenton A, Viney ME, Lello J: Detecting interspecific macro parasite interactions from ecological data: patterns and process.

    Ecol Lett 2010, 13:606-615. OpenURL

  36. Anderson RM, May RM: Infectious Diseases of Humans: Dynamics and Control. Oxford and New York: Oxford University Press; 1991:467-529. OpenURL

  37. Hoagland KE, Schad GA: Necator americanus and Ancylostoma duodenale: Life history parameters and epidemiological implications of two sympatric hookworms of humans.

    Exp Parasitol 1978, 44:36-49. OpenURL

  38. Udonsi JK, Atata G: Necator americanus: Temperature, pH, light, and larval development, longevity, and desiccation tolerance.

    Exp Parasitol 1987, 63(2):136-142. OpenURL

  39. Woolhouse MEJ, Dye C, Etard J-F, Smith T, Charlwood JD, Garnett GP, Hagan P, Hii JLK, Ndhlovu PD, Quinnell RJ, et al.: Heterogeneities in the transmission of infectious agents: Implications for the design of control programs.

    Proc Natl Acad Sci 1997, 94:338-342. OpenURL

  40. Hudson PJ, Dobson AP: Macroparasites: Observed Patterns. In Ecology of Infectious Diseases in Natural Populations. Edited by Grenfell BT, Dobson AP. Cambridge: Cambridge University Press; 1995:144-176. OpenURL

  41. Keymer AE: Density-dependent mechanisms in the regulation of intestinal helminth populations.

    Parasitology 1982, 84:573-587. OpenURL

  42. Chylinski C, Boag B, Stear MJ, Cattadori IM: Effects of host characteristics and parasite intensity on growth and fecundity of Trichostrongylus retortaeformis infections in rabbits.

    Parasitology 2009, 136(01):117-123. OpenURL

  43. Paterson S, Viney ME: Host immune responses are necessary for density dependence in nematode infections.

    Parasitology 2002, 125:283-292. OpenURL

  44. Grenfell BT, Wilson K, Isham VS, Boyd HEG, Dietz K: Modelling patterns of parasite aggregation in natural populations: trichostrongylid nematode-ruminant interactions as a case study.

    Parasitology 1995, 111(Suppl):S135-S151. OpenURL

  45. Quinnell RJ, Keymer AE, Medley GF: The regulation of gastrointestinal helminth populations.

    Phil Trans Roy Soc B 1990, 330:191-201. OpenURL

  46. Galvani AP: Immunity, antigenic heterogeneity, and aggregation of helminth parasites.

    J Parasitol 2003, 89:232-241. OpenURL

  47. Dash KM: Interaction between Oesophagostomum columbianum and Oesophagostomum venulosum in sheep.

    Int J Parasitol 1981, 11:201-207. OpenURL

  48. Holland C: Interactions between Moniliformis (Acanthocephala) and Nippostrongylus (Nematoda) in the small intestine of laboratory rats.

    Parasitology 1984, 88:303-315. OpenURL

  49. Poulin R: Interactions between species and the structure of helminth communities.

    Parasitology 2001, 122:S3-S11. OpenURL

  50. Howard SC, Donnelly CA, Chan MS: Methods for estimation of associations between multiple species parasite infections.

    Parasitology 2001, 122:233-251. OpenURL

  51. Poulin R: Detection of interspecific competition in parasite communities.

    J Parasitol 2005, 91(5):1232-1235. OpenURL

  52. Behnke JM, Ali NMH, Jenkins SN: Survival to patency of low-level infections with Trichuris muris in mice concurrently infected with Nematospiroides dubius.

    Annals Trop Med Hyg 1984, 78:509-517. OpenURL

  53. Yoshida A, Maruyama H, Yabu Y, Amano T, Kobayakawa T, Ohta N: Immune responses against protozoal and nematodal infection in mice with underlying Schistosoma mansoni infection.

    Parasitol Int 1999, 48:73-79. OpenURL

  54. Holmes JC: Effects of concurrent infections on Hymenolepis diminta (Cestoda) and Moniliformis dubius (Acanthocephala). II. Effects on growth.

    J Parasitol 1962, 48:87-96. OpenURL

  55. Holmes JC: Effects of concurrent infections on Hymenolepis diminta (Cestoda) and Moniliformis dubius (Acanthocephala). II. Effects in hamsters.

    J Parasitol 1962, 48:97-100. OpenURL

  56. Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL, Daszak P: Global trends in emerging infectious diseases.

    Nature 2008, 451:990-993. OpenURL

  57. Basáñez M-G, French MD, Walker M, Churcher TS: Paradigm lost: how parasite control may alter pattern and process in human helminthiases.

    Trends Parasitol 2012, 28:161-171. OpenURL