- Open Access
Slaving and release in co-infection control
Parasites & Vectorsvolume 6, Article number: 157 (2013)
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.
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.
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.
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.
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 . 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 . Co-infections have important consequences to parasite ecology  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. 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 . 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 . 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 . 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.
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 . 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 . 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.
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:
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 . 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 . 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):
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:
The hats over Î and Ê indicate that they are at their endemic equilibrium level. These can then be incorporated into the outbreak threshold calculation:
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.
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.
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 . 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.
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.
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 . 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 .
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 . 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.
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.
Petney TN, Andrews RH: Multiparasite communities in animals and humans: frequency, structure and pathogenic significance. Int J Parasitol. 1998, 28: 377-393.
Chan MS: The global burden of intestinal nematodes infections – fifty years on. Parasitol Today. 1997, 13: 438-443.
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.
Pullan RBS: The global limits and population at risk of soil-transmitted helminth infections in 2010. Parasit Vectors. 2012, 5 (1): 81-
Cox FEG: Concomitant infections, parasites and immune responses. Parasitology. 2001, 122 (Suppl): S23-S38.
Graham AL: Ecological rules governing helminth–micro parasite co infection. Proc Natl Acad Sci. 2008, 105: 566-570.
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.
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.
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-
Brooker S, Clements ACA: Spatial heterogeneity of parasite co-infection: determinants and geostatistical prediction at regional scales. Int J Parasitol. 2009, 39: 591-597.
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.
Utzinger J, De Savigny D: Control of neglected tropical diseases: integrated chemotherapy and beyond. PLoS Med. 2006, 3: e112-
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-
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.
Besier B: New anthelmintics for livestock: the time is right. Trends Parasitol. 2007, 23: 21-24.
Waller PJ: Anthelmintic resistance. Vet Parasitol. 1997, 72: 391-412.
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.
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.
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.
Anderson RM, May RM: Helminth Infections of Humans: Mathematical Models, Population Dynamics, and Control. Adv Parasitol. 1985, 24: 1-101.
Churcher TS, Filipe JAN, Basáñez M-G: Density dependence and the control of helminth parasites. J Anim Ecol. 2006, 75: 1313-1320.
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.
Behnke JM: Structure in parasite component communities in wild rodents: predictability, stability, associations and interactions or pure randomness?. Parasitology. 2008, 135: 751-766.
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.
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.
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.
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.
Anderson RM, May RM: Regulation and stability of host-parasite population interactions: I. Regulatory processes. J Anim Ecol. 1978, 47: 219-247.
Dobson AP: The population-dynamics of competition between parasites. Parasitology. 1985, 91: 317-347.
Roberts MG, Dobson AP: The population-dynamics of communities of parasitic helminths. Math Biosci. 1995, 126: 191-214.
Gatto M, DeLeo GA: Interspecific competition among macro parasites in a density-dependent host population. J Math Biol. 1998, 37: 467-490.
Pugliese A: Coexistence of macro parasites without direct interactions. Theor Popul Biol. 2000, 57: 145-165.
Fenton A: Worms and germs: the population dynamic consequences of microparasite-macroparasite co-infection. Parasitology. 2008, 135: 1545-1560.
Bottomley C, Isham VS, Basáñez M-G: Population biology of multispecies helminth infection: interspecific interactions and parasite distribution. Parasitology. 2005, 131: 417-433.
Fenton A, Viney ME, Lello J: Detecting interspecific macro parasite interactions from ecological data: patterns and process. Ecol Lett. 2010, 13: 606-615.
Anderson RM, May RM: Infectious Diseases of Humans: Dynamics and Control. 1991, Oxford and New York: Oxford University Press, 467-529.
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.
Udonsi JK, Atata G: Necator americanus: Temperature, pH, light, and larval development, longevity, and desiccation tolerance. Exp Parasitol. 1987, 63 (2): 136-142.
Woolhouse MEJ, Dye C, Etard J-F, Smith T, Charlwood JD, Garnett GP, Hagan P, Hii JLK, Ndhlovu PD, Quinnell RJ: Heterogeneities in the transmission of infectious agents: Implications for the design of control programs. Proc Natl Acad Sci. 1997, 94: 338-342.
Hudson PJ, Dobson AP: Macroparasites: Observed Patterns. Ecology of Infectious Diseases in Natural Populations. Edited by: Grenfell BT, Dobson AP. 1995, Cambridge: Cambridge University Press, 144-176.
Keymer AE: Density-dependent mechanisms in the regulation of intestinal helminth populations. Parasitology. 1982, 84: 573-587.
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.
Paterson S, Viney ME: Host immune responses are necessary for density dependence in nematode infections. Parasitology. 2002, 125: 283-292.
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.
Quinnell RJ, Keymer AE, Medley GF: The regulation of gastrointestinal helminth populations. Phil Trans Roy Soc B. 1990, 330: 191-201.
Galvani AP: Immunity, antigenic heterogeneity, and aggregation of helminth parasites. J Parasitol. 2003, 89: 232-241.
Dash KM: Interaction between Oesophagostomum columbianum and Oesophagostomum venulosum in sheep. Int J Parasitol. 1981, 11: 201-207.
Holland C: Interactions between Moniliformis (Acanthocephala) and Nippostrongylus (Nematoda) in the small intestine of laboratory rats. Parasitology. 1984, 88: 303-315.
Poulin R: Interactions between species and the structure of helminth communities. Parasitology. 2001, 122: S3-S11.
Howard SC, Donnelly CA, Chan MS: Methods for estimation of associations between multiple species parasite infections. Parasitology. 2001, 122: 233-251.
Poulin R: Detection of interspecific competition in parasite communities. J Parasitol. 2005, 91 (5): 1232-1235.
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.
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.
Holmes JC: Effects of concurrent infections on Hymenolepis diminta (Cestoda) and Moniliformis dubius (Acanthocephala). II. Effects on growth. J Parasitol. 1962, 48: 87-96.
Holmes JC: Effects of concurrent infections on Hymenolepis diminta (Cestoda) and Moniliformis dubius (Acanthocephala). II. Effects in hamsters. J Parasitol. 1962, 48: 97-100.
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.
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.
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.
The authors declare that they have no competing interests.
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.