Skip to content


You're viewing the new version of our site. Please leave us feedback.

Learn more

BMC Infectious Diseases

Open Access
Open Peer Review

This article has Open Peer Review reports available.

How does Open Peer Review work?

Variation in loss of immunity shapes influenza epidemics and the impact of vaccination

BMC Infectious DiseasesBMC series – open, inclusive and trusted201717:632

Received: 3 November 2016

Accepted: 5 September 2017

Published: 19 September 2017



Protective antibody immunity against the influenza A virus wanes in 2–7 years due to antigenic drift of the virus’ surface proteins. The duration of immune protection is highly variable because antigenic evolution of the virus is irregular. Currently, the variable nature of the duration of immunity has had little attention in analyses of the impact of vaccination, including cost-effectiveness studies.


We developed a range of mathematical transmission models to investigate the effect of variable duration of immunity on the size of seasonal epidemics. The models range from simple conceptual to more realistic, by distinguishing between infection- versus vaccination-induced immunity, by inclusion of primary vaccine failure, by assuming a leaky vaccine, and by the inclusion of age-dependent contact patterns.


We show that annual variation in the duration of immunity causes large variation in the size of epidemics, and affects the effectiveness of vaccination. Accumulation of susceptible individuals in one or more mild seasons results in a disproportionately large outbreak in a subsequent season. Importantly, variation in the duration of immunity increases the average infection attack rate when the vaccination coverage is around the outbreak threshold. Specifically, in a tailored age-stratified model with a realistic reproduction number (R 0 = 1.4) and vaccination coverage of 25%, we find that the attack rate in unvaccinated children (<10 years old) is negligible if the duration of immunity is constant, while on average 2.8% (2.5–97.5% percentiles: 1.8–4.1%) of the children are infected if the duration of immunity is variable. These findings stem from the buildup of susceptibility over multiple seasons by waning of immunity, and the nonlinear relation between susceptibility and infection attack rates.


The models illustrate that variation in the duration of immunity impacts the long-term effectiveness of vaccination, and that vaccine effectiveness cannot be judged for each year in isolation. Our findings have implications for vaccination strategies that aim to maximize the vaccination coverage while extending the age range of persons eligible for vaccination.


The influenza virus is responsible for a substantial health burden with annual infection attack rates of 5–10% in adults and 20–30% in children [1]. Annual vaccination can protect against influenza infection and therefore vaccination is recommended for elderly and other individuals at high-risk of developing serious disease [24]. Despite efforts to improve vaccines against influenza, vaccination still provides only partial protection, and antibodies generated after infection or vaccination protect a person only for a limited number of years. This is because the evolution of influenza’s surface proteins prevents neutralization of evolved viruses by preexisting antibodies. This so-called antigenic drift is irregular [5], with relative conserved periods that last 2–7 years [68], and is only modestly predictable [9, 10]. The variable duration of immunity is a cause of substantial year-to-year differences of estimates of vaccine efficacy and vaccine effectiveness [1116], illustrating the need for robust models to assess vaccination policies.

Most epidemiological and cost-effectiveness analyses of influenza vaccination assume that the rate at which immunity is lost does not vary between years [1733]. In this were the case, each influenza season could be analyzed in isolation. In the more realistic scenario of a variable rate at which immunity is lost due to antigenic drift, the fraction of the population that is susceptible (henceforth susceptibility of the population) will vary between years. In this case, it is necessary to track the number of infections that happened in previous seasons, and that still give protection against the currently circulating virus. The importance of linking the infection history of successive seasons is recognized [34, 35], but difficult to implement in prospective studies [36], since the level of preexisting immunity in the population as well as the antigenic evolution of the virus are hard to predict.

Here we study the impact of variation in the duration of immunity on the infection attack rates and the epidemic peaks in a disease transmission framework. To focus on the key elements of acquired immunity, vaccination, and influenza evolution, we formulate an idealized model in which the susceptibility in each season depends on the cumulative effect of all infections and vaccinations in the past. We show that variation in the duration of immunity results in alternating mild and severe seasons, and in an increase in the average infection attack rate and peak prevalence. Especially around the critical vaccination coverage, we find that variation in the duration of immunity increases the number of infections. We also show that these findings are not restricted to highly simplified models, but are also observed in more realistic influenza transmission models.



For each season we calculate population susceptibility using the accumulation of influenza infections and influenza vaccinations over previous years [37]. Susceptibility increases due to waning of immunity and demographic turnover, and decreases by vaccination and infection. Throughout, we focus on populations in temperate zones where epidemics generally last for 4–12 weeks, and where there is hardly influenza circulation in the remainder of the year [38]. In our models, vaccination takes place just before the epidemic, while loss of immunity occurs after the epidemic when there is no virus circulation. Demographic turnover is not separately modeled, as births increase the number of susceptible persons in a similar fashion as loss of immunity does if population size is constant. We map susceptibility at the start of the season to the attack rate using the reproduction number, which is defined as the average number of infections caused by one infected individual in a fully susceptible population (as described in more detail below). Throughout, we compare scenarios with constant and variable duration of immunity. In the variable immunity model, the fractions of the immune individuals that stays immune in the next season are drawn from a Beta distribution that mimics the punctuated antigenic drift of the H3N2 subtype observed using hemagglutination inhibition tests [6, 7].

Below we formulate a basic model to show that a variable duration of immunity increases the epidemic peaks and infection attack rates. We then adapt the basic model so that vaccination provides protection for one season. Finally, we consider three extensions with (i) a ‘leaky’ vaccine, (ii) an age-structured population model with age-specific contact patterns, and (iii) realistic latent and infectious period distributions in an age-structured population with leaky vaccination [35]. The extensions show that the phenomena observed in the basic model are robust, and do not depend on specific model assumptions.

Model structure

