Introduction

SARS-CoV-2 has created a pandemic in which morbidity and mortality have been partially mitigated in many areas by widespread vaccination. COVID-19 vaccines have been extremely effective at preventing severe disease (vaccine efficacy, VE > 90%1), while also reducing susceptibility to infection (VES) and risk of onward transmission (VEI). In spite of these reductions, so-called vaccine breakthrough infections and subsequent transmission have been widely documented2, and have increased dramatically with the emergence of the omicron variant in late 20213,4. These developments raise the question of how to best mitigate transmission in partially vaccinated populations.

Prior to the approval of COVID-19 vaccines, transmission mitigation via regular and repeated screening was shown to be an effective approach to break chains of transmission and decrease the burden of COVID-19 using both RT-PCR5,6,7 and rapid antigen testing7,8. Specifically, screening via testing, which we hereafter refer to simply as screening in most cases, is effective at the community level because it decreases transmission from individuals who are already infected7,9. However, policy proposals in 2021 and early 2022 shifted to focus routine testing requirements on only the unvaccinated population, including an Italian requirement announced in October, 202110 and a U.S. requirement for healthcare workers beginning February, 202211. By reducing rates of transmission from only the unvaccinated population, such policies may be limited by the extent to which transmission is, in fact, driven by the unvaccinated. This raises critical questions about projected policy impacts relative to other non-pharmaceutical interventions (NPIs)12,13—particularly in areas where the unvaccinated population is small.

The role of vaccines in reducing transmission is complex and changing. First, VES and VEI vary depending on which vaccine was administered14. Second, both VES and VEI wane with time since vaccination15,16,17, but can be boosted to higher levels for those receiving an additional dose18. Third, those who have experienced a SARS-CoV-2 infection also show decreased risks of reinfection and subsequent transmission14, providing partial protection to those who are previously infected and remain unvaccinated, but also augmenting protection for those who are both vaccinated and previously infected18,19. Finally, preliminary estimates of VES and VEI, and their prior infection equivalents, are markedly lower for the omicron variant3,20. Thus, the relative estimates of risk reductions due to vaccination, prior infection, or both, as well as the sizes of the populations falling into each category of immunity, will affect transmission dynamics—with or without testing.

In this study, we model the spread of SARS-CoV-2 in populations of mixed vaccination status, focusing on three critical questions. First, how do vaccinated and unvaccinated populations each contribute to community spread and hospitalizations, and how do those contributions vary with rates of vaccination and prior infection? Second, how do testing-based screening programs focused on unvaccinated individuals alone affect community spread and hospitalizations? Third, how effective are delta-era screening strategies likely to be against variants with higher breakthrough and reinfection rates? Our study’s goals are not to make perfectly calibrated predictions but instead to elucidate more general principles of transmission and unvaccinated-only testing in partially vaccinated populations. As such, our analyses consider a wide range of parameters and scenarios.

Results

Unvaccinated-only testing-based screening programs have been discussed and implemented during transmission of both the delta and omicron variants, yet these variants differ in their transmission and disease parameters, particularly among vaccinated or previously infected individuals—the focus of the present study. As such, the analyses that follow incorporate a range of empirically estimated parameters for the delta variant and plausible parameters associated with the omicron variant.

High vaccination rates drive total infections and hospitalizations down, increase the proportions of vaccine breakthroughs, and shift the drivers of transmission

To examine the dynamics of transmission in a population with mixed vaccination status, we first modeled transmission within and between communities of vaccinated (V) and unvaccinated (U) individuals in the absence of a screening program. Based on a standard susceptible exposed infected recovered (SEIR) model, we tracked the four transmission modes by which an infection might spread: U → U, U → V, V → U, and V → V (Fig. 1a). A constant and equivalent fraction of both populations was assumed to have experienced prior SARS-CoV-2 infection, resulting in four categories of imperfect immunity: unprotected (unvaccinated with no prior infection), infection-acquired, vaccine-acquired, and both vaccine- and infection-acquired (so-called “hybrid” immunity). To account for introductions of infection from outside the population, all susceptible individuals were subject to a small, constant rate of exposure, with infection-acquired and vaccine-acquired immunity providing partial protection against subsequent infection.

Fig. 1: Vaccination affects which population drives transmission and dominates infections and hospitalizations.
figure 1

a Diagram of four transmission modes within and between vaccinated and unvaccinated communities, where vaccines and prior infection decrease risks of infection, transmission, and hospitalization. b Number of infected individuals and c new daily hospitalizations over time (solid black), stratified by unvaccinated (dashed gray) and vaccinated (solid gray) populations. d Daily transmission events separated and colored by transmission mode (see legend). e Cumulative infections and f hospitalizations as a percentage of total infections/hospitalizations without vaccination (black), and the percent of each accounted for by vaccine breakthroughs (gray) for varying vaccination rates. g Transmission mode (see legend) as a percentage of cumulative infections for varying vaccination rates. Black arrows in eg indicate vaccination rate at which Reff = 1; vertical dashed lines indicate the lowest vaccination rates for which vaccinated individuals account for the majority of infections, hospitalizations, and transmission as annotated. \({R}_{0}^{{{\mbox{NPI}}}}\,=4\) for all plots, with baseline VE and immunity parameters vs. the delta variant (Materials and Methods, Supplementary Tables 1 and 2); no screening. bd: 58% vaccination rate and 35% rate of prior infection.

