Next Article in Journal
Intrauterine Fetal Demise After Uncomplicated COVID-19: What Can We Learn from the Case?
Next Article in Special Issue
Modeling the Transmission of the SARS-CoV-2 Delta Variant in a Partially Vaccinated Population
Previous Article in Journal
Prediction of SARS-CoV-2 Variant Lineages Using the S1-Encoding Region Sequence Obtained by PacBio Single-Molecule Real-Time Sequencing
Previous Article in Special Issue
Modelling Mutation in Equine Infectious Anemia Virus Infection Suggests a Path to Viral Clearance with Repeated Vaccination
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Influence of Vaccine Administration on COVID-19 Testing Strategies

by
Jonathan E. Forde
1,*,† and
Stanca M. Ciupe
2,†
1
Department of Mathematics and Computer Sciences, Hobart and William Smith Colleges, Geneva, NY 14456, USA
2
Department of Mathematics, Virginia Tech, Blacksburg, VA 24060, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 30 September 2021 / Revised: 13 December 2021 / Accepted: 15 December 2021 / Published: 19 December 2021
(This article belongs to the Special Issue Mathematical Modeling of Viral Infection)

Abstract

:
Vaccination is considered the best strategy for limiting and eliminating the COVID-19 pandemic. The success of this strategy relies on the rate of vaccine deployment and acceptance across the globe. As these efforts are being conducted, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is continuously mutating, which leads to the emergence of variants with increased transmissibility, virulence, and resistance to vaccines. One important question is whether surveillance testing is still needed in order to limit SARS-CoV-2 transmission in a vaccinated population. In this study, we developed a multi-scale mathematical model of SARS-CoV-2 transmission in a vaccinated population and used it to predict the role of testing in an outbreak with variants of increased transmissibility. We found that, for low transmissibility variants, testing was most effective when vaccination levels were low to moderate and its impact was diminished when vaccination levels were high. For high transmissibility variants, widespread vaccination was necessary in order for testing to have a significant impact on preventing outbreaks, with the impact of testing having maximum effects when focused on the non-vaccinated population.

1. Introduction

The emergence in 2019 of the novel acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Wuhan, China [1] has had devastating global consequences including loss of lives [2], strained healthcare systems [3,4,5], and economic recessions. Protective measures, such as masking, physical distancing, testing, contact tracing, quarantine, and isolation, while effective when applied rigorously, have proved insufficient in limiting the spread of SARS-CoV-2 [6,7,8]. Development of COVID-19 vaccines has been the main focus of public health organizations and pharmaceutical companies across the world [9,10,11].
The U.S. Food and Drug Administration (FDA) has approved two mRNA (Pfizer-BioNTech and Moderna) and one viral vector-based (Janssen) COVID-19 vaccines [12], which have consistently demonstrated effectiveness against disease and increased protection against infection [13,14,15,16]. Despite early control in highly vaccinated communities, vaccine shortage in low-access countries and vaccine hesitancy in high-access countries or regions, has led to selection of new variants [17], which might overcome vaccine-induced immunity [18,19].
Delta variant (B.1.617.2) was identified in the United States in March 2021 and accounted for 98.4 % of new infections by 23 September 2021 [17]. It is believed to be more contagious than the ancestral SARS-CoV-2. One study quantified the mean basic reproductive number for the Delta variant to be R 0 = 5.08 , compared to R 0 = 2.79 for the ancestral strain [20]. Another estimated that the Delta variant was 58– 120 % more transmissible than Alpha variant (B.1.1.7) across New England [21]. Recent studies measuring the effectiveness of vaccines against the transmission of the Delta variant found a reduction in the Pfizer-BioNTech effectiveness to 33.5 % and 88 % after one and two doses, compared to 49 % and 94 % effectiveness against the transmission of the Alpha variant [19]. CDC defines vaccine breakthrough infection as the detection of SARS-CoV-2 RNA or antigen in a respiratory specimen collected from a person ≥14 days after receipt of all recommended doses of an FDA-authorized COVID-19 vaccine [22]. The increased transmissibility of the Delta variant and the observed waning in vaccine effectiveness against infection highlight the critical importance of vaccinating an entire community and of offering boosters to individuals vaccinated at least six months prior [23]. With vaccines only recently approved for children aged 5 to 11, it is also important to adhere to rigorous COVID-19 prevention strategies such as masking and testing. One question of importance to the public health authorities, and the goal of this study, is to determine the best strategy for using testing in an increasingly vaccinated population in order to limit transmission.
SARS-CoV-2 diagnostic and surveillance testing is an important intervention tool for controlling transmission. While still widely used in monitoring the non-vaccinated population, little is known about the potential benefit of testing vaccinated population. Given limited resources, it is important to identify who to test in order to most effectively control transmission at minimal cost.
In this study, we developed a multi-scale mathematical model to examine the effect of testing on blocking SARS-CoV-2 infections in populations vaccinated with the Pfizer-BioNTech vaccine [24,25,26]. We expanded our previous work, which connected a within-host model for the time of SARS-CoV-2 infectiousness onset based on an individual virus dynamics and a between-host model for the transmission at the population level [24]. In particular, we incorporated a variable for the population vaccinated with the Pfizer-BioNTech vaccine and two populations of vaccinated individuals that become infected and have asymptomatic or symptomatic disease. We connected the vaccine effectiveness in preventing infection with the age of vaccination and the effect of the vaccine on reported individual virus profile when a breakthrough infection occurs [27]. We were interested in determining how vaccination level combined with surveillance testing can help reduce an infectious event with variants of various infectiousness level.

2. Materials and Methods

