Next Article in Journal
A Miniaturized Microbe-Silicon-Chip Based on Bioluminescent Engineered Escherichia coli for the Evaluation of Water Quality and Safety
Next Article in Special Issue
Detection of Nail Oncometabolite SAICAR in Oral Cancer Patients and Its Molecular Interactions with PKM2 Enzyme
Previous Article in Journal
Fiscal Policy Dilemma in Resolving Agricultural Risks: Evidence from China’s Agricultural Insurance Subsidy Pilot
Previous Article in Special Issue
COVID-19 Vaccine Priority Strategy Using a Heterogenous Transmission Model Based on Maximum Likelihood Estimation in the Republic of Korea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Projecting the Impact of SARS-CoV-2 Variants and the Vaccination Program on the Fourth Wave of the COVID-19 Pandemic in South Korea

Department of Mathematics, Soongsil University, Seoul 06978, Korea
Int. J. Environ. Res. Public Health 2021, 18(14), 7578; https://doi.org/10.3390/ijerph18147578
Submission received: 27 May 2021 / Revised: 13 July 2021 / Accepted: 14 July 2021 / Published: 16 July 2021
(This article belongs to the Special Issue Computational Modeling in Biology and Medicine)

Abstract

:
Vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are currently administered in South Korea; however, vaccine supply is limited. Considering constraints in vaccine supply and the emergence of variant strains, we evaluated the impact of coronavirus disease (COVID-19) vaccination program in reducing incidence, ICU hospitalization, and deaths in South Korea. We developed an age-structured model of SARS-CoV-2 transmission parameterized with Korean demographics and age-specific COVID-19 outcomes. Using our model, we analyzed the impact of the COVID-19 vaccination program during the fourth wave of the pandemic in South Korea in reducing disease burden. We projected that the vaccination program can reduce the overall attack rate to 3.9% from 6.9% without vaccination, over 150 days, starting from 5 July 2021. The highest relative reduction (50%) was observed among individuals aged 50–59 years. Vaccination markedly reduced adverse outcomes, such as ICU hospitalizations and deaths, decreasing them by 45% and 43%, respectively. In the presence of the Delta variant, vaccination is expected to reduce the overall attack rate to 11.9% from 26.9%. Our results indicate that the impact of vaccination can be substantially affected by the emergence of SARS-CoV-2 variants. Furthermore, herd immunity is unlikely to be achieved with the potential emergence of the Delta variant, inconsistent with the blueprint of the South Korean government.

1. Introduction

Despite unprecedented social distancing measures and vaccine development efforts, the rapid spread of coronavirus disease (COVID-19) continues to impose a heavy burden on healthcare at a global scale, resulting in 183 million reported cases and 3.9 million deaths worldwide, as of 4 July 2021 [1]. Since the first case of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was reported in January 2020, the virus has been constantly evolving, and several variants have emerged with mutations that alter the receptor binding domain of the spike protein [1,2]. Certain variants, called variants of concern (VOC), have shown their rapid emergence within populations, with increased potential for transmission or clinical implications, including the Alpha variant first identified in the United Kingdom (UK), the Beta variant in South Africa, the Gamma variant in Brazil, and the Delta variant in India [3]. Among these variants, the SARS-CoV-2 variant Alpha, also known as B.1.1.7, is estimated to have emerged in the UK in September 2020 and has become the dominant circulating SARS-CoV-2 variant in England [3]. The Alpha variant is found to have a selective advantage over the other variants, being 56% more transmissible (95% confidence interval (CI): 50–74%) across three regions, i.e., Denmark, Switzerland, and the United States (US), compared to the pre-existing variants of SARS-CoV-2 [4,5,6]. Furthermore, the Delta variant, also known as B.1.617.2, was first identified in India in December 2020 and has now spread to South Korea. Data from the UK suggest that the Delta variant is 40–60% more transmissible than the Alpha variant.
In South Korea, the detection rate of major virus variants (Alpha, Beta, Gamma, and Delta) in community transmission in June 2021 was over 30%; meanwhile, the cases of Delta variant have surged [7]. Furthermore, on 7 July 2021, it was announced that the fourth wave has hit South Korea, with more than 1200 new daily cases reported on average in July. The Korea Disease Control and Prevention Agency (KDCA) announced that over the first week of July 2021, 26.5% of newly confirmed COVID-19 patients in the capital area were infected with the Delta variant, expressing concern that the Delta variant cases may become prevalent in the country [7].
All authorized COVID-19 vaccines have demonstrated an efficacy of 65–95% against symptomatic, laboratory-confirmed COVID-19 caused by the non-variant strain [2]. However, the efficacy of COVID-19 vaccines might be potentially affected by the emergence of mutant SARS-CoV-2 variants, which are potentially more virulent and have increased infectivity. It has been suggested that the currently authorized COVID-19 vaccines may provide some protection against the various strains, including the Alpha variant [8]. However, reduced antibody neutralization and efficacy have been observed for the Delta variant.
Specifically, there is no evidence that the efficacy of Moderna or Pfizer COVID-19 vaccines is reduced against the Alpha variant [8]. In the UK, when the Alpha variant was the predominant circulating strain, two doses of the Pfizer COVID-19 vaccine demonstrated an efficacy of 85–86% against SARS-CoV-2 infection and symptomatic COVID-19 [2]. Novavax announced that its vaccine showed 85% efficacy against the Alpha variant versus 89% with the non-variant strain, while the AstraZeneca-Oxford vaccine reportedly offered 74% efficacy against the Alpha variant and 84% against the non-variant strain [6]. However, the efficacy of vaccines against the Delta variant is potentially reduced. The clinical trials of the Pfizer and AstraZeneca-Oxford COVID-19 vaccines showed reduced sera neutralization for the Delta variant, with a protective efficacy of 79% and 60% against this variant, respectively [8,9]. However, the protective efficacy of the Pfizer and AstraZeneca-Oxford vaccine was lower with a single dose, with 36% and 30% efficacy against symptomatic diseases by the Delta variant, respectively [10].
As SARS-CoV-2 variants continue to spread globally, it is imperative to assess their impact on vaccination programs at the population level. In this study, we constructed an age-structured, two-strain model of SARS-CoV-2 transmission and vaccination, and applied our model to the South Korean population, where domestic transmission of the Delta variant has been increasingly reported. We examined the impact of a novel variant on vaccination campaigns at the population level during the fourth wave of the pandemic, under various scenarios related to the transmission potential and the emergence of novel variants. Based on our simulation results, we examined whether herd immunity could be achieved with current vaccination programs in the presence of the Delta variant.

