#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Spatial Structure, Transmission Modes and the Evolution of Viral Exploitation Strategies


Why do some viruses use latent infection strategies and in which environments can such prudent host exploitation strategies evolve? Theory of mathematical epidemiology and social evolution predicts that prudent pathogens can evolve when they cluster in space and share the mutual benefits of prudent exploitation of the host population. Here we extend the theory to study the effect of spatial structure on the evolution of pathogens that can transmit both horizontally and vertically. We explore these effects experimentally by competing the latent bacteriophage λ and the virulent mutant λcI857 in spatially structured epidemics and gradually erode spatial structure. We show that the latent bacteriophage λ only wins in a spatially structured epidemic. Yet, a breakdown of the epidemic structure by long-range transmission reduces the benefit of latency by 500 fold. This demonstrates that long-rage transmission, as for example by air travel, might select for much more virulent types of previously latent viruses.


Published in the journal: . PLoS Pathog 11(4): e32767. doi:10.1371/journal.ppat.1004810
Category: Research Article
doi: https://doi.org/10.1371/journal.ppat.1004810

Summary

Why do some viruses use latent infection strategies and in which environments can such prudent host exploitation strategies evolve? Theory of mathematical epidemiology and social evolution predicts that prudent pathogens can evolve when they cluster in space and share the mutual benefits of prudent exploitation of the host population. Here we extend the theory to study the effect of spatial structure on the evolution of pathogens that can transmit both horizontally and vertically. We explore these effects experimentally by competing the latent bacteriophage λ and the virulent mutant λcI857 in spatially structured epidemics and gradually erode spatial structure. We show that the latent bacteriophage λ only wins in a spatially structured epidemic. Yet, a breakdown of the epidemic structure by long-range transmission reduces the benefit of latency by 500 fold. This demonstrates that long-rage transmission, as for example by air travel, might select for much more virulent types of previously latent viruses.

Introduction

When individuals compete for a common resource natural selection often favours more aggressive exploitation strategies. This may lead to resource exhaustion and, consequently, to population extinction, a process known as the 'tragedy of the commons'[1]. Because pathogens compete for a common resource (the host population), the same process may select for extreme exploitation strategies. Yet, two main ecological factors can alter this evolutionary outcome and promote the evolution of less aggressive exploitation strategies.

First, epidemic spread reduces the density of susceptible hosts and can feed back on the selective pressure that favours intermediate exploitation strategies [25]. Indeed, decreased availability of susceptible hosts weakens selection for higher transmission rates [68]. Ultimately, this may select for intermediate evolutionarily stable strategies, balancing the benefit (transmission) and the cost (virulence: induced host mortality) of host exploitation. We recently tested this idea using the bacteriophage λ in experimental epidemics spreading through well-mixed environments [9]. We confirmed that more virulent strains are indeed selected for during early epidemics, when uninfected hosts are abundant, but also that natural selection favours latent strains of the virus as disease prevalence increases.

Second, theoretical studies on spatially structured populations indicate that localized transmission often favours more prudent host exploitation strategies [1021]. Indeed, when epidemics are spatially structured, extreme strategies of host exploitation may lead to the over-exploitation of the local host supply and therefore fail to invade. The intensity of local competition for susceptible hosts depends, however, on the precise genetic and epidemiological structure of the pathogen population [21]. The effect of spatial structure on the evolution of pathogens has had some experimental support [2223] which is, however, based on experiments with obligate killing and strictly horizontally transmitting pathogens. Previous experimental work has shown that vertical transmission selects for lower virulence [2425] but the effect of spatial structure on vertically transmitted pathogens remains an open question. In the present study we explore the effect of spatial structure on the evolution of pathogens that can be transmitted both vertically and horizontally. We first develop and analyse a theoretical model to understand the evolutionary dynamics of pathogens in spatially structured epidemics (Fig 1). This model allows infected hosts to reproduce and to transmit their pathogens vertically. We use this model to study the competition between different pathogen genotypes under various levels of mixing (see supporting information in S1 Text). Then, we confront these predictions with experiments performed with the horizontally and vertically transmitting temperate bacteriophage λ. Bacteriophage λ is an avirulent pathogen that controls its own replication by the virulence repressor cI and integrates into the host genome as a prophage without killing the host cell—a process known as lysogenization. The virulence repressor cI furthermore controls prophage λ reactivation and excludes superinfection by a second viral particle [26]. Yet, spontaneous mutations in cI can markedly increase virulence [27]. We studied competition between temperate bacteriophage λ and its virulent mutant λcI857. Mutant λcI857 carries an unstable virulence repressor protein (cI857) that biases the viral life cycle towards horizontal transmission: whereas λ lysogenizes about 60% of newly infected cells killing only about 40% of cell, λcI857 kills 96% of host cells upon first infection (see Fig S2B in S1 Text). In addition λcI857 shows a higher rate of prophage reactivation, sacrificing infected host cells for increased production of free viral particles (Fig S2A in S1 Text). Furthermore, λcI857 is less efficient in the exclusion of superinfection (see Fig S2B in S1 Text). Whereas the temperate λ wildtype shows an avirulent host exploitation strategy the virulent mutant λcI857 readily exploits and kills its host cell. To experimentally investigate the effect of spatial structure on the success of these opposing strategies we competed λ and λcI857 in experimental epidemics with different degrees of spatial structure (Manipulation of spatial structure is described in Material & Methods and [28]).

