New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
 Agricultural Sciences
 Anthropology
 Applied Biological Sciences
 Biochemistry
 Biophysics and Computational Biology
 Cell Biology
 Developmental Biology
 Ecology
 Environmental Sciences
 Evolution
 Genetics
 Immunology and Inflammation
 Medical Sciences
 Microbiology
 Neuroscience
 Pharmacology
 Physiology
 Plant Biology
 Population Biology
 Psychological and Cognitive Sciences
 Sustainability Science
 Systems Biology
Dynamical resonance can account for seasonality of influenza epidemics

Contributed by Simon A. Levin, October 6, 2004
Abstract
Influenza incidence exhibits strong seasonal fluctuations in temperate regions throughout the world, concentrating the mortality and morbidity burden of the disease into a few months each year. The cause of influenza's seasonality has remained elusive. Here we show that the large oscillations in incidence may be caused by undetectably small seasonal changes in the influenza transmission rate that are amplified by dynamical resonance.
The underlying cause of seasonal oscillations in influenza incidence (1) remains unclear despite at least 80 years of investigation (2). These oscillations are presumably due to some mechanism that causes seasonal changes in the effective transmission rate of the virus from person to person (often called seasonal forcing). But laboratory experiments and epidemiological studies have failed to establish whether these transmission changes are due to direct effects of temperature and humidity on transmission (3), to changes in mixing patterns [e.g., school terms (4) or simply more time spent indoors], or to other factors, such as increased viral production under winter conditions (evidenced by higher febrile reaction to identical influenza exposures; ref. 5). In fact, as we show here, it may be impossible to establish the underlying cause of seasonality in influenza epidemics, because the large observed oscillations in incidence can be generated by seasonal changes in the transmission rate that are too small to measure.
For an infectious disease that induces permanent immunity, simple deterministic models indicate that disease incidence has an intrinsic tendency to oscillate, even in the absence of exogenous driving factors. But these intrinsic oscillations are always damped (6). Such models are often called susceptible, infectious, recovered (SIR) models, referring to the transition of an individual between susceptible, infectious, and recovered states.
For influenza, immunity is not permanent. After recovering from one antigenic variant of the virus, a person is usually at least partially susceptible to new variants within a few years, due to gradual evolution of the virus (7). We model this process of antigenic drift by simply allowing people to lose their resistance to the circulating virus and hence to move back to the susceptible class after a few years [i.e., a SIRsusceptible (SIRS) model]. For simplicity, we do not explicitly model the exposed population but instead include people infected but not yet infectious in the “I” box. Including an exposed class yields similar results.
In a SIRS model, the intrinsic period of oscillation is approximately , where D is the mean infectious period, and L is the average duration of immunity (4, 6). The basic reproductive number is defined as the expected number of secondary infections that a single infectious individual will cause in a wholly susceptible population.
If the period of endogenous SIRS oscillations (T) is near the period of seasonal forcing (1 year), then these two factors may resonate to produce greatly amplified oscillations in incidence (8). Appropriate parameter values for influenza (9–11), 6–10 days for D, 4–8 years for L, and 4–16 for , yield an endogenous oscillatory period, T, between 0.4 and 1.5 years. Therefore, the intrinsic period of oscillation for influenza is potentially close to the (annual) period of seasonal forcing, in contrast to human viral diseases that engender permanent effective immunity and typically have intrinsic interepidemic periods of 2–5 years (4).
Whereas deterministic models treat the number of people in each state as a continuous variable, in reality, individuals are discrete, and the process is more appropriately modeled as a stochastic Markov chain. The resulting demographic stochasticity causes fluctuations in disease incidence to persist (12). Demographic stochastisticity is more important if the average duration of immunity is long compared with the duration of infection (13). In the case of influenza, the ratio between these two time scales is on the order of a few years to 1–2 weeks, or ≈100.
Methods
The red curves in Fig. 1 represent solutions to the deterministic SIRS epidemic model defined by the equations where t is time in years, N is the total population size, and NSI gives the number of resistant (immune) individuals. The contact rate β(t) varies sinusoidally according to the formula In this formulation, equals Dβ_{0}. The parameter values corresponding to Fig. 1a are n = 500,000, L = 4 yr, D = 0.02 yr, β_{0} = 500 per year, and β_{1} = 0.02. The parameter values in Fig. 1b are n = 500,000, L = 8 yr, D = 0.025 yr, β_{0} = 400 per year, and β_{1} = 0.02. The blue curves in Fig. 1 represent solutions to the discretestate continuoustime Markov chain in which events occur at the exponential rates given in Table 1.
The Markov chain was simulated by using the standard Gillespie algorithm (14).
Results
Fig. 1 shows the effect of resonance, using simple deterministic and stochastic models of influenza dynamics in a population of 500,000 people, with very small sinusoidal forcing of the transmission rate. If resonance is weak (i.e., the endogenous period T is far from 1 year), then the resulting oscillations in incidence are much smaller than are actually observed for influenza. If resonance is strong, however, the resulting dynamics exhibit regular annual epidemics in which incidence increases by a factor of three or more from winter to summer. Fig. 2 shows the endogenous period and the resulting size of oscillations for parameters chosen at random from the plausible range for influenza. Due to resonance, large oscillations (average peaktotrough ratio >5) occur for about onefifth of parameter sets, even when the transmission rate varies by only a few percent.
Researchers have worked to identify and measure a seasonal component of influenza transmission with the goal of explaining large annual fluctuations in incidence. But, as we have seen here using simple models, these large fluctuations may be caused by exogenous seasonal changes in transmission that are too small to detect, amplified by the endogenous population dynamics of the host–pathogen system.
Acknowledgments
We are grateful for valuable discussions with Viggo Andreasen, Tom Reichert, and Lone Simonsen and helpful comments from Herbert Hethcote and an anonymous reviewer. This work was supported by the National Institutes of Health (Grant 1R01GM6072901, to S.A.L.) and the Canadian Institutes of Health Research (to D.J.D.E.).
Footnotes

