Molecular, Immunological, and Clinical Characterization of SARS-CoV-2 Infection in Breakthrough-Infected Healthcare Workers During the Last Wave of The Pandemic in Shiraz, Iran

authors:

avatar Sepideh Saeb ORCID 1 , avatar Gholam Reza Talei ORCID 2 , avatar Azam Ghaziasadi ORCID 3 , avatar Arash Letafati ORCID 3 , avatar Sayyad Khanizadeh ORCID 2 , avatar Marzieh Jamalidoust ORCID 4 , *

Department of Virology, School of Medicine, Lorestan University of Medical Sciences, Khorramabad, Iran
Department of Virology, Hepatitis Research Center, Medical School, Lorestan University of Medical Sciences, Khorramabad, Iran
Department of Virology, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran
Department of Virology, Alborzi Clinical Microbiology Research Center, Namazi Hospital, Shiraz University of Medical Sciences, Shiraz, Iran

how to cite: Saeb S, Talei G R, Ghaziasadi A, Letafati A, Khanizadeh S, et al. Molecular, Immunological, and Clinical Characterization of SARS-CoV-2 Infection in Breakthrough-Infected Healthcare Workers During the Last Wave of The Pandemic in Shiraz, Iran. Jundishapur J Microbiol. 2023;16(11):e142644. https://doi.org/10.5812/jjm-142644.

Abstract

Background:

The coronavirus disease 2019 (COVID-19) poses a threat to the global economy and public health. Mutations in the spike protein of the virus can impact the functional characteristics and effectiveness of vaccines.

Objectives:

This study aimed to analyze mutations in three key regions of the spike protein, investigate clinical presentations, and assess protective immune levels in healthcare workers (HCWs) who had received full vaccine doses but were re-infected.

Methods:

In this cross-sectional study, 49 breakthrough-infected HCWs were included from November 2021 to May 2022. After confirming COVID-19 reinfection, the SARS-CoV-2 genome was extracted using the ROJE extraction kit, and genetic variation analysis was performed through Sanger sequencing. Blood samples were collected with prior consent, and ELISA tests were performed to determine antibody levels. Clinical presentations were recorded, and independent t-tests revealed no significant gender-based differences in the C-terminal point mutation.

Results:

In this study, 75.5% (37/49) of the included HCWs were female (P-value < 0.05), and 25 were qualified for PCR and gene sequencing. Mutations were observed in 25, 10, and 1 sequence(s) of the C-terminal domain, N-terminal domain, and RBD regions of the S1 gene, respectively. The mutations had no significant correlation with the patient's gender, age, or occupation, but they were significantly more prevalent in those with underlying diseases. 63.3% (31/49) of patients had high or very high IgG levels, and none had undetectable antibody levels at the time of reinfection. Loss of the sense of smell (69.4% - 34/49), sore throat (65.3% - 32/49), headache (59.2% - 29/49), and cough (57.1% - 28/49) were the most prevalent clinical manifestations (P > 0.05) in breakthrough-infected HCWs, aligning with the pattern of symptoms seen in the Omicron wave. However, the loss of the sense of taste showed significant results concerning clinical manifestation (P < 0.05). The examination of mutations revealed the presence of the Omicron variant in the majority of individuals. Point mutations in the C-terminal region did not significantly vary based on age, gender, or vaccine type. No significant difference was observed between vaccine types and clinical symptoms.

Conclusions:

In conclusion, this study identified spike protein mutations in reinfections among vaccinated healthcare workers. While mutations were prevalent, no significant correlations were found with demographics or vaccine types. Symptoms resembled the Omicron variant, notably with the loss of the sense of taste as a significant marker. Detectable antibody levels post-reinfection suggest that vaccine-induced immunity remains robust. Continuous monitoring of virus variants is crucial for optimizing vaccination strategies in the face of evolving strains.

1. Background

In December 2019, a new lung-related disease was reported in Wuhan, China, by Health Organizations, and it was officially named Coronavirus disease 2019 (COVID-19) (1). During the COVID-19 pandemic, healthcare workers (HCWs) were at the frontline of care for patients suspected or officially diagnosed with COVID-19 despite insufficient training or personal protective equipment (PPE). The risk among HCWs was high during the initial months of the pandemic (2-4). The pandemic is characterized by the rapid development of vaccines using conventional and new platforms, with immediate use to immunize the world population (5). Although different vaccines can significantly reduce hospitalization and mortality, breakthrough infections have been reported with all the vaccines used (6-8). The symptoms of the disease in breakthrough-infected individuals or those who have been re-infected can vary from mild to severe (9, 10).