Fig. 1. Schematic representation of the pathogen life cycle.
Schematic representation of the pathogen life cycle.
In our model we assume that hosts can either be uninfected or infected by the wild type (in green) or a mutant genotype (in red). Transmission can take two different routes. Horizontal transmission occurs between a susceptible host and an infected host with genotype i with rate βi[I|S] (see supporting information Theory in S1 Text). Vertical transmission occurs when and infected host reproduces in an empty site with rate bI δq0/i, where bI is the reproduction rate of infected hosts, δ is the fidelity of vertical transmission and q0/i measures the average local density of empty sites around hosts infected with genotype i. Among hosts infected by the same genotype we distinguish hosts that got infected horizontally (noted IiH) or vertically (noted IiV). Both types of hosts can reproduce and have an increased mortality rate αi (i.e. virulence) due to the infection with genotype i.

Results

Theory

We study a spatial version of the compartmental model depicted in Fig 1. This allows us to make predictions on the epidemiological and evolutionary consequences of spatial structure. First, Fig 2A shows that lower levels of mixing slow down the spread of the infection because in a spatially structured environment, pathogens experience a lower density of susceptible hosts. Second, the change in frequency of a virulent mutant strain is (see Supporting Information Theory in S1 Text):

The first term measures the contribution of the horizontal transmission route and depends on the local density of susceptible hosts in the neighbourhood of a host infected by a mutant or wild-type parasite The second term measures the contribution of vertical transmission and depends on the local density of empty sites experienced by hosts infected by each type of parasite. Indeed, a vertical transmission event is conditional upon a reproduction event and thus on the availability of empty sites in the neighbourhood. Finally, the last term measures the cost of virulence.

Fig. 2. Effect of spatial structure on epidemiology and evolution in the model.
Effect of spatial structure on epidemiology and evolution in the model.
The top figures (A) show simulation results for the density of susceptible (dashed line) and infected hosts (full line) for different level of mixing (left: no mixing (gH = gP = 0), middle: intermediate (gH = gP = 0.5), right: full mixing (gH = gP = 1)). The figures in the middle (B) show simulation results for the change in mutant frequency in the total pathogen population (black), in the horizontally infected hosts (red) or in vertically infected hosts (blue) for different levels of mixing. The figures at the bottom (C) indicate the values of the different components of the selection coefficient identified in Eq (1): the horizontal transmission term (red), the vertical transmission component (blue) and the sum of the previous two terms (black) and the magnitude of the cost of virulence (dashed black). Note that the global mutant frequency increases only when the sum of the transmission terms (full black line) is higher than the cost of virulence (dashed line). Our simulations are the result of the numerical integration of the deterministic approximation of the full spatial model, using Improved Pair Approximation (IPA; [31]) with parameters ϕ = 1/6 and θ = 2/5, which correspond to a triangular lattice (each site has 6 neighbours, see [31] for details). Parameters: αw = 0.1, αm = 1, βw = 2, βm = 3.5, δ = 0.99, d = 0.01. Initial conditions: S(0) = 0.8, IwH(0)=ImH(0)=0.01,IwV(0)=ImV(0)=10−5.

