The role of spatial structure in the evolution of viral innate immunity evasion: A diffusion-reaction cellular automaton model

Most viruses have evolved strategies for preventing interferon (IFN) secretion and evading innate immunity. Recent work has shown that viral shutdown of IFN secretion can be viewed as a social trait, since the ability of a given virus to evade IFN-mediated immunity depends on the phenotype of neighbor viruses. Following this idea, we investigate the role of spatial structure in the evolution of innate immunity evasion. For this, we model IFN signaling and viral spread using a spatially explicit approximation that combines a diffusion-reaction model and cellular automaton. Our results indicate that the benefits of preventing IFN secretion for a virus are strongly determined by spatial structure through paracrine IFN signaling. Therefore, innate immunity evasion can evolve as a cooperative or even altruistic trait based on indirect fitness effects that IFN shutdown exerts on other members of the viral population. We identify key factors determining whether evasion from IFN-mediated immunity should evolve, such as population bottlenecks occurring during viral transmission, the relative speed of cellular infection and IFN secretion, and the diffusion properties of the medium.

Most viruses have evolved strategies for preventing interferon (IFN) secretion and evading innate immunity. Recent work has shown that viral shutdown of IFN secretion can be viewed as a social trait, since the ability of a given virus to evade IFN-mediated immunity depends on the phenotype of neighbor viruses. Following this idea, we investigate the role of spatial structure in the evolution of innate immunity evasion. For this, we model IFN signaling and viral spread using a spatially explicit approximation that combines a diffusion-reaction model and cellular automaton. Our results indicate that the benefits of preventing IFN secretion for a virus are strongly determined by spatial structure through paracrine IFN signaling. Therefore, innate immunity evasion can evolve as a cooperative or even altruistic trait based on indirect fitness effects that IFN shutdown exerts on other members of the viral population. We identify key factors determining whether evasion from IFN-mediated immunity should evolve, such as population bottlenecks occurring during viral transmission, the relative speed of cellular infection and IFN secretion, and the diffusion properties of the medium.

Author summary
Innate immunity signaled by interferons is an important defense system against pathogens. Thus, it is often assumed that most pathogens can rapidly evolve strategies to avoid immune recognition or block interferon signaling. However, the presence of some interferon-stimulating pathogen variants in a given host may be sufficient for triggering innate immune responses, negatively impacting the fitness of interferon-blocking pathogens and posing a major hurdle for the evolution of innate immunity evasion. Here, we simulate the spread of viral infections and interferon-mediated innate immune responses in cell populations. We find that, by allowing for the insulation of interferon-stimulating and interferon-blocking virus variants, spatial structure should play a central role in the evolution of innate immunity evasion. In turn, we identify features determining spatial structure, such as the number of founder viruses, infection dynamics, medium viscosity, and virion size. Our results shed light on the evolution of a key process involved in pathogenesis and virulence.

