Influenza A gradual and epochal evolution: insights from simple models

Mathematical modeling for estimating influenza vaccine efficacy: A case study of the Valencian Community, Spain

Vaccine efficacy and its quantification is a crucial concept for the proper design of public health vaccination policies. In this work we proposed a mathematical model to estimate the efficacy of the influenza vaccine in a real-word scenario. In particular, our model is a SEIR-type epidemiological model, which distinguishes vaccinated and unvaccinated populations. Mathematically, its dynamics is governed by a nonlinear system of ordinary differential equations, where the non-linearity arises from the effective contacts between susceptible and infected individuals. Two key aspects of this study is that we use a vaccine distribution over time that is based on real data specific to the elderly people in the Valencian Community and the calibration process takes into account that over one influenza season a specific proportion of the population becomes infected with influenza. To consider the effectiveness of the vaccine, the model incorporates a parameter, the vaccine attenuation factor, which is related with the vaccine efficacy against the influenza virus. With this framework, in order to calibrate the model parameters and to obtain an influenza vaccine efficacy estimation, we considered the 2016-2017 influenza season in the Valencian Community, Spain, using the influenza reported cases of vaccinated and unvaccinated. In order to ensure the model identifiability, we choose to deterministically calibrate the parameters for different scenarios and we find the one with the minimum error in order to determine the vaccine efficacy. The calibration results suggest that the influenza vaccine developed for 2016-2017 influenza season has an efficacy of approximately 76.7%, and that the risk of becoming infected is five times higher for an unvaccinated individual in comparison with a vaccinated one. This estimation partially agrees with some previous studies related to the influenza vaccine. This study presents a new integrated mathematical approach to study the influenza vaccine efficacy and gives further insight into this important public health topic.

On the probability of strain invasion in endemic settings: Accounting for individual heterogeneity and control in multi-strain dynamics

•

Endemic infection can insulate host populations from invasion by mutant variants.

•

The timing of control implementation strongly influences its efficacy.

•

Controls that exacerbate host heterogeneity outperform those that curtail it.

•

Differential control can facilitate strain invasion and eventual replacement.

Pathogen evolution is an imminent threat to global health that has warranted, and duly received, considerable attention within the medical, microbiological and modelling communities. Outbreaks of new pathogens are often ignited by the emergence and transmission of mutant variants descended from wild-type strains circulating in the community. In this work we investigate the stochastic dynamics of the emergence of a novel disease strain, introduced into a population in which it must compete with an existing endemic strain. In analogy with past work on single-strain epidemic outbreaks, we apply a branching process approximation to calculate the probability that the new strain becomes established. As expected, a critical determinant of the survival prospects of any invading strain is the magnitude of its reproduction number relative to that of the background endemic strain. Whilst in most circumstances this ratio must exceed unity in order for invasion to be viable, we show that differential control scenarios can lead to less-fit novel strains invading populations hosting a fitter endemic one. This analysis and the accompanying findings will inform our understanding of the mechanisms that have led to past instances of successful strain invasion, and provide valuable lessons for thwarting future drug-resistant strain incursions.

Phylodynamic theory of persistence, extinction and speciation of rapidly adapting pathogens

Rapidly evolving pathogens like influenza viruses can persist by changing their antigenic properties fast enough to evade the adaptive immunity, yet they rarely split into diverging lineages. By mapping the multi-strain Susceptible-Infected-Recovered model onto the traveling wave model of adapting populations, we demonstrate that persistence of a rapidly evolving, Red-Queen-like state of the pathogen population requires long-ranged cross-immunity and sufficiently large population sizes. This state is unstable and the population goes extinct or 'speciates' into two pathogen strains with antigenic divergence beyond the range of cross-inhibition. However, in a certain range of evolutionary parameters, a single cross-inhibiting population can exist for times long compared to the time to the most recent common ancestor ([Formula: see text]) and gives rise to phylogenetic patterns typical of influenza virus. We demonstrate that the rate of speciation is related to fluctuations of [Formula: see text] and construct a 'phase diagram' identifying different phylodynamic regimes as a function of evolutionary parameters.

Chimera states in multi-strain epidemic models with temporary immunity

We investigate a time-delayed epidemic model for multi-strain diseases with temporary immunity. In the absence of cross-immunity between strains, dynamics of each individual strain exhibit emergence and annihilation of limit cycles due to a Hopf bifurcation of the endemic equilibrium, and a saddle-node bifurcation of limit cycles depending on the time delay associated with duration of temporary immunity. Effects of all-to-all and non-local coupling topologies are systematically investigated by means of numerical simulations, and they suggest that cross-immunity is able to induce a diverse range of complex dynamical behaviors and synchronization patterns, including discrete traveling waves, solitary states, and amplitude chimeras. Interestingly, chimera states are observed for narrower cross-immunity kernels, which can have profound implications for understanding the dynamics of multi-strain diseases.