We model the interaction between a susceptible class S ( t ) , vaccinated class v ( t , α ) , infected class of asymptomatic individuals i a ( τ , t ) , infected class of symptomatic individuals i s ( τ , t ) , infected class of vaccinated asymptomatic individuals i v a ( τ , t ) , infected class of vaccinated symptomatic individuals i v s ( τ , t ) , individuals recovered from natural infection R ( t ) , individuals vaccinated after natural infection R v 1 ( t ) and individuals recovered after being vaccinated and then infected R v 2 ( t ) . The independent variables are the age of infection in an individual τ , the age of vaccination in an individual α , and the time-since-outbreak in the population t. The parameters are the transmission rate β , the infection weighting functions λ j , the birth rate b, the death rate μ , the disease-induced mortality rate m, the vaccination rate ν , and the degree of protection after vaccine administration η . Moreover, we consider testing in both vaccinated and non-vaccinated populations, with a testing capacity of C tests per day, leading to case detection rates r j , with j { a , s , v a , v s } .
As in our earlier study [24], we assume that an individual’s infection status is given by its virus profile at age of infection τ . In particular, given virus profiles for infected individuals, we link test positivity to the ages of infection during which virus load is above the sensitivity threshold. Similarly, we determine the ages of infection during which the virus load is high enough to allow transmission. We define
τ 1 j = age   for   onset   of   virus   detectability   by   PCR   test , τ 2 j = age   for   onset   of   infectiousness , τ 3 j = age   for   end   of   infectiousness , τ 4 j = age   for   end   of   virus   detectability   by   PCR   test ,
where j { a , s , v a , v s } . The infectivity weighting functions λ j are
λ a ( τ ) = γ , for τ [ τ 2 a , τ 3 a ] 0 , otherwise , λ s ( τ ) = 1 , for τ [ τ 2 s , τ 3 s ] 0 , otherwise , λ v a ( τ ) = γ , for τ [ τ 2 v a , τ 3 v a ] 0 , otherwise , λ v s ( τ ) = 1 , for τ [ τ 2 v s , τ 3 v s ] 0 , otherwise ,
where parameter 0 < γ < 1 represents the relative infectiousness of asymptomatic carriers, in comparison with symptomatic carriers [22]. For all infected classes, 0 τ τ q . For τ > τ q infections are considered resolved, and recovered individual are not susceptible to reinfection.
Daily testing rate. As before [24], we define a daily per capita testing rate, ρ , corresponding to an overall testing capacity of C tests per day at time t in a population of size N ( t ) to be given by
ρ = ln 1 C N ( t ) ,
with C < N ( t ) . When testing is administered to only non-vaccinated individuals,
N ( t ) = S ( t ) + R ( t ) + 0 τ q i a ( τ , t ) + i s ( τ , t ) d τ ,
and the case detection rate functions r j ( τ , t ) become
r j ( τ , t ) = ρ ( t ) , for t 0 and τ 1 j τ τ 4 j 0 , otherwise ,
for j { a , s } , and r v a ( τ , t ) = r v s ( τ , t ) = 0 for all τ .
When testing is administered to both vaccinated and non-vaccinated individuals,
N ( t ) = S ( t ) + R ( t ) + R v 1 ( t ) + R v 2 ( t ) + 0 τ q i a ( τ , t ) + i s ( τ , t ) + i v a ( τ , t ) + i v s ( τ , t ) d τ ,
and the case detection rate functions r j ( τ , t ) become
r j ( τ , t ) = ρ ( t ) , for t 0 and τ 1 j τ τ 4 j 0 , otherwise ,
where j { a , s , v a , v s } . We assume a test return delay of days and that individuals who receive a positive test result are isolated, and can no longer transmit the virus. Lastly, we ignore the possibility of reinfection.
Model equations. On the domain t 0 , 0 τ τ q , and α 0 , the between-host model equations under combined vaccination and testing are
d S d t = b μ S ν S β S ( t ) Λ ( t ) , v α + v t = μ v ( α , t ) β ( 1 η ( α ) ) v ( α , t ) Λ ( t ) , i a τ + i a t = ( μ + m ) i a ( τ , t ) r a ( τ , t ) i a ( τ , t ) e ( μ + m ) , i s τ + i s t = ( μ + m ) i s ( τ , t ) r s ( τ , t ) i s ( τ , t ) e ( μ + m ) , i v a τ + i v a t = ( μ + m ) i v a ( τ , t ) r v a ( τ , t ) i v a ( τ , t ) e ( μ + m ) , i v s τ + i v s t = ( μ + m ) i v s ( τ , t ) r v s ( τ , t ) i v s ( τ , t ) e ( μ + m ) , d R d t = μ R ν R + i a ( τ q , t ) + i s ( τ q , t ) + 0 τ q r a ( τ , t ) i a ( τ , t ) e ( μ + m ) + r s ( τ , t ) i s ( τ , t ) e ( μ + m ) d τ , d R v 1 d t = μ R v 1 + ν R , d R v 2 d t = μ R v 2 + i v a ( τ q , t ) + i v s ( τ q , t ) + 0 τ q r v a ( τ , t ) i v a ( τ , t ) e ( μ + m ) + r v s ( τ , t ) i v s ( τ , t ) e ( μ + m ) d τ .
where
Λ ( t ) = 0 τ q λ a ( τ ) i a ( τ , t ) + λ s ( τ ) i s ( τ , t ) + λ v a ( τ ) i v a ( τ , t ) + λ v s ( τ ) i v s ( τ , t ) d τ
is the weighted infectious population. The subscript represents a delay of days, so τ = τ , t = t .
The initial conditions are
S ( 0 ) = S 0 = b μ V 0 I 0 , v ( α , 0 ) = V 0 δ ( α α f ) , i a ( τ , 0 ) = ( 1 f ) I 0 δ ( τ ) , i s ( τ , 0 ) = f I 0 δ ( τ ) , i v a ( τ , 0 ) = 0 , i v s ( τ , 0 ) = 0 , R ( 0 ) = 0 , R v 1 ( 0 ) = 0 , R v 2 ( 0 ) = 0 .
The boundary conditions are
v ( 0 , t ) = ν S ( t ) , i a ( 0 , t ) = ( 1 f ) β Λ ( t ) S ( t ) + 0 α v 1 η ( α ) v ( α , t ) d α , i s ( 0 , t ) = f β Λ ( t ) S ( t ) + 0 α v 1 η ( α ) v ( α , t ) d α , i v a ( 0 , t ) = ( 1 f ) β Λ ( t ) α v 1 η ( α ) v ( α , t ) d α , i v s ( 0 , t ) = f β Λ ( t ) α v 1 η ( α ) v ( α , t ) d α ,
where f is the fraction of infections that are symptomatic, α v is the age of vaccination after which an individual’s virus load is reduced if infected, and α f is the age of vaccination after which vaccines provide full protection. For α < α v , vaccinated individuals are equivalent to susceptible individuals from the point of view of their virus profile, but they are less likely to get infected. Parameters { β , μ , m , f , α v } are taken from literature, and δ is the Dirac delta function.

2.1. Cumulative Statistics

In order to compare model results with commonly tabulated public health data, we define several cumulative population statistics derived from the model state variables.
The cumulative number of cases up to time t, Σ ( t ) , is given by the equation
Σ ( t ) = I 0 + 0 t i a ( 0 , s ) + i s ( 0 , s ) + i v a ( 0 , s ) + i v s ( 0 , s ) d s .
The cumulative number of positive case detections at time t, P ( t ) , is given by the equation
d P d t = 0 r a ( τ , t ) i a ( τ , t ) e ( μ + m ) + r s ( τ , t ) i s ( τ , t ) e ( μ + m ) + r v a ( τ , t ) i v a ( τ , t ) e ( μ + m ) + r v s ( τ , t ) i v s ( τ , t ) e ( μ + m ) d τ ,
with initial condition P ( 0 ) = 0 .
The cumulative number of SARS-CoV-2 naive individuals who have reached full vaccination status by time t is given by
F ( t ) = V 0 + 0 t v ( α f , s ) d s ,
and the cumulative number of previously infected individuals reaching full vaccination status by time t is given by
C V R ( t ) = 0 t ν R ( s α f ) e μ α f d s .
Thus, the cumulative number of fully vaccinated individuals at time t is
T ( t ) = F ( t ) + C V R ( t ) .
Breakthrough cases are infections of fully vaccinated individuals, so the cumulative number of breakthrough cases up to time t, B ( t ) , is the number of new infections occurring in individuals with α α f . B ( t ) is given by the equation
d B d t = β Λ ( t ) α f 1 η ( α ) v ( α , t ) d α ,
with initial condition B ( 0 ) = 0 .

2.2. Parameter Choices

We parametrize our model based on efficacy data from Pfizer-BioNTech vaccine against the Alpha and Delta variants [19,28]. Therefore, we set the rate 0 η ( α ) 1 describing the degree of protection against the Alpha variant α days after vaccine administration to
η ( α ) = 0 , for α < 7 , 0.49 , for 7 α < α f , 0.94 , for α α f ,
and against Delta variant, α days after vaccine administration, to
η ( α ) = 0 , for α < 7 , 0.33 , for 7 α < α f , 0.88 , for α α f .
The CDC defines a fully vaccinated individual as one who is ≥14 days past receipt of all recommended doses of an FDA-authorized vaccine [22]. In our model, this corresponds to α α f = 35 days. Levine et al. have shown that, for the first 11 days following the Pfizer-BioNTech vaccination, the cycle threshold (CT) values in infected vaccinated individuals do not change compared to those of infected yet non-vaccinated individuals [27]. We therefore set α v = 11 days. Afterwards, virus levels decrease by four-fold in vaccinated individuals [27]. The reduction in virus load leads to both shorter infectivity period, ranging between τ 2 j = 2.8 and τ 3 j = 9.8 days post infection, and shorter time for detection by tests, ranging between τ 1 j = 0.97 and τ 4 j = 10.22 days post infection ( j = { v a , v s } ) for RT-PCR (see Figure 1). We assume a PCR test return delay of = 1 days.
We assume that the Alpha variant transmissibility rate is β = 0.25 , as in our previous study [24], and the Delta variant transmissibility rate is 1.6-times higher, β = 0.4 , in agreement with recent studies [20,21]. The other parameters { μ , m , f , τ q } are as in our previous study [24]. A summary of parameters is given in Table 1.