Introduction
Innate immune signaling during early infection constitutes an important line of defense against viruses. Infected cells recognize pathogen-associated molecular patterns and secrete type-I interferons (IFNs), which act in autocrine and paracrine manners to halt the spread of the infection locally [1,2]. Autocrine signaling triggers antiviral responses in the secretor cell including gene expression arrest and apoptosis, whereas paracrine signaling induces a virusresistant state in neighbor cells. Innate immunity imposes a strong selective pressure on viruses, which have consequently evolved a variety of evasion mechanisms at the level of viral sensing, signal transduction, and/or gene expression [2][3][4]. If successful, these mechanisms lead to avoidance or shutdown of IFN secretion.
We have recently shown that a mutant of vesicular stomatitis virus (VSV) that fails to block IFN secretion reduces the fitness of neighbor, IFN-blocking VSV by triggering local antiviral responses [5]. This interference defines IFN shutdown as a social trait, since the fitness of a given virus infecting a given cell (focal virus) depends on other members of the viral population (neighborhood). A neighborhood constituted by viruses that block or do not stimulate IFN provides a fitness advantage relative to being located near IFN-stimulating viruses. This fitness effect is indirect because it is determined by cells not infected by the focal virus, is exerted through paracrine signaling, and depends on the spatial distribution of infected cells and immunized cells.
Blocking IFN should also have a direct effect on the focal virus determined by autocrine signaling. In principle, since autocrine signaling is aimed at reducing viral progeny production, IFN blockade should directly benefit the focal virus. However, if the viral infection cycle proceeds faster than autocrine antiviral responses, this direct benefit may not be realized. In addition, directing viral proteins to prevent IFN secretion (by avoiding recognition or blocking production) might have a direct cost for the focal virus. For instance, in VSV, the matrix protein M prevents IFN secretion by inhibiting host gene expression [6], but M is also an essential structural component of the virion. Hence, using M proteins for blocking gene expression in the nucleus might be costly for virion morphogenesis. Another potential source of direct costs to IFN blockade is that anti-IFN viral proteins could downregulate cellular gene expression or trigger premature apoptosis, reducing the availability of cellular resources for viral progeny production. The existence of such costs is suggested by our previous results with VSV [5].
Therefore, IFN shutdown should have direct and indirect effects on viral fitness. Because indirect effects depend on other members of the viral and cellular populations, tackling the evolution of innate immunity evasion requires considering the spatial structure and the dynamics of viral spread and IFN-mediated responses. In general, the term spatial structure refers to a non-random arrangement of individuals in space, and can originate from physical barriers, limited dispersal, a tendency to aggregate, demographical history, and other processes related to population dynamics [7]. For instance, in viral infections, host-to-host transmission and intra-host dissemination is limited by anatomical barriers. Yet, even in the absence of such barriers, most viruses infecting solid tissues exhibit spatial structure, as evidenced by the formation of infection foci. Such foci can result from a simple diffusion-reaction process involving adsorption of viral particles to cells, production of viral progeny, and diffusion of progeny particles in the medium before reaching a new cell [8][9][10][11]. Analogously, the spatial structure of immune responses mediated by IFN and other cytokines can be determined by the diffusion of signaling molecules and the effect of these signals in receptor cells [12]. In addition to infection and immunity, diffusion-reaction models have been used for addressing a variety of eco-evolutionary questions such as, for instance, how bacterial CRISPR limits viral spread [13] and how spatial obstacles perturb population expansion fronts and promote random genetic drift [14].
Another useful approach to investigating viral spread is provided by cell automaton models, in which individual cells are simulated as elements of a grid, making spatial structure explicit [15][16][17][18][19][20][21]. This approach has been previously used for studying how the interaction between viruses and the immune system is influenced by space. For instance, it has been shown that spatial clustering of infected cells (which is a consequence of viral spread in foci) tends to reduce the ability of cytotoxic T-lymphocytes (CTLs) to clear the infection, but that clustering actually increases clearing efficacy if chemotaxis allows CTLs to perform a non-random search for infected cells [22]. Another work explored how spatial structure influences the evolution of CTL escape mutants in HIV-1 [23]. CTL escape mutant should be able to avoid CTL-mediated lysis regardless of the presence of non-mutant viruses in the same host because CTL-mediated lysis is epitope-specific. Nevertheless, CTLs can also secrete cytokines capable of repressing infection in an epitope-independent manner, for instance by reducing the overall susceptibility of cells to viral infection. CTL escape mutants might be sensitive to this mode of action and, consequently, their fitness might be dependent on the presence of CTL-stimulating viruses in the neighborhood.
Here, we investigate how natural selection determines the evolution of viral innate immunity evasion in the presence of spatial structure. For this, we use a cell automaton in which infection and IFN-mediated immunity occur as diffusion-reaction processes. From these simulations, we infer direct and indirect fitness effects associated to IFN shutdown, as well as quantitative descriptors of spatial structure. Within this framework, we analyze how viral demography, the timing of infection and immunity, and the physical properties of the virions and the medium should determine the evolution of IFN evasion.

