- Research article
- Open Access
- Open Peer Review
Let the sun shine in: effects of ultraviolet radiation on invasive pneumococcal disease risk in Philadelphia, Pennsylvania
BMC Infectious Diseasesvolume 9, Article number: 196 (2009)
Streptococcus pneumoniae is a common cause of community acquired pneumonia and bacteremia. Excess wintertime mortality related to pneumonia has been noted for over a century, but the seasonality of invasive pneumococcal disease (IPD) has been described relatively recently and is poorly understood. Improved understanding of environmental influence on disease seasonality has taken on new urgency due to global climate change.
We evaluated 602 cases of IPD reported in Philadelphia County, Pennsylvania, from 2002 to 2007. Poisson regression models incorporating seasonal smoothers were used to identify associations between weekly weather patterns and case counts. Associations between acute (day-to-day) environmental fluctuations and IPD occurrence were evaluated using a case-crossover approach. Effect modification across age and sex strata was explored, and meta-regression models were created using stratum-specific estimates for effect.
IPD incidence was greatest in the wintertime, and spectral decomposition revealed a peak at 51.0 weeks, consistent with annual periodicity. After adjustment for seasonality, yearly increases in reporting, and temperature, weekly incidence was found to be associated with clear-sky UV index (IRR per unit increase in index: 0.70 [95% CI 0.54-0.91]). The effect of UV index was highest among young strata and decreased with age. At shorter time scales, only an association with increases in ambient sulphur oxides was linked to disease risk (OR for highest tertile of exposure 0.75, 95% CI 0.60 to 0.93).
We confirmed the wintertime predominance of IPD in a major urban center. The major predictor of IPD in Philadelphia is extended periods of low UV radiation, which may explain observed wintertime seasonality. The mechanism of action of diminished light exposure on disease occurrence may be due to direct effects on pathogen survival or host immune function via altered 1,25-(OH)2-vitamin-D metabolism. These findings may suggest less diminution in future IPD risk with climate change than would be expected if wintertime seasonality was driven by temperature.
Many infectious diseases of public health importance exhibit predictable periodicity, with major increases in incidence during a specific season of the year . Empirical evidence of such "seasonality" has been noted by physicians for centuries, and has been prominent enough to become a part of our vernacular (e.g., "cold and flu season") [1, 2]. Despite this wealth of experiential evidence, the mechanisms underlying seasonality are poorly understood, especially diseases characterized by person-to-person transmission .
Invasive bacterial disease due to Streptococcus pneumoniae and other respiratory pathogens exhibits striking seasonality in its occurrence [3–7]. Pneumococcal infections are a common cause of severe, community-acquired illnesses, including community-acquired pneumonia requiring hospitalization, bacteremia, and meningitis . While the introduction of antibiotics dramatically reduced the case fatality rate (CFR) for pneumococcal disease, the current CFR for bacteremic pneumococcal disease is still estimated at 5-10% in the United States, and may be twice as high among the elderly and in cases of meningitis [9, 10]. The emergence of antimicrobial resistance to beta-lactam agents, macrolides, and other antibiotic classes is an important clinical concern . Although the introduction of conjugate pneumococcal vaccines has been associated with a reduction in disease incidence , the recent increase in invasive infection by non-vaccine serotypes , which may be highly resistant to commonly-used antimicrobials , suggests that this microorganism will persist in challenging both the medical and public health communities.
The incidence of IPD peaks in the winter months, with annual periodicity [4, 5], but the forces that drive this characteristic seasonality are unknown. Wintertime seasonality of communicable respiratory diseases are often assumed to be driven by seasonal changes in environmental conditions (e.g., diminished ultraviolet radiation (UV) exposure, decreased temperature) but such associations may be confounded by other seasonally varying factors, including population behaviours (e.g., clustering indoors), co-occurrence of other infections (e.g., influenza) , and frequency of laboratory testing .
A more thorough understanding of the effect of environmental factors on seasonal IPD incidence could offer significant insight into pathogenesis, improve disease forecasting, and help determine the likely direction of pneumococcal disease incidence in the face of global climate change . Our objective was to investigate how environmental factors influence IPD occurrence in Philadelphia County. We used both traditional analytic methods (i.e., Poisson regression with seasonal smoothers) and a novel case-crossover method to examine the effects of acute weather fluctuations on IPD occurrence. Both methods reduce confounding by environmental, behavioural, and infectious influences that might otherwise distort the observed magnitude of environmental effects on disease risk.
Philadelphia County encompasses an area of 369 km2 in south-eastern Pennsylvania, and is coterminous with the City of Philadelphia (population 1,517,550 in the year 2000 ). The population receives public-health services from the Philadelphia Department of Public Health (PDPH). IPD has been a notifiable condition in the Commonwealth of Pennsylvania since 2002; Pennsylvania uses the uniform case definition endorsed by the National Notifiable Diseases Surveillance System . A case is considered "confirmed" when a consistent clinical syndrome occurs in association with the isolation of S. pneumoniae from a normally sterile site (e.g. blood, cerebrospinal or pleural fluid). Data on IPD case occurrence in Philadelphia was obtained from PDPH records, and included date of onset, age, sex, race and ethnicity of the patient, isolation site, and fatal outcome (if known).
Environmental Exposure Data
Meteorological data including temperature, relative humidity, wind speed, atmospheric pressure, and precipitation for the period from 2002 to 2006 was obtained from the weather station at Philadelphia International Airport, located eight kilometres southwest of Philadelphia's city center . Information pertaining to air quality in Philadelphia County during the years of interest-including concentrations of lead, ozone, sulphur oxides and particulate matter-was obtained from the Environmental Protection Agency . Because daily readings were taken at various locations throughout the region, the arithmetic means of the air quality values were used as exposure variables. UV index forecast estimates for Philadelphia during the same period were retrieved from the National Weather Service Climate Prediction Center . Clear-sky UV indices represent an integral of measured UV radiation levels weighted by the ability of the different UV wavelengths to cause skin erythema. The issued UV index is a similar measure, which accounts for the effect of clouds on radiation transmission; because of inconsistencies in cloud measurement during the study period, we used the clear-sky UV index as our exposure variable.
Rates of invasive pneumococcal disease were calculated using demographic data for Philadelphia County from the year 2000 US Census, as well as 2006 population estimates from the Bureau of the Census, with linear interpolation and extrapolation used to generate estimates for population by age and sex in other years . We evaluated the seasonality of disease occurrence through construction of periodograms and autocorrelograms [15, 22] for weekly case counts. As yearly periodicity was observed, we estimated seasonal and year-on-year trends in IPD occurrence using Poisson regression models that incorporated sine and cosine oscillators, with 52 week (annual) frequencies (i.e., incorporated fast Fourier transforms) [7, 22].
Using these parameters, the expression for the expected number of case counts for a given week, E (Y) is given by:
where cases is an autoregressive model term reflecting the cumulative case count in the month prior to case occurrence, (i.e., cases = ).
The phase-shift of the composite waveform generated by combining sine and cosine components of the above equation can be approximated as the arctangent of β 2/β 3, and can be used to estimate the timing of peak disease incidence . We also included model terms that controlled for longer term trends in incidence of invasive pneumococcal disease, which may have reflected the initiation of surveillance, the introduction of public funding for conjugate pneumococcal vaccination , changes in medical diagnostic practices, or other long-term changes in real or apparent pneumococcal epidemiology. As year-on-year trends in disease occurrence reflected a non-linear increase in disease risk, we used models that incorporated both linear and quadratic yearly terms.
The quadratic model term was statistically significant, but is difficult to interpret, thus we present our final Poisson model with separate linear yearly terms for the period from 2002 to 2003, and the period from 2004 to 2007 . We evaluated the impact of environmental exposures on weekly IPD incidence by incorporating exposure variables into Poisson models both individually, and using a backwards elimination algorithm (with variables retained for P < 0.2 ).
To explore the possibility of effect modification by subject characteristics, we evaluated stratum-specific estimates of effect for age categories and genders. Heterogeneity of effects across strata was assessed using meta-analytic techniques, including both graphical inspection and calculation of meta-analytic Q-statistics . We further explored the sources of between-stratum heterogeneity through construction of meta-regression models that estimate the contribution of group-level covariates to between-stratum variation in effects .
We used a case-crossover approach to evaluate acute (i.e., daily) associations between environmental exposures and IPD occurrence. This approach provides a means for evaluating the association between brief, transient exposures and rare outcomes. The design is characterized by "self matching", in that cases serve as their own controls. In the context of environmental epidemiology, a "case" is a day on which a case occurred, while a "control" is an appropriately selected day on which a case did not occur . We used a time-stratified 2:1 matched case-crossover design in which hazard periods were defined as the reported date of IPD onset from Philadelphia County public health. Beginning on January 1, 2002 the person-time at risk was divided into three-week time strata. Control periods were chosen by matching the hazard period by day of the week within the stratum, and could both precede, both follow, or straddle the hazard period [26, 27]. Random directionality of control selection was used in order to avoid biases that can occur with unidirectional or uniform bidirectional control selection . The 1-3 day incubation period of S. pneumoniae was used to estimate the lag days between acute environmental occurrence and case onset, or plausible effect period . We also evaluated effects during the period immediately preceding incubation (i.e., 4-6 day lags) to evaluate the possibility that environmental conditions might affect risk via enhanced transmission of S. pneumoniae. We evaluated the effects of both raw environmental exposures, and quantile ranks within time strata through construction of conditional logistic regression models . Analyses were performed using SAS version 8.01 (SAS Institute, Cary, NC) and Stata version 9.1 (Stata Corporation, College Station, TX).
Between January 1, 2002 and April 30, 2007 there were 602 reported cases of IPD in Philadelphia County, for a crude annual incidence of 6.84 cases per 100,000 [95% CI 6.32 to 7.42]. Incidence of IPD was highest at the age extremes, with a slight male:female predominance (incidence rate ratio (IRR) 1.16, 95% CI 0.99 to 1.36); the case fatality rate was approximately 10% (Table 1).
Both spectral decomposition and construction of autocorrelograms identified annual periodicity of infection, with peak incidence in mid-February (phase = 6 weeks) (Figures 1 and 2). Strong statistical evidence for seasonal oscillation was obtained from Poisson regression models (P for seasonal oscillation < 0.001). A significant annual increase in incidence was seen throughout the study period, though this was more marked prior to 2004 (IRR per year 1.34, 95% CI 1.08 to 1.66) than subsequently (IRR 1.22, per year 95% 1.16 to 1.34) (Figure 3). We found no clear trends in the incidence of IPD or case-fatality in individual age groups, and no significant heterogeneity was detected between age groups with respect to year-on-year trends in incidence or case-fatality (see Additional File 1 and Additional File 2).
Weekly Environmental Effects
In univariable models the risk of IPD increased with several seasonally oscillating environmental exposures, including temperature, humidity, pressure, air pollution, and UV radiation, as shown in Table 2. Risk of IPD increased with average weekly barometric pressure, sulphur and nitrous oxides, and decreased with average weekly temperature, relative humidity, and UV index. However, after controlling for seasonal oscillation and longer term temporal trends, only cooling-degree days (i.e., average number of degrees above 18°C), maximum temperature, and clear-sky UV index were independently associated with case occurrencein a final multivariable model (Table 2). Risk increased with increasing maximal temperature, but decreased with cooling degree days, indicating a threshold of 18°C above which the qualitative effect of temperature changed.
There was significant heterogeneity in the effect of clear-sky UV index in individuals aged less than 15 years (IRR 0.26, 95% CI 0.11 to 0.61) as compared to older individuals (IRR 0.81, 95% CI 0.64 to 1.03), (Q-statistic 6.31 on 1 d.f., P = 0.012). No heterogeneity in UV effects was found when evaluated in individuals aged over 64 as compared to younger individuals (Q-statistic 0.48 on 1 d.f., P = 0.523), or when evaluated in males as compared to females (Q-statistic 1.68 on 1 d.f., P for heterogeneity = 0.194). Meta-regression models constructed using stratum-specific estimates of effect for age groups, by 10-year age increments, identified a log-linear trend in the effect of UV radiation on disease risk by age stratum, with strongest effects in youngest strata, and effects diminishing with increasing age (change in natural logarithm of IRR = 0.13 per decade, 95% CI 0.04 to 0.23, P = 0.007) (Figure 4).
Acute Environmental Effects
Evaluating associations between environmental and meteorological exposures and IPD risk using a case-crossover approach, we identified an inverse association between ambient levels of sulphur oxides and disease risk during the likely incubation period (Table 3). No other significant associations between environmental exposures and risk were identified either during the incubation period, or in the period immediately prior to incubation. In particular, occurrence of IPD was not affected by daily changes in clear-sky UV index, temperature or cooling-degree days, in contrast to associations on longer time scales described above.
Notwithstanding the existence of vaccination and effective antibiotic therapy, invasive pneumococcal disease remains an important source of population morbidity and mortality. The seasonality of IPD is well recognized, but poorly understood. Epidemiological mechanisms invoked to explain this pattern have included co-occurrence of other infectious diseases , wintertime social gatherings , and seasonal oscillation in immune function . However, concurrent seasonal changes in a variety of environmental, behavioural, and epidemiogical exposures make identification of causal associations particularly challenging . We attempted to address this challenge by using analytic approaches that should control for seasonal confounders, known and unknown, at different time scales. At a weekly time scale we found increases in UV radiation to be most strongly associated with decreased numbers of invasive pneumococcal disease cases, though average temperatures also appeared to influence disease risk. At short time scales, fluctuations in ambient air quality, as manifested by differences in concentrations of sulphur oxides, were associated with changes in risk. The direction of this association was at variance with existing models relating air pollution to pneumonia occurrence [3, 31, 32].
The casual (as opposed to causal) association between low UV radiation in the winter and surges in respiratory disease has been noted previously , and has been proposed as an important driver of influenza seasonality, but has not to our knowledge been evaluated in a way that accounts for coincident seasonal changes in other seasonally oscillatory factors. The degree to which such seasonal oscillation can result in "just so" stories that lead to misattribution of causation to non-specific seasonal exposures is highlighted in the univariable analyses we conducted without including seasonal oscillators. In these models, a variety of environmental conditions, including weather variables and air quality indices, were strongly associated with IPD risk. However, after controlling for non-specific seasonality, only UV radiation (and, more weakly, temperature) were associated with disease risk; indeed, the apparent protective effect of UV radiation was actually strengthened after controlling for seasonal oscillation. The interpretation of such a model is that increases in UV radiation reduce IPD risk, even after accounting for the fact that IPD risk is maximal during low-UV periods of the year.
An important consideration is whether changes in UV radiation, operating at a weekly time scale, constitute a biologically plausible mechanism that explains seasonal oscillation in pneumococcal disease risk. Indeed, there are several mechanisms that may have substantial biological plausibility. Modulation of risk may occur through direct effects of UV radiation on host immune function: Dowell reviewed a variety of immunological changes associated with diminished UV radiation exposure in experimental settings, and noted that in granulocyte and monocyte function were reduced during periods of short light exposure . UV radiation also influences the production of 1,25-(OH)2-vitamin D, which has important immunomodulatory functions . Namely, enhanced maturation of macrophages, macrophage secretion of bactericidal substances such as lysozomal enzyme phosphatase and hydrogen peroxide [33, 35], and secretion of antimicrobial peptides (including cathelicidins and defensins) by both immune cells, and respiratory tract epithelium . Vitamin D deficiency is associated with a marked increase in the risk of pneumonia , and most human vitamin D is acquired via sun exposure. Thus, extended periods of low UV light could result in an increased susceptibility to S. pneumoniae infection resulting from a lack of vitamin D production, though the week-to-week fluctuations in risk described in this paper may be too rapid to represent a vitamin D effect.
An alternative mechanism of action of UV radiation in reducing IPD risk could be direct effects of radiation on pathogen survival. Many bacterial respiratory pathogens, including pneumococcus, are transmitted in respiratory secretions over short distances (i.e., via "large droplet transmission"), and thus encounter the physical environment directly during transmission events. The bactericidal effects of UV-B radiation have been well known for decades; such radiation inactivates bacteria by causing harmful genetic mutations through creation of pyrimidine dimers . A model of UV effect via diminished transmissibility, rather than decreased host susceptibility, is supported by our finding that UV effect is strongest in the youngest individuals in the population (i.e., toddlers), where disease risk is likely to be driven by mobility, contact with peers, prolonged carriage and carelessness with respiratory secretions. We found very little protective effect against IPD in oldest individuals, whose risk may be more strongly linked to immune senescence than to high rates of contact with infectious contemporaries [39, 40].
An unexpected finding in our case-crossover analysis was the identification of increased levels of ambient sulphur oxides with diminished risk of IPD. This may represent a chance association, due to the established association between sulphur dioxide and adverse respiratory outcomes , as well as the previously described correlation between ambient sulphur dioxide and pneumonia risk . Nonetheless, this association may warrant further exploration; for example, ambient air pollutants might have adverse effects on respiratory tract pathogens as well as hosts.
Other findings in this study are also worthy of comment; in Philadelphia, the wintertime peak in pneumococcal incidence occurred about six weeks later than previously been described in U.S. adults . In addition, a gradual increase in cases was observed during the 5-year study period. We suspect that this increase is less likely to represent a true surge in disease rates, which have actually been falling in the U.S. with the introduction of 7-valent conjugate pneumococcal vaccine ; rather, we suspect that this increase, which is more attenuated after 2003, represents the gradual increase in reporting that commonly follows implementation of new infectious disease reporting requirements .
Like any observational study, ours has several limitations. First among these is our heavy reliance on public health surveillance data which is known to suffer from under-reporting of notifiable infectious diseases . Thus our data set may be incomplete, consisting of only a subset of the cases of IPD in Philadelphia during the study period. However, selection bias will only be introduced if the notification of infectious diseases to public health was correlated with meteorological patterns (e.g., cases which occur during days with higher UV index have a greater likelihood of being diagnosed and reported), which seems unlikely . Second, we obtained weather data from a single site at Philadelphia International Airport, and air pollution and UV data were averaged over several locations throughout the county, which may not represent the true exposure status of individual cases. This is probable non-differential misclassification, and will bias the results towards the null hypothesis. The effects seen here are most likely conservative, weakening the strength of observed associations and suggesting our estimates tend towards the lower bound .
In summary, we described the occurrence of IPD in a major U.S. urban center and found that incidence was associated with marked wintertime seasonality that may be partly explained by diminished exposure to UV-B radiation in winter months. Further study is needed, but this result is consistent with observed patterns of respiratory infectious disease, and would be consistent with several biologically plausible models of effect. As the nature of future changes in UV radiation related to climate change are more difficult to predict than general changes in temperatures or precipitation patterns, the implications of these findings for future pneumococcal disease epidemiology are unclear . Nonetheless, we believe this observation goes some way towards explaining the notable seasonality of IPD, and could conceivably lead to novel disease control strategies through improved understanding of this common and virulent infectious disease.
Dowell SF, Ho MS: Seasonality of infectious diseases and severe acute respiratory syndrome-what we don't know can hurt us. The Lancet Infectious Diseases. 2004, 4: 704-8. 10.1016/S1473-3099(04)01177-6.
Singer C: A Short History of Medicine. Ancient Greece. 1928, New York: Oxford University Press, 1-26. 1
Kim PE, Musher DM, Glezen WP, Rodriguez-Barradas MC, Nahm WK, Wright CE: Association of invasive pneumococcal disease with season, atmospheric conditions, air pollution, and the isolation of respiratory viruses. Clin Infect Dis. 1996, 22: 100-6.
Talbot TR, Poehling KA, Hartert TV, Arbogast PG, Halasa NB, Edwards KM, Schaffner W, Craig AS, Griffin MR: Seasonality of invasive pneumococcal disease: temporal relation to documented influenza and respiratory syncytial viral circulation. American Journal of Medicine. 2005, 118: 285-91. 10.1016/j.amjmed.2004.09.016.
Dowell SF, Whitney CG, Wright C, Rose CE, Schuchat A: Seasonal patterns of invasive pneumococcal disease. Emerging Infectious Diseases. 2003, 9: 573-9.
Fisman DN, Lim S, Wellenius GA, Johnson C, Britz P, Gaskins M, Maher J, Mittleman MA, Spain CV, Haas CN, et al: It's not the heat, it's the humidity: wet weather increases legionellosis risk in the greater Philadelphia metropolitan area. Journal of Infectious Diseases. 2005, 192: 2066-73. 10.1086/498248.
Jensen ES, Lundbye-Christensen S, Pedersen L, Sorensen HT, Schonheyder HC: Seasonal variation in meningococcal disease in Denmark: relation to age and meningococcal phenotype. Scandinavian Journal of Infectious Diseases. 2003, 35: 226-9. 10.1080/00365540310010921.
Austrian R: The pneumococcus at the millennium: not down, not out. Journal of Infectious Diseases. 1999, 179 (Suppl 2): S338-41. 10.1086/513841.
Robinson KA, Baughman W, Rothrock G, Barrett NL, Pass M, Lexau C, Damaske B, Stefonek K, Barnes B, Patterson J, et al: Epidemiology of invasive Streptococcus pneumoniae infections in the United States, 1995-1998: Opportunities for prevention in the conjugate vaccine era. JAMA. 2001, 285: 1729-35. 10.1001/jama.285.13.1729.
Ortqvist A, Hedlund J, Kalin M: Streptococcus pneumoniae: epidemiology, risk factors, and clinical features. Seminars in Respiratory & Critical Care Medicine. 2005, 26: 563-74.
Daneman N, McGeer A, Green K, Low DE: Macrolide resistance in bacteremic pneumococcal disease: implications for patient management. Clin Infect Dis. 2006, 43: 432-8. 10.1086/505871.
Whitney CG, Farley MM, Hadler J, Harrison LH, Bennett NM, Lynfield R, Reingold A, Cieslak PR, Pilishvili T, Jackson D, et al: Decline in invasive pneumococcal disease after the introduction of protein-polysaccharide conjugate vaccine. N Engl J Med. 2003, 348: 1737-46. 10.1056/NEJMoa022823.
Emergence of antimicrobial-resistant serotype 19A Streptococcus pneumoniae--Massachusetts, 2001-2006. MMWR Morb Mortal Wkly Rep. 2007, 56: 1077-80.
Pichichero ME, Casey JR: Emergence of a multiresistant serotype 19A pneumococcal strain not included in the 7-valent conjugate vaccine as an otopathogen in children. Jama. 2007, 298: 1772-8. 10.1001/jama.298.15.1772.
Fisman DN: Seasonality of infectious diseases. Annu Rev Public Health. 2007, 28: 127-43. 10.1146/annurev.publhealth.28.021406.144128.
Intergovernmental Panel on Climate Change: Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. 2007, Cambridge, UK: Cambridge University Press
Population Division, US Census Bureau: Population estimates. 2007, City: U.S. Department of Commerce, Last accessed January 29, 2008., [http://www.census.gov/popest/estimates.php]
Case definitions for infectious conditions under public health surveillance. Centers for Disease Control and Prevention. MMWR Recomm Rep. 1997, 46: 1-55.
National Climate Data Center, NOAA Satellite and Information Service: Online climate data directory. 2007, City: National Oceanic and Atmospheric Administration, U.S. Department of Commerce, Last accessed January 29, 2008., [http://lwf.ncdc.noaa.gov/oa/climate/climatedata.html#surface]
U.S. Environmental Protection Agency: Air Quality System. 2007, City: U.S. Environmental Protection Agency, Last accessed January 29, 2008., [http://www.epa.gov/ttn/airs/airsaqs/aqsweb/]
Climate Prediction Center, NW Service: UV Index Bulletins Archive. 2007, City: National Oceanic and Atmospheric Administration, U.S. Department of Commerce, Last accessed January 29, 2008., [http://www.cpc.ncep.noaa.gov/products/stratosphere/uv_index/uv_archive.shtml]
Diggle P: Time Series: A Biostatistical Introduction. 1990, Oxford, UK: Oxford University Press
Woodward M: Modelling quantitative outcome variables. Epidemiology: Study Design and Data Analysis. 2005, Boca Raton, FL: Chapman & Hall/CRC, 2
Hosman D, Lemeshow S: Applied Logistic Regression. 1989, New York: John Wiley and Sons, 82-134. 1
Deeks J, Altman D: Effect measures for meta-analysis of trials with binary outcomes. Systematic Reviews in Health Care: Meta-Analysis in Context. Edited by: Egger M, Davey Smith G, Altman D. 1995, London, UK: BMJ Publishing Group, 313-335.
Levy D, Lumley T, Sheppard L, Kaufman J, Checkoway H: Referent selection in case-crossover analyses of acute health effects of air pollution. Epidemiology. 2001, 12: 186-92. 10.1097/00001648-200103000-00010.
Lumley T, Levy D: Bias in the case crossover design: implications for studies of air pollution. Environmetrics. 2000, 689-704. 10.1002/1099-095X(200011/12)11:6<689::AID-ENV439>3.0.CO;2-N.
Control of Communicable Diseases Manual. 2000, Washington, DC: American Public Health Association, 17
Dowell SF: Seasonal variation in host susceptibility and cycles of certain infectious diseases. Emerging Infectious Diseases. 2001, 7: 369-74.
Dushoff J, Plotkin JB, Levin SA, Earn DJ: Dynamical resonance can account for seasonality of influenza epidemics. Proc Natl Acad Sci USA. 2004, 101: 16915-6. 10.1073/pnas.0407293101.
Schwartz J: PM10, ozone, and hospital admissions for the elderly in Minneapolis-St. Paul, Minnesota. Arch Environ Health. 1994, 49: 366-74. 10.1080/00039896.1994.9954989.
Health effects of outdoor air pollution. Committee of the Environmental and Occupational Health Assembly of the American Thoracic Society. Am J Respir Crit Care Med. 1996, 153: 3-50.
Cannell JJ, Vieth R, Umhau JC, Holick MF, Grant WB, Madronich S, Garland CF, Giovannucci E: Epidemic influenza and vitamin D. Epidemiology & Infection. 2006, 134: 1129-40.
Dowell SF: Seasonal variation in host susceptibility and cycles of certain infectious diseases. Emerg Infect Dis. 2001, 7: 369-74.
van Etten E, Mathieu C: Immunoregulation by 1,25-dihydroxyvitamin D3: basic concepts. Journal of Steroid Biochemistry & Molecular Biology. 2005, 97: 93-101.
Schutte BC, McCray PB: [beta]-defensins in lung host defense. Annual Review of Physiology. 2002, 64: 709-48. 10.1146/annurev.physiol.64.081501.134340.
Muhe L, Lulseged S, Mason KE, Simoes EA: Case-control study of the role of nutritional rickets in the risk of developing pneumonia in Ethiopian children. Lancet. 1997, 349: 1801-4. 10.1016/S0140-6736(96)12098-5.
Brickner PW, Vincent RL, First M, Nardell E, Murray M, Kaufman W: The application of ultraviolet germicidal irradiation to control transmission of airborne disease: bioterrorism countermeasure. Public Health Rep. 2003, 118: 99-114. 10.1016/S0033-3549(04)50225-X.
Hogberg L, Geli P, Ringberg H, Melander E, Lipsitch M, Ekdahl K: Age- and serogroup-related differences in observed durations of nasopharyngeal carriage of penicillin-resistant pneumococci. J Clin Microbiol. 2007, 45: 948-52. 10.1128/JCM.01913-06.
Smithson SL, Kolibab K, Shriner AK, Srivastava N, Khuder S, Westerink MA: Immune response to pneumococcal polysaccharides 4 and 14 in elderly and young adults: analysis of the variable light chain repertoire. Infect Immun. 2005, 73: 7477-84. 10.1128/IAI.73.11.7477-7484.2005.
Lexau CA, Lynfield R, Danila R, Pilishvili T, Facklam R, Farley MM, Harrison LH, Schaffner W, Reingold A, Bennett NM, et al: Changing epidemiology of invasive pneumococcal disease among older adults in the era of pediatric pneumococcal conjugate vaccine. Jama. 2005, 294: 2043-51. 10.1001/jama.294.16.2043.
Janes GR, Hutwagner LC, Cates W, Stroup DF, Williamson GD: Descriptive Epidemiology: Analyzing and Interpreting Surveillance Data. Principles and Practice of Public Health Surveillance. Edited by: Steven REC, Teutsch M. 2000, New York: Oxford University Press, 116-168. 2
Doyle TJ, Glynn MK, Groseclose SL: Completeness of notifiable infectious disease reporting in the United States: an analytical literature review. American Journal of Epidemiology. 2002, 155: 866-74. 10.1093/aje/155.9.866.
Grimes DA, Schulz KF: Bias and causal associations in observational research. Lancet. 2002, 359: 248-52. 10.1016/S0140-6736(02)07451-2.
Talaas P, Kaurola J, Kylling A, Shindell D, Sausen R, Dameris M, Grewe V, Herman J, Damski J, Steil B: The impact of greenhouse gases and halogenated species on future solar UV radiation doses. Geophysical Research Letters. 2000, 27: 1127-30. 10.1029/1999GL010886.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/9/196/prepub
The authors declare that they have no competing interests.
AW performed the majority of the analysis and interpretation of the data, and drafted the manuscript. CJ and CS acquired the case and exposure data, and each performed revisions of the manuscript. LK analysed and interpreted data, and was involved in drafting and revising the manuscript. VN adapted the analysis tools for the data, assisted in its interpretation and was involved in critically reviewing the manuscript. DF performed data analysis, interpretation, and was involved in both the drafting and critical review process of the manuscript. All authors have given approval for the final version to be published.