Influenza emergence in the face of evolutionary constraints

Waning immunity drives respiratory virus evolution and reinfection

Reinfections with respiratory viruses such as influenza viruses and coronaviruses are thought to be driven by ongoing antigenic immune escape in the viral population. However, this does not explain why antigenic variation is frequently observed in these viruses relative to viruses such as measles that undergo systemic replication. Here, we suggest that the rapid rate of waning immunity in the respiratory tract is the key driver of antigenic evolution in respiratory viruses. Waning immunity results in hosts with immunity levels that protect against homologous reinfection but are insufficient to protect against infection with a heterologous, antigenically different strain. As such, when partially immune hosts are present at a high enough density, an immune escape variant can invade the viral population even though that variant cannot infect fully immune hosts. Invasion can occur even when the variant's immune escape mutation incurs a fitness cost, and we expect the expanding mutant population will evolve compensatory mutations that mitigate this cost. Thus the mutant lineage should replace the wild-type, and as immunity to it builds, the process will repeat. Our model provides a new explanation for the pattern of successive emergence and replacement of antigenic variants that has been observed in many respiratory viruses. We discuss our model relative to others for understanding the drivers of antigenic evolution in these and other respiratory viruses.

A speed limit on serial strain replacement from original antigenic sin

Many pathogens evolve to escape immunity, yet it remains difficult to predict whether immune pressure will lead to diversification, serial replacement of one variant by another, or more complex patterns. Pathogen strain dynamics are mediated by cross-protective immunity, whereby exposure to one strain partially protects against infection by antigenically diverged strains. There is growing evidence that this protection is influenced by early exposures, a phenomenon referred to as original antigenic sin (OAS) or imprinting. In this paper, we derive constraints on the emergence of the pattern of successive strain replacements demonstrated by influenza, SARS-CoV-2, seasonal coronaviruses, and other pathogens. We find that OAS implies that the limited diversity found with successive strain replacement can only be maintained if [Formula: see text] is less than a threshold set by the characteristic antigenic distances for cross-protection and for the creation of new immune memory. This bound implies a "speed limit" on the evolution of new strains and a minimum variance of the distribution of infecting strains in antigenic space at any time. To carry out this analysis, we develop a theoretical model of pathogen evolution in antigenic space that implements OAS by decoupling the antigenic distances required for protection from infection and strain-specific memory creation. Our results demonstrate that OAS can play an integral role in the emergence of strain structure from host immune dynamics, preventing highly transmissible pathogens from maintaining serial strain replacement without diversification.

A speed limit on serial strain replacement from original antigenic sin

Many pathogens evolve to escape immunity, yet it remains difficult to predict whether immune pressure will lead to diversification, serial replacement of one variant by another, or more complex patterns. Pathogen strain dynamics are mediated by cross-protective immunity, whereby exposure to one strain partially protects against infection by antigenically diverged strains. There is growing evidence that this protection is influenced by early exposures, a phenomenon referred to as original antigenic sin (OAS) or imprinting. In this paper, we derive new constraints on the emergence of the pattern of successive strain replacements demonstrated by influenza, SARS-CoV-2, seasonal coronaviruses, and other pathogens. We find that OAS implies that the limited diversity found with successive strain replacement can only be maintained if R 0 is less than a threshold set by the characteristic antigenic distances for cross-protection and for the creation of new immune memory. This bound implies a "speed limit" on the evolution of new strains and a minimum variance of the distribution of infecting strains in antigenic space at any time. To carry out this analysis, we develop a theoretical model of pathogen evolution in antigenic space that implements OAS by decoupling the antigenic distances required for protection from infection and strain-specific memory creation. Our results demonstrate that OAS can play an integral role in the emergence of strain structure from host immune dynamics, preventing highly transmissible pathogens from maintaining serial strain replacement without diversification.

Non-pharmaceutical interventions and the emergence of pathogen variants

Non-pharmaceutical interventions (NPIs), such as social distancing and contact tracing, are important public health measures that can reduce pathogen transmission. In addition to playing a crucial role in suppressing transmission, NPIs influence pathogen evolution by mediating mutation supply, restricting the availability of susceptible hosts, and altering the strength of selection for novel variants. Yet it is unclear how NPIs might affect the emergence of novel variants that are able to escape pre-existing immunity (partially or fully), are more transmissible or cause greater mortality. We analyse a stochastic two-strain epidemiological model to determine how the strength and timing of NPIs affect the emergence of variants with similar or contrasting life-history characteristics to the wild type. We show that, while stronger and timelier NPIs generally reduce the likelihood of variant emergence, it is possible for more transmissible variants with high cross-immunity to have a greater probability of emerging at intermediate levels of NPIs. This is because intermediate levels of NPIs allow an epidemic of the wild type that is neither too small (facilitating high mutation supply), nor too large (leaving a large pool of susceptible hosts), to prevent a novel variant from becoming established in the host population. However, since one cannot predict the characteristics of a variant, the best strategy to prevent emergence is likely to be an implementation of strong, timely NPIs.