SARS-CoV-2 is an enveloped virus with a positive-sense, single-stranded (SS), and non-segmented RNA genome that exhibits an ever-evolving nature, resulting in the emergence of multiple viral variants with alterations in genomic or amino acid profiles. Among the most important variants are those occurring in the spike protein region of the virus, which can lead to evasion of the host immune system or increased binding affinity of the receptor-binding domain (RBD) of the virus to its host angiotensin-converting enzyme II (ACE2) receptors. However, due to ongoing changes in the form of single nucleotide polymorphisms (SNPs) and viral mutations, controlling and predicting its spread is crucial. Some mutations directly affect the pathogenicity of COVID-19 and the transmission rate (11).

Important mutations in the RBD of the SARS-CoV-2 virus can lead to changes in the pathogenicity of the virus and disease severity. The formation of new variants of the SARS-CoV-2 genome, which contains the S protein necessary for ACE2 receptor binding and fusion, is crucial. The large structure S glycoprotein of SARS-CoV-2 consists of two subunits: A receptor-binding fragment (S1) and a fusion fragment (S2). S1 consists of three domains: N-terminal domain (NTD), ACE2 receptor domain, and C-terminal domain (CTD) (12-14).

In a study conducted in Iran, researchers investigated variations occurring in the spike protein of SARS-CoV-2, particularly in the receptor-binding domain (RBD). By analyzing samples from Iranian patients, they identified mutations like A475V, L452R, V483A, and F490L within this region. These mutations potentially alter the virus's functionality and its interaction with antibodies (15). Additionally, Eslami et al. investigated 176 high-quality S glycoprotein sequences of Iranian SARS-CoV-2 from public databases between April 2020 and May 2021. They identified various mutations and discussed their potential impact on infectivity and response to neutralizing antibodies. Notably, the Delta variant showed higher transmissibility in Iran and was predicted to become dominant, with some instances carrying an additional mutation (E1202Q) potentially aiding viral-cell fusion. Common mutations like N-terminal domain (NTD) deletion at position I210 and P863H in the fusion peptide-heptad repeat 1 span region were also observed in Iranian SARS-CoV-2. These mutations hold significance in predicting disease spread and influencing vaccine and drug design to combat the virus (16).

The study of mutations occurring in the S protein is of great importance for conducting scientific research on therapeutic strategies considering the immune status induced by the vaccine and the occurrence of clinical symptoms with different severity. By analyzing mutations in critical regions of the S gene of the virus, we can examine the level of genetic diversity within this group of viruses and investigate whether vaccination has been able to control genetic variations. Furthermore, by examining the extent of mutations in different regions of the S gene, we can assess whether these mutations have any impact on disease severity and the occurrence of serious complications.

2. Objectives

This study was conducted to investigate the levels of mutations in the RBD, CTD, and NTD regions of the S gene of the SARS-CoV-2 virus in breakthrough-infected healthcare workers (HCWs) during the Omicron wave in Iran. Additionally, we compared the severity of clinical symptoms with mutations occurring at key sites of the S virus gene and monitored the antibody response, which can provide insights into infection rates among HCWs.

3. Methods

3.1. Study Design

This is an observational cross-sectional study conducted from November 2021 to May 2022. The study population consisted of healthcare workers (HCWs) from Shiraz University of Medical Sciences (SUMS) who had received the COVID-19 vaccination but were reinfected. The inclusion criteria for this study were HCWs aged 18 - 60 who had received a minimum of 2 doses of the COVID-19 vaccine and tested positive for SARS-CoV-2 using the rt-qPCR assay. A total of 49 HCWs were included in the study, and the sampling was done using the census method. Written informed consent was obtained from all study participants or their legally authorized representatives.

3.2. Sample Collection

Nasal and nasopharyngeal swabs were collected from the HCWs, and RNA extraction was performed using the ROJE extraction kit (ROJE, Iran) following the manufacturer's instructions. The presence of the virus in their throat and nasal samples was tested using real-time reverse transcription-polymerase chain reaction (rt-qPCR). Positive samples were stored in a freezer at -70°C until testing. A 5 mL blood sample was taken from the HCWs to assess their IgG antibody levels. Genetic evaluation of different regions of the S gene, including the RBD, CTD, and NTD regions, was conducted using the Nested PCR method. The primers used for different regions of the S gene are listed in Table 1.

Table 1.

Primary and Nested Primers and Their Sequences Used in This Study