In our basic model, individuals are either fully susceptible to infection or fully immune, either due to natural infection or vaccination. The fraction of susceptible individuals (susceptibility) is determined by the cumulated immunity during outbreaks and vaccination in past seasons. Three processes link the susceptibility s t+1 in season t+1 to the susceptibility s t in the previous season, i) vaccination, ii) the yearly influenza epidemic, and iii) the loss of immunity during the inter-epidemic period by virus evolution and demographic turnover. As these processes largely take place sequentially and on different timescales, we model these processes sequentially. In the model, vaccination is applied before the epidemic, and loss of immunity takes place after the epidemic in the inter-epidemic period. Schematically, the model looks as follows:
$$ s_{t} \xrightarrow{\text{vaccination}} s{\prime}_{t} \xrightarrow{\text{epidemic}} s{\prime}{\prime}_{t} \xrightarrow{\text{waning immunity}} s_{t+1}~. $$
In the first step, a fraction of the population is vaccinated, such that the fraction remaining susceptible after vaccination is given by,
$$ s{\prime}_{t}=\left(1-v_{c}v_{e}\right)s_{t}, $$

where v c and v e are the vaccination coverage (the proportion of individuals who are vaccinated) and vaccine efficacy (the proportion of vaccinated individuals that is protected against infection), respectively. Thus, a fraction v c v e is protected against infection. Throughout, we vary v c and take v e =0.5. To keep the model simple and make our arguments as transparent as possible, we assume that each year vaccination is given to a random set of individuals.

Next, part of the population is infected during the epidemic, and will be immune afterward:
$$s{\prime}{\prime}_{t}=s{\prime}_{t}-z_{t}, $$
where z t :=z t (s) is the attack rate (the fraction infected in season t). In a fully susceptible population (s=1) the attack rate is given by z t =1− exp(−z t R 0), with R 0 representing the basic reproduction number, defined as the expected number of infections caused by one infected individual in a completely susceptible population. Estimates of the reproduction number for influenza (R 0) range between 1-2 [38, 39]. Here, we take R 0=1.4 as default value.
Due to immunity acquired in previous years, only a fraction s is susceptible, such that
$$ z_{t}=s{\prime}(1-\exp(- z_{t} R_{0})). $$
In the last step, immunity is lost, so that
$$s_{t+1}=1-\gamma_{t}(1-s{\prime}{\prime}_{t}), $$
where γ t is the fraction of the population that retains its immunity. Loss of immunity is caused by the reduction of cross-protection between the current strain and strains circulating previous seasons. Throughout, we explore the impact of annual variation in the rate at which immunity is lost due to antigenic drift of the virus (i.e. the effect of variability in γ t ). In the deterministic scenario we use \(\gamma =\frac {5}{7}\), and in the stochastic scenario we draw γ from a realistic Beta(5,2) distribution that has a mean of 5/7. In this manner, the deterministic scenario arises naturally as a special case of the stochastic scenario when the parameters of the beta distribution tend to infinity. In both scenarios, immunity lasts 3.5 years on average in both scenarios.

Combining the above, we obtain the following set of equations that map the susceptibility s t in year t to the susceptibility s t+1 in year t+1:

$$\begin{array}{*{20}l} z_{t} &= \left(1-v_{c}v_{e}\right) s_{t} \left(1-\exp\left(- z_{t} R_{0}\right)\right), \end{array} $$
$$\begin{array}{*{20}l} s_{t+1} &= 1-\gamma_{t} \left(1-\left(1-v_{c}v_{e}\right) s_{t} + z_{t}\right). \end{array} $$

These equations can be solved numerically.

The map in Eq. 3 is used to calculate the infection attack rate, but does not yield the infection prevalence over time (the epidemic curve). For this, we use the dynamical system corresponding to the above model. The ordinary differential equations (ODEs) for the fractions of the population in the susceptible (S), infected and infectious (I), and recovered (R) compartments during the epidemic are given by
$$\begin{array}{*{20}l} \frac{dS}{dt} & = - \beta S I, \\ \frac{dI}{dt} & = \beta S I - \nu I, \\ \frac{dR}{dt} & = \nu I. \end{array} $$
Throughout, we fix ν and set β=ν R 0, such that the infection attack rate is independent of ν. In Fig. 1 we take ν=2 (day−1) and S(0)=s′, I(0)=0.000001, and R(0)=1−S(0)−I(0), where s′ is the susceptibility after vaccination as given in Eq. 1. The differential equations are solved with a standard ODE solver using Python.
Fig. 1

Variation in the duration of immunity increases the height of epidemic peaks. In the absence of vaccination, a pattern of regular epidemics is observed if the duration of immunity is constant (blue line), while an irregular pattern with alternating small and large epidemics is observed if the duration of immunity is variable (red line). After the introduction of vaccination (vertical gray line), the pattern of regular epidemics persists in the constant immune duration scenario, albeit at a lower level (blue line). If the duration of immunity is variable, the irregular pre-vaccination pattern of alternating small and large epidemics is exacerbated (red line). The mean of the prevalence peaks is substantially higher in the variable immunity scenario than in the constant immune duration scenario (dotted lines). The means of the prevalence peaks are calculated over 1,000 seasons using a burn-in period of 50 years. See Methods for details

Short-term protection of vaccination

Natural immunity (i.e. immunity acquired after infection) is likely to last longer than immunity provided by vaccination. In the basic model presented above (Eq. 3), we assume that vaccination gives the same duration of protection as natural infection. Here we present the other extreme that vaccination protects individuals for one season only. In this case the the model equations are given by
$$z_{t} = \left(1-v_{c}v_{e}\right) s_{t} \left(1-\exp\left(- z_{t} R_{0}\right)\right), $$
$$s_{t+1} = 1-\gamma_{t} \left(1 - s_{t} + z_{t}\right), $$
where only the attack rate depends on vaccination.

Leaky vaccine model

It is known that influenza vaccines provide partial protection against infection and disease [1214]. We therefore extend the basic model to include the possibility that the vaccine gives partial protection while infection still gives full protection. In this model, protection after infection and vaccination are assumed to wane at identical rates, while a vaccinated individual who is subsequently infected gains full protection.

In an all-or-nothing vaccine model, v e is the proportion of vaccinated individuals that is protected, while with a leaky vaccine, we define v r to be the probability that vaccinated individuals are protected in a single exposure that would have led to transmission to an unvaccinated individual. Unless stated otherwise, we take v e =1.0, and v r =0.5.