3. Results

3.1. Alpha Variant Dynamics in the Absence of Testing

To simulate an initially undetected outbreak in a partially vaccinated population, we assumed that at time t 0 = 0 , when a total of V ( 0 ) = V 0 = 30 % of individuals have been fully vaccinated, I 0 = 5 % of the population was infected with the Alpha variant, f = 70 % with symptomatic and 1 f = 30 % with asymptomatic disease. The vaccine efficacy was given by (11), ν = 1 % additional vaccines were administered daily starting with day t 1 = 20 , and no testing was considered. Daily symptomatic, asymptomatic, breakthrough asymptomatic and breakthrough symptomatic cases peaked at 10.6 % , 4.6 % , 0.31 % and 0.13 % of the entire population at days 14, 14, 47 and 47, respectively (see Figure 2A, left panel, red, blue, magenta and cyan lines). At day 100, a cumulative Σ n o T e s t s ( 100 ) = 38.4 % of the population had been infected when not fully vaccinated (see Figure 2A, middle panel, magenta line) and a cumulative B ( 100 ) = 1.69 % of the population had been infected while fully vaccinated (see Figure 2A, middle panel, green line). Lastly, at day 100 a cumulative F ( 100 ) = 42 % of the naive population has been fully vaccinated (see Figure 2A, right panel, cyan line), a cumulative C V R ( 100 ) = 9.8 % of the population had been fully vaccinated after recovering from natural infection (see Figure 2A, right panel, magenta line), and a cumulative T ( 100 ) = 52 % of the population had been vaccinated (see Figure 2, right panel, black line). A quantity of interest is the percent breakthrough number, B ( 100 ) / F ( 100 ) , defined as the percent of cumulative naive fully vaccinated individuals who get infected divided by the cumulative fully vaccinated population. In this case, we obtained a percent breakthrough case rate B ( 100 ) / F ( 100 ) = 4 % at day 100.
To more closely determine the relationship between the percent breakthrough cases at day 100, B ( 100 ) / F ( 100 ) , the vaccination level at the start of the outbreak, V 0 , and daily vaccination level starting at day 20, ν , we derived a heat map for the percent breakthrough cases for smaller V 0 and ν increments (see Figure 3A). We observed that the percent breakthrough cases range between 0.72 % for V 0 = 80 % and ν = 5 % and 6.56 % for V 0 = 10 % and ν = 1 % (see Figure 3A). Having a large percent of the population vaccinated at the time of the outbreak resulted in a decrease in the proportion of fully vaccinated individuals experiencing breakthrough cases.

3.2. Alpha Variant Model Outcomes in the Presence of Testing

We investigated the effect of fixed daily testing with capacity C, administered beginning at day t 1 = 20 in two cases. Case 1: only the non-vaccinated group was tested; Case 2: both the vaccinated and non-vaccinated groups were tested. In particular, we quantified the effects of RT-PCR tests with return delay of one day in reducing total infections 100 days after the outbreak, Σ n o T e s t ( 100 ) Σ T e s t ( 100 ) .
For low vaccinated population, increased testing resulted in increased reduction in cumulative cases, regardless of the testing strategy. In particular, for V 0 = 10 % vaccinated population and no testing, the outbreak resulted in Σ n o T e s t s ( 100 ) = 63 % infections by day 100 (see Figure 4, grey heatmaps). When testing only non-vaccinated individuals, this number was reduced by 5.9 , 16.4 and 28.3 percentage points for C = 0.03 , C = 0.1 and C = 0.5 , respectively (see Figure 4A, colored heatmaps). When testing everyone, this number was reduced by 5.1 , 14.8 and 28.2 percentage points for C = 0.03 , C = 0.1 and C = 0.5 , respectively (see Figure 4B, colored heatmaps).
For high vaccinated population, increased testing had no effect on cumulative cases, regardless of the testing strategy. In particular, for V 0 = 80 % vaccinated population and no testing, the outbreak resulted in Σ n o T e s t s ( 100 ) = 7.6 % infections by day 100 (see Figure 4, grey heatmaps). Testing affected this number by 0.1 percentage points for C 0.03 when testing non-vaccinated individuals and for C 0.3 when testing everyone, respectively (see Figure 4, colored heatmaps).
Note that, for both testing scenarios, the amount of reduction decreased with the increase of V 0 (for a fixed ν ) and increased with the increase of testing capacity C (for a fixed V 0 ). Moreover, the highest decrease happened for the highest testing in the least vaccinated population. Lastly, in all cases, the model predicted that testing only non-vaccinated results in improved outcomes.

3.3. Delta Variant Dynamics in the Absence of Testing

We examined the effect of the Pfizer–BioNTech vaccine on blocking infections with the Delta variant. To account for the increased infectiousness of Delta, we increased the transmissibility rate to β = 0.4 , 1.6 -times higher than in the case of the Alpha variant. Moreover, we modified the Pfizer–BioNTech effectiveness to (12).
As before, we simulated an initially undetected outbreak in a partially vaccinated population, with V 0 = 30 % of the population fully vaccinated at the time of the outbreak, ν = 1 % additional vaccines were administered daily starting with day t 1 = 20 , and no testing was considered. Daily symptomatic, asymptomatic, breakthrough symptomatic and breakthrough asymptomatic cases peaked at 21.7 % , 9.2 % , 2.2 % and 0.9 % of the entire population at days 22, 22, 29 and 29, respectively (see Figure 2B, left panel, red, blue, magenta and cyan lines). At day 100, a cumulative Σ n o T e s t s ( 100 ) = 66.7 % of the population had been infected when not fully vaccinated (see Figure 2B, middle panel, magenta line) and a cumulative B ( 100 ) = 6.44 % of the population has been infected while fully vaccinated (see Figure 2B, middle panel, green line). Lastly, at day 100, the model predicted that a cumulative F ( 100 ) = 33.5 % of the naive population had been fully vaccinated (see Figure 2B, right panel, cyan line), a cumulative C V R ( 100 ) = 17.2 % of the population had been fully vaccinated after recovering from natural infection (see Figure 2B, right panel, magenta line), and a cumulative T ( 100 ) = 51 % of the population had been vaccinated (see Figure 2B, right panel, black line). The percent breakthrough case rate at day 100 became B ( 100 ) / F ( 100 ) = 19.22 % .
When we expanded our analysis to include smaller V 0 and ν increments, the model prediction for breakthrough cases at day 100 increased to B ( 100 ) / F ( 100 ) = 21.3 % (compared to 6.56 % for the Alpha variant) for V 0 = 10 % and ν = 1 % and to B ( 100 ) / F ( 100 ) = 3.94 % (compared to 0.72 % for the Alpha variant) for V 0 = 80 % and ν = 5 % (see Figure 3B). For this more contagious variant, having a large percent of the population vaccinated at the time of the outbreak significantly decreased the breakthrough cases.

3.4. Delta Variant Model Outcomes in the Presence of Testing