Traditionally, the basic reproductive number R0 is defined as the number of secondary infections generated by a typical infector in an entirely susceptible population, i.e., a population without any NPIs. At the time of writing, NPIs such as masking, ventilation and physical distancing are commonplace in many areas, so we hereafter define \({R}_{0}^{{{\mbox{NPI}}}}\) to be the expected number of secondary infections generated by a typical infector in a population with possible NPIs, but prior to any impacts of population immunity. Furthermore, because precise estimates of \({R}_{0}^{{{\mbox{NPI}}}}\) vary by context, variant, and over time, we consider a range of \({R}_{0}^{{{\mbox{NPI}}}}\) values from 4 to 6. In our baseline modeling scenario, vaccines were assumed to reduce susceptibility to infection by VES = 65%, the likelihood of transmission to others by VEI = 35%, and the likelihood of disease progression to hospitalization conditioned on infection by VEP = 86%, values which land within plausible literature estimates for the effectiveness of two doses of mRNA vaccine against the delta variant in the absence of dramatic waning and without boosting14,18,21,22. Though less often studied in the literature, we assumed that prior SARS-CoV-2 infection would lead to 63% and 13% decreases in risks of infection and transmission, respectively, based on a statistical model relating immunity to neutralization18, and that hybrid immunity would be superior to either vaccination or prior SARS-CoV-2 infection alone. We further assumed an additional 54% decrease in risk of hospitalization conditioned on infection for individuals with prior SARS-CoV-2 infection23, and that individuals with hybrid immunity receive the greater of vaccinated or prior immunity’s protection against hospitalization. See Materials and Methods and Supplementary Tables 1 and 2 for a complete description of the model and parameters.

In a modeled population of N = 20,000 with 58% vaccination rate (corresponding to U.S. estimates as of Nov. 4, 202124) and 35% past infection rate, outbreaks still occurred, despite assuming a partially mitigated delta variant (\({R}_{0}^{{{\mbox{NPI}}}}\,=4\)). During the ensuing outbreak, 59% of total infections and 89% of hospitalizations occurred in unvaccinated individuals (Fig. 1b, c), despite making up only 42% of the population. Furthermore, the peak burden of disease occurred first in the unvaccinated community and then 1 week later in the vaccinated community (Fig. 1b), a known consequence of disease dynamics in populations with heterogeneous susceptibility and transmissibility25,26. By categorizing transmission events into four distinct modes (Fig. 1a), we observe that infections during a delta outbreak in both communities were driven predominantly and consistently by the unvaccinated community (U → U, U → V; Fig. 1d), but that there was nevertheless some transmission from the vaccinated community (breakthrough transmission). These differences occurred despite a “well-mixed” modeling assumption—namely, that an individual with a given vaccination status is no more or less likely to associate with a member of their own group vs. the other group.

Vaccination and past infection rates vary widely across the U.S.24 and the world27 due to impacts of both vaccine availability27 and refusal28, as well as the success or failure of transmission mitigation policies. We therefore asked how a population’s vaccination and past infection rates would affect our observations about infections, hospitalizations, and the relative impacts of the four modes of transmission. This analysis revealed three important points.

First, our results reinforce the fact that increased vaccination rates lead to decreased total infections and hospitalizations, both before and after the herd immunity threshold at Reff = 1 (Fig. 1e, f). Moreover, when large proportions of the population are also partially protected by immunity from prior infection, the vaccination levels required to reach Reff = 1 decrease considerably (Fig. 2a). For instance, increasing prior infection rates from 35 to 50% decreases the required vaccination rate for Reff = 1 from 87 to 80% under baseline modeling parameters. Combinations of immunity from past infection and vaccination thus have the potential to create a herd immunity frontier, beyond which transmission is no longer self-sustaining even in the absence of screening. We caution that although total infections and hospitalizations may appear equal along lines of constant Reff (Fig. 2a, b), both actually decrease as vaccination rates increase, due to vaccines’ superior protection, relative to prior infection, against infection and hospitalization.

Fig. 2: Vaccination and prior infection rates affect epidemic potential, vaccine breakthroughs, and drivers of transmission.
figure 2

Heatmaps show a the total number of infections, b the total number of hospitalizations, c the percentage of total infections occurring in the unvaccinated population, d the percentage of total infections caused by the unvaccinated population, and e the percentage of total hospitalizations occurring in the unvaccinated population for simulated epidemics (see text). White annotation curves show (a, b) isoclines of the effective reproductive number Reff calculated at t = 0, and the line of parameters along which (c) 50% of infections were breakthroughs, d 50% of transmission was due to breakthrough infections, and e 50% of hospitalizations were breakthroughs. N = 20,000 and \({R}_{0}^{{{\mbox{NPI}}}}\,=4\) for all plots, with baseline VE and immunity parameters vs. the delta variant (Materials and Methods, Supplementary Tables 1 and 2); no testing. See Supplementary Fig. 2 for \({R}_{0}^{{{\mbox{NPI}}}}\,=6\).

Second, as vaccination rates increased, the fraction of infections classified as vaccine breakthroughs increased (Fig. 1e), creating a transition point such that when 68% of the population was vaccinated, 50% of all infections were breakthrough infections under our baseline modeling conditions for the delta variant. To determine whether this transition point of 68% was sensitive to the precise fraction of the population with immunity from past infection (35%, Fig. 1), we varied the fraction with infection-acquired immunity between 0 and 100%, finding that the 50/50 breakthrough infection transition occurred between 63 and 75% vaccine coverage (Fig. 2c). Because vaccines provide an additional level of protection against hospitalization VEP, the 50/50 breakthrough hospitalization transition occurs at rates of vaccination of 90–96% (Figs. 1f and 2e). Thus, our results set the expectation that increasing vaccination rates will decrease total infections and hospitalizations, yet a higher proportion of both will be breakthroughs, irrespective of levels of immunity due to prior infection.

Third, as vaccination rates increased, the unvaccinated community ceased to be the primary driver of transmission. Under our baseline modeling conditions (\({R}_{0}^{{{\mbox{NPI}}}}\,=4\), 35% with infection-acquired immunity), this transition occurred when 80% or more of the population was vaccinated (Fig. 1g). When we varied the fraction of the population with infection-acquired immunity between 0 and 100%, this transition point varied from 76 to 82% (Fig. 2d). Thus, while COVID-19 hospitalizations remain concentrated primarily in unvaccinated populations, only a minority of infections will occur in, and be driven by, the unvaccinated community when vaccine coverage is sufficiently high. Note that this implies that unvaccinated individuals living in highly vaccinated communities will still be exposed to SARS-CoV-2 and thus remain at risk of infection and severe disease.

