Global and Local Mutations in Bangladeshi SARS-CoV-2 Genomes

Corona Virus Disease-2019 (COVID-19) warrants comprehensive investigations of publicly available Severe Acute Respiratory Syndrome-CoronaVirus-2 (SARS-CoV-2) genomes to gain new insight about their epidemiology, mutations and pathogenesis. Nearly 0.4 million mutations were identified so far in ∼60,000 SARS-CoV-2 genomic sequences. In this study, we compared 207 of SARS-CoV-2 genomes reported from different parts of Bangladesh and their comparison with 467 globally reported sequences to understand the origin of viruses, possible patterns of mutations, availability of unique mutations, and their apparent impact on pathogenicity of the virus in victims of Bangladeshi population. Phylogenetic analyses indicates that in Bangladesh, SARS-CoV-2 viruses might arrived through infected travelers from European countries, and the GR clade was found as predominant in this region. We found 2602 mutations including 1602 missense mutations, 612 synonymous mutations, 36 insertions and deletions with 352 other mutations types. In line with the global trend, D614G mutation in spike glycoprotein was predominantly high (95.6%) in Bangladeshi isolates. Interestingly, we found the average number of mutations in ORF1ab, S, ORF3a, M and N of genomes, having nucleotide shift at G614 (n=459), were significantly higher (p≤0.001) than those having mutation at D614 (n=215). Previously reported frequent mutations such as P4715L, D614G, R203K, G204R and I300F were also prevalent in Bangladeshi isolates. Additionally, 87 unique amino acid changes were revealed and were categorized as originating from different cities of Bangladesh. The analyses would increase our understanding of variations in virus genomes circulating in Bangladesh and elsewhere and help develop novel therapeutic targets against SARS-CoV-2.