S Gene Segment and Oligo-NucleotideSequencesAmplicon Size BPLocation of PrimerTm
N-terminal domain
SB5 Primary F5'-AACCAGAACTCAATTACCCCC21619-2163959º
SB6 Primary R5'-TTTGAAATTACCCTGTTTTCC50522103-2212354º
SB7 Nested F5'-TCAATTACCCCCTGCATACAC21628-2164859º
SB8 Nested R5'-ATTACCCTGTTTTCCTTCAAG49022097-2211755º
ACE2 receptor
SS1 Primary F5'-TGTGTTGCTGATTATTCTGTC22643-2266355º
Binding domain
SS2 Primary R5'-AAAGTACTACTACTCTGTATG46023082-2310254º
SS3 Nested F5'-ATTCTGTCCTATATAATTCCG22656-2267654º
SS4 Nested R5'-TACTCTGTATGGTTGGTAACC43723072-2309257º
C-terminal domain
VF1 Primary F5'-AATCATTACTACAGACAACAC24901-2492154º
VF2 Primary R5'-CAATCAAGCCAGCTATAAAAC33825218-2523855º
VF3 Nested F5'-AGACAACACATTTGTGTCTGG24913-2493357º
VF4 Nested R5'-GCTATAAAACCTAGCCAAATG31525207-2522755º

3.3. Genomic Analysis

We conducted a genomic analysis of the re-infected healthcare workers (HCWs) referred to as breakthrough-infected individuals. Initially, positive samples were screened using the SARS-CoV-2 nucleic acid kit. Subsequently, the positive samples' RNA was converted into cDNA using the Viragen kit. For cDNA synthesis, 20 µL reaction tubes were prepared. After thorough vortexing and spinning, the reaction tubes underwent thermal cycling in the thermocycler: First, a 10-minute incubation at 25°C, followed by a 60-minute incubation at 47°C, and finally, a 5-minute incubation at 85°C to conclude the RNA-to-cDNA synthesis. After completion, 1 µL of the synthesized cDNA was utilized for PCR amplification. The nested PCR method was used for evaluation, and gradient PCR was performed to determine the most suitable annealing temperature based on the type of primers used in each PCR reaction. Subsequently, 10 microliters of the second-round PCR products were electrophoresed on a 2% agarose gel, and the PCR products were sent to the virology laboratory at SUMS for sequencing to validate the study results.

3.4. Enzyme-Linked Immunoassay

An enzyme-linked immunoassay (ELISA) assay was conducted to measure anti-SARS-CoV-2 IgG levels in the serum of breakthrough-infected patients. The ELISA test employed the Pishtaz Teb kit, manufactured in Iran, to detect antibodies against the SARS-CoV-2 virus in the samples.

3.5. Mutational Analysis

Amplified fragments of the CTD gene were sent to the virology laboratory at REACTED for sequence determination. Sanger sequencing was employed to determine the fragment using the ABI3500 genetic analyzer. The obtained sequences were analyzed using Blast, and their accuracy was evaluated on NCBI. Validated sequences were then registered in the gene bank. BioEdit software was utilized for sequence analysis. Initially, low-quality nucleotides at the beginning and end of chromatograms were removed, and the entire sequences were corrected. Subsequently, both sequences were assembled, and any mismatches were identified. Finally, the corrected sequences were saved as FASTA files.

3.6. Clinical Presentation

Demographic characteristics, initial symptoms, clinical signs, and 60-day mortality rates were recorded and followed by a general physician.

3.7. Statistical Analysis

The collected data were entered into the statistical software SPSS version 21 and analyzed using the analysis of variance (ANOVA) method. In this study, a significance level of a P-value less than 0.05 was considered statistically significant.

4. Results

Based on the collected information, the nursing group comprised 27 individuals, representing 55.1% of the sample size. The medical group consisted of 6 individuals, accounting for 12.2% of the sample, and the service group comprised 16 individuals, representing 32.7% of the sample.

Additional information is presented in Table 2.

Table 2.

Demographic and Immunological Characteristics of Breakthrough-Infected HCWs (N = 49)

VariablesNo. (%)
Age, y
Less than 258 (16.3)
26 - 3518 (36.7)
Above 3623 (46.9)
Sex
Male12 (24.5)
Female37 (75.5)
Type of vaccine
Sinopharm17 (34.7)
AstraZeneca19 (38.8)
Sputnik V9 (18.4)
Bharat2 (4.1)
PastoCovac2 (4.1)
Result of antibody test
Below 1208 (16.3)
121 - 17010 (20.4)
171 - 22030 (61.2)
Above 2211 (2)