These findings are driven by reductions in susceptibility, disease severity, and infectiousness arising from vaccination, prior SARS-CoV-2 infection, or both. However, quantitative estimates of those reductions vary depending on which vaccine was administered17, time since vaccination or SARS-CoV-2 infection15,16,17, whether an additional “booster” dose was given18, and the variant circulating at the time of the study29,30. We therefore sought to determine how our findings might change under different sets of assumptions about vaccine effectiveness by comparing our baseline scenario (VES = 0.65, VEI = 0.35, VEP = 0.86) with a waning/low immunity scenario (VES = 0.5, VEI = 0.1, VEP = 0.80) and a boosted/high immunity scenario (VES = 0.8, VEI = 0.6, VEP = 0.90), as well as a scenario reflecting plausible VE values based on early observations for the omicron variant (VES = 0.35, VEI = 0.05, VEP = 0.7722,31). To explore the impact of these changes in vaccine effectiveness, we simulated outbreaks for all combinations of vaccination and infection-acquired immunity rates under the four VE scenarios. Across simulations, total infections and hospitalizations were well predicted by calculating Reff at the start of each simulation (Eq. (3); Methods). In particular, outbreaks were small when vaccination or past infection rates crossed the herd immunity threshold (Reff < 1). When Reff > 1, total infections monotonically increased as Reff increased (Supplementary Fig. 1). The herd immunity threshold was impossible to cross with vaccination alone in the waning or omicron VE scenarios with partially mitigated transmission (\({R}_{0}^{{{\mbox{NPI}}}}\,=4\), Fig. 3a, d and Supplementary Fig. 1), and in waning, baseline, and omicron VE scenarios with unmitigated transmission (\({R}_{0}^{{{\mbox{NPI}}}}\,=6\); Supplementary Fig. 1), as evidenced by the fact that the Reff = 1 curves either fail to intersect the vaccination rate axis or appear at all.

Fig. 3: Transition points for breakthrough infections, hospitalizations, and transmission.
figure 3

All panels show curves representing the vaccination and prior infection rates infections (pink), transmission (green), and hospitalizations (purple) are composed of equal numbers of vaccinated an unvaccinated individuals, with majority-breakthrough regions to the right of each line as indicated, for a waning/low, b baseline, and c boosted/high VE vs. the delta variant, and d plausible VE vs. the omicron variant. Black lines indicate Reff = 1 isoclines, which do not appear in d due to Reff > 1. See Supplementary Table 2 for immunity parameter values. \({R}_{0}^{{{\mbox{NPI}}}}\,=4\) in all panels.

Waning, boosting, or omicron-specific assumptions altered the proportions of infections and hospitalizations occurring in, and transmission from, the unvaccinated vs. vaccinated communities. All else being equal, waning or omicron VE led to increased fractions of breakthrough infections and hospitalizations, and increased transmission from the vaccinated community, while boosted VE led to decreases of all three. In turn, the population vaccination rates at which the majority of infections or hospitalizations were breakthroughs shifted down for waning or omicron VE (Fig. 3a, d), while the vaccination rate at which the majority of transmission was driven by vaccinated individuals shifted up for boosted VE (Fig. 3c).

Among the four transition points identified in transmission dynamics, we observe that, in each VE scenario, Reff is driven by both vaccination and past infection rates, as evidenced by curvature in Reff = 1 isoclines (Fig. 3, black lines). In contrast, isoclines representing the transition points between majority-unvaccinated vs. majority-breakthrough infections (Fig. 3, pink lines), between majority-unvaccinated vs. majority-breakthrough hospitalizations (Fig. 3, purple lines), and between majority-unvaccinated vs. majority-breakthrough transmission (Fig. 3, green lines) are relatively insensitive to variation in rates of infection-acquired immunity, as evidenced by vertical or near-vertical isoclines. These findings suggest that the relative proportions of breakthrough infections, hospitalizations, and transmission are driven more by vaccination rates and VE, but not by rates of past infection or proximity to herd immunity; indeed, after the herd immunity threshold, all three isoclines show essentially no variation. These observations suggest that unvaccinated-only screening programs, which decrease rates of U → U and U → V transmission, may be highly effective only in regimes where transmission is driven by the unvaccinated (i.e., to the “left” of green isoclines, Fig. 3), an intuition we now explore in detail.

The impacts of unvaccinated-only screening depend on population immunity, compliance, and VE

To explore the impact of unvaccinated-only screening on population transmission, we modified our simulations so that a positive test would result in an unvaccinated individual isolating to avoid infecting others7,32. We considered test sensitivity equivalent to RT-PCR with a 1-day delay between sample collection and diagnosis under three screening paradigms: weekly testing with 50% compliance—a value which reflects observed compliance with a weekly testing mandate in a university setting5—weekly testing with 99% compliance, and, specifically for the omicron variant, twice-weekly testing with 99% compliance.

Our analysis shows that the benefits of an unvaccinated-only screening program fall into one of three categories, depending on the population vaccination rate and transmission dynamics. These categories align with three distinct regions in parameter space, denoted in Fig. 4 as regions I, II and III. In region I, screening is insufficient to fully control transmission, yet nevertheless markedly reduces the peak number of total infections, colloquially “flattening the curve” (Fig. 4a). In region II, screening successfully brings transmission under control (Fig. 4b). In region III, screening has little impact on transmission due to the fact that outbreaks are already mitigated by population immunity and other control measures (Fig. 4c). Unvaccinated-only screening is therefore impactful in the first two regions, sufficient for transmission control in only the second region, and largely inconsequential to transmission in the third.

Fig. 4: The impact of unvaccinated-only screening corresponds to three distinct parameter regions.
figure 4