Asynchrony between virus diversity and antibody selection limits influenza virus evolution

Seasonal influenza viruses create a persistent global disease burden by evolving to escape immunity induced by prior infections and vaccinations. New antigenic variants have a substantial selective advantage at the population level, but these variants are rarely selected within-host, even in previously immune individuals. Using a mathematical model, we show that the temporal asynchrony between within-host virus exponential growth and antibody-mediated selection could limit within-host antigenic evolution. If selection for new antigenic variants acts principally at the point of initial virus inoculation, where small virus populations encounter well-matched mucosal antibodies in previously-infected individuals, there can exist protection against reinfection that does not regularly produce observable new antigenic variants within individual infected hosts. Our results provide a theoretical explanation for how virus antigenic evolution can be highly selective at the global level but nearly neutral within-host. They also suggest new avenues for improving influenza control.

Dynamic interactions of influenza viruses in Hong Kong during 1998-2018

Influenza epidemics cause substantial morbidity and mortality every year worldwide. Currently, two influenza A subtypes, A(H1N1) and A(H3N2), and type B viruses co-circulate in humans and infection with one type/subtype could provide cross-protection against the others. However, it remains unclear how such ecologic competition via cross-immunity and antigenic mutations that allow immune escape impact influenza epidemic dynamics at the population level. Here we develop a comprehensive model-inference system and apply it to study the evolutionary and epidemiological dynamics of the three influenza types/subtypes in Hong Kong, a city of global public health significance for influenza epidemic and pandemic control. Utilizing long-term influenza surveillance data since 1998, we are able to estimate the strength of cross-immunity between each virus-pairs, the timing and frequency of punctuated changes in population immunity in response to antigenic mutations in influenza viruses, and key epidemiological parameters over the last 20 years including the 2009 pandemic. We find evidence of cross-immunity in all types/subtypes, with strongest cross-immunity from A(H1N1) against A(H3N2). Our results also suggest that A(H3N2) may undergo antigenic mutations in both summers and winters and thus monitoring the virus in both seasons may be important for vaccine development. Overall, our study reveals intricate epidemiological interactions and underscores the importance of simultaneous monitoring of population immunity, incidence rates, and viral genetic and antigenic changes.

Early prediction of antigenic transitions for influenza A/H3N2

Influenza A/H3N2 is a rapidly evolving virus which experiences major antigenic transitions every two to eight years. Anticipating the timing and outcome of transitions is critical to developing effective seasonal influenza vaccines. Using a published phylodynamic model of influenza transmission, we identified indicators of future evolutionary success for an emerging antigenic cluster and quantified fundamental trade-offs in our ability to make such predictions. The eventual fate of a new cluster depends on its initial epidemiological growth rate––which is a function of mutational load and population susceptibility to the cluster––along with the variance in growth rate across co-circulating viruses. Logistic regression can predict whether a cluster at 5% relative frequency will eventually succeed with ~80% sensitivity, providing up to eight months advance warning. As a cluster expands, the predictions improve while the lead-time for vaccine development and other interventions decreases. However, attempts to make comparable predictions from 12 years of empirical influenza surveillance data, which are far sparser and more coarse-grained, achieve only 56% sensitivity. By expanding influenza surveillance to obtain more granular estimates of the frequencies of and population-wide susceptibility to emerging viruses, we can better anticipate major antigenic transitions. This provides added incentives for accelerating the vaccine production cycle to reduce the lead time required for strain selection.

The efficacy of annual seasonal influenza vaccines depends on selecting the strain that best matches circulating viruses. This selection takes place 9–12 months prior to the influenza season. To advise this decision, we used an influenza A/H3N2 phylodynamic simulation to explore how reliably and how far in advance can we identify strains that will dominate future influenza seasons? What data should we collect to accelerate and improve the accuracy of such forecasts? And importantly, what is the gap between the theoretical limit of prediction and prediction based on current influenza surveillance? Our results suggest that even with detailed virological information, the tight race between the antigenic turnover dynamics and the vaccine development timeline limits early detection of emerging viruses. Predictions based on current influenza surveillance do not achieve the theoretical limit and thus our results provide impetus for denser sampling and the development of rapid methods for estimating viral fitness.

