Omicron variant genome evolution and phylogenetics

Abstract

Following the discovery of the SARS-CoV-2 Omicron variant (B.1.1.529), the global COVID-19 outbreak has resurfaced after appearing to be relentlessly spreading over the past 2 years. This new variant showed marked degree of mutation, compared with the previous SARS-CoV-2 variants. This study investigates the evolutionary links between Omicron variant and recently emerged SARS-CoV-2 variants. The entire genome sequences of SARS-CoV-2 variants were obtained, aligned using Clustal Omega, pairwise comparison was computed, differences, identity percent, gaps, and mutations were noted, and the identity matrix was generated. The phylogenetics of Omicron variants were determined using a variety of evolutionary substitution models. The ultrametric and metric clustering methods, such as UPGMA and neighbor-joining (NJ), using nucleotide substitution models that allowed the inclusion of nucleotide transitions and transversions as Kimura 80 models, revealed that the Omicron variant forms a new monophyletic clade that is distant from other SARS-CoV-2 variants. In contrast, the NJ method using a basic nucleotide substitution model such as Jukes–Cantor revealed a close relationship between the Omicron variant and the recently evolved Alpha variant. Based on the percentage of sequence identity, the closest variants were in the following order: Omicron, Alpha, Gamma, Delta, Beta, Mu, and then the SARS-CoV-2 USA isolate. A genome alignment with other variants indicated the greatest number of gaps in the Omicron variant's genome ranging from 43 to 63 gaps. It is possible, given their close relationship to the Alpha variety, that Omicron has been around for much longer than predicted, even though they created a separate monophyletic group. Sequencing initiatives in a systematic and comprehensive manner is highly recommended to study the evolution and mutations of the virus.