Public health interventions for epidemics: implications for multiple infection waves
© Wessel et al; licensee BioMed Central Ltd. 2011
Published: 25 February 2011
Epidemics with multiple infection waves have been documented for some human diseases, most notably during past influenza pandemics. While pathogen evolution, co-infection, and behavioural changes have been proposed as possible mechanisms for the occurrence of subsequent outbreaks, the effect of public health interventions remains undetermined.
We develop mean-field and stochastic epidemiological models for disease transmission, and perform simulations to show how control measures, such as drug treatment and isolation of ill individuals, can influence the epidemic profile and generate sequences of infection waves with different characteristics.
We demonstrate the impact of parameters representing the effectiveness and adverse consequences of intervention measures, such as treatment and emergence of drug resistance, on the spread of a pathogen in the population. If pathogen resistant strains evolve under drug pressure, multiple outbreaks are possible with variability in their characteristics, magnitude, and timing. In this context, the level of drug use and isolation capacity play an important role in the occurrence of subsequent outbreaks. Our simulations for influenza infection as a case study indicate that the intensive use of these interventions during the early stages of the epidemic could delay the spread of disease, but it may also result in later infection waves with possibly larger magnitudes.
The findings highlight the importance of intervention parameters in the process of public health decision-making, and in evaluating control measures when facing substantial uncertainty regarding the epidemiological characteristics of an emerging infectious pathogen. Critical factors that influence population health including evolutionary responses of the pathogen under the pressure of different intervention measures during an epidemic should be considered for the design of effective strategies that address short-term targets compatible with long-term disease outcomes.
Epidemics of infectious diseases have been observed throughout history, with substantial variability in their dynamical patterns. The 1918 influenza pandemic is a notorious case documented as the most devastating epidemic with over 50 million deaths and multiple outbreaks in many geographic areas worldwide [1, 2]. Distinct pandemic infection waves were recorded with an 8 to 15 week interval; the latter were more severe than the first and were associated with the majority of deaths [2, 3]. Although several factors may be involved, such as the effect of seasonal changes, demographics, and evolution of the virus, the true mechanism by which subsequent waves occur is not fully understood. Nor is it clearly understood how different control measures and strategies for deployment of limited health resources may interfere with disease dynamics and the occurrence of later infection waves.
Recent epidemiological and modelling studies have attempted to provide explanatory theories for the mechanisms of multiple outbreaks of an infectious pathogen capable of establishing an epidemic [2, 4–10]. Spontaneous behavioural changes (e.g., a change in the number of contacts due to modified behaviour of susceptible individuals) have been shown to affect the course of infection events and produce subsequent outbreaks in an epidemic episode . This has been further investigated through modelling “concerned awareness” of individuals that may result in contagion dynamics of fear and disease , and the implementation of public health control measures (e.g., social distancing) that may interfere with the individuals’ contact patterns during the epidemic . Co-infection has also been suggested as a possible explanation for multiple infection outbreaks as a result of increased transmissibility in co-infected individuals and non-synchronicity in the time course of the two co-circulating infections . Other possible mechanisms include transient post-infection immunity and evolutionary changes that may occur in the characteristics of the infectious pathogens [2, 4, 10].
In this study, we consider the occurrence of multiple infection waves of a pathogen from a public health perspective, and develop mathematical models to investigate how intervention measures may affect the transmission dynamics in a population. Specifically, we are interested in exploring the impact of changes in policy-relevant parameters on the patterns of disease spread during the course of an epidemic. These parameters may reflect the effectiveness of intervention strategies (e.g., treatment or isolation of infected cases) in reducing disease transmission, or their epidemiological consequences (e.g., emergence of drug resistance), and may therefore play an important role in determining the outcome of disease control activities. The significance of this work thus relates to the process of public health decision-making, in particular when confronting the emergence of a novel infectious disease with substantial uncertainty regarding the epidemiological characteristics of the invading pathogen.
For the purpose of this investigation, we develop both mean-field and stochastic epidemiological models that describe the transmission dynamics of a disease in the population, and incorporate treatment and isolation of infected cases as control measures. We parameterize these models to simulate the spread of influenza as a case study, and determine the impact of control parameters on disease dynamics. We illustrate the occurrence of multiple infection waves associated with different treatment levels and the development of drug resistance in the population under the scenario of limited capacity for treatment and isolation of infectious individuals. We compare the results obtained by simulating the mean-field model with those observed in the stochastic model, and discuss our findings in the context of epidemiology and public health.
The model structure
To formulate the models for describing disease epidemic, we assume that the population is initially entirely susceptible to the infectious pathogen. It is assumed that the infection can be treated with drugs, but the pathogen may develop resistance during the course of treatment with potential for transmission. Since resistance emergence may impose fitness cost on pathogen replication and transmission , we assume that the drug-resistant pathogen is less transmissible than the drug-sensitive pathogen. Treatment is assumed to reduce transmissibility of the drug-sensitive infection, but remains ineffective against drug-resistant infection. We also assume that the recovery from infection confers immunity to re-infection with either drug-sensitive or resistant pathogens. Considering epidemics with relatively short time-courses, we ignore the effect of recruitment, natural death, and other demographic variables of the population.
Baseline transmission rate of the sensitive infection
Expose period (duration between exposure and start of infectiousness)
Duration of infectiousness
Rate of developing drug resistance following start of treatment
Relative transmissibility of the sensitive infection following treatment
Relative transmissibility of the resistant infection
Fraction of infectious individuals diagnosed for treatment or isolation
Fraction of treated infectious individuals without isolation
Capacity for treatment and isolation of infectious individuals
Basic reproduction number of the sensitive infection
Details of the model in its stochastic form are provided in the Appendix.
which highlights the importance of isolation for controlling the spread of resistant infection.
Simulations and results
To simulate the models, we considered influenza infection as a case study, for which emergence and spread of drug-resistance during an outbreak can result from treatment of infected individuals. We assumed that the epidemic is triggered by a drug-sensitive influenza virus, and investigated the role of several key model parameters in changing the epidemic patterns and generating multiple waves of infection. These parameters include the fractions of infected individuals identified for treatment or isolation, and the basic reproduction number of disease which varies within the estimated range published in the literature (Table 1). Since public health resources may be limited during an epidemic, we also defined a parameter (T c ) as the capacity for treatment of infected individuals including those who are isolated (i.e., the percentage of the total population that can be treated). To illustrate various scenarios, we initially seeded a susceptible population of size S0 = 10,000 with E0 = 10 individuals exposed to the sensitive virus, and assumed that treatment can result in the emergence of resistance with the relative transmissibility δ r = 0.8 during the outbreak. Other parameter values are given in Table 1.
As the reproduction number of the sensitive infection increases (Figures 2c-2d), higher treatment levels are required for the resistant infection to prevail and cause a significant outbreak . For a reduced level of treatment and a higher transmissibility of the sensitive virus, corresponding to the epidemic profile in Figure 2c, we observed two infection waves, both of which are caused by the spread of the sensitive virus, with generation of very few cases of resistant infection. In this scenario, run-out occurs before epidemic is contained, and a second infection wave takes place. Similar dynamics can occur with two subsequent waves of resistant infections for a significantly higher treatment level (Figure 2d). However, the second wave that occurs after the treatment capacity is fully dispensed (run-out scenario) leads to a major reduction in susceptibility of the population; thereby ending the epidemic. These simulations indicate that multiple infection waves could occur due to limited resources for treatment/isolation of infected cases, the ways that such resources are deployed during the outbreak, the evolutionary responses of the pathogen to control measures (e.g., emergence of drug resistance), or a combination thereof. We performed further experiments with small changes in these parameters, and observed significant influences on the epidemic dynamics that can be associated with the elimination or creation of an infection wave. It is worth noting that the above scenarios can take place even for sufficient drug stockpiles for which run-out does not occur, if a policy for adaptation (e.g., reduction) of treatment at the population level is implemented due to wide-spread drug-resistance .
For comparison purposes, we simulated the stochastic version of the model using a Markov Chain Monte Carlo method and observed sequences of infection waves for different sets of parameter values (see Appendix). Consistent with previous observations , the stochastic model displays a later peak time of infection waves (with lower magnitudes) than the homogeneous mean-field model. This depends not only on the treatment level, but also on other parameters involved in the spread of sensitive and resistant infections, such as the reduction in the potentially infectious contacts and the fitness of resistance. Furthermore, stochastic effects can play a significant role in determining disease dynamics even during the outbreak well past the initial establishment phase of the epidemic. This is illustrated in Figure 4c of the Appendix that the epidemic dies out after the first outbreak in the stochastic model; whereas a second wave of infection takes place in the mean-field model with a larger magnitude compared to the first outbreak.
Stellar advances in the prevention and management of infectious diseases have been achieved since the great influenza pandemic of 1918. Yet, emerging pathogens often inflict incalculable devastation to humanity. The global mobilization with rapid international transportation between populations makes the impact of such diseases even more dramatic with potential socioeconomic upheaval. This was recognized in 2003 with the appearance of severe acute respiratory syndrome (SARS) as the first major infectious disease threat of the 21st century , and was recently experienced with the worldwide spread of a swine-origin influenza A virus H1N1, that led the World Health Organization to declare this virus as the cause of an influenza pandemic on June 11, 2009 . Public health responses to the emergence of new diseases often involve difficult decisions on optimal use of health resources over very short timelines. Such decisions are further confounded by substantial uncertainties regarding the epidemiological characteristics of the novel infectious pathogen, the effectiveness of public health intervention strategies, and the evolutionary responses of the pathogen under the pressure of control measures . From a population health perspective, it is therefore imperative to look beyond short-term targets and account for long-term disease outcomes in strategy development and implementation. This is particularly important for preventing multiple infection outbreaks that may result from imprudent use of resources or unintended adverse consequences of disease containment strategies.
Given the historical evidence for the occurrence of multiple infection waves [2, 3, 7], several modelling studies have attempted to provide explanatory theories for these events in a single epidemic course [2, 5, 6, 8–10]. In this study, we developed mean-field and stochastic models to investigate possible causes of sequential outbreaks from a public health perspective. Our results show that epidemic dynamics can be substantially affected by factors that influence policy design and implementation (e.g., treatment level or isolation of infected individuals), and parameters that determine the effectiveness and consequences of control measures (e.g., reduction in infectiousness due to treatment or emergence of drug-resistance). Furthermore, the initial number of infections can influence disease outcomes. While mean-field and stochastic models may exhibit similar epidemic behaviour, we also observed differences in their predictions in terms of the speed with which disease spreads through the population (with further delay in the peak time of outbreaks in the stochastic model); the magnitudes of infection outbreaks; and more importantly, the occurrence of infection waves (see Appendix). The latter is particularly influenced by stochastic effects, in addition to the structure of contact patterns and heterogeneity in population interactions . Previous work [4, 24] provides a solid foundation for extension of this study through the development of network dynamical models of disease transmission in which heterogeneous contacts between individuals are accounted for.
In this study, we simplified the models and included compartments corresponding to some possible stages of a disease; yet we understand that different pathogens may cause infections with different clinical manifestations and infectiousness periods. For example, influenza is known to have a short latent period of less than 2 days before becoming infectious , followed by a pre-symptomatic infection during which disease can be transmitted without showing clinical symptoms; however, the latent period of SARS is estimated to be longer and may be comparable to the duration of a complete course of influenza infection . It is also well-documented that influenza can be transmitted in asymptomatic form without developing clinical symptoms ; while evidence for asymptomatic transmission of SARS is rather scant. These discrepancies in infection stages of human diseases, combined with the ability of the pathogens to overcome the pressures that are applied to limit their replication and spread, can profoundly impact not only the feasibility and effectiveness of control measures, but also the dynamics of disease over the course of an epidemic. Our study highlights these considerations for further investigation, while demonstrating possible mechanisms for the occurrence of multiple infection waves in a single epidemic. Future research in this direction should address some limitations of the present study, including a systematic exploration of parameter space to characterize which intervention parameter regimes are more likely to give rise to sequences of infection outbreaks, and to determine the sensitivity of model outputs (epidemic dynamics) on parameter changes.
Although models considered here are simulated for influenza infection as a case study, understanding the interplay between intervention parameters, evolutionary responses of the pathogens, and epidemic dynamics remains a critical objective of public health for many diseases , including HIV, tuberculosis, malaria, and several bacterial infections. Such diseases often share common features, including the emergence and prevalence of drug resistant pathogens under the pressure of drug treatment. The initial rise of resistance is generally associated with fitness costs that make the resistant pathogen less capable of competing with the sensitive pathogen (as the dominant competitor) in a given host population . However, evolutionary mechanisms (e.g., compensatory mutations ) may improve the fitness of resistant pathogens, and therefore intervention measures may result in further selection of resistance, as has been documented for the global spread of seasonal influenza drug resistance that appears to be associated with fitness enhancement processes . This suggests that future modelling efforts should integrate factors that govern pathogen-host interactions with the mechanisms of disease epidemiology to guide public health in devising novel and effective means of infection control.
Appendix: Stochastic model
Possible transitions between model compartments that can occur in Δt units of time
transition during Δt
infection of a susceptible
S(t) → S(t) – 1
β(I + δ T I T + δ r (I r + I T,r ))S(Δt) + o(Δt)
increase in exposure (sensitive infection)
E(t) → E(t) + 1
β(I + δ T I T )S(Δt) + o(Δt)
decrease in exposure (sensitive infection)
E(t) → E(t) – 1
µ E E(Δt) + o(Δt)
increase in exposure (resistant infection)
E r (t) → E r (t) + 1
δ r β(I r + I T,r )S(Δt) + o(Δt)
decrease in exposure (resistant infection)
E r (t) → E r (t) – 1
µ E E r (Δt) + o(Δt)
increase in untreated sensitive infection
I(t) → I(t) + 1
(I – p)µ E E(Δt) + o(Δt)
recovery from treated infection
I(t) → I(t) – 1
γI(Δt) + o(Δt)
increase in untreated resistant infection
I r (t) → I r (t) + 1
(I – p)µ E E r (Δt) + o(Δt)
recovery from untreated resistant infection
Ir(t) → I r (t) – 1
γI r (Δt) + o(Δt)
increase in treated sensitive infection
I T (t) → I T (t) + 1
pqµ E E(Δt) + o(Δt)
recovery from treated sensitive infection
I T (t) → I T (t) – 1
(α + γ)I T (Δt) + o(Δt)
increase in treated resistant infection
I T,r (t) → I T,r (t) + 1
pqµ E E(Δt) + αIΔt + o(Δt)
recovery from treated resistant infection
I T,r (t) → I T,r (t) – 1
γI T,r (Δt) + o(Δt)
increase in isolated infection
J(t) → J(t) + 1
p(1 – q)µ E (E + E r )(Δt) + o(Δt)
recovery from disease in isolation
J(t) → J(t) – 1
γJ(Δt) + o(Δt)
Algorithm for stochastic simulations
For simulating the stochastic model, we used the Markov Chain Monte Carlo method, with an initial E(0) = 10 exposed individuals to sensitive infection in a population of S0 = 10, 000 susceptibles. A key parameter in these simulations is the step-size of the Monte Carlo method. Using a fixed step-size requires a large number of steps to guarantee that the transitions between subpopulations take place and disease transmission can occur, which is computationally very demanding in terms of both timing and resources. To reduce such a computational load, we implemented an adaptive step-size method  to estimate the transition time to the next event (Δt) by calculating the sum of the frequencies of all possible events, given by η = β(I + δ T I T )S(t) + δ r β(I r + I T,r )S(t) + (1 – p)µ E (E + E r ) + pqµ E (E + E r ) + αI T + p(1 – q)µ E (E + E r ) + γ(I + I r + I T,r + J + I T ). Then, by choosing Δt = U1/η, where U1 is uniform distribution in the interval [0,1], we ordered all possible events as an increasing fraction of η and generated another uniform deviate (U2 ∈ [0,1]) to determine the nature of the next event. For the convergence of the results, we ran these simulations for 1000 samples, and considered the average of sample realizations of the stochastic process to generate infection curves.
This research was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC), and the Mathematics of Information Technology and Complex Systems (MITACS).
This article has been published as part of BMC Public Health Volume 11 Supplement 1, 2011: Mathematical Modelling of Influenza. The full contents of the supplement are available online at http://0-www.biomedcentral.com.brum.beds.ac.uk/1471-2458/11?issue=S1.
- Cox NJ, Tamblyn SE, Tam T: Influenza pandemic planning. Vaccine. 2003, 21: 1801-1803. 10.1016/S0264-410X(03)00076-8.View ArticlePubMedGoogle Scholar
- Taubenberger JK, Morens DM: 1918 influenza: the mother of all pandemics. Emerg Infect Dis. 2006, 12: 15-22.PubMed CentralView ArticlePubMedGoogle Scholar
- Gani R, Hughes H, Fleming D, Griffin T, Medlock J, Leach S: Potential impact of antiviral drug use during influenza pandemic. Emerg Infect Dis. 2005, 9 (11): 1355-1362.View ArticleGoogle Scholar
- Alexander ME, Dietrich SM, Hua Y, Moghadas SM: A comparative evaluation of modelling strategies for the effect of treatment and host interactions on the spread of drug resistance. J Theor Biol. 2009, 259: 253-263. 10.1016/j.jtbi.2009.03.029.View ArticlePubMedGoogle Scholar
- Caley P, Philp DJ, McCracken K: Quantifying social distancing arising from pandemic influenza. J R Soc Interface. 2008, 5: 631-639. 10.1098/rsif.2007.1197.PubMed CentralView ArticlePubMedGoogle Scholar
- Epstein JM, Parker J, Cummings D, Hammond RA: Coupled contagion dynamics of fear and disease: mathematical and computational explorations. PLoS ONE. 2008, 3: e3955-10.1371/journal.pone.0003955. doi:10.1371/journal.pone.0003955PubMed CentralView ArticlePubMedGoogle Scholar
- Hsieh Y.-H, Chen CWS: Turning points, reproduction number, and impact of climatological events for multi-wave dengue outbreaks. Trop Med Int Health. 2009, 14 (6): 628-638. 10.1111/j.1365-3156.2009.02277.x.View ArticlePubMedGoogle Scholar
- Merler S, Poletti P, Ajelli M, Caprile B, Manfredi P: Coinfection can trigger multiple pandemic waves. J Theor Biol. 2008, 254: 499-507. 10.1016/j.jtbi.2008.06.004.View ArticlePubMedGoogle Scholar
- Poletti P, Caprile B, Ajelli M, Pugliese A, Merler S: Spontaneous behavioural changes in response to epidemics. J Theor Biol. 2009, 260: 31-40. 10.1016/j.jtbi.2009.04.029.View ArticlePubMedGoogle Scholar
- Rios-Doria D, Chowell G: Qualitative analysis of the level of cross-protection between epidemic waves of the 1918-1919 influenza pandemic. J Theor Biol. 2009, 261: 584-92. 10.1016/j.jtbi.2009.08.020.View ArticlePubMedGoogle Scholar
- Domingo E, Holland JJ: RNA viruses mutations and fitness for survival. Annu Rev Microbiol. 1997, 51: 151-178. 10.1146/annurev.micro.51.1.151.View ArticlePubMedGoogle Scholar
- Diekmann O, Heesterbeek JAP: Mathematical Epidemiology of Infectious Diseases. 2000, Chichester: WileyGoogle Scholar
- Brauer F, Castillo-Chavez C: Mathematical Models in Population Biology and Epidemiology. Text in Applied Mathematics. 2001, New York: Sprringer-VerlagGoogle Scholar
- van den Driessche P, Watmough J: Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2000, 180: 29-48. 10.1016/S0025-5564(02)00108-6.View ArticleGoogle Scholar
- Arino J, Bowman CS, Moghadas SM: Antiviral resistance during pandemic influenza: implications for stockpiling and drug use. BMC Infect Dis. 2009, DOI:10.1186/1471-2334-9-8Google Scholar
- Moghadas SM: Management of drug resistance in the population: influenza as a case study. Proc R Soc B. 2008, 275: 1163-1169. 10.1098/rspb.2008.0016.PubMed CentralView ArticlePubMedGoogle Scholar
- Lessler J, Reich NG, Brookmeyer R, Perl TM, Nelson KE, Cummings DA: Incubation periods of acute respiratory viral infections: a systematic review. Lancet Infect Dis. 2009, 9: 291-300. 10.1016/S1473-3099(09)70069-6.PubMed CentralView ArticlePubMedGoogle Scholar
- Moghadas SM, Bowman CS, Röst G, Wu J: Population-wide emergence of antiviral resistance during pandemic influenza. PLoS ONE. 2008, 3: e1839-10.1371/journal.pone.0001839.PubMed CentralView ArticlePubMedGoogle Scholar
- Longini IM, Nizam A, Xu S, Ungchusak K, Hanshaoworakul W, Cummings DAT, Halloran ME: Containing pandemic influenza at the source. Science. 2005, 309: 1083-1087. 10.1126/science.1115717.View ArticlePubMedGoogle Scholar
- Ferguson NM, Cummings DAT, Cauchemez S, Fraser C, Riley S, Meeyai A, Iamsirithaworn S, Burke DS: Strategies for containing an emerging influenza pandemic in Southeast Asia. Nature. 2005, 437: 209-214. 10.1038/nature04017.View ArticlePubMedGoogle Scholar
- Gumel AB, Ruan S, Day T, Watmough J, Brauer F, van den Driessche P, Gabrielson D, Bowman C, Alexander ME, Ardal A, Wu J, Sahai BM: Modelling strategies for controlling SARS outbreaks. Proc R Soc B. 2004, 271: 2223-2232. 10.1098/rspb.2004.2800.PubMed CentralView ArticlePubMedGoogle Scholar
- World now at the start of 2009 influenza pandemic [press release]. Geneva (Switzerland): World Health Organization, http://www.who.int/mediacentre/news/statements/2009/h1n1_pandemic_phase6_20090611/en/index.html; accessed April 7, 2010
- Moghadas SM, Day T, Bauch C, Driedger SM, Brauer F, Greer AL, Yan P, Wu J, Pizzi NJ, Fisman DN: Modelling of pandemic influenza: a guide for the perplexed. CMAJ. 2009, 181: 171-173. 10.1503/cmaj.090885.View ArticleGoogle Scholar
- Débarre F, Bonhoeffer S, Regoes RR: The effect of population structure on the emergence of rug resistance during influenza pandemics. J R Soc Interface. 2007, 4: 893-906. 10.1098/rsif.2007.1126.PubMed CentralView ArticlePubMedGoogle Scholar
- Fraser C, Riley S, Anderson RM, Ferguson NM: Factors that make an infectious disease outbreak controllable. Proc Natl Acad Sci USA. 2004, 101: 6146-6151. 10.1073/pnas.0307506101.PubMed CentralView ArticlePubMedGoogle Scholar
- Moghadas SM: Gaining insights into human viral diseases through mathematics. Eur J Epidemiol. 2006, 21: 337-342. 10.1007/s10654-006-9007-z.View ArticlePubMedGoogle Scholar
- Handel A, Regoes RR, Antia R: The role of compensatory mutations in the emergence of drug resistance. PLoS Comput Biol. 2006 Oct;2(10):e137. 2006, 2 (10): 10.1371/journal.pcbi.0020137.Google Scholar
- Simonsen L, Viboud C, Grenfell BT, Dushoff J, Jennings L, Smit M, Macken C, Hata M, Gog J, Miller MA, Holmes EC: The genesis and spread of reassortment human influenza A/H3N2 viruses conferring adamantane resistance. Mol Biol Evol. 2007, 24: 1811-1820. 10.1093/molbev/msm103.View ArticlePubMedGoogle Scholar
- Rohani P, Keeling MJ, Grenfell BT: The interplay between deterministic and stochasticity in childhood diseases. Am Nat. 2002, 159: 469-481. 10.1086/339467.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.