Model
We considered N S cells susceptible to infection by a virus (uninfected), N E infected cells that are not yet virion-producers (eclipse phase), N P virus-producer cells, and N D cells killed by the virus (Fig 1A). Killed cells were not replaced by new cells and hence became equivalent to empty positions in the grid. We considered two virus variants, W and D. For simplicity, we ignored cells co-infected with both variants. The W variant blocked or did not stimulate IFN production (actively preventing IFN production or avoiding recognition by the innate immune system were equivalent in the context of this model). In contrast, cells infected with the D variant detected the virus and became IFN producers. The model allowed the relative speed of infection and IFN production to vary, such that eclipse-phase cells could be primed for IFN production (N Eε ) or be IFN producers (N Eπ ), and the same applied to D-producing cells (N Pε and N Pπ , respectively).
Based on known dynamics [24][25][26], though, we initially assumed that virion release preceded IFN secretion (N Eπ � 0). We also incorporated the observation that uninfected cells respond to IFN in a dose-dependent manner [5,27] and become immunized (N R ), but do not produce IFN themselves [28]. Infection (N S !N E ) and immunization (N S !N R ) were simulated as Poisson stochastic processes occurring with probability for each cell and simulation time unit Δt (0.5 min). For infection, λ = k V VN S , where k V is the infection rate (infectivity), V the local virion concentration and N S equals 1 or 0. For immunization, λ = k R βN S , where k R is the immunization rate, β the local IFN concentration, and N S equals 1 or 0. All other cellular state transitions were modeled as random processes occurring with cumulative probability P = 0 for t < τ, P = 0.5 for t = τ, and P = 1 for t > τ, where τ is the half time of the corresponding transition (e.g. τ EP for N E ! N P ; Fig 1A). Viral dynamics obeyed a diffusion-reaction process in an orthogonally divided space, the partial derivative of the virion concentration with time in every grid position being The first two terms correspond to the reaction part, where r v is the virion production rate of infected cells. For the D virus, the N P term was replaced with N Pε + N Pπ (that is, all D-producing cells regardless of IFN status). Wherein needed, we used different virion production rates for W and D, r v (W) and r v (D), respectively. The δ V parameter is the virion degradation/outflow rate. We ignored the loss of virions due to adsorption to cells. The last term describes diffusion, where D is the Stokes-Einstein diffusion coefficient, Δx is the grid size and ΔV the virion concentration difference between grid subunits. To calculate it, we considered only  adjacent cells for each grid position (x, y), such that: In the limits of the grid, we used a continuous-system approach, meaning that the neighbors of the rightmost cells were the leftmost cells, and the neighbors of the uppermost cells were the bottommost cells. For reasons of computational efficiency simulations were performed using a time unit Δt = 0.5 min, but a finer time resolution was needed to calculate the diffusion process. We thus used shorter time units for the diffusion part only and assumed interim quasisteady states for all other variables in the system. Analogously, for IFN: Viral spread occurred exclusively through the above diffusion-reaction process (for illustration, S1 Fig shows the diffusion process of virions and IFN alone, without the reaction term). We did not consider other types of spread, such as cell-to-cell spread, active transport of viruses, or cell mobility. Cellular proliferation was also ignored, since the time scale of the infection was shorter than the typical cell division time of normal (non-tumoral) cells, with some exceptions (e.g activated lymphocytes). For simplicity, we also ignored details of the intracellular infection and immunization processes, and we assumed that, in each cell, after the eclipse phase both IFN and virion production was linear with time. In nature, viruses can follow different replication mechanisms [29] and, hence, virion release is not necessarily linear with time. Also, IFN production is regulated by positive feedback loops in the infected cell [30], which we did not consider.
Initially, we assumed that IFN had no effect on already infected cells but we later relaxed this assumption to include autocrine effects as well as post-infection paracrine effects. For this, we allowed infected cells to respond to IFN by undergoing apoptosis (S2 Fig; N A cells). We modeled IFN-triggered cell death (N E !N A , N P !N A ) as a Poisson stochastic process as above, with λ = k A βN P for virion-producing cells and λ = k A βN E for cells in eclipse phase, where k A is the IFN-induced apoptosis rate of infected cells.
The cellular automaton diffusion-reaction simulations were performed using MATLAB R2018b scripts (S1 File).
We calculated the growth rate of each virus variant as where N I(0) is the initial number of infected cells and N I the number of cells infected at time t including both producer and eclipse-phase cells (N I = N E + N P for the W variant; N I = N Eε + N Eπ + N Pε + N Pπ for the D variant). We used log 10 instead of the natural logarithm to facilitate visualization of the results.
To obtain a time-integrated metric of viral fitness, we calculated where N � is the cumulative number of cells infected throughout the simulation and t f the final time point (36 h). Following our previous work [5], we defined the fitness (f) of a focal virus conditional to the presence of other viruses in the same population (neighborhood; Fig 1B). For this, we used the pure W infection (i.e. W virus in a W neighborhood) as reference and we expressed fitness relative to it, such that: where R W|W is the growth rate of a pure W infection, or for time-integrated fitness. Therefore, by definition pure W infections had fitness f = 0. Let us then consider the extreme situation in which, due to the presence of D neighbors, the spread of the W virus would be hampered by as many IFN-immunized cells as in a pure D infection. For this scenario, we defined the fitness of W as Hence, −b is the cost of being in a D neighborhood and, reversely, b is the benefit of being in a W neighborhood. Although b depends on paracrine signaling (that is, the benefit of avoiding IFN secretion), it may also be indirectly influenced by other features of the neighborhood such as, for instance, local virion abundance, which determines competition for cellular resources. Next, let us consider a situation in which the spread of a D virus would not be hampered by the immunization of neighbor cells (D virus in W neighborhood). We defined the fitness of D infections accordingly as f D|W = c. Here, c can be interpreted as the direct effect of blocking IFN for the focal virus, independent of neighborhood. If c > 0, blocking IFN has a direct cost, whereas if c < 0 blocking IFN has a direct benefit. Finally, we defined the fitness of a D virus in a D neighborhood (pure D infection) as f D|D = c−b, which is simply the sum of the independent direct and indirect terms.
Having considered the four extreme cases in which the neighborhood was purely D or purely W, we then allowed for intermediate scenarios (Fig 1B). For this, we defined the overall fitness of W as where r W is a parameter that quantifies how strongly W viruses are influenced by W neighbors.
It follows that where r D quantifies how strongly D viruses are influenced by D neighbors. Thus, f D = c−r D b. Natural selection will favor the W variant and, therefore, will promote the evolution of IFN evasion if f W −f D >0. Consequently: The quantity r W +r D −1 describes the difference between the neighborhoods of W and D. More precisely, r W +r D −1 is the difference between the indirect IFN-mediated effects experienced by W and D viruses. By denoting r = r W +r D −1, the condition for the evolution of IFN shutdown is: This expression corresponds to Hamilton´s rule [31]. Hence, the condition for natural selection to favor viruses that block or avoid IFN signaling can be written in terms of classical social evolution theory. In the social evolution field, r has been often defined as the genetic relatedness between interacting partners for the relevant trait. In our context, such relatedness depends strictly on the spatial structure of infection and immunity. If the two variants are fully insulated, the benefits/costs of preventing/stimulating paracrine signaling will be felt only by viruses of the same kind as the focal virus and, hence, r = 1 (Fig 1B). In contrast, if the innate immune response triggered by the IFN-stimulating variant affects both variants equally, r = 0. If IFN-blocking and IFN-stimulating viruses are spatially segregated to some extent (r > 0), selection may indirectly favor IFN-blocking viruses because they take a greater share of the benefits of preventing paracrine signaling. Notice that all parameters considered (b, c, r W , r D , and r) are time-dependent because fitness and spatial structure vary as the infection and the immune response progress. Also, these parameters are defined at the population level, and are not inferred from an analysis of each individual cell in the grid but, instead, from total cell counts in a given population. An Excel spreadsheet for the calculation of these parameters is available upon request.

