- Open Access
Seasonality of respiratory viruses and bacterial pathogens
Antimicrobial Resistance & Infection Controlvolume 8, Article number: 125 (2019)
Seasonal variation has been observed for various bacterial and viral infections. We aimed to further study seasonality of respiratory viruses and bacterial pathogens in relation to antibiotic use, as well as meteorological parameters.
An ecologic study of antibiotic exposure, meteorological parameters, detection of respiratory viruses and clinical isolates of Clostridioides difficile, Methicillin-resistant Staphylococcus aureus (MRSA), Streptococcus pneumoniae, and Escherichia coli and Klebsiella pneumoniae (grouped together as gram-negative bacteria; GNB) in Rhode Island from 2012 to 2016.
Peak detection of C. difficile occurred 3 months after the peak in antibiotic prescriptions filled (OR = 1.24, 95% CI, 1.07–1.43; P = 0.006). Peak MRSA detection was noted 7 months after the peak in antibiotic prescriptions filled (OR = 1.69, 95% CI, 1.21–2.35; P = 0.003) and 10 months after the peak in respiratory virus detection (OR = 1.04, 95% CI, 1.01–1.06; P = 0.003). Peak GNB detection was noted 2 months after the peak mean monthly ambient temperature (OR = 1.69, 95% C.I., 1.20–2.39; P = 0.004). Peak detection of S. pneumoniae was noted at the same time as the peak in detection of respiratory viruses (OR = 1.01, 95% C.I., 1.00–1.01; P = 0.015).
We identified distinct seasonal variation in detection of respiratory viruses and bacterial pathogens. C. difficile seasonality may, in part, be related to antibiotic prescriptions filled; GNB seasonality may be related to ambient temperature and S. pneumoniae may be related to concurrent respiratory viral infections.
Many host-related risk factors have been established for bacterial infections; however, non-intrinsic factors such as respiratory viral infections, weather conditions, and latitude are among other factors that lead to seasonality of such infections [1,2,3,4,5,6]. In addition, Clostridioides difficile infection is temporally associated with respiratory tract infections [6,7,8], likely due to antibiotic overprescribing for respiratory viral infections.
A better understanding of seasonality of various pathogens, as well as an association with antibiotic prescriptions may allow us to devise and implement future public health interventions (e.g., vaccination and antimicrobial stewardship), particularly in the community setting. The purpose of this study was to assess seasonality of respiratory viruses and bacterial pathogens in relation to antibiotic use and meteorological parameters in one locale.
This was an ecological study involving detection of bacterial pathogens and respiratory viruses in Rhode Island. Microbiologic data included specimens sent to the Lifespan Microbiology Laboratory from inpatients and outpatients assessed in the Lifespan healthcare system (including Rhode Island Hospital/Hasbro Children’s Hospital, Miriam Hospital, Newport Hospital, Bradley Hospital, and Gateway Healthcare). The system is licensed for 1,165 beds, with 62,100 patient discharges, 257,081 emergency department visits, and 757,380 outpatient visits , as well as patient specimens from other outpatient facilities throughout Rhode Island between January 1, 2012 and December 31, 2016. We also obtained data regarding antibiotics prescriptions, total precipitation and average temperature.
Data source for detection of respiratory viruses and bacterial pathogens
Rhode Island Hospital infection control software system (TheraDoc; Premier, Charlotte, NC) was queried to find positive microbiology laboratory identification of respiratory viruses identified by the respiratory viral panel (Luminex, Austin, TX), rapid influenza test (Xpert; Cepheid, Sunnyvale, CA) or rapid Respiratory Syncytial Virus test (Xpert; Cepheid, Sunnyvale, CA) from an outpatient’s or inpatient’s nasopharyngeal swab or bronchoscopic lavage specimen. The target population included all ages regardless of insurance status. The respiratory viral panel included testing for Influenza A/B, respiratory syncytial virus (RSV) A/B, Coronavirus, Parainfluenza virus, human Metapneumovirus, Adenovirus, Rhinovirus/Enterovirus (the panel does not differentiate Rhinovirus and Enterovirus). Clostridioides difficile toxin was identified by PCR testing of stool specimens (BD GeneOhm; Xpert); growth in clinical cultures identified cases involving Methicillin-resistant Staphylococcus aureus (MRSA), Streptococcus pneumoniae, and Escherichia coli and Klebsiella pneumoniae (grouped together as gram-negative bacteria; GNB). Only the first encounter with one of the above-noted pathogens was included if a patient had multiple episodes of positive laboratory testing for that pathogen within 30 days. We excluded MRSA identified only from nares screening. Incidence reflected counts with no denominator based on the assumption that there was no significant year-to-year variability between practices in obtaining respiratory specimens (Additional file 1: Figure S1).
Antibiotics and meteorological data source
Antibiotic prescription data from January 2012 through December 2016 was used to assess antibiotic exposure. Data was obtained from the Rhode Island All Payer Claim Database (ACPD), provided by Rhode Island Department of Health. The ACPD collects and stores payer enrollment data, medical claims, pharmacy claims, and provider data on a monthly basis from commercial insurers, Medicare and Medicaid regarding antibiotic prescriptions filled . The denominator for prescriptions per 1,000 population was derived from Rhode Island census data . Antibiotic classes were: cephalosporins, clindamycin, macrolides, metronidazole, nitrofurantoin, penicillins, quinolones, tetracyclines, and trimethoprim-sulfamethoxazole.
Monthly aggregated data on average temperature and total precipitation from January 2012 through December 2016 in Rhode Island was retrieved from the National Centers for Environmental Information .
The null hypothesis was no seasonal cross-correlation between respiratory viral activity, antibiotic prescriptions filled, weather variables and detection of different bacterial pathogens. We constructed a longitudinal incidence rate using a monthly dataset which included respiratory viruses detected, antibiotic prescriptions filled, and meteorological parameters (total precipitation and average temperature), as well as detected bacterial pathogens (C. difficile, MRSA, GNB, and S. pneumoniae) (Fig. 1). A seasonal trend decomposition procedure, based on Locally Weighted Scatterplot Smoothing (STL), was conducted for each bacterial pathogen to assess for seasonality and trends associated with detection of respiratory viruses, antibiotic prescription, and meteorological parameters . An additive decomposition model was used. To assess for a correlation between a time series and a given a number of lags, we measured the cross-correlation of Xt and Yt + k for each month. Separate cross-correlation functions were applied to determine the specific bacterial pathogens and their highest correlation with detection of respiratory viruses, antibiotic prescriptions filled, or meteorological parameters on defined time lags. We calculated odds ratios (OR) with 95% confidence intervals, estimating risk of elevated incidence compared to annual average incidence. Due to the small numbers and multiple testing, the p-values should only be seen as descriptive measures. Analyses were performed using R (ver. 3.4.3; R Development Core Team, Vienna, Austria). Packages used were: forecast, TSA (time series analysis) , and ASTA (applied statistical time series analysis) .
From January 2012 through December 2016, 6,857 respiratory viruses were detected and 3,065,789 antibiotic prescriptions were filled with distinct seasonality from December through February (Fig. 2-A). Peak precipitation and temperature were noted from June through August annually (Fig. 2-B). There were 16,419 bacterial pathogens identified: 2,354 C. difficile, 4,026 MRSA, 9,767 GNB, and 272 S. pneumoniae (Fig. 2-C). The level of detection of bacterial pathogens varied in month of onset, peak prevalence and duration of detection.
The peak in antibiotic prescriptions filled occurred 3 months before the peak in detection of C. difficile (OR = 1.24, 95% C.I., 1.07–1.43; P = 0.006; Table 1). Peak detection of MRSA occurred 7 months after the peak in antibiotic prescriptions filled (OR = 1.69, 95% C.I., 1.21–2.35; P = 0.003) and 10 months after the peak in respiratory virus detection (OR = 1.04, 95% C.I., 1.01–1.06; P = 0.003). The peak GNB detection occurred 2 months after the peak mean monthly temperature (OR = 1.69, 95% C.I., 1.20–2.39; P = 0.004). The peak in detection of S. pneumoniae occurred at the peak in detection of respiratory viruses (OR = 1.01, 95% C.I., 1.00–1.01; P = 0.015).
Peak detection of C. difficile occurred a few months after the peak in antibiotic prescriptions during winter months, as previously described [7, 8]. Previous reports have also correlated increased C. difficile infections 1–2 months after peak macrolide and fluoroquinolone use .
We previously reported peaks of MRSA infections during the summer and autumn seasons . This seasonality may, in part, be associated with increased antibiotic prescriptions some months earlier, as previously reported . The 7-month lag that we detected is longer than found in another study which detected a temporal relationship between fluoroquinolone prescriptions and ciprofloxacin-resistant MRSA after a 1-month lag period . A recent study revealed that antibiotics lacking MRSA activity, including fluoroquinolones, promote MRSA overgrowth in the nares and this may be related to finding increased MRSA detection in clinical isolates after exposure to such antibiotics . Not surprisingly, a national antibiotic stewardship program has been associated with decreased MRSA prevalence . We also found the circulation of respiratory viruses is likely associated with bacterial co-infection or superinfection due to S. pneumoniae (0-month lag with OR = 1.01) as previously reported .
Our findings are similar to a study of 132 US hospitals that found higher outdoor temperatures associated with an increased frequency of bloodstream infections caused by GNB . A review of the published literature also revealed that increased temperature correlated with bloodstream infections due to GNB . A study from Brazil found that average temperature on the day of diagnosis was associated with Klebsiella spp. (OR 1.19; 95% CI 1.07–1.33) and A. baumannii (OR 1.20; 95% CI 1.07–1.34) . However, this finding may not be replicated in different geographic locations with varying weather conditions . The reasons for this association are unclear but likely, in part, reflects some advantage in environmental growth conditions with warmer temperatures among various GNB. It is possible that seasonal changes in the human microbiome plays a role .
There are limitations of this study. We collected de-identified data using infection control software rather than detailed medical record data. Misclassification bias may be present because we could not exclude clinically-insignificant microbial colonizers or bystanders which may have overestimated the number of true infections. We may have had selection bias since the microbiologic data was collected within our hospital system and outpatient facilities served by our microbiology laboratory and this data may not be reflective of all Rhode Islanders. Respiratory viral infections in children could have a different seasonal pattern compared to adults. There may be detection bias since healthcare providers may not have thought of ordering tests for respiratory virus detection outside of the winter season. In addition, there may have been a change in the denominator population captured over the time course that could potentially influence these findings. Since this is an ecological study, the temporal association can be easily confounded with other factors that are seasonal such as duration of antibiotics and the level of humidity. Viruses and bacteria in this study were detected in patients in one of our three acute care hospitals which comprises 55% of non-Veteran’s Administration medical center acute care beds in Rhode Island, and at various outpatient settings in Rhode Island. Antibiotic prescription data was derived from a database that included people throughout Rhode Island, some of whom were not cared for in our healthcare system. Thus, we cannot associate antibiotic prescriptions filled and microbiologic test results on an individual level. Despite these potential limitations, we are not aware of other studies that assessed the relationships among bacterial pathogens, respiratory viruses, antibiotic exposure, and meteorological parameters in the same geographic region.
In conclusion, we attempted to study the complex ecology of bacterial and viral infections, and assess the potential impact of antibiotic exposure and climactic fluctuations. The temporal association between antibiotic prescriptions filled with detection of C. difficile highlights the potential impact of robust antibiotic stewardship on reducing C. difficile risk. Prevention of respiratory viral infections may reduce risk of S. pneumoniae infections as previously described. A better understanding of the seasonality of GNB and MRSA may help to focus future infection prevention efforts. As such, continued research is needed to better understand why these infections follow seasonal patterns.
Availability of data and materials
All Payer Claim Database
Applied statistical time series analysis
Methicillin-resistant Staphylococcus aureus
Respiratory syncytial virus
Locally Weighted Scatterplot Smoothing
Time series analysis
Rucinski SL, Murphy MP, Kies KD, Cunningham SA, Schuetz AN, Patel R. Eight years of clinical Legionella PCR testing illustrates a seasonal pattern. J Infect Dis. 2018;218:669–70.
Fisman D, Patrozou E, Carmeli Y, Perencevich E, Tuite AR, Mermel LA. Geographical variability of bacteremia study group. Geographical variability in the likelihood of bloodstream infections due to gram-negative bacteria: correlation with proximity to the equator and health care expenditure PLoS One. 2014;9:e114548.
Weinberger DM, Grant LR, Steiner CA, et al. Seasonal drivers of pneumococcal disease incidence: impact of bacterial carriage and viral activity. Clin Infect Dis. 2014;58:188–94.
Richet H. Seasonality in gram-negative and healthcare-associated infections. Clin Microbiol Infect. 2012;18:934–40.
Mermel LA, Machan JT, Parenteau S. Seasonality of MRSA infections. PLoS One. 2011;6:e17925.
Paul M. Seasonality in infectious diseases: does it exist for all pathogens? Clin Microbiol Infect. 2012;18:925–6.
Brown KA, Daneman N, Arora P, Moineddin R, Fisman DN. The co-seasonality of pneumonia and influenza with Clostridium difficile infection in the United States, 1993-2008. Am J Epidemiol. 2013;178:118–25.
Polgreen PM, Yang M, Bohnett LC, Cavanaugh JE. A time-series analysis of Clostridium difficile and its seasonal association with influenza. Infect Control Hosp Epidemiol. 2010;31:382–27.
Lifespan. 2017 Annual Report. https://www.lifespan.org/about-lifespan/lifespan-reports. Accessed 1 Sept 2018.
Rhode Island Department of Health. Rhode Island All-Payer Claims Database. https://www.riapcd-optout.com. Accessed 2 Sept 2018.
Rhode Island Department of Labor and Training. Labor Market Information. http://www.dlt.ri.gov/lmi/census.htm. Accessed 1 Sept 2018.
National Centers for Environmental Information. https://www.ncdc.noaa.gov. Accessed 1 Sept 2018.
Cleveland R, Cleveland McRae WJ, Terpenning I. STL: a seasonal-trend decomposition procedure based on loess. J Off Stat. 1990;6:3–73.
TSA: Time Series Analysis. https://cran.r-project.org/web/packages/TSA/index.html. Accessed 4 Feb 2019.
ASTSA: Applied Statistical Time Series Analysis. https://cran.r-project.org/web/packages/astsa/index.html. Accessed 4 Feb 2019.
Gilca R, Fortin E, Frenette C, Longtin Y, Gourdeau M. Seasonal variations in Clostridium difficile infections are associated with influenza and respiratory syncytial virus activity independently of antibiotic prescriptions: a time series analysis in Quebec, Canada. Antimicrob Agents Chemother. 2012;56:639–46.
Monnet DL, MacKenzie FM, López-Lozano JM, et al. Antimicrobial drug use and methicillin-resistant Staphylococcus aureus, Aberdeen, 1996-2000. Emerg Infect Dis. 2004;10:1432–41.
Sun L, Klein EY, Laxminarayan R. Seasonality and temporal correlation between community antibiotic use and resistance in the United States. Clin Infect Dis. 2012;55:687–94.
Kanwar A, Cadnum JL, Jencson AL, Donskey CJ. Impact of antibiotic treatment on the burden of nasal Staphylococcus aureus among hospitalized patients. Antimicrob Agents Chemother. 2018;62:e00609–18.
Lawes T, Lopez-Lozano JM, Nebot CA, et al. Effects of national antibiotic stewardship and infection control strategies on hospital-associated and community-associated methicillin-resistant Staphylococcus aureus infections across a region of Scotland: a non-linear time-series study. Lancet Infect Dis. 2015;15:1438–49.
Eber MR, Shardell M, Schweizer ML, Laxminarayan R, Perencevich EN. Seasonal and temperature-associated increases in gram-negative bacterial bloodstream infections among hospitalized patients. PLoS One. 2011;6:e25298.
Rodrigues FS. Clemente de Luca FA, Ribeiro da Cunha A, Fortaleza CMCB. Season, weather and predictors of healthcare-associated Gram-negative bloodstream infections: a case-only study. J Hosp Infect. 2019;101:134–41.
Davenport ER, Mizrahi-Man O, Michelini K, Barreiro LB, Ober C, Gilad Y. Seasonal variation in human gut microbiome composition. PLoS One. 2014;9(3):e90731. https://journals.plos.org/plosone/article/file?id=10.1371%2Fjournal.pone.0090731&type=printable.
The authors appreciate the generous support from Seth Peters, MPH, Rebecca Reece, MD, and Rhode Island Department of Health, Division of Preparedness, Response, Infectious Disease and Emergency Medical Services who generously gave of their time to gather and assimilate the APCD data. We also thank Jason Machan, PhD for his helpful suggestions.
Ethics approval and consent to participate
The study was approved by the Rhode Island Hospital Institutional Review Board.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Trend of positive respiratory viral panel (RVP) testing over time. (PNG 30 kb)