Fig 2B shows numerical simulations based on the pair approximation [2931] under different levels of mixing (see supporting information Theory in S1 Text). Note that, under all level of mixing, the frequency of the virulent mutant initially increases. Indeed, early on during the epidemics, parasites have access to a large density of susceptible hosts, which favours the spread of more transmissible mutants. This benefit, however, is only transient. As the pool of susceptible hosts is exhausted, more prudent pathogens end up being selected under all mixing treatments. Spatial structure, however, has a dramatic impact on the virulent mutant: less mixing drastically reduces the transient benefit of virulence. To understand this effect, we represent on Fig 2C the effect of spatial structure on the different components of the selection coefficient in Eq (1). Spatial structure acts through two different effects: (i) space reduces the horizontal transmission benefit (red curve in Fig 2C) and (ii) space generates benefit for the virulence mutant through vertical transmission (blue curve in Fig 2C). The first effect is the classical effect of space often discussed as a self-shading or kin selection argument [14,21,32]. The second effect is due to the accumulation of empty sites near virulence mutants, which is caused by the death of hosts infected by related (and thus virulent) genotypes. In other words, virulence mutant have more opportunities to transmit vertically than avirulent genotypes.

In addition, to better grasp the effect of space on pathogen evolution, we track the change in mutant frequency over two types of infected hosts: hosts that acquired the infection horizontally and hosts that acquired the infection vertically. Fig 2B reveals the effect of mixing on the change in frequency of the mutant in these two compartments. In a well-mixed environment the frequency of the mutant is always higher in the horizontally infected hosts. In contrast, in a structured environment, the frequency of the virulent mutant can be higher in vertically infected hosts. This is due to the effect of spatial structure on the opportunities for vertical transmission. As explained above the virulent mutant has more opportunities for vertical transmission because the local density of empty sites is higher (see Fig 2B). Another way to see this effect is to track the ratio WH / WV which provides a measure of the relative contribution of the two routes of transmission to mutant fitness. When this ratio is above 1 the change in mutant frequency is mainly governed by the horizontal transmission route. In contrast, when it is below one the vertical transmission route contributes more to the mutant fitness. Fig 3 shows that during the early stage of the epidemics horizontal transmission drives the evolution of the virulent mutant under all mixing treatments but that the contribution of vertical transmission is always increased by spatial structure.

Fig. 3. Effect of mixing on the relative contribution of horizontal (WH) and vertical transmission (WV) to the fitness of the virulent mutant in the model.
Effect of mixing on the relative contribution of horizontal (W<sub>H</sub>) and vertical transmission (W<sub>V</sub>) to the fitness of the virulent mutant in the model.
Each curve corresponds to different levels of mixing: no mixing gH = gP = 0 (dotted line), intermediate mixing gH = gP = 0.5 (dashed line) and full mixing with gH = gP = 1 (full line). The contribution of horizontal transmission is maximal in the early stage of the epidemic because many susceptible hosts are still available. More mixing increases the contribution of horizontal transmission to the fitness of the mutant throughout the epidemics.

Experiment

The above model is an attempt to capture the interplay between the effects of spatial structure and vertical transmission on the evolution of pathogen virulence. To test these predictions experimentally we followed competition of the avirulent λ and virulent λcI857 in spatial epidemics after an initial inoculation of a 1:1 ratio of cells infected by λ or λcI857 onto a biofilm of susceptible cells. Thereupon spatial structure was either left undisturbed or disturbed for 30s or 24h in the absence of liquid, or 24h in the presence of liquid (24h-wet). By engineering fluorescent protein expression cassettes (CFP and YFP) into the viral genome we tagged infected cells by a fluorescent colour. This enabled us to directly visualize pathogen-host structure and spatial prevalence in these four environments (see Fig 4). In the undisturbed treatment the epidemic spreads in a circular front segregating by viral types (λCFP or λYFP) but large areas remain uninfected (see Fig 4A). In the 30s disturbance treatment the epidemic spreads through the entire biofilm and creates epidemic clusters dominated by a single type of virus (Fig 4B). After 24h disturbance epidemic clusters disappear, but the biofilm remains structured (Fig 4C). In the 24h-wet treatment epidemic structure is completely homogenized (Fig 4D). Thus, as expected from our model simulation (Fig 2A), higher mixing is associated with a significant increase in overall prevalence of infection (χ12 = 36.52, p<0.0001, Fig S3 in S1 Text).

Fig. 4. Competition of fluorescently marked virus (λCFP and λYFP—blue and green areas) during invasion into uninfected cells (black areas) at different degrees of mixing (undisturbed, 30s, 24h, 24h-wet).
Competition of fluorescently marked virus (λCFP and λYFP—blue and green areas) during invasion into uninfected cells (black areas) at different degrees of mixing (<i>undisturbed</i>, <i>30s</i>, <i>24h</i>, <i>24h-wet</i>).
(A) Undisturbed environment—Circular epidemic front segregates into single-virus sectors (blue and green sectors, white arrow depicts the centre of epidemic inoculation). (B) 30s disturbance—Epidemic clusters that are dominated by a single type of virus appear. (C,D) Longer periods of disturbance (24h) and a liquid surface layer (24h-wet) progressively homogenize the spatial genetic structure of virus types.