Selection for IFN shutdown is determined by spatial structure
We investigated the dynamics of viral spread and innate immune responses in a square grid containing 47,961 cells. Our model considered virion and IFN diffusion as well as cell infection and immunization (Fig 1A). For this, we used the empirically-determined size and degradation rate of VSV particles and the measured diffusion coefficient of chicken IFN. Generally, parameter values corresponded to a rapidly growing lytic virus and adherent IFN-producing cells. Yet, the model is generally applicable to different types of viruses. Details of parameter values are provided in Table 1. The simulated infections progressed as foci as a result of the diffusion-reaction process, reproducing the typical spread mode of many viruses (Fig 2A). In simulations containing only one type of virus variant, IFN-blocking (W) and IFN-stimulating (D) viruses initially spread at similar rates. However, D infections subsequently became halted by innate immunity, whereas W infections progressed and invaded the entire cell population. In mixed infections initiated with an equal input of W and D, the growth of both variants was halted as the immune response was deployed (Fig 2A). Hence, in mixed infections the fitness of W and D relative to a pure W infection (f W and f D , respectively) decreased with time.
The fitness difference between W and D in mixed infections can be expressed as f W −f D = rb −c (see Model; Fig 1B) where b is the indirect benefit of being in a W, IFN suppressive neighborhood, c is the direct fitness effect of suppressing IFN on the focal virus, and r is the difference between the neighborhoods of W and D determined by spatial structure. In the absence of spatial structure, r = 0 and hence f W −f D = -c. Hence, to infer c we performed simulations of mixed infections in which spatial structure was disrupted by equalizing virion concentration in the grid at every time unit Δt (0.5 min). Knowing c, we obtained b by comparing the fitness of a pure D infection relative to a pure W infection, f D|D = c−b (Fig 2B). In these simulations, Spatial structure in viral innate immunity evasion IFN had no effect on already infected cells and, hence, autocrine signaling was not allowed. Consequently, we obtained c = 0. In contrast, b increased with time as the cell population became immunized.
In mixed infections, initially r was close to 1 because each virus was only weakly influenced by neighbors from other infection foci. Yet, as the infection progressed and IFN-mediated immunity expanded throughout the cell population, r-values decreased (Fig 2B). We calculated the parameters describing spatial structure associated to each variant as r W = 1+f W /b, r D = (c−f D )/b, and r = r W +r D −1, where r W describes to what extent the fitness of W viruses is determined by neighbor viruses of the same type W and, analogously, r D describes to what extent the fitness of D viruses is determined by neighbor viruses of the same type D (see Model for details; Fig 1B). We found that r W decayed with time faster than r D , indicating that the neighborhood of both variants gradually became more similar to a pure D infection than to a pure W infection. This shows that IFN functions as a harmful diffusible molecule for the virus, akin to a pollutant or poison. Thus, we define IFN as a "public bad", by contraposition to wellknown diffusible public goods (e.g. microbial secreted enzymes).
We ended simulations at 36 h, a time point at which pure D infections died out, and calculated a time-integrated fitness value using the cumulative number of infected cells (N � ), from which we obtained the corresponding time-integrated values of each parameter in our model ( Table 2). Despite the inhibitory effects of IFN, W remained fitter than D. Because c = 0, this was strictly due to spatial structure which, in turn, emerged from the diffusion-reaction infection and immunity processes and allowed W to be on average less adversely affected by IFNmediated immunity than D. These results suggest that spatial structure is a key determinant of the evolution of IFN suppression in viruses.

