When does pathogen evolution maximize the basic reproductive number in well-mixed host-pathogen systems?

The role of a programmatic immune response on the evolution of pathogen traits

In modelling pathogen evolution during epidemics, it is important to understand the interactions between within-host infection dynamics and between-host pathogen transmission. Multiscale models often assume an immune response that is highly responsive to pathogen dynamics. Empirical evidence, however, suggests that the immune response in acute infections is triggered and programmatic. This leads to somewhat more predictable infection trajectories where transition times and, consequently, the infectious window are non-exponentially distributed. Here, we develop a within-host model where the immune response is triggered by pathogen growth but otherwise develops independently, and use this to obtain analytic expressions for the infectious period and peak pathogen load. This allows us to model the basic reproductive number in terms of explicit functional relationships among within-host traits including the growth rate of the pathogen. We find that the dependence of pathogen load and the infectious window on within-host parameters constrains the evolution of the pathogen growth rate. At low growth rate, selection favours a higher pathogen load and therefore increasing pathogen growth rate. At high growth rates, selection for a longer infectious window trades off against selection against the effects of virulence. At intermediate growth rates the basic reproductive number is relatively insensitive to changes in the growth rate. The resulting "flat" region of the pathogen fitness landscape is due to the stability of the programmatic immune response in clearing the infection.

A Mathematical Framework for Predicting Lifestyles of Viral Pathogens

Despite being similar in structure, functioning, and size, viral pathogens enjoy very different, usually well-defined ways of life. They occupy their hosts for a few days (influenza), for a few weeks (measles), or even lifelong (HCV), which manifests in acute or chronic infections. The various transmission routes (airborne, via direct physical contact, etc.), degrees of infectiousness (referring to the viral load required for transmission), antigenic variation/immune escape and virulence define further aspects of pathogenic lifestyles. To survive, pathogens must infect new hosts; the success determines their fitness. Infection happens with a certain likelihood during contact of hosts, where contact can also be mediated by vectors. Besides structural aspects of the host-contact network, three parameters appear to be key: the contact rate and the infectiousness during contact, which encode the mode of transmission, and third the immunity of susceptible hosts. On these grounds, what can be said about the reproductive success of viral pathogens? This is the biological question addressed in this paper. The answer extends earlier results of the author and makes explicit connection to another basic work on the evolution of pathogens. A mathematical framework is presented that models intra- and inter-host dynamics in a minimalistic but unified fashion covering a broad spectrum of viral pathogens, including those that cause flu-like infections, childhood diseases, and sexually transmitted infections. These pathogens turn out as local maxima of numerically simulated fitness landscapes. The models involve differential and integral equations, agent-based simulation, networks, and probability.

Why is sterility virulence most common in sexually transmitted infections? Examining the role of epidemiology

Sterility virulence, or the reduction in host fecundity due to infection, occurs in many host-pathogen systems. Notably, sterility virulence is more common for sexually transmitted infections (STIs) than for directly transmitted pathogens, while other forms of virulence tend to be limited in STIs. This has led to the suggestion that sterility virulence may have an adaptive explanation. By focusing upon finite population models, we show that the observed patterns of sterility virulence can be explained by consideration of the epidemiological differences between STIs and directly transmitted pathogens. In particular, when pathogen transmission is predominantly density invariant (as for STIs), and mortality is density dependent, sterility virulence can be favored by demographic stochasticity, whereas if pathogen transmission is predominantly density dependent, as is common for most directly transmitted pathogens, sterility virulence is disfavored. We show these conclusions can hold even if there is a weak selective advantage to sterilizing.

Viral invasion fitness across a continuum from lysis to latency

The prevailing paradigm in ecological studies of viruses and their microbial hosts is that the reproductive success of viruses depends on the proliferation of the 'predator', that is, the virus particle. Yet, viruses are obligate intracellular parasites, and the virus genome-the actual unit of selection-can persist and proliferate from one cell generation to the next without lysis or the production of new virus particles. Here, we propose a theoretical framework to quantify the invasion fitness of viruses using an epidemiological cell-centric metric that focuses on the proliferation of viral genomes inside cells instead of virus particles outside cells. This cell-centric metric enables direct comparison of viral strategies characterized by obligate killing of hosts (e.g. via lysis), persistence of viral genomes inside hosts (e.g. via lysogeny), and strategies along a continuum between these extremes (e.g. via chronic infections). As a result, we can identify environmental drivers, life history traits, and key feedbacks that govern variation in viral propagation in nonlinear population models. For example, we identify threshold conditions given relatively low densities of susceptible cells and relatively high growth rates of infected cells in which lysogenic and other chronic strategies have higher potential viral reproduction than lytic strategies. Altogether, the theoretical framework helps unify the ongoing study of eco-evolutionary drivers of viral strategies in natural environments.

Pathogen evolution in finite populations: slow and steady spreads the best

The theory of life-history evolution provides a powerful framework to understand the evolutionary dynamics of pathogens. It assumes, however, that host populations are large and that one can neglect the effects of demographic stochasticity. Here, we expand the theory to account for the effects of finite population size on the evolution of pathogen virulence. We show that demographic stochasticity introduces additional evolutionary forces that can qualitatively affect the dynamics and the evolutionary outcome. We discuss the importance of the shape of the pathogen fitness landscape on the balance between mutation, selection and genetic drift. This analysis reconciles Adaptive Dynamics with population genetics in finite populations and provides a new theoretical toolbox to study life-history evolution in realistic ecological scenarios.