Total number of infections with no screening (black) and weekly testing with 99% compliance (pink) are shown for a 25%, b 75%, and c 95% population vaccination rates. d Effective reproductive number over various population vaccination rates, where Reff = 1 is denoted by gray dashed line. The impacts of screening fall into three categories (see text) depending on whether vaccination rate falls into region I, II, or III, as annotated. \({R}_{0}^{{{\mbox{NPI}}}}\,=4\) and 35% rate of prior infection with baseline immunity parameters (Materials and Methods, Supplementary Table 2).

The three regions that correspond to different impacts of screening on transmission are separated by boundaries which can be estimated from two analytical calculations of Reff—one which includes the effects of screening and one which does not (Eq. (3), Methods). The boundary separating regions I and II is given by those parameters for which Reff = 1 with screening, while the boundary separating regions II and III is given by those parameters for which Reff = 1 without screening (Fig. 4d). Thus, the value of a screening testing program in reducing infections can be evaluated based on which of three regions the current vaccination rate, prior infection rate, and VE fall into.

To illustrate the value of this Reff-based analysis, we considered vaccination rates and prior infection rates ranging from 0–100% and varied VE between waning, baseline, and boosted scenarios for the delta variant. Across scenarios, dramatic relative reductions in cumulative infections are concentrated within the envelope between the boundaries of Reff = 1 with and without screening, i.e., region II (Fig. 5). Outside of this effective screening envelope, percent reductions in cumulative infections decreased markedly, either because unvaccinated-only screening flattened the infection curve but had little impact on cumulative infections (region I), or because existing population immunity prevented large outbreaks in the first place (region III). Assuming a 35% past infection rate and \({R}_{0}^{{{\mbox{NPI}}}}\,=4\), region III appeared only for baseline and boosted vaccine effectiveness assumptions, and only when vaccination rates were ~90% or greater (baseline VE) or 75% or greater (boosted VE). Sensitivity analyses show that increasing \({R}_{0}^{{{\mbox{NPI}}}}\) to 6, potentially representing pre-pandemic contact rates and the SARS-CoV-2 delta variant, cause region III to shrink further (Supplementary Fig. 5). Thus continued screening for SARS-CoV-2 among the unvaccinated may be of limited value when vaccination rates are sufficiently high.

Fig. 5: The impacts of unvaccinated-only screening depend on population immunity, compliance, and vaccine effectiveness.
figure 5

Percent reduction in cumulative infections due to screening over various population vaccination rates assuming waning/low (a, d), baseline (b, e), and boosted/high (c, f) VE with once-weekly screening at 50% (top row) and 99% (bottom row) compliance. White lines indicate the population immunity rate at which Reff = 1 with screening (solid) and without screening (dashed), which divide the space into three regions, labeled I, II and III. See Supplementary Tables 1 and 2 for parameter values. \({R}_{0}^{{{\mbox{NPI}}}}\,=4\) in all panels; see Fig. 3 for \({R}_{0}^{{{\mbox{NPI}}}}\,=6\).

The role of compliance—the fraction of scheduled tests that are actually taken—can also be clarified by examining the three regions of screening testing impact. Both the simulations and equations for Reff show that increasing compliance from 50% (Fig. 5, row 1) to 99% (Fig. 5, row 2) causes the lower boundary of the effective screening envelope to shift to lower vaccination and prior infection rates, decreasing the size of region I and increasing the size of region II. Moreover, increased compliance increases the magnitude of infection reductions within both regions, visible as an intensification of color in the infection reduction heatmaps (Fig. 5). As a result of these observations, we conclude that, in addition to test sensitivity, frequency, and turnaround time7, high participation in screening programs is critical to expanding the impact of unvaccinated-only screening testing programs. However, we also note that compliance had little effect in region III where Reff < 1, a result which parallels analysis of universal screening programs9.

This Reff-based analysis of transmission is not restricted to unvaccinated-only screening programs. To illustrate this, we considered an identical set of simulations as in Fig. 5 but with universal screening, i.e., screening via testing of the vaccinated and unvaccinated populations alike. Universal screening caused the boundary between regions I and II (Reff = 1 with screening) to shift, expanding the size of the effective screening envelope (Supplementary Fig. 6). While we present these results here for completeness, universal testing in mixed vaccination status populations have been investigated elsewhere prior to the present work9.

The impact of screening on hospitalizations is also predicted well by the Reff-based effective screening envelope. While hospitalizations were not identical across all equal-Reff combinations of vaccination and prior infection rates, dramatic relative reductions in cumulative hospitalizations were nevertheless clearly concentrated within region II, with decreasing relative reductions in regions I and III (Supplementary Fig. 6). We therefore find that analysis based only on the effective screening envelope, calculated via Eq. (3) (Methods), is useful in predicting the impact of screening regimens—both unvaccinated-only and universal—on reductions in cumulative infections and hospitalizations alike.

The omicron variant’s rapid spread, and in particular its increased rates of reinfection and vaccine breakthrough, raise key questions about the role of unvaccinated-only screening programs and whether populations considering such policies may fall into region I, II, or III defined above. Because estimates of omicron-specific immunity parameters remain either limited or nonexistent at the time of writing, we assumed plausible lower values of each based on extant data22,31, alongside lower infection hospitalization rates for omicron in general (Supplementary Tables 1 and 2). Under such assumptions, weekly unvaccinated-only screening with 50% compliance was ineffective at reducing cumulative infections even though screening reduced the peak magnitude of infections (Fig. 6d). Regardless of compliance, all prior infection and vaccination rate combinations with a weekly screening policy were in region I, indicating that the magnitude of peak infections can be mitigated but the impact on cumulative infections is low (Fig. 6a, b). Doubling the frequency of screening to twice weekly with 99% compliance creates a bifurcated landscape, with highly effective screening only in settings with 18–40% vaccination rates (Fig. 6c). For vaccination rates above 50%, even twice-weekly unvaccinated-only screening programs with near-perfect compliance are unlikely to dramatically impact the spread of the omicron variant (region I). Universal screening showed comparatively higher impact, yet, nevertheless, only twice-weekly testing regimens created broad region II regimes in which community spread was controlled (Supplementary Fig. 6).