In the leaky vaccine model we need to track the fractions of the population that are susceptible (s) and vaccinated (p, partially immune). Following the same steps as before, we obtain
$$\begin{array}{*{20}l} {}s_{t} \!\rightarrow\! s{\prime}_{t}&\,=\,(\!1\,-\,v_{c} v_{e})s_{t} \!\rightarrow \!s{\prime}{\prime}_{t}\!=s{\prime}_{t} \,-\, z^{s}_{t} \!\rightarrow \!s_{t+1}\,=\,1-\gamma_{t}(1-s{\prime}{\prime}_{t}),\\ {}p_{t} \!\rightarrow \!p{\prime}_{t}&\,=\,p_{t}\,+\,v_{c} v_{e} s_{t} \rightarrow\! \!p{\prime}{\prime}_{t}\,=\,p{\prime}_{t} \,-\, z^{p}_{t}\! \rightarrow\!\! p_{t+1}=\gamma_{t} p{\prime}{\prime}_{t}~, \end{array} $$

where \(z^{s}_{t}\) and \(z^{p}_{t}\) are the attack rates of susceptible and vaccinated individuals, and \(z_{t}=z^{s}_{t}+z^{p}_{t}\) is the overall attack rate. As in the basic model, each year a fraction 1−γ t of the vaccinated and partially immune individuals lose immunity.

Combining the above, we obtain the following map linking s t and p t to s t+1 and p t+1:
$$\begin{array}{*{20}l} s_{t+1} & = 1-\gamma_{t} \left(1 - \left(1-v_{c} v_{e}\right)s_{t} + z^{s}_{t} \right)\\ p_{t+1}&=\gamma_{t} \left(p_{t} + v_{c} v_{e} s_{t}-z_{t}^{p} \right) \\ z^{s}_{t} & = s'_{t} \left(1 - \exp\left(-R_{0} \left(z^{s}_{t} + z^{p}_{t}\right)\right)\right)\\ z^{p}_{t} & = p'_{t} \left(1 - \exp\left(-R_{0} \left(1-v_{r}\right) \left(z^{s}_{t} + z^{p}_{t}\right) \right)\right)~. \end{array} $$

As before, the above equations are readily solved numerically.

In a similar vein as before, we obtain a system of ODEs for the fractions of the population in the susceptible (S), partially protected (P), infected and infectious (I), and protected (R) compartments:
$$\begin{array}{*{20}l} \frac{dS}{dt} & = - \beta S I \\ \frac{dP}{dt} & = - \beta (1-v_{r}) P I \\ \frac{dI}{dt} & = \beta (S + (1-v_{r}) P) I - \nu I \\ \frac{dR}{dt} & = \nu I~, \end{array} $$

with initial conditions S(0)=s , P(0)=p , I(0)=0.000001, and R(0)=1−S(0)−P(0)−I(0).

Age-stratified model

We extend the main model (Eq. 3) in another direction to include age structure. We consider three age classes: children (≤ 10 years), adults (> 10 and ≤ 60 years), and elderly (> 60 years). These broad classes are chosen to reflect the fact that young children are the drivers of transmission, while serious disease occurs mostly in elderly. Moreover, in many countries healthy elderly are eligible for vaccination (usually > 60 or > 65 years). The attack rates, with i {children, adults, elderly}, are given by
$$ z_{t,i} = s{\prime}_{t,i} \left(1 - \exp \left(- R_{0} \sum_{j} g_{ij} z_{t,j} \right) \right). $$
In the above, the next generation matrix G with elements g ij is given by
$$g_{ij}=q f_{j} c_{ij}, $$
where c ij =c ji is the (symmetric) contact rate between individuals in age groups i and j, and f j is the proportion of the population in age group j [40]. The proportionality parameter q is used to scale G such that the dominant eigenvalue equals 1 and R 0 in Eq. 5 corresponds to the basic reproduction number.
For simplicity, we assume a uniform population demography and lifespan of 80 years, such that f=(1/8,5/8,2/8). Contacts between children are intense, while children and elderly have fewer contacts. In our analyses we take the following contact matrix
$$\boldsymbol{C}=\left(\begin{array}{ccc} 9 & 3 & 1\\ 3 & 3 & 2 \\ 1 & 2 & 2 \end{array}\right)~, $$
which is roughly in agreement with the contact structure as observed in representative contact survey studies in the Netherlands [40], and elsewhere [41].
The susceptibilities s t+1,i in the next season are given by
$$\begin{array}{*{20}l} s_{t+1,i} & = 1-\gamma_{t} \left(1 - \left(1-v_{c}v_{e}\right) s_{t,i} + z_{t,i} \right), \end{array} $$

with z t,i as in Eq. 5.

To calculate the peak prevalence in the age-stratified model (Eq. 5), we consider the age-stratified SIR-model with next-generation matrix G, i.e.
$$\begin{array}{*{20}l} \frac{dS_{i}}{dt} & = - \beta S_{i} \sum_{j} g_{ij} I_{j} \\ \frac{dI_{i}}{dt} & = \beta S_{i} \sum_{j} g_{ij} I_{j} - \nu I_{i} \\ \frac{dR_{i}}{dt} & = \nu I_{i}, \end{array} $$

where i,j {children, adults, elderly}.

Tailored influenza transmission model

As a final extension, we consider a tailored transmission model that incorporates both age-dependent contact patterns and a partially immune compartment. In addition, we include compartments for latently infected infected individuals (i.e. those who are infected but not yet infectious), and assume that both the latent and infectious periods are gamma (Erlang) distributed. Such dynamical systems form the backbone of many influenza transmission models [35, 42]. Denoting the forces of infection by λ i (λ i {children, adults, elderly}), the model dynamics during the epidemic are specified by the following system of ODEs:
$$\begin{array}{*{20}l} \frac{dS_{i}}{dt} & = - \lambda_{i} S_{i}, \\ \frac{dP_i}{dt} &= -(1-v_r) \lambda_{i} P_{i},\\ \\ \frac{dE^{(1)}_{i}}{dt} & = \lambda_{i} (S_{i} + (1-v_{r})P_{i}) - \nu E^{(1)}_{i},\\ \frac{dE^{(2)}_{i}}{dt} & = \nu E^{(1)}_{i} - \nu E^{(2)}_{i},\\ \\ \frac{dI^{(1)}_{i}}{dt} & = \nu E^{(2)}_{i} - \mu I^{(1)}_{i},\\ \frac{dI^{(2)}_{i}}{dt} & = \mu I^{(1)}_{i} - \mu I^{(2)}_{i},\\ \\ \frac{dR_{i}}{dt} & = \mu I^{(2)}_{i}, \end{array} $$
with forces of infection λ i given by,
$$\lambda_{i} = \beta \sum_{j} g_{ij} \left(I^{(1)}_{j} + I^{(2)}_{j}\right), $$
and i,j {children, adults, elderly}. Following [35], we take ν=2.5 (day −1) and μ=1.1 (day −1). Furthermore, we take β=μ R 0/2, and other parameters are as specified earlier.


