ABSTRACT
Coinfections pose serious threats to health and exacerbate parasite burden. If coinfection is detrimental, then what within-host factors facilitate it? Equally importantly, what hinders it, say via exclusion or priority effects? Such interactions ought to stem from their within-host environment (‘niche’), i.e., resources that parasites steal from hosts and immune cells that kill them. Yet, despite two decades of empirical focus on within-host infection dynamics, we lack a mechanistic framework to understand why coinfection arises and the diverse range of its’ consequences. Hence, we construct a trait-based niche framework, one that illustrates general principles that govern parasite competition for a resource and apparent competition for immune cells. We show that coinfection requires a competition-resistance tradeoff and that each parasite most impacts the niche factor to which it is most sensitive. These predictions then provide mechanistic interpretation for infection outcomes seen in a variety of extant experiments: Why does nutrient supplementation shift relative frequencies of coinfecting parasites? When and how does sequence of parasite invasion allow only early invading parasites to win? How does intrinsic variation in immune response shape coinfection burden? Together, this mechanistic framework of parasite competition offers new perspectives to better predict within-host infection dynamics and improve individual health.
INTRODUCTION
Coinfections pose a serious threat to health both at the individual and population scale for hosts. Broadly, coinfection is the successful concurrent infestation of a host with two or more parasite species. Coinfection worsens human health (76% of 2000 publications in a meta-analysis) and exacerbates infection burden (57%; Griffiths et al 2011). For instance, bacterial co-infections cause two-fold increase in mortality in COVID patients (Shah et al 2023). Similarly, coinfecting hghhhelminth species increased the odds of anaemia five-to-eight-fold among children (Ezeamama et al 2005). Coinfection at the individual scale can also influence population-scale disease dynamics (Ezenwa & Jolles 2011). Populations having higher frequency of coinfection experience larger epidemics than those with epidemics of one parasite (Susi et al 2015). Yet, despite virulent costs, coinfections are pervasive amongst hosts, including humans, wildlife, and livestoc/ agriculture (Vogels et al 2018; Ezenwa et al 2010; Halliday, Penczykowski et al 2020). Therefore, we need to better understand why and how parasites coinfect their hosts.
Nonetheless, coinfection represents just one of numerous outcomes of concurrent pathogenesis. That fact poses fundamental questions about coinfection. First, if coinfection is indeed detrimental, then what within-host factors facilitate it? Equally importantly, what prevents it? Within-host parasite competition can lead to coinfection (i.e., within-host parasite coexistence), singly infection (through exclusion or priority effects), or clearance (no infection) from hosts (Vogel et al 2018). Yet, we lack theoretical predictions to elucidate those outcomes of competition. For instance, irrespective of timing, sequential exposure to some parasites always lead to coinfection, yet to others, only early infecting parasites prevail (Fig. 1A; Clay et al 2019; Devevey et al 2015). Second, when they coinfect why do some parasites become more abundant than others in relative and absolute senses? For instance, nutrient supply or resources can shift this ‘community structure’ of parasites within hosts favouring one parasite over the other (Fig. 1B; Fellous & Koella 2009, Budischak et al 2015). Finally, not all coinfected hosts are equal. Some hosts present higher coinfection burden than others, hinting at variation in immunological resistance (Fig. 1C; Fuess et al 2021, Halliday et al 2018). How does such variation govern coinfection burden? Presently, these disparate bits about coinfection remain unconnected, lacking a synthetic glue. Here, we seek to catalyse creation of a framework linking genesis of coinfection to its consequences (Lively et al 2014; Restif & Graham 2015).
A mechanistic, within-host framework of parasite competition provides a start. To illustrate, we develop and evaluate within-host models of coinfection, synthesizing insights through an ecological lens (following Ramesh & Hall 2023). Broadly, coinfection (i.e., within-host coexistence of parasites) occurs when each parasite species competes more strongly with themselves (intraspecifically) than with each other (interspecifically; Fig. 1D). Those competitive differences could be measured phenomenologically. Even better, they can be quantified from trait-based niche models for parasites competing within hosts (Graham et al 2008, Cressler et al 2014). In such models, parasites interact within host ‘ecosystems’ for shared resources while facing attack by energetically costly immune cells. Simple ecological rules, then, govern within-host parasite competition. To illustrate, we first borrow from old and new niche models to outline how shared energy (resources) and/or immune cells could govern divergent coinfection outcomes (Fig. 2). Second, using a case study of a two parasite – immune cells – energy model (2PIE) niche model, we link outcomes of infection to competitive abilities, nullclines (niches), feedback, and key traits (Figs. 3-7; Table S1,2). Third, our approach scripts how to link a priori predictions to interpretation of experimental outcomes using these types of models (Fig. 8-9). Finally, we offer suggestions for future theory and experiments (Fig. 10, Table 1). Taken together, we lay a general mechanistic framework to understand coinfection using within-host niche-based competition model.
Section I: A phenomenological approach to coinfection
To visualize these interactions, we label the direct effect of species j on growth rate of species i, yielding (hereafter) interspecific positive effects (black arrow) and negative effects (red arrow) and intraspecific self-limitation (red curve; Fig. 1,2). Then, we can write the two-species Lotka–Volterra competition model in terms of intra (aii)- and inter (aij)-specific per capita competition coefficients: where Pi is the density of parasite species i, and ri is its intrinsic (maximum per capita) growth rate. This model assumes that per capita growth rates (fitness) decrease linearly with density of each species with slopes aii and aij (Table 1A). Successful coinfection satisfies two conditions. First, each single-species (‘boundary’) carrying capacity (1/aii) must be invasible by the other parasite when rare. That condition is met when each species reduces its own growth rate more than that of its competitor (i.e., a11 > a21 and a22 > a12). Second, joint intraspecific (a11 a22) must exceed joint interspecific competition (a12 a21). When both are met, the feasible coexistence (‘interior’) equilibrium is stable with negative feedback (Table 1A).
Coinfection is hindered in one of two ways. First, strong positive feedback at the interior generates priority effects when joint inter-exceeds intra-specific competition (a12 a21 > a11 a22). Priority effects also require mutually uninvasible boundary equilibria and a feasible interior (which are both met when a11 < a21 and a22 < a12). Here, high enough initial density of parasites determines the winner. Second, competitive exclusion arises with asymmetric competition. Specifically, if P2 cannot invade P1’s boundary (a11 < a21) but P1 can invade P2’s boundary (a22 > a12), then P1 excludes P2. Conversely, P2 excludes P1 when a11 > a21 but a22 < a12. Hence, mutually or asymmetrically strong interspecific competition can hinder coinfection.
Empirically, estimates of competition coefficients can come from densities and growth rates of parasites. Various methods exist, including a response surface design. In the context of macroparasites, this design involves growing species either alone or at the lowest densities that lead to infection. It also involves growing species with increasing density of intra- and inter-specific competitors. Through this design, aii and aij can be estimated, then used to predict if species can successfully infect when rare (Freckleton & Watkinson 2000; reviewed in Hart et al 2018; see Table 1A). However, this approach remains phenomenological, measuring coefficients only in a specific environment. That context-dependence limits predictions of infection outcomes across gradients of traits, nutrient supply, or in allocation to immunity. Instead, predictions in these other contexts require a mechanistic framework of competition between parasites involving host immunity and/or host resources.
Section II: Within-host feedbacks driving divergent infection outcomes
To move toward creation of such a framework, we outline six models of within-host competition between parasites involving immunity, resources, and/or direct interference (Fig. 2). All potentially contain mechanisms for coinfection, single infection due to exclusion or priority effects (aka, alternative states or founder control), and no infection. Although details can vary, such divergent outcomes in these models follow general principles of species competition as outlined above. In essence, coinfection occurs when strong intraspecific competition (mediated directly or looped indirectly via ‘niche’ factors such as resources or immune cells) outweighs strong interspecific competition. Stronger intraspecific competition leads to net negative feedback on the feasible (’interior’) stable equilibrium, facilitating coinfection.
Two parasite – shared energy (2PE) or – shared immune cells (2PI)
Parasite species that share a resource (energy) or suffer attack from a shared immune system engage in exploitative or apparent competition, respectively (Fig. 2A,B). In exploitative competition (Fig. 2A), a parasite that can survive at the lower equilibrium energy (resource) level (E*) can outcompete the other. For instance, among two clones of a rodent malaria (P. chabaudi), the superior competitor for red blood cells (RBC) always excludes the inferior clone (De Roode et al 2005). Competition for a single, shared resource can lead to a competitive hierarchy among multiple parasite species. For instance, in a malarial-hookworm system, one Plasmodium species (P. falciparum) outcompetes hookworms which themselves outcompete another Plasmodium (P. vivax) for red blood cells (Budischak et al 2018). Analogous rules can characterize host systems with a shared immune system (Fig. 2B). In models, when two exponentially growing parasites share immune cells, the species that can withstand highest density of immune cells (i.e., highest I*) wins via apparent competition (Fenton & Perkins 2010). In a possibly analogous experiment, immune-mediated interactions led to competitive suppression of an avirulent malarial clone, allowing the virulent clone to dominate (Råberg et al 2006). Therefore, exploitative competition for a single resource or apparent competition involving shared immune cells can lead to competitive exclusion alone.
Three aspects of within-host parasite biology can introduce coinfection or priority effects to pure resource or apparent competition. First, if some mechanism imposes self-limitation on both parasites, then they can coinfect even if they share only one immune niche. For instance, donor-controlled reproduction can create strong self-limitation on macroparasite density (Fenton & Perkins 2010). Second, non-linearity in immune clearance rates can generate self-limitation, particularly when clearance accelerates at higher parasite density (Fenton & Perkins 2010). Third, interference competition between parasites can facilitate priority effects or coinfection (Amarasekare 2002; Fig. 2A-B with mutually inhibitory effects between P1 and P2 [grey shading]). If the inferior energy competitor is superior at interference, priority effects ensue. Then, at high enough initial densities, the inferior energy competitor can win instead of being excluded.
However, if the interference also confers a benefit to the interacting species (e.g., if killing competitor larvae and consuming it increases per capita growth rate), then coinfection ensues (following Amarasekare 2002). These general principles for joint interference and exploitative competition (Fig. 2A) likely also apply to immune-mediated apparent competition, a possibility for future exploration (Fig. 2B).
Two parasite – two energy (2P2E; Resource ratio) or Two parasite – two immune cells (2P2I; Immune ratio)
When parasites (Pj) simultaneously compete for two resource or energy sources (Ei) via exploitative competition, infection outcomes can include coinfection and priority effects (Tilman 1982). Consider two parasite species competing for two substitutable resources (like in the resource ratio model: Fig. 2C). Here, coinfection minimally requires each species to trade off their requirements for each energy resource. Then sufficiently intermediate supply ratios must permit each single-species (‘boundary’) equilibria to fall within each competitors’ niches (enabling mutual invasibility). With both conditions met, coinfection arises if each parasite has larger impact on the resource to which its fitness is most sensitive (see Appendix Section 2). In testing the resource ratio model, Lacroix et al. observed a distinct competitive hierarchy, where one plant virus altered success of a competing cereal dwarf virus under nitrogen: phosphorus (N:P) supply. In co-inoculated plants, the cereal dwarf virus infection rates decreased with elevated P supply rate, while the addition of N significantly increased its’ infection rate (Lacroix et al 2014; Smith 2014). This suggests that, at the very least, N:P shapes how each parasite species differentially impacts their resources, setting the stage for coinfection.
Analogous rules likely apply to hosts with two types of immune response (2P2I, the immune ratio model; Fig. 2D). If parasites act as substitutable resources for two generalist immune responses, coinfection arises if each parasite also has larger impact on the immune response to which its’ fitness is most sensitive (Appendix Section 2; derived for case with parasite self-limitation). Modification of this 2P2I framework captures well-studied within-host interactions. First, the host can mount an independent, specialist immune response to each parasite (with no links between Pi – Ij, Fenton & Perkins 2010). Therefore, the separate dynamics of each immune-parasite pair determines coinfection. Similar results emerge for parasites separated spatially (e.g., in different host tissues: Cervi et al 2004; Karvonen et al 2006). Second, two parasites can each be attacked by specialist immune responses that inhibit each other. For instance, in mice, competing parasites interact via T-helper cells, where Th1 attacks intra-cellular malaria, Th2 attacks intestinal nematodes, but Th1 and Th2 inhibit each other (Griffiths et al 2015). Such interference creates negative feedback that can further enable coinfection (e.g., links between I1 and I2 [grey shading]; Fenton & Perkins 2010).
Two parasite – immune cell – energy (2PIE)
Parasites can also engage in simultaneous exploitative and immune-mediated apparent competition (Cressler et al 2014; Ramesh & Hall 2023). Consider common variations of two competing parasites (Pj) that steal energy (E) from hosts and immune cells (I) that kill them (2PIE; Fig. 2E-F; Table S1). In one, only parasites induce production of immune cells (2PIEi, induced immunity). In the other, energy is continuously allocated to maintain baseline immune function, even without parasites (2PIEc, constitutive immunity). This structure resembles a food web in which two prey species share a resource while attacked by a generalist predator (Holt et al 1994; Leibold 1996). This keystone predation (diamond) model anticipates the tradeoffs and niche dimensions that would govern exploitative and apparent competition parasites (Ramesh & Hall 2023). When combined, these forms of competition enable coinfection, priority effects, or competitive exclusion (see below Fig. 3; Ramesh & Hall 2023). Constitutive immunity (Fig. 2F) reduces opportunities (parameter space) for coinfection, lowering parasite burden while maintaining higher energy for other metabolic work (see below; Fig. 9).
The remainder of this review will focus on joint immune and energy (resource) competition in the 2PIE model. Most organisms have some type of (costly) immune defences to fight parasites that steal host energy (reviewed in Lochmiller & Deerenberg 2000; Zuk & Stoehr 2002). Furthermore, principles that govern strength of intra- v inter-specific competition in 2PIE models should apply to others. With the 2PIE model, then, we ask three questions: (a) What within-host feedback drives divergent outcomes of infection? (b) How does this framework offer new ways to resolve results in previous coinfection experiments (Fig. 1A-C)? (c) How can this framework guide future theory and experiments?
Section III: A mechanistic framework linking feedbacks, trait ratios, and infection outcomes
In this section, we construct a framework for within-host competition of parasites that connects trade-offs, nullclines (niches), feedbacks, and traits. As we will learn, the 2PIE model follows the familiar tale of intra- vs inter-specific competition governing coexistence vs priority effects (given a feasible interior equilibrium; Fig. 1D). Although specifics between 2PIEi and the others (Fig. 2A-F) surely differ, this approach lays a general script to generate mechanistic a priori predictions and glean better understanding of coinfection dynamics from within-host niche models. Many of the mathematical details are found elsewhere (Appendix Section 1).
The Model (see also S1, Table S1)
Growth rate of immune cells I (Eq. 2A): In 2PIEi, immune cells (I) increase after attack on parasite j (Pj) at rate , inducing consumption of of energy, E, with conversion efficiency eI of energy into an immune cell after debiting loss rate mI. In 2PIEc, I also increases with allocation of E at baseline rate ab.
Growth rate of parasite Pj (Eq. 2B): The two parasites consume host energy with feeding rate and energy per parasite conversion, . The parasites are lost due to attack by immune cells and die at (shared) background rate mP.
Growth rate of energy, E (Eq. 2C): The host consumes a resource S via a Monod function with maximal assimilation rate fEand half-saturation constant h. (This function merely pays homage to non-linear feeding behaviour of hosts following Cressler et al 2014). That resource, converted to energy within the host (E), is lost at fixed rate r for use by hosts (for metabolic needs). Its net production, then, is f(S) – r E. Additionally, host energy is consumed via induction from immune attack on parasites (proportional to the triple product Pj E I), and from consumption (theft) by parasites themselves. Hence, parasite j ‘consumes’ energy both indirect and directly described, a sum grouped below as fj(I). In 2PIEc, energy is also lost via allocation to baseline immune function at rate ab. The 2PIEi,c models are thus:
Outcomes of competition
Depending on their traits, competing parasites can successfully coinfect, singly infect (via competitive exclusion or priority effects), or produce no infection in a host. Such outcomes can be captured along gradients of, say, supply of resources to hosts, S, and feeding rate of parasites on energy within hosts, (Fig. 3A, Table S1). The lines within the 2D bifurcation diagram represent qualitative shifts of within-host composition determined by minimal requirements. We highlight two such shifts or ‘transcriticals’ (T1 and T2) governing coinfection, but we build to it. Very low S prevents persistence of either parasite, resulting in no infection (a nutritional clearance; Fig. 3A, yellow). Increasing S, now supports only single infection by parasite 1 (P1) because it excludes the other via resource competition (P2, at E-P1 [lighter] or E-P1-I [darker] ‘boundary equilibria’). Now, intermediate S, such that T2 < S < T1, allows parasite 2 (P2) to invade that P1 boundary and P1 invasion into the P2’s boundary equilibrium (E-P2-I). Hence, that region represents mutual invasibility and coinfection (Fig. 3A, orange), situated between regions of competitive exclusion (‘single infection’). At higher S, P1 is excluded allowing P2 to win (E-P2-I; Fig. 3A, dark purple). Notice, increasing both and S flips the geometry of those transcriticals. Now, higher S, such that T1 < S < T2, does not allow P1 invasion into the P2’s boundary equilibrium and vice-versa. Hence, this region represents mutual uninvasibility and priority effects aka founder control (Fig. 3A, grey; see Appendix S1 for details). The offsets show sample dynamics when hosts are exposed to concurrent infection, started at one or two initial densities, in each of these regions (Fig. 3B – F). In the no infection region, insufficient energy caused failed infection, even at high initial densities of parasites (Fig. 3D). In the single infection regions, the winner always excludes the other, even with initially low densities (Fig. 3A, E). In the coinfection region, all trajectories head toward a stable (interior) equilibrium (Fig. 3C). With priority effects, the parasite starting more densely wins and excludes the other (i.e., top: P1 wins; bottom: P2 wins; Fig. 3F).
Mechanism of competition: A within-host “assembly rules” and nullcline approach
These predictions hinge on intra- vs inter-specific competition. Intraspecific competition, mediated indirectly via host energy and immune cells, produces negative feedback loops while interspecific competition yields positive ones (Fig. 1D). When strength of intra-exceeds inter-specific competition, net negative feedback on the feasible interior equilibrium can enable coinfection (coexistence). Conversely, greater inter-than intra-specific competition generates net positive feedback that can produce priority effects. This section actualizes these concepts for 2PIEi using niches and nullclines.
Assembly rules
Broadly, coinfection requires a competition-resistance trade-off: the species that competes superiorly for energy (without immune cell) must remain more vulnerable to immune attack. It also requires intermediate resource supply, S, as consumed by hosts, f(S), falling in between two thresholds (described below: T2 < f(S) < T1). The assembly of a within-host ‘community’ of parasites can be understood via hosts exposed to (i) single species (Fig. 4A, 5A dotted and dashed lines) and (ii) both species (Fig. 4B, 5B solid lines). The single species ‘boundary’ equilibria provide insights into their minimum resource requirements of parasites for energy (Ej*, akin to R* sensu Tilman 1982) and of immune cells for the product of energy and parasites (EPj*). They also highlight maximal immune cell densities supported by each parasite (Ij*, akin to P* sensu Holt 1977). These quantities provide a start for assembly of coinfection.
First, the route to single infection follows resource- and immune-based assembly rules along a nutrient supply gradient, S (Figs. 4B, 5B). Very low nutrient supply (S) meets neither parasite’s minimum energy requirements (E* < Ej*). Hence, very low S prevents infection (i.e., clearance via starvation; yellow). However, with high enough S, E* increases, eventually meeting that minimum (Ej*). Then, Pj can invade and increase with S (pinning E at Ej*; E-Pj regions in lighter). Eventually, Pj’s density crosses a minimal threshold (EPj*) that induces immune activation (darker colors). With higher S still, immune attack reduces density of Pj freeing up more energy (where energy with immune cells is denoted by Ej,I*). While both parasite species follow similar assembly rules, they differ in a key trade-off: P1 is the superior competitor without immune cells (E1* < E2*) but is more vulnerable to immune attack . That competition-(resistance trade-off provides a first prerequisite for coinfection.
Comparison of those key energy and immune quantities yields more conditions for coinfection (Fig. 4B). Along the nutrient gradient, S, the more vulnerable parasite, P1 invades (E-P1). It outcompetes (excludes) the more resistant parasite, P2, via resource competition (since E1* < E2*). At some S following immune activation (EP1*), P1 becomes the superior apparent competitor (I1* > I2*) but inferior energy competitor (E2,I* < E1,I*) – a shift in competitive dominance. Then, with further S still, the E-I niche environment set by P1 enables P2 to invade (at T2: orange).
In a window of coinfection, the within-host niche environment remains constant while parasite densities shift (Fig. 4B). Here, energy and immune cells stay at Eint* and Iint*, respectively. However, with higher S, the structure of this community shifts, with P2 increasing while P1 declines. With this shift, the energy and parasite density needed to keep the immune system activated, Eint*(P1*+ P2*), increases until P2 could support immune activation alone (at EP2*; Fig. 4D). Broadly, in this region, the more vulnerable P1 remains the superior apparent competitor (Fig. 4D, I1* > I2*) but weaker competitor for energy (Fig. 4C, E2* < E1,I* shifts to E2,I* < E1,I *). Conversely, the more resistant P2 is the superior energy competitor but weaker apparent competitor. This biology means that the more vulnerable parasite species also produces more immune cells, the niche factor to which it is most sensitive (see below). That combination imposes a brake on its own growth rate (negative feedback), thus promoting coinfection. With S past this threshold (T1), P2 excludes P1 (via resource competition: E2,I* < E1,I*) but more immune cells lower P2’s density (purple, Figs. 5B-E).
Priority effects can arise with different trait combinations for the two parasites (here: at higher ; Fig. 5B). Along a gradient of nutrient supply to hosts, S, the region of priority effects is approached similarly to the coinfection case: first the more vulnerable parasite invades (when its E1* is met), then the immune system is activated (when its EP1* is reached). However, P1 always remains the superior energy competitor, unlike the flip seen before the coinfection window for coinfection (Fig. 5C). Past T1, more resistant P2 cannot invade until T2 is crossed; with higher S, more resistant P2 excludes P1 via apparent competition (since I2* > I1*; Fig. 5D). Biologically, this means that the more resistant species produces more immune cells, the niche factor to which their competitors are most sensitive (see below), creating positive feedback and thus priority effects. With decreasing S, P1 could not invade until T1 was crossed. Hence, the interior equilibrium in the priority effects window is unstable (a saddle), separating two stable states, dominance by P1 (at its E-P1-I boundary equilibrium) or by P2 (at E-P2-I).
A nullcline (niche) interpretation
Nullclines help to characterize the niche environment that enables it coinfection. At a given nutrient supply S, nullclines denote conditions in which a parasite (Pj) or niche component (energy, immune cells) neither grows nor declines. Combinations of energy (E) and immune (I) cells that fall to the right of a parasite’s nullcline (higher energy, lower immune cells) sit within its fundamental niche (Fig. 6A; P1, blue shading; P2, purple). Then, the impact of each parasite on their within-host niche determines I and E nullclines. The nullcline for E is the combination of P1 and P2 that ‘consumes’ all production of energy at the interior (Eint*). Below this line, energy grows (dE/dt > 0, yellow in Fig. 6B). The Pj-axis intercepts note the density of each parasite, when alone, that consumes all net energy production. The nullcline for I is the combination of P1 and P2 that meets the minimal EP* needed to maintain its activation. All points above the line increase production of immune cells (dI/dt > 0; green in Fig. 6B). At each intercept , the parasite alone supports immune activation (at EPj*) given Eint*.
Coinfection (Fig. 6C - E; black, closed circle): Co-infection requires that parasite nullclines cross, niche nullclines cross, and nutrient supply is intermediate. The nullclines of the parasites cross in E-I space because of their competition-resistance trade-off (Appendix section S1). The shallower slope of less resistant P1’s nullcline means it is more sensitive to I while more resistant P2 is more sensitive to E. When they coinfect, both parasites set their relevant single-species (boundary) equilibria within the fundamental niche of the other. Hence, each single-parasite environment can be invaded by the other (Fig. 6C; P1 equilibrium [blue dot] falls within P2’s niche [purple shaded] and vice-versa). Now P1 is the superior apparent competitor (I1* > I2*) while P2 is the superior energy competitor (i.e., has lower E*: E1,I * > E2,I * [as shown] or E2*; Fig. 4C). Meanwhile, the I nullcline is steeper in P1-P2 space (Fig. 6E). Hence P1 has larger impact on immune cells while P2 has larger impact on energy. Therefore, each parasite exerts largest impact on the niche factor to which it is most sensitive. Based on the nullcline intercepts, coinfection requires that the interior niche would support surplus but deficient than needed to meet immune system’s EP* requirement (Fig. 6E; when ). The coinfection equilibrium, then, combines densities of the superior apparent (P1) and energy competitor (P2) that consumes all net energy production while meeting minimal requirements of immune cells – eliminating those surpluses and deficits.
Priority effects (Fig. 6F-H; white, open circle): During priority effects, those three elements change. First, while parasite nullclines still cross (due to the competition-resistance trade-off), now both parasites (Pj) alone create uninvasible immune (I) – energy (E) niche environments (Fig. 6F; P1 blue dot outside P2 purple shaded niche and vice-versa). Furthermore, in this region P1 is always the superior energy competitor, without (E1 * < E2 *) and with immune cells (E1,I * < E2,I*), while P2 becomes the superior apparent competitor (I2* > I1*; Fig. 6F). Second, the E nullcline is now steeper than the I nullcline (Fig. 6J). Hence, impacts flip: less resistant P1 has larger impact on E while more resistant P2 has larger impact on I. Therefore, each parasite has a larger impact on the niche factor to which its competitor is most sensitive. Thus, the system would provide excess of the more resistant parasite and a deficiency of the vulnerable one enabling priority effects.
Factors enabling coinfection: Feedback loops and trade-offs of traits
Outcomes of coinfection v priority effects depend on feedback and traits that govern them. At the most fundamental scale, co-infection must be ‘feasible’ (set by an environment supporting both species) at an equilibrium having net negative feedback. Such net negative feedback emerges when intraspecific competition exceeds interspecific competition. That competition depends on loops connecting parasites to their niche (Fig. 7A, S1). Loops are chains of interactions linking increased density of a species to growth rates of others, eventually returning to that species (Puccia & Levins 1991). For instance, interspecific competition can be traced starting with an increase in P1. (i) An increase in P1’s density can ‘fuel immune cells’ that suppress its competitor, freeing up energy resources, thereby increasing its own growth rate via higher ‘births’. (ii) P1 can also ‘starve immune cells’ by reducing P2’s density via consumption of E. Through both chains of interactions, P1 indirectly benefits from an increase in its own density via gains in energy, hence birth rate, or lower death rate, respectively. Those two positive (destabilizing) loops then push against two negative loops, (iii and iv) “Pi is attacked, Pj eats”. These latter loops add the stabilizing, consumer-resource-like, intraspecific competitive interactions within which each parasite is enmeshed. For instance, (iii) a small increase in P1 increases immune attack, thereby reducing P1 (P1-I loop) while P2 is slowed by resource consumption (P2-E loop). Then, (iv) those roles reverse, i.e., resources brake P1 while the immune system slows parasite P2. Summed, those two loops (iii and iv) jointly determine the amount of intraspecific competition (negative feedback). Then, the 2PIE model follows the recognizable script of intra- vs inter-specific competition governing coexistence vs priority effects. If the strength of intraspecific competition loops (iii + iv: negative) exceeds that of interspecific loops (i + ii: positive), net negative feedback leads to coinfection. If, instead, interspecific exceed intraspecific loop strength, net positive feedback generates priority effects.
With some rearrangement, these intra- and inter-specific loops governing feedback can more intimately connect to the niche environment and assembly rules revealed above (Fig. 7B, S1). The first involves the difference in ratios of sensitivity of each parasite to immune cells, I, vs. energy, E (denoted by αj; following Pásztor et al 2016). The second involves the ratio of impacts of each parasite on I vs. E (εj). Both ratios can be written as proportional to ratios of traits or traits and minimal requirements (Fig. 7D and E, respectively; see details in Appendix Section S1). Coinfection requires that each parasite has a larger impact on the niche factor to which it is most sensitive. If the more vulnerable P1 is most sensitive to and has highest impact on I, then α1 > α2 and ε1 < ε2. Such symmetry in ratios generates net negative feedback. For priority effects, the competition-resistance trade-off still yields α1 > α2. However, because P1 has highest impact on E (to which P2 is most sensitive), the impact ratios flip, ε1 < ε2. This asymmetry in sensitivity vs. impact ratios produces net positive feedback.
Section IV: A focus on the traits and quantities to measure in the future
Presently, most within-host competition experiments observe patterns of coinfection, then infer mechanisms. The next phase of experimentation should shift toward feedback, minimal requirements, and traits of mechanistic niche models (like 2PIEi). Above, we learned that coinfection requires that each parasite has greater impact on the niche factor to which it is most sensitive (generating net negative feedback); with priority effects, each parasite has greatest impact on the niche factor to which its competitor is most sensitive (yielding net positive feedback). Future research on coinfection dynamics could aim to quantify the traits and minima’s underlying these sensitivity and impact frameworks.
Quantifying traits: The (a)symmetry of the sensitivity-impacts ratios governing coinfection and priority effects can be better understood by delineating the combination of traits (Fig. 7D, Appendix Section S1). For instance, the sensitivity ratios are governed by killing rate of immune cells on parasites , feeding rate of parasites , and the cost of producing a new parasite . On the other hand, the impact ratios depend on how attacked parasites induce energy allocation to immune cells and the per capita ‘consumption’ of energy by parasite j directly and indirectly .
Quantifying key traits and minimal requirements: Experimental tests of sensitivity-impact ratios could also centre on measurement of key quantities such as the minimal requirement of each parasite (E*) for energy and of the immune system for activation (EP*; Fig. 7E, Table 1C). For instance, the difference in sensitivity ratios boils down to killing rate of immune cells on parasites, , and each parasite’s minimal energy requirement without immune cells, Ej* Similarly, the difference in impact ratios depend on per capita ‘consumption’ of energy by parasites, , and the immune cells’ minimal energy-parasite requirement EPj* (Appendix S1). These traits and quantities could be estimated via experiments (in principle).
Despite this call for a focus on traits and minimal requirements, we acknowledge inherent challenges in measuring them. To move forward, experimentalists can leverage interdisciplinary methodologies. For instance, measurement of state variables like parasite (Pj*), energy (Ej*), or immune densities (Ij*) over time can allow estimation of parameters of the model. Alternatively, traits could be measured with radioisotope labelling of resources to track parasite feeding rate (Gomez-Amaro et al 2015), live imaging to track immune killing rate (Galli et al 2021), ICP-MS to measure energetic content of infected host (Cassat et al 2018), etc. As a pay-off, the combination of trait measurements and models generates a priori predictions of infection outcomes and the feedback mechanisms that govern them (notably demonstrated in Budischak et al 2015; Griffiths et al 2015). Ultimately, such approaches may catalyse a new wave of theory-grounded, niche-mechanistic combinations of modelling and experimentation.
Section V: Explanations for divergent infection outcomes using 2PIE
The model of exploitative and apparent competition between parasites, 2PIEi,c, makes predictions (Fig. 8A, 9) that can help to contextualize and interpret previous experiments (Fig. 8B-F). In this next section, we envision how 2PIEi could produce these various outcomes, backing-out mechanism post hoc. Thus, 2PIEi offers a way to potentially resolve otherwise seemingly inconsistent experimental outcomes. Then, it can guide creation of new coinfection models tested in future experiments. Ideally, such efforts would start with parameterized predictions that test divergence of infection outcomes a priori.
Coinfected or not? What within-host factors facilitate or inhibit coinfection? (Fig. 8B, C)
Priority effects
Theory: Stronger inter-than intra-specific competition ensures priority effects (via net positive feedback) where the parasite with sufficiently high initial dose wins (Fig. 8A, B; Table 1B).
Empirical evidence: In a rare, unequivocal demonstration of priority effects, the order of arrival of bacterial (B. burgdorferi) strains determined which won and excluded the other within a host mouse (Fig. 8B). Immune response likely did not explain the priority effects (Devevey et al 2015), so future work can pinpoint mechanisms (e.g., resource competition, interference, etc.) facilitating them (Fig. 2B; Table 1B).
Guiding future experiments: Varying initial parasite dose (dose-dependent assays) provide on possible way to delineate coinfection vs priority effects. When parasites coinfect, dose does not change outcome, despite varying initial P2 and constant initial P1. With priority effects, sufficiently high initial density determines winner (Table 1B). Such an experiment would reveal whether competitive exclusion, or priority effects led to single infection. Even better, measurement of key quantities involved in sensitivity-impact ratios can delineate mechanism (e.g., superior resource v resistant competitor as in Fig. 5).
Competitive exclusion
Theory: Differences in parasite traits can separate exclusion from coinfection. For example, for a given nutrient supply to hosts, a parasite with higher feeding rate (all else equal), can exclude its competitor resulting in single infection (Fig. 8A, C). In contrast, one with lower feeding rate may successfully coinfect. Empirical evidence: Mice infected with malarial strain P1 or P3 (blue) were allowed to compete pairwise against a common malarial strain P2 (purple; Fig. 8C; De Roode et al 2005). In the P3 - P2 pairing, regardless of order or arrival or delay between infections, P3 competitively excluded P2 (Fig. 8C).
Guiding future experiments: Higher competitive ability of P3 may reflect higher feeding rate on host energy. In contrast, if strain P1 had lower feeding rate (hence less competitive), it could stably coinfect with P2 (Fig. 8C). Such hypotheses about resource acquisition traits could be tested in the future (Table 1B).
Coinfection hierarchy and community structure (Fig. 8D-F)
Theory: When parasites coinfect, nutrient supply to hosts (S) can mediate densities of the within-host energy and immunity, and hence relative and absolute densities of parasites (Fig. 8D). Thus, nutrient supply can shift this ‘community structure’ of parasites within hosts. For instance, increasing nutrient supply favours the more resistant P2 over the less resistant P1.
Empirical evidence: This prediction could explain shifts in community structure of parasites in two different systems (Fig. 8E). The first arose with mice infected by two species of gastrointestinal worms (case I; Budischak et al 2015); the second used mosquito larvae infected by a microsporidian and protozoan parasite (case II; Fellous & Koella 2009). Both demonstrate how changing nutrient supply (S) favours one species (purple; positive slope) over the other (blue; negative slope; Fig. 8D). In both cases, relative abundance of coinfecting species can shift (or not) depending on S, too. However, future trait measurements would need to establish that higher S favoured the more resistant parasite as 2PIE predicts (Table 1C).
Guiding future experiments: Much may be learned about coinfection and nutrient supply with focus on individual hosts (reviewed in Ezenwa 2021). Hosts can exhibit high intra- and inter-individual variation in infection outcomes (Merill & Cáceres 2018). To illustrate, parasite densities in the mosquito case ranged among all individuals from coinfection to single infection via exclusion / priority effects to no infection (see reference to shaded regions in Fig. 8D v 8F; replotting of Fellous & Koella 2009). Such variation among hosts could arise from individual differences in resource acquisition. If so, individuals nominally fed the same amount (in a treatment) may fall functionally along different supply points (like in Fig. 8D). Individuals consuming less resources would favour the less resistant parasite (P1) while those eating much more would favour the more resistant one (P2) - yielding exclusion in both cases (Fig. 8D, F). The remaining hosts became coinfected, with parasites reaching different densities among hosts. Future, experimental tests can test for the full range of possible infection outcomes (and hence ‘assembly rules’ leading to coinfection outcomes) at the individual scale, say, along broader nutrient gradients (Table 1C).
Furthermore, niche-based insights can guide more predictive experiments at the within-host scale which can then be scaled to the population linking within-to between-host dynamics. For instance, fluctuation in nutrient supply to host could shift competitive outcomes within hosts that then alters multi-parasite outbreaks at the population scale (Hite & Cressler 2018; Ezenwa & Jolles 2011) or ecosystem processes (Kendig et al 2020).
Coinfection burden (Fig. 9)
Finally, some hosts present higher coinfection burden (total density) than others, hinting at intrinsic host resistance via immune clearance. How does variation in immune response govern coinfection burden?
Theory: In comparison to the model with only immune induction (2PIEi), baseline energy allocation to immunity squeezes parameter space for coinfection and reduces parasite burden (2PIEc, constitutive immunity; Fig. 9A-B). Despite the higher allocation (via ab), reduced burden allows the infected host to maintain more energy for other metabolic work (Fig. 9A-B). At lower nutrient supply (S), ‘no infection’ by parasite shifts from nutritional clearance (E, yellow; Fig, 9A) to combination of immune and nutritional form of clearance (E-I, green; Fig, 9B). The 2PIEc model also predicts that increasing baseline allocation favours the more resistant P2 over P1, eventually excluding it (Fig. 9E).
Empirical evidence: Consistent with these predictions, plant hosts treated with immune-signalling hormone experienced a lower prevalence of a less aggressive parasite, increased burden of infection by a more aggressive parasite, and experienced fewer co-infections (Halliday et al 2018). Thus, hosts with higher allocation to immunity can resist infection more than those with lower allocation, leading to variation among hosts in infection burden (all else equal).
Guiding future experiments: First, one could test such predictions using strains with immune suppression of specific genes or metabolites involved in constitutive immunity (like in Chen et al 2005). That could shift competition from a 2PIEc framework (with ab > 0) to a more 2PIEi like one (with ab = 0 at the extreme). Another test could involve measuring competitive outcomes in host genotypes differing allocation to baseline immunity (Fuess et al 2021; Table 1D).
Section VI: Conclusion
Why do divergent infection dynamics arise within a host? Competing parasites within a host can coinfect (coexist), singly infect (through exclusion, or priority effects), or become cleared (no infection; Fig. 1, 3). Those differing outcomes significantly affect health of individual hosts and even alter population-level disease outbreaks (Mideo et al 2008). Such divergent infection outcomes ought to stem from niche interactions, i.e., with the resources that parasites steal from hosts and the immune cells that kill them (Cressler et al 2014, Graham 2008). Yet, despite two decades of empirical focus on infection dynamics, many of these studies present a collection of disparate results without a synthetic glue. Clearly then, we need a comprehensive framework that explains why parasites coinfect and why they might not.
Here we build a within-host framework of parasite competition based on ecological theory (Fig. 2). Using a two parasite – immune cells – energy model (2PIE), we illustrate general principles governing parasite competition via join exploitative and immune-mediated apparent competition (Ramesh & Hall 2023). Those forms of competition underlie a within-host framework that connects infection outcomes to competitive abilities, nullclines (niches), feedbacks, and traits (Figs. 3 – 7; see details in Appendix Section S1). Notably, we delineate how the interplay of three quantities – minimum resource requirements of (1) parasites for energy (Ej*, akin to R*), and of (2) the product of energy and parasites (EPj*) for immune cells, and (3) maximal immune cells supported by each parasite (Ij*) – provide a start for assembly of coinfection.
We show that coinfection requires that each parasite exerts greater impact on its more sensitive niche factor; that arrangement introduces net stabilizing, negative feedback. With priority effects, each parasite more strongly impacts the factor to which its competitor is most sensitive, leading to net destablizing, positive feedback (Fig. 3, 4).
That framework also provides mechanistic insights into and explanations for some experimental results (Fig. 8-9). First, it explains two ways to hinder coinfection: priority effects can favour early invaders, while competitive exclusion always inhibits one parasite (arising via e.g., fast feeders at low nutrient supply; Fig. 5, 8B-C). Second, increasing nutrient supply in the coinfection region can favour the more resistant over the less resistant parasite, shifting community structure along nutrient gradients (Fig. 8D). Finally, greater investment in constitutive immunity squeezes opportunity (parameter space) for coinfection and reduces parasite burden, freeing up energy for metabolic needs (Fig. 9). Together, the 2PIE models provide interpretation for otherwise puzzling outcomes from various experiments. Although specifics will vary, these principles should apply other niche models involving, say, two energy or immune sources (Fig. 2, 2P2E, 2P2I; see Appendix Section S2). This synthesis of experiments with models underscores need for parametrized trait-based experiments that better predict within-host infection outcomes (Table 1).
Unpacking such within-host mechanisms can improve our understanding of individual health. For instance, deworming trials of hosts coinfected with malaria and gastrointestinal worms show that increased availability of RBCs allow malaria to proliferate within the host, ultimately making them sicker (Budischak et al 2018). Disentangling the mechanisms can then allow for correct course of treatment plan to improve individual health (e.g., malarial drugs or vaccines followed by deworming). Together, our work provides a synthesis of within-host niche-based frameworks, laying a theoretical groundwork for a mechanistic understanding of competition outcomes. We then used model predictions to contextualize key findings from the past two decades of experimentation. With this synthesis, we aim to catalyse a new wave of theory-grounded, niche-mechanistic combinations of modelling and experimentation. Ultimately, such a new wave could help to mitigate the severity and comorbidities of coinfections but also advance the development of preventive drug therapeutics and vaccines, ultimately enhancing individual health.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Competing Interest
Authors declare no competing interests
Acknowledgements
AR was supported by the IU Center for the Integrative Study of Animal Behavior, NSF DEB (1655656), and a MSU Presidential Postdoctoral Fellowship. We thank J Lau and the Klausmeier-Litchman lab for helpful comments on the manuscript; D Brisson, J de Roode, S Budischak, S Fellous for generously sharing their data. AR credits J Coltrane and A Ali Khan for artistic inspiration.
Footnotes
↵# Co-author Email: sprhall{at}iu.edu