Explaining the geographical origins of seasonal influenza A (H3N2)

Most antigenically novel and evolutionarily successful strains of seasonal influenza A (H3N2) originate in East, South and Southeast Asia. To understand this pattern, we simulated the ecological and evolutionary dynamics of influenza in a host metapopulation representing the temperate north, tropics and temperate south. Although seasonality and air traffic are frequently used to explain global migratory patterns of influenza, we find that other factors may have a comparable or greater impact. Notably, a region's basic reproductive number (R0) strongly affects the antigenic evolution of its viral population and the probability that its strains will spread and fix globally: a 17-28% higher R0 in one region can explain the observed patterns. Seasonality, in contrast, increases the probability that a tropical (less seasonal) population will export evolutionarily successful strains but alone does not predict that these strains will be antigenically advanced. The relative sizes of different host populations, their birth and death rates, and the region in which H3N2 first appears affect influenza's phylogeography in different but relatively minor ways. These results suggest general principles that dictate the spatial dynamics of antigenically evolving pathogens and offer predictions for how changes in human ecology might affect influenza evolution.

Seasonality in risk of pandemic influenza emergence

Influenza pandemics can emerge unexpectedly and wreak global devastation. However, each of the six pandemics since 1889 emerged in the Northern Hemisphere just after the flu season, suggesting that pandemic timing may be predictable. Using a stochastic model fit to seasonal flu surveillance data from the United States, we find that seasonal flu leaves a transient wake of heterosubtypic immunity that impedes the emergence of novel flu viruses. This refractory period provides a simple explanation for not only the spring-summer timing of historical pandemics, but also early increases in pandemic severity and multiple waves of transmission. Thus, pandemic risk may be seasonal and predictable, with the accuracy of pre-pandemic and real-time risk assessments hinging on reliable seasonal influenza surveillance and precise estimates of the breadth and duration of heterosubtypic immunity.

Self-enforcing regional vaccination agreements

In a highly interconnected world, immunizing infections are a transboundary problem, and their control and elimination require international cooperation and coordination. In the absence of a global or regional body that can impose a universal vaccination strategy, each individual country sets its own strategy. Mobility of populations across borders can promote free-riding, because a country can benefit from the vaccination efforts of its neighbours, which can result in vaccination coverage lower than the global optimum. Here we explore whether voluntary coalitions that reward countries that join by cooperatively increasing vaccination coverage can solve this problem. We use dynamic epidemiological models embedded in a game-theoretic framework in order to identify conditions in which coalitions are self-enforcing and therefore stable, and thus successful at promoting a cooperative vaccination strategy. We find that countries can achieve significantly greater vaccination coverage at a lower cost by forming coalitions than when acting independently, provided a coalition has the tools to deter free-riding. Furthermore, when economically or epidemiologically asymmetric countries form coalitions, realized coverage is regionally more consistent than in the absence of coalitions.

Prediction, dynamics, and visualization of antigenic phenotypes of seasonal influenza viruses

Human seasonal influenza viruses evolve rapidly, enabling the virus population to evade immunity and reinfect previously infected individuals. Antigenic properties are largely determined by the surface glycoprotein hemagglutinin (HA), and amino acid substitutions at exposed epitope sites in HA mediate loss of recognition by antibodies. Here, we show that antigenic differences measured through serological assay data are well described by a sum of antigenic changes along the path connecting viruses in a phylogenetic tree. This mapping onto the tree allows prediction of antigenicity from HA sequence data alone. The mapping can further be used to make predictions about the makeup of the future A(H3N2) seasonal influenza virus population, and we compare predictions between models with serological and sequence data. To make timely model output readily available, we developed a web browser-based application that visualizes antigenic data on a continuously updated phylogeny.

Modelling the role of immunity in reversion of viral antigenic sites