2. Materials and Methods

2.1. Model Structure

We developed an age-structured model of COVID-19 transmission and vaccination considering two co-circulating strains (Figure 1). Our model assumed that the mutant strain would be more transmissible on contact with an infectious person than the pre-existing, wild-type strain, i.e., υ > 1 . The model classified individuals according to their vaccination status and the natural history of COVID-19, considering the two strains. Additionally, as unvaccinated individuals are susceptible to both the wild and mutant strains, they were classified as follows: susceptible ( S k ); exposed ( E j , k ), asymptomatic ( A j , k ); pre-symptomatic ( P j , k ); mildly symptomatic ( M j , k ); severely symptomatic ( I j , k ); and recovered ( R j , k ). Similarly, partially (or fully) vaccinated individuals are divided into susceptible ( S V k   or S W k ), exposed ( E V j , k or E W j , k ), asymptomatic ( A V j , k or A W j , k ), pre-symptomatic ( P V j , k or P W j , k ), mildly symptomatic ( M V j , k or M W j , k ), severely symptomatic ( I V j , k or I W j , k ), and recovered ( R V j , k or R W j , k ) categories. Here, the subscript k refers to the age group of individuals (k = 1, 2, …, 8), while the subscript j indicates the original strain (j = 1) or variant alpha (j = 2). Model population was stratified into eight age groups of 0–9, 10–19, 20–29, 30–39, 40–49, 50–59, 60–69, and 70+ years based on the Korean census data.
Vaccinated individuals followed the analogous natural history of COVID-19, although they had a lower probability of both clinical and subclinical infection. However, it was assumed that vaccinated individuals who develop clinical or subclinical infection would be as infectious as unvaccinated individuals with clinical or subclinical infection. We parameterized our model with daily contacts within and between age groups using synthetic country-specific contact matrices that were constructed based on the POLYMOD study [11,12]. Our simulations were run with and without the introduction of SARS-CoV-2 variants. For cases with SARS-CoV-2 variants, we assumed that Delta would be introduced on day 25 (i.e., 30 July 2021), when school holidays and summer vacation would begin in South Korea. Our mathematical model with two strains of SARS-CoV-2 and vaccination is given by the following equations:
S k = ( θ k ( t ) + λ 1 , k ( t ) + λ 2 , k ( t ) )   S k , E 1 , k = λ 1 , k ( t ) S k γ E E 1 , k , A 1 , k = ( 1 h k ) γ E E 1 , k γ A A 1 , k , P 1 , k = h k γ E E 1 , k γ P P 1 , k , M 1 , k = ρ k γ P P 1 , k γ M M 1 , k , I 1 , k = ( 1 ρ k ) γ P P 1 , k γ I I 1 , k , R 1 , k = γ A A 1 , k + γ M M 1 , k + γ I I 1 , k σ λ 2 , k ( t ) R 1 , E 2 , k = λ 2 , k S k + σ λ 2 , k ( t ) R 1 γ E E 2 , k , A 2 , k = ( 1 h k ) γ E E 2 , k γ A A 2 , k , P 2 , k = h k γ E E 2 , k γ P P 2 , k , M 2 , k = ρ k γ P P 2 , k γ M M 2 , k , I 2 , k = ( 1 ρ k ) γ P P 2 , k γ I I 2 , k , R 2 , k = γ A A 2 , k + γ M M 2 , k + γ I I 2 , k , S V k = θ k ( t ) S k { ( 1 η 1 , k ) λ 1 , k ( t ) + ( 1 η 2 , k ) λ 2 , k ( t ) } S V k ψ k ( t )   S V k , E V 1 , k = ( 1 η 1 , k ) λ 1 , k S V k { ( 1 h k ) + ( 1 δ 1 , k ) h k } γ E E V 1 , k , A V 1 , k = ( 1 h k ) γ E E V 1 , k γ A A V 1 , k , P V 1 , k = ( 1 δ 1 , k ) h k γ E E V 1 , K γ P P V 1 , K , M V 1 , k = ρ k γ P P V 1 , k γ M M V 1 , k , I V 1 , k = ( 1 ρ k ) γ P P V 1 , k γ I I V 1 , k , R V 1 , k = γ A A V 1 , k + γ M M V 1 , k + γ I I V 1 , k , E V 2 , k = ( 1 η 2 , k ) λ 2 , k ( t ) S V k { ( 1 h k ) + ( 1 δ 1 , k ) h k } γ E E V 2 , k , A V 2 , k = ( 1 h k ) γ E E V 2 , k γ A A V 2 , k , P V 2 , k = ( 1 δ 2 , k ) h k γ E E V 2 , k γ P P V 2 , k , M V 2 , k = ρ k γ P P V 2 , k γ M M V 2 , k , I V 2 , k = ( 1 ρ k ) γ P P V 2 , k γ I I V 2 , k , R V 2 , k = γ A A V 2 , k + γ I I V 2 , k S W k = ψ k ( t )   S V k { ( 1 α 1 , k ) λ 1 , k ( t ) + ( 1 α 2 , k ) λ 2 , k ( t ) } S W k , E W 1 , k = ( 1 α 1 , k ) λ 1 , k S W k { ( 1 h k ) + ( 1 χ 1 , k ) h k } γ E E W 1 , k , A W 1 , k = ( 1 h k ) γ E E W 1 , k γ A A W 1 , k , P W 1 , k = ( 1 χ 1 , k ) h k γ E E W 1 , K γ P P W 1 , K , M W 1 , k = ρ k γ P P W 1 , k γ M M W 1 , k , I W 1 , k = ( 1 ρ k ) γ P P W 1 , k γ I I W 1 , k , R W 1 , k = γ A A W 1 , k + + γ M M W 1 , k + γ I I W 1 , k , E W 2 , k = ( 1 α 2 ) λ 2 , k S W k ( 1 χ 2 , k ) γ E E W 2 , k , A W 2 , k = ( 1 h k ) γ E E W 2 , k γ A A W 2 , k , P W 2 , k = ( 1 χ 2 , k ) h k γ E E W 2 , K γ P P W 2 , K , M W 2 , k = ρ k γ P P W 2 , k γ M M W 2 , k , I W 2 , k = ( 1 ρ k ) γ P P W 2 , k γ I I W 2 , k , R W 2 , k = γ A A W 2 , k + + γ M M W 2 , k + γ I I W 2 , k ,
where
λ 1 , k = β k j = 1 16 c k j ( q A ( A 1 , j + A V 1 , j + A W 1 , j ) + ( P 1 , j + P V 1 , j + P W 1 , j ) + q M ( M 1 , j + M V 1 , j + M W 1 , j ) + q I ( I 1 , j + I V 1 , j + I W 1 , j ) N j ) and   λ 2 , k = υ λ 1 , k .