The results presented in Table 2 indicate that among the age distribution groups, the highest percentage, 46.9% (23/49), was observed in individuals aged 36 years and above. Regarding gender, women constituted the majority, accounting for 75.5% (37/49) of the cases. In terms of vaccine type, AstraZeneca had the highest usage rate at 38.8% (19/49).

As shown in Table 2, 16% of the breakthrough-infected HCWs had low IgG levels (≤120 IU), while none of them had IgG levels below 50 IU. It is noteworthy that the majority of them had high IgG antibody levels (63%). The most prevalent symptoms among the breakthrough-infected subjects were loss of smell sense, which was observed in 69.4% (34/49) of cases, followed by sore throat in 65.3% (32/49), headache in 59.2% (29/49), and cough in 57.1% (28/49). Shortness of breath was reported in 30.6% (15/49) of the breakthrough-infected HCWs (Table 3). No deaths were observed within 60 days after SARS-CoV-2 re-infection.

Table 3.

Clinical Presentation of SARS-CoV-2 Infection in Breakthrough-Infected HCWs, Shiraz, Iran (N = 49) a

VariablesNegativePositiveP-Value
Fever26 (53.1)23 (46.9)≥ 0.05
Headache20 (40.8)29 (59.2)≥ 0.05
Losing smell sense15 (30.6)34 (69.4)≥ 0.05
Sore throat17 (34.7)32 (65.3)≥ 0.05
Shortness of breath34 (69.4)15 (30.6)≥ 0.05
General fatigue31 (63.3)18 (36.7)≥ 0.05
Runny nose33 (67.3)16 (32.7)≥ 0.05
Nausea36 (73.5)13 (26.5)≥ 0.05
Cough21 (42.9)28 (57.1)≥ 0.05
Losing taste sense38 (77.6)11 (22.4)≤ 0.05 b
Underlying disease≤ 0.05 b
Obesity46 (93.9)3 (6.1)
Diabetes47 (95.9)2 (4.1)
Asthma46 (93.9)3 (6.1)

The regression analysis revealed a significant positive relationship (P = 0.001, t = 6.196) between the impact of point mutations in the C-terminal region and the antibody levels in the healthcare staff. Table 4 displays the differences in point mutations in the C-terminal region based on demographic and background variables.

Table 4.

Some Spike Protein Mutations in the N-terminal, RBD, and C-terminal Regions

Sample CodeN-ter (SB)RBD (SS)C-ter (VF)Mutations
5××E65Q
6××No mutation
8×No mutation
9×E65Q
19×M36I, L37I, K45E, R51X, S61X
31×No mutation
32×No mutation
33××S61X, Y33X, S35X
34×No mutation
35×No mutation
38×S61X,
39×No mutation
41×D26N, R29S, K38N, N61K, S98N, E105A

Mutations were observed in the 25, 10, and 1 sequence of the C-terminal domain (CTD), N-terminal domain (NTD), and RBD (receptor binding domain) regions of the S gene, respectively. Amino acid mutations in the CTD included E65Q, M36I, K45E, R51X, and S61X, while those in the N-terminal region included S35X, Y33X, and S61X. Additionally, amino acid mutations in the RBD region included D26N, R29S, K38N, N61K, S98N, and E105A. Most of the mutations that occurred in the S region of the genome were related to the CTD and can be analyzed and evaluated statistically. Protein sequence analysis of the S protein did not reveal any significant mutations in the virus. Detailed information regarding the mutations is provided in Table 4. The nucleotide sequence data were deposited in the GenBank database under the accession numbers “OR073457.1, OR073473.1, OR073458.1, OR073474.1, OR073475.1, OR073459.1, OR073476.1, OR073460.1, OR073477.1, OR073461.1, OR073478.1, OR073462.1, OR073463.1, OR073479.1, OR073464.1, OR073465.1, OR073468.1, OR073480.1, OR073470.1, OR073481.1, OR073471, OR073472.1, OR073482.1”.

The independent t-test showed that the point mutation in the C-terminal region did not significantly differ based on gender. Furthermore, the ANOVA indicated that the point mutation in the C-terminal region did not significantly differ based on age and vaccine type. According to the conducted ANOVA, there was no significant difference between vaccine types in clinical symptoms such as fever, chills, headache, loss of smell, sore throat, diarrhea, fatigue, shortness of breath, body aches, runny nose, nausea, vomiting, cough, and loss of taste.

4.1. Phylogenetic Tree

The phylogenetic results of our study on the S gene of SARS-CoV-2 demonstrated similarity with Epi delta-oq414654 corona Iran-om858820 Omicron-NC 045512.2 Wuhan (Figure 1).