Variation in the duration of immunity increases the height of epidemic peaks

To investigate the impact of variation in the duration of immunity on the epidemic dynamics, we compare scenarios with and without variation in the duration of immunity. For the moment, we consider the case without vaccination. If the duration of immunity is constant, the number of individuals who are infected over an epidemic is balanced exactly by the number of individuals added to the susceptible pool during the inter-epidemic period by waning of immunity. As a consequence, we observe regular annual epidemics in this scenario (Fig. 1).

When we take the more realistic assumption that the duration of immunity is variable, and vary the fraction of the population that loses its protection in each year, the situation is different. In this case, years with no or hardly any influenza activity are interspersed with years of large influenza epidemics (Fig. 1). Susceptible individuals who escape infection in mild seasons (i.e. seasons in which by chance only a small fraction of the population loses immunity) continue to be susceptible in the next season, resulting in an accumulated pool of susceptible individuals. Importantly, the mean peak of the prevalence is larger in the variable case, while on average the same fraction of the population becomes susceptible by losing protection against the virus each season. The disproportionally large outbreaks are the result of indirect effects in the infection dynamics, namely that each additional infection increases the exposure and the probability of infection of susceptible individuals. In other words, a positive feedback of infections exists during the epidemics that cause many infections in short time when the pool of susceptible individuals is large.

Variation in the duration of immunity reduces the impact of vaccination

Next, we incorporate vaccination into the model. We consider two types of individuals (susceptible and immune to infection), and assume a vaccination coverage (i.e. per person probability of vaccination) of 10% and vaccine efficacy of 0.5 (Methods). Hence, a vaccinated person is fully immune with probability 0.5 and fully susceptible with probability 0.5, i.e. the probability of primary vaccine failure is 0.5. This assumption is relaxed later by adding partial protection to the model.

Even in this idealized context, we find large variation in the size of epidemics in the presence of vaccination, due to variation in the duration of immunity (Fig. 1). Although the mean of the epidemic peaks is lower post-vaccination than pre-vaccination, and large outbreaks occur less frequently, vaccination cannot prevent the irregular pattern of years with small or no outbreaks interspersed with occasional large outbreaks. In contrast, if the duration of protection is constant, vaccination has a substantial protective effect in the first years after its introduction (the honeymoon period [43]), and a regular pattern of (much smaller) yearly epidemics reappears in subsequent seasons.

Comparing the average attack rates in the scenarios with constant and variable duration of immunity, we find a significant difference in the post-vaccination era (Fig. 2 b). No such difference is observed without vaccination (Fig. 2 a). This is surprising, as on average immunity is lost at identical rates in both scenarios. Here, individuals who would have escaped from infection in the scenario with constant duration of immunity can be infected when the duration of immunity is short in the scenario with variable duration of immunity.
Fig. 2

Attack rates are increased in models with variable duration of immunity. A comparison of the attack rates with a constant (blue) and a variable (red) duration of immunity without vaccination (a) and with vaccination (with a coverage of 20%) (b). With vaccination the attack rates are significantly higher in the scenario with variable duration of immunity, while in both cases the average duration of immunity is identical. Shown are median values of 1,000 simulations. Black error bars denote the 2.5−97.5 percentiles. For each scenario, the mean infection attack rate is calculated over a period of 100 years after a burn-in period of 50 years

Elevated attack rates near the critical vaccination coverage

For influenza, the vaccination coverage varies substantially over time and between countries/regions. Therefore, we investigate how the probability of an epidemic, the infection attack rate, and peak prevalence are molded by vaccine uptake (Fig. 3). As expected, in the scenario with constant duration of immunity there is a sharp transition at a critical vaccination coverage: yearly epidemics occur when the vaccination coverage is below the outbreak threshold, and no epidemics occur when the vaccination coverage is above the threshold. In the latter case, the increase in the number of immune individuals by vaccination is such that it is enough to prevent a critical accretion of susceptible individuals by waning of immunity and demographic turnover (Fig. 3 a). This does not hold anymore in the scenario with variable immunity. In this scenario, epidemics may occur even when the vaccination coverage is high, and the influx of susceptible individuals by waning immunity and demographic turnover is on average more than counterbalanced by vaccination. At low vaccination coverage no outbreaks occur in some of the years, but the average infection attack rates are comparable with those obtained in the scenario with constant immunity. In a similar vein as before (Fig. 1), we find that the peak prevalence is increased in the case of variable duration of immunity at any vaccination coverage (Fig. 3 b). Here the (absolute) difference is most pronounced if the vaccination coverage is low.
Fig. 3

Infection attack rates and epidemic peaks are increased in models with variable duration of immunity. Shown are (a) the yearly probability that an epidemic will occur, (b) the peak prevalence in an epidemic, and (c) the infection attack rate, all as a function of the vaccination coverage. The figures show these quantities with constant (blue lines) and variable duration of immunity (red lines). Notice that if the duration of immunity is variable, the infection attack rate and the epidemic peak are always higher than if the duration of immunity is constant. Also notice that the difference between the infection attack rates in both scenarios is most pronounced at a vaccination coverage around the outbreak threshold (approximately at a vaccination coverage of 21%). The red line and the shaded area indicate the median and the 2.5-97.5 percentiles of 1,000 simulations, respectively. Each simulation lasts 150 years. Mean values are calculated over the last 100 years

A further comparison of scenarios with variable and constant immunity shows that the attack rate is higher in the scenario with variable immunity, especially when the vaccination coverage is close to the outbreak threshold (Fig. 3 c). Here the accumulated number of susceptible individuals determines to a large extent whether an outbreak occurs, and the attack rates do not depend solely on the loss of immunity since the last epidemic. This is illustrated in Fig. 4, which shows that the attack rate correlates not only positively with the susceptible fraction at the start of the current season, but also with the susceptible fraction at the start of the previous season (Fig. 4 a). In this regime, i.e. above a vaccination coverage of 15% (Fig. 4 b), epidemics occur if enough susceptible individuals have accumulated over multiple seasons. Notice that in the absence of vaccination, the correlation between the infection attack rate and susceptibility in the previous season is negative (Fig. 4 a), because in this case a high susceptibility at the start of the influenza season is likely to result in an outbreak in the same season, making the occurrence of an outbreak in the next season less likely.
Fig. 4

