A central problem in infection biology is understanding why two individuals exposed to identical infections have different outcomes. We have developed an experimental model where genetically identical, co-housed Drosophila given identical systemic infections experience different outcomes, with some individuals succumbing to acute infection while others control the pathogen as an asymptomatic persistent infection. We found that differences in bacterial burden at the time of death did not explain the two outcomes of infection. Inter-individual variation in survival stems from variation in within-host bacterial growth, which is determined by the immune response. We developed a model that captures bacterial growth dynamics and identifies key factors that predict the infection outcome: the rate of bacterial proliferation and the time required for the host to establish an effective immunological control. Our results provide a framework for studying the individual host-pathogen parameters governing the progression of infection and lead ultimately to life or death.https://doi.org/10.7554/eLife.28298.001
Sick individuals do not all respond to an infection in the same way. One individual may experience mild symptoms and recover easily, while another may suffer devastating illness or even death. A number of factors are often assumed to account for these differences, including the sex, age and genes of the individuals, and differences in the environments the individuals have been exposed to. However, random variations in how an individual’s immune system interacts with the infection could also play an important role in recovery.
Duneau et al. have now studied how genetically identical fruit flies who were raised in the same environment respond to different bacterial infections. This enabled them to develop a mathematical model that describes how a bacterial infection develops in an individual. In an initial phase, the bacteria proliferate freely. If the immune defenses activate in time to control the infection, the number of bacteria in the fly decreases to a constant level and the infection enters a long-term, or chronic, phase. The sooner this happens the more likely it is that the fly will survive. If the immune control happens too late, the infection enters a terminal phase and the fly will die once the number of bacteria increases to a certain level.
The model therefore reveals that the precise time at which the immune system takes control of the bacterial population – termed the “Time to Control” – determines the outcome of the infection. Duneau et al. confirmed this by injecting bacteria into identical flies. A small variation in the Time to Control was sometimes the difference between a fly living or dying. Understanding what controls this apparently random variation is key to understanding infection and potentially developing more efficient treatments for a wide range of diseases – not just those caused by bacteria, but also those caused by viruses and parasites, like HIV and malaria.https://doi.org/10.7554/eLife.28298.002
Two hosts exposed to apparently identical infections may have dramatically different outcomes. In one case, the host may experience only mild symptoms and recover easily, while in another the host may suffer devastating illness or even death. In many pathogenic infections, only a fraction of the infected hosts die from the infection or carry a high pathogen burden, while other hosts may strongly limit parasite proliferation. For example, mosquitoes infected with malaria experience high and low pathogen load (Molina-Cruz et al., 2008), rats infected with Enterococcus faecalis experience bimodal pathogen burden (Frank et al., 2015), and even the abundance of gut microbes exhibits bimodal variation (Lahti et al., 2014). Genetically identical hosts may still experience drastically different infection outcomes (e.g. clonal population of Daphnia magna infected by one strain of bacteria [Duneau et al., 2012]). These observations suggest that the process of infection could itself be intrinsically variable. While such differences can sometimes stem from genetic or environmental variation among hosts, in other cases the difference in outcome appears to be more arbitrary or stochastic, and the relative contributions of genetic, environmental and random individual variation remain unclear.
Variation in within-host disease processes can affect pathogen evolution by modulating the within-host evolutionary dynamic during the course of the infection (Alizon et al., 2011; Mideo et al., 2011). Theory predicts that variation in traits such as bacterial growth rate or host response time to an infection can affect the outcome of the infection and thus the epidemiology of the disease (Alizon et al., 2011; Fenton et al., 2006). To better understand the selection pressures that shape host-parasite interactions, it is important to empirically open the ‘within-host black box’ and characterize the different phases of an infection and to define the evolutionarily relevant parameters in the host and parasite (Alizon et al., 2013).
Drosophila is a prime model to study host-microbe interactions (Buchon et al., 2014). Drosophila relies on both humoral and cellular responses to fight infection (Lemaitre and Hoffmann, 2007). The humoral response consists of the secretion of antimicrobial peptides (AMPs) into the hemolymph (insect blood) and the activation of the melanization cascade, both of which limit microbial proliferation. The cellular responses include encapsulation of foreign bodies and phagocytosis by dedicated immune cells. Antimicrobial peptides are constitutively expressed in a restricted set of barrier tissues at low levels but are massively upregulated in the fat body in response to infection. Two major signaling pathways, the Toll and Imd pathways, drive the induction of AMPs in the fat body. Importantly, we and other groups have observed strong inter-individual variation in the outcome of Drosophila infected with a wide range of bacteria and viruses (Chambers et al., 2014; Clemmons et al., 2015; Ferreira et al., 2014; Howick and Lazzaro, 2014; Kutzer and Armitage, 2016), which suggests that Drosophila is a good model to study the stochastic nature of infection.
In the present work, we aim to identify key parameters of the intra-host dynamic that underlie the inter-individual variation in the outcome of infection. We first demonstrate that infection of Drosophila melanogaster by a series of bacterial pathogens leads to highly variable infection outcomes. Notably, infections with the Gram-negative Providencia rettgeri and the Gram-positive Enterococcus faecalis result in one of two binary outcomes: some individuals die with a high pathogen burden, whereas others survive indefinitely with fairly asymptomatic and low-burden persistent infections. Different individuals of the same age and genotype may experience either outcome, even when co-housed in the same rearing vial. We find that bacterial burden upon host death did not correlate with the time post-injection at which death occurs, and that the lethal burden varies across bacterial species and strains. Surviving flies sustain a much lower pathogen burden, which also varies across bacterial species and can depend on the initial inoculum and on the immune response. Using a variety of experimental manipulations, we determined that inter-individual variation occurs as a consequence of the interplay between bacterial growth rate and the activation of the host immune response. We developed a theoretical model that describes the probabilistic nature of bacterial growth and thus host survival based on three key factors: the net rate of bacterial proliferation (µ), the timing of effective immunological control (Tc), and an inferred threshold pathogen density (ntip) that must be reached to enable the bacteria to switch to a lethally acute infection instead of a chronically persistent one. Our model predicts infection dynamics accurately and suggests that inter-host variation in survival must originate in the ability of the pathogen to reach the ntip threshold before effective immune control is established. We observe empirically that the quantitative activation of the immune response varies among individual hosts, suggesting that inter-individual stochasticity early in infection might lead to differences in the probability of ultimate survival. Altogether, our results illustrate the mechanisms underlying the variable nature of infection outcome and provide a framework for studying the distinct parameters that govern the progression of infection and lead ultimately to life or death.
Infection with different pathogens produces dramatically different outcomes that range from the certain death of all host individuals to the nearly certain survival of all host individuals. In order to define the range of outcomes, we infected male Drosophila melanogaster with a panel of seven bacteria by injection of a controlled dose into the abdominal cavity (Figure 1A). Infection with Providencia alcalifaciens and Serratia marcescens Db11 was lethal to the flies in a short timeframe, killing all individuals within one day of infection. Conversely, nearly all individuals survived infection with E. coli and Erwinia carotovora Ecc15 (formally known as Pectinobacterium carotovorum carotovorum 15) (Figure 1A). Infections with different natural pathogens of Drosophila including P. rettgeri, E. faecalis, P. burhodogranariea resulted in the death of only a fraction of the infected population (Figure 1A). It is commonly assumed that inter-individual differences can stem from genetic differences or environmental variation, including nutrition, gut microbes or rearing conditions. We therefore controlled for these sources of variation by infecting same-aged, genetically identical hosts (isogenic lines of the Drosophila Genome Resource Panel and w1118), reared in a common environment with a controlled microbiota and given experimentally identical infections (Figure 1B). Even under experimentally identical conditions, we continued to see dual outcomes of infection. In addition, the relative probability of the infection achieving either outcome varied among genotypes. Inter-individual variation could also stem from small variation in fly handling at the time of infection. For instance, differences in exposure to the CO2 used for anesthesia during the experiments could influence the probability of death upon infection (Helenius et al., 2009). However, we found that even 15 min of exposure to CO2 prior to or after injection did not affect the chance of survival relative to 2 min of exposure (Figure 1—figure supplement 1). These results demonstrate stochastic variation in the outcome of infection (life or death) even among genetically identical individuals, although host and bacterial genotype combine to determine the probability of surviving infection.
One could reasonably predict that individuals that survive their infections may do so by having cleared the infecting pathogen. To test this, we quantified bacterial loads in individual flies that were still alive seven days after injection with E. coli, E. carotovora, E. faecalis, and P. rettgeri. At this point, mortality curves had plateaued and little to no further death was observed (Figure 1A). We found that none of the bacteria had been entirely cleared from the host and nearly all individual hosts still had measurable loads of bacteria (Figure 1C). In the case of the most benign pathogens, which did not lead to the death of any host individuals, a few host individuals either cleared the infection or sustained loads below our detection threshold of 10 bacteria per fly (~22% for E. coli and ~8% for E. carotovora). However, most surviving hosts developed persistent, chronic infections. In addition, the bacterial load sustained during the persistent phase of infection was constant but distinct for each bacterium (e.g., Figure 1—figure supplement 2A). We termed this parameter the ‘Set Point Bacterial Load’ (SPBL), in analogy to the commonly used ‘Set Point Viral Load’ for chronic viral infections (Fraser et al., 2014). Interestingly the SPBL for E. coli and E. carotovora were significantly lower than those from E. faecalis and P. rettgeri infections, which also differed between each other (Figure 1C). The SPBL for E. coli and E. carotovora, our two least virulent bacteria, was a few hundred bacteria per infected host. In contrast, the SPBL for pathogens of intermediate virulence, such as P. rettgeri and E. faecalis, centered around a few thousand bacteria, and no single fly ever cleared the infection. In addition, we found the SPBL for P. rettgeri to be significantly different among three D. melanogaster host genotypes (Figure 1C), increasing with the size of the initial inoculum (Figure 1—figure supplement 2B).
We next examined the individual hosts that died from their infections. In the case of P. rettgeri, only a subset of infected hosts die. Mortality begins as soon as 16 hr post-injection and continues over the course of a few days before survivorship plateaus (Figure 1A,B). We quantified the Bacterial Load Upon Death (BLUD) by sampling individual flies within 30 min after their death, before the confounding effect of post-mortem bacterial proliferation. We compared those loads to pathogen burdens of still-living individuals sampled from the population all along the course of infection (Figure 2Ai) or at a unique time point where 5% of the population died at once (Figure 2Aii). The P. rettgeri loads in recently-dead individuals were 223.4 ± 21.6 (mean ±sd) bacteria/fly, regardless of the time post-injection at which the host died (Figure 2B and Figure 2—figure supplement 1A). These loads were consistently 100-fold higher than those of living individuals (215.8 ± 23.8 bacteria/fly; Figure 2A). The narrow range of bacterial load upon death and the fact that live flies carry fewer bacteria than recently dead flies (Figure 2Aii) suggests that the BLUD is a lethal load and cannot be sustained in living individuals. We found that there is a BLUD associated with every bacterium and that the BLUD differed significantly between bacteria (Figure 2C), indicating that this threshold is a property of the pathogen. Surprisingly, the BLUD did not reflect the time it took various pathogens to kill the host. For example, the BLUD for S. marcescens, which killed wild-type flies within 8 hr of injection, was 220.2 ± 20.7 bacteria/fly. This is lower than that of P. alcalifaciens (224.3 ± 20.5), which killed within 24 hr, but higher than for S. aureus (218.1 ± 20.7), which killed within 28 hr (Figure 2C). Because the BLUD varies across bacteria but is constant for a given bacterial strain in a given environment, we propose that it provides a novel measure of bacterial pathogenicity or toxicity to the host.
Reciprocally, we reasoned that the BLUD also represents a measure of host disease tolerance corresponding to the maximal bacterial load that can be sustained before the individual dies from the infection. As a measure of disease tolerance, the BLUD might vary among hosts. We therefore next measured both inter- and intra-specific variation in the BLUD after infection with P. rettgeri. We first assessed the BLUD of P. rettgeri-infected flies belonging to four different Drosophila species (D. sechellia, D. erecta, D. mojavensis and D. ananassae) and to four different D. melanogaster wild-derived isogenic lines (RAL-714, RAL-359, RAL-138 and RAL-882). We found that the BLUD varied significantly across Drosophila species (Figure 2D) and across different D. melanogaster genotypes (Figure 2E). To our surprise, however, the BLUD did not depend on the initial inoculum starting the infection (Figure 2—figure supplement 1B). We tested this by challenging wild-type D. melanogaster lines with P. rettgeri infection at five different doses ranging from ~300 (OD600 = 0.01) to ~150,000 (OD600 = 5) bacteria injected. Flies injected with higher doses died earlier than flies injected with lower doses (Figure 2—figure supplement 1C–D) but BLUD did not correlate with the dose injected (Figure 2—figure supplement 1B). These data affirm that hosts die at a lethal burden that does not depend on their initial infection dose or on their time to death.
Having described the BLUD as a measure of disease tolerance, we sought to determine whether the immune response, the main mechanism of bacterial elimination (i.e. resistance), affected the BLUD. We did not detect significant differences in the BLUD between any of the immune-deficient hosts and their corresponding wild-type genotypes (Figure 2F). This result supports the idea that the BLUD is an exclusive measure of tolerance and not resistance. Collectively, our results indicate that bacterial infections of Drosophila result in binary outcomes. Either the host dies from acute infection with a stereotypical bacterial burden that is independent of dose, the BLUD, or the host survives but sustains a chronic infection at low and also stereotypical bacterial burden that does depend on initial dose, the SPBL. These binary outcomes are general to D. melanogaster and are observed in all wild-type host genotypes and species that we tested, although the relative proportion of hosts entering either state may vary with pathogen identity and host genotype (Figure 1A and B). Although the BLUD is a key parameter to describe bacterial virulence and host disease tolerance, it does not explain the observed inter-individual variation in the outcome of infection (death versus chronic persistence).
We have shown that individual hosts may either survive or die from their infections, and we hypothesized that among-individual variation in survival might result from differences in bacterial growth dynamics at early stages of infection. In order to test this hypothesis, we initiated an in-depth study of bacterial growth in response to three categories of pathogens: those that quickly kill all infected hosts, those that kill none, and those that kill only a fraction of the infected population. We measured the bacterial load in individual flies at regular one-hour intervals for the first 16 hr after injection (Figure 3). Upon infection with two avirulent bacteria, E. coli and E. carotovora, bacterial burdens inside individual flies decreased soon after injection (Figure 3A), suggesting that these bacteria are controlled within a couple of hours by the host, even though they are not completely cleared. After injection with bacteria that kill all infected hosts in few hours (P. alcalifaciens or S. marcescens), the pathogen burdens increased monotonically until all individuals were dead (Figure 3B), suggesting that the host fails to control the infection. The pattern was very different for pathogens killing only a fraction of the infected host population: E. faecalis, P. burhodogranariea, P. rettgeri (Figure 3C–E). For these, we found that bacteria proliferated homogenously in all individuals during the early phase of infection (up to around 6–8 hr), but that subsequently two groups of hosts could be distinguished. Hosts of one group carried high and increasing bacterial loads, and we hypothesized that these individuals are fated to die. Hosts in the other group exhibited a smaller bacterial load that stabilized over time, and we hypothesize that these are the individuals who will sustain chronic infections. Remarkably, the in vivo proliferation of P. rettgeri in the early phase of infection is equivalent to the rate of in vitro growth in LB medium (Figure 3—figure supplement 1), suggesting that the bacteria are growing largely without constraint in this early phase. The same is true of the always-lethal P. alcalifaciens and S. marcescens (Figure 3—figure supplement 1). Proliferation continues at this same rate in the host group with the higher bacterial load up until the BLUD is reached, at which point the host dies. These results are consistent with the hypothesis that host variation in the outcome of infection might stem from differences in the ability of the fly to control infection in the early stage.
Although perhaps unlikely, one possible explanation for the observed data would be that a subset of hosts infected with P. rettgeri, E. faecalis, P. burhodogranariea do not survive their infections because a mutation that provides resistance to the host immune system occurs in the bacterial population, allowing these bacteria to continue to grow instead of being controlled. If this were true, the bacteria that kill the fly would have evolved hyper-virulence and be able to outcompete their wild-type counterparts. Conversely, the flies that survive their infections might conceivably do so because the bacteria infecting them sustained loss-of-virulence mutations. To test these hypotheses, we sampled bacteria from flies bearing either low (i.e. close to the SPBL) or high bacterial load (i.e. close to the BLUD; Figure 4Ai) at 12 hr post-injection and used them to infect new cohorts of flies. In both cases, we found that bacteria isolated from infected hosts were as likely or more likely than the original bacterial clone to kill new hosts, although we see no difference in mortality upon infection with bacteria from high-load versus low-load flies (Figure 4Aii). These data collectively provide no evidence that the probability of a host controlling or failing to control bacterial infection results from heritable changes in the bacterial population during infection.
We next hypothesized that the host immune response might be primarily responsible for regulating pathogen burden, and that variation among individual hosts in immune activity might determine whether infections proceed to acute lethality or settle into chronic persistence. To test the hypothesis that the existence of a binary outcome to infection depends on the immune system, we studied P. rettgeri growth during an in vivo infection and the probability of death in hosts deficient for either the cellular or humoral immune responses.
Engulfment by macrophage-like cells, called phagocytes, is considered the first line of immune defense once the parasite is inside an insect. We thus tested whether phagocyte activity predicts inter-individual variation in bacterial growth. We infected hosts in which phagocytes were genetically ablated with P. rettgeri and compared them to infected control hosts (Figure 5A). Individual hosts separated into high-load and low-load groups over the first 16 hr of infection both in the presence and in the absence of phagocytes. We therefore concluded that phagocytosis is not a determinant of the switch between lethal and persistent P. rettgeri infection under our experimental conditions.
We next investigated the contribution of the humoral response to inter-individual variation in bacterial control. Gram-negative bacteria are mainly controlled by the Imd pathway, a major regulator of AMP production (Buchon et al., 2014). We tested whether flies deficient for the Imd pathway (DreddEP142 mutant, Figure 5B) show inter-individual variation in P. rettgeri proliferation. No individual DreddEP142 mutant fly could control infection and all mutant hosts died with a high pathogen burden (Figure 5B). This demonstrates that the Imd pathway is required to control P. rettgeri infection and shift the infection to a persistent state in a subset of hosts.
The Toll pathway provides the primary response against Gram-positive bacteria but can also be activated by host damage and pathogen-derived virulence factors (El Chamy et al., 2008). We infected Toll-deficient (spzrm7 mutants) and wild-type hosts with P. rettgeri and observed the emergence of discrete high-load and low-load groups of flies of both genotypes during the first 16 hr of infection (Figure 5C). This demonstrates that Toll-deficient flies are sometimes able to control P. rettgeri infection. We next tested the role of the Toll pathway upon infection with a Gram-positive bacterium, E. faecalis, which also results in either acute lethality or chronic persistence in wild-type flies (Figure 3C). Without exception, flies deficient for the Toll pathway (spzrm7 mutants) failed to control E. faecalis infection and all hosts died with high bacterial loads (Figure 5D). These data collectively demonstrate that the specific and appropriate arm of the humoral immune response is critical for the control of bacterial infection and establishment of persistent infection.
Altogether, our data indicate that the capacity to control infection versus allowing it to proceed to the lethal state is determined by the action of the specific and appropriate arm of the immune response, but that the lethal pathogen burden is independent of immunity. This reinforces the idea that the BLUD is a readout of disease tolerance, but that the probability of the ultimate outcome of infection is mechanistically distinct and determined by host immunological control.
We have empirically demonstrated that even genetically identical hosts differ in their ability to control infection, which ultimately leads to binary outcomes of death at high pathogen burden or survival with chronically persistent infection. In order to quantitatively define how the kinetics of bacterial proliferation and host response determine ultimate infection outcome, we developed a model to estimate bacterial growth within each individual host (illustrated in Figure 6A). Our model assumes that in the early phase of infection, the bacteria grow unchecked within the host at rate μ. Empirical support for this component of the model comes from the observation that bacteria proliferate at the same rate as in in vitro culture over the first few hours of infection in a wild-type host and until death in an immune deficient host (as shown in Figure 3—figure supplement 1). We assume that there is a time post-injection where the immune response of an individual fly becomes sufficiently active to restrict bacterial proliferation, which we term the ‘time to control’ (Tc), and we define t̄c as the average time to control of a given population. This is based on the empirical observation that an intact immune response is required for establishment of a persistent infection (see Figure 5B) and that the timing of our observed divergence in host trajectories corresponds roughly to the time at which the D. melanogaster humoral immune response begins to become active (Lemaitre and Hoffmann, 2007). We also observed that infections that kill only a fraction of the hosts reach a state where flies either have a high bacterial load, in which case they die, or a low bacterial load, in which case they survive. We model this divergence between the two types of bacterial kinetics by hypothesizing a threshold bacterial number (ntip) that, if reached prior to control, results deterministically in ultimate host death at high bacterial load (the BLUD).
By applying our model to the empirical data, we can obtain for each individual host a probability corresponding to the likelihood that control is effective (see Figures 3 and 5, where white dots are infections that are most likely uncontrolled while red dots are most likely to be controlled). This probability is computed as the probability that bacterial load does reach ntip before control becomes effective. Using this computation, we estimated the overall probability that an individual host controls an infection by a given bacterium. We found that the estimated probability of control in our model, which is obtained from the measurements of bacterial load in different flies at different times along the infection (i.e. there was no repeated measurement on the host), correlates well with the empirically observed proportion of hosts that are killed (Spearman rank correlation test: S = 3.02, r = 0.97, p-value=8.17e-6; Figure 6B). This is both an indication that the model adequately describes our data and a confirmation that the probability of death is determined by the early phase of infection as described by the model.
One parameter that our model allows to estimate from experimental data on early growth is the average time to control (t̄c; Figure 6C). We hypothesized that ‘priming’ the immune system prior to infection should dramatically reduce this parameter and therefore yield much higher probability of controlling the infection. To test that, we infected transgenic flies that had constitutive activation of the Imd pathway with P. rettgeri and monitored bacterial growth. We found that bacterial loads in hosts with an immune system that was genetically activated prior to infection began to decrease immediately after injection, demonstrating effective control at the moment of infection (t̄cc0) (Figure 7A). Surprisingly, however, the bacteria grew to acute lethality in a small proportion of individuals with constitutive Imd pathway activation (Figure 7A), suggesting that the Imd pathway alone is not always sufficient to control P. rettgeri.
We compared the average time to control (t̄c) estimated from the model with experimental measurements of the time at which host immunity impacts bacterial load after infection with the same bacteria. We quantified proliferation of E. carotovora and P. rettgeri in wild-type and Imd pathway mutant hosts (Figure 7B) and the growth of E. faecalis in wild-type and Toll pathway mutant hosts (Figure 5D), using the logic that bacterial loads should become measurably different between wild-type and immune deficient flies soon after control by the immune response occurs in the wild-type host. Using this method, we found that the effect of the immune response is detectable starting at 2 hr post-injection for E. carotovora, at 8 hr post-injection for P. rettgeri, and at 2 hr post-injection for E. faecalis (Figures 7B and 5D). These empirical measurements of t̄c ranges overlap with the confidence intervals of t̄c that our model predicts for these infections (Figure 6C).
We found variation in t̄c across bacterial infections, and we next hypothesized that this variation, because it impacts the probability of control (Pc), could partly determine the outcome of infection. In our model, the probability of control also depends on parameters that determine early bacterial growth. To analyze the contribution of t̄c and bacterial growth rate to the outcome of infection, we compared the likelihood (i.e. how well the model fits the data) of a complete model fitting the data through the estimation of 11 parameter values for each type of infections to the likelihood of models where some of the parameters are held constant for all infection types. For example, to estimate the role of the parameter t̄c, we compared likelihood of our full model, where each type of infection has its own value of t̄c, with a model fitting our data by estimating a common t̄c value for all infection types. A significantly better fit of the full model reveals a significant role of the parameter which was held constant in the second model. We performed this analysis for all the parameters involved in modeling the control of infection (Control model, i.e. t̄c, Vc, ntip) or the parameters involved in modeling early bacterial growth (Growth model, i.e. tlag, µ, σb) to test their relative impact on how accurately the model predicts survivorship (Table 2). We found that both sets of parameters involved in modeling the control of infection and early bacterial growth both have a strong impact on adequately predicting bacterial load (Table 2). We found that variation in survivorship is best explained by variation in t̄c when only bacteria displaying intermediate survival are considered (i.e. when excluding E. coli, E. carotovora which have very low growth and mortality rates, and S. marcescens and P. alcalifaciens which have high growth and mortality rates, see Figure 6B). Conversely, when all bacteria were included in the analysis, differences in survivorship are best explained by differences in early bacterial growth (Table 2). These results suggest that the probability of surviving infection is affected by the average time at which control occurs (t̄c), but that this effect can be masked by large differences in growth rates (µ).
Our analyses so far indicate that variation in infection outcome is strongly dependent on t̄c, which should be determined by both the kinetics of the host immune response and the bacterial sensitivity to host immunity. For infection by a given bacterium, with a known and constant t̄c, our model predicts that inter-individual variation in infection outcome must stem from variance in time to control (Tc), quantified by Vc. Setting this variance to zero for all infections indeed produces a probability of control which is either zero, with all individuals dying from infection, or one, with all infections being controlled. Since we see empirical variability in control of some infections, we hypothesized this might arise from biological variance in Tc.
We empirically tested the hypothesis that the speed and magnitude of Imd pathway induction might vary among individual flies, potentially resulting in variation in Tc that leads to variation in ultimate infection outcome. To test this, we first determined the kinetics of Imd pathway activation after infection with P. rettgeri by quantifying mRNA levels of the Diptericin gene, which encodes an antimicrobial peptide and provides an excellent read-out of Imd pathway activity (Buchon et al., 2014; Lemaitre et al., 1997). We began by measuring Diptericin expression in pools of 20 flies Canton S. Diptericin expression was extremely low at the time of injection (Figure 7C, left), but a measurable induction was detected at 2 hr post-injection and its level continuously increased during the first 8 hr post-injection (up to ~100 fold induction compared to the start of the infection). Thus, the Imd pathway is activated transcriptionally very early upon infection, well before we detected any impact on bacterial proliferation. To test whether inducibility is variable among individual hosts, we quantified Diptericin expression level in single flies at 4 hr post-injection (Figure 7C, right). Time points equal or earlier than 8 hr were chosen because they occur before t̄c and before we begin to see inflated variance in the bacterial load of individual flies. At later times of infection, higher bacterial number leads to increased immune system activation in a positive feedback loop, so variation in AMP gene expression level becomes an indicator of the differences in bacterial loads instead of an estimation of the variance of induction in response to a unique bacterial load. At 4 hr post-injection, however, immune induction can be detected (Figure 7C) but immune control is not yet effective (Figure 7B) and no divergence in bacterial load is yet observed (e.g. Figure 3D and E). At 4 hr post-injection, we detected strong among-individual variation in Diptericin expression, ranging from little induction to strong immune response, with approximately 5-fold difference between extremes (Figure 7C), which may result from experimentally undetectable differences in handling, infection, or developmental history. We then did paired measurements of AMP expression and bacterial loads from single flies at 8 hr post inoculation, which is just before the predicted t̄c. The Imd pathway is activated and Diptericin is expressed primarily in the D. melanogaster fat body, most of which is located in the abdomen, and we determined that the bacterial burden in the head of a fly is tightly correlated with the load estimated from the thorax and abdomen: (r = 0.84, p=8.8e-9; Figure 7—figure supplement 1), so we tested whether pathogen burden estimated from the heads of individual flies could be predicted by the Diptericin expression measured from abdominal RNA extractions. We found no significant correlation between bacterial load and the level of Diptericin transcripts 8 hr after injection (Spearman correlation: r = 0.06, S = 4225.88, p=0.75). This indicates that empirical variability in Imd activation is independent of variation in bacterial load at the early stage of infection. The disconnect may allow the infections of some hosts to grow through the critical threshold ntip prior to t̄c, ultimately leading to host death, while the seemingly identical infections of other hosts are effectively controlled.
We evaluated infection of Drosophila by multiple bacteria that cover a range of pathogenicity, ranging from completely lethal infections that kill all hosts to fairly benign infections where most or all hosts survive. Most interestingly, we demonstrate binary outcome of infection with pathogens such as P. rettgeri or E. faecalis. Even among hosts that are identical in genotype, sex and age and are raised in a common environment, a fraction of hosts die with high bacterial burden while the remainder survive indefinitely with persistent low pathogen load. We identify a set of key parameters that determine which of these outcomes occurs and we propose a general model framework to describe the dynamic interplay between host and pathogen. We find that host death is associated with a bacterial load, the BLUD, which varies with host genotype and pathogen strain but does not vary with inoculum dose. The BLUD does not correlate with the time a host will die from infection. Instead, we suggest that inter-individual variation in the probability of death stems from variation in the time it takes for each individual fly to establish effective immunological control of the infection (Tc), with the understanding that this must occur before the pathogen reaches a critical density (ntip) if control is to be achieved. Hosts that survive their infections do so carrying a fixed bacterial load, the SPBL.
Based on the quantification of bacterial load in the host, we have built a model that describes inter-individual variability in bacterial growth (Figure 6A). In an early phase of the infection by pathogens, bacterial proliferation occurs at a rate similar to that of in vitro growth, suggesting that growth is essentially unconstrained in the first hours of infection. After a brief delay, the host immune response becomes sufficiently active so as to limit or perhaps even reverse bacterial growth. This time, which we parameterize as t̄c in our model, will certainly depend on both the kinetics of immune induction by the host and the sensitivity of the pathogen to the host immune response. In a following phase of our model, the infection resolves in one of two ways. The bacteria may continue to proliferate at the unconstrained rate until the host dies at a lethal load, the BLUD, or the bacteria may cease proliferation and settle into a phase of long-term persistent infection with a burden defined as the SPBL. In principle, the SPBL could be complete elimination of the bacteria, although we very rarely see that with D. melanogaster hosts. We infer that the probability of an infection entering into the lethal versus persistent state is defined by the probability that control is effective before bacteria grow above the critical threshold density, ntip. Thus, the probability of lethal infection is most influenced and therefore predicted by the rate of bacterial proliferation (µ) and the average time required to mount an effective immune response, t̄c.
The simplicity of our statistical model yields some crucial insights and implications. First, a complex and apparently stochastic organism-level phenotype (survival versus death) can be satisfactorily predicted by modeling the probability of transitions between a few discrete states of infection (early growth, uncontrollable growth, and decrease/stabilization after control). In the course of infection, each may occur in discrete temporal phases (termed ‘early phase’, ‘resolution phase’, ‘terminal phase’ or ‘chronic phase’ in Figure 6A) that have distinct and possibly deterministic outcomes, as those described in (Duneau et al., 2011; Ebert et al., 2016; Hall et al., 2017; Levin and Antia, 2001; Schmid-Hempel and Ebert, 2003). Furthermore, each of these states and the parameters that determine them has a distinct mechanistic underpinning and may be differentially subjected to natural selection. For example, host evolution of shortened time to control through increased inducibility of the immune response would have a different genetic basis than one that involved reduced bacterial proliferation through host sequestration of nutrients (e.g. metal ions or carbon). Both mechanisms might result in apparent ‘resistance’ to infection, both would contribute to increased probability of survival from infection, and each would carry distinct costs or tradeoffs. From the pathogen’s perspective, changes in its growth rate or sensitivity to host immune responses would also alter the probability that they become lethal versus persistent, also with distinct mechanistic bases and costs or tradeoffs. Each predictive parameter in the model can be experimentally manipulated, and the entire model can be theoretically evaluated to identify evolutionary optima in different environments and under different transmission requirements. We note that this model can be easily extended to infection by other classes of pathogen, to specific tissues such as the gut, or to hosts that display secondary immune responses (e.g., antibody-mediated acquired immune responses).
Interestingly, the two endpoint parameters of infection, bacterial load in survivors (SPBL) and in hosts that succumb (BLUD), were remarkably invariant among individual hosts of a same genotype even though they varied substantially across infections by different pathogens, host species and genotypes. The mechanisms underlying the death of a fly upon infection have been curiously overlooked, and the nature of damage to D. melanogaster from bacterial infection remains elusive (Dionne and Schneider, 2008). We find that flies die at a specific bacterial load that does not depend on initial infection dose or eventual time to death. This implies that death is unlikely to be a direct consequence of accumulated infectious damage, since longer term exposure to low pathogen burdens might trigger as much damage as short-term exposure to high bacterial burden. Instead, our results are more compatible with a model for death occurring as a consequence of sepsis or multiple organ failure driven by pathogen load (Baue, 1975). Of course, the lethal load of different pathogens would vary depending on the toxicity of the particular pathogen. In that sense, we can argue that the BLUD is a measure of pathogenicity of the bacterium, as it represents the minimal density of bacteria that kills a host. We hesitate to use the word ‘virulence’ though, because the BLUD does not correlate with the speed at which pathogens kill or with the probability of host death. Reciprocally, the BLUD also represents a measure of host disease tolerance (Soares et al., 2017) as it corresponds to the maximal bacterial load that can be sustained in the host without dying from the infection.
Our results are somewhat inconsistent with those of Haine et al., 2008, which documented rapid cellular elimination of bacteria injected into Tenebrio molitor and posited an only minor role for the humoral response. In contrast, our study indicates that the timing and intensity of AMP production is the most crucial factor for controlling bacterial infection in Drosophila. One intriguing hypothesis for reconciling the two observations is that adult Drosophila have only a small number of hemocytes due to loss of hemocytes upon aging (Guillou et al., 2016), while this might not be the case in Tenebrio. Thus, the relative importance of cellular versus humoral immunity may differ between the two insects. Additionally, different bacterial pathogens were employed in the two studies, and different microbes may be differentially sensitive to alternative arms of the immune response so generalizations should be made with caution.
We found that complete bacterial clearance rarely happens in Drosophila. With all seven pathogens in this study, in no case were the bacteria fully cleared from all hosts even a week after the initial infection, even in cases where 100% of hosts survived the infection. Persistent infections stabilized at a fixed load for each pathogen, the SPBL. On rare occasion, individual flies in our study would carry persistent infections for several days before suddenly dying with a pathogen burden that had reached the BLUD. This suggests that there are some conditions under which a persistent infection can re-emerge as an acute infection. This is reminiscent of chronic infections by multiple human pathogens, including HIV, eukaryotic parasites such as Toxoplasma gondii, and bacteria such as Salmonella enterica, Mycobacterium tuberculosis or Streptococcus pneumoniae. These parasites are considered to be ‘specialists’ at chronic infection, yet our results suggest that qualitatively similar phenomena can be obtained from infection with a broad range of generalist bacteria. As the SPBL varies with both host and pathogen genotype, it is presumably subject to selection in both. In that context, it is worth noting that forcing an infection into the persistent state is effectively a disease tolerance strategy from the host perspective, yet the pathogen burden borne may still carry lifelong fitness cost. From the pathogen’s perspective, either persistence or acute lethality may be favored depending on the conditions of transmission. Moreover, selection for high growth rate due to competition between bacterial genotypes within a host (Alizon et al., 2009; 2013) may trade off against the capacity to establish persistent infection.
Levin and Antia (2001) suggested that ‘us and other living organisms are little more than soft, thin-walled flasks of culture media’. With the present work, we argue that hosts are slightly more complex and that potentially minute variations in host or pathogen physiology during the early phase of infection can have dramatic effects on the ultimate outcome. We have defined three parameters (t̄c, μ, nTip) that are sufficient to predict ultimate infection outcome, although we still do not know whether variation in those parameters is due to micro-environmental variation, uncontrolled plasticity in developmental history, somatic mutations, or other factors that are uncontrolled or uncontrollable in experiments (e.g. depth of penetration of the needle during injection, site at which bacteria accumulate etc.) and in nature (e.g. time since last meal or mating, psychological status of the fly etc.). Nevertheless, what is becoming clear is that small differences in pathogen infectivity and host immunological control, especially in the early stages of infection, may manifest as large differences in the outcome of infection.
Drosophila melanogaster were reared on glucose-yeast medium (82 g/L yeast, 82 g/L glucose, 1% Drosophila agar, supplemented with 2.5 mg/L methylparaben and 10 mL of a solution of phosphoric and propionic acid: 41.5 ml Phosphoric Acid + 418 mL Propionic Acid + 540.5 mL distilled water). At day two after eclosion, adults were isolated in groups of five males and five females. All experiments were conducted with mated males 5 to 8 day post-eclosion. Rearing and experiments were conducted at 25°C (±1°C) with a 12 hr light/dark cycle. Canton S, Oregon R and w1118 were used as wild-type laboratory strains. Eight fully isogenic lines from the Drosophila Genome Reference Panel (DGRP: RAL-85, RAL-138, RAL-359, RAL-882, RAL-594, RAL-707, RAL-712, RAL-714) were randomly chosen and used to assay the role of inter-individual genetic variation on several infection parameters (Mackay et al., 2012). All wild-type stocks and isogenic lines were made axenic then re-associated with a controlled mixture of 5 bacterial species that compose the core microbiota of Drosophila (Wong et al., 2011) to ensure homogeneity of associated microbiological communities. Briefly, eggs of these stocks were sterilized by bleaching and reassociated with a mixture of 5 common gut microbes of Drosophila (Acetobacter pomorum, Acetobacter tropicalis, Lactobacillus plantarum, Lactobacillus brevis, Lactobacillus fructivorans) as previously published. Immune mutant stocks have been described previously and include mutants of the Imd pathway (PGRP-LE112, PGRP-LCΔE double mutants (Takehana et al., 2004), DreddEP1412 (RRID:BDSC_10456, Leulier et al., 2000) and RelishE20 mutants (RRID:BDSC_55714, Hedengren et al., 1999), mutants of the Toll pathway (spzrm7, [Lemaitre et al., 1996]), and mutants of the melanization cascade (double mutant PPO1Δ, PPO2Δ, [Binggeli et al., 2014]). We used the Gal80TS; daughterless-Gal4 driver in combination with UAS-Imd to ubiquitously induce the Imd pathway at the adult stage. To generate phagoless adult flies, we ablated phagocytes by inducing the pro-apoptotic gene Bax in hemocytes using the hemocyte-specific driver Hml-Gal4 (Hml-Gal4 >UAS GFP [RRID:BDSC_30140]; UAS-Bax, Gal80TS) in adult flies at 2 days post-eclosion (Defaye et al., 2009).
The bacterial strains used in this study include the Gram-negative Providencia rettgeri (strain Dmel) (Juneja and Lazzaro, 2009), P. burhodogranariea (Juneja and Lazzaro, 2009), Serratia marcescens (strain Db11) (Kurz et al., 2003), E. coli (Type strain), Pectinobacterium carotovorum carotovorum 15 (strain Ecc15-GFP; formerly genus Erwinia, [Basset et al., 2000]), Pseudomonas entomophila (Vodovar et al., 2005), and Providencia alcalifaciens (Juneja and Lazzaro, 2009), as well as the Gram-positive Enterococcus faecalis (isolated from wild-caught D. melanogaster by B. Lazzaro), and Staphylococcus aureus (strain PIG1 [Liu et al., 2005]). Cultures were grown to saturation overnight at 37°C (29°C for E. carotovora) in LB liquid medium (LB broth, Miller, VWR). Saturated cultures were suspended and diluted in Phosphate Buffered Saline (PBS, pH 7.4) to the desired optical density (OD600 = 1 or as otherwise indicated). We injected 23 nL of bacterial suspension (c.a. 30,000 bacteria for OD600 = 1) into each fly abdomen using a Nanoject II (Drummond) (Khalil et al., 2015). Flies were anesthetized with light CO2 for about five minutes during the injection procedure and were observed shortly after injection to confirm recovery from manipulations. We measured host survivorship post-injection in groups of 20 or 50 males kept with ad libitum access to food. Differences in survivorship were tested with Cox regression models (R package Survival).
To characterize the dynamics of within-host bacterial loads, after being dipped in Ethanol 70% and washed in PBS to limit contamination by external bacteria, individual flies were homogenized in 500 µl of sterile PBS with an HT homogenizer (OPS Diagnostics) at each timepoint post-injection. The homogenate was then diluted to 1:100 or 1:1000 in PBS to ensure that plate counts remained within the limits of resolution of the plating system. We plated 70 µl onto LB agar using a WASP II Autoplate spiral plater (Microbiology International). Plates were incubated overnight at 37°C (29°C for E. carotovora) and bacterial colonies were counted using an EZ-Count Automated Colony Counter (Microbiology International) to estimate the number of viable bacteria per fly. To estimate the Bacterial Load Upon Death (BLUD), infected hosts were checked every 30 min and newly dead flies (flies not moving and on their side or back) were collected and homogenized, with bacterial load quantified as described above.
Total RNA was extracted from either 20 flies or single flies with TRIzol reagent (Invitrogen). Template RNA (1 µg) was used to generate cDNA by reverse transcription using the SuperScript II cDNA synthesis kit (Promega) and then analyzed by quantitative polymerase chain reaction (qPCR) using the PerfeCTa SYBR Green SuperMix (Quantabio). Expression values were normalized to RpL32. Primer sequences are available in method supplementary 1. At least three independent repeats were done.
We used non-parametric tests to compare bacterial loads between groups of flies: Wilcoxon for two groups comparison and Kruskal-Wallis for more than two groups. The only exception was the comparison of loads between immune deficient mutants and wild-type controls, where, to take into account that we used more than one genotype to test the same biological hypothesis, we used a linear regression with genotype as a random effect. We used Spearman correlations to assess the relationship between bacterial loads and several quantitative traits (initial inoculum, time to death, and Dpt expression). All data analysis was performed using R version 3.3.1 (R Core Team , 2017).
Our results demonstrate that experimentally identical hosts exposed to the same pathogenic infections can either control the infection or succumb to acute bacterial proliferation. In addition, we see that bacterial growth in the absence of immune control has a similar rate to that of in vitro LB cultures (see Results). We therefore built a model that consists of a mixture of two demographic models, one describing within-host bacterial growth and one modeling bacterial elimination (see Figure 6A for a representation of the conceptual model this statistical model is based upon). More precisely, we considered that each bacterial load estimated on a single fly (n) was a log-normal random variable with the mean computed from a Baranyi model (Baranyi and Roberts, 1994) when bacteria grow in the host uncontrolled by the fly's immune system with
or from an exponential decrease model, when bacteria are controlled by the host, with
With increasing time, t, n begins at n0 and approaches nmax in the Baranyi model, while it decreases from nc to asymptotically reach zero in the exponential decrease model. Actual bacterial load was assumed to follow a log-normal distribution, with average given by either of the two previous equations and variance fixed at σ2b for the Baranyi model and σ2d for the exponential decrease.
To reflect the probabilistic nature of the infection process, we developed a model to compute the probability that bacteria are either controlled (and thus have an exponential decrease) or are not (and thus have a Baranyi growth dynamic). We considered that the time required for the host immune defense to efficiently control bacterial infection follows a Gamma distribution with a fixed average tlag + t̄c and a variance Vc. t̄c is the average time host defenses take to control bacteria after they have started to grow (Tc), and Vc quantifies inter-individual variation in Tc. We further assumed that the host cannot control infections after the bacterial load reaches a fixed threshold, which we termed the tipping point (ntip). The eleven parameters of this model are all described in Table 1.
Under the model, the probability that a host survives infection is the probability that the control of bacterial growth occurs before the tipping point is reached. More precisely, if Pc is the probability of survival, ttip is the time at which the tipping point is reached, and Tc is a random variable corresponding to the time of control observed for a given host, we can write
The probability of survival Pc therefore depends on parameters that determine how fast bacteria grow initially (n0, nmax, tlag, µ), on the tipping point itself (ntip), and on parameters, t̄c and Vc, that govern the distribution of time to control Tc.
All these parameters can be estimated by adjusting the model on experimental data consisting of bacterial loads measured over time in different individuals. This was performed by computing the likelihood of the model for each bacterial load x observed at time t as
where is the likelihood under the Baranyi model and the likelihood under the exponential decrease model. The global log-likelihood of the model is then obtained by summing log-likelihoods computed for each point of the datasets. This was possible because each point of our data set consists of a single individual host, and all measurements at all sampling times are therefore somewhat independent from each other (i.e. the measurements are not repeated on the same fly). We illustrated simulated outcomes of this model of within-host dynamics in Figure 6—figure supplement 1.
We estimated the 11 parameters of this mixture model by maximizing log-likelihood, using the Optim procedure in R. Confidence intervals for parameter estimates were then obtained by bootstrapping data for each experimental time point. This model, once adjusted, allows estimation of the probability that a host controls an infection by computing the probability that effective control occurs before ntip is reached. Probability of control therefore depends on the parameters of the Baranyi model (n0, nmax, tlag, µ) and of those describing the distribution of time to control (t̄c, Vc, ntip).
Virulence evolution and the trade-off hypothesis: history, current state of affairs and the futureJournal of Evolutionary Biology 22:245–259.https://doi.org/10.1111/j.1420-9101.2008.01658.x
Epidemiological and clinical consequences of within-host evolutionTrends in Microbiology 19:24–32.https://doi.org/10.1016/j.tim.2010.09.005
Multiple infections and the evolution of virulenceEcology Letters 16:556–567.https://doi.org/10.1111/ele.12076
A dynamic approach to predicting bacterial growth in foodInternational Journal of Food Microbiology 23:277–294.https://doi.org/10.1016/0168-1605(94)90157-0
Immunity in Drosophila melanogaster--from microbial recognition to whole-organism physiologyNature Reviews Immunology 14:796–810.https://doi.org/10.1038/nri3763
Thorax injury lowers resistance to infection in Drosophila melanogasterInfection and Immunity 82:4380–4389.https://doi.org/10.1128/IAI.02415-14
R: A Language and Environment for Statistical ComputingR Foundation for Statistical Computing, Vienna, Austria.
Models of infectious diseases in the fruit fly Drosophila melanogasterDisease Models and Mechanisms 1:43–49.https://doi.org/10.1242/dmm.000307
Pathogen responses to host immunity: the impact of time delays and memory on the evolution of virulenceProceedings of the Royal Society B: Biological Sciences 273:2083–2090.https://doi.org/10.1098/rspb.2006.3552
The evolutionary consequences of stepwise infection processesTrends in Ecology & Evolution 32:612–623.https://doi.org/10.1016/j.tree.2017.05.009
Providencia sneebia sp. nov. and Providencia burhodogranariea sp. nov., isolated from wild Drosophila melanogasterInternational Journal Of Systematic And Evolutionary Microbiology 59:1108–1111.https://doi.org/10.1099/ijs.0.000117-0
Systemic bacterial infection and immune defense phenotypes in Drosophila melanogasterJournal of Visualized Experiments, 10.3791/52613, 25992475.
Tipping elements in the human intestinal ecosystemNature Communications 5:4344.https://doi.org/10.1038/ncomms5344
The host defense of Drosophila melanogasterAnnual Review of Immunology 25:697–743.https://doi.org/10.1146/annurev.immunol.25.022106.141615
Staphylococcus aureus golden pigment impairs neutrophil killing and promotes virulence through its antioxidant activityThe Journal of Experimental Medicine 202:209–215.https://doi.org/10.1084/jem.20050846
Reactive oxygen species modulate Anopheles gambiae immunity against bacteria and PlasmodiumJournal of Biological Chemistry 283:3217–3223.https://doi.org/10.1074/jbc.M705873200
On the evolutionary ecology of specific immune defenceTrends in Ecology & Evolution 18:27–32.https://doi.org/10.1016/S0169-5347(02)00013-7
Disease tolerance and immunity in host protection against infectionNature Reviews Immunology 17:83–96.https://doi.org/10.1038/nri.2016.136
Low-diversity bacterial community in the gut of the fruitfly Drosophila melanogasterEnvironmental Microbiology 13:1889–1900.https://doi.org/10.1111/j.1462-2920.2011.02511.x
Bruno LemaîtreReviewing Editor; Ecole Polytechnique Fédérale de Lausanne, Switzerland
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Stochastic variation in the initial phase of bacterial infection predicts the probability of survival in D. melanogaster" for consideration by eLife. Your article has been favorably evaluated by Wendy Garrett (Senior Editor) and three reviewers, one of whom, Bruno Lemaître (Reviewer #1), is a member of our Board of Reviewing Editors. The following individual involved in review of your submission has agreed to reveal their identity: Andrea Graham (Reviewer #2).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
This elegant and thorough study provides rigorous quantitative dissection, under controlled environmental conditions, of the contributions of host and pathogen genetics plus random individual variation – in this case, stochasticity in the infection process itself – to the dynamics and outcome of infection. The authors make excellent use of the tools available for Drosophila, including isogenic lines and immune knockouts, as well as a wide array of bacterial species. They are therefore able to quantify variation among pathogen species in odds of killing the host, in the set point bacterial load, and in the consistent bacterial burden at which death occurs (regardless of the time at which that threshold burden is reached, interestingly). Likewise, their extensive empirical data show variation among host species and isogenic lines in bacterial load upon death, which appear independent of inoculum size and immune deficiency. The authors go on to show that stochastic variation within many of their host-pathogen pairings leads to bimodal outcomes: death or chronic infection. The authors explain this bifurcation as a race between rates of bacterial proliferation and immune response induction. The evolutionary implications of this are clear – e.g., selection on rates of immune response induction or bacterial proliferation. Indeed, evolutionary biologists, ecologists, microbiologists, physiologists, and immunologists working in all sorts of biological systems will be intrigued and inspired by this array of findings.
There are, however, a number of important issues to clarify before considering this paper for publication.
1) Impact of initial conditions:
A) In particular, the overall conclusion is that variations in the timing of initiating the humoral immune response (tc) determine the ultimate fate of infected animal – either death or survival with chronic infection. One trivial possibility that is not directly addressed is that (small) differences in initial inoculation (with little or no variation in the tc) are actually responsible. This would be consistent with data present in Figure 2—figure supplement 1, and the mathematical model presented, without postulating variation in tc. It seems that empirically determined initial dose is only reported in one figure panel (Figure 5D). Given the variation in the pulled capillary needle and the injection device, how can the authors be certain small variations don't tip the balance?
1B) While it looks like that all individual flies were treated identically, another parameter that is clearly different for each individual fly is the exposure to CO2 during manipulations/injection. It is probable that those flies that were injected first (had lower exposure to CO2) in the end were survivors, while those injected last (had longer exposure to CO2) succumbed. This should be experimentally addressed, considering that CO2 is known to increase Drosophila sensitivity to infections (Helenius et al., 2009).
Collectively, the authors should provide more information to support their statements that this is not small variation in the initial condition that explain the stochasticity. While the influence of CO2 might be easily tested, the impact of differences in the initial dose could be better characterized by adding more measurement of the initial dose for several bacteria. Another way could be to monitor stochasticity in infected flies with two related doses in the same setting (Dose 1 being always slightly superior of Dose 2) and see the existence of flies that die faster with Dose 2 compared to Dose 1).
2) Analysis of Immune response and bacterial growth:
The analysis of Imd signaling (Dpt induction) on single flies is a key part of this paper, yet is massively under-developed. Only one time point is shown, and given that transcript is measured, when mature AMP is the functional output, it is not at all clear that this one time point is relevant. Then, in the next panel, a different time point is examined to argue that Dipt levels do not correlate with bacterial loads. It can be argued that an earlier time point, such as 2 hours, is really the most relevant, as the model argues that the timing of the initiation of the response is deterministic, and perhaps at the initiating time point a correlation between load and response would be observed. Did the authors try to analyze immune response and bacterial growth in live flies using adequate fluorescent reporter (Dpt-mcherry, P. rettgeri GFP) ? While being less quantitative, this could bring another set of data that support their main conclusions.
On this line, the notion that the methods used here demonstrate bacterial growth in an individual fly (see language in Figure 3 legend) is grossly overstated. Clearly, as a population, the evidence clearly demonstrates bacterial growth in vivo. However, the bacterial growth within any individual fly cannot be assessed by these methods.
3) Clarity of the text and relevant information:
Many parameters affect infection such a temperature, sex, injection methods…. This information should be in the Results section not hidden in the Materials and methods. It is not clear what each survival graph shows, is it a representative result? How many times was it repeated? How many flies? Males or females? For bacterial load graphs, were all individual flies taken from one experiment (one biological replicate) or is it pooled from several replicates? It is often difficult to know how the experiments were done. All these details should be provided. Moreover, we could not exclude batch effect. This is actually suggested by Figure 4B with the new bacteria appearing more virulent.
4) All the data should be made publicly available to allow other scientists to use them. I would suggest to include all the raw data in excel sheets (bacterial count, survival data; biological repeats…).
5) The idea that the bacteria proliferate at the same rate in LB medium and flies should be documented for other strains.
6) Description of the model: more information on the model should be included to the main text. As it is written it is hard to understand this section of the paper. Moreover, even as written it is not fully and clearly explained to a non-expert in the field of modeling. For example, the meaning and import of mini graph of "Frequency vs. tc was unclear to one reviewer, and the relationship of Pc to the area below "the" curve was likewise unclear, as there are 2 curves presented. Some parameters are introduced but described later: Vc is shown in the figure but mentioned only in the Discussion. The confusion on the description of the model stem partly because of the high density at which series of hypotheses, tests and results are presented, and different sets of empirical results invoked (e.g., the third paragraph of the subsection “A mixture model to capture within-host bacterial growth dynamics”, which covers huge terrain very rapidly). Perhaps greater use of paragraph breaks or subheadings would help the reader navigate.
Globally, a longer paragraph with the model, the definition of each parameter could improve the impact of this paper.
[Editors' note: further revisions were requested prior to acceptance, as described below.]
Thank you for resubmitting your work entitled "Stochastic variation in the initial phase of bacterial infection predicts the probability of survival in D. melanogaster" for further consideration at eLife. Your revised article has been favorably evaluated by Wendy Garrett (Senior Editor), a Reviewing Editor, and one reviewer.
The manuscript has been improved but there are some remaining issues that need to be addressed before acceptance, as outlined below:
Both the external reviewer and the reviewing editor found that the authors have adequately answered to most of reviewers’ comments. As you can see from the reviewer 3, there is regret that the paper does not include in vivo monitoring of antimicrobial peptide gene expression to back up the model. Nevertheless, it is true that GFP reporters are slow and might not be the appropriate tool to monitor the immune dynamic at early time point. To address this in particular, we would recommend that you adjust your text, in both the Results and the Discussion, to make it explicitly clear that you are inferring a variance in the initiation of response, but do not (yet) have tools available to demonstrate it.
This remains one of the most thought provoking papers I have read in recent years. I think it will be a landmark in the field.
That being said, I continue to struggle with the big issue, which is not directly addressed by the revisions. In short, the authors have inferred that the differences in the timing or strength of the initial immune response following infection determine the outcome (death at BLUD or chronic infection at SPBL). Surprisingly, they write in the rebuttal that they have been unable to measure AMP induction via fluorescent protein reporters at these early time points, at least in intact flies. Ideally, I would love some type of data that can probe this inference.
Another more significant issue, is (still) the presentation of the mathematical model. Lots is made of the probability calculation in Figure 6 and beyond, but the way this probability is determine is not explained in the Results section. It is presented in detail in the Materials and methods, but frankly the math is over my head. I would prefer some discussion of this computation directly in the Results section so as to make sense of Figure 6B, for example.
Finally, the rationale for the time points used in 7B (4 hour) and 7C (8 hours) is presented, and argued compellingly, in the rebuttal but this logic is still missing from the actual text of the article. I would recommend including.https://doi.org/10.7554/eLife.28298.032
- Nicolas Buchon
- Brian P Lazzaro
- Nicolas Buchon
- David Duneau
- David Duneau
- David Duneau
- Nicolas Buchon
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank Patrícia Beldade, Gary Blissard, Jeffrey Hogdson, Marine Cambon, Kathleen Gordon, Vanika Gupta, and Joo Hyun Im for comments on the manuscript, Hidetoshi Inamine for thoughtful discussions. We thank Bruno Lemaitre and Luis Teixeira for fly stocks. This work was funded by NSF 1354421, NSF 1656118 and funding from the Empire State Stem Cell Fund through New York State Department of Health NYSTEM contract C029542 for NB, NIH grant R01 AI083932 to BPL and a fellowship from the Swiss National Foundation (P300P3_147874) and from the French Laboratory of Excellence project ‘TULIP’ (ANR-10-LABX-41; ANR-11-IDEX-0002–02) to DD.
- Bruno Lemaître, Reviewing Editor, Ecole Polytechnique Fédérale de Lausanne, Switzerland
© 2017, Duneau et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.