The mixing treatment strongly affects competitive fitness W of λcI857 (χ12 = 79.15, p<0.0001, Fig 5). In the undisturbed environment λcI857 fitness decreases by ~100 fold after the first transfer day, in contrast in the 24h-wet environment λcI857 fitness increases by ~5 fold. As expected, intermediate spatial disturbance treatments (30s and 24h) yield intermediate levels of λcI857 fitness (Fig 5). This demonstrates that erosion of the spatial structure increases the fitness of the virulent λcI857 by ~500 fold. Note that additional transfers did not affect this pattern except for the treatment 30s where λcI857 fitness increases with time (t = 5.08, p<10–4).

Fig. 5. Mixing selects for virulence.
Mixing selects for virulence.
A 1:1 mixture of temperate λ and its virulent mutant λcI857 was used to seed spatial epidemics. Spatial structure was disturbed to different degree (undisturbed, 30s, 24h, 24h-wet) and transferred onto a fresh biofilm of susceptible hosts for 3 consecutive days. Disturbance of the epidemic structure has a strong effect on the competitive fitness of the virulent λcI857. The undisturbed environment strongly selects against virulence whereas in a well-mixed environment (24h-wet) virulence is beneficial. As expected, intermediate levels of mixing (30s, 24h) lead to intermediate fitness of λcI857. Symbols represent 6 replications for each mixing treatment (with 3 independent replicates for each marker-mutant combination).

To disentangle the contribution of horizontal and vertical transmission to the fitness of λcI857 we repeated the experiment for the first round of spatial transfer with uninfected cells marked by the red fluorescent protein mCherry. This way we identified new infections from horizontal transmission as doubly coloured cells (mCherry and colour of the new infecting virus) and vertical transmission from the initial lysogen as singly coloured cells (colour of the originally integrated virus). From these data we calculated the competitive fitness of λcI857 among all new infections (horizontal fitness WH) and fitness among all vertical transmission events (vertical fitness WV) in each spatial environment. Fig 6 presents the effect of mixing on the ratio WH / WV between horizontal and vertical fitness components of λcI857. In the absence of mixing both routes of transmission contribute to λcI857 fitness. Interestingly, in agreement with our theoretical prediction (Fig 2C), we found that higher levels of mixing promote horizontal transmission (χ12 = 14,82, p = 0.002). For instance, in the highly mixed environment (24h-wet treatment), 3 out of 4 new infections result from the horizontal transmission route (Fig 6).

Fig. 6. Effect of mixing on the relative contribution of horizontal (WH) and vertical transmission (WV) to the fitness of the virulent λcI857.
Effect of mixing on the relative contribution of horizontal (W<sub>H</sub>) and vertical transmission (W<sub>V</sub>) to the fitness of the virulent λcI857.
More mixing (from undisturbed, 30s, 24h to 24h-wet) increases the fitness contribution of horizontal transmission to the competitive ability of λcI857 (WH) relative to the fitness contribution of vertical transmission (WV) (see Fig 3). Each box represents the median, the first and the third quartiles. Dashed lines delimit 1.5 times the inter-quartile range, above which individual counts are considered outliers and marked as empty circles.

Discussion

