Characterization of altered genomic landscape of SARS-CoV-2 variants isolated in Saudi Arabia in a comparative in silico study

SARS-CoV-2 has become one of the unprecedented global health challenge for human population. Genomic signature studies of SARS-CoV-2 reveals relation between geographical location of the isolates and genetic diversity. The present work is an in silico, cross sectional study aimed to determine the genetic heterogeneity of SARS-CoV-2 variants isolated in Saudi Arabia compared to the first isolated strain NC_045512 (reference sequence). Each sequence was aligned against the reference sequence using local alignment search tool (NCBI) Nucleotide-BLAST. A total of 58 SARS-CoV-2 genomes were isolated in KSA and retrieved from NCBI. Our study shows that KSA variants demonstrated homology ranging between 99.96 and 99.98 % compared to the reference strain. There are 89 nucleotide changes that have been identified among the KSA variants; the most common nucleotide change was C: T accounting for 50.6% (45/89). These nucleotides changes resulted in 53.9% (48/89) missense mutations and 42.7% (38/89) silent mutations; while the majority of mutations- 48.3% (43/89) occurred in ORF1ab gene. All structural genes displayed mutations; N gene harbored 16.9% (15/89) mutations, S gene displayed 15.7% (14/89) mutations, M gene exhibited 2.2% (2/89) mutations and E gene showed only 1 mutation which was silent. The most frequently changed nucleotide was C3037T (silent mutation) and A23403G (D614G), each of which occurred in 57 variants out of 58 followed by C14408T (P4715L) and C241T (5′UTR) which were found in 56 and 55 variants respectively. The Phylogenetic trees showed that SARS-CoV-2 variants isolated in Saudi Arabia clustered together closely.