Phylogenetic tree based on sequencing and analysis of the S region
Phylogenetic tree based on sequencing and analysis of the S region

5. Discussion

In this study, specific mutations in the breakthrough-infected HCWs have been examined, but no significant mutations were observed in most patients who had received two doses of the vaccine before being infected with COVID-19. These results indicate that the vaccines these individuals received were effective in reducing the development of key mutations in the S-critical region gene, particularly in RBD. On the other hand, vaccination can be effective in reducing the transmission of the Omicron variant and improving the condition of vaccinated individuals who have acquired this variant. Several studies have shown reductions in hospital admissions, mortality, and ICU admissions, in addition to reducing disease severity in people with breakthrough infections (17, 18).

In the current study, certain factors played a role in determining the severity of COVID-19 disease. The most important of these were the use of a ventilator, an oxygen reservoir, the use of vasopressors, and treatment with IVIG or corticosteroids. None of our patients had severe disease, and only 31% of them had mild disease. This study revealed that the severity of the breakthrough infection among HCWs during the COVID-19 Omicron wave was low, which is consistent with other studies (19). Based on the results of this study, it can be said that the vaccines developed against the Omicron variant can still maintain their effectiveness in preventing the disease and reducing its severity in cases of infection.

Several studies have investigated the mutation rate in different regions of the S gene of the SARS-CoV-2 virus among vaccinated individuals. However, there is a lack of specific research examining the mutation rate in the RBD, CTD, and NTD regions of the S gene in vaccinated individuals. One study demonstrated that the mutation rate of SARS-CoV-2 is estimated to be around 1.12 × 10-3 mutations per site per year, which is similar to those observed in the influenza virus and HIV (20). Another study focused on the mutation rate of SARS-CoV-2 during experimental evolution and revealed that the basic mutation rates ranged from 10-6 to 10-3 for various positive-sense RNA viruses (21). Based on the results of this study, it appears that there is a significant and positive relationship between point mutations in the C-terminal region and antibodies in HCWs.

This finding indicates that point mutations in this region may have an impact on the functionality of antibodies against the Omicron variant. It can be suggested that point mutations in the C-terminal region could introduce changes in the protein structure of antibodies, which may lead to a decrease in the capacity of the antibodies to detect and neutralize the Omicron variant (22, 23). In other words, these mutations can potentially reduce the effectiveness of antibodies against the Omicron variant, making individuals with these mutations more susceptible to Omicron infection (24). However, further studies are required to confirm this interpretation and fully elucidate the effect of point mutations in the C-terminal region on antibodies and their responses to the Omicron variant. Additionally, it should be noted that other variables may also influence the effectiveness of antibodies against the Omicron variant, such as antibody levels, individual genetic diversity, and other environmental and educational factors that may play a role in antibody performance (24).

The results of the current study showed that there was no significant difference in the impact of point mutations in the C-terminal region between men and women. This finding indicates that the effect of point mutations in this region on the performance of antibodies against the Omicron variant is comparable between genders. Additionally, the results demonstrated that there was no difference in the impact of point mutations in the C-terminal region based on age and vaccine type. This result suggests that point mutations in this region affect antibody functionality against the Omicron variant similarly across individuals of different ages and vaccine types.

It should be noted that the differences between groups may be influenced by a combination of various factors (25). For example, individuals in different age groups may have received different types of vaccines, which could have an impact on the results (26). Additionally, the presence of underlying diseases in some groups may also contribute to the differences. However, the results showed that there was a difference in the impact of point mutations in the C-terminal region based on the presence of underlying diseases. This finding may indicate that the presence of underlying diseases can influence the severity and functionality of antibodies against the Omicron variant. Furthermore, it is worth mentioning that other factors, such as concurrent infection with other strains of the COVID-19 virus, hygiene conditions, and environmental factors, may also play a role in the impact of point mutations in the C-terminal region (27).

The results of the current study showed that there was no difference in clinical symptoms such as fever, chills, headache, loss of taste and smell, sore throat, diarrhea, fatigue, shortness of breath, body aches, nausea, vomiting, cough, and runny nose among the five types of vaccines: Sinovac, AstraZeneca, Sputnik V, Bharat, and Pfizer. This finding indicates that all five types of vaccines have performed similarly in reducing clinical symptoms of COVID-19, and there is no significant difference among them in this regard. This result is important because it demonstrates that each of these vaccine types is capable of controlling and reducing clinical symptoms of COVID-19.