Antigenic sites in viral pathogens exhibit distinctive evolutionary dynamics due to their role in evading recognition by host immunity. Antigenic selection is known to drive higher rates of non-synonymous substitution; less well understood is why differences are observed between viruses in their propensity to mutate to a novel or previously encountered amino acid. Here, we present a model to explain patterns of antigenic reversion and forward substitution in terms of the epidemiological and molecular processes of the viral population. We develop an analytical three-strain model and extend the analysis to a multi-site model to predict characteristics of observed sequence samples. Our model provides insight into how the balance between selection to escape immunity and to maintain viability is affected by the rate of mutational input. We also show that while low probabilities of reversion may be due to either a low cost of immune escape or slowly decaying host immunity, these two scenarios can be differentiated by the frequency patterns at antigenic sites. Comparison between frequency patterns of human influenza A (H3N2) and human RSV-A suggests that the increased rates of antigenic reversion in RSV-A is due to faster decaying immunity and not higher costs of escape.

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.

Environmental dependence of genetic constraint

The epistatic interactions that underlie evolutionary constraint have mainly been studied for constant external conditions. However, environmental changes may modulate epistasis and hence affect genetic constraints. Here we investigate genetic constraints in the adaptive evolution of a novel regulatory function in variable environments, using the lac repressor, LacI, as a model system. We have systematically reconstructed mutational trajectories from wild type LacI to three different variants that each exhibit an inverse response to the inducing ligand IPTG, and analyzed the higher-order interactions between genetic and environmental changes. We find epistasis to depend strongly on the environment. As a result, mutational steps essential to inversion but inaccessible by positive selection in one environment, become accessible in another. We present a graphical method to analyze the observed complex higher-order interactions between multiple mutations and environmental change, and show how the interactions can be explained by a combination of mutational effects on allostery and thermodynamic stability. This dependency of genetic constraint on the environment should fundamentally affect evolutionary dynamics and affects the interpretation of phylogenetic data.

A deterministic model for influenza infection with multiple strains and antigenic drift

We describe a multiple strain Susceptible Infected Recovered deterministic model for the spread of an influenza subtype within a population. The model incorporates appearance of new strains due to antigenic drift, and partial immunity to reinfection with related circulating strains. It also includes optional seasonal forcing of the transmission rate of the virus, which allows for comparison between temperate zones and the tropics. Our model is capable of reproducing observed qualitative patterns such as the overall annual outbreaks in the temperate region, a reduced magnitude and an increased frequency of outbreaks in the tropics, and the herald wave phenomenon. Our approach to modelling antigenic drift is novel and further modifications of this model may help improve the understanding of complex influenza dynamics.

The antigenic evolution of influenza: drift or thrift?

It is commonly assumed that antibody responses against the influenza virus are polarized in the following manner: strong antibody responses are directed at highly variable antigenic epitopes, which consequently undergo 'antigenic drift', while weak antibody responses develop against conserved epitopes. As the highly variable epitopes are in a constant state of flux, current antibody-based vaccine strategies are focused on the conserved epitopes in the expectation that they will provide some level of clinical protection after appropriate boosting. Here, we use a theoretical model to suggest the existence of epitopes of low variability, which elicit a high degree of both clinical and transmission-blocking immunity. We show that several epidemiological features of influenza and its serological and molecular profiles are consistent with this model of 'antigenic thrift', and that identifying the protective epitopes of low variability predicted by this model could offer a more viable alternative to regularly update the influenza vaccine than exploiting responses to weakly immunogenic conserved regions.

Canalization of the evolutionary trajectory of the human influenza virus

BACKGROUND

Since its emergence in 1968, influenza A (H3N2) has evolved extensively in genotype and antigenic phenotype. However, despite strong pressure to evolve away from human immunity and to diversify in antigenic phenotype, H3N2 influenza shows paradoxically limited genetic and antigenic diversity present at any one time. Here, we propose a simple model of antigenic evolution in the influenza virus that accounts for this apparent discrepancy.

RESULTS

In this model, antigenic phenotype is represented by a N-dimensional vector, and virus mutations perturb phenotype within this continuous Euclidean space. We implement this model in a large-scale individual-based simulation, and in doing so, we find a remarkable correspondence between model behavior and observed influenza dynamics. This model displays rapid evolution but low standing diversity and simultaneously accounts for the epidemiological, genetic, antigenic, and geographical patterns displayed by the virus. We find that evolution away from existing human immunity results in rapid population turnover in the influenza virus and that this population turnover occurs primarily along a single antigenic axis.

CONCLUSIONS

Selective dynamics induce a canalized evolutionary trajectory, in which the evolutionary fate of the influenza population is surprisingly repeatable. In the model, the influenza population shows a 1- to 2-year timescale of repeatability, suggesting a window in which evolutionary dynamics could be, in theory, predictable.