Factors promoting spatial insulation of IFN-suppressing and IFNstimulating variants
As shown above, the fitness of IFN-suppressing virus variants depends on their spatial insulation from IFN-stimulating variants. A straightforward factor determining such segregation is the initial fraction of infected cells. To illustrate this, we performed simulations in which we varied the initial number of infected cells from N I(0) = 2 to N I(0) = 1000, keeping constant the Spatial structure in viral innate immunity evasion initial frequency of each variant at 50% as well as the initial number of susceptible cells. We found that for N I(0) = 2, the W variant was largely unaffected by the presence of the D virus, whereas for N I(0) = 1000 the interference was such that both variants showed similarly low fitness (Fig 3). The relative fitness of the D variant also decayed as N I(0) increased due to selfinterference (except for N I(0) = 1000, due to saturation of W spread). The indirect benefit of suppressing IFN secretion (b) tended to increase with N I(0) because larger N I(0) values produced more immunized cells in the absence of IFN suppression. However, spatial structure dropped from r � 1 for N I(0) = 2 to r � 0 for N I(0) = 1000. This effect was driven by a gradual populating of the W neighborhood with IFN-immunized cells, as shown by the decreasing r W values. Spatial structure determined the fitness advantage of W over D through the rb term, such that selection for IFN suppression was stronger for lower N I(0) values. These results suggest that IFN evasion should more easily evolve in viruses that experience severe population bottlenecks during transmission and hence produce few, isolated initial infection foci. Another factor that should determine selection for IFN suppression is medium viscosity, since viscosity determines the diffusion coefficient of IFN and virions. We found that the spatial structure parameter r increased with medium viscosity. The more local nature of infection and immunization at higher viscosities allowed the W variant to be less adversely affected by D viruses, increasing the fitness difference between W and D and thus making selection for IFN suppression stronger (S3 Fig). We also set out to test whether selection for IFN suppression depended on virion size. For this, we performed simulations in which we varied virion size from 50 nm to 500 nm. We found that the spatial structure parameter r tended to increase with larger virion sizes, but that this effect was weak (S4 Fig) because IFN

