Variability in transmissibility of the 2009 H1N1 pandemic in Canadian communities
© Mostaço-Guidolin et al.; licensee BioMed Central Ltd. 2011
Received: 9 August 2011
Accepted: 13 December 2011
Published: 13 December 2011
The prevalence and severity of the 2009 H1N1 pandemic appeared to vary significantly across populations and geographic regions. We sought to investigate the variability in transmissibility of H1N1 pandemic in different health regions (including urban centres and remote, isolated communities) in the province of Manitoba, Canada.
The Richards model was used to fit to the daily number of laboratory-confirmed cases and estimate transmissibility (referred to as the basic reproduction number, R0), doubling times, and turning points of outbreaks in both spring and fall waves of the H1N1 pandemic in several health regions.
We observed considerable variation in R0 estimates ranging from 1.55 to 2.24, with confidence intervals ranging from 1.45 to 2.88, for an average generation time of 2.9 days, and shorter doubling times in some remote and isolated communities compared to urban centres, suggesting a more rapid spread of disease in these communities during the first wave. For the second wave, R e , the effective reproduction number, is estimated to be lower for remote and isolated communities; however, outbreaks appear to have been driven by somewhat higher transmissibility in urban centres.
There was considerable geographic variation in transmissibility of the 2009 pandemic outbreaks. While highlighting the importance of estimating R0 for informing health responses, the findings indicate that projecting the transmissibility for large-scale epidemics may not faithfully characterize the early spread of disease in remote and isolated communities.
Although the overall Canadian experience of the 2009 H1N1 pandemic influenza was perceived to be relatively mild for most individuals, the disease disproportionately affected several vulnerable and underserved communities, including First Nation reserves in northern Manitoba, remote and isolated communities in Nunavut, and Aboriginal communities on Vancouver Island [1–3]. The variability in disease outcomes is often described by factors pertinent to the population demographics, environmental characteristics, underlying health conditions, and healthcare access [3–6]. A direct effect of these factors appears on the transmissibility of disease in different population settings, which may contribute to larger outbreaks in vulnerable communities as a result of differential prevalence .
Transmissibility can be described as the basic reproduction number (R0), which in the epidemiological context, is defined as the average number of new cases generated by an infectious individual in a fully susceptible population . This quantity, once estimated for a new disease, can be used to identify target strategies for mitigating disease outcomes. Following the emergence of the H1N1 pandemic strain, attempts were made to estimate R0, largely for urban population settings [9–12]. These estimates suggest that the novel H1N1 strain was less transmissible (range: 1.2-1.7) than the 1918-1919 pandemic virus (range: 1.8-3.0) . However, H1N1 virus transmission appeared to vary significantly in remote and isolated communities, and among vulnerable population groups [1, 3], and therefore prior estimates of R0 for the urban population settings may not be applicable to devising targeted measures for future planning in these communities.
Given the experience of the 2009 pandemic in the province of Manitoba, Canada, with severe outbreaks in different health regions during both spring and fall waves, we sought to investigate the variability in transmissibility of the disease using the data for laboratory confirmed cases of H1N1 infection. The Richards model  was applied and fit to cumulative incidence of confirmed cases to estimate the transmissibility and turning points of the outbreaks in different population settings, and compared the results for the first and second waves. This modelling approach allows us to predict changes in the course of a single outbreak, determine the cumulative incidence at the turning point of the epidemic, and project the final number of identified cases (the cumulative incidence of infection) at the end of epidemic. In practical terms, determining the turning point of an outbreak will provide important public health information for timely implementation of effective intervention strategies. The results are compared with those of previous studies on transmissibility of pandemic H1N1 influenza. Finally, the findings are placed in the context of public health practice for future planning regarding the prevention and control of emerging infectious diseases, and the impact of constraints imposed by data and the model on the results presented in this study are discussed.
Materials and methods
To estimate plausible ranges for R0, the Richards model  was employed to fit to data collected for laboratory-confirmed cases of H1N1 infection in a single outbreak. This type of modelling, which describes the dynamics of cumulative incidence of infection, has been widely used in the study of plant disease epidemics [14, 15], and more recently for the epidemics in human populations [16–19].
was used, where τ is constant of integration related to the time for change in the growth rate of the cumulative cases. By fitting this solution to the cumulative case counts, the rate of infection spread (r) at the onset of each outbreak was obtained. Assuming that the generation time (the average time taken for the primary case to infect secondary cases ) follows a gamma distribution with shape parameter α and scale parameter β[23, 24], the reproduction number can be estimated as R0 = (1 + rβ) α using the moment generation method . The terminology of effective reproduction number, R e , is used for the second wave, given immunity from the first wave in the population. Turning points of the outbreaks (t m ) were also estimated by determining the inflection points at which I(t m ) = K/(1 + α)1/α. We should point out that previous work  has referred to the parameter τ as the turning point, which is theoretically incorrect; although in numerical terms it may lead to approximations that are close to the actual turning point. We also estimated the doubling time (the time it takes for the number of cases to double) at the early stages of the outbreak as t d = T ln 2(R0 - 1) .
Based on epidemiological studies conducted in the United States for the infector-infectee pairs in households , an average generation time of 2.9 days was considered with a standard deviation of 1.4 days for the Gamma distribution with α = 4.2 and β = 0.68 . We estimated R0 and R e using this generation time and also a shorter average generation time of 2.5 days with standard deviation of 0.9 days  that follows a Gamma distribution with shape and scale parameters α = 7.72 and β = 0.324, respectively.
A laboratory-confirmed case was defined as an individual with influenza-like illness or severe respiratory illness who tested positive for pandemic H1N1 influenza A virus by real-time reverse-transcriptase PCR (RT-PCR) or viral culture. The first case of H1N1 infection in Manitoba was identified (tested positive) on May 2, 2009, and data were reported by the earliest date of symptom onset, initial care, specimen collection, hospital admission, and/or intensive care unit admission. For the present study, data use was approved by the Human Research Ethics Board of the University of Manitoba (H2009:339) and Health Information Privacy Committee of Manitoba (2009/2010-40).
The solution of the Richards model was fit to the cumulative number of confirmed cases by applying the nonlinear least squares method using the R-language and environment for statistical computing. All fits were performed starting from the day on which the first laboratory confirmed case was identified in the associated health region. The turning point (t m ) of the outbreaks was determined using the Newton-Raphson method at the inflection point.
Estimates of fitting parameters and the basic reproduction number for the first wave of the 2009 pandemic in different health regions in the province of Manitoba
r (95% CI)
Estimates of fitting parameters and the effective reproduction number for the second wave of the 2009 pandemic in different health regions in the province of Manitoba
r (95% CI)
R e (95% CI)
Estimates for R0 and R e with 95% confidence intervals with the average generation time of 2.5 days
R e (95% CI)
During the second wave, Central, Brandon, and Assiniboine were the most impacted health regions in the southern part of the province of Manitoba. The mean effective reproduction number R e was estimated to be 1.63 (95% CI: 1.56-1.74) for Central, 1.74 (95% CI: 1.57-2.07) for Brandon, and 1.49 (95% CI: 1.46-1.52) for Assiniboine health regions with the average generation time of 2.9 days (Table 2). For comparison purposes, R e was also estimated for Interlake, Winnipeg, and for the entire province. We observed that R e estimates for Interlake and the entire province remained virtually the same for both waves, but slightly increased for Winnipeg in the second wave. Turning points of the outbreaks were projected between November 7 and 12. The results of parameter estimates are reported in Table 2, and best fits for Central, Winnipeg, and Manitoba are illustrated in Figure 2d-f.
The occurrence of pandemic outbreaks appeared to be associated with a geographic shift from northern Manitoba (e.g., Burntwood, North Eastman, and Norman) in the first wave to more southern parts of the province (e.g., Assiniboine, Central, and Brandon) in the second wave (Figure 1). Estimates for the first wave indicate that the plausible range of R0 for the entire province of Manitoba is virtually identical to that obtained for its largest urban centre (the Winnipeg health region). In the second wave, transmissibility appears to have been lower for the entire province than the Winnipeg health region. In comparing timelines of epidemic spread, our estimates of turning points show that the time interval for the pace of faster epidemic growth (the period before turning point) was comparable to that of the first wave, and occurred about 5 to 6 weeks after the onset of outbreaks in Manitoba. This is notwithstanding higher levels of pre-existing immunity and wider use of antiviral drugs in the second wave before vaccine distribution (which were expected to synergistically decelerate the rate of disease spread).
Discussion and conclusions
In the event of an emerging disease, estimates of transmissibility and other parameters pertinent to the nature of the infection provide critical information for guiding public health in optimizing health policy decisions and mitigating the impact of disease on affected populations . In the case of the 2009 pandemic, estimates of R0 helped characterize the epidemiological patterns of early disease spread in urban centres with large populations [9, 10, 12], and adapt preparedness measures to address the urgent situation in these settings [28, 29]. However, such estimates were not available for remote and isolated communities and other residential areas with small populations (e.g., Indigenous peoples on reserves), which experienced disproportionate rates of illness with poor outcomes, in particular during outbreaks in the spring wave. In Canada, severe outbreaks were observed in several locations in the province of Manitoba during the first wave, most notably in communities in the Burntwood, Norman and North Eastman health regions.
In contrast to previous work on estimating R0[9, 10, 27], the results of this study suggest a considerable geographic variation in transmissibility of the 2009 pandemic for the first wave, with estimates comparable to those of the 1918 (H1N1-Spanish) pandemic strain between 1.8 and 3.0 . For the second wave in Manitoba, our estimates indicate a reduction in R e with ranges that are considered to be of moderate transmissibility and comparable to those of the 1957 (H2N2) (range: 1.5-1.7) and 1968 (H3N2) pandemics (range: 1.6-2.2), yet significantly higher than the transmissibility of seasonal influenza (range: 1.1-1.4) [12, 30]. The higher estimates of R0 in the first wave of the 2009 pandemic are linked to outbreaks in remote and isolated communities, and populations of relatively small size; although estimates of the urban centre of the Winnipeg health region are close to the upper bound of previous estimates . Compared to estimates in the province of Ontario (range: 1.25-1.38) , and national Canadian estimates (range: 1.12-1.54) , the results of this study show that projecting the transmissibility using data for large-scale epidemics may not faithfully characterize the transmission pattern of early spread in rural areas and isolated communities.
The results of this study have several important implications for future planning to meet the threat of novel influenza viruses. Of greatest concern to public health are the geographic spread, severity, time course of an outbreak, and determining the most effective utilization of available social, preventive, and therapeutic resources to mitigate the impacts of an emerging disease. Understanding the relationship between differential severity and differential transmissibility of a disease in distinct populations can inform possible scenarios for allocation and optimal distribution of health resources prior to and during the spread of an infection to reduce vulnerability of the populations and alleviate disease outcomes. Disease management is also concerned with the transmissibility that is closely related to the final size of an epidemic in the long-term disease dynamics. The parameter R0 has a clear biological significance: it determines how fast the infection will spread through a population that has been previously unexposed to the disease, its meaning is independent of the model used to estimate its value, and it is dimensionless and thus directly comparable across populations and geographic regions. Our investigation in this study provides compelling evidence for the existence of strong links between transmissibility and differential severity of the disease, in particular in communities with limited access to healthcare. In small communities with strong social ties (e.g., Indigenous populations), frequent interactions between individuals are considerably higher than those in large urban centres, which can result in repeated exposure to infection and may account for a higher between-household transmission rate comparable to the rates of secondary household transmission observed at more southern latitudes. Given factors such as multigenerational households, high between-household social interactions, environmental characteristics, limited access to healthcare, and differential prevalence of predisposing health conditions and other types of health disparities, disease control strategies for remote and isolated communities may be significantly different from those applied to urban centres. Early real-time estimates of transmissibility will therefore be crucial for identifying the type and intensity of intervention measures that are effective in disease containment. In the absence of such information during the early stages of the 2009 pandemic, the existing recommendations for the use of antiviral drugs were modified to include the treatment of moderate to severe cases and individuals with pre-existing conditions or those at risk of developing poor outcomes, but to exclude treatment of mild cases or contacts . For future planning, a less conservative strategy may be beneficial for remote and isolated communities that are more likely to suffer from increased rates of disease transmission. A less conservative strategy might include a broader coverage of antiviral drugs with possible extension to prophylaxis of close contacts in these "high-risk" populations. Although not as cost-effective as vaccination, antiviral use for treatment and post-exposure prophylaxis is far more economical than hospitalization or intensive care. Furthermore, without a virus-specific vaccine (which, for example, could have been the case in the fall wave of the H1N1 pandemic had the herald wave not struck in the spring of 2009), antiviral medication may be the only pharmaceutical option and therefore availability and strategic use of drugs is crucial for mitigating disease in these vulnerable communities. With unique population characteristics that place some remote and isolated communities at increased risk for adverse health consequences , it is imperative for public health officials to identify transmission characteristics of early spread for the implementation of effective, feasible, and more economical health responses [32, 33].
This study has several limitations that warrant further investigation. For estimating the transmissibility of different health regions, we used data for laboratory confirmed cases of H1N1 infection, which may introduce some bias due to differential rates of testing across age groups, over time, and possibly between different health regions. It is difficult to assess the magnitude and direction of such biases, but we understand that they influence estimates and comparative analysis of R0 within and between health regions. While our results demonstrate the variability in diseases transmissibility in different population settings, identification of factors responsible for such variation remains an important objective for future work, and in the context of the 2009 pandemic, understanding the role of these factors merit further population-specific research. For example, we used estimates of generation times from two previous study, which appear to be higher than the early estimate of 1.9 days for generation times based on the La Gloria data from initial outbreaks in Mexico . These estimates may have been slightly inflated by the inclusion of a small number of imported cases. Nevertheless, for a given generation time, estimates of R0 for remote and isolated communities will remain higher than those for more urban regions. However, the generation time is a key parameter in our study, and we note that it may vary between different populations with distinct demographic variables, which will in turn affect estimates of the reproduction number. Furthermore, R0 can be calculated as the product of the frequency of effective contact (leading to transmission) and the average infectious period. It is possible that in remote communities, due to the limited access to health services and medical care, an average infectious individual may experience a longer contagious period than an average infectious person in urban centres, hence contributing to a greater reproduction number without having to make more frequent contact per unit time.
The model employed here does not include the immune status of the populations; yet we understand that estimates of R0 for the first wave may be affected by pre-existing immunity [34, 35], and therefore the actual value of R0 (when the entire population is susceptible) may be greater than estimates reported here. For the second wave, the effect of pre-existing immunity becomes more pronounced by the reduced level of susceptibility due to infection in the spring wave in northern communities. Furthermore, unlike the first wave, antiviral drugs were used earlier and in a much wider scale . There were also some limitations in data for laboratory confirmed cases of the second wave. For instance, in order to maintain the laboratory capacity and responsiveness, laboratory testing of non-severe cases was suspended on November 20, 2009, unless infected individuals met indications outlined in the Manitoba Health Clinical Care Guidelines for Pandemic H1N1. Although this modification was made well beyond the establishment phase of the second wave, it may have resulted in a significant under-reporting at the peak of pandemic outbreaks, and therefore estimates of the final size of the cumulative number of confirmed cases and other fitting parameters may be influenced. Despite these limitations, the findings of this study underscore the importance of rapid identification of transmission characteristics in the early stages of an emerging disease to identify effective and population specific intervention strategies, and optimize health responses.
The work was supported by the Canadian Institutes of Health Research Catalyst Grant for Pandemic Preparedness (Grant No.: 104047), the Mathematics of Information Technology and Complex Systems, and the Natural Sciences and Engineering Research Council of Canada. The funders had no role in study design, data collection and analysis, preparation of the manuscript, or decision to publish. The authors would like to thank Dr. Shelly Deeks and Steven Johnson for their contributions to the manuscript and the reviewers for their insightful comments that have improved this paper.
- Zarychanski R, Stuart TL, Kumar A, Doucette S, Elliott L, Kettner J, Plummer F: Correlates of severe disease in patients with 2009 pandemic influenza (H1N1) virus infection. CMAJ. 2010, 182 (3): 257-64. 10.1503/cmaj.091884.PubMedPubMed CentralView ArticleGoogle Scholar
- Kondro W: Dispensing antivirals in underserved communities. CMAJ. 2009, 181 (9): E199-200. 10.1503/cmaj.109-3055.PubMedPubMed CentralView ArticleGoogle Scholar
- Campbell A, Rodin R, Kropp R, Mao Y, Hong Z, Vachon J, Spika J, Pelletier L: Risk of severe outcomes among patients admitted to hospital with pandemic (H1N1) influenza. CMAJ. 2010, 182 (4): 349-55. 10.1503/cmaj.091823.PubMedPubMed CentralView ArticleGoogle Scholar
- To K, Wong S, Li I, Hung I, Tse H, Woo P, Chan KH, Yuen KY: Concurrent comparison of epidemiology, clinical presentation and outcome between adult patients suffering from the pandemic influenza A (H1N1) 2009 virus and the seasonal influenza A virus infection. Postgrad Med J. 2010, 86 (1019): 515-21. 10.1136/pgmj.2009.096206.PubMedView ArticleGoogle Scholar
- Bautista E, Chotpitayasunondh T, Gao Z, Harper Sa, Shaw M, Uyeki TM, Zaki SR, Hayden FG, Hui DS, Kettner JD, Kumar A, Lim M, Shindo N, Penn C, Nicholson KG: Clinical aspects of pandemic 2009 influenza A (H1N1) virus infection. New Engl J Med. 2010, 362 (18): 1708-19.PubMedView ArticleGoogle Scholar
- La Ruche G, Tarantola A, Barboza P, Vaillant L, Gueguen J, Gastellu-Etchegorry M: The 2009 pandemic H1N1 influenza and indigenous populations of the Americas and the Pacific. Euro Surveill. 2009, 14 (42): 1-6.Google Scholar
- World Health Organization: [http://www.who.int/environmental_health_emergencies/disease_outbreaks/communicable_diseases/en/index.html]
- Heesterbeek OD, P JA: Mathematical Epidemiology of Infectious Diseases. John Wiley and Sons Ltd. 2000Google Scholar
- Fraser C, Donnelly CA, Cauchemez S, Hanage WP, Van Kerkhove MD, Hollingsworth TD, Griffin J, Baggaley RF, Jenkins HE, Lyons EJ, Jombart T, Hinsley WR, Grassly NC, Balloux F, Ghani AC, Ferguson NM, Rambaut A, Pybus OG, Lopez-Gatell H, Alpuche-Aranda CM, Chapela IB, Zavala EP, Guevara DME, Checchi F, Garcia E, Hugonnet S, Roth C: Pandemic potential of a strain of influenza A (H1N1): early findings. Science. 2009, 324 (5934): 1557-61. 10.1126/science.1176062.PubMedPubMed CentralView ArticleGoogle Scholar
- Pourbohloul B, Ahued A, Davoudi B, Meza R, Meyers La, Skowronski DM, Villaseñor I, Galván F, Cravioto P, Earn DJD, Dushoff J, Fisman D, Edmunds WJ, Hupert N, Scarpino SV, Trujillo J, Lutzow M, Morales J, Contreras A, Chávez C, Patrick DM, Brunham RC: Initial human transmission dynamics of the pandemic (H1N1) 2009 virus in North America. Influenza Other Respi Viruses. 2009, 3 (5): 215-22. 10.1111/j.1750-2659.2009.00100.x.PubMed CentralView ArticleGoogle Scholar
- White LF, 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 2009 influenza A/H1N1 pandemic in the USA. Influenza Other Respi Viruses. 2009, 3 (6): 267-76. 10.1111/j.1750-2659.2009.00106.x.PubMed CentralView ArticleGoogle Scholar
- Yang Y, Sugimoto JD, Halloran ME, Basta NE, Chao DL, Matrajt L, Potter G, Kenah E, Longini IM: The transmissibility and control of pandemic influenza A (H1N1) virus. Science. 2009, 326 (5953): 729-33. 10.1126/science.1177373.PubMedPubMed CentralView ArticleGoogle Scholar
- Turner ME, Bradley EL, Kirk Ka, Pruitt KM: A theory of growth. Math Biosci. 1976, 29 (3-4): 367-373. 10.1016/0025-5564(76)90112-7.View ArticleGoogle Scholar
- Segarra J, Jeger MJ, van den Bosch F: Epidemic dynamics and patterns of plant diseases. Phytopathology. 2001, 91 (10): 1001-10. 10.1094/PHYTO.2001.91.10.1001.PubMedView ArticleGoogle Scholar
- Hau B, Kosman E: Comparative analysis of flexible two-parameter models of plant disease epidemics. Phytopathology. 2007, 97 (10): 1231-44. 10.1094/PHYTO-97-10-1231.PubMedView ArticleGoogle Scholar
- Zhou G, Yan G: Severe acute respiratory syndrome epidemic in Asia. Emerg Infect Dis. 2003, 9 (12): 1608-[http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3034341/]PubMedPubMed CentralGoogle Scholar
- Hsieh YH, Chen CWS: Turning points, reproduction number, and impact of climatological events for multi-wave dengue outbreaks. Trop Med Int Health. 2009, 14 (6): 628-38. 10.1111/j.1365-3156.2009.02277.x.PubMedView ArticleGoogle Scholar
- Hsieh YH, Ma S: Intervention measures, turning point, and reproduction number for dengue, Singapore, 2005. Am J Trop Med Hyg. 2009, 80: 66-71. [http://www.ncbi.nlm.nih.gov/pubmed/19141842]PubMedGoogle Scholar
- Hsieh YH, Fisman DN, Wu J: On epidemic modeling in real time: an application to the 2009 Novel A (H1N1) influenza outbreak in Canada. BMC Res Notes. 2010, 3: 283-10.1186/1756-0500-3-283.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsoularis a, Wallace J: Analysis of logistic growth models. Math Biosci. 2001, 179: 21-55.View ArticleGoogle Scholar
- Hsieh YH: Pandemic influenza A (H1N1) during winter influenza season in the southern hemisphere. Influenza Other Respi Viruses. 2010, 4 (4): 187-97. 10.1111/j.1750-2659.2010.00147.x.View ArticleGoogle Scholar
- Anderson RM, Fraser C, Ghani AC, Donnelly Ca, Riley S, Ferguson NM, Leung GM, Lam TH, Hedley AJ: Epidemiology, transmission dynamics and control of SARS: the 2002-2003 epidemic. Philos Trans R Soc Lond Series B Biol Sci. 2004, 359 (1447): 1091-105. 10.1098/rstb.2004.1490. [http://www.ncbi.nlm.nih.gov/pubmed/15306395]View ArticleGoogle Scholar
- Donnelly Ca, Finelli L, Cauchemez S, Olsen SJ, Doshi S, Jackson ML, Kennedy ED, Kamimoto L, Marchbanks TL, Morgan OW, Patel M, Swerdlow DL, Ferguson NM: Serial intervals and the temporal distribution of secondary infections within households of 2009 pandemic influenza A (H1N1): implications for influenza control recommendations. Clin Infect Dis. 2010, 52 (Supplement 1): S123-S130.View ArticleGoogle Scholar
- McBryde E, Bergeri I, Gemert CV: Early transmission characteristics of influenza A (H1N1) v in Australia: Victorian state, 16 May-3 June 2009. Eurosurveillance. 2009, 14 (42): 1-6.Google Scholar
- Wallinga J, Lipsitch M: How generation intervals shape the relationship between growth rates and reproductive numbers. Proc R Soc B: Biol Sci. 2007, 274 (1609): 599-604. 10.1098/rspb.2006.3754.View ArticleGoogle Scholar
- Hahné S, Donker T, Meijer A, Timen A, Van Steenbergen J, Osterhaus A, Van Der Sande M, Koopmans M, Wallinga J, Coutinho R, Others: Epidemiology and control of influenza A(H1N1)v in the Netherlands: the first 115 cases. Euro Surveill. 2009, 14 (27): 1-4.Google Scholar
- Tuite AR, Greer AL, Whelan M, Winter AL, Lee B, Yan P, Wu J, Moghadas S, Buckeridge D, Pourbohloul B, Fisman DN: Estimated epidemiologic parameters and morbidity associated with pandemic H1N1 influenza. CMAJ. 2010, 182 (2): 131-6. 10.1503/cmaj.091807.PubMedPubMed CentralView ArticleGoogle Scholar
- Ajelli M, Merler S, Pugliese a, Rizzo C: Model predictions and evaluation of possible control strategies for the 2009 A/H1N1v influenza pandemic in Italy. Epidemiol Infect. 2010, 139 (June 2009): 1-12.Google Scholar
- Halder N, Kelso JK, Milne GJ: Analysis of the effectiveness of interventions used during the 2009 A/H1N1 influenza pandemic. BMC Public Health. 2010, 10: 168-10.1186/1471-2458-10-168.PubMedPubMed CentralView ArticleGoogle Scholar
- Viboud C, Tam T, Fleming D, Handel A, Miller Ma, Simonsen L: Transmissibility and mortality impact of epidemic and pandemic influenza, with emphasis on the unusually deadly 1951 epidemic. Vaccine. 2006, 24: (44-46):6701-7.View ArticleGoogle Scholar
- Public Health Agency of Canada: Clinical recommendations for patients presenting with respiratory symptoms during the 2009-2010 influenza season. [http://www.phac-aspc.gc.ca/alert-alerte/h1n1/pdf/H1N1/Decision/Tree/oct23e.pdf]
- Moghadas SM, Pizzi NJ, Yan P, Wu J, Tamblyn DF SE: Canada in the face of the 2009 H1N1 pandemic. Influenza Sother Respi Viruses. 2011, 5: 83-88. 10.1111/j.1750-2659.2010.00184.x.View ArticleGoogle Scholar
- Laskowski M, Mostaço-Guidolin LC, Greer AL, Wu J, Moghadas SM: The impact of demographic variables on disease spread: influenza in remote communities. Nat: Sci Rep. 2011, 1 (i): 1-7.Google Scholar
- Greenbaum Ja, Kotturi MF, Kim Y, Oseroff C, Vaughan K, Salimi N, Vita R, Ponomarenko J, Scheuermann RH, Sette A, Peters B: Pre-existing immunity against swine-origin H1N1 influenza viruses in the general human population. Proc Natl Acad Sci USA. 2009, 106 (48): 20365-70. 10.1073/pnas.0911580106.PubMedPubMed CentralView ArticleGoogle Scholar
- Duvvuri VRSK, Moghadas SM, Guo H, Duvvuri B, Heffernan JM, Fisman DN, Wu GE, Wu J: Highly conserved cross-reactive CD4+ T-cell HA-epitopes of seasonal and the 2009 pandemic influenza viruses. Influenza Other Respi Viruses. 2010, 4 (5): 249-58. 10.1111/j.1750-2659.2010.00161.x.View ArticleGoogle Scholar