↵ ‡ To whom correspondence should be addressed. Email: dushoff{at}eno.princeton.edu.

Author contributions: J.D., D.J.D.E., and J.B.P. designed research; J.D. performed research; and J.D., J.B.P., S.A.L., and D.J.D.E. wrote the paper.

Abbreviations: SIR, susceptible, infectious, recovered; SIRS, SIRsusceptible.
 Copyright © 2004, The National Academy of Sciences
References
 ↵
 ↵
 ↵

↵
Anderson, R. M. & May, R. M. (1991) Infectious Diseases of Humans: Dynamics and Control (Oxford Science, Oxford).

↵
Nelson, R. J., Demas, G. E., Klein, S. L. & Kriegsfeld, L. J. (2002) Seasonal Patterns of Stress, Immune Function and Disease (Cambridge Univ. Press, Cambridge).

↵
Wilson, E. B. & Worcester, J. (1945) Proc. Natl. Acad. Sci. USA 31 , 294298. pmid:16578169

↵
Cox, N. J. & Bender, C. A. (1995) Semin. Virol. 6 , 359370.

↵
Dietz, K. (1976) in Lecture Notes in Biomathematics: Mathematical Models in Medicine, eds. Berger, J., Bühler, W., Repges, R. & Tautu, P. (Springer, Berlin), pp. 115.

↵
Douglas, R. G., Jr. (1975) in The Influenza Viruses and Influenza, ed. Kilbourne, E. D. (Academic, New York), pp. 395447.

Fine, P. (1982) in Influenza Models: Prospects for Development and Use, ed. Selby, P. (MTP Press, Lancaster), pp. 1585.
 ↵

↵
Bartlett, M. S. (1957) J. R. Stat. Soc. A120, 4870.
 ↵
 ↵