2.2. Disease Dynamics

In our mathematical model, the risk of infection for people susceptible to COVID-19 depended on contact with individuals infected by pre-existing wild-type or mutant SARS-CoV-2 that could be in asymptomatic, pre-symptomatic, or symptomatic stages of the disease. Our model assumes that individuals with clinical symptoms are more infectious than asymptomatic individuals [13]. Specifically, it was assumed that the infectivity of asymptomatic, mild symptomatic, and severely symptomatic individuals would be 26%, 22% and 44%, respectively, relative to the pre-symptomatic stage, considering the impact of self-isolation of symptomatic patients (Table 1) [14,15,16,17]. For infected individuals, the average incubation period was assumed to be 5.2 days, and an age-specific proportion ( h k ) of infected individuals were assumed to develop symptomatic disease following a pre-symptomatic stage [18,19]. The pre-symptomatic stage and infectious period were assumed to be 2.3 days and 3.2 days on average, respectively [18,19]. Those who did not develop symptomatic disease were assumed to remain asymptomatic until recovery, with an average infectious period of 5 days [20,21]. Symptomatic cases would have an age-specific probability of developing mild or severe illness. We assumed that individuals who recovered from infection by pre-existing original strain of SARS-CoV-2 were partially susceptible to infection by SARS-CoV-2 variant.

2.3. Vaccination