Human mobility networks and persistence of rapidly mutating pathogens

Rapidly mutating pathogens may be able to persist in the population and reach an endemic equilibrium by escaping hosts' acquired immunity. For such diseases, multiple biological, environmental and population-level mechanisms determine the dynamics of the outbreak, including pathogen's epidemiological traits (e.g. transmissibility, infectious period and duration of immunity), seasonality, interaction with other circulating strains and hosts' mixing and spatial fragmentation. Here, we study a susceptible-infected-recovered-susceptible model on a metapopulation where individuals are distributed in sub-populations connected via a network of mobility flows. Through extensive numerical simulations, we explore the phase space of pathogen's persistence and map the dynamical regimes of the pathogen following emergence. Our results show that spatial fragmentation and mobility play a key role in the persistence of the disease whose maximum is reached at intermediate mobility values. We describe the occurrence of different phenomena including local extinction and emergence of epidemic waves, and assess the conditions for large-scale spreading. Findings are highlighted in reference to previous studies and to real scenarios. Our work uncovers the crucial role of hosts' mobility on the ecological dynamics of rapidly mutating pathogens, opening the path for further studies on disease ecology in the presence of a complex and heterogeneous environment.

Host immunity and pathogen diversity: A computational study

The distinctive features of human influenza A phylogeny have inspired many mathematical and computational studies of viral infections spreading in a host population, but our understanding of the mechanisms that shape the coupled evolution of host immunity, disease incidence and viral antigenic properties is far from complete. In this paper we explore the epidemiology and the phylogeny of a rapidly mutating pathogen in a host population with a weak immune response, that allows re-infection by the same strain and provides little cross-immunity. We find that mutation generates explosive diversity and that, as diversity grows, the system is driven to a very high prevalence level. This is in stark contrast with the behavior of similar models where mutation gives rise to a large epidemic followed by disease extinction, under the assumption that infection with a strain provides lifelong immunity. For low mutation rates, the behavior of the system shows the main qualitative features of influenza evolution. Our results highlight the importance of heterogeneity in the human immune response for understanding influenza A phenomenology. They are meant as a first step toward computationally affordable, individual based models including more complex host-pathogen interactions.

A cross-immunization model for the extinction of old influenza strains

Given the frequent mutation of antigenic features, the constancy of genetic and antigenic diversity of influenza within a subtype is surprising. While the emergence of new strains and antigenic features is commonly attributed to selection by the human immune system, the mechanism that ensures the extinction of older strains remains controversial. To replicate this dynamics of replacement current models utilize mechanisms such as short-lived strain-transcending immunity, a direct competition for hosts, stochastic extinction or constrained antigenic evolution. Building on the idea of short-lived immunity we introduce a minimal model that exhibits the aforementioned dynamics of replacement. Our model relies only on competition due to an antigen specific immune-response in an unconstrained antigenic space. Furthermore the model explains the size of typical influenza epidemics as well as the tendency that new epidemics are associated with mutations of old antigens.

The effects of a deleterious mutation load on patterns of influenza A/H3N2's antigenic evolution in humans

Recent phylogenetic analyses indicate that RNA virus populations carry a significant deleterious mutation load. This mutation load has the potential to shape patterns of adaptive evolution via genetic linkage to beneficial mutations. Here, we examine the effect of deleterious mutations on patterns of influenza A subtype H3N2's antigenic evolution in humans. By first analyzing simple models of influenza that incorporate a mutation load, we show that deleterious mutations, as expected, act to slow the virus's rate of antigenic evolution, while making it more punctuated in nature. These models further predict three distinct molecular pathways by which antigenic cluster transitions occur, and we find phylogenetic patterns consistent with each of these pathways in influenza virus sequences. Simulations of a more complex phylodynamic model further indicate that antigenic mutations act in concert with deleterious mutations to reproduce influenza's spindly hemagglutinin phylogeny, co-circulation of antigenic variants, and high annual attack rates.

DOI:http://dx.doi.org/10.7554/eLife.07361.001

