key: cord-0880145-2cgnlitt authors: Grabowski, F.; Kochanczyk, M.; Lipniacki, T. title: Omicron strain spreads with the doubling time of 3.2-3.6 days in South Africa province of Gauteng that achieved herd immunity to Delta variant date: 2021-12-09 journal: nan DOI: 10.1101/2021.12.08.21267494 sha: 9df4b776155c421a9db85255901f57c0cac83f52 doc_id: 880145 cord_uid: 2cgnlitt Omicron, the novel, highly mutated SARS-CoV-2 Variant of Concern (belonging to the Pango lineage B.1.1.529), was first collected on November 8, 2021, in Gauteng province of South Africa. By the end of November 2021 it has spread towards fixation in Gauteng and was detected on all continents. Based on data collected till December 7, 2021, we showed the exponential growth of the Omicron variant over the four-week period in Gauteng (November 8-December 5, 2021) with the doubling time equal 3.38 day [CI 95%: 3.18-3.61 day]. Log-linear regression suggests that the spread began around October 10, 2021, however due to stochasticity in the initial spread this estimate is likely inaccurate. Phylogenetic analysis indicates that the Omicron strain started to diverge in between October 28 and November 5, 2021. This implies that the hidden spread of Omicron before October 10, 2021 (which would suggest slower strain growth) is unlikely. The very short doubling time of Omicron in Gauteng, a province that has reached herd immunity to the Delta variant (implied by the decrease of the weekly number of cases between July and October, 2021, at no significant mobility restrictions), suggests that Omicron will cause abrupt outbreaks of COVID-19 epidemics across the world, and will become the (temporarily) dominant strain. Omicron, the novel SARS-CoV-2 Variant of Concern (Pango lineage B.1.1.529, Nextstrain clade identifier 21K, first collected in Gauteng province of South Africa, November 8, 2021 GISAID sequence accession ID: EPI_ISL_6913995) is characterized by 30 amino acid substitutions, three small deletions and one small insertion in the spike protein compared to the original SARS-CoV-2 virus 1 . Altogether, Omicron has 51 amino-acid mutations, and its closest known epigenetic sister has 15 mutations (GISAID sequence accession ID: EPI_ISL_622806) with only 9 common mutations, implying the distance of 42 mutations from the last common ancestor (based on the phylogenetic tree generated by Nextstrain). The collection date of the Omicron epigenetic sister, September 13, 2020, suggests more than a year of evolution in an isolated niche, possibly in an immunocompromised host, but more data is necessary to rule out or confirm the existence of hidden branches (see Ref. 2 for discussion). The lineage started spreading rapidly in Gauteng province in November 2021, reaching fixation by the end of that month and causing an abrupt epidemic outbreak in the province. Based on GISAID data 3 and Gauteng province data on daily COVID-19 cases 4 , both accessed on December 7, 2021, we estimated the Omicron strain doubling time, and the strain divergence date. The 2 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted December 9, 2021. ; https://doi.org/10.1101/2021.12.08.21267494 doi: medRxiv preprint epidemic outbreak (Fig. 1A) . In Fig. 1C we show the exponential growth of Omicron variant in weeks 45-48, 2021 (November 8-December 5). Omicron cases were estimated as a product of the proportion of Omicron genomes to all collected genomes in a given week and the respective weekly number of total COVID-19 confirmed cases (see Table 1 ). 3 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted December 9, 2021. We have demonstrated the exponential growth of the Omicron strain in the South African province of Gauteng in the four-week period from November 8 to December 5, 2021, with the doubling time equal 3.38 day [CI 95%: 3.18-3.61 day]. Based on the mutation accumulation profile we found that the Omicron strain started diverging between October 28 and November 5, 2021, which is after October 10, the date suggested by log-linear regression of the number of weekly cases. The discrepancy can be explained by the occurrence of a su-4 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted December 9, 2021. ; https://doi.org/10.1101/2021.12.08.21267494 doi: medRxiv preprint perspreading event or events at the beginning of November that could have caused a sharp growth of Omicron infections before the more regular exponential growth has started. This is also suggested by the fact that before the week of November 8-15, in which 61 Omicron genomes were collected worldwide, no Omicron genomes were found. The estimated date of divergence between October 28 and November 5, 2021, in conjunction with the lack of detected Omicron genomes before November 8, implies that circulation of Omicron in October or September 2021 (which could potentially explain the abrupt outbreak) is unlikely. Before the Omicron outbreak, the Delta variant was the dominant strain in Gauteng, and between July and October the COVID-19 epidemic was receding without significant mobility reduction, suggesting that the population of Gauteng has reached herd immunity to the Delta variant. The population-level immunity has been apparently overcome by the Omicron variant. Omicron accumulated more than 30 mutations in its spike protein, with 15 substitutions in the receptor binding domain (RBD, residues 319-541) alone 9 . Many of these RBD mutations are thought to decrease potency of neutralizing antibodies 10 . This is in line with the recent study which demonstrated that the spread of the Omicron variant increased substantially the hazard ratio for reinfection versus primary infection 11 . The Omicron doubling time found here to be in range 3.18-3.61 days for the Gauteng province is larger but comparable to the doubling times during the first COVID-19 pandemic outbreaks in spring 2020 that were found in range between 1.86 and 2.88 for China, Italy, France, Germany, Spain, UK, Switzerland and New York State. Such short Omicron doubling time implies that it will likely outcompete the Delta variant, become (temporarily) the dominant strain and cause epidemic outbreaks across the world. These outbreaks may be hard to control by the current vaccines (as already suggested by abrupt growths of Omi-9 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, ECDC, https://www.ecdc.europa.eu/en/publications-data/covid-19-threatassessment-spread-omicron-first-update (visited on 12/07/2021). 10 Callaway E. & Ledford H. (2021) Nature 600, 197-199. 11 Pulliam J. R. et al. (2021 ) medRxiv, 10.1101 /2021 .11.11.21266068 10.1101 /2021 5 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted December 9, 2021. ; https://doi.org/10.1101 https://doi.org/10. /2021 cron cases in Denmark 12 and UK 13 with 69% and 77% vaccination rates respectively) as well as by lockdowns due to increasing lockdown fatigue. All used data were retrieved on December 7, 2021. The Omicron strain divergence date and mutation accumulation rate were determined by the linear regression assuming Poisson distribution of the number of mutations in a give time point. The Omicron phylogenetic tree was accessed from Nextstrain (https://nextstrain.org). The Omicron weekly cases where estimated based on GISAID data and the Gauteng province COVID-19 cumulative cases (DSFSI at University of Pretoria, https://github.com/dsfsi/covid19za/tree/ master/data) corrected from 8099 to 605 for November 23, 2021 using the information from South African National Institute for Communicable Diseases (NICD) (https://www. nicd.ac.za/latest-confirmed-cases-of-covid-19-in-south-africa-23-november-2021). According to the information provided by NICD the difference results from retrospective addition of 7494 antigen tests. The mobility in Gauteng province was accessed from COVID-19 Community Mobility Reports from Google (https://www.google.com/covid19/mobility). The weekly average was calculated based on workday days. We analyzed publicly available genome datasets retrieved from GISAID, https://www.gisaid.org and Community Mobility Reports https://www.google.com/covid19/mobility/. Cases data for Gauteng can be found in daily reports from the National Institute for Communicable Diseases https: //www.nicd.ac.za as aggregated by University of Pretoria https://github.com/dsfsi/covid19za. 12 Status på omikron-varianten (B.1.1.529) pr. 08.12.21, https://www.ssi.dk/aktuelt/nyheder/2021/statuspa-omikron-varianten-b11529-pr-071221. 13 Variant of concern: Omicron, VOC-21NOV-01 (B.1.1.529), https : / / assets . publishing . service . gov . uk / government/uploads/system/uploads/attachment_data/file/1038404/Technical_Briefing_30.pdf. 6 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted December 9, 2021. ; https://doi.org/10.1101/2021.12.08.21267494 doi: medRxiv preprint Data curation and Software, Frederic Grabowski; Investigation, Frederic Grabowski and Marek Kochańczyk; Supervision, Marek Kochańczyk and Tomasz Lipniacki; Visualization, Frederic Grabowski; Writing -original draft We are very grateful to the GISAID Initiative and all its data contributors, i.e., the Authors from the originating laboratories responsible for obtaining the specimens and the submitting laboratories where genetic sequence data were generated and shared via the GISAID Initiative, on which this research is based. This study was supported by the National Science Centre Poland OPUS grant 2018/29/B/NZ2/00668.The funding agency had no role in study design, data collection, and analysis, decision to publish, or preparation of the manuscript. The authors have no competing interests.