Fig. 6: Unvaccinated-only screening during omicron transmission cannot achieve Reff < 1 except in low-vaccination and high-frequency regimes.
figure 6

Percent reduction in cumulative infections due to screening over various population vaccination rates assuming plausible parameters for immunity vs. the omicron variant, with a once-weekly screening at 50% compliance, b once-weekly screening at 99% compliance, and c twice-weekly screening at 99% compliance. d Number of individuals infected over time, under screening scenarios denoted A, B, C, compared with no screening (black) with 58% vaccination rate and 35% rate of prior infection. Solid white line indicates Reff = 1 with screening; Reff = 1 is not achievable without screening. See Supplementary Tables 1 and 2 for parameter values. \({R}_{0}^{{{\mbox{NPI}}}}\,=4\) in all panels; see Supplementary Fig. 4 for \({R}_{0}^{{{\mbox{NPI}}}}\,=6\) and Supplementary Fig. 5 for universal testing.

Unvaccinated-only screening shifts the balance of unvaccinated vs. breakthrough transmission but not infection or hospitalization

By reducing transmission from unvaccinated individuals, screening programs specifically mitigate U → U and U → V transmission modes, thus diminishing the role of the unvaccinated population in transmission dynamics and amplifying the relative role of vaccine breakthrough transmission. As a consequence, we observe that in the presence of screening, the vaccination rates at which the unvaccinated cease to drive a majority of transmission decrease by up to 15 percentage points (Fig. 7b), with the largest decreases for 99% compliance and waning VE vs. delta, and the smallest decreases for 50% compliance and boosted VE vs. delta, or in all screening scenarios vs. the omicron variant. Under waning/low, baseline and omicron VE scenarios, unvaccinated-only screening programs shrink the regime in which the unvaccinated population drives outbreaks.

Fig. 7: Screening and vaccine effectiveness affect transition points to majority-breakthrough regimes.
figure 7

The vaccination rates at which the vaccinated population makes up the majority of a infections, b transmission, and c hospitalizations for low, moderate, and high vaccine effectiveness scenarios. Minimum (filled circle) and maximum (open circle) endpoints show the variation in transition points over all combinations of vaccination and prior infection rates for no screening (black), 50% compliance (purple), 99% compliance (pink) over all possible values for past infection rates. \({R}_{0}^{{{\mbox{NPI}}}}\,=4\) for all plots; see Supplementary Fig. 7 for \({R}_{0}^{{{\mbox{NPI}}}}\,=6\).

In contrast, unvaccinated-only screening programs had little effect on the percentage of infections or hospitalizations that were vaccine breakthroughs. Instead, majority-breakthrough regimes remained primarily dependent on vaccination rates and vaccine effectiveness (Fig. 7a, c), with transitions to majority-breakthrough infection regimes beginning at 55 to 67% vaccination rates (waning VE, delta), 63 to 75% vaccination rates (baseline VE, delta), 83 to 84% vaccination rates (boosted VE, delta), and 50 to 55% vaccination rates (omicron). Transitions to majority-breakthrough hospitalizations occurred at 83 to 88% (waning VE, delta), 91 to 96% (baseline VE, delta), 96 to 99% vaccination rates (boosted VE, delta), and 58 to 83% vaccination rates (omicron), regardless of screening. We therefore conclude that unvaccinated-only screening programs do not markedly alter the expectations of majority-breakthrough infections or hospitalizations at high vaccination levels, particularly if VE is low or waning.

Discussion

In this analysis, we find that in communities with mixed vaccination status, routine SARS-CoV-2 screening programs focused only on the unvaccinated may reduce infections and hospitalizations, but in a manner dependent on two conditions. First, effective screening via testing requires high participation to be most impactful, reinforcing the need for mechanisms to encourage or enforce high participation. Second, when immunity from vaccination and past infection are high enough to curtail transmission on their own, or in concert with effective NPIs12,13, testing the remaining unvaccinated population averts few infections and hospitalizations in both relative and absolute terms. On the other hand, when transmission due to reinfection and/or vaccine breakthrough is sufficiently high, unvaccinated-only screening will at best “flatten the curve” of infections, but cannot adequately control infections and hospitalizations except when testing twice weekly with near-perfect participation in low-vaccination communities. Once communities reach vaccination rates of ~40% or more, even twice-weekly unvaccinated-only screening, with near-perfect compliance and isolation adherence, cannot control the omicron variant. Thus, targeted unvaccinated screening programs are highly effective only in a restricted region of epidemiological parameter space, results echoed by similar work analyzing universal screening programs9.

Key to understanding our study are three observations and findings. First, an unvaccinated-only screening program simply tests fewer and fewer individuals as vaccination rates increase. Second, the relative role of the unvaccinated population in driving transmission decreases as vaccination rates increase, regardless of vaccine effectiveness. Third, when vaccine effectiveness against infection and transmission wanes, unvaccinated-only screening programs decrease in impact. As a consequence, our work broadly suggests that unvaccinated-only screening is most beneficial when a population is undervaccinated and is close to, but has not yet achieved, herd immunity—region II in our analyses—leading to the recommendation that such testing programs be used in conjunction with other NPIs. Indeed, our work finds that unvaccinated-only screening alone is generally insufficient to markedly reduce infections and hospitalizations when (1) the population is far from the herd immunity threshold, inclusive of existing NPIs, in either direction, (2) vaccination rates are high, or (3) testing is weekly and/or compliance is low. Against a backdrop of waning immunity and continued emergence of antigenically distinct variants, even herd immunity may be, at best, temporary. In this context, weekly unvaccinated-only screening programs alone are an insufficient countermeasure for the omicron variant.