IFN suppression as an altruistic trait
The above results indicate that IFN shutdown can be favored by natural selection acting on indirect fitness effects, as long as infection and immunity are spatially structured. This suggests that IFN shutdown could evolve even if it imposes a direct fitness cost to the actor (c > 0). From a social evolution perspective, this implies that innate immunity evasion could evolve as an altruistic trait. To illustrate this, we made IFN shutdown costly by assigning a direct advantage to D in terms of virion production rate (higher r v and thus, higher number of virions produced per cell). Initially, D expanded faster than W throughout the cell population. However, as cells became immunized, the spread of D was halted whereas W continued to expand ( S5  Fig). Using the cumulative number of infected cells to calculate fitness as above, we found that W remained fitter than D for a virion production cost of up to twofold, i.e. r v (D) = 2r v (W) (Fig  4). As expected, c-values increased as the advantage of D in terms of virion productivity became larger. Less intuitively, b also increased despite the fact that parameters controlling paracrine innate immunity were not changed. This occurred because D neighborhoods contained more virions than W neighborhoods and hence experienced stronger local competition for cells. Hence, b captured fitness effects associated to spatial structure including but not limited to paracrine signaling. To show this, we repeated the above simulations switching innate immunity off (k R = 0). Under these conditions, D outcompeted W owing to its higher virion production rate. Yet, again, b increased as the excess virion production of D became larger, even if no paracrine response was possible, confirming that b captured the effect of local competition for cells (S6 Fig). Spatial structure in viral innate immunity evasion

Contribution of autocrine signaling to direct fitness
To explore the situation in which IFN blockade provides a direct benefit (c < 0), we allowed for autocrine IFN signaling. To implement this, infected cells responded to IFN by undergoing apoptosis (S2 Fig). This also allowed infected cells to respond to IFN secreted by other cells (post-infection paracrine signaling). Without changing the rate at which IFN made noninfected cells immune (k R ), we increased the rate at which IFN triggered apoptosis (k A ). This showed that, even for relatively strong post-infection IFN effects, autocrine signaling had a minimal impact on the direct fitness component c, which was always close to zero and much smaller than the indirect components b or rb (S7 Fig). The reason for this is that virion production preceded IFN responses, such that IFN signaling on already infected cells was relatively unimportant. This is in line with our recent experimental results, which indicated a low efficacy of IFN (autocrine or paracrine) in already infected cells [5].
In contrast, when we allowed IFN to be released before virions (5 h versus 6 h) and IFNtriggered apoptosis was strong (k A >k R ), b decreased and-c increased by the same amount ( Fig  5). Hence, the indirect benefits of blocking IFN secretion partially became direct benefits. Despite this drop in b, the rb fitness component was largely insensitive to changes in k A . This occurred because, as k A increased, D infections became more localized and thus produced fewer immunized cells. This increased spatial structure and allowed the W virus to be less severely affected by the presence of D, as shown by the greater r W and r values. As a result, the rb and-c terms were similarly large, indicating that spatial structure-dependent indirect fitness effects contributed significantly to the benefits of IFN suppression even in the presence of strong post-infection IFN-induced apoptosis. Spatial structure in viral innate immunity evasion