During an epidemic, disease prevalence increases, and the opportunities to find uninfected hosts dwindle. This may select for prudent exploitation of the remaining host population. In a previous study, we demonstrated the impact of such epidemiological feedback on the evolution of bacteriophage λ [9]. Here, we use the same experimental system to study the effect of spatial structure on the evolution of bacteriophage λ. To account for latency and vertical transmission in the life cycle of bacteriophage λ, we extend previous theory on the effect of spatial structure on pathogen evolution by allowing for vertical transmission of the pathogen. Our model tries to capture all the details of the bacteriophage λ life-cycle, but for tractability we do not explicitly model the free living stage of the virus. Our central result is that the fitness of a virulent strain relative to a more prudent strain decreases as the environment becomes increasingly structured. A general insight of previous theory is that the outcome of selection on parasite life-history traits depends on the balance between the genetic and epidemiological structures of the parasite population [21]. At the genetic level, parasites infecting a cluster of hosts tend to be related when transmission occurs locally. Hence, competition for the local host supply tends to take place among related pathogens (kin competition) [18,21,32]. At the epidemiological level, local pathogen spread reduces the local density of uninfected hosts, thereby degrading the parasite's local environment and diminishing future opportunities for horizontal transmission [5,14]. From models without vertical transmission, it has been shown that more prudent strategies of host exploitation are selected for if the genetic structure of the parasite population is above a threshold determined by the epidemiological structure of the host population [33]. As interactions become more localized, genetic structure increases and the balance is tilted at the advantage of more prudent strains [21]. We show that vertical transmission does not alter the general theoretical expectation that spatial structure should select for more prudent host exploitation strategies. Yet our analysis reveals that in a spatially structured environment vertical transmission can affect evolutionary dynamics via two main effects. First, vertical transmission affects opportunities of horizontal transmission through a modification of the local density of susceptible hosts (the first term in Eq (1)). Second, vertical transmission generates a force (the second term in Eq (1)) that may counteract the effect of spatial structure. During the early stage of the epidemic this effect remains very low but further work is required to investigate the impact of the effects of vertical transmission on evolutionary stable virulence strategies. Our model is useful to decouple the change in mutant frequency in horizontally and vertically infected hosts. This distinction between different host types is artificial but it helps to grasp the interplay between spatial structure and transmission routes. As in other situations where the pathogen can appear in different states [9] or in different hosts [8, 3436], tracking the change in frequency between different compartments enables us to generate a prediction regarding the effect of spatial structure on the relative contribution of the two transmission routes on pathogen fitness (Fig 2C).

This new theoretical model generates predictions that we tested experimentally with bacteriophage λ. Unlike our previous experimental study in chemostats [9], the present experiment has been realized on agar plates to manipulate the amount of spatial structuring. Yet, on a normal agar plate bacterial growth eventually enters stationary phase and the epidemiological and evolutionary dynamics halt. Consequently, our experiments allow us to test predictions during the initial phase of the epidemics only. During this phase we made three qualitative predictions that we verified experimentally. First, we confirmed that lower levels of mixing reduce the speed of the epidemic. Second, we confirmed that a virulent strain is outcompeted by a more prudent strain when the environment becomes increasingly structured. This is our main result and it agrees well with previous experimental studies realized with horizontally transmitted pathogens [2223]. Third, our experimental setup allows monitoring the relative contribution of horizontal and vertical transmission among new infections and confirms that mixing decreases the contribution of vertical transmission to pathogen fitness.

The evolution of host exploitation strategies can often be explained at the level of within-host selection. For parasites competing for within-host resources, more prudent host exploitation should be selected for when the probability of multiple infections decreases [4]. In our model, superinfection exclusion by a latent λ prophage guarantees that only related pathogen individuals compete within the same host cell and, hence, limits selection for increased replication and virulence at the within-host level. Similarly, many covertly infecting and genome-integrating viral pathogens take within- host competition out of the equation by active exclusion of superinfection. In herpes and retroviruses, for example, superinfection exclusion is well described and might be a key requirement for the maintenance of viral latency [3739]. Yet some mutations are known to enable λ to break through superinfection exclusion [26]. We did not observe such mutants in our present experiments, but we were able to extend our theory to explore the effect of spatial structure on the epidemiology and evolution of such virulence mutants that can superinfect already infected bacteria (see supporting information Theory in S1 Text). This analysis shows that mixing does affect the epidemiology and evolution in this alternative scenario: more mixing enhances the spread of the virulence mutant. Hence, this analysis suggests that the effect of spatial structure on the transient dynamics we report in the present study is robust to this modification of the life cycle. We believe that similar predictions may be feasible to understand and predict the epidemiology and evolution of a broad range of pathogens with complex transmission modes.

Methods

Theory

The theoretical analysis is presented in the supporting information (S1 Text). We examine the evolutionary dynamics of a vertically and horizontally transmitted pathogen in a spatially structured environment. The analytical model reveals that the change in frequency can be decomposed in different terms that characterize the local environment of competing variants of the pathogens. These terms can be tracked in numerical simulations (Fig 2). Further understanding of the evolutionary dynamics can be obtained by focusing on the change in frequency of a mutant pathogen in hosts that have acquired the pathogen horizontally (IH) and in hosts that have acquired the pathogen vertically (IV).

Experiments