Indeed, while our analysis focused on a single screening-based intervention in isolation, unvaccinated-only (or universal) screening programs are typically implemented alongside other NPIs5,6,8. These NPIs, including limitations on gatherings, increasing the availability of personal protective equipment, and school or restaurant closures, were estimated to have reduced the effective reproductive number Reff by 0.1–0.2 in 2020, particularly when implemented early12, and by around 10% in 202113. In comparison, an unvaccinated-only weekly screening policy with a realistic (50%) compliance rate5 and 58% vaccination rate, would reduce Reff by an estimated 10.2% (Eq. (3)), decreasing further as vaccination rates increase, and compliance and isolation adherence decrease. Thus, while our analysis ranks vaccinate-or-test policies as potentially competitive with high-impact NPIs12,13, such screening will decrease in impact as vaccination rates inch higher. Because prior work has shown that pandemic countermeasures also vary in their impact depending on time, vaccination, and the presence of other NPIs or behaviors13, an empirical assessment of vaccinate-or-test programs would be valuable. However, just as empirical estimates of NPIs’ impacts on Reff include wide uncertainty12,13, similar estimates for unvaccinated-only screening programs are also likely to be highly uncertain.

Our study elucidates three critical transitions as vaccination rates increase. First, when vaccination rates are sufficiently high, a majority of the albeit reduced number of infections will be vaccine breakthrough infections. This fact should come as no surprise, as this transition must occur at some point for any vaccine below 100% effectiveness; for the delta variant, our modeling estimates it to take place between 63 and 75% vaccine coverage (baseline VE; 55–67% vaccinated with waning VE; 83–84% vaccinated with boosted VE), while for the omicron variant, we estimate it to take place between 55 and 59% vaccine coverage. Second, a transition to majority-breakthrough hospitalizations will occur at some point greater than the transition to majority-breakthrough infections, a natural consequence of VEP > 0. Third, while community spread is driven by the unvaccinated at low-vaccination rates, it is driven by the vaccinated population at high vaccination rates (Fig. 3). These vaccination rate transition points separating majority-unvaccinated transmission and majority-breakthrough transmission are driven lower by unvaccinated-only screening programs (Fig. 7). Taken together, these results suggest that while the overall number of infections during an outbreak decreases as vaccination rates increase, assuming VES > 0, vaccine breakthrough infections and transmission events from vaccinated individuals should not be surprising in highly vaccinated populations—vaccine effectiveness is imperfect. Consequently, in anticipation of continued community transmission even in highly vaccinated communities, those at increased risk of severe COVID-19 should take additional precautions to limit their risk of infection or severe disease.

Our analyses identify two limitations of screening via testing programs in reducing community transmission which generalize beyond the specific scenarios investigated herein. First, the ability of a screening program to prevent community spread is restricted to a limited “envelope” of past infection rate and vaccination rate combinations such that Reff without screening is greater than one, and Reff with screening is less than one. Second, the width of that effective screening envelope, and thus the effectiveness of a screening program to control transmission more broadly, are highly sensitive to compliance and participation: weekly screening with 50% compliance—a rate which reflects observed compliance in a population with a weekly screening mandate5—is likely to be relatively ineffective. However, although our analyses focus on the benefits of testing in reducing transmission, testing also plays an important role in diagnosis and treatment, detection of variants, situational awareness and surveillance, and decreasing pressure on the healthcare system during outbreaks. Furthermore, testing focused on the unvaccinated population may provide additional incentives to get vaccinated and thus avoid regular testing. Our study did not explore the benefits of unvaccinated-only testing mandates for these additional purposes.

Our analysis is limited in at least three different manners. First, our modeling incorporated fixed parameters that are difficult to estimate in practice. For instance, while our analysis considered boosted, baseline, and waning scenarios for vaccines’ reductions in susceptibility VES, infectiousness VEI, and hospitalization given infection VEP based on ranges of estimates in the current literature, few studies are available to guide estimates of similar risk reductions associated with prior SARS-CoV-2 infection, with or without vaccination (but see refs. 14,18). Moreover, real-time estimates for emerging variants of concern such as omicron require observational study and are thus unavailable for prospective policy analyses. Alternative parameter assumptions may be explored via the provided open source code. Second, we assumed perfect isolation after receiving a positive test result. Were this assumption to be violated by imperfect or delayed isolation, we predict a proportional loss of screening impact across all scenarios. Third, our model assumes values of \({R}_{0}^{{{\mbox{NPI}}}}\) and immunity associated with the delta variant and plausible values for the omicron variant, but other emerging variants may dramatically shift the values of these parameters. These limitations affect the exact vaccination and past infection rates at which the three transitions identified in our study occur, and thus our analyses describe fundamental phenomena but do not make projections or predictions for specific communities.

Our analysis also uses a well-mixed SEIR model framework, inheriting two key limitations which merit direct discussion. First, we assume that vaccination and past infection statuses are uncorrelated at the population level. In reality, they may be anticorrelated due to the protective effects of vaccination, or because those with past infection may choose to forgo subsequent vaccination. We similarly assumed no homophily in contact patterns based on vaccination status, following the well-mixed assumption of the SEIR model framework, yet those who choose to be vaccinated may be more likely to be situated in a social network with others who choose to be vaccinated, and vice versa33. Second, compartmental SEIR models such as ours assume uniform infectiousness in the I compartment, contrasting empirical observations34 and more sophisticated models7,32. While our model’s latent and infectious periods are well aligned with other SEIR models9,35,36,37, they nevertheless lead to unrealistically long generation times. Decreasing these periods proportionally to achieve the same reproductive number while aligning more closely with generation time estimates38 would change the time-scale across all simulations, but would not impact the cumulative metrics or dynamics discussed in our key results.