Discussion
In previous work, we showed that IFN shutdown can be considered a social process because the IFN-stimulating status of a given virus modifies the fitness of neighbor viruses [5]. Specifically, mutants that stimulated IFN production had a negative impact on the fitness of IFNstimulating and IFN-suppressing members of the viral population, whereas IFN blockers provided little or no benefits to IFN stimulators. Hence, from the virus perspective, IFN functions as a diffusible "public bad", and the ability of natural selection to promote IFN shutdown depends critically on spatial structure. Here, we have investigated this process quantitatively using a biophysically well-defined system in which a diffusion-reaction model was implemented in cell automaton simulations. We have shown that spatial structure in the form of infection and immunization foci can emerge from IFN and virion concentration gradients, which in turn depend on demography (e.g. transmission bottlenecks) and the properties of the medium (e.g. viscosity). Because IFN is smaller and hence diffuses faster than virions, immunization can reach larger areas than infection, even if cells release virions faster than IFN. Delayed but effective innate immune responses capable of successfully arresting virus progression by producing a ring of immunized cells around infection foci have been experimentally shown previously using IFN-stimulating virus variants [24].
We found that, whenever virion release precedes IFN secretion, the fitness associated to IFN evasion should not be strongly determined by direct effects exerted on the infected cell, but rather by indirect (neighborhood) effects. An important consequence of indirect fitness effects being larger than direct effects is that IFN evasion could be favored by natural selection even if costly for the actor, meaning that IFN evasion could evolve as an altruistic trait. Cooperation and altruism have long been studied by evolutionary biologists [31,[36][37][38][39][40]. In principle, natural selection should not favor traits that diminish individual fitness because cooperators can be invaded by cheaters that benefit from cooperative traits without reciprocating. However, cooperation breakdown can be avoided if cooperators tend to interact preferentially with other cooperators. Spatial structure allows for such assortment in a wide variety of organisms and is hence considered a fundamental factor driving the evolution of cooperation [41][42][43][44]. Spatial structure should actually be the main driver of cooperation in simple organisms for which other mechanisms of assortment such as learning and memory are not possible, as is the case of viruses. Here, we have shown that particle diffusion provides a simple and general physical basis for cooperator assortment in viruses. Within the context of innate immunity evasion, in mixed populations containing both IFN-stimulating and IFN-suppressing virus variants, IFN and virion concentration gradients allow IFN suppressors to be less severely affected by immunization than IFN stimulators.
The role of spatial structure in the evolution of innate immunity evasion should be particularly critical when shutting down IFN entails a direct fitness cost, since IFN suppression then becomes an altruistic trait. In future work, it would be interesting to elucidate the direct and indirect fitness effects of IFN suppression in different viruses, and the mechanisms involved in these effects. Yet, we found that in the absence of such direct costs, and even when IFN suppression has direct benefits, the evolution of IFN shutdown should depend strongly on spatial structure, provided that the ability of IFN to halt viral spread is mainly determined by paracrine signaling. As we have shown, the importance of spatial structure increases if virion release precedes IFN signaling. Currently available empirical evidence supports the view that this is the case for different viruses [24,25]. Although further work is required to establish more general conclusions, virion release prior to IFN signaling should be frequent in viral infections because rapid replication and progeny release tend to be critical for the survival of many viruses.
Spatial structure is pervasive in viral infections. In addition to the basic diffusion-reaction process investigated here, spatial structure can result from other processes such as localized cell-to-cell virus transfer, subversion of cellular structures for enhancing viral spread, use of mobile cells for disseminating the infection, limited virus trafficking among organs or body compartments, inter-host transmission bottlenecks, and so on [45][46][47][48][49][50]. In future work it would be interesting to investigate the validity of our results in these more complex frameworks. Other complex scenarios, such as tridimensional structures and the spread of the virus and cytokines through vessels could be considered. A common qualitative observation, though, is that viral infections often exhibit metapopulation structures in which quasi-isolated demes are initiated by small number of founder viral particles, increasing genetic drift and intra-deme genetic relatedness. This should make the evolution of cooperative traits such as IFN evasion more likely. In some cases, though, spatial structure might be less evident, such as in blood-borne viruses at the intra-host level. Potentially, selection for IFN evasion might be less efficient in these viruses, or these viruses may utilize different IFN evasion strategies, such as antagonizing downstream antiviral responses triggered by IFN signaling in infected cells. green. An eclipse phase between viral sensing (N ε ) and IFN secretion (N π ) is also considered. The relative speed of infection and IFN production can vary. Hence, D-infected cells can be in four possible stages (N Eε , N Pε , N Eπ , N Pπ ). k v : virion infectivity (infection rate). k R : IFN immunization rate. τ EP : viral eclipse half time. τ επ : half time between infection and IFN secretion. τ PD : half time between virion production and cell death. Thus τ PD + τ EP , is the total duration of the infection cycle. r v : virion production rate. K = r v τ PD is thus the number of virions produced per infected cell. δ V : virion degradation/outflow rate. r β : IFN production rate of immunized cells. δ β : IFN degradation/outflow rate. Infected cells also respond to IFN by undergoing apoptosis (N A cells).