We next investigated the effect of fixed daily testing with capacity C, administered beginning at day t 1 = 20 in reducing total infections 100 days after the outbreak, Σ n o T e s t ( 100 ) Σ T e s t ( 100 ) , when only the non-vaccinated and when everyone was tested.
When less than 60 % of the population was vaccinated at the time of the outbreak, increased testing led to increased reduction in cumulative cases, regardless of the testing strategy. In particular, for V 0 = 10 % vaccinated population and no testing, the outbreak resulted in Σ n o T e s t ( 100 ) = 85.6 % of the population being infected by day 100 (see Figure 5, grey heatmaps). When testing non-vaccinated individuals only, this number was reduced by 1.9 , 6.5 and 18.5 percentage points for C = 0.03 , C = 0.1 and C = 0.5 , respectively (see Figure 5A, colored heatmaps). When testing everyone, this number was reduced by 1.7 , 5.9 and 18.5 percentage points for C = 0.03 , C = 0.1 and C = 0.5 , respectively (see Figure 5B, colored heatmaps). For V 0 = 30 % and no testing, the outbreak resulted in Σ n o T e s t ( 100 ) = 66.7 % of the population being infected by day 100 (see Figure 5, grey heatmaps). When testing non-vaccinated individuals only, this number was reduced by 3.5 , 10.5 and 19.9 percentage points for C = 0.03 , C = 0.1 and C = 0.5 , respectively (see Figure 5A, colored heatmaps). When testing everyone, this number was reduced by 2.6 , 8.3 and 20.2 percentage points for C = 0.03 , C = 0.1 and C = 0.5 , respectively (see Figure 5B, colored heatmaps). Finally, for V 0 = 60 % vaccinated population and no testing, the outbreak resulted in Σ n o T e s t ( 100 ) = 32.1 % of the population being infected by day 100 (see Figure 5, grey heatmaps). When testing only non-vaccinated individuals, this number was reduced by 4.2 , 8.1 and 9.9 percentage points for C = 0.03 , C = 0.1 and C = 0.5 , respectively (see Figure 5A, colored heatmaps). When testing everyone, this number was reduced by 2.3 , 5.9 and 10.4 percentage points for C = 0.03 , C = 0.1 and C = 0.5 , respectively (see Figure 5B, colored heatmaps). As with the Alpha variant, the amount of reduction increased with the increase of C (for a fixed V 0 ), regardless of the testing scenarios. Unlike the Alpha variant case, however, the reduction did not decrease monotonically with the increase of V 0 (for fixed C). The maximum reduction occured at different V 0 values for different testing capacities C. In particular, for both testing strategies, the maximum reductions occured at V 0 = 40 % for C = 0.1 and V 0 = 30 % for C = 0.3 .
When more than 60 % of the population was vaccinated, increased testing resulted in increased reduction in cumulative cases when we tested everyone, but not when we only tested the non-vaccinated. That is why, for high vaccinated populations, testing everyone led to a slight increase in case reductions, compared with testing only the non-vaccinated.

3.5. The Role of Testing for Emerging Variant Dynamics

In an ongoing and evolving pandemic, the specific characteristics of the virus and of human behavior are unknown. Hence, the parameters upon which our model relies are uncertain. We therefore investigated the role of testing in reducing cases when we assumed outbreaks with variants of increased infectiousness. First, we considered populations of different vaccination levels, where vaccines maintained high efficacy rates (11), and second, we considered a 50 % vaccinated population where variants induce decreased vaccines efficacy. In both cases we assumed a fixed testing capacity of C = 0.1 with a PCR test with a return delay of one day.
When the vaccine efficacy is unchanged by increased variant infectiousness, testing led to higher case reduction for intermediate vaccination levels. The peak of the reduction decreased and shifted to the right as the infectiousness level increased (see Figure 6). For example, for high infectivity rate β = 0.9 we obtained maximum testing effect for a population with 80 % background vaccination level, while for a lower infectivity rate β = 0.3 we obtained maximum testing effect for a population with 40 % background vaccination level. The percent case reduction was higher when we tested only the non-vaccinated population but the peak reduction occured at the same background vaccination level, regardless of the testing strategy (see Figure 6A,B).
We next investigated the effect of testing in reducing infectious cases when the vaccine efficacy decreased. We fixed the background vaccination level to V 0 = 0.5 and determined the percent case reduction for different β values when the efficacy of the second dose varied between η 2 = 0.35 and η 2 = 0.95 (see Figure 7). For β values for which R 0 < 1 , the effect of testing decreased with increasing vaccine efficacy. For β values for which R 0 > 1 , the effect of testing increased as vaccine efficacy increased. For R 0 > 1 , testing had maximum effect for variants with R 0 close to 1 ( β = 0.2 , 0.3 ), regardless of testing strategy (see Figure 7A,B), and testing non-vaccinated individuals alone only marginally improved the outcomes.

4. Discussion