On the Reproduction Number of a Gut Microbiota Model

A spatially structured linear model of the growth of intestinal bacteria is analysed from two generational viewpoints. Firstly, the basic reproduction number associated with the bacterial population, i.e. the expected number of daughter cells per bacterium, is given explicitly in terms of biological parameters. Secondly, an alternative quantity is introduced based on the number of bacteria produced within the intestine by one bacterium originally in the external media. The latter depends on the parameters in a simpler way and provides more biological insight than the standard reproduction number, allowing the design of experimental procedures. Both quantities coincide and are equal to one at the extinction threshold, below which the bacterial population becomes extinct. Optimal values of both reproduction numbers are derived assuming parameter trade-offs.

Pathogen evolution under host avoidance plasticity

Host resistance consists of defences that limit pathogen burden, and can be classified as either adaptations targeting recovery from infection or those focused upon infection avoidance. Conventional theory treats avoidance as a fixed strategy which does not vary from one interaction to the next. However, there is increasing empirical evidence that many avoidance strategies are triggered by external stimuli, and thus should be treated as phenotypically plastic responses. Here, we consider the implications of avoidance plasticity for host-pathogen coevolution. We uncover a number of predictions challenging current theory. First, in the absence of pathogen trade-offs, plasticity can restrain pathogen evolution; moreover, the pathogen exploits conditions in which the host would otherwise invest less in resistance, causing resistance escalation. Second, when transmission trades off with pathogen-induced mortality, plasticity encourages avirulence, resulting in a superior fitness outcome for both host and pathogen. Third, plasticity ensures the sterilizing effect of pathogens has consequences for pathogen evolution. When pathogens castrate hosts, selection forces them to minimize mortality virulence; moreover, when transmission trades off with sterility alone, resistance plasticity is sufficient to prevent pathogens from evolving to fully castrate.

Characterization of the endemic equilibrium and response to mutant injection in a multi-strain disease model

We explore a model of an antigenically diverse infection whose otherwise identical strains compete through cross-immunity. We assume that individuals may produce upon infection different numbers of antibody types, each of which matches the antigenic configuration of a particular epitope, and that one matching antibody type grants total immunity against a challenging strain. In order to reduce the number of equations involved in the analytic description of the dynamics, we follow the strategy proposed by Kryazhimskiy et al. (2007) and apply a low-order closure reminiscent of a pair approximation. Using this approximation, we go beyond the numerical studies of Kryazhimskiy et al. (2007) and explore the analytic properties of the ensuing model in the absence of mutation. We characterize its endemic equilibrium, comparing with the results of agent based simulations of the full model to assess the performance of the closure assumption. We show that a particular choice of immune response leads to a degenerate endemic equilibrium, where different strain prevalences may exist, breaking the symmetry of the model. Finally we study the behavior of the system under the injection of mutant strains. We find that the build up of diversity from a single founding strain is extremely unlikely for different choices of the population׳s immune response.

Evolution of pathogen virulence across space during an epidemic

We explore pathogen virulence evolution during the spatial expansion of an infectious disease epidemic in the presence of a novel host movement trade-off, using a simple, spatially explicit mathematical model. This work is motivated by empirical observations of the Mycoplasma gallisepticum invasion into North American house finch (Haemorhous mexicanus) populations; however, our results likely have important applications to other emerging infectious diseases in mobile hosts. We assume that infection reduces host movement and survival and that across pathogen strains the severity of these reductions increases with pathogen infectiousness. Assuming these trade-offs between pathogen virulence (host mortality), pathogen transmission, and host movement, we find that pathogen virulence levels near the epidemic front (that maximize wave speed) are lower than those that have a short-term growth rate advantage or that ultimately prevail (i.e., are evolutionarily stable) near the epicenter and where infection becomes endemic (i.e., that maximize the pathogen basic reproductive ratio). We predict that, under these trade-offs, less virulent pathogen strains will dominate the periphery of an epidemic and that more virulent strains will increase in frequency after invasion where disease is endemic. These results have important implications for observing and interpreting spatiotemporal epidemic data and may help explain transient virulence dynamics of emerging infectious diseases.

Evolutionary dynamics of infectious diseases in finite populations

In infectious disease epidemiology the basic reproductive ratio, R0, is defined as the average number of new infections caused by a single infected individual in a fully susceptible population. Many models describing competition for hosts between non-interacting pathogen strains in an infinite population lead to the conclusion that selection favors invasion of new strains if and only if they have higher R0 values than the resident. Here we demonstrate that this picture fails in finite populations. Using a simple stochastic SIS model, we show that in general there is no analogous optimization principle. We find that successive invasions may in some cases lead to strains that infect a smaller fraction of the host population, and that mutually invasible pathogen strains exist. In the limit of weak selection we demonstrate that an optimization principle does exist, although it differs from R0 maximization. For strains with very large R0, we derive an expression for this local fitness function and use it to establish a lower bound for the error caused by neglecting stochastic effects. Furthermore, we apply this weak selection limit to investigate the selection dynamics in the presence of a trade-off between the virulence and the transmission rate of a pathogen.