- Research article
- Open Access
- Open Peer Review
Effect modification of environmental factors on influenza-associated mortality: a time-series study in two Chinese cities
BMC Infectious Diseasesvolume 11, Article number: 342 (2011)
Environmental factors have been associated with transmission and survival of influenza viruses but no studies have ever explored the role of environmental factors on severity of influenza infection.
We applied a Poisson regression model to the mortality data of two Chinese metropolitan cities located within the subtropical zone, to calculate the influenza associated excess mortality risks during the periods with different levels of temperature and humidity.
The results showed that high absolute humidity (measured by vapor pressure) was significantly (p < 0.05) associated with increased risks of all-cause and cardiorespiratory deaths, but not with increased risks of pneumonia and influenza deaths. The association between absolute humidity and mortality risks was found consistent among the two cities. An increasing pattern of influenza associated mortality risks was also found across the strata of low to high relative humidity, but the results were less consistent for temperature.
These findings highlight the need for people with chronic cardiovascular and respiratory diseases to take extra caution against influenza during hot and humid days in the subtropics and tropics.
Influenza used to be considered as a "cold" disease as it usually returns every cold winter in temperate countries. However, recent studies have shown that influenza can be active throughout the year in the warm tropics and subtropics and the disease burden of influenza there can be as heavy as that in temperate climates [1–4]. It has been proposed that influenza seasonality is driven by complicated interactions between antigenic drifts of virus strains, environmental factors, host susceptibility and behavior changes [5–7]. Among these potential factors, environmental factors including temperature and relative humidity have been most thoroughly explored by laboratory and observational studies [8, 9]. Recent studies raised a hypothesis that absolute humidity is one of drivers for influenza seasonality in temperate regions [10, 11], but the mechanism behind various seasonal patterns of influenza outbreaks under different climates remains unclear.
Given the association of environmental factors with both influenza virus activity and mortality, these factors have been adjusted for as confounders in the statistical models for influenza associated mortality burden [2, 12]. However, few of previous laboratory or epidemiological studies tackled the potential modifying role of environmental factors on severity of infection, or mortality risks associated with influenza. Such effect modification is plausible because extreme weather conditions could have a synergistic effect with influenza on mortality burden. Both cold and hot temperatures have been associated with increased mortality risks of respiratory diseases in numerous studies [13–15]. Our previous study in Hong Kong also found a two-peak seasonal variation in mortality burden of influenza, which is similar to the pattern of influenza seasonality, suggesting that environmental factors might also affect the severity of seasonal influenza infection . In this study we applied a Poisson regression model to the data of two subtropical Chinese cities: Guangzhou and Hong Kong, to examine the possible effect modification of environmental factors on the severity of influenza infection measured by the mortality attributable to influenza. These two cities are geographically close, with Guangzhou located at latitude 23°N and Hong Kong at 21°N (Figure 1). Both cities have a typical subtropical climate, but on average Hong Kong has a higher temperature and humidity. Hong Kong also has a larger population than Guangzhou (6.8 million vs. 3.7 million).We examined three environmental factors which have been documented to regulate virus survival and transmission: temperature, relative humidity and absolute humidity.
Influenza surveillance data were obtained from the Hong Kong Department of Health (DH) for the period of 1998-2006 in Hong Kong and the Guang Dong Provincial Centre for Disease Control and Prevention (GDCDC) for 2004-2006 in Guangzhou . In both cities the consultation rates for influenza-like illness (ILI) from sentinel hospitals or doctors were reported on the weekly basis, with diagnosis criteria of fever (≥38°C), cough or sore throat. Nasopharyngeal aspirate or swab was collected from ILI patients and then tested of influenza viruses by immunofluorescence and cell culture. Positive specimens were further typed (subtyped) into A/H3N2, A/H1N1 and B by haemagglutination inhibition test with strain-specific antisera provided by the World Health Organization. During 2004-2006, annual averages of 35,256 and 1,872 specimens were tested by the DH and GDCDC, respectively.
Mortality data for the period of 2004-2006 were obtained from the Hong Kong Census and Statistics Department and the Guangzhou Department of Health. According to Hong Kong Law all deaths from natural causes are required to be registered, so the mortality data in Hong Kong was fairly complete. The registered death data collected by the Guangzhou Department of Health covered the residents from the eight urban districts, but did not include the immigrant population and residents living in the two suburban districts. Therefore, the population denominator of Guangzhou did not change much during the study period. Given the fact that most immigrant workers are young adults and few died in Guangzhou, we can assume that our mortality data was complete and representative of the general population of Guangzhou. We aggregated the weekly numbers of deaths with underlying cause of cardiorespiratory disease, pneumonia and influenza, and all-cause mortality. The corresponding International Classification of Diseases, Tenth Revision (ICD10) codes adopted by Guangzhou and Hong Kong are I00-I99, J00-J99 for cardiorespiratory, J10-J18 for pneumonia and influenza, and A00-R99 for all-cause deaths. We used the accidental deaths (ICD10 codes S00-T989) in Hong Kong as the control disease. Weekly mean temperature and relative humidity were separately derived from the National Meteorological Information Centre of China and Hong Kong Observatory for Guangzhou and Hong Kong respectively. Weekly mean vapor pressure was calculated as a metric for absolute humidity following an equation provided by Basu et al. ,
Vapor pressure will be used for absolute humidity hereafter in this paper.
A generalized additive model (GAM) with a log link function and Poisson error was fitted to the weekly numbers of deaths. The Poisson model has been widely applied to the estimation of influenza associated disease burden and recently been validated by an empirical dataset of laboratory confirmed influenza cases . First, long-term trends and seasonal patterns of cause-specific mortality counts, as well as environmental factors including temperature and relative humidity, were adjusted for as confounders by building a core model:
Here y t denotes the number of deaths at week t. ns(t), ns(temp t ), and ns(humd t ) denote the natural cubic spline smoothing functions of time, weekly average temperature and relative humidity. The natural spline smoothing function was used to remove the small variation while maintaining the major trend of each variable, i.e. to make them smoother. The aim of smoothing was to increase the efficiency in estimating the model coefficients . Given the high correlation between temperature and vapor pressure, the smoothing function of vapor pressure was not added in order to avoid collinearity between the variables. The adequacy of this core model was evaluated by the absence of any obvious pattern in partial autocorrelation functions of its residuals (Additional file 1: Figure S1). Second, the weekly proportions of specimens positive for influenza A or B were then added into the core model as a variable for influenza virus activity to obtain a main effect model.
The main effects of influenza have been presented elsewhere . To explore the effect modification of environmental factors on influenza associated mortality, we added into the core model the product terms of influenza proportion variable and dummy variables for periods of normal and extreme (high or low) weather conditions as interaction terms between virus activity and environmental factors. For example, the interaction model for temperature and influenza mortality was:
where Lowtemp i = 1 for the periods within low temperature ranges and 0 for otherwise, and Midtemp t and Hightemp t are similarly defined as the dummy variables for the middle and high temperature periods. The smoothing function of temperature ns(temp t ) was kept in the model in order to adjust for the association of temperature and mortality. It is reasonable to assume that the cutoff point of extreme weather may differ across cities as people may adapt to prevailing climates. We used the first (25%) and third quartiles (75%) of weekly average temperatures (or humidity) as the cutoff points to define the low, middle and high temperature (or humidity) periods . The presence of effect modification by environmental factors was evaluated using likelihood ratio tests between the interaction and main effect models. To measure the effects of influenza on mortality, we computed the percentage change of mortality counts associated with 1% increases of influenza virus activity for the low, middle and high periods. The formula for the low temperature period is
Here β1 was obtained from the above interaction model. To test whether our results were robust to various definitions of periods, we chose two sets of extra cutoff points: 20th and 80th, 30th and 70th percentiles of weekly average data in each city. Although we removed the autocorrelation and seasonal trends within mortality data, there are still concerns that such adjustment was inadequate and the remaining uncontrolled seasonal factors may cause interaction terms of environmental factors and virus activity to appear significant in our models. To rule out this possibility, we used accidental deaths that were expected to be unrelated to influenza infection as a control mortality group to show that our findings are unlikely the spurious results of under-adjustment of seasonal confounding factors in modeling. Some previous studies used the anomalies to assess the effects of meteorology factors on influenza .We therefore conducted a sensitivity analysis by defining the strata by anomalies, instead of absolute values of meteorology factors. The anomalies were defined as the deviations of observed metrological data from a seasonal curve with a constant and a sinusoidal pair fitted to these observed data. All the analyses were performed using the mgcv package of R software (version 2.5.1.) .
Hong Kong has a larger population than Guangzhou (6.8 million versus 3.7 million) during our study period. On average, Hong Kong has a slightly higher temperature, relative humidity and vapor pressure, and smaller annual variations than Guangzhou (Table 1). The mean of the weekly proportions of specimens positive for influenza A or B was higher in Hong Kong than in Guangzhou. In 2004 and 2006, there were two peaks of virus activity in Hong Kong (one in February/March and another in June/July), but only one broad peak in Guangzhou (Figure 2). The influenza seasonality was similar between these two cities in 2005.
More deaths with underlying causes of cardiorespiratory, pneumonia and influenza or all-cause were recorded in the low temperature (or low vapor pressure) periods for both Guangzhou and Hong Kong (Table 2). For both cities, the mortality counts did not show any obvious difference across levels of relative humidity. High influenza virus activity coincided with high levels of temperature, relative humidity or vapor pressure, with the only exception of Guangzhou which had higher mean proportions when temperature is within the middle range.
Significant interaction between temperature and influenza on the mortality risks was only found in cardiorespiratory mortality in Hong Kong (p < 0.05), and the interaction between relative humidity and influenza was found significant for all-cause mortality in both Guangzhou and Hong Kong, and for cardiorespiratory mortality only in Guangzhou (Table 3). The patterns of influenza impact across the low to high temperature periods were not consistent among the two cities. In Guangzhou, the highest risk of mortality tended to be observed in the middle-temperature period, whereas for Hong Kong the largest changes in risk were found in the periods with high temperature for the three mortality categories (Table 3). An increasing pattern of influenza associated mortality risks could be observed along low-, mid- and high-relative humidity periods in Guangzhou and Hong Kong, but most of estimates for low- and high-periods were not statistically significant (Table 3).
The interaction between vapor pressure levels and virus activity were found to be significant (p < 0.05) for all-cause and cardiorespiratory mortality, but not for pneumonia and influenza mortality (Table 3). Consistently higher mortality risks were found at the high levels of vapor pressure, with only exception of P&I mortality in Hong Kong. The influenza effects on mortality of cardiorespiratory, pneumonia and influenza, and all causes were found significant during the middle- and high-vapor pressure periods (with the only exception of all-cause mortality in Hong Kong), but not significant when vapor pressures remained at the relatively low levels. For the high vapor pressure periods, all-cause excess mortality counts attributable to influenza would increase by 0.35% and 0.26% for per 1% increase of virus activity, and the corresponding increases for cardiorespiratory mortality were 0.54% and 0.49% for Guangzhou and Hong Kong, respectively (Table 3).
We also assessed the effect modification of environmental factors on influenza effects (i.e. interaction of each factor and influenza) for the age group younger than 65 years (< 65) and the elderly aged 65 years or older (≥65). The results of the ≥65 age group were consistent with those for the all-ages group, with an increasing trend over vapor pressure levels observed for all-cause and cardiorespiratory mortality in the two cities (Figure 3). For the < 65 age group, this trend could also be observed in most city-specific disease categories, with the only exception of all-cause mortality in Hong Kong. But the interaction terms were statistically significant (p < 0.05) only in the ≥65 group. The modifying effects of temperature and relative humidity on influenza effects were found quite similar between the all-ages and ≥65 age groups, but not between the all-ages and < 65 age groups (data not shown).
The models with indicators for different cutoff points returned similar estimates for temperature, relative humidity and vapor pressure (Additional file 1: Figure S2). We did not observe any significant interaction between influenza and environmental factors in terms of their effects on the control mortality category of accidental mortality in Hong Kong. Influenza associated accidental mortality risks were also not statistically significant.
To ensure the standardized comparison between Guangzhou and Hong Kong, we decided to apply the same modeling approach to the data of same study period, because the core model would have been slightly different if we used the longer time series of Hong Kong. To check the robustness of our conclusions, we repeated the above analysis using a longer time series data of Hong Kong during 1998-2006 and the estimates were shown in Additional file 1: Table S1. The statistical significance of interaction terms was close to those from the study period of 2004-2006. The estimates for the low temperature periods became larger and statistically significant; those for the high relative humidity periods were smaller and comparable to the middle relative humidity; and for the low vapor pressure periods, the estimates were similar but became significant. The increasing trend across the low, middle and high levels of temperature and relative humidity was less evident.
The results of stratification analysis by anomalies are shown in Additional file 1: Table S2. The estimates were similar to those for the periods defined by absolute values of meteorological factors, in terms of magnitude and changing patterns. But the likelihood ratio tests showed more significant interaction for temperature or vapor pressure, and less for relative humidity.
In this study we quantified influenza associated mortality risks at the various ranges of temperature and humidity, and compared the results between two large subtropical cities. An increasing pattern of influenza-associated mortality risks on all-cause and cardiorespiratory along the low to high periods was found for temperature, relative and absolute humidity in both cities during the period of 2004-2006. The interaction between vapor pressure indicators and virus activity were also consistently significant for relative and absolute humidity. Although lower vapor pressure (or temperature) has been found to facilitate virus transmission and survival in the guinea pig model , our results suggested that higher vapor pressure (or temperature) was associated with a higher mortality risk attributable to influenza. Severity of seasonal influenza epidemics were not only determined by virus transmission efficiency and outdoor weather conditions, but also largely affected by host resistance, indoor living environment and social behavior . The high vapor pressure periods coincided with the summer peaks of influenza in both Guangzhou and Hong Kong, indicating that influenza could pose a higher risk when reaching its peak in seasons with high vapor pressure in subtropical cities.
Our results may help to interpret our previous findings that the effects of influenza were significantly higher in the humid and warm spring/summer period than in the dry and cold winter period in Hong Kong . The extreme low vapor pressure was usually recorded during December-January and the highest appeared during June-July, which coincided with the trough and peak periods of excess risks associated with influenza viruses. Effect modification of environmental factors was only detected in all-cause and cardiorespiratory mortality, but not in pneumonia and influenza, suggesting that the synergistic interaction between high humidity (or temperature) and virus activity may mainly lie in their similar regulation pathways in cardiovascular systems. These results are also in agreement with our previous findings that pneumonia and influenza mortality risks attributable to influenza did not exhibit a seasonal variation . Extreme heat has been documented to increase blood viscosity through evaporation of body fluid and trigger intravascular coagulation through damaging endothelial cells . Influenza infection has a similar pro-thrombotic effect by inducing inflammation around blood vessels and rupturing atherosclerotic plaques . As a result, mortality risks would be dramatically raised by the stress of both extreme weather and influenza infections. The results suggested that we need extra precautionary measures to reduce influenza infections in people with cardiovascular diseases, especially under the frequent hot and humid weather conditions experienced in the tropical and subtropical areas. Although the experiments of influenza virus transmission between guinea pig hosts found the higher transmission rates occurred under dry air (vapor pressure below 10hPA) , our results indicated that the mortality risks associated with influenza under the low vapor pressure environment were lower than the rest of study period. The reason could be the short time of exposure to the very low level of vapor pressure. During our study period there were only 2 and 13 weeks with an average vapor pressure below 10hPA in Hong Kong and Guangzhou, respectively. Most weekly average vapor pressures during the low vapor pressure period were within the range of 10-20hPA, in which the guinea pig experiments showed dramatically reduced transmission rates . In future, we may examine the seasonal variation in influenza effects in other cities to assess whether such a seasonal variation, if common in subtropical and tropical cities, is consistently determined by environmental factors, or by other factors such as host immunity and virus virulence.
The results for the different age groups suggested that the modification effects of environmental factors may mainly lie in the elderly aged over 65 years, as the consistent increasing trend over the low to high vapor pressure periods was only observed in this age group. However, since over 65% of deaths occurred in this age group for both cities (Table 1), the small numbers of weekly death counts in the younger age group (< 65 years) may not have had enough power to allow assessment of effect modification based on the data over the 3 years. A future study with a long study period or a large population may help answer whether young people maybe also expose to higher mortality risks during the hot and humid days.
In this study, we used the quartiles of weekly data in each city, to separately define the periods with normal (middle) and extreme (low and high) weather. Given the difference in weather conditions between these cities, we think that it is not appropriate to use the same cutoff points for temperature or humidity to compare their modification effects on influenza associated mortality, as people living in hot subtropical and tropical regions may adapt well to the year-round hot and humid climate and have a higher threshold for adverse effects of weather. For example, although it is widely accepted that the temperature effects on mortality exhibited a U- or V-shape curve in both temperate and tropical/subtropical areas, the turning point of this curve varied across different cities. A study conducted in 11 cities of the US found that the turning point of temperature for its effects on mortality could range from 18.4°C to 32.4°C . To our best knowledge, so far there are no studies that have ever assessed the effect modification of temperature on influenza effects. Therefore, the commonly adopted cutoff points of city-specific quartiles seem appropriate at this stage [21, 27].
There are several limitations in our study. Firstly, our study is based on 3 years of surveillance data which may not have enough power to allow assessment of exposure-response curves for the effects of environmental factors. Nevertheless, our findings did suggest an increasing trend of influenza associated mortality risks across the periods of low, middle and high vapor pressure, although such findings may be applicable only to the warm climates. Secondly, we only investigated the effect modification of environmental factors through a simple interaction model, but there were other unadjusted factors, including host susceptibility and virulence of influenza strains. These factors are unlikely to work independently with environmental factors. Other environmental factors such as ultraviolet radiation , rainfall  have been proposed to play a role in the regulation of influenza seasonality, although evidence is rather limited compared with the three factors we chose to investigate . Lastly, we did not adjust for the vaccination rate in our model. In 2003, vaccination rate was 191 doses per 1,000 total population in Hong Kong , slightly higher than the rate of 129 doses/1,000 total population in Guangzhou . However, it is not clear when people received vaccination; therefore we were unable to assess the role of vaccination in our study.
This study provides a piece of key evidence to the effect of environmental factors on severity of seasonal influenza under warm climates and helps reveal the mechanism behind global influenza seasonality. It also highlights the need for people with chronic cardiovascular and respiratory conditions to take extra caution against influenza during the hot and humid days in the subtropics.
Chow A, Ma S, Ling AE, Chew SK: Influenza-associated deaths in tropical Singapore. Emerg Infect Dis. 2006, 12: 114-121.
Wong CM, Chan KP, Hedley AJ, Peiris JSM: Influenza-associated mortality in Hong Kong. Clin Infect Dis. 2004, 39: 1611-1617. 10.1086/425315.
Wong CM, Yang L, Chan KP, Leung GM, Chan KH, Guan Y, et al: Influenza-associated hospitalization in a subtropical city. PLoS Med. 2006, 3: e121-10.1371/journal.pmed.0030121.
Viboud C, Alonso WJ, Simonsen L: Influenza in tropical regions. PLoS Med. 2006, 3: e89-10.1371/journal.pmed.0030089.
Dowell SF: Seasonal variation in host susceptibility and cycles of certain infectious diseases. Emerg Infect Dis. 2001, 7: 369-374.
Dowell SF, Ho MS: Seasonality of infectious diseases and severe acute respiratory syndrome-what we don't know can hurt us. Lancet Infect Dis. 2004, 4: 704-708. 10.1016/S1473-3099(04)01177-6.
Lofgren E, Fefferman NH, Naumov YN, Gorski J, Naumova EN: Influenza seasonality: underlying causes and modeling theories. J Virol. 2007, 81: 5429-5436. 10.1128/JVI.01680-06.
Schulman JL, Kilbourne ED: Airborne transmission of influenza virus infection in mice. Nature. 1962, 195: 1129-1130. 10.1038/1951129a0.
Lowen AC, Mubareka S, Steel J, Palese P: Influenza virus transmission is dependent on relative humidity and temperature. PLoS Pathog. 2007, 3: 1470-1476.
Shaman J, Kohn M: Absolute humidity modulates influenza survival, transmission, and seasonality. Proc Natl Acad Sci USA. 2009, 106: 3243-3248. 10.1073/pnas.0806852106.
Shaman J, Pitzer VE, Viboud C, Grenfell BT, Lipsitch M: Absolute Humidity and the Seasonal Onset of Influenza in the Continental United States. PLoS Biol. 2010, 8: e1000316-10.1371/journal.pbio.1000316.
Warren-Gash C, Bhaskaran K, Hayward A, Leung GM, Lo SV, Wong CM, et al: Circulating influenza virus, climatic factors, and acute myocardial infarction: a time series study in England and wales and Hong Kong. J Infect Dis. 2011, 203: 1710-1718. 10.1093/infdis/jir171.
Braga AL, Zanobetti A, Schwartz J: The effect of weather on respiratory and cardiovascular deaths in 12 U.S. cities. Environ Health Perspect. 2002, 110: 859-863.
Hajat S, Kovats RS, Lachowycz K: Heat-related and cold-related deaths in England and Wales: who is at risk?. Occup Environ Med. 2007, 64: 93-100.
Baccini M, Biggeri A, Accetta G, Kosatsky T, Katsouyanni K, Analitis A, et al: Heat effects on mortality in 15 European cities. Epidemiology. 2008, 19: 711-719. 10.1097/EDE.0b013e318176bfcd.
Yang L, Wong CM, Chan KP, Chau PY, Ou CQ, Chan KH, et al: Seasonal effects of influenza on mortality in a subtropical city. BMC Infect Dis. 2009, 9: 133-10.1186/1471-2334-9-133.
Center for Health Protection: Sentinel surveillance. [http://www.chp.gov.hk]
Basu R, Feng WY, Ostro BD: Characterizing temperature and mortality in nine California counties. Epidemiology. 2008, 19: 138-145. 10.1097/EDE.0b013e31815c1da7.
Yang L, Chiu SS, Chan KP, Chan KH, Wong HS, Peiris JSM, et al: Validation of statistical models for estimating hospitalization associated with influenza and other respiratory viruses. PLoS One. 2011, 6: e17882-10.1371/journal.pone.0017882.
Hastie TJ, Tibshirani RJ: Generalized Additive Models. 1990, London: Chapman and Hall
Yang L, Ma S, Chen PY, He JF, Chan KP, Chow A, et al: Influenza associated mortality in the subtropics and tropics: results from three Asian cities. Vaccine. 2011
Ren C, Williams GM, Mengersen K, Morawska L, Tong S: Does temperature modify short-term effects of ozone on total mortality in 60 large eastern US communities?--an assessment using the NMMAPS data. Environ Int. 2008, 34: 451-458. 10.1016/j.envint.2007.10.001.
R Development Core Team. R: A Language and Environment for Statistical Computing. [2.2.1]. 2005, Vienna: R Foundation for Statistical Computing
Bouchama A, Knochel JP: Heat stroke. N Engl J Med. 2002, 346: 1978-1988. 10.1056/NEJMra011089.
Madjid M, Aboshady I, Awan I, Litovsky S, Casscells SW: Influenza and cardiovascular disease: is there a causal relationship?. Tex Heart Inst J. 2004, 31: 4-13.
Curriero FC, Heiner KS, Samet JM, Zeger SL, Strug L, Patz JA: Temperature and mortality in 11 cities of the eastern United States. Am J Epidemiol. 2002, 155: 80-87. 10.1093/aje/155.1.80.
Medina-Ramon M, Zanobetti A, Schwartz J: The effect of ozone and PM10 on hospital admissions for pneumonia and chronic obstructive pulmonary disease: a national multicity study. Am J Epidemiol. 2006, 163: 579-588. 10.1093/aje/kwj078.
Sagripanti JL, Lytle CD: Inactivation of influenza virus by solar radiation. Photochem Photobiol. 2007, 83: 1278-1282. 10.1111/j.1751-1097.2007.00177.x.
Moura FE, Perdigao AC, Siqueira MM: Seasonality of influenza in the tropics: a distinct pattern in northeastern Brazil. Am J Trop Med Hyg. 2009, 81: 180-183.
Tamerius J, Nelson M, Zhou S, Viboud C, Miller M, Alonso W: Global influenza seasonality: reconciling patterns across temperate and tropical regions. Environ Health Perspect. 2011, 119: 439-445.
The Macroepidemiology of Influenza Vaccination Study Group: The macro-epidemiology of influenza vaccination in 56 countries, 1997-2003. Vaccine. 2005, 23: 5133-5143.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/11/342/prepub
We thank Peter Chapman for making the map and proofreading the manuscript. This study was supported by the Research Fund for the Control of Infectious Diseases/Health and Health Services Research Fund ; and the Area of Excellence Scheme of the University Grants Committee of the Hong Kong Special Administrative Region [AoE/M-12/06].
JSMP served as ad hoc consultant to pharmaceutical firms Crucell MV and Sanofi Pasteur. Other authors declare that they have no competing interests.
CMW and JSMP initiated the study; PYC, JFH, CQO and APD collected and cleaned the data; LY and KPC conducted the data analysis; LY and CMW drafted the manuscript. All authors read and approved the final manuscript.