Reaching herd-immunity to SARS-CoV-2 through vaccination and natural infection is made harder by vaccine hesitancy, vaccine shortage around the world, emergence of new variants, and delayed vaccination approval for children aged 5 to 11 years. Therefore, additional public health interventions such as masks, social distancing, and surveillance testing are still needed. In previous work, we have used mathematical models to show that surveillance testing can be an effective public health intervention to reduce a SARS-CoV-2 outbreak in a non-vaccinated population [24]. The significance of surveillance testing needs to be re-evaluated in the context of vaccination, emergence of new variants, and waning immunity. CDC recommends that vaccinated people are tested only when experiencing COVID-19 symptoms or after coming into close contact with someone with COVID-19 [33]. To better determine the role of surveillance tests in a rapidly changing environment, we developed a multi-scale mathematical model of COVID-19 transmission in a mixed vaccinated and non-vaccinated population. Our model investigated how vaccination levels impact the effectiveness of testing. We compared two testing strategies: one in which tests were administered to both vaccinated and non-vaccinated individuals and one in which tests were administered to non-vaccinated individuals alone. Additionally, we considered variants of increased infectiousness.
In the case of variants with low transmission rates, such as the Alpha variant, where vaccines are highly effective in blocking transmission, we found that testing remained an effective intervention when the overall vaccination is low to moderate. For higher vaccination levels, the impact of testing was diminished, even relative to the smaller outbreaks that occured in those scenarios. For any fixed testing capacity, the number of cases prevented by testing decreased with increasing vaccination level.
For variants with higher transmission rates, such as the Delta variant, where vaccine efficacy in blocking transmission was reduced, a more complex pattern of testing effectiveness was apparent. For low vaccination, the impact of testing was low, as testing was not sufficient to overcome the force of infection created by the Delta variant. This difference to the case of the Alpha variant was a result of the magnitude of the outbreak, which was driven by increased transmissibility and the reduced effectiveness of the vaccine against the Delta variant. Interestingly, as the vaccination level increased, the number of cases prevented by testing increased as well, even though the number of cases that would occur in the absence of testing declined. As vaccination level increased further, the effectiveness of a fixed testing capacity declined again, due to the significant reduction in the number of cases expected in the absence of testing, as in the case of the Alpha variant. Thus, when considering the Delta variant, the impact of a fixed testing capacity was highest for moderate vaccination level and lower for low and high vaccination level.
We investigated differences in testing strategies. In the United States, it is common for surveillance testing to focus on the non-vaccinated population [33]. Here we compared the differences in outcomes for testing strategies that included only the non-vaccinated and those that included both vaccinated and non-vaccinated populations. We found that testing strategies that focused on the non-vaccinated population were generally more effective than broad testing strategies. For the Alpha variant, this was true for all cases considered. For the Delta variant, a broad testing strategy was preferable when the testing capacity exceeded the non-vaccinated population, a result that has been reported in other modeling studies [34]. This indicated that the most effective strategy should focus first on testing the non-vaccinated population, then using excess capacity in the non-vaccinated population.
As expected, for fixed vaccination level, increasing testing capacity increased the number of cases prevented. However, for high vaccination level, the impact of additional testing eventually saturated, and further testing had little or no effect.
We have also studied the prevalence of infection within the fully vaccinated population, so-called “breakthrough” cases. For the Alpha variant, the prevalence of breakthrough cases was uniformly low (ranging between 0.7 % and 6.45 % ), and decreased as the population vaccination level increased. These results are similar to those from clinical studies which reported an Alpha variant incidence of breakthrough infections of 0.5 % in health care workers in US who received either the Pfizer–BioNTech or the Moderna vaccine [35], of 2.6 % in health care workers in Israel who received both doses of Pfizer–BioNTech vaccine [36], and of 0.08 % in New York metropolitan area vaccinated with either Pfizer–BioNTech, Moderna, or Janssen vaccines [37]. For the Delta variant, our model predicted that breakthrough cases were much more prevalent in the vaccinated population and ranged between 3.94 % and 21.3 % . These numbers are comparable with the 8.4 % breakthrough cases reported in Houston hospitals [38] and lower than the 28 % reported in the Washington, DC area [39]. While the increased transmissibility of the Delta variant and the decreased effectiveness of vaccine were necessary conditions for this increase in breakthrough cases, the primary driving factor was the extent of the outbreak in the non-vaccinated population. As the vaccinated population became larger, both the total number of breakthrough cases and their prevalence decreased.
The results for the Alpha and Delta variants presented in this study were based on uncertain parameters. To better describe the dependence of our model predictions on evolving viral parameters, we quantified percent case reductions due to testing for variants of increased infectiousness and for variants that make vaccines less effective in preventing infection. We found that increased variant transmission rates (which may be due to the evolution of inherent viral characteristics influencing infectiousness or leading to higher viral load, or due to changing compliance with public health interventions) resulted in decreased testing efficacy. Interestingly, as variant transmission rates increased, testing was most useful when applied to populations with high vaccination background. The decrease in testing effectiveness can be compensated by increasing the testing frequency and decreasing the delay in results return (not shown). Applying these additional measures is especially important when variants of higher infectiousness become dominant, as they lead to larger and faster spreading outbreaks.
Similarly, we found that increased variant transmission rates resulted in decreased testing efficacy, regardless of the efficacy of the second dose of vaccine. The impact of second dose vaccine efficacy was less pronounced than the impact of the background vaccination rate. This indicates that even if vaccine efficacy wanes over time in individuals, maintaining a high vaccination rate in the population is an important factor in ensuring the effectiveness of testing strategies for outbreak reduction. As with the Alpha and Delta scenarios, focusing testing on non-vaccinated individuals resulted in improved outcomes, regardless of transmission rate or vaccine efficacy.
Our modeling approach includes several simplifying assumptions, some of which can be relaxed to generalize our results in a variety of ways. First, it is important to note that, in this study, the measure for the effectiveness of testing is the prevention of transmission, not the prevention of disease, hospitalization, or death. The true extent of transmission is generally not known in an ongoing outbreak, even when testing is widespread. Public health outcomes and vaccine effectiveness are generally expressed in terms of preventing disease, with all three vaccines approved for use in the United States being highly effective in preventing hospitalization and deaths regardless of the variant [40]. While preventing disease is the immediate goal, we argue that preventing transmission even when a large fraction of the population has been vaccinated (through surveillance testing) is a worthy long-term goal that may allow emergent variant strains resistant to vaccination to go extinct before becoming the next dominant strain [41]. Secondly, we assumed that a positive test return results in the complete isolation of the infected individual, who can no longer transmit the virus. This is a best case scenario, and further work is needed to incorporate other metrics, such as compliance rate.
We have started the outbreak with a large number of infectious individuals. When we decreased that initial burden from 5 % to 1 % , the epidemic curve flattened and the time to the peak lengthened. As a result, the effect of testing in lower vaccinated populations increased significantly (not shown). As before, the effect of testing in the highly vaccinated population saturated. This indicates that early detection of nascent outbreaks is essential for the effectiveness of testing as a public health intervention.
Our study is limited to PCR testing with a return delay of one day. In the past, we have investigated the trade-off between employing faster, cheaper, yet less sensitive antigen tests at a more frequent rate. When the same testing capacity was used, the antigen tests underperformed the PCR tests with one day delay (not shown). This was due to decreased interval of detection for the antigen tests ( τ 1 v = 3.2 and τ 4 v = 6.7 days post infection). If the frequency of antigen test administration was increased, or in places where the PCR return is long, antigen test can present a reliable alternative for surveillance testing.
Lastly, our model assumed that the individuals that received their first vaccine dose and then get infected within the first α v = 11 days have the same virus profile as a non-vaccinated individual. The results of this study are not sensitive to the value of α v (not shown). In the model, such individuals may then complete vaccination following recovery, which may differ from current public health practices.

5. Conclusions

In summary, we have developed a multi-scale model of SARS-CoV-2 transmission in a vaccinated population. We found that, when variants with low transmission rates and for which vaccines are highly efficacious (such as the Alpha variant) are dominant, testing is effective when vaccination levels are low to moderate and its impact is diminished when vaccination levels are high. When variants with increased transmission rates and for which vaccines are less efficacious (such as the Delta variant) are dominant, widespread vaccination is necessary in order to prevent significant outbreaks. When only moderate vaccination level can be achieved, frequent testing can significantly reduce the cumulative size of the outbreak, and the impact of testing is greatest when it is focused on the remaining non-vaccinated population.

Author Contributions