As of 5 July 2021, two COVID-19 vaccines (i.e., AstraZeneca-Oxford and Pfizer COVID-19 vaccines) have been widely used in South Korea: a viral vector vaccine with lower efficacy against symptomatic disease and an mRNA vaccine that demonstrated high efficacy against symptomatic disease [22,23,24]. Specifically, as of 5 July 2021, 19,537,712 doses of approved COVID-19 vaccines have been administered in South Korea, including 11,366,200 (58.2%) doses of the AstraZeneca-Oxford vaccine and 7,002,955 (35.8%) doses of the Pfizer vaccine [25].
Our simulation approach is modeled after the vaccination program in South Korea (Table 2). Vaccines were modeled as a two-dose regimen with 19.5% and 10% of the population being partially and fully vaccinated, respectively, on day 1 (5 July 2021) [25]. In baseline vaccination scenarios, vaccination was assumed to continue according to the current rollout plan, and we examined its potential impact on symptomatic disease, hospitalizations, and deaths. Vaccine administration in South Korea is partly constrained due to the limited vaccine supply and the capacity to administer doses. In addition, vaccination against COVID-19 is performed in order of decreasing age, with the minimum age being 18. In our model, vaccines were administered to all adults aged ≥20 years in descending order of age group. As of 5 July 2021, 84% of individuals aged ≥60 years received at least one dose of COVID-19 vaccine. Thus, in our model, unvaccinated individuals aged ≥60 years were vaccinated with the first dose until the coverage level of 87% was reached, following which vaccines were administered to individuals aged 40–59 years and then to those aged 20–39 years starting on August 1 and September 1, respectively. All susceptible individuals were assumed to be eligible for vaccination. The detailed supply schedule used in this analysis was based on public announcements [26]. Specifically, the constant daily rate of vaccine administration in younger (aged 20–39) and older adults (aged 40–59) was scaled to meet the target vaccine coverage levels of 86% and 76% by day 90 (25 September), respectively. With this rollout plan, we predicted to achieve 70% coverage level of the entire population by 25 September 2021. This rate corresponds to the goal of vaccination of 36 million individuals with the first dose by the end of September 2021, as outlined by the KDCA [26]. In addition, we assumed that the vaccination campaign with the second dose would start on day 95, 110, and 120 for individuals aged ≥60, 40–59 years, and 20–39 years, respectively. The constant daily rate of vaccine administration was scaled to achieve 65% of second-dose coverage level of the entire population by 30 November 2021. Infection dynamics continued during the simulation, and outcomes were evaluated for 150 days.
We assumed “leaky” vaccine protection in which all vaccinated individuals were partially protected against infection and subject to some residual risk of infection and symptomatic disease. This contrasts with an “all-or-none” vaccine efficacy, in which a proportion of vaccinated individuals would have complete protection. The baseline values of vaccine efficacy used in the model are listed in Table 1. The overall vaccine effectiveness against disease (denoted by VE) was modeled as a function of the risk of infection and risk of symptomatic disease, conditional to infection [27]. Specifically, the vaccine effectiveness ( V E j , k ) against strain j (j = 1 for wild-type and j = 2 for SARS-CoV-2 variant) among age group k is represented as
V E j , k = 1 ( 1 η j , k ) ( 1 δ j , k )
where η j , k represents vaccine efficacy against infection and δ j , k represents the efficacy against symptomatic disease, conditional to infection. First-dose and second-dose effectiveness values (Table 1) were based on estimates from prior studies where extended dose interval strategy was employed [28,29]. We assumed that vaccine effectiveness after the first and the second dose in preventing symptomatic disease due to infection with the wild type, V E 1 , k (k = 3, …, 6), was 70% and 80%, respectively, reflecting the efficacy of the AstraZeneca or Pfizer vaccine [22,23,24,28,29,30]. An additional scenario with the emergence of the Delta variant was evaluated, considering the efficacy after receiving one dose and two doses of the AstraZeneca or Pfizer vaccine (Table 1) [22,24]. Based on a prior study, we assumed that the vaccine efficacy against the Delta variant, V E 2 , k (k = 3, …, 6), after receiving a single dose and two doses of the AstraZeneca or Pfizer vaccine is 32% and 65%, respectively (Table 1). For all scenarios, vaccine efficacy was assumed to be reduced by a factor of 30% in vaccinated individuals aged ≥60 years, based on observed reduction in vaccine effectiveness among older individuals [31]. In the absence of data for vaccine efficacy against infection or transmission, we assumed that vaccination provides protection against symptoms of both the pre-existing original strain and the Delta variant with 50% lower efficacy than its efficacy against infection (Table 1). We assumed that the vaccine efficacy would not wane during the relatively short time span used for this analysis (150 days).
Table 1. Description of the model parameters and their estimates.
Table 1. Description of the model parameters and their estimates.
NotationDescriptionEstimatesReferences
1 / γ E Incubation period (days)5.2 days[18]
1 / γ A Asymptomatic period (days)5.0 days[16,20,21]
1 / γ P Pre-symptomatic period (days)2.3 days[14,16,19]
1 / γ M Infectious period from onset of mild symptoms (days)3.2 days[16,20]
1 / γ S Infectious period from onset of severe symptoms (days)3.2 days[16,20]
1 σ Reduction of infection by variant strain due to cross immunity0.7[32]
υ Relative transmission potential of strain 2 compared to strain 12[10]
q A Relative infectivity of asymptomatic individuals compared to the pre-symptomatic stage0.26[15]
q M Relative infectivity of mildly symptomatic individuals compared to the pre-symptomatic stage, considering isolation0.22[14,17,33]
q I Relative infectivity of severely symptomatic individuals compared to the pre-symptomatic stage, considering isolation0.44[14,17,33]
h k Proportion of infections that are symptomatic0.29 for k = 1; 0.21 for
k = 2;
0.27 for k = 3; 0.33 for
k = 4;
0.40 for k = 5; 0.49 for
k = 6;
0.63 for k = 7; 0.69 for
k = 8
[13]
ρ k Proportion of symptomatic cases that exhibit mild symptoms0.90 for k = 1 and 2;
0.85 for k = 3, 4, and 5;
0.60 for k = 6 and 7;
0.20 for k = 8
[34,35]
η 1 , k Vaccine efficacy, after the first dose to before the second dose, against infection by strain 1 among the age group k0.53 for k = 1, 2, …, 6;
0.37 for k = 7 and 8
[22,23,24,31]
η 2 , k Vaccine efficacy, after the first dose to before the second dose, against infection by strain 2 among the age group k0.21 for k = 1, 2, …, 6;
0.14 for k = 7 and 8
[10,29,36,37]
δ 1 , k Vaccine efficacy, after the first dose to before the second dose, against symptomatic disease by strain 1 among the age group k0.58 for k = 1, 2, …, 6;
0.41 for k = 7 and 8
[22,23,24,31]
δ 2 , k Vaccine efficacy, after the first dose to before the second dose, against symptomatic disease by strain 2 among the age group k0.14 for k = 1, 2, …, 6;
0.10 for k = 7 and 8
[10,29,36,37]
α 1 , k Vaccine efficacy, after the second dose, against infection by strain 1 among the age group k0.64 for k = 1, 2, …, 6;
0.37 for k = 7 and 8
[38,39]
α 2 , k Vaccine efficacy, after the second dose, against infection by strain 2 among the age group k0.48 for k = 1, 2, …, 6;
0.33 for k = 7 and 8
[8,9]
χ 1 , k Vaccine efficacy, after the second dose, against symptomatic disease by strain 1 among the age group k0.37 for k = 1, 2, …, 6;
0.26 for k = 7 and 8
[22,23,24,31]
χ 2 , k Vaccine efficacy, after the second dose, against symptomatic disease by strain 2 among the age group k0.33 for k = 1, 2, …, 6;
0.23 for k = 7 and 8
[8,9]

2.4. Model Implementation

For model calibration, an effective reproduction number of 1.2 was used to account for the effect of current non-pharmaceutical COVID-19 interventions in South Korea [25]. To incorporate the age distribution of pre-existing immunity in the population, the information on case distribution as well as vaccinated proportions in age groups were used as the initial conditions (Table 2). We simulated the model across a time span of 150 days, assuming pre-existing natural and vaccine-induced immunity at the start of vaccination.
Table 2. Cumulative percentage of people who have received a COVID-19 vaccine in South Korea by vaccination status and age group as of 5 July 2021 [25].
Table 2. Cumulative percentage of people who have received a COVID-19 vaccine in South Korea by vaccination status and age group as of 5 July 2021 [25].
Partially Vaccinated (%)Fully Vaccinated (%)
Total19.510.4
Age group20–296.44.1
30–3916.414.4
40–499.04.7
50–598.63.6
60–6979.93.3
70 and above33.651.1

3. Results

The transmission probability per contact was calibrated to an effective reproduction number Re = 1.2. For the base-case scenario without the emergence of a SARS-CoV-2 variant, the attack rate was projected to be 6.9% on day 150 without further vaccination.

3.1. Attack Rates

Without the emergence of a SARS-CoV-2 variant, we found that the vaccination program against COVID-19 in South Korea would avert 45% of symptomatic infections with a baseline scenario (Figure 2). The highest relative reduction (55–56%) was observed among individuals aged 40–59 years. We found that with an increased transmission potential, i.e., Re = 1.3, the vaccination program would reduce the attack rate from 9.2% to 4.9%, averting 47% of symptomatic infections.
We found that, in the absence of vaccination and the emergence of the Delta variant, the attack rate would be 26.9% over 150 days (Figure 3), compared to 6.9% in the absence of a SARS-CoV-2 variant. Furthermore, the variant would become a dominant strain within 65 days after introduction. With the emergence of a SARS-CoV-2 variant in at the end of July, assuming a baseline transmission potential of the Delta variant, the vaccination program in South Korea would reduce cumulative symptomatic infections by 56% over the same period. As vaccine coverage level increases over time, the larger benefits of a vaccination program in averting symptomatic infections would be achieved especially in the presence a SARS-CoV-2 variant; however, herd immunity is unlikely to be achieved.

3.2. Hospitalization and Deaths

In the absence of vaccination, total intensive care unit (ICU) hospitalizations and deaths were projected to be 12.4 and 5.1 per 10,000 persons, respectively, on day 150 without the emergence of variant strains. Vaccination with reduced efficacy in elderly individuals still markedly reduced hospitalizations and deaths (Figure 4). Vaccination reduced infection outcomes, such as ICU hospitalizations and deaths, decreasing it by 45% and 43%, respectively. When we consider the situation after emergence of the Delta variant, further reduction in hospitalizations and deaths can be expected. In this scenario, total ICU hospitalizations and deaths were projected to be 77.6 and 35.1 per 10,000 persons, respectively, on day 150 without vaccines. Vaccination would reduce total ICU hospitalizations and deaths to 31.3 and 14.5 per 10,000 persons, respectively.

4. Discussion

Vaccination is expected to have a substantial impact in mitigating COVID-19 outbreaks [40]. However, owing to challenges such as vaccine supply and rollout, coupled with the potential emergence of highly transmissible SARS-CoV-2 variants, it is debatable whether the current vaccination strategy would generate herd immunity. It has been also suggested that available vaccines should be used for a single-dose COVID-19 vaccination strategy rather than a two-dose regimen, to double vaccine coverage. A few countries, including the UK and Canada, have approved guidelines to defer the second dose by up to 12 and 16 weeks, respectively [16]. The South Korean government announced plans to defer the second dose by up to 12 weeks, and aimed to vaccinate 36 million people with their first doses by September 2021. In this study, we evaluated the impact of COVID-19 vaccination program in South Korea considering the potential emergence of SARS-CoV-2 variants and its impact on the vaccination program.
Our results indicate that a SARS-CoV-2 variant with higher transmissibility compared to the currently circulating strain would become dominant within two months of being introduced. Our findings are consistent with prior studies that indicated that the variant alpha would become the predominant variant in the US, 2–3 months after introduction [3,41]. Even in the presence of SARS-CoV-2 variants, a vaccination program was shown to have potential in significantly reducing the disease burden associated with COVID-19. We project that fully vaccinating 33 million people against COVID-19 would avert at least 45% of symptomatic infections, 45% of hospitalizations, and 43% of deaths. Nevertheless, we found that herd immunity is unlikely to be achieved by the end of November 2021 with the emergence of the Delta variant, in contrast to the country’s vaccination blueprint with a goal to achieve herd immunity by November 2021.
Our study has some limitations. First, we did not explicitly model high-risk individuals including healthcare workers and comorbid individuals. Inclusion of those individuals may improve the accuracy of model predictions. Second, our parameter values for transmissibility of SARS-CoV-2 variants and the vaccine efficacy against it are based on early reports of the Delta variant in the UK; however, these estimates may vary across different populations or settings.

5. Conclusions

Our study unequivocally suggests a substantial impact of vaccination program on mitigating the severity of COVID-19 infection at the population level in South Korea. However, this strategy is unlikely to help achieve herd immunity with the potential emergence of more contagious SARS-CoV-2 variants, as the current pace of vaccine rollout is insufficient to prevent the exacerbation of cases, hospitalizations, and deaths expected.

Funding

This work was supported by the National Research Foundation of Korea grant funded by the Korean government (MSIT) [No. 2018R1C1B6001723].

Institutional Review Board Statement

Ethical review and approval were waived for this study as only aggregated publicly available data were used.

Informed Consent Statement

Not applicable.

Data Availability Statement

The daily number of confirmed cases and deaths associated with COVID-19 in South Korea was obtained from publicly available sources, available at https://www.kdca.go.kr (accessed on 5 July 2021).

Conflicts of Interest

The author does not have commercial or other association that might pose a conflict of interest.

References

  1. WHO COVID-19 Weekly Epidemiological Update. Available online: https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports (accessed on 5 July 2021).
  2. Christie, A.; Mbaeyi, S.A.; Walensky, R.P. CDC Interim Recommendations for Fully Vaccinated People: An Important First Step. JAMA J. Am. Med. Assoc. 2021, 325, 1501–1502. [Google Scholar] [CrossRef]
  3. Galloway, S.E.; Paul, P.; MacCannell, D.R.; Johansson, M.A.; Brooks, J.T.; MacNeil, A.; Slayton, R.B.; Tong, S.; Silk, B.J.; Armstrong, G.L. Emergence of SARS-CoV-2 B.1.1.7 Lineage—United States, December 29, 2020-January 12, 2021. Mmwr. Morb. Mortal. Wkly. Rep. 2021, 70, 95–99. [Google Scholar] [CrossRef]
  4. Davies, N.G.; Abbott, S.; Barnard, R.C.; Jarvis, C.I.; Kucharski, A.J.; Munday, J.D.; Pearson, C.A.B.; Russell, T.W.; Tully, D.C.; Washburne, A.D. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Science 2021, 372, eabg3055. [Google Scholar] [CrossRef]
  5. Challen, R.; Brooks-Pollock, E.; Read, J.M.; Dyson, L.; Tsaneva-Atanasova, K.; Danon, L. Risk of mortality in patients infected with SARS-CoV-2 variant of concern 202012/1: Matched cohort study. BMJ 2021, 372, n579. [Google Scholar] [CrossRef]
  6. Gómez, C.E.; Perdiguero, B.; Esteban, M. Emerging SARS-CoV-2 Variants and Impact in Global Vaccination Programs against SARS-CoV-2/COVID-19. Vaccines 2021, 9, 243. [Google Scholar] [CrossRef] [PubMed]
  7. Several Local Gov’ts Waging All-out Battle Against Coronavirus Variants. Available online: https://en.yna.co.kr/view/AEN20210511008900315?section=search (accessed on 11 May 2021).
  8. Mahase, E. Covid-19: Where are we on vaccines and variants? BMJ 2021, 372, n597. [Google Scholar] [CrossRef] [PubMed]
  9. Sheikh, A.; McMenamin, J.; Taylor, B.; Robertson, C. SARS-CoV-2 Delta VOC in Scotland: Demographics, risk of hospital admission, and vaccine effectiveness. Lancet 2021, 397, 2461–2462. [Google Scholar] [CrossRef]
  10. Callaway, E. Delta coronavirus variant: Scientists brace for impact. Nature 2021, 595, 17–18. [Google Scholar] [CrossRef] [PubMed]
  11. Prem, K.; Cook, A.R.; Jit, M. Projecting social contact matrices in 152 countries using contact surveys and demographic data. PLoS Comput. Biol. 2017, 13, e1005697. [Google Scholar] [CrossRef]
  12. Mossong, J.; Hens, N.; Jit, M.; Beutels, P.; Auranen, K.; Mikolajczyk, R.; Massari, M.; Salmaso, S.; Tomba, G.S.; Wallinga, J. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008, 5, e74. [Google Scholar] [CrossRef] [PubMed]
  13. Davies, N.G.; Klepac, P.; Liu, Y.; Prem, K.; Jit, M.; Eggo, R.M. Age-dependent effects in the transmission and control of COVID-19 epidemics. Nat. Med. 2020, 26, 1205–1211. [Google Scholar] [CrossRef]
  14. Moghadas, S.M.; Fitzpatrick, M.C.; Sah, P.; Pandey, A.; Shoukat, A.; Singer, B.H.; Galvani, A.P. The implications of silent transmission for the control of COVID-19 outbreaks. Proc. Natl. Acad. Sci. USA 2020, 117, 17513–17515. [Google Scholar] [CrossRef] [PubMed]
  15. Sayampanathan, A.A.; Heng, C.S.; Pin, P.H.; Pang, J.; Leong, T.Y.; Lee, V.J. Infectivity of asymptomatic versus symptomatic COVID-19. Lancet 2021, 397, 93–94. [Google Scholar] [CrossRef]
  16. Moghadas, S.M.; Vilches, T.N.; Zhang, K.; Nourbakhsh, S.; Sah, P.; Fitzpatrick, M.C.; Galvani, A.P. Evaluation of COVID-19 vaccination strategies with a delayed second dose. PLoS Biol. 2021, 19, e3001211. [Google Scholar] [CrossRef] [PubMed]
  17. Jarvis, C.I.; Van Zandvoort, K.; Gimma, A.; Prem, K.; Klepac, P.; Rubin, G.J.; Edmunds, W.J. Quantifying the impact of physical distance measures on the transmission of COVID-19 in the UK. BMC Med. 2020, 18, 124. [Google Scholar] [CrossRef]
  18. Li, Q.; Guan, X.; Wu, P.; Wang, X.; Zhou, L.; Tong, Y.; Ren, R.; Leung, K.S.M.; Lau, E.H.Y.; Wong, J.Y. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N. Engl. J. Med. 2020, 382, 1199–1207. [Google Scholar] [CrossRef]
  19. He, X.; Lau, E.H.Y.; Wu, P.; Deng, X.; Wang, J.; Hao, X.; Lau, Y.C.; Wong, J.Y.; Guan, Y.; Tan, X. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nat. Med. 2020, 26, 672–675. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Li, R.; Pei, S.; Chen, B.; Song, Y.; Zhang, T.; Yang, W.; Shaman, J. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (COVID-19). Science 2020, 368, 489–493. [Google Scholar] [CrossRef] [Green Version]
  21. Gatto, M.; Bertuzzo, E.; Mari, L.; Miccoli, S.; Carraro, L.; Casagrandi, R.; Rinaldo, A. Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures. Proc. Natl. Acad. Sci. USA 2020, 117, 10484–10491. [Google Scholar] [CrossRef] [Green Version]
  22. Baden, L.R.; El Sahly, H.M.; Essink, B.; Kotloff, K.; Frey, S.; Novak, R.; Diemert, D.; Spector, S.A.; Rouphael, N.; Creech, C.B. Efficacy and Safety of the mRNA-1273 SARS-CoV-2 Vaccine. N. Engl. J. Med. 2020, 384, 403–416. [Google Scholar] [CrossRef] [PubMed]
  23. Polack, F.P.; Thomas, S.J.; Kitchin, N.; Absalon, J.; Gurtman, A.; Lockhart, S.; Perez, J.L.; Perez Marc, G.; Moreira, E.D.; Zerbini, C. Safety and Efficacy of the BNT162b2 mRNA Covid-19 Vaccine. N. Engl. J. Med. 2020, 383, 2603–2615. [Google Scholar] [CrossRef] [PubMed]
  24. Voysey, M.; Clemens, S.A.C.; Madhi, S.A.; Weckx, L.Y.; Folegatti, P.M.; Aley, P.K.; Angus, B.; Baillie, V.L.; Barnabas, S.L.; Bhorat, Q.E. Safety and efficacy of the ChAdOx1 nCoV-19 vaccine (AZD1222) against SARS-CoV-2: An interim analysis of four randomised controlled trials in Brazil, South Africa, and the UK. Lancet 2021, 397, 99–111. [Google Scholar] [CrossRef]
  25. KCDC. The Updates of COVID-19 in Republic of Korea. In Centers for Disease Control and Prevention Korea: 2021. Available online: http://www.kdca.go.kr/board/board.es?mid=a20501010000&bid=0015 (accessed on 5 July 2021).
  26. Kim, A. Korea to vaccinate 12 million against COVID-19 by June: Minister. The Korea Herald. 15 March 2021. [Google Scholar]
  27. Matrajt, L.; Eaton, J.; Leung, T.; Dimitrov, D.; Schiffer, J.T.; Swan, D.A.; Janes, H. Optimizing vaccine allocation for COVID-19 vaccines shows the potential role of single-dose vaccination. Nat. Commun. 2021, 12, 3449. [Google Scholar] [CrossRef]
  28. Lumley, S.F.; Rodger, G.; Constantinides, B.; Sanderson, N.; Chau, K.K.; Street, T.L.; O’Donnell, D.; Howarth, A.; Hatch, S.B.; Marsden, B.D.; et al. An observational cohort study on the incidence of SARS-CoV-2 infection and B.1.1.7 variant infection in healthcare workers by antibody and vaccination status. Clin. Infect. Dis. 2021, ciab608. [Google Scholar] [CrossRef] [PubMed]
  29. Bernal, J.; Andrews, N.; Gower, C.; Stowe, J.; Robertson, C.; Tessier, E.; Simmons, R.; Cottrell, S.; Roberts, R.; O’Doherty, M. Early Effectiveness of COVID-19 Vaccination with BNT162b2 mRNA Vaccine and ChAdOx1 Adenovirus Vector Vaccine on Symptomatic Disease, Hospitalisations and Mortality in Older Adults in England. Available online: https://www.medrxiv.org/content/10.1101/2021.03.01.21252652v1 (accessed on 30 March 2021).
  30. Iacobucci, G. Covid-19: Infections fell by 65% after first dose of AstraZeneca or Pfizer vaccine, data show. BMJ 2021, 373, n1068. [Google Scholar] [CrossRef]
  31. Soiza, R.L.; Scicluna, C.; Thomson, E.C. Efficacy and safety of COVID-19 vaccines in older people. Age Ageing 2021, 50, 279–283. [Google Scholar] [CrossRef]
  32. Kissler, S.M.; Tedijanto, C.; Goldstein, E.; Grad, Y.H.; Lipsitch, M. Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. Science 2020, 368, 860–868. [Google Scholar] [CrossRef] [PubMed]
  33. Ferretti, L.; Wymant, C.; Kendall, M.; Zhao, L.; Nurtay, A.; Abeler-Dorner, L.; Parker, M.; Bonsall, D.; Fraser, C. Quantifying SARS-CoV-2 transmission suggests epidemic control with digital contact tracing. Science 2020, 368, 6491. [Google Scholar] [CrossRef] [Green Version]
  34. Shoukat, A.; Wells, C.R.; Langley, J.M.; Singer, B.H.; Galvani, A.P.; Moghadas, S.M. Projecting demand for critical care beds during COVID-19 outbreaks in Canada. CMAJ 2020, 192, E489–E496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Moghadas, S.M.; Shoukat, A.; Fitzpatrick, M.C.; Wells, C.R.; Sah, P.; Pandey, A.; Sachs, J.D.; Wang, Z.; Meyers, L.A.; Singer, B.H. Projecting hospital utilization during the COVID-19 outbreaks in the United States. Proc. Natl. Acad. Sci. USA 2020, 117, 9122–9126. [Google Scholar] [CrossRef] [Green Version]
  36. Dagan, N.; Barda, N.; Kepten, E.; Miron, O.; Perchik, S.; Katz, M.A.; Hernan, M.A.; Lipsitch, M.; Reis, B.; Balicer, R.D. BNT162b2 mRNA Covid-19 Vaccine in a Nationwide Mass Vaccination Setting. N. Engl. J. Med. 2021, 384, 1412–1423. [Google Scholar] [CrossRef]
  37. Hall, V.J.; Foulkes, S.; Saei, A.; Andrews, N.; Oguti, B.; Charlett, A.; Wellington, E.; Stowe, J.; Gillson, N.; Atti, A.; et al. COVID-19 vaccine coverage in health-care workers in England and effectiveness of BNT162b2 mRNA vaccine against infection (SIREN): A prospective, multicentre, cohort study. Lancet 2021, 397, 1725–1735. [Google Scholar] [CrossRef]
  38. Pawlowski, C.; Lenehan, P.; Puranik, A.; Agarwal, V.; Venkatakrishnan, A.J.; Niesen, M.J.M.; O’Horo, J.C.; Virk, A.; Swift, M.D.; Badley, A.D. FDA-authorized mRNA COVID-19 vaccines are effective per real-world evidence synthesized across a multi-state health system. Med. (N. Y.) 2021. [Google Scholar] [CrossRef]
  39. Shrotri, M.; Krutikov, M.; Palmer, T.; Giddings, R.; Azmi, B.; Subbarao, S.; Fuller, C.; Irwin-Singer, A.; Davies, D.; Tut, G. Vaccine effectiveness of the first dose of ChAdOx1 nCoV-19 and BNT162b2 against SARS-CoV-2 infection in residents of long-term care facilities in England (VIVALDI): A prospective cohort study. Lancet Infect. Dis. 2021. [Google Scholar] [CrossRef]
  40. Moghadas, S.M.; Vilches, T.N.; Zhang, K.; Wells, C.R.; Shoukat, A.; Singer, B.H.; Meyers, L.A.; Neuzil, K.M.; Langley, J.M.; Fitzpatrick, M.C. The impact of vaccination on COVID-19 outbreaks in the United States. Clin. Infect. Dis. Off. Publ. Infect. Dis. Soc. Am. 2021. [Google Scholar] [CrossRef]
  41. Parpia, A.S.; Martinez, I.; El-Sayed, A.M.; Wells, C.R.; Myers, L.; Duncan, J.; Collins, J.; Fitzpatrick, M.C.; Galvani, A.P.; Pandey, A. Racial disparities in COVID-19 mortality across Michigan, United States. EClinicalMedicine 2021, 33, 100761. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Model diagram of COVID-19 transmission and vaccination with two co-circulating strains, denoted by n (n = 1 or 2). All individuals are stratified by age, although age indices have been omitted for clarity.