Correlations between the infection attack rate and susceptibility in the current and previous year. (a) The Pearson correlation between the number of susceptible individuals at the start of one influenza season and the infection attack rate in the next influenza season is negative without vaccination (\(\hat {\rho }=-0.24\); p < 0.00001), as a high susceptibility in one year implies a high probability of an epidemic in that year, and a low probability in the next year. At higher vaccination coverage, however, the accumulation of susceptible persons over multiple years becomes more important, and the correlation between susceptibility in one year and the infection attack rate in the next year is reversed (\(\hat {\rho }=0.13\); p < 0.00001). The switch in correlation from negative to positive is seen in (b) as a function of vaccination coverage (blue line). The correlation between the number of susceptible individuals at the start of the influenza season and the infection attack rate in the same year is positive, with correlations decreasing with increasing vaccination coverage (cyan line). Parameters are as in Fig. 3; in the right panel of (a) a vaccination coverage of 25% is used. For each vaccination coverage, we analyze time series of 10,000 years

Thus, the analyses show that a variable duration of immunity changes the epidemiological dynamics in two ways. First, epidemics cease to occur every year, while the size of epidemic peaks becomes highly variable and on average higher than in the constant immunity scenario. Second, the accumulation of susceptible individuals over multiple seasons in the variable immunity scenario is able to lift susceptibility over the outbreak threshold where this would not be possible in the corresponding scenario with fixed duration of immunity. Notice that these phenomena only occur in a multi-season framework, and are a consequence of the fact that the final size function is convex near the outbreak threshold (Additional file 1: Figure S1). We argue that these findings hold in general, and do not depend on specific modeling assumptions (Additional file 1: Text S1).

The phenomenon of increased attack rates in models with variable duration of immunity disappears in models in which vaccination induces protection for at most one season. Such scenarios may not be wholly unrealistic [44]. In this case, the mean epidemic infection attack rates are identical in the models with fixed and variable duration of immunity over the full range of vaccination coverages (Additional file 1: Figure S3c), the reason being that such a large fraction of the population is vaccinated around the critical vaccination coverage that no matter how much immunity is lost, an outbreak remains unlikely.

Tailored epidemic models

More detailed models are commonly used to study influenza transmission dynamics [18, 20, 23, 33, 35]. We therefore investigate an epidemic model with additional layers of complexity to study the robustness of the above results (Methods). First, we assume Erlang (Gamma) distributed latent and infectious periods [35]. Second, we stratify the population into three age classes: children (less than 10 years old), adults (10–60 years), and elderly (older than 60 years). This extension captures the essence of heterogeneity in contact patterns by age [4042], in particular the observation that contacts are intense among children, intermediate among adults and between adults and children, and low with and among the elderly. To be able to compare results with the earlier models without age-structure, we assume indiscriminate vaccination, i.e. vaccination coverage is equal in all age groups. Third, since influenza vaccines provide partial protection to infection [1214], we extend the model to account for the leakiness of the vaccine. With this assumption, heavily exposed vaccinated individuals are more likely to be infected than vaccinated individuals who are sporadically exposed.

In the tailored model, we find that the epidemic peak as well as the attack rates are still increased when the duration of immunity is variable (Fig. 5). This is true for all age groups, and for the unvaccinated as well as vaccinated subpopulations. In line with the heterogeneity in contact patterns, and in agreement with incidence data [1], we find that the attack rates are highest in children, followed by adults and the elderly. When the vaccination coverage is such that the population is close to the outbreak threshold in the constant immunity scenario (20–30%), attack rates are higher with than without a variable duration of immunity, like we observed in the basic model (Fig. 3). Specifically, at a vaccination coverage of 25% the attack rate in unvaccinated children is only 0.5% if the duration of immunity is constant, while it is 2.8% (95% CI: 1.8%, 4.1%) if the duration of immunity is variable (Fig. 5). Similarly, at 25% vaccination coverage only 0.3% of the vaccinated children are infected in the model with a constant duration of immunity, while 1.5% (95% CI: 1.0%, 2.2%) of the vaccinated children are infected if the duration of immunity is variable. Similar results are obtained when the basic model is extended with either a leaky vaccine, or an age-dependent transmission framework (Additional file 1: Figure S4). Together, the extensions of the basic model illustrate that the finding of increased attack rates due to a variable duration of immunity is generic, and is expected to be found in even more complex models.
Fig. 5

Increased attack rates in an age-specific epidemic model. Variation in duration of immunity increases the epidemic peak (a) and the attack rates in the unvaccinated (b) as well as in the vaccinated (c) subpopulations. In this model we assume a leaky vaccine, use age-specific contact patterns, and take realistic latent and infectious periods (using two latent compartments and two infectious compartments). The time series shows periodic behavior with periods of 2 and 3 years (at vaccination coverages 0-13% and 14-15%, respectively) for the scenario with constant duration of immunity. Shown are the averages of the periodic epidemics, explaining the irregularity of the blue curves. Percentile bands are not drawn in b and c


The effectiveness of humoral immune responses developed during past influenza infections to current strains is determined by the antigenic evolution of the virus [68]. Typically, the effective duration of immunity is relatively short but longer than a single influenza season, so that past influenza epidemics influence the epidemiology of the present epidemic. It is also known that the period of immunity varies substantially, and is driven by antigenic evolution of the virus. Here we have shown that linking the year-to-year epidemiology and including variation in the duration of immunity in transmission models together are expected to have a profound impact on the transmission dynamics of influenza. Specifically, the total size of epidemics becomes highly irregular, with increased mean peak prevalence and reduced impact of vaccination on infection attack rates.

These phenomena are not only observed in overly simple influenza models, but also in more realistic ones that take important aspects of influenza epidemiology into account. For instance, our findings are preserved in models that distinguish between infection- and vaccination-induced immunity, include both all-or-nothing and leaky immunity, and take age-structure and the essence of human contact patterns into account. In fact, the results stem from the non-linear relation between population susceptibility and infection attack rate (Additional file 1: Text S1), a hallmark of the transmission dynamics of infectious diseases. Hence, it is expected that the above findings will persist in still more complex models. Such extensions could for instance include an explicit distinction between influenza A and B subtypes, heterogeneity in contact rates, or explicit strain dynamics resulting in buildup and loss of immunity [45].