More broadly, our work is situated within a family of research which uses mathematical modeling to estimate the impact of targeted countermeasures or strategies in populations with heterogeneous susceptibility, transmissibility, and/or contact rates. Other areas of focus include the allocation of scarce personal protective equipment to reduce transmission39, the prioritization of vaccines by subpopulation40,41,42, proactive screening programs in specific workplace structures43 or contact networks9, immunity “passport” programs32, or immune shielding strategies 44. While our analyses are directed at SARS-CoV-2, this work illustrates contributes general principles to this literature by showing that screening programs focused on testing the unvaccinated may be less effective than hoped in the face of high vaccination rates, waning vaccine effectiveness, or low compliance with testing.

Methods

SEIR model

Our analyses are based on a continuous time ordinary differential equation compartmental model with susceptible, exposed, infectious, and recovered (SEIR) compartments, stratified into vaccinated V and unvaccinated U groups. In addition to tracking infections among these two groups separately, we also tracked infections from both groups separately, enabling us to investigate four modes of transmission: from U to U, from U to V, from V to U, and from V to V. In all simulations, we used a constant total population size of N = 20,000 and denoted the vaccinated fraction of the population with ϕ.

To incorporate the possibility that individuals may have experienced prior infections, we further subdivided U and V into SARS-CoV-2 naive and SARS-CoV-2 experienced subpopulations, such that a fraction ψ of each was assumed to be previously infected and 1 − ψ remains naive. For notation, we denote the subpopulations of U to be u (unvaccinated, naive) and x (unvaccinated, experienced/prior infection), and the subpopulations of V to be v (vaccinated, naive) and h (vaccinated, experienced). We assumed that vaccination and SARS-CoV-2 experience statuses were fixed at the start of each simulation and immutable throughout, such that there was no ongoing vaccination, and individuals who were infected and recovered during each simulation were not reassigned to SARS-CoV-2 experienced status45.

We denote the protective effects of immunity as XE, VE, HE, expressed as reductions in risk due to prior infection alone (x), vaccination alone (v), or prior infection and vaccination (i.e., so-called “hybrid” immunity; h), respectively. Immunity was modeled to (1) decrease the risk of infection upon exposure, (2) decrease the risk of transmission upon infection, placing our vaccine and immunity model in the broader category of leaky models46, and (3) decrease the risk of disease progression (i.e., hospitalization) upon infection. Reductions in the risk of infection upon exposure (XES, VES, HES), reductions in the risk of transmission when infected (XEI, VEI, HEI), and reductions in the risk of hospitalization when infected (XEP, VEP, HEP) were parameterized separately, based on ranges of estimates from the literature. Note that VEH, the reduction in risk of hospitalization due to vaccination, is more commonly reported in the literature than VEP, the reduction in risk of hospitalization due to vaccination conditional on infection. So, we used the formula \({{{\mbox{VE}}}}_{P}=1-\frac{1-{{{\mbox{VE}}}}_{H}}{1-{{{\mbox{VE}}}}_{S}}\) to estimate values for VEP. We used the same relationship to estimate XEP. See Supplementary Table 2. Due to broad uncertainty in these effects over time since exposure14,45 or vaccination14,15,16, by vaccine manufacturer and schedule17,18,47,48, by context30,49, and by variant18, our analyses intentionally consider a range of values. We assumed that hybrid immunity against infection, HES, and transmission, HEI, would always be superior to either vaccination alone or prior infection alone, via the simple formula HE = (1 − XE)VE + XE and hybrid immunity against hospitalization given infection HEP = max{VEP, XEP}.

Supplementary Fig. 8 shows a model schematic diagram for the SEIR model used in the manuscript, where solid and dashed lines denote movement and transmission between classes, respectively. In lieu of including hospitalization as a model compartment, we computed total hospitalizations in each immunity class via multiplication of total infections, variant-specific infection hospitalization rates and (1 − RRP), where RRP is the risk reduction against progression to hospitalization given infection (i.e., 0, VEP, XEP, or HEP) (See Supplementary Tables 1 and 2).

To model a community with open boundaries, we included a uniform risk of exposure to infection from an external source at a rate of N−1 per person per day. For instance, in a completely naive population, Su/N individuals would be infected per day. After including the protective effects of vaccination and past infection this resulted in importation of infections at per-capita rates of (1 − VES)N−1, (1 − HES)N−1, (1 − XES)N−1, and N−1 new infections per day in the v, h, x, and u groups respectively.

All simulations were run for 270 days, and all individuals were initially in one of the susceptible compartments Su, Sx, Sv, or Sh in proportions (1 − ϕ)(1 − ψ), (1 − ϕ)ψ, ϕ(1 − ψ), and ϕψ, respectively. Model equations were solved using lsoda solver from the package deSolve, R version 4.1.0.

Incorporation of community testing

Screening the unvaccinated population via community testing, and subsequent isolation of those testing positive, was modeled by increasing the rate at which infected individuals were removed from the unvaccinated Iu and Ix compartments. Similarly, universal screening regardless of vaccination status was modeled by increasing the rate at which infected individuals were removed from all I compartments, Iu, Ix, Iv and Ih. The effectiveness of screening tests was assumed to be equal for vaccinated and unvaccinated individuals. We estimated increased rates of removal taking into account (1) the calibrated trajectories of viral loads within individual infection50, (2) the relationship between viral load and infectiousness7, (3) the frequency of testing, (4) the test’s analytical sensitivity (i.e., limit of detection) and turnaround time32, and (5) screening compliance and valid sample rates, i.e., the fraction of scheduled or mandated tests which actually produce a valid sample5. In particular, our adaptation takes a previous model7,32 and updates viral load dynamics for the delta variant of SARS-CoV-251,52, the dominant variant at the time of the present analysis. To incorporate the effectiveness of screening θ, we reduce the duration of infectiousness 1/γ by a factor (1 − θ). Parameter values for θ are found in Supplementary Table 1, and are based on weekly PCR testing with a 1-day turnaround, analytical limit of detection of 103 RNA copies per ml sample, and compliance rates of 50% (as in ref. 5) or 99% (as in ref. 8). These values assume that individuals immediately and successfully isolate upon receiving a positive diagnosis. We note that estimated effects of rapid antigen tests (with higher analytical limits of detection, but zero turnaround time) are highly similar to PCR testing under the assumptions above, provided that the community testing program frequencies and compliance rates are identical7.