Host exploitation strategies of λ and λcI857 were quantified by three independent assays. (1) Virus production (PFU/mL): Lysogens were grown to OD600nm = 0.6 at 30°C and shifted to 35°C for 2h until lysis occurred. Viral titers were determined by qPCR (Roche LightCycler480 primers F:5’AATGAAGGCAGGAAGTA3’ R:5’GCTTTCCATTCCATCGG3’) on filtered lysates (Millipore 96-well filtration plate Multiscreen HTS). CP values were calibrated to PFU/ml by top-agar plating a dilution series of a λvir (3x109 pfu). (2) Growth of infected hosts (CFU/mL): Lysogens were diluted to OD600nm = 0.07 and grown for 6h at 35°C in eight replicates in 96-well plates (900 rpm, Titramax shaker (Heidolph, Germany)). Final OD600nm was measured in an Infinity200 microplate reader (Tecan, Austria). OD600nm values were calibrated to CFU/ml by plating a dilution curve. (3) Lysogenization rate was determined by challenging non-infected E.coli MG1655 of OD600nm = 0.1 with 108 PFU/mL free virus particles of λCFP, λYFP, λcI857CFP and λcI857YFP for 24h. After 24h, the proportion of lysogenized (fluorescent) cells was determined by flow cytometry (Fig S2 in S1 Text).

Fluorescently labelled phage were constructed as described in [9]. Spatial epidemics were started by spreading a 108 cfu/ml suspension of uninfected E.coli MG1655 (starved 24h in 0.2%Maltose, 10mM MgSO4) onto agar plates until dried. Thereafter a 1:1 mixture of cells infected by λ and λcI857 was inoculated in localized spots by an array of 0.4mm stainless steel tattooing needles at about 4mm distance. Epidemic structure was disturbed by rolling 4mm sterile stainless steel beads over the dry agar surface for 30s or 24h or over agar surface humidified by sterile saline solution (24h-wet) or kept undisturbed. Agar plates were incubated at 35°C in agitation (250rpm), but only in the 24h and 24h-wet treatment steel beads remained on the agar surface. To control for marker effects all competitions were carried out with reciprocal marker-mutant combinations (λYFP vs. λcI857CFP and λCFP vs. λcI857YFP). Spatial structure was photographed on an Olympus BH-2 RCF microscope with a 2x fluorite objective and filter cube #69380 (Chroma Thechnology Corp., VT, USA). For the daily transfer of the epidemic structure the needle array described above was poked into plates that had been incubated for 24 hours and subsequently the needle array was poked into agar plates containing a layer of uninfected cells as described above. For flowcytometry the bacterial layer was washed off with 2 ml saline and analysed on a B&D Fortessa flowcytometer (Excitation CFP, YFP, mCherry at 405, 488 and 561 nm).

Horizontal and vertical transmission was estimated by repeating the spatial epidemics for a single cycle with lysogens of λCFP, λYFP, λcI857CFP and λcI857YFP carrying pBAD18 and uninfected E.coli MG1655 carrying plasmid pRSET-mCherry. 20μg/ml Ampicillin was added to the medium to prevent spontaneous plasmid loss. Horizontal transmission events were scored by double coloration (mCherry+CFP and mCherry+YFP). Horizontal and vertical fitness were calculated as mentioned above by only considering mCherry+ or mCherry- cells.

The relative fitness of λcI857 was calculated as W = (ft /(1- ft))/ (f0 /(1- f0)) where ft measures the frequency at time t. In general the frequency was calculated over the whole host population. In Fig 6 we distinguish horizontal fitness, WH, and vertical fitness, WV, by focusing either on the frequency in newly infected bacteria (bacteria carrying the plasmid pRSET-mCherry) or on the frequency in previously infected bacteria (bacteria that do not carry the plasmid pRSET-mCherry), respectively. In Fig 4 we present the ratio WH / WV for different mixing treatments.

Statistical analysis

Analyses were carried out using the R statistical package (version 3.0.2). The fitness W of λcI857 (Fig 5, S1 Data) was log transformed prior to fitting in a mixed effects model (lme, nlme package). Mixing (undisturbed, 30 s, 24h and 24h-wet) and transfer day (1, 2 and 3) were fitted as fixed effects while population was nested within the marker effect (the two reciprocal marker-mutant combinations) to account for the repeated measurements at different points in time. The ratio WH / WV (Fig 6, S1 Data) was fitted in a mixed effect model with mixing as a fixed factor and marker as a random effect. The prevalence of the infection (Fig S3 in S1 Text and S1 Data) was fitted in a mixed effect model with mixing as a fixed factor and marker colour as a random effect. Maximal models, including all higher-order interactions, were simplified by sequentially eliminating non-significant terms and interactions to establish a minimal model. The significance of explanatory variables in mixed effect models was established using a likelihood ratio test which is approximately distributed as a chi-square distribution [40]. A posteriori contrasts were carried out at the final time point of the experiments in Figs 5, 6 and S3 (in S1 Text) by aggregating factor levels together and by testing the fit of the simplified model using a likelihood ratio test [41].

