New Research In
The impact of COVID-19 nonpharmaceutical interventions on the future dynamics of endemic infections
Edited by Alan Hastings, University of California, Davis, CA, and approved October 13, 2020 (received for review June 24, 2020)
Significance
Nonpharmaceutical interventions (NPIs), such as social distancing, reduce not only COVID-19 cases but also other circulating infections such as influenza and RSV. The susceptible population for these infections will increase while NPIs are in place. Using models fit to historic cases of RSV and influenza, we project large future outbreaks of both diseases may occur following a period of extended NPIs. These outbreaks, which may reach peak numbers in the winter, could increase the burden to healthcare systems.
Abstract
Nonpharmaceutical interventions (NPIs) have been employed to reduce the transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), yet these measures are already having similar effects on other directly transmitted, endemic diseases. Disruptions to the seasonal transmission patterns of these diseases may have consequences for the timing and severity of future outbreaks. Here we consider the implications of SARS-CoV-2 NPIs for two endemic infections circulating in the United States of America: respiratory syncytial virus (RSV) and seasonal influenza. Using laboratory surveillance data from 2020, we estimate that RSV transmission declined by at least 20% in the United States at the start of the NPI period. We simulate future trajectories of both RSV and influenza, using an epidemic model. As susceptibility increases over the NPI period, we find that substantial outbreaks of RSV may occur in future years, with peak outbreaks likely occurring in the winter of 2021–2022. Longer NPIs, in general, lead to larger future outbreaks although they may display complex interactions with baseline seasonality. Results for influenza broadly echo this picture, but are more uncertain; future outbreaks are likely dependent on the transmissibility and evolutionary dynamics of circulating strains.
Nonpharmaceutical interventions (NPIs) have proven effective in reducing the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in many contexts (1⇓⇓⇓–5). Policy measures including social distancing, school closures, travel restrictions, and the use of masks in public spaces have been implemented to reduce the transmission of the virus. In addition to SARS-CoV-2, NPIs may also reduce the transmission of other directly transmitted, respiratory infections (6, 7). Understanding the possible influence of a SARS-CoV-2 NPI period on the incidence of these infections remains a key question for the broader public health impact of the pandemic. Furthermore, the implications of relaxing NPIs for future outbreaks of these other infections have not been fully considered.
Many endemic, directly transmitted, respiratory infections exhibit distinct seasonal and longer-term cycles in incidence (8⇓–10). While climate may drive the seasonality of these diseases in some cases (11⇓⇓–14), other directly transmitted infections, such as measles, are driven primarily by seasonal cycles of population aggregation such as the timing of school semesters (15, 16). Secular changes in susceptible recruitment, for instance, due to vaccination campaigns or declines in birth rates, can disrupt long-run patterns of infection dynamics (17, 18). Similarly, human movement via either displacement or migration has also been shown to alter patterns of infection (19). While there has been less work to identify the polymicrobial implications of nonpharmaceutical control measures, evidence from the 1918 influenza pandemic suggests that NPIs may have reduced measles transmission by 38% (20).
Two important directly transmitted, viral respiratory diseases circulating in the US population are seasonal influenza and respiratory syncytial virus (RSV). Seasonal influenza accounts for significant annual mortality, with the ongoing evolution of the virus’ antigenic sites leading to evasion of the host immune system (21, 22). Epidemics of seasonal influenza at higher latitudes are driven largely by variations in climate (12, 13). While there is some evidence of herd immunity, a complex interaction between alternating subtypes and antigenic drift determines year-to-year variation in susceptibility and corresponding outbreak size (10, 23).
RSV causes lower respiratory tract infections in young infants, and contributes to approximately 5% of under-five deaths globally (24), with no vaccine currently available. Previous models show RSV epidemics exhibit limit cycle behavior, tuned by climate-driven seasonality (Materials and Methods) (11, 25). In most regions in the United States, RSV and influenza exhibit peak incidence in the winter months, coinciding with cold, dry climatic conditions (11, 13).
Here we consider the impact of nonpharmaceutical control measures on the incidence of these two respiratory infections. We focus primarily on RSV, with the simpler limit cycle dynamics presenting an opportunity to probe interactions with NPIs. High interannual irregularity in influenza dynamics, driven by variation in circulating strains and subtypes, makes it harder to attribute the impact of NPIs. We first evaluate the influence of control measures targeting SARS-CoV-2 using influenza and RSV surveillance data. Since changes to physician visits for both viruses could be driven by behavioral responses to control measures, we look at the percent positive tests for both viruses as reported from laboratory surveillance data (see Discussion for the limitations of this measure).
Results
Fig. 1 shows the percent positive tests for RSV (Fig. 1A) and influenza (Fig. 1B) for 2019–2020 (highlighted) and four preceding years, for four states (RSV data with at least 2 y of observations were not available for other states). A national emergency in response to the COVID-19 pandemic was declared on March 13, 2020 in the United States, shown with the dashed line. Following the declaration, many states put in place control measures to limit the spread of SARS-CoV-2. Despite the declaration occurring after the typical seasonal peak in cases, a decline in prevalence is observed beyond mean seasonal levels. In Florida, where RSV cases tend to persist throughout the year (11, 25), observed RSV prevalence is reduced to near zero in March 2020. A similar pattern is visible in Hawaii for influenza, where cases are normally persistent. In Fig. 1C, we show the 2019–2020 change in percentage positive influenza tests relative to weekly mean over the previous four seasons. The 2019–2020 influenza season appears to have been more severe than average, with a relative increase in prevalence prior to March 2020 possibly driven by increased circulation of influenza subtype B (SI Appendix, Figs. S1 and S3). However, following the declaration of emergency, declines to below average levels can be observed across almost all reporting states.
Reduction in RSV and Influenza cases since March 2020. The percent positive laboratory tests for (A) RSV and (B) influenza across four US states. Data from 2020 are highlighted in red (RSV) and light blue (influenza). Data from previous seasons (2016–2019) are highlighted in gray. (C) The 2020 change relative to seasonal mean for influenza for all available US states (RSV surveillance data are only available for select states). Dashed line shows timing of the declaration of national emergency.
To explore the possible implications of control (i.e., NPI) periods for the future dynamics of influenza and RSV, we use epidemiological models and consider a range of possible scenarios for the length and intensity of control measures. Given the current uncertainty in the future course of the COVID-19 pandemic, and how responses might change over time, we cannot make precise predictions of future outcomes. For RSV, we use the time series Susceptible–Infected–Recovered (SIR) model (26, 27), fitted to historic US case data described in previous work (11). Specifically, we evaluate how NPI perturbations impact the epidemic limit cycles of RSV. We first consider a range of control period lengths and percent reduction in transmission based on Florida and Texas seasonality. We consider these two states for two reasons. First, the states exhibit distinct dynamics patterns of RSV incidence, with Florida having persistent cycles and an earlier summertime outbreak, and Texas having annual, wintertime outbreaks. Second, in terms of data availability, both these locations have recent surveillance data as well as historic case data, enabling the estimation of possible changes to transmission in 2020. Minnesota was not considered here because the percent positive surveillance data do not accurately reveal the biennial dynamics in this location (11).
Fig. 2 A and B shows the impact of these varied controls on peak incidence (I/N), peak proportion susceptible (S/N), and timing of peak I/N. Major dynamic effects are caused by a buildup of susceptible individuals as NPIs reduce transmission. Longer controls, with a greater reduction in transmission, generally lead to a greater increase in susceptibility and larger resulting outbreaks. For Florida, these outbreaks tend to occur in the summer months, but can occur throughout the year. For Texas, where seasonal transmission peaks in the winter, peak outbreaks occur only in the winter months, with the earliest outbreak in 2022.
RSV simulations for Florida and Texas. (A and B) Surface plots show the change in peak incidence per capita and peak susceptibility per capita, relative to pre-2020 maxima, for varied lengths of control (weeks) and percent reduction in transmission. Black dashed line in the first plot row shows the region above which minimum incidence drops below 1, that is, where local extinction is possible. The lower surface plot shows the timing of peak incidence in this period. Results for (A) Florida and (B) Texas are shown. (C) Simulations of future RSV epidemics, assuming a control period of 1 y and a 20% reduction in transmission, are shown for Florida and Texas. Gray block represents the NPI period, red line is proportion infected (I/N), and blue dashed line is proportion susceptible (S/N).
We use RSV laboratory surveillance data for Florida and Texas to parameterize the actual reduction in transmission caused by SARS-CoV-2 NPIs. We find that a reduction in transmission of 20% is able to conservatively capture the decline in prevalence recently observed in the surveillance data (Fig. 2C and SI Appendix, Fig. S4). Using this model parameterization, we run simulations with a control period of 1 y. Results from Florida and Texas, shown in Fig. 2C, indicate an increased likelihood of severe RSV outbreaks after the control period has ended.
We run simulations to investigate the potential impact of control measures on RSV for over 300 US counties and Mexican states using the time series SIR model fitted to historic county-level case data (Fig. 3) (11). Counties with short time series (less than 5 y of data) and sparse numbers of cases (under 10 at peak) are removed. We compare the impact of two periods of control: lasting 6 mo (Fig. 3 A and C) and lasting 1 y (Fig. 3 B and D). Although the 6-mo control period occurs outside of the peak season of the virus, substantial RSV outbreaks are still projected as a lagged response to the SARS-CoV-2 NPIs. In general, the longer, 1-y control period results in larger RSV outbreaks, yet complex interactions with seasonality arise. For New York County, the shorter control period results in a large outbreak in the following winter (2021–2022), but the longer control period results in a more persistent but less intense outbreak. In contrast, a large RSV outbreak is observed in Miami after a year of control measures. In most cases, simulated dynamics eventually return to the pre-NPI attractor. For Boulder County, by contrast, control periods have complex interactions with the seasonal biennial epidemics of the disease. In these deterministic simulations, a longer control period in Boulder County causes the epidemic trajectory to shift to a separate coexisting attractor (Fig. 3E and SI Appendix, Fig. S8). In general, the timing and size of future outbreaks will depend on the interaction between the dynamics of susceptibility and the seasonality of transmission.
RSV simulations for US counties and Mexican states. Simulations for four US counties with either (A) 6 mo or (B) 1 y of controls. Simulations for all US counties (with population > 500,000) and Mexican states in data with (C) 6-mo or (D) 1-y control period, where max incidence prior to the control period is set to 1. (E) Susceptible–Infected phase plane plot for Boulder, CO, showing epidemic trajectory with incidence time series above. The epidemic settles on a coexisting attractor postcontrol shown by the distinct precontrol (dark blue) and postcontrol (dark red) stable trajectories.
Compared to RSV, influenza epidemics exhibit a less uniform seasonal pattern. Gradual evolution of the influenza virus’ antigenic sites (antigenic drift) means population susceptibility changes over time (21), and different subtypes may circulate each year with different levels of transmission (28). In our preliminary analysis, we therefore focus on the overall dynamics of susceptibility, ignoring year-to-year differences in circulating strains. We simulate influenza using a Susceptible–Infected–Recovered–Susceptible (SIRS) model, developed in previous analyses to explore influenza seasonality in the United States, where
Fig. 4 shows the results using the influenza model under two control scenarios (6 mo and 1 y with a 20% reduction in transmission) and two transmission scenarios (high
Influenza simulations for New York County. Simulations using a (A and B) 6-mo and (C and D) 1-y control period for both (A and C) high (
Discussion
NPIs put in place to limit the spread of SARS-CoV-2 are already beginning to affect the transmission of other directly transmitted, endemic diseases. Our results suggest that a buildup of susceptibility during these control periods may result in large outbreaks in the coming years. Shorter NPIs may occur outside of the peak season of the disease but still lead to elevated future outbreaks. Longer NPIs may overlap with seasonal peak forcing, resulting in larger future outbreaks on average, but with complex transient effects (a more detailed analysis of NPI timing is shown in SI Appendix, Figs. S9 and S10). Results for RSV in the United States suggest that these outbreaks may reach their peak in the winter of 2021–2022. This finding appears robust even when we account for possible imported cases. Following perturbation, RSV generally returns to the endemic attractor, but more complex behavior is possible (Fig. 3E).
Preliminary results for influenza suggest outbreaks may occur outside of the typical season, coinciding with the end of the control period. However, we do not address complex features of the influenza virus such as circulating subtypes or the implications of global NPIs for antigenic drift. The latter may prove significant, for example, if the evolution rate first declines with NPIs, then rebounds (29). More broadly, our results suggest that healthcare systems may need to prepare for future outbreaks of non–COVID-19 infections, as NPIs are relaxed. These outbreaks may occur several years after initial NPIs were put into place.
There are several caveats to these results. First, we are at the early stages of understanding the implications of SARS-CoV-2 NPIs for endemic infections. In our model, we used a fixed reduction in transmission; however, this may not capture heterogeneities in NPIs across locations and over time. As more surveillance data become available, tracking further changes to endemic disease prevalence will be important. Serological surveys, currently used to measure exposure and potential immunity to SARS-CoV-2, could similarly be employed to monitor these polymicrobial responses (30⇓–32).
Second, and importantly, an influx of COVID-19 cases could artificially lower the percent positive test data we use to calibrate reductions in RSV and influenza transmission. For RSV, this seems unlikely, as the mean age of infection is much lower than for COVID-19, and cases are unlikely to overlap. Most individuals are infected with RSV before the age of 2 y (18). In contrast, infants appear less likely to be infected with SARS-CoV-2, and severe respiratory presentations are unusual (33, 34). We also note that total tests for RSV have declined since March (SI Appendix, Fig. S2). For influenza, it is more plausible that an influx of COVID-19 cases could be biasing results; however, the sharp decrease observed across states right after the national emergency declaration suggests that it is unlikely to be primarily driven by this factor. The total influenza specimens collected have also declined (SI Appendix, Fig. S1). As a further check, we create an incidence proxy, following ref. 35, by multiplying percent influenza-like illness by percent positive tests (SI Appendix, Fig. S11). Substantial declines in this measure are still observed following the emergency declaration.
Possible biases could inflate our measure of the percent reduction in transmission due to NPIs (estimated at 20%). However, similarly, the relatively limited time horizon of our data, and the fact that current NPIs are concurrent with the seasonal lower transmission period, could mean we are underestimating this reduction. Recent evidence from Hong Kong estimated a 33 to 44% reduction in influenza transmission due to NPIs (7). We therefore stress that our results are uncertain, and we cannot yet make precise predictions. Unbiased polymicrobial surveillance would enable better attribution of NPI effects and is a crucial area for future development. In addition, interactions between the SARS-CoV-2 virus and endemic viruses may be more complex than described here. Immunological relationships between viruses, both competitive and cooperative, may have broad-scale implications for future infection dynamics (36). The impact of NPIs on strain structure of RSV (37) is an important area for future work.
Finally, although we have primarily focused on the United States, outcomes may be more severe in Southern Hemisphere locations where NPI timing aligns with the peak season for seasonal wintertime diseases. Our results also illustrate the potential for COVID-19 NPI to impact the dynamics and persistence of a much wider range of infections. Increased surveillance, serological surveys, and local modeling efforts will help determine the future dynamics and risk from these circulating infections.
Materials and Methods
Data.
Recent (2016–2020) disease data based on laboratory results from either antigen or PCR tests for RSV are obtained from the corresponding government websites for each state: Florida (unspecified) (38), Minnesota (antigen) (39), Oregon (antigen and PCR) (40), and Texas (antigen) (41). Even though a few other states report RSV surveillance, we do not include them in our analysis, as they do not provide information on RSV circulation from previous years. Some RSV data are extracted from the graphs of the state surveillance reports, as raw values are unavailable. Influenza surveillance data are obtained from Centers for Disease Control and Prevention FluView Interactive (42). Historic RSV data (pre-2010) used to train the RSV model come from hospitalizations data originally obtained from the State Inpatient Databases of the Healthcare Cost and Utilization Project maintained by the Agency for Healthcare Research and Quality (43). Population data for the United States are obtained from publicly available combined files of US Census Bureau data available via the National Bureau of Economic Research (44). US birth data are downloaded from the Centers for Disease Control (45). Transmission in the influenza model relies on specific humidity data taken from NASA’s Modern-Era Retrospective analysis for Research and Applications version 2 dataset (46, 47).
Models.
We first calculate location-specific seasonal transmission rates using the time series Susceptible–Infected–Recovered model (TSIR), a discrete time adaptation of the SIR model (26, 48). County-level transmission rates were calculated for a previous study (11). The TSIR model describes the number of infected and susceptible individuals as a set of difference equations. The number of susceptible individuals is given by
Using this formulation, it is shown that a linear regression of cumulative births on cumulative cases gives
We generate forward simulations using county seasonal transmission rates,
The percentage reduction in transmission is estimated by comparing model simulations to laboratory RSV data from 2020 for Texas and Florida. State-level data from Minnesota are not used because the laboratory data does not capture the biennial cycles of incidence that exist in this state (11). Other states do not have multiple years of available data to compare present reductions in prevalence. Laboratory test data are scaled to the model projection using the 2016–2020 mean, that is,
For influenza we use a climate-driven SIRS model (13, 14, 28). Antigenic drift of the influenza virus results in a seasonal return to susceptibility, meaning TSIR methods are not appropriate for this infection. The model is described by a series of differential equations (51),
Data Availability.
Data for this study come from publicly available datasets. Code and collated data to recreate the main results are available via GitHub at https://github.com/rebaker64/NPIs.
Acknowledgments
R.E.B. is supported by the Cooperative Institute for Modeling the Earth System. This study is supported by the Princeton Institute for International and Regional Studies and the Princeton Environmental Institute.
Footnotes
- ↵1To whom correspondence may be addressed. Email: racheleb@princeton.edu.
Author contributions: R.E.B., S.W.P., C.J.E.M., and B.T.G. designed research; R.E.B. and S.W.P. performed research; R.E.B., W.Y., and G.A.V. analyzed data; and R.E.B., S.W.P., W.Y., G.A.V., C.J.E.M., and B.T.G. wrote the paper.
The authors declare no competing interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.2013182117/-/DCSupplemental.
Published under the PNAS license.
References
- ↵
- S. Pei,
- S. Kandula,
- J. Shaman
- ↵
- H. Tian et al.
- ↵
- S. Lai et al.
- ↵
- M. Chinazzi et al.
- ↵
- S. Hsiang et al.
- ↵
- N. H. Leung et al.
- ↵
- B. J. Cowling et al.
- ↵
- M. E. Martinez
- ↵
- C. Viboud,
- W. J. Alonso,
- L. Simonsen
- ↵
- C. Viboud et al.
- ↵
- R. E. Baker et al.
- ↵
- J. Shaman,
- M. Kohn
- ↵
- J. Shaman,
- C. Viboud,
- V. E. Pitzer,
- B. T. Grenfell,
- M. Lipsitch
- ↵
- R. E. Baker,
- W. Yang,
- G. A. Vecchi,
- C. J. E. Metcalf,
- B. T. Grenfell
- ↵
- B. T. Grenfell,
- O. N. Bjørnstad,
- J. Kappey
- ↵
- M. J. Ferrari et al.
- ↵
- B. Bolker,
- B. T. Grenfell
- ↵
- V. E. Pitzer et al.
- ↵
- A. Wesolowski et al.
- ↵
- A. D. Becker,
- A. Wesolowski,
- O. N. Bjørnstad,
- B. T. Grenfell
- ↵
- R. Webster,
- W. Laver,
- G. Air,
- G. Schild
- ↵
- V. N. Petrova,
- C. A. Russell
- ↵
- D. J. Smith et al.
- ↵
- H. Wang et al.
- ↵
- V. E. Pitzer et al.
- ↵
- B. F. Finkenstädt,
- B. T. Grenfell
- ↵
- B. T. Grenfell,
- O. N. Bjørnstad,
- B. F. Finkenstädt
- ↵
- W. Yang,
- M. Lipsitch,
- J. Shaman
- ↵
- N. Arinaminpathy et al.
- ↵
- A. K. Winter et al.
- ↵
- M. J. Mina et al.
- ↵
- M. J. Mina et al.
- ↵
- J. B. Dowd et al.
- ↵
- H. Hong,
- Y. Wang,
- H. T. Chung,
- C. J. Chen
- ↵
- E. Goldstein,
- S. Cobey,
- S. Takahashi,
- J. C. Miller,
- M. Lipsitch
- ↵
- S. Nickbakhsh et al.
- ↵
- L. White et al.
- ↵
- Florida Department of Health
- ↵
- Minnesota Department of Health, Weekly influenza & respiratory activity: Statistics
- ↵
- Oregon Health Authority
- ↵
- Texas Health and Human Services
- ↵
- Centers for Disease Control and Prevention
- ↵
- Agency for Healthcare Research and Quality
- ↵
- National Bureau of Economic Research
- ↵
- Centers for Disease Control and Prevention
- ↵
- R. Gelaro et al.
- ↵
- National Aeronautics and Space Administration
- ↵
- O. N. Bjørnstad,
- B. F. Finkenstädt,
- B. T. Grenfell
- ↵
- K. Glass,
- Y. Xia,
- B. T. Grenfell
- ↵
- A. D. Becker,
- B. T. Grenfell
- ↵
- J. Dushoff,
- J. B. Plotkin,
- S. A. Levin,
- D. J. Earn
We recommend
-
Inference of seasonal and pandemic influenza transmission dynamics
Wan Yang et al., Proc Natl Acad Sci U S A, 2015
-
Forecasting the spatial transmission of influenza in the United States
Sen Pei et al., Proc Natl Acad Sci U S A, 2018
-
The seasonality of nonpolio enteroviruses in the United States: Patterns and drivers
Margarita Pons-Salort et al., Proc Natl Acad Sci U S A, 2018
-
Moving beyond metagenomics to find the next pandemic virus.
Vincent Racaniello, Proc Natl Acad Sci U S A, 2016
-
SARS-like WIV1-CoV poised for human emergence
Vineet D. Menachery et al., Proc Natl Acad Sci U S A, 2016
-
Large, delayed outbreaks of endemic diseases possible following COVID-19 controls
MedicalXpress, 2020
-
Qiagen Obtains CE Mark for Multiplex Respiratory Infection Test, Saliva-Based SARS-CoV-2 Assay
Ben Butkus et al., GenomeWeb
-
Eurobio Gets CE Mark for Multiplexed SARS-CoV-2, Influenza A/B, RSV Test
staff reporter, 360Dx, 2020
-
Will the COVID-19 virus become endemic?
MedicalXpress, 2020
-
Management of Renal Carcinoma During Pregnancy
The Doctor's Channel, 2019
Sign up for Article Alerts
Article Classifications
- Biological Sciences
- Ecology