Each year, up to 15% of the world's population experience symptoms of an influenza infection, also commonly known as flu. The most common culprit is a strain of the virus called influenza type A subtype H3N2. One reason that so many people become infected each year is that this virus evolves rapidly. Within a few years, proteins on the surface of the virus known as antigens become less recognizable to the immune system of a person who has been previously infected. This means that the person can become ill with the virus again because their immune system cannot mount an effective response to the evolved virus strain.

Influenza virus strains evolve rapidly because their genetic material accumulates mutations quickly. Although some of these mutations are beneficial to the virus, other mutations are harmful and reduce the ability of the virus to spread. Sometimes beneficial mutations may occur alongside harmful ones, but it is not known how the harmful mutations affect the evolution of the virus.

Here, Koelle and Rasmussen used computer models of H3N2 influenza to examine the effect of harmful mutations on the evolution of this virus population. The models show that harmful mutations limit how quickly the antigens can evolve. Also, the presence of these harmful mutations effectively acts as a sieve: they allow only large changes in the antigens to establish in the virus population.

The models suggest that there are three routes by which large changes in the antigens on H3N2 viruses may occur. The first is by a single mutation that has a big effect on the antigens in viruses that only carry a few harmful mutations, but these large mutations would not happen very often. Another route may be through more common mutations that have only a small or moderate benefit, which would allow the virus to become more common in the population before it acquires a beneficial mutation with a much greater effect. The third possibility is that a large beneficial mutation may arise in viruses that have many harmful mutations. These harmful mutations may initially limit the ability of the virus to spread, but over time, some of these harmful mutations may then be lost.

Koelle and Rasmussen found that the computer models could recreate the patterns of virus evolution that have been observed in real strains of H3N2. Researchers use predictions of influenza evolution to help them decide which virus strains should be included in flu vaccines each year. Koelle and Rasmussen findings indicate that harmful mutations should be considered when making these predictions.

DOI:http://dx.doi.org/10.7554/eLife.07361.002

The impact of aggregating serogroups in dynamic models of Neisseria meningitidis transmission

Background

Neisseria meningitidis (Nm) is a pathogen of multiple serogroups that is highly prevalent in many populations. Serogroups associated with invasive meningococcal disease (IMD) in Canada, for example, include A, B, C, W-135, X and Y. IMD is a rare but serious outcome of Nm infection, and can be prevented with vaccines that target certain serogroups. This has stimulated the development of dynamic models to evaluate vaccine impact. However, these models typically aggregate the various Nm serogroups into a small number of combined groups, instead of modelling each serogroup individually. The impact of aggregation on dynamic Nm model predictions is poorly understood. Our objective was to explore the impact of aggregation on dynamic model predictions.

Methods

We developed two age-structured agent-based models--a 2-strain model and a 4-strain model--to simulate vaccination programs in the Canadian setting. The 2-strain model was used to explore two different groupings: C, versus all other serogroups combined; and B, versus all other serogroups combined. The 4-strain model used the four groupings: C, B, Neisseria lactamica, versus all other serogroups combined. We compared the predicted impact of monovalent C vaccine, quadrivalent ACWY vaccine (MCV-4), and monovalent B vaccine (4CMenB) on the prevalence of serogroup carriage under these different models.

Results

The 2-strain and 4-strain models predicted similar overall impacts of vaccines on carriage prevalence, especially with respect to the vaccine-targeted serogroups. However, there were some significant quantitative and qualitative differences. Declines in vaccine-targeted serogroups were more rapid in the 2-strain model than the 4-strain model, for both the C and the 4CMenB vaccines. Sustained oscillations, and evidence for multiple attractors (i.e., different types of dynamics for the same model parameters but different initial conditions), occurred in the 4-strain model but not the 2-strain model. Strain replacement was also more pronounced in the 4-strain model, on account of the 4-strain model spreading prevalence more thinly across groups and thus enhancing competitive interactions.

Conclusions

Simplifying assumptions like aggregation of serogroups can have significant impacts on dynamic model predictions. Modellers should carefully weigh the advantages and disadvantages of aggregation when formulating models for multi-strain pathogens.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-015-1015-8) contains supplementary material, which is available to authorized users.

Estimating dynamic transmission model parameters for seasonal influenza by fitting to age and season-specific influenza-like illness incidence

