Emergence of Omicron FN.1 a descendent of BQ.1.1 in Botswana

Botswana, like the rest of the world, has been significantly impacted by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In December 2022, we detected a monophyletic cluster of genomes comprising a sublineage of the Omicron variant of concern (VOC) designated as B.1.1.529.5.3.1.1.1.1.1.1.74.1 (alias FN.1, clade 22E). These genomes were sourced from both epidemiologically linked and unlinked samples collected in three close locations within the district of Greater Gaborone. In this study, we assessed the worldwide prevalence of the FN.1 lineage, evaluated its mutational profile, and conducted a phylogeographic analysis to reveal its global dispersal dynamics. Among approximately 16 million publicly available SARS-CoV-2 sequences generated by 30 September 2023, only 87 were of the FN.1 lineage, including 22 from Botswana, 6 from South Africa, and 59 from the UK. The estimated time to the most recent common ancestor of the 87 FN.1 sequences was 22 October 2022 [95% highest posterior density: 2 September 2022-24 November 2022], with the earliest of the 22 Botswana sequences having been sampled on 7 December 2022. Discrete trait reconstruction of FN.1 identified Botswana as the most probable place of origin. The FN.1 lineage is derived from the BQ.1.1 lineage and carries two missense variants in the spike protein, S:K182E in NTD and S:T478R in RDB. Among the over 90 SARS-CoV-2 lineages circulating in Botswana between September 2020 and July 2023, FN.1 was most closely related to BQ.1.1.74 based on maximum likelihood phylogenetic inference, differing only by the S:K182E mutation found in FN.1. Given the early detection of numerous novel variants from Botswana and its neighbouring countries, our study underscores the necessity of continuous surveillance to monitor the emergence of potential VOCs, integrating molecular and spatial data to identify dissemination patterns enhancing preparedness efforts.