Supporting Information

Attachment 1

Attachment 2


Zdroje

1. Hardin G. (1968). The Tragedy of the Commons. The population problem has no technical solution; It Requires a Fundamental Extension in Morality. Science 162(3859): 1243–1248. 17756331

2. Gandon S. and Day T. (2009). Evolutionary Epidemiology and the Dynamics of Adaptation. Evolution 63(4): 826–838. doi: 10.1111/j.1558-5646.2009.00609.x 19154387

3. Anderson R.M., and May R.M. (1982). Coevolution of hosts and parasites. Parasitology 85: 411–426. 6755367

4. Frank S. A. (1996). Models of parasite virulence. The Q. Rev. Biol. 71(1): 37–78. 8919665

5. Dieckmann U., Metz J.A.J., Sabelis M.W., and Sigmund K. (2002). Adaptive Dynamics of Infectious Diseases: In Pursuit of Virulence Management. Cambridge University Press, Cambridge MA.

6. Lenski R.E., and May R.M. (1994). The evolution of virulence in parasites and pathogens: Reconciliation between two competing hypotheses. J. Theor. Biol. 169(3): 253–265. 7967617

7. Day T., and Proulx S.R. (2004). A general theory for the evolutionary dynamics of virulence. The American Naturalist 163(4): E40–63. 15122509

8. Day T., and Gandon S. (2007). Applying population-genetic models in theoretical evolutionary epidemiology. Ecol. Letters 10(10): 876–888. 17845288

9. Berngruber T.W., Froissart R., Choisy M., and Gandon S (2013). Evolution of virulence in emerging epidemics. PLoS Pathogens 9(3): e1003209. doi: 10.1371/journal.ppat.1003209 23516359

10. Boerlijst M.C., Lamers M.E., and Hogeweg P. (1993). Evolutionary consequences of spiral waves in a host parasitoid system. Proc. R. Soc. B. 253: 15–18.

11. Sato K., Matsuda H., and Sasaki A.(1994). Pathogen invasion and host extinction in lattice structured populations. J. Math. Biol. 32(3): 251–268. 8182357

12. van Baalen M., and Sabelis M.W.(1995). The milker-killer dilemma in spatially structured predator-prey interactions. Oikos 74(3): 391–400.

13. Rand D. A., Keeling M., and Wilson H.B. (1995). Invasion, stability and evolution to criticality in spatially extended, artificial host-pathogen ecologies. Proc. R. Soc. B.: 259(1354): 55–63.

14. Boots M, and Sasaki A. (1999). Small worlds’ and the evolution of virulence: Infection occurs locally and at a distance. Proc. R. Soc. B.: 266(1432): 1933–1938. 10584335

15. Keeling M. J. (1999). The effects of local spatial structure on epidemiological invasions. Proc. R. Soc. B. 266(1421): 859–867. 10343409

16. Haraguchi Y., and Sasaki A. (2000). The evolution of parasite virulence and transmission rate in a spatially structured population. J. Theor. Biol. 203: 85–96. 10704294

17. Johnson C.R., and Boerlijst M.C. (2002). Selection at the level of the community: The importance of spatial structure. TREE 17(2): 83–90.

18. van Baalen M. (2002). Contact networks and the evolution of virulence, in Dieckmann U., Metz J. A. J., Sabelis M. W. and Sigmund K. (eds) Adaptive Dynamics of Infectious Diseases. Cambridge University Press, Cambridge.

19. Read J. M., and Keeling M. J. (2003). Disease evolution on networks: The role of contact structure. Proc. R. Soc. B.: 270(1516): 699–708. 12713743

20. Lion S., and van Baalen M. (2008). Self-structuring in spatial evolutionary ecology. Ecol. Letters 11(3): 277–295. 18070102

21. Lion S., and Boots M. (2010). Are parasites prudent in space?” Ecol. Letters 13(10): 1245–1255. doi: 10.1111/j.1461-0248.2010.01516.x 20727004

22. Kerr B., Neuhauser C., Bohannan B. J. M., and Dean A. M. (2006). Local migration promotes competitive restraint in a host-pathogen ‘Tragedy of the Commons’. Nature 442(7098): 75–78. 16823452