Transmission modes and forces of infection

Inclusive of all effects introduced above, the forces of infection are given by:

$${\lambda }_{u}=\, \alpha \Big(\frac{{I}_{u}}{{N}_{u}}{c}_{u\to u}+[1-{{{\mbox{XE}}}}_{I}]\frac{{I}_{x}}{{N}_{x}}{c}_{x\to u}+[1-{{{\mbox{VE}}}}_{I}]\frac{{I}_{v}}{{N}_{v}}{c}_{v\to u}\\ +[1-{{{\mbox{HE}}}}_{I}]\frac{{I}_{h}}{{N}_{h}}{c}_{h\to u}\Big)+\frac{1}{N}$$
(1)
$${\lambda }_{i}=\, \Big[\alpha \Big(\frac{{I}_{u}}{{N}_{u}}{c}_{u\to i}+[1-{{{\mbox{XE}}}}_{I}]\frac{{I}_{x}}{{N}_{x}}{c}_{x\to i}+[1-{{{\mbox{VE}}}}_{I}]\frac{{I}_{v}}{{N}_{v}}{c}_{v\to i}\\ + [1-{{{\mbox{HE}}}}_{I}]\frac{{I}_{h}}{{N}_{h}}{c}_{h\to i}\Big)+\frac{1}{N}\Big][1-{({{{{{\rm{R}}}}}}{{{{{{\rm{R}}}}}}}_{S})}_{i}]\,,$$
(2)

where \(i=\left\{x,v,h\right\}\), and reductions in susceptibility due to immunity are given by \({({{{{{\rm{R}}}}}}{{{{{{\rm{R}}}}}}}_{S})}_{i}=\left\{{{{\mbox{XE}}}}_{S},{{{\mbox{VE}}}}_{S},{{{\mbox{HE}}}}_{S}\right\}\), correspondingly. The parameter α is the probability an unvaccinated, SARS-CoV-2 naive individual is infected by an infectious contact, tuned to achieve the desired \({R}_{0}^{{{\mbox{NPI}}}}\), ci → j is the number of times an individual in group j is contacted by individuals from group i per day, and Nj is a convenience variable representing the number of people in subpopulation j.

To produce counts of how many infections were caused by each of the transmission modes U → U, U → V, V → U, and V → V, we integrated the appropriate terms of Eqs. (1) and (2) over the duration of each simulation. For instance, the cumulative number of vaccinated infections caused by the unvaccinated population is given by integrating over the forces of infection from u and x to v and h:

$$U\to V=\, \alpha \int\nolimits_{0}^{270}\Big[\frac{{I}_{u}(t)}{{N}_{u}}\left({c}_{u\to v}{S}_{v}(t)[1-{{{\mbox{VE}}}}_{S}]+{c}_{u\to h}{S}_{h}(t)[1-{{{\mbox{HE}}}}_{S}]\right)\,\cdots \\ +[1-{{{\mbox{XE}}}}_{I}]\frac{{I}_{x}(t)}{{N}_{x}}\left({c}_{x\to v}{S}_{v}(t)[1-{{{\mbox{VE}}}}_{S}]+{c}_{x\to h}{S}_{h}(t)[1-{{{\mbox{HE}}}}_{S}]\right)\Big]dt$$

Reproductive number

The basic reproductive number R0 is defined as the expected number of secondary infections created by a typical infector in an entirely susceptible population, in the absence of any NPIs. Given the variety of environments in which SARS-CoV-2 spreads, and the presence of various permanent or semi-permanent NPIs, we use \({R}_{0}^{{{\mbox{NPI}}}}\) in this work to denote the reproductive number in an entirely susceptible population inclusive of varying levels of now-baseline NPIs for the delta, omicron, and future variants. We consider \({R}_{0}^{{{\mbox{NPI}}}}\,=4\) (Main and Supplementary Figures) and \({R}_{0}^{{{\mbox{NPI}}}}\,=6\) (Supplementary Figures). For unvaccinated-only screening programs, this model’s effective reproductive number is given by:

$${R}_{{{\mbox{eff}}}}={R}_{0}^{{{\mbox{NPI}}}}\,\left[{f}_{u}(1-\theta )+{f}_{x}{r}_{x}(1-\theta )+{f}_{v}{r}_{v}+{f}_{h}{r}_{h}\right],$$
(3)

where fu = (1 − ψ)(1 − ϕ), fx = ψ(1 − ϕ), fv = (1 − ψ)ϕ, and fh = ϕψ represent the fractions of the population in the unvaccinated, experienced, vaccinated, and hybrid immunity groups, respectively, and rx = (1 − XEI)(1 − XES), rv = (1 − VEI)(1 − VES), and rh = (1 − HEI)(1 − HES) are the cumulative impacts of immunity on each group. Setting the above equation equal to a constant produces isoclines shown in plots throughout the paper. The reduction in Reff due to screening is given by:

$${R}_{{{{{{{{\rm{noscreening}}}}}}}}}-{R}_{{{{{{{{\rm{screening}}}}}}}}}={R}_{0}^{{{\mbox{NPI}}}}\,\theta (1-\phi )\left[1-\psi (1-{r}_{x})\right],$$
(4)

a function linear in each of its variables which goes to zero as the vaccination rate ϕ approaches 1.

For universal screening programs, similar calculations yield:

$${R}_{{{\mbox{eff}}}}^{{{{{{{{\rm{universal}}}}}}}}}={R}_{0}^{{{\mbox{NPI}}}}\,(1-\theta )\left[{f}_{u}+{f}_{x}{r}_{x}+{f}_{v}{r}_{v}+{f}_{h}{r}_{h}\right],$$
(5)

differing from Eq. (3) only in the terms to which (1 − θ) applies. For a complete derivation of these equations, see Supplementary Materials.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.