Dynamic transmission models are essential to design and evaluate control strategies for airborne infections. Our objective was to develop a dynamic transmission model for seasonal influenza allowing to evaluate the impact of vaccinating specific age groups on the incidence of infection, disease and mortality. Projections based on such models heavily rely on assumed 'input' parameter values. In previous seasonal influenza models, these parameter values were commonly chosen ad hoc, ignoring between-season variability and without formal model validation or sensitivity analyses. We propose to directly estimate the parameters by fitting the model to age-specific influenza-like illness (ILI) incidence data over multiple influenza seasons. We used a weighted least squares (WLS) criterion to assess model fit and applied our method to Belgian ILI data over six influenza seasons. After exploring parameter importance using symbolic regression, we evaluated a set of candidate models of differing complexity according to the number of season-specific parameters. The transmission parameters (average R0, seasonal amplitude and timing of the seasonal peak), waning rates and the scale factor used for WLS optimization, influenced the fit to the observed ILI incidence the most. Our results demonstrate the importance of between-season variability in influenza transmission and our estimates are in line with the classification of influenza seasons according to intensity and vaccine matching.

Capturing the dynamics of pathogens with many strains

Pathogens that consist of multiple antigenic variants are a serious public health concern. These infections, which include dengue virus, influenza and malaria, generate substantial morbidity and mortality. However, there are considerable theoretical challenges involved in modelling such infections. As well as describing the interaction between strains that occurs as a result cross-immunity and evolution, models must balance biological realism with mathematical and computational tractability. Here we review different modelling approaches, and suggest a number of biological problems that are potential candidates for study with these methods. We provide a comprehensive outline of the benefits and disadvantages of available frameworks, and describe what biological information is preserved and lost under different modelling assumptions. We also consider the emergence of new disease strains, and discuss how models of pathogens with multiple strains could be developed further in future. This includes extending the flexibility and biological realism of current approaches, as well as interface with data.

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.

Spreading of multiple epidemics with cross immunization

Pathogen-host relationships are the result of an ongoing coevolutionary race where the immune system of the host attempts to eliminate the pathogen, while the successful pathogen mutates to become invisible for the host's immune system. We here propose a minimal pathogen-host evolution model that takes into account cross immunization and allows for evolution of a spatially heterogeneous immune status of a population of hosts. With only the mutation rate as a determining parameter, the model allows us to produce an evolutionary tree of diseases which is highly branched, but hardly ever splits into separate long-lived trunks. Side branches remain short lived and seldom diverge to the extent of losing all cross immunizations.

The impact of past epidemics on future disease dynamics

Many pathogens spread primarily via direct contact between infected and susceptible hosts. Thus, the patterns of contacts or contact network of a population fundamentally shape the course of epidemics. While there is a robust and growing theory for the dynamics of single epidemics in networks, we know little about the impacts of network structure on long-term epidemic or endemic transmission. For seasonal diseases like influenza, pathogens repeatedly return to populations with complex and changing patterns of susceptibility and immunity acquired through prior infection. Here, we develop two mathematical approaches for modeling consecutive seasonal outbreaks of a partially-immunizing infection in a population with contact heterogeneity. Using methods from percolation theory we consider both leaky immunity, where all previously infected individuals gain partial immunity, and polarized immunity, where a fraction of previously infected individuals are fully immune. By restructuring the epidemiologically active portion of their host population, such diseases limit the potential of future outbreaks. We speculate that these dynamics can result in evolutionary pressure to increase infectiousness.

Role of environmental persistence in pathogen transmission: a mathematical modeling approach

Although diseases such as influenza, tuberculosis and SARS are transmitted through an environmentally mediated mechanism, most modeling work on these topics is based on the concepts of infectious contact and direct transmission. In this paper we use a paradigm model to show that environmental transmission appears like direct transmission in the case where the pathogen persists little time in the environment. Furthermore, we formulate conditions for the validity of this modeling approximation and we illustrate them numerically for the cases of cholera and influenza. According to our results based on recently published parameter estimates, the direct transmission approximation fails for both cholera and influenza. While environmental transmission is typically chosen over direct transmission in modeling cholera, this is not the case for influenza.

Explaining rapid reinfections in multiple-wave influenza outbreaks: Tristan da Cunha 1971 epidemic as a case study

Influenza usually spreads through the human population in multiple-wave outbreaks. Successive reinfection of individuals over a short time interval has been explicitly reported during past pandemics. However, the causes of rapid reinfection and the role of reinfection in driving multiple-wave outbreaks remain poorly understood. To investigate these issues, we focus on a two-wave influenza A/H3N2 epidemic that occurred on the remote island of Tristan da Cunha in 1971. Over 59 days, 273 (96%) of 284 islanders experienced at least one attack and 92 (32%) experienced two attacks. We formulate six mathematical models invoking a variety of antigenic and immunological reinfection mechanisms. Using a maximum-likelihood analysis to confront model predictions with the reported incidence time series, we demonstrate that only two mechanisms can be retained: some hosts with either a delayed or deficient humoral immune response to the primary influenza infection were reinfected by the same strain, thus initiating the second epidemic wave. Both mechanisms are supported by previous empirical studies and may arise from a combination of genetic and ecological causes. We advocate that a better understanding and account of heterogeneity in the human immune response are essential to analysis of multiple-wave influenza outbreaks and pandemic planning.

