- Research article
- Open Access
The role of cellular immunity in Influenza H1N1 population dynamics
BMC Infectious Diseases volume 12, Article number: 329 (2012)
Pre-existing cellular immunity has been recognized as one of the key factors in determining the outcome of influenza infection by reducing the likelihood of clinical disease and mitigates illness. Whether, and to what extent, the effect of this self-protective mechanism can be captured in the population dynamics of an influenza epidemic has not been addressed.
We applied previous findings regarding T-cell cross-reactivity between the 2009 pandemic H1N1 strain and seasonal H1N1 strains to investigate the possible changes in the magnitude and peak time of the epidemic. Continuous Monte-Carlo Markov Chain (MCMC) model was employed to simulate the role of pre-existing immunity on the dynamical behavior of epidemic peak.
From the MCMC model simulations, we observed that, as the size of subpopulation with partially effective pre-existing immunity increases, the mean magnitude of the epidemic peak decreases, while the mean time to reach the peak increases. However, the corresponding ranges of these variations are relatively small.
Our study concludes that the effective role of pre-existing immunity in alleviating disease outcomes (e.g., hospitalization) of novel influenza virus remains largely undetectable in population dynamics of an epidemic. The model outcome suggests that rapid clinical investigations on T-cell assays remain crucial for determining the protection level conferred by pre-existing cellular responses in the face of an emerging influenza virus.
In spring 2009, a new triple reassortant strain of the influenza A/H1N1 virus emerged. Exhibiting its unique genome composition with rapid global spread , this new strain caused the first pandemic of the 21st century. The genome combination of the 2009 H1N1 strain was antigenically distinct from the circulating seasonal influenza subtypes, seasonal H1N1 (sH1N1) and H3N2 [2, 3]. With the exception of older individuals (> 60 years of age), all other age-groups did not confer any effective cross-antibody protection against this strain [4, 5]. The lack of cross-reactive neutralizing antibodies in a large fraction of the population made the global spread of the virus readily possible. However, compared to the severity of past pandemics of 1918 H1N1, 1957 H2N2 and 1968 H3N2 , the 2009 H1N1 pandemic appeared relatively mild with a death toll of ≈ 300,000 resulting from respiratory and cardiovascular complications globally . The unexpected mild nature and observed prolonged incubation period for this strain [8–10], received an increased attention towards the role of pre-existing cellular immunity. Population-wide studies suggest the existence of an immune response induced from conserved epitopes between the 2009 H1N1 and seasonal influenza strains [11–13]. An immunoinformatics study on the conservancy analysis and influenza epitope prediction  revealed the high level existence of Hemagglutinin (HA) CD4+ T cell epitope conservancy (≈ 95%−100%) between the sH1N1 and 2009 H1N1 isolates. However, no considerable conservancy between subtypes H3N2 and 2009 H1N1 was observed. This investigation hypothesized that the availability of cross-conserved epitopes and Major Histocompatibility Complex class II, HLA-DRB1 alleles essential to the activation of the T cell repertoire, may have played a role in reducing the severity of the 2009 H1N1 infection despite the lack of virus-specific antibody titres. A recent human experimental study demonstrated that conserved and common epitope-specific pre-existing CD4+ T cell immunity, but not CD8+, plays a critical role in limiting viral shedding and severity of infection in the absence of antibody titers . Other studies have also shown recall responses of CD4+T helper memory cells to these shared conserved epitopes [14, 15].
In this paper, we developed a mathematical model to explore the interplay between the individual cellular cross-reactivity and the population spread of disease, linking micro-dynamics with macro-dynamics within an immuno-epidemiology framework. The goal is to determine whether the cross-reactivity reported in our previous work  can be observed at the population level infection dynamics, i.e., to measure the relative change in the magnitude and timing of the epidemic peak.
Our methodology is based on the development of a continuous time Monte-Carlo Markov Chain (MCMC) model to simulate the transmission dynamics of a novel influenza virus by taking into account pre-existing immunity as a result of prior exposure to seasonal influenza strains.
We divided the population into classes of individuals susceptible to the new virus with no prior exposure (S s ) and with prior exposure (S p ) having some level of pre-existing T-cell immunity. Upon exposure to the virus, these individuals move to the exposed classes (E s , E p ), and after the exposed period has elapsed move to the asymptomatically infected classes (A s , A p ) or symptomatically infected classes (I s , I p ). Finally, infected individuals move to the recovered class (R) upon recovery, gaining immunity against re-infection. A schematic of our model is provided in Figure 1. We defined ⋀ s and ⋀ p as the force of infection for the S s and S p populations, given by
where β is the baseline transmission rate, 0 < δ a < 1 denotes a reduction in transmissibility of asymptomatic infection, and 0 < δ p < 1 represents a possible reduction in the susceptibility of individuals in the S p subpopulation. We note that in the absence of specific antibodies, the reduction in susceptibility to infection will likely be very low (δ p ≈ 1), even in the presence of T-cell immunity. However, this cellular immunity can reduce infectiousness (i.e., that corresponds to a reduction in transmissibility in the model), and mitigate illness (i.e., which corresponds to a higher probability of undergoing asymptomatic infection) .
The stochastic model for the transmission dynamics of infection is described in Table 1. The MCMC simulation moves forward in time through the determination of transition probabilities which, in turn, determine the amount of time that elapses between events. We obtained transition probabilities in the MCMC simulation using the rates defined in Table 1. These transition probabilities were then compared to random numbers generated in the simulation and this determines which event will occur. The classes of subpopulations were subsequently adjusted to reflect this event and new transition probabilities were obtained. This procedure continued until infection was wiped out from the population (i.e., when the exposed and infected classes were equal to zero).
Parameterization of the model
We considered a population of 2000 susceptible individuals, and divided it into two groups based on their infection history with seasonal influenza A virus strains, H3N2 and sH1N1. Data from 2001 to 2009 indicate that subtype H3N2 has been the dominant circulating strain, making up over 80% of the viruses typed and subtyped in Europe [16, 17]. Thus, we assumed an upper bound of 20% of individuals considered susceptible to the new virus with prior exposure to the sH1N1 strain in the pre-pandemic era, and considered a lower bound of 80% of the population fully susceptible with no prior exposure to sH1N1.
Cross-reactive pre-existing T-cell memory responses can provide protection by reducing the susceptibility of individuals S p (this may be very low reduction), and decreasing the severity of disease and its transmissibility [13, 18–20]. We assumed that the 2009 H1N1 susceptibility of individuals with prior exposure to sH1N1 is reduced due to T-cell cross-reactivity of conserved epitopes (by a factor δ p ), but remains unchanged for those who had no prior exposure to sH1N1 (see [11, 12]). Cross-reactivity was assigned at a level of 52% for S p , based on conserved T-cell epitopes [11, 12]. To reflect the reduction in the severity of illness and transmission, we first assumed that the probability of an exposed individual in E p developing symptomatic infection was reduced compared to that of an exposed individual in E s , i.e., 0 < p p < p s < 1. We also assumed that the transmission of infection from asymptomatically infectious individuals is reduced (by a factor δ a ) compared to symptomatically infectious cases. In our model, we let the latent period 1/γ p (for those with prior exposure to sH1N1) vary between the length of the latent period 1/γ s and the total period of infection 1/γ s + 1/α s (for those with no prior exposure to sH1N1). Furthermore, it was assumed that the total illness period following exposure is the same for all infected individuals (i.e., symptomatic and/or asymptomatic). Summarizing the above assumptions, we may highlight two main points: individuals with prior exposure to sH1N1 have (i) lower infectiousness and therefore lower transmissibility if infectious; and (ii) lower probability of developing symptomatic infection if exposed; compared to those with no prior exposure.
Note that existing literature does not provide any information on the total period of infection (i.e., 1/γ + 1/α) following exposure. While a prolonged incubation period has been reported as a result of pre-existing immunity , clinical experiments suggest a reduced period of illness . Due to the lack of sufficient data and information, we have steered clear of changing the period of infection for those with prior exposure, and assumed that 1/γ p + 1/α p = 1/γ s + 1/α s . Parameter values for our model described in Table 2 are derived from the above assumptions and the published literature for seasonal and the 2009 H1N1 pandemic. We seeded simulations with initial numbers of infections I s (0) = 4 and I p (0) = 2. These initial values were chosen so that infection would progress in the population (i.e., if these were too small many simulations would clear the infection due to stochasticity).
Basic reproduction number
The basic reproduction number (denoted by R 0) is defined as the the average number of secondary cases generated by a single infectious case introduced into an entirely susceptible population [30, 31]. According to this definition, it is expected that an epidemic will occur if R 0 > 1, and die out if R 0 < 1. For our model, R 0 was determined to be:
where N 0, , correspond to the initial sizes of the total population, the population of susceptibles with no cross-reactivity, and the population of susceptibles with 52% cross reactivity (considered as the level of immunity), respectively. The reader may consult  for a review of different methods of R 0 derivation. Other parameters in the expression of R 0are described in the model diagram illustrated in Figure 1. We used estimates of R 0 to calculate the transmission rate β for a given set of parameter values from Table 2. Reported R 0 values for initial spread of the 2009 H1N1 virus lie in the range 1.4−1.6 in Mexico, 1.7−1.8 in the United States, and 1.25−1.38 in Ontario, Canada [10, 21–24]. Note that, R 0 from these studies is, in fact, a measurement of the effective reproduction number (denoted by R e ) since prior immunity exists in the initial population . To be consistent with [10, 21–24], we refer to these measurements as R 0.
Results and discussion
Since immune cross-reactivity is expected to prolong the incubation period and reduce the severity of illness [13, 18–20], we simulated the stochastic dynamical model to investigate the effect of these factors in the epidemic profile. Simulations were run and mean and standard deviation (for 1000 runs) were recorded, when the initial population with cross-protection varies from 5% to 20% of the total population, the infectious period decreases from 5 to 0.5 days, and the proportion of the exposed class with cross protection E p becoming symptomatic I p varies from 0% to 60%. Figure 2 shows the variation in the peak magnitude (left column) and time to this peak (right column) for the infectious subpopulation I s + A s when p p = 0 (top), 0.05p s (middle) and p s (bottom). In each scenario, the exposed period lasts for on average 1.5 (black), 2.6 (red), 4.3 (green), 6 (blue) days, and the initial size of the population with cross protection varies from 5% to 20%. The figure shows the mean of the simulation runs and two standard errors from the mean, where (standard error = SD/(number of simulations)). These simulations demonstrate that a change in the peak magnitude (Figure 2, left column) can result from a change in the initial size of the population with pre-existing immunity, but the change in infectious period and the fraction of E p moving to the I p class have virtually no effect. The time of the peak magnitude was not affected by any change in these model parameters (Figure 2, right column).
To determine the effect of cross protection δ p on the peak magnitude and time to peak, we also determined these data points for the case where cross reactivity provided no protection to individuals previously exposed to sH1N1 (δ p = 1), and only provided reduced illness and transmission. Again, we observed that, as the initial size of the population with prior exposure to sH1N1 increases, the peak magnitude of A s + I s subpopulation decreases, but the time to peak was not affected. We also observed that variation in the exposed period from 1.5−6 days and the fraction of the E p class moving to the infectious symptomatic class I p from 0 < p p < p s played no significant effect in the model outcomes. The results from these simulations differed very little from the case when cross protection was included in the simulations (when δ p = 0.48; simulations not shown). Further simulations demonstrate that our results hold true for larger population sizes of susceptible individuals.
During an epidemic, infected individuals who can be identified (either as influenza-like illness or laboratory confirmed case) constitute a fraction of the subpopulation I s with symptoms. In our simulations, the peak magnitude of this subpopulation changes very little, and the time to peak is unchanged from that found above. Hence it is unlikely that the impact of pre-existing cellular immunity would be observed in passive surveillance data associated with case reporting.
The 2009 H1N1 pandemic has been characterized as the mildest pandemic on record [8–10]. Typically, disease was mild in all age groups, which is not common to all pandemics . Recent studies of the 2009 H1N1 pandemic suggest that this may be attributed to the host counter responses to infection, induced from the pre-existing CD4+ and CD8+T cell responses to the conserved epitopes [18–20]. Our previous work found that there is a high CD4+ T-cell epitope conservancy and Major Histocompatibility Complex (MHC) class II, HLA-DRB1 promiscuity between sH1N1 (which has been circulating for some 33 years), and the 2009 H1N1 pandemic strain . We hypothesized that the consistent circulation and infective ability of the sH1N1 strain may be a critical factor to the atypical mild nature of 2009 H1N1, because both strains shared common and conserved epitopes, which remain intact in evolutionary mechanisms. In the current work, we extended this epitope conservancy study to explore the interplay between individual cross-reactivity and the population dynamics of disease spread using the MCMC simulation approach.
We developed a simulation model to determine the effect of cross-reactivity on the dynamics of H1N1 infection in a population. The results were used to compare infection curves in the presence and absence of cross-reactivity, and when the initial size of the population with pre-existing immunity to the new strain (resulted from previous exposure to sH1N1) changes. We found that the infection curves of the asymptomatically and symptomatically infected populations undergo a very small change in peak magnitude and peak time. While pre-existing immunity can reduce the likelihood of clinical disease and mitigate illness, our findings suggest that the effects of T-cell immunity would be unlikely to be observed through only surveillance data collected for influenza-like illnesses or laboratory confirmed cases. This has important implications for public health planning and the development of targeted strategies (e.g., vaccination) for mitigating the impact of disease in the population. First, the incidence of an influenza epidemic may not provide a reliable measurement on the novelty of the virus or the level of pre-existing immunity in the population. Second, the effect of pre-existing cellular immunity in alleviating disease outcomes (e.g., hospitalization) remains largely undetectable in population dynamics of an epidemic. These suggest that rapid clinical investigations on T-cell assays remain crucial for determining the protection level conferred by pre-existing cellular responses in the face of an emerging influenza virus.
Our study focused on a fixed measurement of R 0, varying the transmission rate (β) and the initial proportion of the population (). A similar interesting question related to our current study is to determine how peak magnitude and time change if β is fixed and is varied, which will change R 0. This is a course for future work.
Our study focused on the effect of T-cell immunity developed from sH1N1 in providing partial protection against infection caused by a novel influenza strain. In the context of the 2009 H1N1 pandemic, previous work in primates has shown sH1N1 primed animals cleared the infection rapidly, suggesting the role of cross-reactive T-cell responses . In the absence of specific protective antibodies, epitope-specific CD4+ T cell immunity can help reduce viral shedding and mitigate illness, with evidence of cytolytic activity . In future work, combining the effects of both (cellular and humoral) types of immunity will be an important step in determining how immune system dynamics will impact the epidemic course, which may provide critical information that can be used to project plausible patterns of disease spread in the population at the early stages of disease outset.
Shinde V, Bridges C, Uyeki T, Shu B, Balish A, Xu X, Lindstrom Sea: Triple-reassortant swine influenza A (H1) in humans in the United States, 2005-2009. N Engl J Med. 2009, 360: 2616-2625. 10.1056/NEJMoa0903812.
Zimmer S, Burke D: Historical perspective–emergence of influenza A (H1N1) viruses. N Engl J Med. 2009, 361: 279-285. 10.1056/NEJMra0904322.
Itoh Y, Shinya K, Kiso M, Watanabe T, Sakoda Y, Hatta M, Muramoto Y, et al: In vitro and in vivo characterization of new swine-origin H1N1 influenza viruses. Nature. 2009, 460: 1021-1025.
Center for Disease Control: Update: Swine influenza A (H1N1) infections–california and texas april 2009. MMWR: Morbidity and Mortality Weekly Report. 2009, 58: 435-437.
Hancock K, Veguilla V, Lu X, Zhong W, Butler E, Sun H, Liu F, et al: Cross-reactive antibody responses to the 2009 pandemic H1N1 influenza virus. N Engl J Med. 2009, 361: 1945-1952. 10.1056/NEJMoa0906453.
Taubenberger J, Morens D: 1918 influenza: The mother of all pandemics. Emerg Infect Dis. 2006, 12: 15-22.
Dawood F, Iuliano A, Reed C, Meltzer M, Shay D, Cheng P, Bandaranayake D, Breiman R, Brooks W, Buchy P, Feikin D, Fowler K, Gordon A, Hien N, Horby P, Huang Q, Katz M, Krishnan A, Lal R, Montgomery J, Mïlbak K, Pebody R, Presanis A, Razuri H, Steens A, Tinoco Y, Wallinga J, Yu H, Vong S, Bresee J, Widdowson M: Estimated global mortality associated with the first 12 months of pandemic influenza A H1N1 virus circulation: a modelling study. Lancet Infect Dis. 2012, 12: 687-695. 10.1016/S1473-3099(12)70121-4.
Reed C, Angulo F, Swerdlow D, Lipsitch M, Meltzer M, Jernigan D, Finelli L: Estimates of the prevalence of pandemic (H1N1) 2009, united states, April-July 2009. Emerg Infect Dis. 2009, 15: 2004-2007. 10.3201/eid1512.091413.
Perez-Padilla R, de la Rosa-Zamboni D, Ponce de Leon S, Hernandez M, Quinones-Falconi F, Bautista E, Ramirez-Venegas A, et al: Pneumonia and respiratory failure from swine-origin influenza A (H1N1) in Mexico. N Engl J Med. 2009, 361: 680-689. 10.1056/NEJMoa0904252.
Tuite A, Greer A, Whelan M, Winter A, Lee B, Yan P, Wu J, Moghadas S, Buckeridge D, Pourbohloul B, Fisman D: Estimated epidemiologic parameters and morbidity associated with pandemic H1N1 influenza. CMAJ. 2010, 182: 131-136. 10.1503/cmaj.091807.
Duvvuri V, Moghadas S, Guo H, Duvvuri B, Heffernan J, Fisman D, Wu G, Wu J: Highly conserved cross-reactive CD4+ T-cell HA-epitopes of seasonal and the 2009 pandemic influenza viruses. Influenza ORV. 2010, 4: 249-258.
Greenbaum J, Kotturi M, Kim Y, Oseroff C, Vaughan K, Salimi N, Vita R, et al: Pre-existing immunity against swine-origin H1N1 influenza viruses in the general human population. PNAS. 2009, 106: 20365-20370. 10.1073/pnas.0911580106.
Wilkinson T, Li C, Chui C, Huang A, Perkins M, Liebner J, Lambkin-Williams R, Gilbert A, Oxford J, Nicholas B, Staples K, Dong T, Douek D, McMichael A, Xu X: Preexisting influenza-specific CD4+ T cells correlate with disease protection against influenza challenge in humans. Nat Med. 2012, 18: 274-280. 10.1038/nm.2612.
Ge X, Tan V, Bollyky P, Standifer N, James E, Kwok W: Assessment of seasonal influenza A virus-specific CD4 T-cell responses to 2009 pandemic H1N1 swine-origin influenza A virus. J. Virology. 2009, 84: 3312-3319.
Schanen B, De Groot A, Moise L, Ardito M, McClaine E, Martin W, Wittman V, Warren W, Drake D: Coupling sensitive in vitro and in silico techniques to assess cross-reactive CD4(+) T cells against the swine-origin H1N1 influenza virus. Vaccine. 2011, 29: 3299-3309. 10.1016/j.vaccine.2011.02.019.
EISS: Weekly Electronic Bulletin Week 2 :05/01/2009-11/01/2009 16. 2009, [http://www.euroflu.org/cgi-files/bulletin.cgi], January , Issue No 2879
McMichael A, Gotch F, Noble G, Beare P: Cytotoxic T-cell immunity to influenza. N Engl J Med. 1983, 309: 13-17. 10.1056/NEJM198307073090103.
Webby R, Andreansky S, Stambas J, et al: Protection and compensation in the influenza virus-specific CD8+ T cell response. PNAS. 2003, 100: 7235-7240. 10.1073/pnas.1232449100.
McElhaney J, Xie D, Hager W, et al: T cell responses are better correlates of vaccine protection in the elderly. J Immunol. 2006, 176: 6333-6339.
Pourbohloul B, Ahued A, Davoudi B, Meza R, Meyers L, Skowronski D, Villasenor I, Galvan F, Cravioto P, Earn D, Dushoff J, Fisman D, Edmunds W, Hupert N, Scarpino S, Trujillo J, Lutzow M, Morales J, Contreras A, Chavez C, Patrick D, Brunham R: Initial human transmission dynamics of the pandemic (H1N1) 2009 virus in North America. Influenza ORV. 2009, 3: 215-222.
Fraser C, Donnelly C, Cauchemez S, Hanage W, Van Kerkhove M, Hollingsworth T, Griffin J, Baggaley R, Jenkins H, Lyons E, Jombart T, Hinsley W, Grassly N, Balloux F, Ghani A, Ferguson N, Rambaut A, Pybus O, Lopez-Gatell H, Alpuche-Aranda C, Chapela I, Zavala E, Guevara D, Checchi F, Garcia E, Hugonnet S, Roth C: Pandemic potential of a strain of influenza A (H1N1): early findings. Science. 2009, 324: 1557-1561. 10.1126/science.1176062.
White L, Wallinga J, Finelli L, Reed C, Riley S, Lipsitch M, Pagano M: Estimation of the reproductive number and the serial interval in early phase of the influenza A/H1N1 pandemic in the USA. Influenza ORV. 2009, 3: 267-276.
Yang Y, Sugimoto J, Halloran M, Basta N, Chao D, Matrajt L, Potter G, Kenah E, Longini IJ: The transmissibility and control of pandemic influenza A (H1N1) virus. Science. 2009, 326: 729-733. 10.1126/science.1177373.
Longini IJ, Nizam A, Xu S, Ungchusak K, Hanshaoworakul W, Cummings D, Halloran M: Containing pandemic influenza at the source. Science. 2005, 309: 1083-1087. 10.1126/science.1115717.
Ferguson N, Cummings D, Cauchemez S, Fraser C, Riley S, Meeyai A, Iamsirithaworn S, Burke D: Strategies for containing an emerging influenza pandemic in Southeast Asia. Nature. 2005, 437: 209-214. 10.1038/nature04017.
Lessler J, Reich N, Brookmeyer R, Perl T, Nelson K, Cummings D: Incubation periods of acute respiratory viral infections: a systematic review. Lancet Infect Dis. 2009, 9: 291-300. 10.1016/S1473-3099(09)70069-6.
Alexander M, Bowman C, Feng Z, Gardam M, Moghadas S, Rost G, Wu J, Yan P: Emergence of drug-resistance: implications for antiviral control of pandemic influenza. Proc R Soc B. 2007, 274: 1675-1684. 10.1098/rspb.2007.0422.
Moghadas S, Bowman C, Rost G, Wu J: Population-wide emergence of antiviral resistance during pandemic influenza. PLoS One. 2008, 3: e1839-10.1371/journal.pone.0001839.
Heffernan J, Smith R, Wahl L: Perspectives on the basic reproductive ratio. J.R. Soc Interface. 2005, 2: 281-293. 10.1098/rsif.2005.0042.
van den Driessche P, Watmough J: Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2000, 180: 29-48.
McCaw J, McVernon J, McBryde E, Mathews J: Influenza: Accounting for Prior Immunity. Science. 2009, 325 (5944): 1071-[http://www.sciencemag.org/content/325/5944/1071.1.short]
Weinfurter J, Brunner K, Capuano S, Li C, Broman K, et al: Cross-Reactive T Cells Are Involved in Rapid Clearance of 2009 Pandemic H1N1 Influenza Virus in Nonhuman Primates. PLoS Pathogen. 2011, 7: e1002381-10.1371/journal.ppat.1002381.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/12/329/prepub
We thank James McCaw and Cameron Wellard for their insightful comments. This research was funded by NSERC (JMH, SMM, JW), MITACS (JMH, SMM, DNF, JW), MPRIME (JMH, SMM, JW), Early Researchers Award (JMH), OGS (BD), President Susan Mann Dissertation Scholarship (BD), CIHR (DNF,SMM) and the Arthritis Society (GW).
DF has received past research and educational funding from Sanofi Pasteur, Glaxo Smith Kline, and Novartis, all of which manufacture vaccines against influenza.
VD, SMM, JMH and DF were involved in the development of the model. VD, HG and JMH carried out the simulations. VD, JMH and SMM carried out the model analysis. BD participated in model analysis. VD, JMH, SMM wrote the mansucript. DF, BD, HG, JW and GW editted the mansucript. All authors read and approved the final manuscript.
About this article
Cite this article
Duvvuri, V.R., Heffernan, J.M., Moghadas, S.M. et al. The role of cellular immunity in Influenza H1N1 population dynamics. BMC Infect Dis 12, 329 (2012). https://doi.org/10.1186/1471-2334-12-329
- Prior Exposure
- Basic Reproduction Number
- H1N1 Pandemic
- Peak Magnitude