Individuals who are concerned about these clinical symptoms during vaccination can use any of these vaccine types without having particular concerns about differences in clinical symptoms (28, 29). Gaborieau et al.'s study reported on the prevalence and clinical symptoms of COVID-19. In this study, people with different underlying diseases such as asthma, obesity, and immunodeficiency were also examined, and overall, the most common clinical symptoms in patients were fever, runny nose, digestive problems, and respiratory distress (30). In another study in California, cough, runny nose, and general fatigue were reported as symptoms of patients with COVID-19 (31). Compared to our study. We observed loss of smell sense, sore throat, cough, and fever among the most presented clinical manifestations. Our study includes some limitations, such as the presence of underlying diseases, which may subsequently affect the results. Further investigations on a higher population scale and in other groups within society are needed.

5.1. Conclusions

The study discovered that vaccines developed for the Omicron variant are effective in preventing infections and reducing disease severity. Nevertheless, mutations in the C-terminal region may potentially weaken the antibody response to Omicron. Further research is required to gain a comprehensive understanding of this phenomenon. Gender and vaccine type did not appear to significantly influence these mutations, but underlying health conditions could potentially impact the effectiveness of antibodies. Additionally, all five vaccine types demonstrated similar efficacy in reducing COVID-19 symptoms, implying that any of them can be used effectively against Omicron.

