An in silico analysis of early SARS-CoV-2 variant B.1.1.529 (Omicron) genomic sequences and their epidemiological correlates ============================================================================================================================ * Ashutosh Kumar * Adil Asghar * Himanshu N. Singh * Muneeb A. Faiq * Sujeet Kumar * Ravi K. Narayan * Gopichand Kumar * Prakhar Dwivedi * Chetan Sahni * Rakesh K. Jha * Maheswari Kulandhasamy * Pranav Prasoon * Kishore Sesham * Kamla Kant * Sada N. Pandey ## Abstract **Background** A newly emerged SARS-CoV-2 variant B.1.1.529 has worried the health policy makers worldwide due to the presence of a large number of mutations in its genomic sequence, especially in the spike protein region. World Health Organization (WHO) has designated it as a global variant of concern (VOC) and has named as ‘Omicron’. A surge in new COVID-19 cases have been reported from certain geographical locations, primarily in South Africa (SA) following the emergence of Omicron. **Materials and methods** We performed an ***in silico*** analysis of the complete genomic sequences of Omicron available on GISAID (until 2021-12-10) to predict the functional impact of the mutations present in this variant on virus-host interactions in terms of viral transmissibility, virulence/lethality, and immune escape. The mutations present at the receptor binding domain (RBD) of the variants were assessed using an open analysis pipeline which integrates a yeast-display platform with deep mutational scanning. Further, we performed a correlation analysis of the relative proportion of the genomic sequences of specific SARS-CoV-2 variants (in the period of 01 Oct-10th Dec, 2021) with the current epidemiological data (new COVID-19 cases and deaths) from SA to understand whether the Omicron has an epidemiological advantage over existing variants. **Results** Compared to the current list of global VOCs/VOIs (as per WHO) Omicron bears more sequence variation, specifically in the spike protein and host receptor-binding motif (RBM). Omicron showed the closest nucleotide and protein sequence homology with Alpha variant for the complete sequence as well as for RBM. The mutations were found primarily condensed in spike region (28-48) of the virus. Further, the mutational analysis showed enrichment for the mutations decreasing ACE2-binding affinity and RBD protein expression, in contrast, increasing the propensity of immune escape. An inverse correlation of Omicron with Delta variant was noted (r=-0.99, p< .001, 95% CI: -0.99 to - 0.97) in the sequences reported from SA post-emergence of the new variant, later showing a decrease. There has been a steep rise in the new COVID-19 cases in parallel with increase in the proportion of Omicron since the first case (74-100%), on the contrary, the incidences of new deaths have not been increased (r=-0.04, p>0.05, 95% CI =-0.52 to 0.58). **Conclusions** Omicron may have greater immune escape ability than the existing VOCs/VOIs. However, there are no clear indications coming out from the predictive mutational analysis that the Omicron may have higher virulence/lethality than other variants, including Delta. The higher ability for immune escape may be a likely reason for the recent surge in Omicron cases in SA. **Highlights** * Higher immune escape ability than the existing VOCs/VOIs * No clear indications of increased affinity for ACE2 binding * Driving a new COVID-19 wave in South Africa * Outcompeting Delta variant * Currently, no clear evidence for increased virulence/lethality Key words * SARS-CoV-2 variant * COVID-19 * pandemic * immune escape * transmissibility * virulence ## Introduction A new variant of SARS-CoV-2 (lineage B.1.1.529) has been reported from Botswana, South Africa (SA) and multiple other countries1. World Health Organization (WHO) has designated it as a global variant of concern (VOC) and has named as ‘Omicron’2. The new variant has been given a pango-lineage as BA.1. Currently very little is known about epidemiological characteristics of this variant. The presence of a large number of mutations in its genomic sequence, especially in the spike protein region including in the host receptor-binding domain (RBD) has raised the speculations that Omicron can prove a serious epidemiological threat and can be a reason for the next COVID-19 waves globally3. More recently a sub-lineage (pango-lineage BA.2) has appeared in Omicron with slightly varying set of mutations. We performed an ***in silico*** analysis of the complete genomic sequences, of Omicron available on GISAID, to predict the functional impact of the mutations, present in this variant, on virus-host interactions in terms of viral transmissibility, virulence, and the immune escape capabilities. Further, we made an assesment of the relative proportion of the genomic sequences of the SARS-CoV-2 variants and correlated that with the rise in new COVID-19 cases in a most affected global geographical location by Omicron, as to understand whether the new variant has an epidemiological advantage in terms of transmissibility and virulence/lethality over the existing variants. The findings of this study provides important insights about origin and the epidemiological attributes of this variant paving way for the future research. ## Methods ### Data collection #### SARS-CoV-2 genomic sequences The SARS-CoV-2 genomic sequence for Omicron variant and other global variants of concern/interest (VOCs/VOIs) were downloaded from EpiCoVTM database of Global Initiative on Sharing All Influenza Data (GISAID) ([https://www.gisaid.org/](https://www.gisaid.org/)) using automatic search function feeding information for geographical location, SARS-CoV-2 lineage and sample collection and sequence reporting dates (up to 2021-12-10). The optimum length and coverage of the downloaded sequences (used for variant comparisons) were achieved with selecting ‘complete sequence’ and ‘high coverage’ in the search function. #### Epidemiological data The global epidemiological data (daily new cases and deaths) for COVID-19 for the period of 1st Oct-10th Dec, 2021 were accessed from Worldometer ([https://www.worldometers.info/coronavirus](https://www.worldometers.info/coronavirus)). #### Data analysis The mutational analysis on the genomic sequences was performed and 3D structure of spike with amino acid changes in Omicron was generated using CoVsurver app by GISAID ([https://www.gisaid.org/](https://www.gisaid.org/)) employing hCoV-19/Wuhan/WIV04/2019 as a reference strain. Further, a comparative mutational analysis of Omicron with existing global VOCs/VOIs (as per WHO)4 was generated using ‘Compare lineages’ function at Outbreak.info ([https://outbreak.info/](https://outbreak.info/)) using GISAID as the source of genomic sequence data. The protein sequence translation from the viral genomic sequences were done using Expasy by Swiss Institute of Bioinformatics ([https://www.expasy.org/](https://www.expasy.org/)). A comparative assessment of the Omicron nucleotide and protein sequences with existing global VOCs/VOIs were made using NCBI Blast tool ([https://blast.ncbi.nlm.nih.gov/Blast.cgi?](https://blast.ncbi.nlm.nih.gov/Blast.cgi?)). Furthermore, the functional impact of the mutations present at the receptor binding domain (RBD) of the variants were assessed using an open analysis pipeline developed by Starr *et al*, 2020 which integrates a yeast-display platform with deep mutational scanning to determine how all amino acid mutations to measure how all possible RBD amino acid mutations affect ACE2-binding affinity and protein expression (a correlate of protein folding stability) as compared to the wild-type (WT) SARS-CoV-2 strain ([https://jbloomlab.github.io/SARS-CoV-2-RBD_DMS/](https://jbloomlab.github.io/SARS-CoV-2-RBD_DMS/)). The epidemiological correlates of the Omicron variant were assessed based on the comparative analysis of the genomic sequences from GISAID and current epidemiological data (daily new cases and deaths) from one of the most affected region by this variant (last date of collection: 2021-12-10). The number of sequences for each SARS-CoV-2 variants was retrieved using automatic search function feeding information for the lineage and collection dates in EpiCoVTM database of GISAID for the period of 1st Oct, 2021 to 10th Dec, 2021. A 3 day sum of the total number of sequences was noted for each variant and their relative proportions were calculated (in percentage). Data was tabulated, and the distribution of each variant was charted against the COVID-19 epidemiological data (3 day sum of new case and deaths) and statistically analyzed for appreciating the changes in relationship between the variables before and after the emergence of Omicron. #### Statistics An expect (E) value ≤ 0 was considered significant for the sequence homology match through NCBI Blast. An *E*-value close to zero or below together with the higher Max score signified a higher sequence homology ranking. (Further details of the statistical methods predicting significance can be consulted at [https://www.ncbi.nlm.nih.gov/BLAST/tutorial/Altschul-1.html](https://www.ncbi.nlm.nih.gov/BLAST/tutorial/Altschul-1.html)). For the mutational analysis only the mutations present in atleast 75% of sequenced samples were considered for functional characterisation. For the analysis of epidemiological data the statistical tests were performed to evaluate inter-group differences between the SARS-CoV-2 Variants with the help of MS Excel 2019 and XLSTAT package. Normality of the data was examined using Shapiro-Wilk test. Pearson’s correlation and Spearman’s rho tests were performed for the normally distributed and skewed data respectively. A correlation matrix was generated and a linear regression analysis was performed between the comparing variables [presented as r values= 0 to 1 at confidence interval (CI) of 95%]. Results were considered statistically significant at p-value ≤ 0.05. Graphs were plotted to visualize the data trends. #### Ethical statement An ethical clearance was precluded as *in silico* methods and open access/globally shared databases were used to conduct this study. ## Results A total of 3604 genomic sequences of Omicron from 54 countries were uploaded on the GISAID until 2021-12-10 (Fig. S1) which were analyzed for the mutational characteristics. The highest number of Omicron genomic sequences (81.5% of the total) were reported from SA where a rise in the COVID-19 cases has been observed post-emergence of the new variant ([https://www.worldometers.info/coronavirus/country/south-africa/dated2021-12-10](https://www.worldometers.info/coronavirus/country/south-africa/dated2021-12-10)). The mutations found were primarily condensed in spike region (28-48) of the virus, however, frequent non-spike mutations were also noted (20-26). In addition, eight sequences of BA.2 was reported from four countries which included newer mutations in all viral parts including spike protein (Figs. S2-3). In this study, we have focused on analyzing the genomic sequences of currently most prevalent sublineage of Omicron (BA.1). ### Sequence homology of Omicron (BA.1) with WT strains Compared to the current list of global VOCs/VOIs (as per WHO) Omicron showed more sequence variation, more specifically in the spike protein (n21563-25384: 1-1273aa) including RBM (n22869-23089: 438-508aa), where the variations were most prominent (Table S1). The homology of the Omicron to the reference strain (hCoV-19/Wuhan/WIV04/2019) for the spike protein sequence varied from 96.23 to 97.8% (28-48 mutations) in the analyzed sequences (Fig. 1). ![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/22/2021.12.18.21267908/F1.medium.gif) [Figure 1](http://medrxiv.org/content/early/2021/12/22/2021.12.18.21267908/F1) Figure 1 3-D structure of Omicron (BA.1) spike glycoprotein in interaction of human ACE2 showing key amino acid substitutions. (Data source: The figure has been generated from the genomic sequences of Omicron using CoVsurver app from [Gisaid.org](http://Gisaid.org).) ### Sequence homology of Omicron (BA.1) with existing SARS-CoV-2 VOCs/VOIs The analysis of genomic and protein sequence homology of Omicron with the reference strain and current global VOCs/VOIs (as per WHO) showed the highest similarity of Omicron with the Alpha variant for the complete sequence as well as RBM, however, the highest similarity for the complete nucleotide and protein sequences for the spike protein were noted with Beta and Delta variants, respectively (Table S1). ### Mutational analysis The multiple clusters of closely spaced mutations were noted across the sequence, which were most densely placed in the spike protein region, particularly in its S1 subunit, including the host receptor-binding motif (RBM) (Fig. 2, Table 2). Many of the mutations in Omicron are common with the current global VOCs/VOIs (Fig. 3). Table 1 provides a list of the mutations in Omicron (BA.1) (present in at least 75% of sequences) and their functional characteristics based on the existing literature evidence. Based on the available evidence these mutations in pango-lineage BA.1 can be broadly categorized into four major groups: 1. Host receptor binding (10), 2. Host-adaptibility (3), 3. Immune escape (20), and 4. Virus replication (18) (Table 1). ![Figure 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/22/2021.12.18.21267908/F2.medium.gif) [Figure 2](http://medrxiv.org/content/early/2021/12/22/2021.12.18.21267908/F2) Figure 2 The mutational landscape in SARS-CoV-2 variant B.1.1.529 (Omicron, sublineage: BA.1). Abbreviations: RBD— receptor binding domain, RBM—receptor binding motif, ORF—open reading frame. ![Figure 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/22/2021.12.18.21267908/F3.medium.gif) [Figure 3](http://medrxiv.org/content/early/2021/12/22/2021.12.18.21267908/F3) Figure 3 Lineage comparison between Omicron and other global variants of concerns/interest. Showing the mutations with > 75% prevalence in at least one lineage. (Data source: The figure has been generated from outbreak.info, based on the SARS-CoV-2 genomic sequences uploaded in GISAID until 2021-12-6.) View this table: [Table 1](http://medrxiv.org/content/early/2021/12/22/2021.12.18.21267908/T1) Table 1 List of mutations in SARS-CoV-2 variant B.1.1.529 (Omicron, sublineage: BA.1). View this table: [Table 2](http://medrxiv.org/content/early/2021/12/22/2021.12.18.21267908/T2) Table 2 The predicted impact of receptor binding motif (RBM) variations in SARS-CoV-2 variant B.1.1.529 The analysis of the mutations present at the RBD using deep mutational scanning pipeline by Starr *et al*, 2020 reflected prominent ACE2-binding affinity and protein expression changes (Table 2). Notably, mutations decresing the ACE2-binding affinity and protein expression were found significantly greater in number. ### Epidemiological correlates A total (N) of 4224 SARS-CoV-2 genomic sequences (NDelta = 999, NOmicron = 2937, and NOthers = 288) were uploaded on GISAID from SA in the period of study. For the complete duration of the study Delta correlated negatively with the new COVID-19 cases (r=-0.567, p=0.004, 95% CI=-0.79 to - 0.21), however positively with the new deaths (r=0.38, p=0.065, 95% CI= -0.025 to 0.68). The differential analysis of the SARS-CoV-2 genomic sequences from SA prior and after emrgence of the first case of Omicron (dated 2021-11-05, EPI_ISL_7456440) reflected sharp change in the dominance of the variant from Delta to Omicron (Fig. 4). An inverse correlation of Omicron with Delta variant was noted (r=-0.99, p< .001, 95% CI: -0.99 to -0.97) in the period of study. There has been a steep rise in the new COVID-19 cases in parrallel with the increase in proportion of Omicron since the first case of Omicron (74-100% of total genomic sequences after 15-17th Nov, 2021), however, no parrallel increase has been observed in the death cases, which otherwise showed a reverse trend (r=-0.04, p>0.05, 95% CI =-0.52 to 0.58) (Fig. 4). ![Figure 4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/22/2021.12.18.21267908/F4.medium.gif) [Figure 4](http://medrxiv.org/content/early/2021/12/22/2021.12.18.21267908/F4) Figure 4 Epidemiological correlates of Omicron and Delta variants genomic sequences reported on GISAID from South Africa for the period of 1st Oct-10th Dec, 2021. The proportion of Delta and Omicron variants in the total SARS-CoV-2 genomic sequences were correlated with the new COVID-19 cases and deaths in the period of study (3-day sum of each variable). A sharp change in the dominance from Delta to Omicron is observable since the first case (5th Nov, 2021). The rise of Omicron cases has parrelled with the rise in the new COVID-19 cases. In comparison, the Delta variant is showing a fall in the same time period. Notably, there has been no increase in the number of deaths post emergence of Omicron. (Data source: [gisaid.org](http://gisaid.org) and worldometer.info.) ## Discussion Our analysis of the SARS-CoV-2 genomic sequences coupled with epidemiological data from SA unravels multiple observations regarding host-virus interactions which may help predicting the epidemiological potentials of the Omicron. ### Sequence variation with WT strains and existing SARS-CoV-2 VOCs/VOIs Our analyses showed that Omicron has significantly greater variantions than current VOCs and VOIs in comparisons with the wild type strains (Wuhan strain and B.1). The sequence variation was greatest in the spike protein region, more particularly in the receptor-binding motif (RBM). Among the existing VOCs and VOIs, the highest homology of the complete sequence and RBM (nucleotide as well as protein) of Omicron was observed with Alpha variant (Table S1). Interestingly, similar to Alpha variant S-gene target failure (SGTF) in PCR based detection is proving a sensitive method for the detection of Omicron in the clinical samples23. As Omicron bears key mutations from multiple existing VOCs/VOIs, with approximate sequence homology variation, rather than a direct descendency, the multiple recombination events between the variants inside the hosts can be a more plausible explanation for its origin. Notably, compared with the WT strains, Omicron has significantly higher sequence variation than existing VOCs/VOIs, more specifically in spike protein region (Table S1). Further, in the spike protein region, the highest sequence variations is noted in RBM (Table S1). It will be pertinent to explore the evolutionary mechanisms involved in accumulation of such a large number of mutations in Omicron. The speculations have also been raised that the long term persistence of SARS-CoV-2 infection in an immunocompromised subject can be a probable mechanism behind the Omicron origin. Existing literature24–27 supports strong plausibility of this mechanism. Avanzato *et al*24 and Choi *et al*26 reported case studies of persistence of infection and accumulation of novel mutations in SARS-CoV-2 spike gene and RBD in a chronically ill and immunocompromised COVID-19 patients. Another such case has been recently reported by Karim *et al*25. Authors documented long persistence of SARS-CoV-2 infection (more than six month) in a patient with advanced HIV and antiretroviral treatment failure. Through whole genome sequencing for SARS-CoV-2, performed at multiple time-points from the patient samples, authors demonstrated early emergence of the E484K substitution followed by N501Y, K417T, and multiple others mutation (including some novel mutations) in the spike gene and RBD. Increase in the genomic diversity reflecting intra-host evolution of SARS-CoV-2 during prolonged infection was also noted in a recent cohort study by Voloch *et al*27. ### Predicted effect on virus-host interactions The new variant has accumulated multiple closely spaced mutations at receptor-binding motif (RBM) with ACE2, specifically a cluster of mutations have appeared at 477-501aa (Fig. 2). Notably, it has many of the mutations common with the earlier variants of concern (Fig. 3), many of which have been shown to enhance RBD:ACE2 binding in comparison to the WT strain28, as has also been reflected in our analysis (Table 1-2). The selective mutations which are either present at or in the vicinity of the RBM (N440K, S477N, T478K, and N501Y) in majority of the Omicron sequences are believed to stabilize binding with ACE2 (Table 2). D614G, a key mutation characteristically present in all B.1 descendents28, is known to stabilize the trimeric structure and create a more open conformation of the RBD allowing a stronger binding with ACE228. In paradox, our analysis shows that the majority of the novel or rare spike mutations (<0.2% prevalence in the total sequenced samples, Table 1) in Omicron may have deleterious effect on the host interactions owing to their presence at the constrained RBD regions in terms of ACE2 binding (10/15) and/or RBD expression (8/15) (Table 2). Notably, most of the spike mutations which predicted a favorable effect on ACE2 binding or RBD expression or both are present in current VOCs, primarily Delta (T478K), and Alpha (N501Y), Beta (K417N) variants. Further, a set of mutations in Omicron which are present inside (P681H) or in vicinity (D614G, H655Y) of the furin cleavage site (FCS) of SARS-CoV-2 spike protein — a small stretch of peptide (PRRAR) inserted at the intersection of spike segments S1 and S2 (681–685 aa residues)— can enhance proteolytic cleavage of spike protein by a host protease-furin, which is said to improve its fusion to the host-cell membrane29. P681H is characteristically present in multiple VOCs/VOIs such as B.1.1.7, P.1, Q.1 and B.1.621 lineage variants30. A mutation at the same location P681R has been present in Delta variant and its emerging sub-lineages31. The characterization of the individual mutations present on RBM gives a speculation that Omicron may not have more efficient interactions with the host than existing VOCs/VOIs, more specifically Delta. A further assessment of the allosteric influence and dynamic interactions of the mutations present at RBD and other regions of spike protein, and examining in *in situ*/*vivo* studies will be necessary to understand their exact impact on host-receptor binding and its clinical correlates. The clinical data on severity of the illness in Omicron is currently very limited, although the evaluation of early cases has indicated a milder illness compared to the existing VOCs32. ### Viral replication Many of the mutations, especially in the non-spike regions, are linked with viral oligomerization and synthesis and packaging of the ribonucleic acid core (Table 1) and it is likely that they may have a role in replication of the virus inside the host cells33. Particular to notice has been NSP12 P323L, a mutation present in *RNA-dependent RNA polymerase (RdRp)* coding region (Fig. 1, Table 1). Of note this has been a frequently observed mutation in the earlier variants (96.69%) (Table 1). However, whether these mutations will have a positive or a negative impact on the viral replication is little understood yet. Interestingly, the results of a preliminary study34 which employed *ex-vivo* cultures of Omicron isolated from the respiratory tract of the infected patients, and compared that with WT and Delta variant strains. The authors observed that after 24 hours of incubation Omicron replicated 70 times faster than WT and Delta variant strains in the human bronchus. In contrast, it replicated less efficiently (>10 times lower) in the human lung tissue than WT strain, and also lower than Delta variant. ### Immune escape Majority of the spike mutations (18/32) in Omicron have occurred at the known antibody recognition sites (Table 1). Existing studies have established role of these mutations in immune escape against convalescent sera, vaccine acquired antibodies, and therapeutically used monoclonal antibodies (Table 1). The emerging preliminary evidence from *in situ* studies indicate potential immune escape by Omicron against convalescent sera, vaccine acquired antibodies, and therapeutically used monoclonal antibodies 23,35,36. Interestingly, Omicron contains K417N and E484A, which are present in multiple existing variants, and are believed to contribute to immune escape28. Of note, K417 locus is a known epitope for CB26—a therapeutically used monoclonal antibody in COVID-1928. A larger number of mutations in Omicron spike protein, more specifically in RBD may be an evolutionary gain in this variant providing it higher immune escape. A support for this notion comes from a recent study by Nabel *et al*37 who demonstrated that SARS-CoV-2 pseudotypes containing up to seven mutations, as opposed to the one to three found in earlier variants of concern, are more resistant to neutralization by therapeutic antibodies and serum from vaccine recipients. A non-spike mutation (N) P13L present in Omicron (Table 1) is shown to cause complete loss of recognition by epitope-specific (B*27:05-restricted CD8+ nucleocapsid epitope QRNAPRITF9-17) T cells in a recent cell-line based *in situ* study22, however, no such evidence in human samples are currently available. In a recent preliminary study Redd *et al* examined peripheral blood mononuclear cell (PBMC) samples from PCR-confirmed, recovered/ convalescent COVID-19 cases (n=30) for their anti-SARS-CoV-2 CD8+ T-cell responses with Omicron. Authors noted that only one low-prevalence (found in 7%) epitope (GVYFASTEK, restricted to HLA*A03:01 and HLA*A11:01) from the spike protein (T95I) region was mutated in Omicron38. Presence of these mutations raises concerns about escaping T cell immunity by Omicron39 and hence are required to be explored in detail. Currently *in situ/vivo* studies are scarce which have examined Omicron response against current COVID-19 vaccines in use, although, emerging preliminary evidence is supporting the prediction of immune escape for this variant23,35,36,40. Cele *et al*23 tested the ability of plasma from 14 BNT162b2 vaccinated study participants to neutralize Omicron versus WT D614G virus in a live virus neutralization assay.The authors observed that Omicron showed a 41-fold decline in FRNT50 Geometric mean titer (GMT) compared to WT D614G virus in the subjects without previous infection (6/14). Interestingly, those with previous infection showed relatively higher neutralization titers with Omicron (6/14), which indicated that the previous infection, followed by vaccination or booster may increase the levels of neutralization and protection from severe disease in cases of Omicron infection. ### Epidemiological correlates: Omicron vs. Delta variants The analysis of the SARS-CoV-2 genomic sequences from SA indicates the Omicron has gained an advantage in terms of transmissibility over Delta variant (Fig. 4). A third COVID-19 wave driven by Delta variant has recently occurred in SA41 hence current epidemiological characteristics of Delta and newly emerged variant in the local population should be analyzed in this backdrop. Our analyses showed that before the arrival of Omicron, Delta variant was dominant locally, on the contrary, at present majority of the new sequences are from Omicron (Omicron vs. Delta: r=-0.99, p< .001, 95% CI: -0.99 to -0.97) (Fig. 4). The steep rise in the new COVID-19 cases in SA are seem to be driven by Omicron, in contrast, Delta variant linked cases are seeing a decline (Fig. 4). The rapid rise in the new COVID-19 cases linked with the emergence of a new SARS-CoV-2 variant is strongly indicating the commencement of a new COVID-19 wave in SA42. Further, the death, which is considered a strong indicator of the virulence/lethality, is showing a negative correlation (r=-0.04, p>0.05, 95% CI =-0.52 to 0.58) (Fig. 4) with rise in Omicron, however, it correlates positively with Delta variant in the period post-emergence (r=0.38, p=0.065, 95% CI= - 0.025 to 0.68) in the complete duration of the study. The pattern indicates that current deaths are primarily linked with Delta rather than Omicron. At present, the studies are scarce which have examined transmissibility advantage for Omicron. A ∼2.4 (2.0-2.7) times higher transmissibility with Omicron compared to Delta variant in Gauteng, SA was indicated in a recent conference presentation by South African COVID-19 Modelling Consortium43. Another preliminary estimate from UK has indicated that with Omicron risk of spreading infection to the members in a household is three times higher compared to Delta variant44. Based on the current epidemiological patterns observed in SA, an epidemiological advantage to Omicron in comparison to Delta can be inferred in terms of transmissibility, however, our analysis has not indicated for increased lethality with Omicron compared to Delta and other variants currently circulating in SA population. However, it will be too early to assess the virulence/lethality with Omicron and a study of epidemiological patterns over a longer time period will be mandatory to develop a conclusive view on this. It is noteworthy that the presence of an immunological barrier in the population imparted by recent COVID-19 wave mediated by Delta variant may be a likely reason for fall in new cases by this variant. A continuous fall in Delta cases is also noticeable in the period prior to emergence of Omicron (Fig. 4), which is further substantiating this notion. The data records show that a significant proportion of the local population in SA has been fully vaccinated by now (25.2%)45. Notably, the high number of immune escape mutations in Delta has contributed to lowered efficacy of the vaccines, immunity from the natural infections, and therapeutically used antibodies28. As Omicron contains much higher number of immune escape related mutations, including many shared with Delta (Fig. 3), it is likely that Omicron may have added potential for the vaccine breakthorugh infections and reinfections. Any supportive evidence in literature which examined immune escape capability of Omicron in epidemiological data is currently scarce, however, similar speculations have been presented by few other authors and the global health regulatory bodies2,32,46. Interestingly, a higher risk of reinfections with Omicron has been indicated in a recent preprint article by Pulliam *et al*47. The authors of this study have performed a retrospective analysis of routine epidemiological surveillance data to examine whether SARS-CoV-2 reinfection risk has changed through time in SA, in the context of the emergence of the consecutive variants: Beta, Delta, and Omicron. The authors noted that as compared to the first wave driven by WT strains subsequent waves by Beta and Delta variant had lower estimated hazard ratio for reinfection versus primary infection [relative hazard ratio for wave 2 vs. wave 1: 0.75 (CI95: 0.59-0.97); for wave 3 vs. wave 1:0.71 (CI95: 0.56-0.92)] in comparison to Omicron [Omicron surge vs. wave 1: 2.39 (95%CI: 1.88-3.11)]. ### Conclusion ***In silico*** analysis of the viral genomic sequences suggest that the Omicron may have greater immune escape ability than the existing VOCs/VOIs including Delta. There are no clear indications that Omicron may have higher virulence/lethality than other, so far reported variants. Therefore, the higher ability for immune escape can be a likely reason for the recent surge in Omicron cases in SA. ## Limitations The limitations of this study needs to be considered while interpreting its clinical and epidemiological significance. Firstly, as of now very limited number of genomic sequences and epidemiological data from the geographical regions affected with Omicron is available for the analysis, which may have an impact on the results. It is likely that the relative frequency of certain lineage characterising mutations in the Omicron variant may vary in future, as more number of sequences are accumulated over time. The epidemiological correlates for the new variant are only indicative of current preliminary trends and requires to be followed for a longer period, more specifically the deaths. ## Supporting information Figs. S1-3, and Table S1 [[supplements/267908_file02.pdf]](pending:yes) ## Data Availability The data generated in this study has been included with the submission. ## Data sharing Primary data used for this study are publicly available on: SARS-CoV-2 genomic sequence—GISAID database: [https://www.gisaid.org/](https://www.gisaid.org/); Epidemiological data—Worldometer: [https://www.worldometers.info/coronavirus/coronavirus/country/india](https://www.worldometers.info/coronavirus/coronavirus/country/india)). The categorized data for the study period can be availed from the corresponding author on reasonable request. ## Author (s) contributions AK, GK, and PD collected samples and analyzed data. AK wrote first draft. AA and HS performed statistical analysis. MF, SK, RKN, RKJ, CS, MK, PP, KS, KK, and SNP reviewed and edited the paper. All authors consented for submitting final draft. ## Ethical statement An ethical clearance was precluded as *in silico* methods and open access/globally shared databases were used to conduct this study. ## Financial support None ## Conflict of Interest None ## Acknowledgements The study has used SARS-CoV-2 genomic sequence and epidemiological data from GISAID ([https://www.gisaid.org/](https://www.gisaid.org/)), outbreak.info ([http://outbreak.info/](http://outbreak.info/)), and Worldometer ([https://www.worldometers.info/coronavirus/coronavirus/country/india](https://www.worldometers.info/coronavirus/coronavirus/country/india)), respectively. * Received December 18, 2021. * Revision received December 18, 2021. * Accepted December 22, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.GISAID - hCov19 Variants. Accessed July 20, 2021. [https://www.gisaid.org/hcov19-variants/](https://www.gisaid.org/hcov19-variants/) 2. 2.Update on Omicron. Accessed December 4, 2021. [https://www.who.int/news/item/28-11-2021-update-on-omicron](https://www.who.int/news/item/28-11-2021-update-on-omicron) 3. 3.Callaway E. Heavily mutated coronavirus variant puts scientists on alert. Nature. Published online November 25, 2021. doi:10.1038/D41586-021-03552-W [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/D41586-021-03552-W&link_type=DOI) 4. 4.Tracking SARS-CoV-2 variants. Accessed December 16, 2021. [https://www.who.int/en/activities/tracking-SARS-CoV-2-variants/](https://www.who.int/en/activities/tracking-SARS-CoV-2-variants/) 5. 5.Li Q, Wu J, Nie J, et al. The Impact of Mutations in SARS-CoV-2 Spike on Viral Infectivity and Antigenicity. Cell. 2020;182(5):1284–1294.e9. doi:10.1016/j.cell.2020.07.012 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.07.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32730807&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F22%2F2021.12.18.21267908.atom) 6. 6.Shen L, Bard JD, Triche TJ, Judkins AR, Biegel JA, Gai X. Rapidly emerging SARS-CoV-2 B.1.1.7 sub-lineage in the United States of America with spike protein D178H and membrane protein V70L mutations. Emerg Microbes Infect. 2021;10(1):1293–1299. doi:10.1080/22221751.2021.1943540 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/22221751.2021.1943540&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34125658&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F22%2F2021.12.18.21267908.atom) 7. 7.Starr TN, Greaney AJ, Hilton SK, et al. Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding. Cell. 2020;182(5):1295–1310.e20. doi:10.1016/j.cell.2020.08.012 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.08.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32841599&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F22%2F2021.12.18.21267908.atom) 8. 8.Greaney AJ, Starr TN, Gilchuk P, et al. Complete Mapping of Mutations to the SARS-CoV-2 Spike Receptor-Binding Domain that Escape Antibody Recognition. Cell Host Microbe. 2021;29(1):44–57.e9. doi:10.1016/J.CHOM.2020.11.007 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.CHOM.2020.11.007&link_type=DOI) 9. 9.Ho D, Wang P, Liu L, et al. Increased Resistance of SARS-CoV-2 Variants B.1.351 and B.1.1.7 to Antibody Neutralization. Res Sq. Published online January 29, 2021. doi:10.21203/RS.3.RS-155394/V1 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.21203/RS.3.RS-155394/V1&link_type=DOI) 10. 10.Wu J, Zhang L, Zhang Y, et al. The Antigenicity of Epidemic SARS-CoV-2 Variants in the United Kingdom. Front Immunol. 2021;12. doi:10.3389/FIMMU.2021.687869 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/FIMMU.2021.687869&link_type=DOI) 11. 11.Greaney AJ, Loes AN, Crawford KHD, et al. Comprehensive mapping of mutations in the SARS-CoV-2 receptor-binding domain that affect recognition by polyclonal human plasma antibodies. Cell Host Microbe. 2021;29(3):463–476.e6. doi:10.1016/J.CHOM.2021.02.003 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.CHOM.2021.02.003&link_type=DOI) 12. 12.Wu K, Peng G, Wilken M, Geraghty RJ, Li F. Mechanisms of host receptor adaptation by severe acute respiratory syndrome coronavirus. J Biol Chem. 2012;287(12):8904–8911. doi:10.1074/JBC.M111.325803 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjExOiIyODcvMTIvODkwNCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzIyLzIwMjEuMTIuMTguMjEyNjc5MDguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 13. 13.Wibmer CK, Ayres F, Hermanus T, et al. SARS-CoV-2 501Y.V2 escapes neutralization by South African COVID-19 donor plasma. Nat Med. 2021;27(4):622–625. doi:10.1038/S41591-021-01285-X [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/S41591-021-01285-X&link_type=DOI) 14. 14.Korber B, Fischer WM, Gnanakaran S, et al. Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus. Cell. 2020;182(4):812–827.e19. doi:10.1016/j.cell.2020.06.043 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.06.043&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32697968&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F22%2F2021.12.18.21267908.atom) 15. 15.Braun KM, Moreno GK, Halfmann PJ, et al. Transmission of SARS-CoV-2 in domestic cats imposes a narrow bottleneck. bioRxiv Prepr Serv Biol. 2021;17(2):e1009373. doi:10.1101/2020.11.16.384917 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4xMS4xNi4zODQ5MTd2MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzIyLzIwMjEuMTIuMTguMjEyNjc5MDguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 16. 16.Dieterle ME, Haslwanter D, Bortz RH, et al. A Replication-Competent Vesicular Stomatitis Virus for Studies of SARS-CoV-2 Spike-Mediated Cell Entry and Its Inhibition. Cell Host Microbe. 2020;28(3):486–496.e6. doi:10.1016/J.CHOM.2020.06.020 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.CHOM.2020.06.020&link_type=DOI) 17. 17.Baum A, Fulton BO, Wloga E, et al. Antibody cocktail to SARS-CoV-2 spike protein prevents rapid mutational escape seen with individual antibodies. Science. 2020;369(6506):1014–1018. doi:10.1126/SCIENCE.ABD0831 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUwNi8xMDE0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMjIvMjAyMS4xMi4xOC4yMTI2NzkwOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 18. 18.Sabir DK. Analysis of SARS-COV2 spike protein variants among Iraqi isolates. Gene reports. 2022;26:101420. doi:10.1016/J.GENREP.2021.101420 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.GENREP.2021.101420&link_type=DOI) 19. 19.Voss C, Esmail S, Liu X, et al. Epitope-specific antibody responses differentiate COVID-19 outcomes and variants of concern. JCI insight. 2021;6(13). doi:10.1172/JCI.INSIGHT.148855 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI.INSIGHT.148855&link_type=DOI) 20. 20.Lubinski B, Fernandes MH V., Frazier L, et al. Functional evaluation of proteolytic activation for the SARS-CoV-2 variant B.1.1.7: role of the P681H mutation. bioRxiv Prepr Serv Biol. Published online July 29, 2021. doi:10.1101/2021.04.06.438731 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMS4wNC4wNi40Mzg3MzF2MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzIyLzIwMjEuMTIuMTguMjEyNjc5MDguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 21. 21.Ramirez S, Fernandez-Antunez C, Galli A, et al. Overcoming Culture Restriction for SARS-CoV-2 in Human Cells Facilitates the Screening of Compounds Inhibiting Viral Replication. Antimicrob Agents Chemother. 2021;65(7). doi:10.1128/AAC.00097-21 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjE0OiJBQUMuMDAwOTctMjF2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzIyLzIwMjEuMTIuMTguMjEyNjc5MDguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 22. 22.de Silva TI, Liu G, Lindsey BB, et al. The impact of viral mutations on recognition by SARS-CoV-2 specific T cells. iScience. 2021;24(11):103353. doi:10.1016/J.ISCI.2021.103353 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.ISCI.2021.103353&link_type=DOI) 23. 23.SARS-Cov-2 Omicron Has Extensive But Incomplete Escape Of Pfizer BNT162b2 Elicited Neutralization And Requires ACE2 For Infection - allAfrica.com. Accessed December 8, 2021. [https://allafrica.com/view/resource/main/main/id/00130592.html](https://allafrica.com/view/resource/main/main/id/00130592.html) 24. 24.Avanzato VA, Matson MJ, Seifert SN, et al. Case Study: Prolonged Infectious SARS-CoV-2 Shedding from an Asymptomatic Immunocompromised Individual with Cancer. Cell. 2020;183(7):1901–1912.e9. doi:10.1016/J.CELL.2020.10.049 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.CELL.2020.10.049&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F22%2F2021.12.18.21267908.atom) 25. 25.Karim F, Moosa M, Gosnell B, et al. Persistent SARS-CoV-2 infection and intra-host evolution in association with advanced HIV infection. medRxiv. Published online June 4, 2021:2021.06.03.21258228. doi:10.1101/2021.06.03.21258228 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wNi4wMy4yMTI1ODIyOHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMjIvMjAyMS4xMi4xOC4yMTI2NzkwOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 26. 26.Choi B, Choudhary MC, Regan J, et al. Persistence and Evolution of SARS-CoV-2 in an Immunocompromised Host. N Engl J Med. 2020;383(23):2291–2293. doi:10.1056/NEJMC2031364 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2031364&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F22%2F2021.12.18.21267908.atom) 27. 27.Voloch CM, da Silva Francisco Jr R, P de Almeida LG, et al. Intra-host evolution during SARS-CoV-2 prolonged infection. Virus Evol. 2021;7(2):1–11. doi:10.1093/VE/VEAB078 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/VE/VEAB078&link_type=DOI) 28. 28.Kumar A, Rakesh Parashar |, Kumar S, et al. Emerging SARS-CoV-2 variants can potentially break set epidemiological barriers in COVID-19. J Med Virol. Published online November 29, 2021. doi:10.1002/JMV.27467 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/JMV.27467&link_type=DOI) 29. 29.Kumar A, Prasoon P, Kumari C, et al. SARS-CoV-2-specific virulence factors in COVID-19. J Med Virol. Published online 2020. doi:10.1002/jmv.26615 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.26615&link_type=DOI) 30. 30. Alaa Abdel Latif, Julia L. Mullen, et al. P681H mutation Report. outbreak.info. Accessed December 8, 2021. [https://outbreak.info/situation-reports?muts=S%3AP681H](https://outbreak.info/situation-reports?muts=S%3AP681H) 31. 31. Alaa Abdel Latif, Julia L. Mullen, et al. P681R mutation Report. outbreak.info, available at outbreak.info. Accessed December 8, 2021. [https://outbreak.info/situation-reports?muts=S%3AP681R](https://outbreak.info/situation-reports?muts=S%3AP681R) 32. 32.Threat Assessment Brief: Implications of the further emergence and spread of the SARS CoV 2 B.1.1.529 variant of concern (Omicron) for the EU/EEA first update. Accessed December 5, 2021. [https://www.ecdc.europa.eu/en/publications-data/covid-19-threat-assessment-spread-omicron-first-update](https://www.ecdc.europa.eu/en/publications-data/covid-19-threat-assessment-spread-omicron-first-update) 33. 33.Naqvi AAT, Fatima K, Mohammad T, et al. Insights into SARS-CoV-2 genome, structure, evolution, pathogenesis and therapies: Structural genomics approach. Biochim Biophys Acta - Mol Basis Dis. 2020;1866(10):165878. doi:10.1016/j.bbadis.2020.165878 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbadis.2020.165878&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32544429&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F22%2F2021.12.18.21267908.atom) 34. 34.HKUMed finds Omicron SARS-CoV-2 can infect faster and better than Delta in human bronchus but with less severe infection in lung. [https://www.med.hku.hk/en/news/press/20211215-omicron-sars-cov-2-infection?](https://www.med.hku.hk/en/news/press/20211215-omicron-sars-cov-2-infection?) 35. 35.Wilhelm A, Widera M, Grikscheit K, et al. Reduced Neutralization of SARS-CoV-2 Omicron Variant by Vaccine Sera and monoclonal antibodies. medRxiv. Published online December 8, 2021:2021.12.07.21267432. doi:10.1101/2021.12.07.21267432 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4xMi4wNy4yMTI2NzQzMnY0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMjIvMjAyMS4xMi4xOC4yMTI2NzkwOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 36. 36.South African COVID-19 Modelling Consortium. ShewardEtAl.pdf - Google Drive. Accessed December 10, 2021. [https://drive.google.com/file/d/1CuxmNYj5cpIuxWXhjjVmuDqntxXwlfXQ/view](https://drive.google.com/file/d/1CuxmNYj5cpIuxWXhjjVmuDqntxXwlfXQ/view) 37. 37.Nabel KG, Clark SA, Shankar S, et al. Structural basis for continued antibody evasion by the SARS-CoV-2 receptor binding domain. doi:10.1126/SCIENCE.ABL6251 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/SCIENCE.ABL6251&link_type=DOI) 38. 38.Redd AD, Nardin A, Kared H, et al. Minimal cross-over between mutations associated with Omicron variant of SARS-CoV-2 and CD8+ T cell epitopes identified in COVID-19 convalescent individuals. bioRxiv. Published online December 9, 2021:2021.12.06.471446. doi:10.1101/2021.12.06.471446 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMS4xMi4wNi40NzE0NDZ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzIyLzIwMjEuMTIuMTguMjEyNjc5MDguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 39. 39.Grifoni A, Sidney J, Vita R, et al. SARS-CoV-2 human T cell epitopes: Adaptive immune response against COVID-19. Cell Host Microbe. 2021;29(7):1076–1092. doi:10.1016/J.CHOM.2021.05.010 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.CHOM.2021.05.010&link_type=DOI) 40. 40.Public Health England (PHE). SARS-CoV-2 variants of concern and variants under investigation in England Variant of concern: Omicron, VOC21NOV-01 (B.1.1.529) Technical briefing 30. [https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment\_data/file/1038404/Technical\_Briefing\_30.pdf](https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/1038404/Technical_Briefing_30.pdf) 41. 41.Tegally H, Wilkinson E, Althaus CL, et al. Rapid replacement of the Beta variant by the Delta variant in South Africa. medRxiv. Published online September 27, 2021:2021.09.23.21264018. doi:10.1101/2021.09.23.21264018 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wOS4yMy4yMTI2NDAxOHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMjIvMjAyMS4xMi4xOC4yMTI2NzkwOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 42. 42.Kumar A, Dwivedi P, Kumar G, et al. Second wave of COVID-19 in India could be predicted with genomic surveillance of SARS-CoV-2 variants coupled with epidemiological data: A tool for future. medRxiv. Published online June 13, 2021:2021.06.09.21258612. doi:10.1101/2021.06.09.21258612 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wNi4wOS4yMTI1ODYxMnYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMjIvMjAyMS4xMi4xOC4yMTI2NzkwOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 43. 43.Pearson\_Epidemics8\_Omicron.pdf - Google Drive. Accessed December 4, 2021. [https://drive.google.com/file/d/1hA6Mec2Gq3LGqTEOj35RqSeAb\_SmXpbI/view](https://drive.google.com/file/d/1hA6Mec2Gq3LGqTEOj35RqSeAb_SmXpbI/view) 44. 44.Public Health England (PHE). SARSCoV2 variants of concern and variants under investigation in England. Technical briefing 31, 10 December 2021. [https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment\_data/file/1040076/Technical\_Briefing\_31.pdf](https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/1040076/Technical_Briefing_31.pdf) 45. 45.Coronavirus (COVID-19) Vaccinations - Statistics and Research - Our World in Data. Accessed December 6, 2021. [https://ourworldindata.org/covid-vaccinations?country=~ZAF](https://ourworldindata.org/covid-vaccinations?country=%E2%88%BCZAF)) 46. 46.Science Brief: Omicron (B.1.1.529) Variant | CDC. Accessed December 5, 2021. [https://www.cdc.gov/coronavirus/2019-ncov/science/science-briefs/scientific-brief-omicron-variant.html#6](https://www.cdc.gov/coronavirus/2019-ncov/science/science-briefs/scientific-brief-omicron-variant.html#6) 47. 47.Pulliam JRC, Schalkwyk C van, Govender N, et al. Increased risk of SARS-CoV-2 reinfection associated with emergence of the Omicron variant in South Africa. medRxiv. Published online December 2, 2021:2021.11.11.21266068. doi:10.1101/2021.11.11.21266068 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4xMS4xMS4yMTI2NjA2OHYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMjIvMjAyMS4xMi4xOC4yMTI2NzkwOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)