The impact of seasonal and year-round transmission regimes on the evolution of influenza A virus

Punctuated antigenic change is believed to be a key element in the evolution of influenza A; clusters of antigenically similar strains predominate worldwide for several years until an antigenically distant mutant emerges and instigates a selective sweep. It is thought that a region of East-Southeast Asia with year-round transmission acts as a source of antigenic diversity for influenza A and seasonal epidemics in temperate regions make little contribution to antigenic evolution. We use a mathematical model to examine how different transmission regimes affect the evolutionary dynamics of influenza over the lifespan of an antigenic cluster. Our model indicates that, in non-seasonal regions, mutants that cause significant outbreaks appear before the peak of the wild-type epidemic. A relatively large proportion of these mutants spread globally. In seasonal regions, mutants that cause significant local outbreaks appear each year before the seasonal peak of the wild-type epidemic, but only a small proportion spread globally. The potential for global spread is strongly influenced by the intensity of non-seasonal circulation and coupling between non-seasonal and seasonal regions. Results are similar if mutations are neutral, or confer a weak to moderate antigenic advantage. However, there is a threshold antigenic advantage, depending on the non-seasonal transmission intensity, beyond which mutants can escape herd immunity in the non-seasonal region and there is a global explosion in diversity. We conclude that non-seasonal transmission regions are fundamental to the generation and maintenance of influenza diversity owing to their epidemiology. More extensive sampling of viral diversity in such regions could facilitate earlier identification of antigenically novel strains and extend the critical window for vaccine development.

Consequences of host heterogeneity, epitope immunodominance, and immune breadth for strain competition

Consumer-resource dynamics of hosts with their pathogens are modulated by complex interactions between various branches of hosts' immune systems and the imperfectly perceived pathogen. Multistrain SIR models tend to sweep competitive interaction terms between different pathogen strains into a single parameter representing cross-immunity. After reviewing several hypotheses about the generation of immune responses, we look into the consequences of assuming that hosts with identical immune repertoires respond to new pathogens identically. In particular, we vary the breadth of the typical immune response, or the average number of pathogen epitopes a host perceives, and the probability of perceiving a particular epitope. The latter quantity in our model is equivalent both to the degree of diversity in host responses at the population level and the relative immunodominance of different epitopes. We find that a sharp transition to strain coexistence occurs as host responses become narrow or skewed toward one epitope. Increasing the breadth of the immune response and the immunogenicity of different epitopes typically increases the range of cross-immunity values in which chaotic strain dynamics and competitive exclusion occur. Models attempting to predict the outcomes of strain competition should thus consider the potential diversity and specificity of hosts' responses to infection.

A two-tiered model for simulating the ecological and evolutionary dynamics of rapidly evolving viruses, with an application to influenza

Understanding the epidemiological and evolutionary dynamics of rapidly evolving pathogens is one of the most challenging problems facing disease ecologists today. To date, many mathematical and individual-based models have provided key insights into the factors that may regulate these dynamics. However, in many of these models, abstractions have been made to the simulated sequences that limit an effective interface with empirical data. This is especially the case for rapidly evolving viruses in which de novo mutations result in antigenically novel variants. With this focus, we present a simple two-tiered 'phylodynamic' model whose purpose is to simulate, along with case data, sequence data that will allow for a more quantitative interface with observed sequence data. The model differs from previous approaches in that it separates the simulation of the epidemiological dynamics (tier 1) from the molecular evolution of the virus's dominant antigenic protein (tier 2). This separation of phenotypic dynamics from genetic dynamics results in a modular model that is computationally simpler and allows sequences to be simulated with specifications such as sequence length, nucleotide composition and molecular constraints. To illustrate its use, we apply the model to influenza A (H3N2) dynamics in humans, influenza B dynamics in humans and influenza A (H3N8) dynamics in equine hosts. In all three of these illustrative examples, we show that the model can simulate sequences that are quantitatively similar in pattern to those empirically observed. Future work should focus on statistical estimation of model parameters for these examples as well as the possibility of applying this model, or variants thereof, to other host-virus systems.