Several limitations and simplifications deserve scrutiny. First, in our analyses we assumed that vaccine efficacy (i.e. primary vaccine failure in our model) is a fixed quantity that is uncorrelated with other model parameters. This is unlikely to be true in practice. For instance, it seems plausible that vaccine efficacy is correlated with the rate of waning immunity, because a highly drifted virus will result in a poor vaccine match. The implication is that the tendency for overshooting in models with variable immunity may be increased further. As a consequence, our results should be conservative, in the sense that they provide a lower bound for the differences in attack rates in models with fixed and variable duration of immunity.

Another potential limitation is that we did not include blunting of antibody immune responses, which has been reported when using one vaccine strain for multiple seasons ([46] and references therein). This is done partly because this finding is still somewhat contentious, while the potential quantitative impact remains uncertain. Nevertheless, it may well be that variation of influenza epidemics between years may be due in part to this phenomenon. Hence, immune blunting could contribute to decreasing vaccine effectiveness, not only directly but also indirectly by increasing variation in infection attack rates. In fact, we argue that any mechanism that increases variation in infection attack rates is expected to reduce the impact of vaccination.

A final simplification worth mentioning is that we assumed that vaccination, the influenza epidemic, and demographic turnover and virus evolution take place on different timescales. Although this seems a reasonable assumption, timescales are not entirely separated. In particular, it is known that influenza evolves during the influenza season, and that the effectiveness of vaccination may wane during the influenza season [44, 47]. Unfortunately, it is not possible to assess the potential quantitative impact of this in the current model, which is built around the assumption that timescales can be separated.

Despite the importance to use a multi-year framework with a variable duration of immunity to describe the epidemiology of influenza and estimate the impact of vaccination, epidemiological studies invariably treat influenza epidemics as if they are independent of one another (reviewed in [15]). This has potential implications for evaluations of vaccine effectiveness, if only because the effectiveness is expected to be lower if by chance the epidemic in a certain year is larger than average, thereby increasing differences in the probability of infection in the unvaccinated and vaccinated subpopulations. Again, this is caused by the non-linearity in the infection dynamics. For instance, vaccinated persons who escape infection in a mild influenza season affecting a small fraction of the population may well be infected in a severe influenza epidemic affecting a larger fraction of the population, even though their level of immunity and match of the vaccine with the virus is identical in both cases. Thus, the epidemiological history in the population needs to be taken into account when estimating the vaccine effectiveness, for example by including the size of past epidemics as covariate in statistical analyses.

With regard to public health implications it is noteworthy that some countries have started vaccination of healthy children to reduce influenza circulation [2, 3], relying on predictions by models with a fixed duration of immunity [18, 20, 23, 33, 35]. Especially in children, however, the long-term effects of vaccination are unknown. Moreover, by increasing the vaccination coverage, it is conceivable that a situation will be created in which the pool of susceptible persons increases gradually over the years until a large epidemic occurs that affects an unexpectedly large fraction of the population.


Our models show that variation in the duration of immunity negatively affects the effectiveness of vaccination of epidemic pathogens. Such variation is well-documented for influenza A, and is caused by the virus’ irregular evolution. Hence, we call for a multi-year perspective in analyses of the transmission dynamics and vaccine effectiveness of influenza.



We thank Rob J. de Boer, Jaap T. van Dissel, Christiaan H. van Dorp, and Marianne A.B. van der Sande for helpful discussions.


Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO grant 645.000.002).

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files. Computer code is publicly available [37].

Authors’ contributions

RGW, JW, and MVB designed and conceptualized the study. RGW wrote computer code and performed model analyses. RGW, JW, and MVB contributed to interpretation of the results. RGW drafted the initial manuscript. RGW, JW, and MVB contributed to critical revisions of the manuscript. All authors read and approved the final manuscript.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

