ABSTRACT
Objective Using a pandemic influenza model modified for COVID-19, this study investigated the degree of control over pre-symptomatic transmission that common non-pharmaceutical interventions (NPIs) would require to reduce the spread in long-term care facilities.
Methods We created a stochastic compartmental SEIR model with Poisson-distributed transition states that compared the effect of R0, common NPIs, and isolation rates of pre-symptomatic carriers primarily on attack rate, peak cases, and timing in a 200-resident nursing home. Model sensitivity was assessed with 1st order Sobol’ indices.
Results The most rigorous NPIs decreased the peak number of infections by 4.3 and delayed the peak by 9.7 days in the absence of pre-symptomatic controls. Reductions in attack rate were not likely, even with rigorous application of all defined NPIs, unless pre-symptomatic carriers were identified and isolated at rates exceeding 76%. Attack rate was most sensitive to the pre-symptomatic isolation rate (Sobol’ index > 0.7) and secondarily to R0.
Conclusions Common NPIs delayed and reduced epidemic peaks. Reducing attack rates ultimately required efficient isolation of pre-symptomatic cases, including rapid antigen tests on a nearly daily basis. This must be accounted for in testing and contact tracing plans for group living settings.
INTRODUCTION
Despite being less than 1% of the population of the United States, 8% of all confirmed cases and more than 41% of recorded COVID-19 deaths in the US have been linked to nursing homes and other long-term care facilities (LTCFs). Caring for our aging population has seldom been a top priority1,2 even though the proportion of US elderly residing in LTCFs is greater than anywhere else in the world. Nursing homes and residential care communities are staffed by more than 1 million nurses, aides, and social workers, and are home to more than 2.2 million residents.3 Most residents are elderly and thus at a higher risk of infection and mortality from COVID-19.4–6 In this population, the prevalence of advanced age and significant comorbidities such as diabetes and hypertension substantially increases the risk of hospitalization, severe symptoms, and death.7–11
Since the 1918 influenza pandemic there has been extensive discussion about how to best prepare for another pandemic,12–16 with some special considerations for LTCF residents.17–20 In 2018, the World Health Organization (WHO) released the Global Influenza Strategy 2019-2030 with the goals of reducing the burden of seasonal influenza, minimizing the risk of zoonotic influenza, and mitigating the effects of pandemic influenza. These efforts have provided limited guidance on the implementation of non-pharmaceutical interventions (NPIs) for controlling the spread of pandemic infection.16 Controlling exposure pathways through social distancing measures, personal protection measures, and community-based interventions is critically important, uniquely so during the period of time between pandemic pathogen emergence and the development and distribution of therapeutics and vaccines.
In 2008, motivated by mounting concern about pandemic influenza, we produced a mathematical model that anticipated some of the challenges that nursing homes would face in the event of a new influenza pandemic prior to vaccine deployment. The most effective modeled plan required virtually complete facility isolation, complete visitation restriction, rotating staff schedules with a variable isolation period prior to employee reentry, and stringent mitigation protocols imposed on all high-priority services entering a facility.21 Without access to vaccines or other therapeutics, we concluded that the effectiveness of such plans mostly depended on preventing staff from introducing an infection acquired outside of the facility, as mitigation after introduction is rarely effective. The successful containment of SARS/MERS and the relatively low virulence of the 2009 H1N1 influenza virus made the considerable investment and social disruption required for such a pandemic plan difficult to justify.
COVID-19 has dramatically changed that calculation. Though two effective vaccines are on the horizon, control deficits in many places in the US and exponential increases in infection rates in the winter of 2020,22,23 it is critical to identify opportunities and mechanisms that reduce the risk of introducing or re-introducing COVID-19 into LTCFs. While incubation period of influenza viruses are typically 2-3 days, that of COVID-19 appears to be between 5 and 9 days, with a small number of cases taking 12+ days to manifest symptoms.24–28 Our models for influenza recommended the isolation of staff for three or more days at the end of off-periods to allow symptoms to manifest after exposure. However, the longer and uncertain incubation period of COVID-19 would require staff to quarantine for up to 2-week periods at the end of off-duty stints to ensure effectiveness. This is a non-starter. Additionally, serial interval studies provide mounting evidence of a period of pre-symptomatic infectiousness of possibly multiple days, during which an infected individual might unknowingly expose others.29–32 Pre-symptomatic transmission has likely contributed significantly to COVID-19 spread in many LTCFs across the United States and is a significant challenge that NPIs based on influenza transmission dynamics do not address.33–37
In this study, we adapted our previous stochastic SEIR model of a pandemic influenza outbreak in a nursing home to reflect the longer incubation period and significant pre-symptomatic transmission risk of COVID-19. We then characterized the performance of four distinct NPI scenarios across a range of reproductive numbers and isolation rates for pre-symptomatic individuals to identify approaches that most reduced the attack rate as well as the timing and magnitude of the epidemic peak.
METHODS
Nuño et al. 200820 analyzed how effectively four different NPI scenarios prevented the introduction of infection to a nursing home resident population during a hypothetical local outbreak of pandemic influenza. We adapted the stochastic, compartmental Susceptible-Exposed-Infected-Recovered (SEIR) model from this previous study to reflect our current knowledge of progression and transmission parameters for SARS-CoV-2, including significant modifications for residents, to model outcomes in a nursing home with 200 susceptible residents and 83 total staff. Multiple NPI scenarios with varying levels of rigor and incremental pre-symptomatic isolation rates were applied across viable values for R0 to see how the attack rate, mortality rate, number of hospitalizations, and number of ICU admissions responded. Each intervention category described by the original model was restructured and re-parameterized to reflect the current suggested NPI policies for nursing homes.4,38 NPIs modeled include isolation of symptomatic residents and barring symptomatic staff from work; staff and visitor use of PPE when interacting with residents; temperature monitoring for incoming staff and visitors; staff schedules that account for incubation and isolation periods; and restrictions on visitor entry. Continuous outcomes among NPI scenarios were assessed for significance by calculating the difference between two scenarios and applying one-sided t-tests with a null hypothesis set to a difference of 0.
Mathematical Model
A detailed accounting of the structure and parameterization of the model can be found in the Supplemental Digital Content. Four main actors were considered: residents, who could only be exposed in the facility and move to the hospital or ICU; staff, who could be exposed in the community or in the facility; visitors, who had limitations on visiting hours, and could be exposed in the community or the facility; and community members, within which the local outbreak begins. All actors within the facility at a given time were assumed to be mixing randomly, as were all actors in the community. Though permanence of immunity is still unknown, we assume that recovery provides complete immunity for the remainder of the 200-day simulation.
In the SEIR model, where i,j are the actor type (resident, staff, visitor, community member) and their location (facility or community) respectively, all who can be infected begin in the susceptible compartment Sij. When successfully exposed, they move to the asymptomatic/pre-symptomatic silent infection compartment, Eij, where they can infect other susceptible actors. Some develop symptoms, moving to the infected compartment Iij, while others do not and ultimately recover and are removed from the modeling population (complete immunity). Since hospital and ICU capacity are critical concerns in the COVID-19 pandemic, we added two new compartments to the resident outbreak model: infected residents who required hospitalization move to the hospital compartment,, and residents who require intensive care move to the ICU compartment,. Residents could suffer a mortality event in both and, which would remove them from the modeling population. The parameters and definitions presented in Table [1] reflect the known characteristics of the SARS-CoV-2 virus.
A full schematic of the model is available in Supplemental Figure [S1]. State transition rates for residents in the facility are highlighted in Equations [1], below. Exposure rates were related to the infection force, λ, which modifies the transmission rate and transmission reduction from protective measures by taking transient actors into account (Equation [2]).
Baseline Scenario
The Baseline scenario describes the case where no specific preventative action is taken. Symptomatic residents are not isolated, symptomatic staff are allowed to report to work, staff and visitors are not required to use PPE in normal interactions with residents, temperature and other symptom checks are not required for entry, visitors are allowed two hours of visitation time per week, and staff have a normal 5 days on/2 days off work week.
NPI Categories 1-2
Categories 1-2 are defined by isolating symptomatic individuals, light mitigation of spread within the facility, and no restrictions for entrance into the facility by asymptomatic individuals. Isolation of symptomatic residents is required, and symptomatic staff are told to stay home. PPE is required for staff and visitors interacting with residents. Otherwise, entry to the facility does not require symptom checks, staff maintain their standard schedules, and visitation hours remain at baseline levels.
NPI Categories 3-4
Categories 3-4 apply all measures in Categories 1-2 in addition to tighter restrictions on entry to the facility. Temperature checks are done for anyone entering the facility. Staff switch to 12-hour shifts with 3 days off at a time to allow for more time for symptoms to develop before returning to work. Visitor hours are reduced to half of Baseline.
NPI Category 5
Category 5 implements the most restrictive versions of the previous NPIs. Category 1-2 restrictions are still in place and temperature checks are required for entry. Staff switch to a 4 days on/4 days off schedule to further extend time away from residents. Visitors are completely disallowed. Complete specifications and parameterizations for each category are presented in Supplemental Table [S5].
Solving the Model
The literature reports that SARS-CoV-2 has a basic reproduction number, R0, that ranges from 1.5 to 6.47 in the early stages of an outbreak.39 We narrowed our analysis down to an interval of values from 2.0 to 4.0 in increments of 0.2 by varying β from 0.61 to 1.21 in intervals of 0.06. Given the strong evidence that silent infections facilitate the rapid spread of the virus,40–43 we repeated each simulation at 0% to 90% (in increments of 10%) of exposed individuals being identified and isolated, simulating actors isolating themselves after contact with a known positive individual regardless of symptom status or an increasingly stringent testing regimen.44 Since the model facility has a finite population and the outcomes of interest were based on counts, stochasticity was introduced by applying a Poisson distribution to all state transitions.45 100 iterations were run for each combination of inputs to quantify the uncertainty. Simulations were executed in MATLAB (vR2020a, The MathWorks, Inc. 2020).
Sensitivity Analysis
The package ODEsensitivity (v1.1.2; Weber, Theers, and Surmann, 2019) for R (v3.6.3; R Core Team 2020) was used to calculate first-order Sobol’ indices for each primary model output with respect to the two main inputs of concern, R0 and asymptomatic isolation rate \cite{wu2013sensitivity}46. The Sobol’ method is a generalized sensitivity analysis (GSA) that provides robust variance comparisons independent of the linearity or monotonicity of a model; outputs used here are the first-order Sobol’ index, which is the ratio of a model’s variability from the parameter alone to the total variability in the full model.
RESULTS
Four different combinations of increasingly restrictive NPIs (Baseline, Categories 1-2, Categories 3-4, and Category 5) were simulated with incremental isolation rates for asymptomatic cases (0%, 50%, 60%, 70%, 80%, and 90%), each with a R0 ranging from 2.0 to 4.0 in intervals of 0.2. For each scenario, 8 outcomes were calculated: total resident infections and attack rate, total resident hospitalizations, total resident ICU admissions, total resident deaths and mortality rate, number of symptomatic infections at peak, and the timing of peak infections.
Figure [1] displays attack and mortality rates as well as hospital and ICU admissions observed as model outputs across R0 values for the four NPI scenarios at 90% pre-symptomatic isolation. Similarly, Figure [2] also display attack rates but at a high-resolution intervals of 2% for pre-symptomatic isolation rates between 70% and 90%. A complete accounting of mean values and 95% confidence intervals (95% CIs) are presented in Supplementary Tables [S2-S4]. All parenthetical ranges subsequently reported represent 95% CIs. Notably, when the pre-symptomatic isolation rate was 0% (i.e., no isolation), the four NPI scenarios were indistinguishable in terms of the total infections, hospitalizations, ICU admissions, and deaths (Supplementary Figure [S3] and Supplementary Table [S2]).
Attack Rate
In the baseline scenario, total resident infections remained close to 190, or a 95% attack rate, regardless of R0: 189.6 (188-196) at R0 = 2.0 and 190.8 (189-195) at R0 = 4.0. At a 0% pre-symptomatic isolation rate, no NPI scenario showed significant reductions in the attack rate. At 50% isolation, Categories 1-2, Categories 3-4, and Category 5 scenarios all led to lower total infections for R0 < 3.0, but rapidly deteriorated for higher R0 values (Figure [2] and Supplementary Table [S3]). All scenarios showed large reductions in attack rate at 90% isolation (Figure [1] and Supplementary Table [S2]). Hospitalizations, ICU admissions, and mortality rate were all reduced as attack rate fell. Where differences were observable among these four scenarios, the Category 5 scenario outperformed all predecessors. Figure [2] shows that such large reductions only occurred after higher isolation rates were reached. In the worst-case scenario of R0 = 4.0, for instance, reducing attack rate to 80% required isolation rates of 76% for Category 5, 80% for Categories 3-4, and 84% for Categories 1-2. The lowest attack rate achieved by our model during the worst-case scenario of R0 =4.0 was 43.2% (38.3% - 59.0%) at 90% isolation using Category 5 interventions. See Supplementary Table [S5] for the full range of calculated values.
Peak Infections and Peak Times
At 0% isolation all three NPI scenarios showed significant reductions in peak symptomatic infections compared to Baseline. Across all values for R0, peak numbers of infections were reduced by an average of 3.2 (0.8-5.6), 4.0 (1.5-6.6), and 4.3 (2.3-6.3) infections for Categories 1-2, Categories 3-4, and Category 5, respectively. When compared to baseline across all values of R0, peak times were delayed by an average of 3.6 (1.0-6.2), 4.8 (1.3-8.4), and 9.7 (4.7-14.6) days for Categories 1-2, Categories 3-4, and Category 5, respectively. While Categories 1-2 and Categories 3-4 were not significantly different from each other, Category 5 consistently showed longer delays than the other two (p = 0.0004 and 0.001). See Supplementary Figure [5] for more details.
As isolation rates increased incrementally, while Categories 1-2 and Categories 3-4 frequently did not show significant differences from each other (α < 0.05). Category 5 consistently had the largest decreases in peak size and delays to peak time at 90% isolation rate: the peak size had been reduced by an average of 52.9 (49.8 - 56.0) symptomatic infections, and the peak time delay increased to 32.6 (25.5 - 39.6) days. Except for the difference in peak time delay between 80% and 90% isolation (p = 0.496) the performance of every NPI Category was significantly improved over previous isolation rates (α < 0.05).
Sensitivity Analysis
Results of the Sobol’ method sensitivity analysis can be found in Supplementary Figure [S2]. First-order Sobol’ indices were calculated for R0 and the asymptomatic isolation rate with respect to attack rate. Variability in attack rate was more sensitive to the isolation rate for the first 125 days of the outbreak, with a maximum value of over 0.70 at day 27. Similarly, variability in hospitalizations was more sensitive to isolation rate for the first 155 days of the outbreak, with a maximum value of almost 0.75 before day 50. Variability in ICU admissions followed a pattern similar to hospitalizations. Variance in mortality rate, however, was more sensitive to changes in R0 than isolation rate, except for during the early days of an outbreak where contributions were roughly equal.
LIMITATIONS
There are several limitations that are worth considering with the proposed modeling approach presented in this study. Our model assumes that coronavirus transmission occurs from direct human-to-human contact; however, studies have shown that indirect transmission through surfaces and shared air circulation systems can also occur. While we incorporated stochastic fluctuation in the process of transmission, our model does not consider differential rates of transmission due to changes in infection control measures or individual-based variation and contact structure. Neither age nor comorbidity status was included as predictors of transmission, hospitalization, or death, so agents only differed by access to specific compartments. Thus, results should be read in terms of proportional change and not absolute values. Our mortality rates were generally lower than what has been observed in nursing home outbreaks across the US.
DISCUSSION
Community outbreaks were almost always followed by an outbreak in our model, indicating that even the most stringent of NPI scenarios are unlikely to prevent introduction altogether. Baseline conditions invariably lead to an attack rate ≥ 95%, neither varying with NPI scenario nor with the pre-symptomatic isolation rate as defined. Our NPI scenarios, which were parameterized to reflect commonly endorsed control methods in nursing homes and other group living situations, all failed to consistently prevent the virus from entering the facility (average attack rates significantly > 0% in all simulations). Even the most restrictive intervention plan, Category 5, did not make significant differences in attack or mortality rate until 76% of asymptomatic cases were successfully isolated, where a minimal reduction in attack rate to 80% was achieved in the worst case scenario. Thus, even the most stringent social controls modeled are undermined by a long period of pre-symptomatic infectivity. This is consistent with assessments made on many individual nursing homes, which highlighted empirical evidence that pre-symptomatic transmission appears to be the largest contributor to their failure to control an outbreak.33–37
The above is not to be taken as an argument that these NPI scenarios should be jettisoned. While it is true they were unable to reduce the size of the outbreak without fastidious control of pre-symptomatic carriers over the 200 days of simulation, they significantly modified epidemic dynamics in terms of timing and shape of the epidemic curve. In fact, even with no ability to identify or modify pre-symptomatic isolation, all scenarios significantly reduced and delayed the infection peak compared to Baseline. Since each Category expands upon the previous lower-ordinal category, it is not a surprise that the model demonstrated increasing effectiveness with increasing rigor. However, it is notable that not all amplifications resulted in differences that were significant. Comparing the scenarios to each other, the quantitative difference in impact between Categories 3-4 and Category 5 was larger than that between Categories 1-2 and Categories 3-4. Absent the effects of isolation rate and R0, longer periods of off-duty time provided to staff, greater detection of symptomatic individuals, and the blanket ban on visitors remain important contributors to mitigation. Even with an isolation rate of 0%, Category 5 bought the facility about a week of time delay for the epidemic peak. In the real-world, an extra week to act once an outbreak has started allows for additional, possibly novel reduction measures to be implemented, and reduces the immediate pressure on local hospital systems, even if the total number of hospitalizations is not reduced. NPIs, as they are currently conceived and implemented, remain important to our ability to “flatten the curve.” Quantification of such benefits is beyond the scope of this model and has been addressed elsewhere.13,47–49
Inclusion of pre-symptomatic isolation rate as an input unbound to the NPI categories allowed us to explore variations in outcome independent of other NPI elements. Our scenarios would likely have some inherent isolation rate in real world applications. However, the high isolation rate required for any of the NPI scenarios to effectively reduce the attack rate demonstrates that successful mitigation strategies must include a method for identification of exposed individuals of sufficient sensitivity. Frequent and rapid testing for the presence of viral antigen using a technology and implementation capable of sensitivity at or above the indicated level of isolation is essential. This level of identification must be accompanied by facility and staff support to accommodate isolation periods for individuals who have had contact with known infected persons. Without the ability to implement testing and contact tracing in a way that upwards of three-quarters of potential asymptomatic, infectious carriers are continuously discovered, even current best-practices will see only small reductions on the total impact on the facility. Thus, proactive actions to isolate potential cases and identify asymptomatic cases, combined with rigorous control measures appear to be required to prevent the rapid spread of the SARS-CoV-2 virus in group settings.36,40,50
Isolation Rate and Rapid Antigen Tests
Our model demonstrates the importance of planning a tailored pre-symptomatic isolation rate for the level of rigor at which other NPIs are implemented. A Category 5 approach in the model should aim to achieve an absolute minimum of 76% isolation rate of pre-symptomatic individuals; if a less rigorous approach (lower Category) is all that is available, or a higher reduction in attack rate is desired, the burden of isolation becomes substantially higher. A perfectly sensitive test would require testing the target percentage of all person-days at the facility (followed by perfect isolation procedures); however, the four rapid antigen tests approved for use in the United States by the FDA under the Emergency Use Act (EUA) do not offer perfect sensitivity: 97.1% (85.1%-99.9%; BinaxNOW COVID-19 Ag Card), 97.6% (91.6%-99.3%; LumiraDx SARS-CoV-2 Ag Test), 96.7% (83.3%-99.4%; Sofia SARS Antigen FIA), and 84% (67% - 93%; BD Veritor System) compared to PCR tests. Before rigorous evaluation of baseline measurements of success can be done (such as positive predictive value, or PPV), a conservative sampling plan should use the lower bounds of possible sensitivities. Generally, there is a simple relationship between the target isolation rate, Itarget, the lower bound test sensitivity seLB, and the proportion of person-days that testing needs to be done to achieve that rate, PPD (Equation [3]).
If PPD =1, every resident and staff must be tested daily to reach the target isolation rate with no room for error; PPD > 1 means the pre-test probability of false negatives likely prevents that test from achieving the desired pre-symptomatic isolation rate; and PPD<1 indicates that testing all person-days may not be necessary to achieve the desired isolation rate. The lower bounds of sensitivities for the available tests are displayed in Figure [2] along with a heat map of the average attack rate at different isolation rates; the region to the left of a given line highlights the isolation rates achievable with that test in a conservative plan. Note that a 76% isolation rate already rejects the use of the BD Veritor System as insufficiently sensitive, while the LumiraDx test would reliably achieve desired isolation rates up to the 90\% investigated in this model. Utilizing a less rigorous approach would increase the isolation rate required to achieve the same reduction in attack rate.
The specificity of tests on offer appears generally quite high. Falsely identifying an LTCF resident as positive would be of little concern if the resident were merely isolated. If this information were used to transfer such residents to a consolidated facility such as a ward or separate building, then the matter becomes more consequential.
CONCLUSIONS
Though the absolute attack rates discovered in this model should not be used as hard targets for choosing which interventions to apply, isolation rates and their relationships to model outputs can be found in Figure [2]; and these represent the form, if not the precise detail of implementation of models such as these.
These results stand in stark comparison with the results of our foundational model on pandemic influenza, which found significant reductions in impact from even Category 1-2 responses.20 Adapting this model to COVID-19, and especially incorporating the extended pre-symptomatic period, has brought to light the requirements that all of high levels of detection, isolation of asymptomatic carriers, and rigorous applications of suggested NPIs (Category 5) are necessary for control. Results here should be taken as a demonstration of relative importance of NPIs and isolation of pre-symptomatic cases but not as predictions for absolute magnitude for real-world facilities. However, these results are consistent with empirical reports from facilities in high-risk areas, in which the absence of control measures and most particularly, merely identifying and isolating only symptomatic cases, a COVID-19 outbreak will likely spread rapidly within all long term care facilities.6,33,36,37
Data Availability
The full MATLAB code and all required materials to replicate the simulations from this mathematical modeling paper are publicly available at https://doi.org/10.5281/zenodo.4073575
Acknowledgements
We thank Claire and Richard Griffith for generously providing real-world input, suggestions, and discussions of the work presented here. We also thank Baltazar Espinoza Cortes for providing analytical support and advice.
Footnotes
Conflict of Interest: None declared.
Source of Funding: This research was supported by a Seed Fund Award 2020-0000000022 from CITRIS and the Banatao Institute at the University of California.
Data Access: The full MATLAB code and all required materials to replicate the simulations from this mathematical modeling paper are publicly available at https://doi.org/10.5281/zenodo.4073575
LIST OF ABREVIATIONS
- COVID-19
- Coronavirus Disease 2019
- EUA
- Emergency Use Act
- FDA
- Food and Drug Administration
- GSA
- Generalized Sensitivity Analysis
- LTCF
- Long-Term Care Facility
- NPI
- Non-Pharmaceutical Intervention
- PCR
- Polymerase Chain Reaction
- PPE
- Personal Protective Equipment
- SARS-CoV-2
- Severe Acute Respiratory Syndrome Coronavirus 2
- SEIR
- Susceptible Exposed Infected Recovered
- WHO
- World Health Organization