J.E.F. and S.M.C. conceived the study. J.E.F. wrote the code. J.E.F. and S.M.C. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Science Foundation grants No. DMS-1813011 and DMS-2051820.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The code generated will be available upon acceptance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A novel coronavirus from patients with pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727–733. [Google Scholar] [CrossRef] [PubMed]
  2. WHO Cases Dashboard. Available online: https://covid19.who.int/ (accessed on 29 August 2021).
  3. Legido-Quigley, H.; Asgari, N.; Teo, Y.Y.; Leung, G.M.; Oshitani, H.; Fukuda, K.; Cook, A.R.; Hsu, L.Y.; Shibuya, K.; Heymann, D. Are high-performing health systems resilient against the COVID-19 epidemic? Lancet 2020, 395, 848–850. [Google Scholar] [CrossRef] [Green Version]
  4. Rocha, R.; Atun, R.; Massuda, A.; Rache, B.; Spinola, P.; Nunes, L.; Lago, M.; Castro, M.C. Effect of socioeconomic inequalities and vulnerabilities on health-system preparedness and response to COVID-19 in Brazil: A comprehensive analysis. Lancet Glob. Health 2021, 9, e782–e792. [Google Scholar] [CrossRef]
  5. Haldane, V.; De Foo, C.; Abdalla, S.M.; Jung, A.S.; Tan, M.; Wu, S.; Chua, A.; Verma, M.; Shrestha, P.; Singh, S.; et al. Health systems resilience in managing the COVID-19 pandemic: Lessons from 28 countries. Nat. Med. 2021, 27, 964–980. [Google Scholar] [CrossRef]
  6. Ingram, C.; Downey, V.; Roe, M.; Chen, Y.; Archibald, M.; Kallas, K.A.; Kumar, J.; Naughton, P.; Uteh, C.O.; Rojas-Chaves, A.; et al. COVID-19 prevention and control measures in workplace settings: A rapid review and meta-analysis. Int. J. Environ. Res. Public Health 2021, 18, 7847. [Google Scholar] [CrossRef]
  7. Perra, N. Non-pharmaceutical interventions during the COVID-19 pandemic: A review. Phys. Rep. 2021, 913, 1–52. [Google Scholar] [CrossRef]
  8. Matrajt, L.; Leung, T. Evaluating the effectiveness of social distancing interventions to delay or flatten the epidemic curve of coronavirus disease. Emerg. Infect. Dis. 2020, 26, 1740. [Google Scholar] [CrossRef]
  9. Kumar, A.; Dowling, W.E.; Román, R.G.; Chaudhari, A.; Gurry, C.; Le, T.T.; Tollefson, S.; Clark, C.E.; Bernasconi, V.; Kristiansen, P.A. Status report on COVID-19 vaccines development. Curr. Infect. Dis. Rep. 2021, 23, 9. [Google Scholar] [CrossRef]
  10. Callaway, E. The race for coronavirus vaccines: A graphical guide. Nature 2020, 580, 576–577. [Google Scholar] [CrossRef]
  11. Polack, F.P.; Thomas, S.J.; Kitchin, N.; Absalon, J.; Gurtman, A.; Lockhart, S.; Perez, J.L.; Marc, G.P.; Moreira, E.D.; Zerbini, C.; et al. Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine. N. Engl. J. Med. 2020, 383, 2603–2615. [Google Scholar] [CrossRef]
  12. U.S. COVID-19 Vaccine Product Information. Available online: https://www.cdc.gov/vaccines/covid-19/info-by-product/index.html (accessed on 29 August 2021).
  13. Tenforde, M.W. Sustained Effectiveness of Pfizer-BioNTech and Moderna Vaccines Against COVID-19 Associated Hospitalizations Among Adults—United States, March–July 2021. Morb. Mortal. Wkly. Rep. 2021, 70, 1156. [Google Scholar] [CrossRef] [PubMed]
  14. Bernal, J.L.; Andrews, N.; Gower, C.; Robertson, C.; Stowe, J.; Tessier, E.; Simmons, R.; Cottrell, S.; Roberts, R.; O’Doherty, M.; et al. Effectiveness of the Pfizer-BioNTech and Oxford-AstraZeneca vaccines on covid-19 related symptoms, hospital admissions, and mortality in older adults in England: Test negative case-control study. BMJ 2021, 373, n1088. [Google Scholar] [CrossRef] [PubMed]
  15. Haas, E.J.; Angulo, F.J.; McLaughlin, J.M.; Anis, E.; Singer, S.R.; Khan, F.; Brooks, N.; Smaja, M.; Mircus, G.; Pan, K.; et al. Impact and effectiveness of mRNA BNT162b2 vaccine against SARS-CoV-2 infections and COVID-19 cases, hospitalisations, and deaths following a nationwide vaccination campaign in Israel: An observational study using national surveillance data. Lancet 2021, 397, 1819–1829. [Google Scholar] [CrossRef]
  16. Tenforde, M.W.; Patel, M.M.; Ginde, A.A.; Douin, D.J.; Talbot, H.K.; Casey, J.D.; Mohr, N.M.; Zepeski, A.; Gaglani, M.; McNeal, T.; et al. Effectiveness of SARS-CoV-2 mRNA vaccines for preventing Covid-19 hospitalizations in the United States. Clin. Infect. Dis. 2021. [Google Scholar] [CrossRef]
  17. CDC Covid-19 Variants Proportion in USA. Available online: https://covid.cdc.gov/covid-data-tracker/#variant-proportions (accessed on 23 September 2021).
  18. Fowlkes, A. Effectiveness of COVID-19 Vaccines in Preventing SARS-CoV-2 Infection Among Frontline Workers Before and During B. 1.617. 2 (Delta) Variant Predominance—Eight US Locations, December 2020–August 2021. MMWR. Morb. Mortal. Wkly. Rep. 2021, 70, 1167. [Google Scholar] [CrossRef]
  19. Lopez Bernal, J.; Andrews, N.; Gower, C.; Gallagher, E.; Simmons, R.; Thelwall, S.; Stowe, J.; Tessier, E.; Groves, N.; Dabrera, G.; et al. Effectiveness of Covid-19 vaccines against the B.1.617.2 (delta) variant. N. Engl. J. Med. 2021, 385, 585–594. [Google Scholar] [CrossRef]
  20. Liu, Y.; Rocklöv, J. The reproductive number of the Delta variant of SARS-CoV-2 is far higher compared to the ancestral SARS-CoV-2 virus. J. Travel Med. 2021, 28, taab124. [Google Scholar] [CrossRef]
  21. Earnest, R.; Uddin, R.; Matluk, N.; Renzette, N.; Siddle, K.J.; Loreth, C.; Adams, G.; Tomkins-Tinch, C.; Petrone, M.E.; Rothman, J.E.; et al. Comparative transmissibility of SARS-CoV-2 variants Delta and Alpha in New England, USA. medRxiv 2021. [Google Scholar] [CrossRef]
  22. Birhane, M.; Bressler, S.; Chang, G.; Clark, T.; Dorough, L.; Fischer, M.; Watkins, L.F.; Goldstein, J.M.; Kugeler, K.; Langley, G.; et al. COVID-19 Vaccine Breakthrough Infections Reported to CDC—United States, January 1–April 30, 2021. Morb. Mortal. Wkly. Rep. 2021, 70, 792. [Google Scholar]
  23. Pfizer-BioNTech Booster Approval. Available online: https://www.fda.gov/news-events/press-announcements/fda-authorizes-booster-dose-pfizer-biontech-covid-19-vaccine-certain-populations (accessed on 3 September 2021).
  24. Forde, J.; Ciupe, S. Quantification of the tradeoff between test sensitivity and test frequency in a COVID-19 epidemic-a multi-scale modeling approach. Viruses 2021, 13, 457. [Google Scholar] [CrossRef]
  25. Nikin-Beers, R.; Blackwood, J.C.; Childs, L.M.; Ciupe, S.M. Unraveling within-host signatures of dengue infection at the population level. J. Theor. Biol. 2018, 446, 79–86. [Google Scholar] [CrossRef]
  26. Dorratoltaj, N.; Nikin-Beers, R.; Ciupe, S.M.; Eubank, S.G.; Abbas, K.M. Multi-scale immunoepidemiological modeling of within-host and between-host HIV dynamics: Systematic review of mathematical models. PeerJ 2017, 5, e3877. [Google Scholar] [CrossRef]
  27. Levine-Tiefenbrun, M.; Yelin, I.; Katz, R.; Herzel, E.; Golan, Z.; Schreiber, L.; Wolf, T.; Nadler, V.; Ben-Tov, A.; Kuint, J.; et al. Initial report of decreased SARS-CoV-2 viral load after inoculation with the BNT162b2 vaccine. Nat. Med. 2021, 75, 790–792. [Google Scholar] [CrossRef] [PubMed]
  28. Pfizer-BioNTech COVID-19 Vaccine EUA. Available online: https://www.fda.gov/media/144414/download (accessed on 3 September 2021).
  29. Ke, R.; Zitzmann, C.; Ribeiro, R.M.; Perelson, A.S. Kinetics of SARS-CoV-2 infection in the human upper and lower respiratory tracts and their relationship with infectiousness. medRxiv 2020. [Google Scholar] [CrossRef]
  30. He, D.; Zhao, S.; Lin, Q.; Zhuang, Z.; Cao, P.; Wang, M.H.; Yang, L. The relative transmissibility of asymptomatic COVID-19 infections among close contacts. Int. J. Infect. Dis. 2020, 94, 145–147. [Google Scholar] [CrossRef] [PubMed]
  31. Wölfel, R.; Corman, V.M.; Guggemos, W.; Seilmaier, M.; Zange, S.; Müller, M.A.; Niemeyer, D.; Jones, T.C.; Vollmar, P.; Rothe, C.; et al. Virological assessment of hospitalized patients with COVID-2019. Nature 2020, 581, 465–469. [Google Scholar] [CrossRef] [Green Version]
  32. He, X.; Lau, E.H.; Wu, P.; Deng, X.; Wang, J.; Hao, X.; Lau, Y.C.; Wong, J.Y.; Guan, Y.; Tan, X.; et al. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nat. Med. 2020, 26, 672–675. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. CDC Guidelines for Vaccinated People. Available online: https://www.cdc.gov/coronavirus/2019-ncov/vaccines/fully-vaccinated-guidance.html (accessed on 3 September 2021).
  34. Motta, F.C.; McGoff, K.A.; Deckard, A.; Wolfe, C.R.; Moody, M.A.; Cavanaugh, K.; Denny, T.N.; Harer, J.; Haase, S.B. Benefits of Surveillance Testing and Quarantine in a SARS-CoV-2 Vaccinated Population of Students on a University Campus. medRxiv 2021. [Google Scholar] [CrossRef]
  35. Hacisuleyman, E.; Hale, C.; Saito, Y.; Blachere, N.E.; Bergh, M.; Conlon, E.G.; Schaefer-Babajew, D.J.; DaSilva, J.; Muecksch, F.; Gaebler, C.; et al. Vaccine breakthrough infections with SARS-CoV-2 variants. N. Engl. J. Med. 2021, 384, 2212–2218. [Google Scholar] [CrossRef]
  36. Bergwerk, M.; Gonen, T.; Lustig, Y.; Amit, S.; Lipsitch, M.; Cohen, C.; Mandelboim, M.; Gal Levin, E.; Rubin, C.; Indenbaum, V.; et al. COVID-19 breakthrough infections in vaccinated health care workers. N. Engl. J. Med. 2021, 385, 1474–1484. [Google Scholar] [CrossRef] [PubMed]
  37. Duerr, R.; Dimartino, D.; Marier, C.; Zappile, P.; Wang, G.; Lighter, J.; Elbel, B.; Troxel, A.B.; Heguy, A. Dominance of Alpha and Iota variants in SARS-CoV-2 vaccine breakthrough infections in New York City. J. Clin. Investig. 2021, 131, e152702. [Google Scholar] [CrossRef] [PubMed]
  38. Musser, J.M.; Christensen, P.A.; Olsen, R.J.; Long, S.W.; Subedi, S.; Davis, J.J.; Hodjat, P.; Walley, D.R.; Kinskey, J.C.; Gollihar, J.D. Delta variants of SARS-CoV-2 cause significantly increased vaccine breakthrough COVID-19 cases in Houston, Texas. medRxiv 2021. [Google Scholar] [CrossRef]
  39. Luo, C.H.; Morris, C.P.; Sachithanandham, J.; Amadi, A.; Gaston, D.; Li, M.; Swanson, N.J.; Schwartz, M.; Klein, E.Y.; Pekosz, A.; et al. Infection with the SARS-CoV-2 Delta Variant is Associated with Higher Infectious Virus Loads Compared to the Alpha Variant in both Unvaccinated and Vaccinated Individuals. medRxiv 2021. [Google Scholar] [CrossRef]
  40. Moline, H.L.; Whitaker, M.; Deng, L.; Rhodes, J.C.; Milucky, J.; Pham, H.; Patel, K.; Anglin, O.; Reingold, A.; Chai, S.J.; et al. Effectiveness of COVID-19 vaccines in preventing hospitalization among adults aged ≥ 65 years—COVID-NET, 13 states, February–April 2021. Morb. Mortal. Wkly. Rep. 2021, 70, 1088. [Google Scholar] [CrossRef]
  41. Rella, S.A.; Kulikova, Y.A.; Dermitzakis, E.T.; Kondrashov, F.A. Rates of SARS-CoV-2 transmission and vaccination impact the fate of vaccine-resistant strains. Sci. Rep. 2021, 11, 15729. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Virus profiles in non-vaccinated and vaccinated individuals. log 10 virus load per swab over time during Alpha variant natural infection (red line) and vaccination (green line) as given by the within-host model in [29]. Non-vaccinated patients are assumed to be infectious from t = 2.5 days till t = 10.5 days (shaded pink region). Vaccinated patients are assumed to be infectious from t = 2.85 days till t = 9.8 days (shaded purple region). Black horizontal lines correspond to RT-PCR test detection threshold (LOD) of log 10 ( V ) = 2 per swab and antigen test detection threshold (LOD) of log 10 ( V ) = 5 per swab.