Theoretical Biology, Utrecht University
National Institute for Public Health and the Environment


  1. World Health Organization. Seasonal epidemics and disease burden. 2014. Accessed 17 Oct 2016.Google Scholar
  2. Centers for Disease Control and Prevention. CDC influenza vaccination recommendations. Accessed 17 Oct 2016.Google Scholar
  3. Joint Committee on Vaccination and Immunisation Department of Health England. JCVI statement on the annual influenza vaccination programme - extension of the programme to children. 2013. Accessed 17 Oct 2016.Google Scholar
  4. Health Council of the Netherlands, Fighting the flu. The Hague: Health Council of the Netherlands. 2014. publication no. 2014/16 isbn 978-94-6281-005-1. Scholar
  5. Koelle K, Cobey S, Grenfell B, Pascual M. Epochal evolution shapes the phylodynamics of interpandemic influenza a (H3N2) in humans. Science. 2006; 314(5807):1898–903. doi:10.1126/science.1132745.View ArticlePubMedGoogle Scholar
  6. Smith DJ, Lapedes AS, Jong JCd, Bestebroer TM, Rimmelzwaan GF, Osterhaus ADME, Fouchier RAM. Mapping the antigenic and genetic evolution of influenza virus. Science. 2004; 305(5682):371–6. doi:10.1126/science.1097211.View ArticlePubMedGoogle Scholar
  7. Bedford T, Suchard MA, Lemey P, Dudas G, Gregory V, Hay AJ, McCauley JW, Russell CA, Smith DJ, Rambaut A. Integrating influenza antigenic dynamics with molecular evolution. eLife. 2014; 3:01914. doi:10.7554/eLife.01914.View ArticleGoogle Scholar
  8. Fonville JM, Wilks SH, James SL, Fox A, Ventresca M, Aban M, Xue L, Jones TC, Le NMH, Pham QT, Tran ND, Wong Y, Mosterin A, Katzelnick LC, Labonte D, Le TT, van der Net G, Skepner E, Russell CA, Kaplan TD, Rimmelzwaan GF, Masurel N, de Jong JC, Palache A, Beyer WEP, Le QM, Nguyen TH, Wertheim HFL, Hurt AC, Osterhaus ADME, Barr IG, Fouchier RAM, Horby PW, Smith DJ. Antibody landscapes after influenza virus infection or vaccination. Science. 2014; 346(6212):996–1000. doi:10.1126/science.1256427.View ArticlePubMedPubMed CentralGoogle Scholar
  9. Luksza M, Lässig M. A predictive fitness model for influenza. Nature. 2014; 507(7490):57–61. doi:10.1038/nature13087.View ArticlePubMedGoogle Scholar
  10. Neher RA, Bedford T, Daniels RS, Russell CA, Shraiman BI. Prediction, dynamics, and visualization of antigenic phenotypes of seasonal influenza viruses. Proc Natl Acad Sci. 2016; 113(12):1701–9. doi:10.1073/pnas.1525578113.View ArticleGoogle Scholar
  11. Jefferson T, Rivetti A, Di Pietrantonj C, Demicheli V, Ferroni E. Vaccines for preventing influenza in healthy children. Cochrane Database Syst Rev. 2012; 8:004879. doi:10.1002/14651858.CD004879.pub4.Google Scholar
  12. Jefferson T, Di Pietrantonj C, Al-Ansary LA, Ferroni E, Thorning S, Thomas RE. Vaccines for preventing influenza in the elderly. Cochrane Database Syst Rev. 2010; 2:004876. doi:10.1002/14651858.CD004876.pub3.Google Scholar
  13. Jefferson T, Rivetti D, Rivetti A, Rudin M, Di Pietrantonj C, Demicheli V. Efficacy and effectiveness of influenza vaccines in elderly people: a systematic review. Lancet. 2005; 366(9492):1165–74. doi:10.1016/S0140-6736(05)67339-4.View ArticlePubMedGoogle Scholar
  14. Jefferson T, Di Pietrantonj C, Rivetti A, Bawazeer GA, Al-Ansary LA, Ferroni E. Vaccines for preventing influenza in healthy adults. Cochrane Database Syst Rev. 2010:001269. doi:10.1002/14651858.CD001269.pub4.
  15. Osterholm MT, Kelley NS, Sommer A, Belongia EA. Efficacy and effectiveness of influenza vaccines: a systematic review and meta-analysis. Lancet Infect Dis. 2012; 12(1):36–44. doi:10.1016/S1473-3099(11)70295-X.View ArticlePubMedGoogle Scholar
  16. Manzoli L, Ioannidis JPA, Flacco ME, De Vito C, Villari P. Effectiveness and harms of seasonal and pandemic influenza vaccines in children, adults and elderly: a critical review and re-analysis of 15 meta-analyses. Human Vaccines Immunotherapeutics. 2012; 8(7):851–62. doi:10.4161/hv.19917.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Carrat F, Lavenu A, Cauchemez S, Deleger S. Repeated influenza vaccination of healthy children and adults: borrow now, pay later?Epidemiol Infect. 2006; 134(1):63–70. doi:10.1017/S0950268805005479.View ArticlePubMedGoogle Scholar
  18. Vynnycky E, Pitman R, Siddiqui R, Gay N, Edmunds WJ. Estimating the impact of childhood influenza vaccination programmes in England and Wales. Vaccine. 2008; 26(41):5321–30. doi:10.1016/j.vaccine.2008.06.101.View ArticlePubMedGoogle Scholar
  19. Smith KJ, Lee BY, Nowalk MP, Raymund M, Zimmerman RK. Cost-effectiveness of dual influenza and pneumococcal vaccination in 50-year-olds. Vaccine. 2010; 28(48):7620–5. doi:10.1016/j.vaccine.2010.09.053.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Fisman DN, Tuite AR. Estimation of the health impact and cost-effectiveness of influenza vaccination with enhanced effectiveness in Canada. PLOS ONE. 2011; 6(11):27420. doi:10.1371/journal.pone.0027420.View ArticleGoogle Scholar
  21. Michaelidis CI, Zimmerman RK, Nowalk MP, Smith KJ. Estimating the cost-effectiveness of a national program to eliminate disparities in influenza vaccination rates among elderly minority groups. Vaccine. 2011; 29(19):3525–30. doi:10.1016/j.vaccine.2011.02.098.View ArticlePubMedGoogle Scholar
  22. Clements KM, Chancellor J, Nichol K, DeLong K, Thompson D. Cost-effectiveness of a recommendation of universal mass vaccination for seasonal influenza in the United States. Value Health. 2011; 14(6):800–11. doi:10.1016/j.jval.2011.03.005.View ArticlePubMedGoogle Scholar
  23. Pitman RJ, White LJ, Sculpher M. Estimating the clinical impact of introducing paediatric influenza vaccination in England and Wales. Vaccine. 2012; 30(6):1208–24. doi:10.1016/j.vaccine.2011.11.106.View ArticlePubMedGoogle Scholar
  24. Giglio N, Gentile A, Lees L, Micone P, Armoni J, Reygrobellet C, Crepey P. Public health and economic benefits of new pediatric influenza vaccination programs in Argentina. Human Vaccines Immunotherapeutics. 2012; 8(3):312–22. doi:10.4161/hv.18569.View ArticlePubMedPubMed CentralGoogle Scholar
  25. Ndeffo Mbah ML, Medlock J, Meyers LA, Galvani AP, Townsend JP. Optimal targeting of seasonal influenza vaccination toward younger ages is robust to parameter uncertainty. Vaccine. 2013; 31(30):3079–89. doi:10.1016/j.vaccine.2013.04.052.View ArticlePubMedPubMed CentralGoogle Scholar
  26. Newall AT, Dehollain JP, Creighton P, Beutels P, Wood JG. Understanding the cost-effectiveness of influenza vaccination in children: methodological choices and seasonal variability. PharmacoEconomics. 2013; 31(8):693–702. doi:10.1007/s40273-013-0060-7.View ArticlePubMedGoogle Scholar
  27. Newall AT, Dehollain JP. The cost-effectiveness of influenza vaccination in elderly Australians: an exploratory analysis of the vaccine efficacy required. Vaccine. 2014; 32(12):1323–5. doi:10.1016/j.vaccine.2014.01.017.View ArticlePubMedGoogle Scholar
  28. Zhou L, Situ S, Feng Z, Atkins CY, Fung IC-H, Xu Z, Huang T, Hu S, Wang X, Meltzer MI. Cost-effectiveness of alternative strategies for annual influenza vaccination among children aged 6 months to 14 years in four provinces in China. PLOS ONE. 2014; 9(1):87590. doi:10.1371/journal.pone.0087590.View ArticleGoogle Scholar
  29. Yamin D, Balicer RD, Galvani AP. Cost-effectiveness of influenza vaccination in prior pneumonia patients in Israel. Vaccine. 2014; 32(33):4198–205. doi:10.1016/j.vaccine.2014.05.015.View ArticlePubMedPubMed CentralGoogle Scholar
  30. Clements KM, Meier G, McGarry LJ, Pruttivarasin N, Misurski DA. Cost-effectiveness analysis of universal influenza vaccination with quadrivalent inactivated vaccine in the United States. Human Vaccines Immunotherapeutics. 2014; 10(5):1171–80. doi:10.4161/hv.28221.View ArticlePubMedPubMed CentralGoogle Scholar
  31. Preaud E, Durand L, Macabeo B, Farkas N, Sloesen B, Palache A, Shupo F, Samson SI, Vaccines Europe influenza working group. Annual public health and economic benefits of seasonal influenza vaccination: a European estimate. BMC Public Health. 2014; 14(1):813. doi:10.1186/1471-2458-14-813.View ArticlePubMedPubMed CentralGoogle Scholar
  32. Blommaert A, Bilcke J, Vandendijck Y, Hanquet G, Hens N, Beutels P. Cost-effectiveness of seasonal influenza vaccination in pregnant women, health care workers and persons with underlying illnesses in Belgium. Vaccine. 2014; 32(46):6075–83. doi:10.1016/j.vaccine.2014.08.085.View ArticlePubMedGoogle Scholar
  33. Rose MA, Damm O, Greiner W, Knuf M, Wutzler P, Liese JG, Krüger H, Wahn U, Schaberg T, Schwehm M, Kochmann TF, Eichner M. The epidemiological impact of childhood influenza vaccination using live-attenuated influenza vaccine (LAIV) in Germany: predictions of a simulation study. BMC Infect Dis. 2014; 14(1):40. doi:10.1186/1471-2334-14-40.View ArticlePubMedPubMed CentralGoogle Scholar
  34. Baguelin M, Jit M, Miller E, Edmunds WJ. Health and economic impact of the seasonal influenza vaccination programme in England. Vaccine. 2012; 30(23):3459–62. doi:10.1016/j.vaccine.2012.03.019.View ArticlePubMedGoogle Scholar
  35. Baguelin M, Flasche S, Camacho A, Demiris N, Miller E, Edmunds WJ. Assessing optimal target populations for influenza vaccination programmes: an evidence synthesis and modelling study. PLoS Med. 2013; 10(10):1001527. doi:10.1371/journal.pmed.1001527.View ArticleGoogle Scholar
  36. Xiao Y, Moghadas SM. Impact of viral drift on vaccination dynamics and patterns of seasonal influenza. BMC Infect Dis. 2013; 13:589. doi:10.1186/1471-2334-13-589.View ArticlePubMedPubMed CentralGoogle Scholar
  37. Computer code is publicly available at 2017. doi:10.5281/zenodo.845465. Accessed 18 Aug 2017.
  38. te Beest DE, van Boven M, Hooiveld M, van den Dool C, Wallinga J. Driving factors of influenza transmission in the Netherlands. Am J Epidemiol. 2013; 178(9):1469–77. doi:10.1093/aje/kwt132.View ArticlePubMedGoogle Scholar
  39. Gran JM, Iversen B, Hungnes O, Aalen OO. Estimating influenza-related excess mortality and reproduction numbers for seasonal influenza in norway, 1975-2004. Epidemiol Infect. 2010; 138(11):1559–68. doi:10.1017/S0950268810000671.View ArticlePubMedGoogle Scholar
  40. Wallinga J, Teunis P, Kretzschmar M. Using data on social contacts to estimate age-specific transmission parameters for respiratory-spread infectious agents. Am J Epidemiol. 2006; 164(10):936–44. doi:10.1093/aje/kwj317.View ArticlePubMedGoogle Scholar
  41. Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, Massari M, Salmaso S, Tomba GS, Wallinga J, Heijne J, Sadkowska-Todys M, Rosinska M, Edmunds WJ. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLOS Med. 2008; 5(3):74. doi:10.1371/journal.pmed.0050074.View ArticleGoogle Scholar
  42. Lugnèr AK, Boven Mv, Vries Rd, Postma MJ, Wallinga J. Cost effectiveness of vaccination against pandemic influenza in European countries: mathematical modelling analysis. BMJ. 2012; 345:4445. doi:10.1136/bmj.e4445.View ArticleGoogle Scholar
  43. McLean AR, Anderson RM. Measles in developing countries. Part II, The predicted impact of mass vaccination. Epidemiol Infect. 1988; 100(3):419–42.View ArticlePubMedPubMed CentralGoogle Scholar
  44. Belongia EA, Sundaram ME, McClure DL, Meece JK, Ferdinands J, VanWormer JJ. Waning vaccine protection against influenza A (H3N2) illness in children and older adults during a single season. Vaccine. 2015; 33(1):246–51. doi:10.1016/j.vaccine.2014.06.052.View ArticlePubMedGoogle Scholar
  45. Wikramaratna PS, Kucharski A, Gupta S, Andreasen V, McLean AR, Gog JR. Five challenges in modelling interacting strain dynamics. Epidemics. 2015; 10:31–4. doi:10.1016/j.epidem.2014.07.005. Challenges in Modelling Infectious DIsease Dynamics Challenges in Modelling Infectious DIsease Dynamics.View ArticlePubMedGoogle Scholar
  46. Belongia EA, Skowronski DM, McLean HQ, Chambers C, Sundaram ME, De Serres G. Repeated annual influenza vaccination and vaccine effectiveness: review of evidence. Expert Rev Vaccines. 2017; 16(7):1–14.View ArticlePubMedGoogle Scholar
  47. Boni MF, Gog JR, Andreasen V, Feldman MW. Epidemic dynamics and antigenic evolution in a single season of influenza a. Proc Biol Sci. 2006; 273(1592):1307–16. doi:10.1098/rspb.2006.3466.View ArticlePubMedPubMed CentralGoogle Scholar


© The Author(s) 2017