References

  • 1.

    Macera M, De Angelis G, Sagnelli C, Coppola N, Vanvitelli C. Clinical presentation of COVID-19: Case series and review of the literature. Int J Environ Res Public Health. 2020;17(14). [PubMed ID: 32674450]. [PubMed Central ID: PMC7399865]. https://doi.org/10.3390/ijerph17145062.

  • 2.

    Trbojevic-Akmacic I, Petrovic T, Lauc G. SARS-CoV-2 S glycoprotein binding to multiple host receptors enables cell entry and infection. Glycoconj J. 2021;38(5):611-23. [PubMed ID: 34542788]. [PubMed Central ID: PMC8450557]. https://doi.org/10.1007/s10719-021-10021-z.

  • 3.

    Hasoksuz M, Kilic S, Sarac F. Coronaviruses and SARS-COV-2. Turk J Med Sci. 2020;50(SI-1):549-56. [PubMed ID: 32293832]. [PubMed Central ID: PMC7195990]. https://doi.org/10.3906/sag-2004-127.

  • 4.

    Shoja E, Aghamohammadi V, Bazyar H, Moghaddam HR, Nasiri K, Dashti M, et al. Covid-19 effects on the workload of Iranian healthcare workers. BMC Public Health. 2020;20(1):1636. [PubMed ID: 33138798]. [PubMed Central ID: PMC7605333]. https://doi.org/10.1186/s12889-020-09743-w.

  • 5.

    Doke PP, Mhaske ST, Oka G, Kulkarni R, Muley V, Mishra AC, et al. SARS-CoV-2 breakthrough infections during the second wave of COVID-19 at Pune, India. Front Public Health. 2022;10:1040012. [PubMed ID: 36711329]. [PubMed Central ID: PMC9877521]. https://doi.org/10.3389/fpubh.2022.1040012.

  • 6.

    Jalil M, Ashkan Z, Gholamnezhad M, Jamalidoust S, Jamalidoust M. Effect of COVID-19 on healthcare workers' morbidity and mortality compared to the general population in Kohgiluyeh and Boyer-Ahmad Province, Iran. Health Sci Rep. 2023;6(1). e961. [PubMed ID: 36519077]. [PubMed Central ID: PMC9742827]. https://doi.org/10.1002/hsr2.961.

  • 7.

    Patalon T, Gazit S, Pitzer VE, Prunas O, Warren JL, Weinberger DM. Odds of testing positive for SARS-CoV-2 following receipt of 3 vs 2 doses of the BNT162b2 mRNA vaccine. JAMA Intern Med. 2022;182(2):179-84. [PubMed ID: 34846533]. [PubMed Central ID: PMC8634151]. https://doi.org/10.1001/jamainternmed.2021.7382.

  • 8.

    Jamalidoust M, Eilami O, Ashkan Z, Ziyaeyan M, Aliabadi N, Habibi M. The rates and symptoms of natural and breakthrough infection pre- and post- Covid-19 non-mRNA vaccination at various peaks amongst Iranian healthcare workers. Virol J. 2023;20(1):182. [PubMed ID: 37596593]. [PubMed Central ID: PMC10436397]. https://doi.org/10.1186/s12985-023-02156-2.

  • 9.

    Rodda LB, Netland J, Shehata L, Pruner KB, Morawski PA, Thouvenel CD, et al. Functional SARS-CoV-2-specific immune memory persists after mild COVID-19. Cell. 2021;184(1):169-183 e17. [PubMed ID: 33296701]. [PubMed Central ID: PMC7682481]. https://doi.org/10.1016/j.cell.2020.11.029.

  • 10.

    Arafa A, Mohammed Z, Mahmoud O, Elshazley M, Ewis A. Depressed, anxious, and stressed: What have healthcare workers on the frontlines in Egypt and Saudi Arabia experienced during the COVID-19 pandemic? J Affect Disord. 2021;278:365-71. [PubMed ID: 33007626]. [PubMed Central ID: PMC7831891]. https://doi.org/10.1016/j.jad.2020.09.080.

  • 11.

    Yin C. Genotyping coronavirus SARS-CoV-2: methods and implications. Genomics. 2020;112(5):3588-96. [PubMed ID: 32353474]. [PubMed Central ID: PMC7184998]. https://doi.org/10.1016/j.ygeno.2020.04.016.

  • 12.

    Zhang S, Qiao S, Yu J, Zeng J, Shan S, Tian L, et al. Bat and pangolin coronavirus spike glycoprotein structures provide insights into SARS-CoV-2 evolution. Nat Commun. 2021;12(1):1607. [PubMed ID: 33707453]. [PubMed Central ID: PMC7952905]. https://doi.org/10.1038/s41467-021-21767-3.

  • 13.

    Durmaz B, Abdulmajed O, Durmaz R. Mutations observed in the SARS-CoV-2 spike glycoprotein and their effects in the interaction of virus with ACE-2 receptor. Medeni Med J. 2020;35(3):253-60. [PubMed ID: 33110678]. [PubMed Central ID: PMC7584268]. https://doi.org/10.5222/MMJ.2020.98048.

  • 14.

    Aliabadi N, Jamaliduost M, Pouladfar G, Marandi NH, Ziyaeyan M. Characterization and phylogenetic analysis of Iranian SARS-CoV-2 genomes: A phylogenomic study. Health Sci Rep. 2023;6(1). e1052. [PubMed ID: 36686884]. [PubMed Central ID: PMC9841325]. https://doi.org/10.1002/hsr2.1052.

  • 15.

    Karami C, Asghariazar V, Pahlavan Y, Mazloum Jalali K, Tajehmiri A, GhorbaniVanan M, et al. Mutation interest in the receptor-binding domain of spike protein of severe acute respiratory syndrome coronavirus 2 in the Northwest of Iran. Jundishapur Journal of Microbiology. 2023;16(8). https://doi.org/10.5812/jjm-137863.

  • 16.

    Eslami S, Glassy MC, Ghafouri-Fard S. A comprehensive overview of identified mutations in SARS CoV-2 spike glycoprotein among Iranian patients. Gene. 2022;813:146113. [PubMed ID: 34896524]. [PubMed Central ID: PMC8656178]. https://doi.org/10.1016/j.gene.2021.146113.

  • 17.

    Arbel R, Hammerman A, Sergienko R, Friger M, Peretz A, Netzer D, et al. BNT162b2 vaccine booster and mortality due to COVID-19. N Engl J Med. 2021;385(26):2413-20. [PubMed ID: 34879190]. [PubMed Central ID: PMC8728797]. https://doi.org/10.1056/NEJMoa2115624.

  • 18.

    CDC Covid- Vaccine Breakthrough Case Investigations Team. COVID-19 vaccine breakthrough infections reported to CDC - United States, January 1-April 30, 2021. MMWR Morb Mortal Wkly Rep. 2021;70(21):792-3. [PubMed ID: 34043615]. [PubMed Central ID: PMC8158893]. https://doi.org/10.15585/mmwr.mm7021e3.

  • 19.

    Sena GR, Lima TPF, Vidal SA, Duarte M, Bezerra PGM, Fonseca Lima EJ, et al. Clinical characteristics and mortality profile of COVID-19 patients aged less than 20 years old in Pernambuco - Brazil. Am J Trop Med Hyg. 2021;104(4):1507-12. [PubMed ID: 33606669]. [PubMed Central ID: PMC8045659]. https://doi.org/10.4269/ajtmh.20-1368.

  • 20.

    McLean G, Kamil J, Lee B, Moore P, Schulz TF, Muik A, et al. The impact of evolving SARS-CoV-2 mutations and variants on COVID-19 vaccines. mBio. 2022;13(2). e0297921. [PubMed ID: 35352979]. [PubMed Central ID: PMC9040821]. https://doi.org/10.1128/mbio.02979-21.

  • 21.

    Amicone M, Borges V, Alves MJ, Isidro J, Ze-Ze L, Duarte S, et al. Mutation rate of SARS-CoV-2 and emergence of mutators during experimental evolution. Evol Med Public Health. 2022;10(1):142-55. [PubMed ID: 35419205]. [PubMed Central ID: PMC8996265]. https://doi.org/10.1093/emph/eoac010.

  • 22.

    Beeraka NM, Sukocheva OA, Lukina E, Liu J, Fan R. Development of antibody resistance in emerging mutant strains of SARS CoV-2: Impediment for COVID-19 vaccines. Rev Med Virol. 2022;32(5). e2346. [PubMed ID: 35416390]. [PubMed Central ID: PMC9111059]. https://doi.org/10.1002/rmv.2346.

  • 23.

    K A, Sharma A, Kumar D, Singh SK, Gupta G, Chellappan DK, et al. Molecular aspects of Omicron, vaccine development, and recombinant strain XE: A review. J Med Virol. 2022;94(10):4628-43. [PubMed ID: 35705439]. [PubMed Central ID: PMC9349635]. https://doi.org/10.1002/jmv.27936.

  • 24.

    van Gils MJ, Lavell A, van der Straten K, Appelman B, Bontjer I, Poniman M, et al. Antibody responses against SARS-CoV-2 variants induced by four different SARS-CoV-2 vaccines in health care workers in the Netherlands: A prospective cohort study. PLoS Med. 2022;19(5). e1003991. [PubMed ID: 35580156]. [PubMed Central ID: PMC9113667]. https://doi.org/10.1371/journal.pmed.1003991.

  • 25.

    Khandia R, Singhal S, Alqahtani T, Kamal MA, El-Shall NA, Nainu F, et al. Emergence of SARS-CoV-2 Omicron (B.1.1.529) variant, salient features, high global health concerns and strategies to counter it amid ongoing COVID-19 pandemic. Environ Res. 2022;209:112816. [PubMed ID: 35093310]. [PubMed Central ID: PMC8798788]. https://doi.org/10.1016/j.envres.2022.112816.

  • 26.

    Cerqueira-Silva T, Oliveira VA, Boaventura VS, Pescarini JM, Junior JB, Machado TM, et al. Influence of age on the effectiveness and duration of protection of Vaxzevria and CoronaVac vaccines: A population-based study. Lancet Reg Health Am. 2022;6:100154. [PubMed ID: 34957437]. [PubMed Central ID: PMC8692070]. https://doi.org/10.1016/j.lana.2021.100154.

  • 27.

    Nguyen KV. Containing the spread of COVID-19 virus facing to its high mutation rate: approach to intervention using a nonspecific way of blocking its entry into the cells. Nucleosides Nucleotides Nucleic Acids. 2022;41(8):778-814. [PubMed ID: 35532338]. https://doi.org/10.1080/15257770.2022.2071937.

  • 28.

    Hodgson SH, Mansatta K, Mallett G, Harris V, Emary KRW, Pollard AJ. What defines an efficacious COVID-19 vaccine? A review of the challenges assessing the clinical efficacy of vaccines against SARS-CoV-2. Lancet Infect Dis. 2021;21(2):e26-35. [PubMed ID: 33125914]. [PubMed Central ID: PMC7837315]. https://doi.org/10.1016/S1473-3099(20)30773-8.

  • 29.

    Chung JY, Thone MN, Kwon YJ. COVID-19 vaccines: The status and perspectives in delivery points of view. Adv Drug Deliv Rev. 2021;170:1-25. [PubMed ID: 33359141]. [PubMed Central ID: PMC7759095]. https://doi.org/10.1016/j.addr.2020.12.011.

  • 30.

    Gaborieau L, Delestrain C, Bensaid P, Vizeneux A, Blanc P, Garraffo A, et al. Epidemiology and clinical presentation of children hospitalized with SARS-CoV-2 infection in suburbs of Paris. J Clin Med. 2020;9(7). [PubMed ID: 32674306]. [PubMed Central ID: PMC7408757]. https://doi.org/10.3390/jcm9072227.

  • 31.

    Lewnard JA, Mora AM, Nkwocha O, Kogut K, Rauch SA, Morga N, et al. Prevalence and clinical profile of SARS-CoV-2 infection among farmworkers in Monterey County, California: June-November, 2020. medRxiv. 2021. https://doi.org/10.1101/2020.12.27.20248894.