Figure 1. Virus profiles in non-vaccinated and vaccinated individuals. log 10 virus load per swab over time during Alpha variant natural infection (red line) and vaccination (green line) as given by the within-host model in [29]. Non-vaccinated patients are assumed to be infectious from t = 2.5 days till t = 10.5 days (shaded pink region). Vaccinated patients are assumed to be infectious from t = 2.85 days till t = 9.8 days (shaded purple region). Black horizontal lines correspond to RT-PCR test detection threshold (LOD) of log 10 ( V ) = 2 per swab and antigen test detection threshold (LOD) of log 10 ( V ) = 5 per swab.
Viruses 13 02546 g001
Figure 2. Dynamics of Alpha and Delta variants infection over time. Left panels: daily asymptomatic (blue), symptomatic (red), breakthrough asymptomatic (cyan), and breakthrough asymptomatic (magenta) infections over time; Middle panels: cumulative cases Σ ( t ) (magenta), cumulative breakthrough cases B ( t ) (green) over time; Right panels: cumulative fully vaccinated F ( t ) (cyan), cumulative vaccinated after infection C V R ( t ) (magenta) and cumulative total vaccinated T ( t ) (black) over time in the absence of testing. Panel (A): Alpha variant; Panel (B): Delta variant. The background vaccination is 30 % and the other parameters and initial conditions are given in Table 1.
Figure 2. Dynamics of Alpha and Delta variants infection over time. Left panels: daily asymptomatic (blue), symptomatic (red), breakthrough asymptomatic (cyan), and breakthrough asymptomatic (magenta) infections over time; Middle panels: cumulative cases Σ ( t ) (magenta), cumulative breakthrough cases B ( t ) (green) over time; Right panels: cumulative fully vaccinated F ( t ) (cyan), cumulative vaccinated after infection C V R ( t ) (magenta) and cumulative total vaccinated T ( t ) (black) over time in the absence of testing. Panel (A): Alpha variant; Panel (B): Delta variant. The background vaccination is 30 % and the other parameters and initial conditions are given in Table 1.
Viruses 13 02546 g002
Figure 3. Percent breakthrough cases at day 100. Heatmaps for the percent breakthrough cases in the vaccinated population at day 100 versus additional daily vaccines, ν , and background vaccination levels, V 0 . Panel (A): Alpha variant; Panel (B): Delta variant. Parameters and initial conditions are given in Table 1.
Figure 3. Percent breakthrough cases at day 100. Heatmaps for the percent breakthrough cases in the vaccinated population at day 100 versus additional daily vaccines, ν , and background vaccination levels, V 0 . Panel (A): Alpha variant; Panel (B): Delta variant. Parameters and initial conditions are given in Table 1.
Viruses 13 02546 g003
Figure 4. Reduction in Alpha variant cases at 100 days. Heatmaps for the reduction in cumulative cases at 100 days after an outbreak with an Alpha variant, Σ n o T e s t s ( 100 ) Σ T e s t s ( 100 ) , as given by model Equation (1) versus RT-PCR testing capacity with a return delay of 1 day, C, and background vaccination levels, V 0 . Panel (A): Test non-vaccinated only; Panel (B): Test everybody. The gray heatmaps represent the cumulative cases at day 100 in the absence of testing, Σ n o T e s t ( 100 ) . Parameters and initial conditions are given in Table 1.
Figure 4. Reduction in Alpha variant cases at 100 days. Heatmaps for the reduction in cumulative cases at 100 days after an outbreak with an Alpha variant, Σ n o T e s t s ( 100 ) Σ T e s t s ( 100 ) , as given by model Equation (1) versus RT-PCR testing capacity with a return delay of 1 day, C, and background vaccination levels, V 0 . Panel (A): Test non-vaccinated only; Panel (B): Test everybody. The gray heatmaps represent the cumulative cases at day 100 in the absence of testing, Σ n o T e s t ( 100 ) . Parameters and initial conditions are given in Table 1.
Viruses 13 02546 g004
Figure 5. Reduction in Delta variant cases at 100 days. Heatmaps for the reduction in cumulative cases at 100 days after an outbreak with a Delta variant, Σ n o T e s t s ( 100 ) Σ T e s t s ( 100 ) , as given by model Equation (1) versus RT-PCR testing capacity with a return delay of 1 day, C, and background vaccination levels, V 0 . Panel (A): Test non-vaccinated only; Panel (B): Test everybody. The gray heatmaps represent the cumulative cases at day 100 in the absence of testing, Σ n o T e s t ( 100 ) . Parameters and initial conditions are given in Table 1.
Figure 5. Reduction in Delta variant cases at 100 days. Heatmaps for the reduction in cumulative cases at 100 days after an outbreak with a Delta variant, Σ n o T e s t s ( 100 ) Σ T e s t s ( 100 ) , as given by model Equation (1) versus RT-PCR testing capacity with a return delay of 1 day, C, and background vaccination levels, V 0 . Panel (A): Test non-vaccinated only; Panel (B): Test everybody. The gray heatmaps represent the cumulative cases at day 100 in the absence of testing, Σ n o T e s t ( 100 ) . Parameters and initial conditions are given in Table 1.
Viruses 13 02546 g005
Figure 6. Percent reduction in variant cases at 100 days as vaccination increases. Percent reduction in cumulative cases at 100 days after an outbreak with a variants of different infectivity rate β , ( Σ n o T e s t s ( 100 ) Σ T e s t s ( 100 ) ) / Σ n o T e s t s ( 100 ) , as given by model Equation (1) versus background vaccination levels, V 0 . Panel (A): Test non-vaccinated only; Panel (B): Test everybody. Parameters and initial conditions are given in Table 1 and we assume RT-PCR testing capacity C = 0.1 with a return delay of 1 day.
Figure 6. Percent reduction in variant cases at 100 days as vaccination increases. Percent reduction in cumulative cases at 100 days after an outbreak with a variants of different infectivity rate β , ( Σ n o T e s t s ( 100 ) Σ T e s t s ( 100 ) ) / Σ n o T e s t s ( 100 ) , as given by model Equation (1) versus background vaccination levels, V 0 . Panel (A): Test non-vaccinated only; Panel (B): Test everybody. Parameters and initial conditions are given in Table 1 and we assume RT-PCR testing capacity C = 0.1 with a return delay of 1 day.
Viruses 13 02546 g006
Figure 7. Percent reduction in variant cases at 100 days as vaccination efficacy decreases. Percent reduction in cumulative cases at 100 days after an outbreak with a variants of different infectivity rate β , ( Σ n o T e s t s ( 100 ) Σ T e s t s ( 100 ) ) / Σ n o T e s t s ( 100 ) , as given by model Equation (1) versus efficacy of the second dose of the vaccine, η 2 . Panel (A): Test non-vaccinated only; Panel (B): Test everybody. Parameters and initial conditions are given in Table 1, we assume background vaccination V 0 = 0.5 and RT-PCR testing capacity C = 0.1 with a return delay of 1 day.
Figure 7. Percent reduction in variant cases at 100 days as vaccination efficacy decreases. Percent reduction in cumulative cases at 100 days after an outbreak with a variants of different infectivity rate β , ( Σ n o T e s t s ( 100 ) Σ T e s t s ( 100 ) ) / Σ n o T e s t s ( 100 ) , as given by model Equation (1) versus efficacy of the second dose of the vaccine, η 2 . Panel (A): Test non-vaccinated only; Panel (B): Test everybody. Parameters and initial conditions are given in Table 1, we assume background vaccination V 0 = 0.5 and RT-PCR testing capacity C = 0.1 with a return delay of 1 day.
Viruses 13 02546 g007
Table 1. Parameter values and initial conditions used in model Equation (1).
Table 1. Parameter values and initial conditions used in model Equation (1).
Fixed ParametersDescriptionValueSource
β Transmission ratevaried
bBirth rate 7 × 10 4 /day[25]
μ Death rate 7 × 10 4 /day[25]
mDisease induced mortality rate 10 3 /day
fFraction of symptomatic infections 0.7 [30]
γ Relative asymp. infectiousness0.7
α v Age of vaccination for reduced virus11 days[27]
α f Age of vaccination corresponding to full protection35 days[12]
τ q Recovery age14 days[24]
PCR test return delay1 day
τ 1 j ( j = { a , s } )Age of onset of virus detectability in non-vaccinated0.554 days
τ 1 j ( j = { v a , v s } )Age of onset of virus detectability in vaccinated0.974 days
τ 2 j ( j = { a , s } )Age of onset of infectiousness in non-vaccinated 2.5 days[31]
τ 2 j ( j = { v a , v s } )Age of onset of infectiousness in vaccinated 2.8 days[27]
τ 3 j ( j = { a , s } )Age of end of infectiousness in non-vaccinated 10.5 days[30,32]
τ 3 j ( j = { v a , v s } )Age of end of infectiousness in vaccinated 9.8 days[27]
τ 4 j ( j = { a , s } )Age of loss of virus detectability in non-vaccinated10.95 days
τ 4 j ( j = { v a , v s } )Age of loss of virus detectability in vaccinated10.22 days
ν Vaccination ratevaried
t 1 Time when additional vaccination is initiated20 days
Ctesting capacityvaried
Initial ConditionsDescriptionValueSource
S ( 0 ) Susceptible population 0.95 V ( 0 )
V ( α , 0 ) = V ( 0 ) = V 0 Vaccination levelvaried
i s ( τ , 0 ) Infected symptomatic population 0.05 f δ ( τ )
i a ( τ , 0 ) Infected asymptomatic population 0.05 ( 1 f ) δ ( τ )
i v s ( τ , 0 ) Infected vaccinated symptomatic population0
i v a ( τ , 0 ) Infected vaccinated asymptomatic population0
R ( 0 ) Recovered population0
R v 1 ( 0 ) Vaccinated after natural infection0
R v 2 ( 0 ) Vaccinated, infected, and recovered0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Forde, J.E.; Ciupe, S.M. Modeling the Influence of Vaccine Administration on COVID-19 Testing Strategies. Viruses 2021, 13, 2546. https://0-doi-org.brum.beds.ac.uk/10.3390/v13122546

AMA Style

Forde JE, Ciupe SM. Modeling the Influence of Vaccine Administration on COVID-19 Testing Strategies. Viruses. 2021; 13(12):2546. https://0-doi-org.brum.beds.ac.uk/10.3390/v13122546

Chicago/Turabian Style

Forde, Jonathan E., and Stanca M. Ciupe. 2021. "Modeling the Influence of Vaccine Administration on COVID-19 Testing Strategies" Viruses 13, no. 12: 2546. https://0-doi-org.brum.beds.ac.uk/10.3390/v13122546

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop