Young-Ran Ha
Hyun-Jeong Kim1
Yoon-Seok Chung- 1Jeju Branch Office, Honam Regional Center for Disease Control and Prevention, Korea Disease Control and Prevention Agency, Jeju, Republic of Korea
- 2Jeju National Quarantine Station, Korea Disease Control and Prevention Agency, Jeju, Republic of Korea
- 3Division of Infectious Disease Diagnosis Control, Honam Regional Center for Disease Control and Prevention, Korea Disease Control and Prevention Agency, Gwangju, Republic of Korea
Introduction: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has been fueled by new variants emerging from circulating strains. Here, we report results from a genomic surveillance study of SARS-CoV-2 on Jeju Island, Republic of Korea, from February 2021 to September 2022.
Methods: A total of 3,585 SARS-CoV-2 positive samples were analyzed by Sanger sequencing of the gene encoding the spike protein before performing phylogenetic analyses.
Results: We found that the Alpha variant (B.1.1.7) was dominant in May 2021 before being replaced by the Delta variant (B.1.617.2) in July 2021, which was dominant until December 2021 before being replaced by the Omicron variant. Mutations in the spike protein, including N440K and G446S, have been proposed to contribute to immune evasion, accelerating the spread of Omicron variants.
Discussion: Our results from Juju Island, Republic of Korea, are consistent with and contribute to global surveillance efforts crucial for identifying new variants of concern of SARS-CoV-2 and for monitoring the transmission dynamics and characteristics of known strains.
Introduction
The coronavirus disease 2019 (COVID-19) is caused by numerous genetic variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that have emerged throughout the world (Kumar et al., 2021). The rapid genetic changes of SARS-CoV-2 have been influenced by various factors, including low levels of existing immunity to each new strain and the high mutation rates caused by error-prone polymerases (Denison et al., 2011). Screening for changes in genetic characterization could help predict vaccine effectiveness and vaccine-induced immune evasion to predict outbreaks (Mengist et al., 2021; Shaibu et al., 2021).
The SARS-CoV-2 genome encodes four structural proteins: nucleoprotein (NP), membrane (M), envelope (E), and spike (S) proteins (Rathnasinghe et al., 2021). SARS-CoV-2 invades host cells via direct contact with the human angiotensin-converting enzyme 2 (hACE2) receptor and the receptor binding domain (RBD) on the spike protein (Walls et al., 2020; Chen et al., 2022; Rashid and Salih, 2022). Spike (S) protein is a surface glycoprotein comprised of two functional subunits (i.e., S1 and S2). The S1 subunit is composed of an N-terminal domain and an RBD. The S2 subunit contains four domains: (1) a fusion peptide (FP) that mediates the binding of viral and host cell membranes; (2) two heptad repeat regions (HR1 and HR2) important for viral fusion; (3) a transmembrane (TM) domain that is role for membrane fusion activity; and (4) a cytoplasmic (CT) domain (Bosch et al., 2003; Corver et al., 2009; Guruprasad, 2021; Paul et al., 2021).
Upon binding to receptors on the cell surface, the spike protein is activated by TM protease serine 2 (TMPRSS2), promoting the entry of the virus into cells and making it an important target for multiple vaccines (Hoffmann et al., 2020). The protein has also been monitoring mutations, including insertions, deletions, and amino acid substitutions (Walls et al., 2020).
The World Health Organization (WHO) has been monitoring changes to the SARS-CoV-2 genome (World Health Organization [WHO], 2021), classifying emerging variants as variants of concern (VOC), variants of interest (VOI), and variants under monitoring (VUM) (World Health Organization [WHO], 2021). Since late November 2021, the Omicron (B.1.1.529) as a VOC has dominated global spread, replacing previously circulating VOCs, notably Alpha (B.1.1.7), Beta (B.1.531), and Gamma (P.1) by March 2022, and Delta (B.1.627.2) by July 2022 (World Health Organization [WHO], 2021). Since February 2022, the WHO began introduced a fourth variant category called “Omicron subvariants under monitoring,” which included the Omicron subvariants BA.5, BA.2.75, BA.4.6, XBB, and BA.2.3.20 by 29 November 2022 (World Health Organization [WHO], 2021). To be classified as an Omicron subvariants under monitoring, the strain must: (i) belong to a currently circulating VOC by phylogenetic analysis; (ii) be more transmissible than other circulating VOC lineages; and (iii) have amino acid substitutions that confer fitness advantages relative to other circulating variants, in epidemiology studies (World Health Organization [WHO], 2021).
In this study, we sequenced and analyzed 3,585 SARS-CoV-2 positive specimens isolated from COVID-19 positive patients on Jeju Island, Republic of Korea, from January 2021 to September 2022. We analyzed spike protein gene sequences and performed phylogenetic analyses that allowed us to follow the dynamics of Alpha, Beta, Delta, and Omicron variants in the region and report mutations related to known mechanisms of immune evasion.
Materials and methods
Epidemiological data analysis
The map of South Korea was obtained from the Statistical Geographic Information Service (SIGS) offered by Statistics Korea1 (Figures 1A, B). The number of confirmed cases of COVID-19 was provided by Korea Disease Control and Prevention Agency (KDCA2) (Figure 1C).
Figure 1. Map of South Korea showing the number of confirmed cases of COVID-19 from January 2020 to September 2022 (A). Jeju Island including Jeju-Si and Seogwipo-Si (B). The monthly number of confirmed cases of SARS-CoV-2 on Jeju Island from January 2020 to September 2022 (C).
RNA isolation and reverse transcription-polymerase chain reaction
Coronavirus disease 2019 (COVID-19) positive specimens were obtained from Seoul Clinical Laboratories (SCL) and from the Jeju Special Self-Governing Province Institute of Environment Research. Specimens were handled in a Class II biosafety cabinet (Thermo Scientific 1300 Series A2) in a biosafety level 2 (BL2) laboratory. RNA was isolated using an RNA extraction kit (Qiagen, Valencia, CA, USA), as per the manufacturer’s instructions. cDNA was synthesized using SuperScript IV First-Strand Synthesis System (Invitrogen, Carlsbad, CA, USA), and PCR was conducted using PrimeSTAR GXL DNA Polymerase (Takara, Shiga, Japan) from January to August 2021. From September 2021 to September 2022, reverse transcription-polymerase chain reaction (RT-PCR) was conducted using Qiagen Onestep RT-PCR kit (Qiagen, Valencia, CA, USA), DiaStar 2 × OneStep RT-PCR premix kit (SolGent, Daejeon, Republic of Korea), and SEQMAX qPCR one step mastermix (Nine Korea, Republic of Korea). Sequences of primers used in RT-PCR are shown in Table 1. To detect mutations of SARS-CoV-2 spike protein, L71/R75, L76/R79, and L80/R84 primer sets were used. RT-PCR was performed under the following conditions: an initial reverse transcription step at 50°C for 30 min followed by a denaturation step at 95°C for 5–15 min. This was followed by 35 cycles of 30 s at 95°C, 30 s at 58°C, 1 min 30 s at 68°C, and a final extension step at 68°C for 7 min (Ha et al., 2022). To distinguish between BA.4 and BA.5, RT-PCR was instead performed on the N (nucleocapsid) gene under the following conditions: an initial reverse transcription step at 50°C for 30 min followed by a denaturation step at 95°C for 5–15 min. This was followed by 35 cycles of 10 s at 98°C, 15 s at 50°C, 2 min 30 s at 68°C, and a final extension step at 68°C for 7 min.
Table 1. List of primer sequences used for RT-PCR analyses of spike protein (Fragment 1, 2, and 3) and N protein (L94 and R95).
DNA sequencing and analyses
DNA sequencing (Cosmogenetech, Seoul, Republic of Korea and SolGent, Daejeon, Republic of Korea) was performed using a standard protocol. To analyze the sequences for spike protein, L71/R75, L76/R79, and L80/R84 and inner primers including L73/R73, L78/R77, and L83/R82 were used until December 2021 (Table 1). From 2022, inner primers were used only L73, L78, and L83. To detect the mutation in the N gene, L94/R95 were used from July and September 2021. The sequence data from COVID-19 positive specimens on Jeju Island were submitted to GenBank through the National Center for Biotechnology Information (NCBI) (Table 2). Full-length sequences of COVID-19 spike proteins of various geographical origins were downloaded in FASTA format from GISAID. The sequence of SARS-CoV-2 reference genome Wuhan Hu-1 (accession number: NC 045512.2) was retrieved from the NCBI databases. The non-coding 3′ and 5′ regions were trimmed using CLC Genomic Workbench 5.0.1 software (CLC bio, Denmark). Multiple sequence alignments were done using a multiple sequence alignment (MAFFT) programs. Sequences of SARS-CoV-2 with specific mutations were searched in the GISAID database.3 Phylogenetic analysis was done using iTOL.4 Branch support was calculated by bootstrap, consisting of 1,000 alignments. SARS-CoV-2 positive specimens from confirmed cases of Jeju Island were classified into clades and lineages using Nextclade Beta.5
Table 2. Spike protein mutations from SARS-CoV-2 confirmed cases of Jeju Island including Alpha, Beta, Delta, and Omicron.
Results
Epidemiological data analysis
Jeju Island is Korean Republic’s biggest Island, located 80 km away from the mainland (Jang et al., 2019) and includes Jeju-Si and Seogwipo-Si (Figures 1A, B). The number of confirmed cases of COVID-19 was shown from January 2020 to September 2022 (Figure 1A). The cumulative number of COVID-19 cases in Jeju was 333,722, accounting for 1.35% of the total number of cases in the Republic of Korea (Figure 1A). From January 2020 to September 2022, the monthly number of confirmed cases of SARS-CoV-2 on Jeju Island was shown in Figure 1C. The Delta (B.1.627.2) variant emerged and began to spread on the island in May 2021, with a sharp increase in confirmed cases appearing later (e.g., from July to August 2021) (Figure 1C). The Omicron (B.1.1.529) variant was first detected in December 2021 and confirmed cases rose beginning shortly after (e.g., January 2022) (Figure 1C).
Distribution of clades
Of the 3,585 spike protein gene sequences analyzed in our study, 3,322 (92.66%) were local and 263 (7.34%) were imported from off-island (Supplementary Figure 1). SARS-CoV-2 variants were classified into lineages based on the sequences of their spike protein genes (Figure 2). From April to September 2021, strains from clade 20A (i.e., B.1.619, B.1.619.1, and B.1.620) were also observed. B.1.619 lineage was first observed in April 2021 and 20.00, 32.43, and 31.52% from May to July 2021. B.1.620 lineage also observed in April 2021 and 4.83 and 18.92% in May and June 2021 (Supplementary Figure 2). The prevalence of Epsilon (B.1.429), one of the previous circulating VOIs, was 4.35 and 12.24% in March and April 2021, respectively (Figure 2).
Figure 2. Changes in the genomic prevalence of SARS-CoV-2 variants characterized by spike protein mutations from February 2021 to September 2022 across Jeju Island.
Alpha variant (B.1.1.7) has several amino acid exchanges and deletions within the viral spike protein including del69–70HV, del144Y, N501Y, A570D, D614G, P681H, T761I, S982A, and D1118H (Socher et al., 2021). In this study, Alpha variant (OP763726∼OP763743) has an additional spike S98F and D138H mutations, which are present in 1.90 and 1.14% of strains deposited in the GISAID database, as of May 2023 (Table 2). Previous study also noted that S98F, D138H mutation was found in in less than 2% of genomes worldwide (Zarate et al., 2022). The S98F and D138H mutations are also found in Alpha variant such as hCoV-19/Spain/AN-CBA-06831/2021 (EPI_ISL_15678545), hCoV-19/Germany/BW-RKI-I-1053239/2021 (EPI_ISL_16147015), hCoV-19/USA/PA-CDC-4197450/2021 (EPI_ISL_17534404), and hCoV-19/Chile/AP-ISPC-86703/2021 (EPI_ISL_17028381).
The previously circulating VOCs, Beta variant (B.1.351) was detected infrequently, its prevalence decreasing from 4.35 to 0.69% between March and May 2021 (Figure 2). Beta variant contains amino acid mutations within the viral spike protein including L18F, D80A, D215G, deletion at positions 242–244 (L242del, A243del, and L244del), K417N, E484K, N501Y, A701V, A27S, and D614G (Edara et al., 2021; Sanches et al., 2021; Velasco et al., 2021). In this study, Beta variant (OP763744∼OP763745) has an additional spike A67V and H49Y mutations, which are present in 1.35 and 0.11% of strains deposited in the GISAID database, as of May 2023 (Table 2). The A67V mutations are also found in Beta variant such as hCoV-19/Canada/UN-63124/2021 (EPI_ISL_12388094), hCoV-19/Bangladesh/icddrb-GAC2-285/2021 (EPI_ISL_13477157), hCoV-19/USA/IN-CDC-LC26926/2021 (EPI_ISL_17532096), and hCoV-19/Spain/GA-CHUVI-19298619c-2/2021 (EPI_ISL_176218 27). The H49Y mutations are also found in Beta variant such as hCoV-19/South_Africa/NICD-K253-002-172/2021 (EPI_ISL_1607 8823), hCoV-19/South_Korea/KDCA2733s/2021 (EPI_ISL_27120 10), hCoV-19/USA/NY-CDC-QDX24758501/2021 (EPI_ISL_436 7140), and hCoV-19/Philippines/PH-PGC-43487/2021 (EPI_ ISL_5569124).
The Delta variant was first observed in May 2021 and began to dominate from July to December 2021, while Omicron emerged and replaced previous strains starting in January 2022 (Figure 2). The spike gene mutations in Delta variant are T19R, L452R, T478K, D614G, P681R, E156G, G142D, and D950N, with deletions at positions 157 and 158 (Shiehzadegan et al., 2021; Dhawan et al., 2022; Focosi et al., 2022; Singh et al., 2022b). In this study, Delta variant (OP763715) has an additional spike G145D mutation, which are present in 0.01% of strains deposited in the GISAID database, as of May 2023 (Table 2). The G145D mutations was also found in Delta variant such as hCoV-19/USA/CA-HLX-STM-3XXQ3BM69/2021 (EPI_ISL_12833436), hCoV-19/India/TN-INSACOG-CBR-S-7676/2021 (EPI_ISL_14677988), and hCoV-19/Russia/MOW-PMVL-DZ-19421141784/2021 (EPI_ ISL_15884308).
The Omicron variant was first observed in the case of overseas imports in December 2021. It became dominant in 2022, before being supplanted by the BA.2 subvariant between January and April 2022 (Figures 2, 3). BA.4 and BA.5 begin to replace the BA.2 variant beginning in June 2022, while the BA.5 subvariant dominated transmission from August to September 2022 (Figure 3). Between August and September 2021, the prevalence of BA.5 increased from 27.84 to 34.05%. The spike gene mutations in Omicron variant are A67V, Δ69-70, T95I, G142D/Δ143-145, Δ211/L212I, ins214EPE, V213G, L212I, G339D, S371L, S371F, S373P, S375F, T376A, D405N, R408S, K417N, N440K, G446S, S477N, T478K, L452R/L452Q, E484A, F486V, Q493R, G496S, Q498R, N501Y, Y505H, T547K, D614G, H655Y, N679K, P681H, S704L, N764K, D796Y, N856K, A879V, Q954H, N969K, and L981F (Muik et al., 2022; Shanmugaraj et al., 2022; Singh et al., 2022a; Takashita et al., 2022; Yoshida et al., 2022). In this study, Omicron variant (OP763688 and OP763695) has an additional spike L179P and Q271R mutations, which are present in 0.02 and 0.002% of strains deposited in the GISAID database, as of May 2023 (Table 2). The Q271R mutation was also found in Omicron variant such as hCoV-19/Sweden/C-03_SE600_000200436976N2/2022 (EPI_ISL_17712759), hCoV-19/Canada/QC-L00612610001/2023 (EPI_ISL_17641597), and hCoV-19/South_Korea/GJ-HERI-A1144/2022 (EPI_ISL_17703223). The Q271R mutation was also found in Omicron variant such as hCoV-19/France/ARA-HCL-723000653001/2023 (EPI_ISL_17657463), hCoV-19/Brazil/PA-LACENPA-150377499/2023 (EPI_ISL_17514895), and hCoV-19/South_Korea/KDCA227335/2023 (EPI_ISL_17473659).
Figure 3. Prevalence of SARS-CoV-2 Omicron subvariants on Jeju Island from December 2021 to September 2022.
Phylogenetic analyses of currently and previously circulating VOCs
A phylogenetic tree of Alpha (B.1.1.7) variants isolated from patients on Jeju Island between May and early June 2021 (OP763726∼OP763643) was generated. The Alpha lineage from samples isolated on Jeju Island was more closely related to each other than the rest of Alpha lineage (Figure 4A). Similarly, we generated phylogenetic trees of samples isolated from patients on Jeju Island for Beta (B.1.531, March to May 2021, OP763744 and OP763645). Notably, the sequence from sample OP763744 closely resembled that of hCoV-19/France/IDF-ALPIGENE-2107210176 2021 (EPI ISL 11449404) (Figure 4B). The Delta variant on Jeju Island was grouped by confirmation date of COVID-19 infection (Figure 4C). The sequences from OP763699 to OP763709 emerged August and September 2021, and Nextclade analyses classified these as belonging to clade 20I. Notably, the sequences from OP763703 to OP763707 closely match that of hCoV-19/South Korea/KDCA29074/2022 (EPI ISL 10046128). The sequences from OP763711 to OP763719 emerged in November 2021, and Nextclade analyses classified these to clade 21A. We also report that in OP763708, an additional T95I mutation was found in Spike protein (Table 2). The genomes sequences of B.1.1.529 (Omicron) on Jeju Island clustered several subvariants including BA.1.1, BA.2, and BA.5 (Figure 4D).
Figure 4. Phylogenetic analysis of SARS-CoV-2 variants including Alpha (A), Beta (B), Delta (C), and Omicron (D) characterized by spike protein mutations from December 2021 to September 2022 across Jeju Island. Red parts indicate the genomes from confirm cases in Jeju Island.
The sequences from OP763667 to OP763671, OP763696, and OP763698 were emerged in March 2022. These sequences belong to BA.1.1 subvariant. The sequences from OP763672 to OP763681 were also emerged on March 2022 and belong to BA.2 subvariant. The sequences from OP763682, OP763683, OP763685, and OP763686 were emerged on March 2022 and belong to BA.2.12.1. The sequences from OP763687 to OP753695 were emerged in August 2022 and belong BA.5 subvariant (Figure 4D). The sequence from OP763672 to OP763680, OP763682, OP763683, OP763685, and OP763686 closed to hCoV-19/Norway 126 2022 (EPI ISL 16100571) and hCoV-19/Canada NL-PHML-519353 2022 (EPI ISL 16100675). The sequence of OP753695 closed to hCoV-19/South Korea/KDCA33410s/2022 (EPI ISL 16052972). The number of mutations in Spike protein of Alpha, Beta, and Delta variant ranged from 11 to 12, while the number of mutations of Omicron exceeded 30 (Table 2).
Discussion
In this study, we performed genomic surveillance of SARS-CoV-2 samples isolated from patients on Jeju Island 1 year after that original infection, from February 2021 to September 2022, a period during which multiple variants were circulating globally. During the period of our spike protein-based genomic surveillance study on Jeju Island, the Alpha, Delta, and Omicron variants (i.e., currently and previously circulating VOCs) became dominant in turn, while the Beta and Gamma variants did not. In addition, the import and prevalence of mutations in Jeju showed different from the mainland of Korea.
On Jeju Island, the Alpha variant was still dominant in May 2021 before decreasing from 56.55 to 6.52% as the Delta variant increased in circulation (Figure 2). In the Republic of Korea, Alpha variant was decreased from 27.40 to 8.1% in May and July 2021 (Kim et al., 2021). Globally, the Alpha variant dominated transmission until the emergence of the Delta variant in the winter of 2020 (Dong et al., 2022). Our data suggests that the Delta variant began to spread on overseas imported case from Jeju on May 2021, with confirmed cases of COVID-19 increasing that month from 82 to 323. After that, confirmed cases of COVID-19 on Jeju Island increased from 468 to 870 from July to August 2021 (Figure 1C and Supplementary Figure 3). In the Republic of Korea, the Delta variant was first confirmed on April 2021 and exceeded half of all cases by July 2021 (Kim et al., 2021, 2022). Globally, the Delta variant has been replaced by the Omicron variant by the end of 2021 (Lista et al., 2022), while on Jeju Island, Delta variant has been dominant over half the year, only being replaced by Omicron in January 2022 (Figure 2).
The Omicron variant was declared as a VOC on 26 November 2021, after being identified earlier that month in South Africa (Gobeil et al., 2022). In the Republic of Korea, the Omicron variant was first confirmed in the end of November 2021 and exceeded half of all cases on December 2021 (Kim et al., 2022). By February 2022, the Omicron variant accounted for 99.1% of all COVID-19 cases in the Republic of Korea (Park et al., 2023). On Jeju Island, confirmed cases of COVID-19 increased significantly in a short period, from 784 in January 2022 to 134,559 in April 2022 (Figure 1C). In addition, the Omicron variant accounted for more than 99.00% of all COVID-19 cases on February 2022 (Figure 2). From this result, the import and prevalence of mutations in Jeju are different from those in Korea.
Conclusion
In conclusion, we used genomic surveillance and phylogenetic analyses to follow the relative abundance of different SARS-CoV-2 variants, from January 2020 to February 2022, on Jeju Island, Republic of Korea. During this period, multiple variants and their subvariants emerged having enhanced human-to-human transmissibility. The dynamics of variants differed on Jeju Island compared with the mainland. Because of the limited entry points compared to land-locked regions, new variants might be introduced more slowly or take different forms due to other Island characteristics. Being a major center for domestic and international travel, Jeju Island might be at risk of continuous influx of new variants. We conclude that monitoring efforts for COVID-19 cases and SARS-CoV-2 variants is essential for monitoring and controlling outbreaks on Jeju Island.
Data availability statement
The datasets presented in this study can be found in the GenBank. The accession number from OP763667 to OP763745 can be found in the article/Supplementary material.
Ethics statement
This study uses strains obtained from Jeju Special Self-Governing Province Research Institute of Public Health and Environment and Seoul Clinical Laboratories (SCL). Korea Disease Control and Prevention Agency (KDCA) did not require the study to be reviewed or approved by an ethics committee because this study belongs to the case including (1) insignificant impact on donors of human derivatives and the public, (2) research in which researchers do not collect or record personal information, and (3) cases where human material is not directly collected. The Institutional Review Board (IRB) of Korea Disease Control and Prevention Agency (KDCA) granted exemption for this study (IRB exemption number: KDCA-2023-06-01-PE-01).
Author contributions
Y-RH contributed to the conception, design, data acquisition, and drafting of the manuscript. H-JK conducted the experiment and helped sequencing data analysis. J-SP advised the experimental procedure and treatment of specimens. Y-SC conceived the entire study and helped draft the manuscript. All authors read and approved the final manuscript.
Funding
This research was supported by a fund (6137-300-210-13) of the Korea Centers for Disease Control and Prevention Agency (KDCA).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1170766/full#supplementary-material
Footnotes
- ^ https://sgis.kostat.go.kr/
- ^ http://ncov.kdca.go.kr
- ^ https://www.gisaid.org/
- ^ https://itol.embl.de
- ^ http://clades.nextstrain.org/
References
Bosch, B. J., Van Der Zee, R., De Haan, C. A., and Rottier, P. J. (2003). The coronavirus spike protein is a class I virus fusion protein: Structural and functional characterization of the fusion core complex. J. Virol. 77, 8801–8811. doi: 10.1128/jvi.77.16.8801-8811.2003
Chen, J., Wang, R., Gilby, N. B., and Wei, G. W. (2022). Omicron variant (B.1.1.529): Infectivity, vaccine breakthrough, and antibody resistance. J. Chem. Inf. Model. 62, 412–422.
Corver, J., Broer, R., Van Kasteren, P., and Spaan, W. (2009). Mutagenesis of the transmembrane domain of the SARS coronavirus spike glycoprotein: Refinement of the requirements for SARS coronavirus cell entry. Virol. J. 6:230. doi: 10.1186/1743-422X-6-230
Denison, M. R., Graham, R. L., Donaldson, E. F., Eckerle, L. D., and Baric, R. S. (2011). Coronaviruses: An RNA proofreading machine regulates replication fidelity and diversity. RNA Biol. 8, 270–279. doi: 10.4161/rna.8.2.15013
Dhawan, M., Sharma, A., Priyanka, Thakur, N., Rajkhowa, T. K., and Choudhary, O. P. (2022). Delta variant (B.1.617.2) of SARS-CoV-2: Mutations, impact, challenges and possible solutions. Hum. Vaccin. Immunother. 18:2068883. doi: 10.1080/21645515.2022.2068883
Dong, R., Hu, T., Zhang, Y., Li, Y., and Zhou, X. H. (2022). Assessing the transmissibility of the new SARS-CoV-2 variants: From delta to omicron. Vaccines 10:496.
Edara, V. V., Norwood, C., Floyd, K., Lai, L., Davis-Gardner, M. E., Hudson, W. H., et al. (2021). Infection- and vaccine-induced antibody binding and neutralization of the B.1.351 SARS-CoV-2 variant. Cell Host Microbe 29, 516–521.e513.
Focosi, D., Maggi, F., Mcconnell, S., and Casadevall, A. (2022). Spike mutations in SARS-CoV-2 AY sublineages of the Delta variant of concern: Implications for the future of the pandemic. Future Microbiol. 17, 219–221. doi: 10.2217/fmb-2021-0286
Gobeil, S. M., Henderson, R., Stalls, V., Janowska, K., Huang, X., May, A., et al. (2022). Structural diversity of the SARS-CoV-2 Omicron spike. Mol. Cell 82, 2050–2068.e2056.
Ha, Y. R., Shin, E. S., Kim, H. J., Hyeon, E. H., Park, J. S., and Chung, Y. S. (2022). Outbreak of Severe Acute Respiratory Syndrome Coronavirus-2 B.1.620 Lineage in the General Hospital of Jeju Island, Republic of Korea. Front. Microbiol. 13:860535. doi: 10.3389/fmicb.2022.860535
Hoffmann, M., Kleine-Weber, H., Schroeder, S., Kruger, N., Herrler, T., Erichsen, S., et al. (2020). SARS-CoV-2 Cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 181, 271–280.e278.
Jang, G., Park, J., and Lee, C. (2019). Complete genome sequences of novel S-deletion variants of porcine epidemic diarrhea virus identified from a recurrent outbreak on Jeju Island, South Korea. Arch. Virol. 164, 2621–2625. doi: 10.1007/s00705-019-04360-4
Kim, I., Park, A., Lee, H., Kim, H., Kim, J., Kim, J., et al. (2022). Status and characteristics of the SARS-CoV-2 variant outbreak in the Republic of Korea in January 2021. J. Public Health Wkly. Rep. 15, 505–510.
Kim, I., Park, A., Lee, H., Kim, J., Kim, D., Kim, J., et al. (2021). July 2021 status and characteristics of the COVID-19 variant virus outbreak in the Republic of Korea. J. Public Health Wkly. Rep. 14, 3388–3396.
Kumar, A., Dowling, W. E., Roman, R. G., Chaudhari, A., Gurry, C., Le, T. T., et al. (2021). Status Report on COVID-19 Vaccines Development. Curr. Infect. Dis. Rep. 23:9.
Lista, M. J., Winstone, H., Wilson, H. D., Dyer, A., Pickering, S., Galao, R. P., et al. (2022). The P681H Mutation in the spike glycoprotein of the alpha variant of SARS-CoV-2 escapes IFITM restriction and is necessary for type I interferon resistance. J. Virol. 96:e0125022. doi: 10.1128/jvi.01250-22
Mengist, H. M., Kombe Kombe, A. J., Mekonnen, D., Abebaw, A., Getachew, M., and Jin, T. (2021). Mutations of SARS-CoV-2 spike protein: Implications on immune evasion and vaccine-induced immunity. Semin. Immunol. 55:101533. doi: 10.1016/j.smim.2021.101533
Muik, A., Lui, B. G., Bacher, M., Wallisch, A. K., Toker, A., Finlayson, A., et al. (2022). Omicron BA.2 breakthrough infection enhances cross-neutralization of BA.2.12.1 and BA.4/BA.5. Sci. Immunol. 7:eade2283.
Park, A. K., Kim, I. H., Lee, C. Y., Kim, J. A., Lee, H., Kim, H. M., et al. (2023). Rapid emergence of the omicron variant of severe acute respiratory syndrome coronavirus 2 in Korea. Ann. Lab. Med. 43, 211–213.
Paul, D., Pyne, N., and Paul, S. (2021). Mutation profile of SARS-CoV-2 spike protein and identification of potential multiple epitopes within spike protein for vaccine development against SARS-CoV-2. Virusdisease 32, 703–726. doi: 10.1007/s13337-021-00747-7
Rashid, P. M. A., and Salih, G. F. (2022). Molecular and computational analysis of spike protein of newly emerged omicron variant in comparison to the delta variant of SARS-CoV-2 in Iraq. Mol. Biol. Rep. 49, 7437–7445. doi: 10.1007/s11033-022-07545-4
Rathnasinghe, R., Jangra, S., Cupic, A., Martinez-Romero, C., Mulder, L. C. F., Kehrer, T., et al. (2021). The N501Y mutation in SARS-CoV-2 spike leads to morbidity in obese and aged mice and is neutralized by convalescent and post-vaccination human sera. Medrxiv [Preprint]. doi: 10.1101/2021.01.19.21249592
Sanches, P. R. S., Charlie-Silva, I., Braz, H. L. B., Bittar, C., Freitas Calmon, M., Rahal, P., et al. (2021). Recent advances in SARS-CoV-2 spike protein and RBD mutations comparison between new variants Alpha (B.1.1.7, United Kingdom), Beta (B.1.351, South Africa), Gamma (P.1, Brazil) and Delta (B.1.617.2, India). J. Virus Erad. 7:100054. doi: 10.1016/j.jve.2021.100054
Shaibu, J. O., Onwuamah, C. K., James, A. B., Okwuraiwe, A. P., Amoo, O. S., Salu, O. B., et al. (2021). Full length genomic sanger sequencing and phylogenetic analysis of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) in Nigeria. PLoS One 16:e0243271. doi: 10.1371/journal.pone.0243271
Shanmugaraj, B., Malla, A., Khorattanakulchai, N., and Phoolcharoen, W. (2022). SARS-CoV-2 omicron variant: Could it be another threat? J. Med. Virol. 94, 1284–1288.
Shiehzadegan, S., Alaghemand, N., Fox, M., and Venketaraman, V. (2021). Analysis of the delta variant B.1.617.2 COVID-19. Clin. Pract. 11, 778–784.
Singh, P., Sharma, K., Shaw, D., Bhargava, A., and Negi, S. S. (2022a). Mutational characterization of Omicron SARS-CoV-2 lineages circulating in Chhattisgarh, a central state of India. Front. Med. 9:1082846. doi: 10.3389/fmed.2022.1082846
Singh, P., Sharma, K., Singh, P., Bhargava, A., Negi, S. S., Sharma, P., et al. (2022b). Genomic characterization unravelling the causative role of SARS-CoV-2 Delta variant of lineage B.1.617.2 in 2nd wave of COVID-19 pandemic in Chhattisgarh India. Microb. Pathog. 164:105404. doi: 10.1016/j.micpath.2022.105404
Socher, E., Conrad, M., Heger, L., Paulsen, F., Sticht, H., Zunke, F., et al. (2021). Mutations in the B.1.1.7 SARS-CoV-2 spike protein reduce receptor-binding affinity and induce a flexible link to the fusion peptide. Biomedicines 9:525. doi: 10.3390/biomedicines9050525
Takashita, E., Yamayoshi, S., Simon, V., Van Bakel, H., Sordillo, E. M., Pekosz, A., et al. (2022). Efficacy of antibodies and antiviral drugs against omicron BA.2.12.1, BA.4, and BA.5 subvariants. N. Engl. J. Med. 387, 468–470. doi: 10.1056/NEJMc2207519
Velasco, J. M., Chinnawirotpisan, P., Valderama, M. T., Joonlasak, K., Manasatienkij, W., Huang, A., et al. (2021). Coding-complete genome sequences of 11 SARS-CoV-2 B.1.1.7 and B.1.351 variants from metro manila, Philippines. Microbiol. Resour. Announc. 10:e0049821. doi: 10.1128/MRA.00498-21
Walls, A. C., Park, Y. J., Tortorici, M. A., Wall, A., Mcguire, A. T., and Veesler, D. (2020). Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell 183:1735.
World Health Organization [WHO] (2021). Tracking SARS-CoV-2 variants [Online]. Geneva: World Health Organization.
Yoshida, A., Okamura, S., Torii, S., Komatsu, S., Miyazato, P., Sasaki, H., et al. (2022). Versatile live-attenuated SARS-CoV-2 vaccine platform applicable to variants induces protective immunity. iScience 25:105412. doi: 10.1016/j.isci.2022.105412
Keywords: COVID-19, SARS-CoV-2, spike protein, Republic of Korea, Jeju Island
Citation: Ha Y-R, Kim H-J, Park J-S and Chung Y-S (2023) Genomic surveillance of genes encoding the SARS-CoV-2 spike protein to monitor for emerging variants on Jeju Island, Republic of Korea. Front. Microbiol. 14:1170766. doi: 10.3389/fmicb.2023.1170766
Received: 21 February 2023; Accepted: 03 July 2023;
Published: 18 July 2023.
Edited by:
Hye-Ra Lee, Korea University, Republic of KoreaReviewed by:
Diego Forni, Eugenio Medea (IRCCS), ItalyUn Yung Choi, Konkuk University, Chungju, Republic of Korea
Copyright © 2023 Ha, Kim, Park and Chung. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yoon-Seok Chung, cm9sbHN0b25lOTNAa29yZWEua3I=