Figure 1. Model diagram of COVID-19 transmission and vaccination with two co-circulating strains, denoted by n (n = 1 or 2). All individuals are stratified by age, although age indices have been omitted for clarity.
Ijerph 18 07578 g001
Figure 2. Impact of vaccination program against COVID-19 during the fourth wave in the absence of SARS-CoV-2 variant. (A) Projected cumulative attack rates over 150 days with and without further vaccination. The shaded areas represent the simulation outputs with lower and higher reproductive ratios of the wild type with Re = 1.1 and 1.3, respectively. (B) Overall and age-specific relative reduction of attack rates with vaccination, as compared to the outbreak scenario without further vaccination over 150 days.
Figure 2. Impact of vaccination program against COVID-19 during the fourth wave in the absence of SARS-CoV-2 variant. (A) Projected cumulative attack rates over 150 days with and without further vaccination. The shaded areas represent the simulation outputs with lower and higher reproductive ratios of the wild type with Re = 1.1 and 1.3, respectively. (B) Overall and age-specific relative reduction of attack rates with vaccination, as compared to the outbreak scenario without further vaccination over 150 days.
Ijerph 18 07578 g002
Figure 3. Impact of vaccination program against COVID-19 during the fourth wave in the presence of SARS-CoV-2 variant. The shaded areas represent the simulation outputs with lower and higher reproductive ratios of the Delta variant with υ = 1.8 and υ = 2.2 , respectively (A) Projected cumulative attack rates over 150 days with and without further vaccination. (B) Overall and age-specific relative reduction of attack rates with vaccination, as compared to the outbreak scenario without further vaccination over 150 days.
Figure 3. Impact of vaccination program against COVID-19 during the fourth wave in the presence of SARS-CoV-2 variant. The shaded areas represent the simulation outputs with lower and higher reproductive ratios of the Delta variant with υ = 1.8 and υ = 2.2 , respectively (A) Projected cumulative attack rates over 150 days with and without further vaccination. (B) Overall and age-specific relative reduction of attack rates with vaccination, as compared to the outbreak scenario without further vaccination over 150 days.
Ijerph 18 07578 g003
Figure 4. Projected total number of (A) ICU hospitalizations and (B) deaths per 10,000 persons over 150 days.
Figure 4. Projected total number of (A) ICU hospitalizations and (B) deaths per 10,000 persons over 150 days.
Ijerph 18 07578 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shim, E. Projecting the Impact of SARS-CoV-2 Variants and the Vaccination Program on the Fourth Wave of the COVID-19 Pandemic in South Korea. Int. J. Environ. Res. Public Health 2021, 18, 7578. https://doi.org/10.3390/ijerph18147578

AMA Style

Shim E. Projecting the Impact of SARS-CoV-2 Variants and the Vaccination Program on the Fourth Wave of the COVID-19 Pandemic in South Korea. International Journal of Environmental Research and Public Health. 2021; 18(14):7578. https://doi.org/10.3390/ijerph18147578

Chicago/Turabian Style

Shim, Eunha. 2021. "Projecting the Impact of SARS-CoV-2 Variants and the Vaccination Program on the Fourth Wave of the COVID-19 Pandemic in South Korea" International Journal of Environmental Research and Public Health 18, no. 14: 7578. https://doi.org/10.3390/ijerph18147578

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