23. Boots M., and Mealor M. (2007). Local interactions select for lower pathogen infectivity. Science 315(5816): 1284–1286. 17332415

24. Messenger S.L., Molineux I.J., and Bull J.J. (1999). Virulence evolution in a virus obeys a trade-off. Proc Biol Sci 266, 397–404. 10097397

25. Pagán I., Montes N., Milgroom M.G., and García-Arenal F. (2014). Vertical transmission selects for reduced virulence in a plant virus and for increased resistance in the host. PLoS Pathog. 10, e1004293. doi: 10.1371/journal.ppat.1004293 25077948

26. Ptashne M. (1992) A Genetic Switch: Phage Lambda and Higher Organisms. Oxford: Blackwell Publishers.

27. Refardt D, Rainey PB Tuning a genetic switch: Experimental evolution and natural variation of prophage induction. Evolution (2010) 64:1086–1097 doi: 10.1111/j.1558-5646.2009.00882.x 19891623

28. Berngruber T.W., Lion S and Gandon S (2013). Evolution of suicide as a defence strategy against pathogens in a spatially structured environment. Ecol. Lett. 16(4):446–53. doi: 10.1111/ele.12064 23331662

29. Matsuda HN, Ogita A, Sasaki A & Sato K (1992). Statistical mechanics of population: The lattice Lotka-Volterra model. Progress in Theoretical Physics 88:1035–1049.

30. Rand DA (1999). Correlation equations and pair approximations for spatial ecologies. In: Advanced Ecological Theory: Advances in Principles and Applications, ed. McGlade JM., p. 100–142. Blackwell, Oxford UK.

31. Van Baalen M (2000). Pair Approximations for Different Spatial Geometries. In: The Geometry of Ecological Interactions: Simplifying Spatial Complexity, eds. Dieckmann U, Law R & Metz JAJ, p. 359–387. Cambridge University Press.

32. Wild G., Gardner A. and West S.A. (2009) Adaptation and the evolution of parasite virulence in a connected world. Nature 459: 983–986. doi: 10.1038/nature08071 19474791

33. Lion S., and Gandon S. (2015). Evolution of spatially structured host—parasite interactions. J. Evol. Biol. 28(1): 10–28. doi: 10.1111/jeb.12551 25439133

34. Day T. & Gandon S. (2006) Insights from Price’s equation into evolutionary epidemiology. In: Disease evolution: models, concepts and data analyses. eds. Feng Z. Dieckmann U.; Levin S. American Mathematical Society, p. 23–43.

35. Gandon S., and Day T. (2007). The evolutionary epidemiology of vaccination. J R Soc Interface 4, 803–817. 17264055

36. Lélu M., Langlais M., Poulle M.-L., Gilot-Fromont E., and Gandon S. (2013). When should a trophically and vertically transmitted parasite manipulate its intermediate host? The case of Toxoplasma gondii. Proceedings of the Royal Society of London B: Biological Sciences 280, 20131143. doi: 10.1098/rspb.2013.1143 23825211

37. Mador N., Panet A., and Steiner I.. 2002. The latency-associated gene of herpes simplex virus type 1 (HSV-1) interferes with superinfection by HSV-1. J. Neurovirol. 8:97–102. 12491159

38. Nethe M., Berkhout B., and van der Kuyl A.. 2005. Retroviral superinfection resistance. Retrovirology 2:52. 16107223

39. Berngruber T.W., Weissing F.J., and Gandon S. (2010). Inhibition of superinfection and the evolution of viral latency. J. Virol. 84(19): 10200–10208. doi: 10.1128/JVI.00865-10 20660193

40. Bolker B. (2008). Ecological Models and Data in R (Princeton University Press).

41. Crawley M.J. (2007). The R Book, John Wiley & Sons, Hoboken, NJ.

Štítky
Hygiena a epidemiologie Infekční lékařství Laboratoř

Článek vyšel v časopise

PLOS Pathogens


2015 Číslo 4
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

plice
INSIGHTS from European Respiratory Congress
nový kurz

Současné pohledy na riziko v parodontologii
Autoři: MUDr. Ladislav Korábek, CSc., MBA

Svět praktické medicíny 3/2024 (znalostní test z časopisu)

Kardiologické projevy hypereozinofilií
Autoři: prof. MUDr. Petr Němec, Ph.D.

Střevní příprava před kolonoskopií
Autoři: MUDr. Klára Kmochová, Ph.D.

Všechny kurzy
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#