Genetic analysis of SARS-CoV-2 isolates collected from Bangladesh: Insights into the origin, mutational spectrum and possible pathomechanism

Evaluation of Genomic Surveillance of SARS-CoV-2 Virus Isolates and Comparison of Mutational Spectrum of Variants in Bangladesh

The SARS-CoV-2-induced disease, COVID-19, remains a worldwide public health concern due to its high rate of transmission, even in vaccinated and previously infected people. In the endemic state, it continues to cause significant pathology. To elu- cidate the viral mutational changes and screen the emergence of new variants of concern, we conducted this study in Bangladesh. The viral RNA genomes extracted from 25 ran- domly collected samples of COVID-19-positive patients from March 2021 to February 2022 were sequenced using Illumina COVID Seq protocol and genomic data processing, as well as evaluations performed in DRAGEN COVID Lineage software. In this study, the percentage of Delta, Omicron, and Mauritius variants identified were 88%, 8%, and 4%, respectively. All of the 25 samples had 23,403 A>G (D614G, S gene), 3037 C>T (nsp3), and 14,408 C>T (nsp12) mutations, where 23,403 A>G was responsible for increased transmis- sion. Omicron had the highest number of unique mutations in the spike protein (i.e., sub- stitutions, deletions, and insertions), which may explain its higher transmissibility and immune-evading ability than Delta. A total of 779 mutations were identified, where 691 substitutions, 85 deletions, and 3 insertion mutations were observed. To sum up, our study will enrich the genomic database of SARS-CoV-2, aiding in treatment strategies along with understanding the virus's preferences in both mutation type and mutation site for predicting newly emerged viruses' survival strategies and thus for preparing to coun- teract them.

SARS-CoV-2 main protease (M-pro) mutational profiling: An insight into mutation coldspots

SARS-CoV-2 and the COVID-19 pandemic have marked a milestone in the history of scientific research worldwide. To ensure that treatments are successful in the mid-long term, it is crucial to characterize SARS-CoV-2 mutations, as they might lead to viral resistance. Data from >5,700,000 SARS-CoV-2 genomes available at GISAID was used to report SARS-CoV-2 mutations. Given the pivotal role of its main protease (M-pro) in virus replication, a detailed analysis of SARS-CoV-2 M-pro mutations was conducted, with particular attention to mutation-resistant residues or mutation coldspots, defined as those residues that have mutated in five or fewer genomes. 32 mutation coldspots were identified, most of which mediate interprotomer interactions or funneling interaction networks from the substrate-binding site towards the dimerization surface and vice versa. Besides, mutation coldspots were virtually conserved in all main proteases from other CoVs. Our results provide valuable information about key residues to M-pro structure that could be useful in rational target-directed drug design and establish a solid groundwork based on mutation analyses for the inhibition of M-pro dimerization, with a potential applicability to future coronavirus outbreaks.

Genomic epidemiology of early SARS-CoV-2 transmission dynamics in Bangladesh

BACKGROUND

Genomic epidemiology has helped reconstruct the global and regional movement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, there is still a lack of understanding of SARS-CoV-2 spread in some of the world's least developed countries (LDCs).

METHODS

To begin to address this disparity, we studied the transmission dynamics of the virus in Bangladesh during the country's first COVID-19 wave by analysing case reports and whole-genome sequences from all eight divisions of the country.

RESULTS

We detected > 50 virus introductions to the country during the period, including during a period of national lockdown. Additionally, through discrete phylogeographic analyses, we identified that geographical distance and population -density and/or -size influenced virus spatial dispersal in Bangladesh.

CONCLUSIONS

Overall, this study expands our knowledge of SARS-CoV-2 genomic epidemiology in Bangladesh, shedding light on crucial transmission characteristics within the country, while also acknowledging resemblances and differences to patterns observed in other nations.

COVID-19 in Bangladesh: Wave-centric assessments and mitigation measures for future pandemics

The ongoing pandemic COVID-19 caused by Severe Acute Respiratory Coronavirus-2 (SARS-CoV-2) has wreaked havoc globally by affecting millions of lives. Although different countries found the implementation of emergency measures useful to combat the viral pandemic, many countries are still experiencing the resurgence of COVID-19 cases with new variants even after following strict containment guidelines. Country-specific lessons learned from the ongoing COVID-19 pandemic can be utilized in commencing a successful battle against the potential future outbreaks. In this article, we analyzed the overall scenario of the COVID-19 pandemic in Bangladesh from Alpha to Omicron variant and discussed the demographic, political, economic, social, and environmental influences on the mitigation strategies employed by the country to combat the pandemic. We also tried to explore the preparedness and precautionary measures taken by the responsible authorities, the choice of strategies implemented, and the effectiveness of the response initiated by the government and relevant agencies. Finally, we discussed the possible strategies that might help Bangladesh to combat future COVID-19 waves and other possible pandemics based on the experiences gathered from the ongoing COVID-19 pandemic.

Transmission Dynamics and Genomic Epidemiology of Emerging Variants of SARS-CoV-2 in Bangladesh

With the progression of the global SARS-CoV-2 pandemic, the new variants have become more infectious and continue spreading at a higher rate than pre-existing ones. Thus, we conducted a study to explore the epidemiology of emerging variants of SARS-CoV-2 that circulated in Bangladesh from December 2020 to September 2021, representing the 2nd and 3rd waves. We collected new cases and deaths per million daily data with the reproduction rate. We retrieved 928 SARS-CoV-2 sequences from GISAID and performed phylogenetic tree construction and mutation analysis. Case counts were lower initially at the end of 2020, during January-February and April-May 2021, whereas the death toll reached the highest value of 1.587 per million on the first week of August and then started to decline. All the variants (α, β, δ, η) were prevalent in the capital city, Dhaka, with dispersion to large cities, such as Sylhet and Chattogram. The B.1.1.25 lineage was prevalent during December 2020, but the B.1.617.2/δ variant was later followed by the B.1.351/β variant. The phylogeny revealed that the various strains found in Bangladesh could be from numerous countries. The intra-cluster and inter-cluster communication began in Bangladesh soon after the virus arrived. The prominent amino acid substitution was D614G from December 2020 to July 2021 (93.5 to 100%). From February-April, one of the VOC's important mutations, N501Y substitution, was also estimated at 51.8%, 76.1%, and 65.1% for the α, β and γ variants, respectively. The γ variant's unique mutation K417T was detected only at 1.8% in February. Another frequent mutation was P681R, a salient feature of the δ variant, detected in June (88.2%) and July (100%). Furthermore, only one γ variant was detected during the entire second and third wave, whereas no η variant was observed in this period. This rapid growth in the number of variants identified across Bangladesh shows virus adaptation and a lack of strict quarantine, prompting periodic genomic surveillance to foresee the spread of new variants, if any, and to take preventive measures as soon as possible.

SARS-CoV-2 and Emerging Foodborne Pathogens: Intriguing Commonalities and Obvious Differences

The coronavirus disease 2019 (COVID-19) has resulted in tremendous human and economic losses around the globe. The pandemic is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a virus that is closely related to SARS-CoV and other human and animal coronaviruses. Although foodborne diseases are rarely of pandemic proportions, some of the causative agents emerge in a manner remarkably similar to what was observed recently with SARS-CoV-2. For example, Shiga toxin-producing Escherichia coli (STEC), the most common cause of hemolytic uremic syndrome, shares evolution, pathogenesis, and immune evasion similarities with SARS-CoV-2. Both agents evolved over time in animal hosts, and during infection, they bind to specific receptors on the host cell's membrane and develop host adaptation mechanisms. Mechanisms such as point mutations and gene loss/genetic acquisition are the main driving forces for the evolution of SARS-CoV-2 and STEC. Both pathogens affect multiple body organs, and the resulting diseases are not completely cured with non-vaccine therapeutics. However, SARS-CoV-2 and STEC obviously differ in the nature of the infectious agent (i.e., virus vs. bacterium), disease epidemiological details (e.g., transmission vehicle and symptoms onset time), and disease severity. SARS-CoV-2 triggered a global pandemic while STEC led to limited, but sometimes serious, disease outbreaks. The current review compares several key aspects of these two pathogenic agents, including the underlying mechanisms of emergence, the driving forces for evolution, pathogenic mechanisms, and the host immune responses. We ask what can be learned from the emergence of both infectious agents in order to alleviate future outbreaks or pandemics.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Acute Cerebellar Inflammation and Related Ataxia: Mechanisms and Pathophysiology

The cerebellum governs motor coordination and motor learning. Infection with external microorganisms, such as viruses, bacteria, and fungi, induces the release and production of inflammatory mediators, which drive acute cerebellar inflammation. The clinical observation of acute cerebellitis is associated with the emergence of cerebellar ataxia. In our animal model of the acute inflammation of the cerebellar cortex, animals did not show any ataxia but hyperexcitability in the cerebellar cortex and depression-like behaviors. In contrast, animal models with neurodegeneration of the cerebellar Purkinje cells and hypoexcitability of the neurons show cerebellar ataxia. The suppression of the Ca2+-activated K+ channels in vivo is associated with a type of ataxia. Therefore, there is a gap in our interpretation between the very early phase of cerebellar inflammation and the emergence of cerebellar ataxia. In this review, we discuss the hypothesized scenario concerning the emergence of cerebellar ataxia. First, compared with genetically induced cerebellar ataxias, we introduce infection and inflammation in the cerebellum via aberrant immunity and glial responses. Especially, we focus on infections with cytomegalovirus, influenza virus, dengue virus, and SARS-CoV-2, potential relevance to mitochondrial DNA, and autoimmunity in infection. Second, we review neurophysiological modulation (intrinsic excitability, excitatory, and inhibitory synaptic transmission) by inflammatory mediators and aberrant immunity. Next, we discuss the cerebellar circuit dysfunction (presumably, via maintaining the homeostatic property). Lastly, we propose the mechanism of the cerebellar ataxia and possible treatments for the ataxia in the cerebellar inflammation.

A Global Mutational Profile of SARS-CoV-2: A Systematic Review and Meta-Analysis of 368,316 COVID-19 Patients

Since its first detection in December 2019, more than 232 million cases of COVID-19, including 4.7 million deaths, have been reported by the WHO. The SARS-CoV-2 viral genomes have evolved rapidly worldwide, causing the emergence of new variants. This systematic review and meta-analysis was conducted to provide a global mutational profile of SARS-CoV-2 from December 2019 to October 2020. The review was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA), and a study protocol was lodged with PROSPERO. Data from 62 eligible studies involving 368,316 SARS-CoV-2 genomes were analyzed. The mutational data analyzed showed most studies detected mutations in the Spike protein (n = 50), Nucleocapsid phosphoprotein (n = 34), ORF1ab gene (n = 29), 5′-UTR (n = 28) and ORF3a (n = 25). Under the random-effects model, pooled prevalence of SARS-CoV-2 variants was estimated at 95.1% (95% CI; 93.3–96.4%; I² = 98.952%; p = 0.000) while subgroup meta-analysis by country showed majority of the studies were conducted ‘Worldwide’ (n = 10), followed by ‘Multiple countries’ (n = 6) and the USA (n = 5). The estimated prevalence indicated a need to continuously monitor the prevalence of new mutations due to their potential influence on disease severity, transmissibility and vaccine effectiveness.

Investigating the possible origin and transmission routes of SARS-CoV-2 genomes and variants of concern in Bangladesh

The COVID-19 pandemic induced by the SARS-CoV-2 virus and its variants has ravaged most countries around the world including Bangladesh. We have analyzed publicly available genomic data to understand the current COVID-19 outbreak scenario as well as the evolutionary origin and transmission routes of SARS-CoV-2 isolates in Bangladesh. All the early isolates as well as recent B.1.1.7 and B.1.351 variants had already spread across the major divisional cities of Bangladesh. A sex biasness towards male COVID-19 patient samples sequencing has been observed over female patient samples in all age-group, that could be the trend in infection rate. Phylogenetic analysis indicated a total of 13 estimated countries, including Italy, India, United Kingdom, Saudi Arabia, United Arab Emirates, Germany, Australia, New Zealand, South Africa, Democratic Republic of the Congo, United States, Russia, and Denmark, could be the possible origin introduced SARS-CoV-2 isolates in Bangladesh because of regional and intercontinental travel. Recent, B.1.1.7 variant could be imported from a total of 7 estimated countries including UK, India, Nigeria, Spain, Ireland, Australia, and Indonesia, while South Africa and the United States are the most likely sources of B.1.351 variant in Bangladesh. Based on these findings, public health strategies could be designed and implemented to reduce the local transmission of the virus.

Temporal landscape of mutational frequencies in SARS-CoV-2 genomes of Bangladesh: possible implications from the ongoing outbreak in Bangladesh

Along with intrinsic evolution, adaptation to selective pressure in new environments might have resulted in the circulatory SARS-CoV-2 strains in response to the geoenvironmental conditions of a country and the demographic profile of its population. With this target, the current study traced the evolutionary route and mutational frequency of 198 Bangladesh-originated SARS-CoV-2 genomic sequences available in the GISAID platform over a period of 13 weeks as of 14 July 2020. The analyses were performed using MEGA X, Swiss Model Repository, Virus Pathogen Resource and Jalview visualization. Our analysis identified that majority of the circulating strains strikingly differ from both the reference genome and the first sequenced genome from Bangladesh. Mutations in nonspecific proteins (NSP2-3, NSP-12(RdRp), NSP-13(Helicase)), S-Spike, ORF3a, and N-Nucleocapsid protein were common in the circulating strains with varying degrees and the most unique mutations (UM) were found in NSP3 (UM-18). But no or limited changes were observed in NSP9, NSP11, Envelope protein (E) and accessory factors (NSP7a, ORF 6, ORF7b) suggesting the possible conserved functions of those proteins in SARS-CoV-2 propagation. However, along with D614G mutation, more than 20 different mutations in the Spike protein were detected basically in the S2 domain. Besides, mutations in SR-rich region of N protein and P323L in RDRP were also present. However, the mutation accumulation showed a significant association (p = 0.003) with sex and age of the COVID-19-positive cases. So, identification of these mutational accumulation patterns may greatly facilitate vaccine development deciphering the age and the sex-dependent differential susceptibility to COVID-19.

Molecular Analysis of SARS-CoV-2 Circulating in Bangladesh during 2020 Revealed Lineage Diversity and Potential Mutations

Virus evolution and mutation analyses are crucial for tracing virus transmission, the potential variants, and other pathogenic determinants. Despite continuing circulation of the SARS-CoV-2, very limited studies have been conducted on genetic evolutionary analysis of the virus in Bangladesh. In this study, a total of 791 complete genome sequences of SARS-CoV-2 from Bangladesh deposited in the GISAID database during March 2020 to January 2021 were analyzed. Phylogenetic analysis revealed circulation of seven GISAID clades G, GH, GR, GRY, L, O, and S or five Nextstrain clades 20A, 20B, 20C, 19A, and 19B in the country during the study period. The GISAID clade GR or the Nextstrain clade 20B or lineage B.1.1.25 is predominant in Bangladesh and closely related to the sequences from India, USA, Canada, UK, and Italy. The GR clade or B.1.1.25 lineage is likely to be responsible for the widespread community transmission of SARS-CoV-2 in the country during the first wave of infection. Significant amino acid diversity was observed among Bangladeshi SARS-CoV-2 isolates, where a total of 1023 mutations were detected. In particular, the D614G mutation in the spike protein (S_D614G) was found in 97% of the sequences. However, the introduction of lineage B.1.1.7 (UK variant/S_N501Y) and S_E484K mutation in lineage B.1.1.25 in a few sequences reported in late December 2020 is of particular concern. The wide genomic diversity indicated multiple introductions of SARS-CoV-2 into Bangladesh through various routes. Therefore, a continuous and extensive genome sequence analysis would be necessary to understand the genomic epidemiology of SARS-CoV-2 in Bangladesh.

Evolution of SARS-CoV-2: Review of Mutations, Role of the Host Immune System

Since the reporting of the first cases of coronavirus in China and the publication of the first sequence of SARS-CoV-2 in December 2019, the virus has undergone numerous mutations. In Europe, the spring outbreak (March-April) was followed by a drop in the number of cases and deaths. The disease may have evolved into a milder form. The increase in PCR-positive cases in late summer 2020 did not lead to the expected increase in hospitalizations, ICU admissions, and deaths, based on the severity of the disease in the spring. This difference in disease severity could be due to factors independent of the virus or to the evolution of the virus. This review attempts to identify the mutations that have appeared since the beginning of the pandemic and their role in the temporal evolution of the pandemic. There are a cell and humoral type cross-reactivity in a large part of the population to common cold coronaviruses (HCoVs) and SARS-CoV-2. Evolutionarily important mutations and deletions have emerged in the SARS-CoV-2 genes encoding proteins that interact with the host immune system. In addition, one of the major mutations (in viral polymerase) is logically associated with a higher frequency of mutations throughout the genome. This frequency fluctuates over time and shows a peak at the time when the epidemic was most active. The rate of mutations in proteins involved in the relationship to the immune system continues to increase after the first outbreak. The cross-reactivity on the 1 hand and the viral mutations observed on the other hand could explain the evolution of the pandemic until the summer of 2020, partly due to the evolution of the virus in relation to the host immune system. The immunization campaign began in December 2020: concerns are emerging about a possible escape of the circulating variants vaccines in early 2021. These variants could also escape immunity acquired through infection with the 2020 strains.

Extensive genetic diversity with novel mutations in spike glycoprotein of severe acute respiratory syndrome coronavirus 2, Bangladesh in late 2020

In Bangladesh, coronavirus disease 2019 (COVID-19) has been highly prevalent during late 2020, with nearly 500 000 confirmed cases. In the present study, the spike (S) protein of severe acute respiratory coronavirus 2 (SARS-CoV-2) circulating in Bangladesh was genetically investigated to elucidate the diversity of mutations and their prevalence. The nucleotide sequence of the S protein gene was determined for 15 SARS-CoV-2 samples collected from eight divisions in Bangladesh, and analysed for mutations compared with the reference strain (hCoV-19/Wuhan/WIV04/2019). All the SARS-CoV-2 S genes were assigned to B.1 lineage in G clade, and individual S proteins had 1-25 mutations causing amino acid substitution/deletion. A total of 133 mutations were detected in 15 samples, with D614G being present in all the samples; 53 were novel mutations as of January 2021. On the receptor-binding domain, 21 substitutions including ten novel mutations were identified. Other novel mutations were located on the N-terminal domain (S1 subunit) and dispersed sites in the S2 subunit, including two substitutions that remove potential N-glycosylation sites. A P681R substitution adjacent to the furin cleavage site was detected in one sample. All the mutations detected were located on positions that are functionally linked to host transition, antigenic drift, host surface receptor binding or antibody recognition sites, and viral oligomerization interfaces, which presumably related to viral transmission and pathogenic capacity.

Genomic exploration light on multiple origin with potential parsimony-informative sites of the severe acute respiratory syndrome coronavirus 2 in Bangladesh

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a new strain of beta coronavirus that has spread worldwide within a short period of time and has been responsible for the current COVID-19 pandemic. This novel virus shows high transmission and adaptability frequency into the host with rapid changes in genomic sequences. In this study, we analyzed the complete genome of 41 strains isolated in Bangladesh to understand the evolutionary route and genetic variations of this rapidly evolving virus. The phylogenetics, parsimony informative sites and mutation analyses were performed using MEGA X, Multiple sequence alignment program (MAFFT), and Virus Pathogen Resource. The phylogenetic analysis of the studied genomes along with the reference genome suggested that the viral strains found in Bangladesh might be coming from multiple countries such as France, Germany, India, the USA, and Brazil. After entering into the country, intra-cluster and inter-cluster began to circulate in the 8 individual divisions of Bangladesh. We also identified 26 parsimony-informative sites along with the 9 most important sites for virus evolution. Genome-wide annotations revealed 256 mutations, of which 10 were novel (NSP3, RdRp, Spike) in Bangladeshi strains where I120F(NSP2), P323L(RdRp), D614G (Spike), R203K, G204R(N) are the most prominent. Most importantly, numerous mutations were flourishing in the N protein gene (67) followed by S (45), RdRp (38), NSP2 (34), NSP3 (20), and ORF8 (6) gene. Moreover, nucleotide deletion analysis found nine deletions throughout the genomes including in ORF7a (8), ORF8 (1) with one insertion (G) at 265 positions in only one genome. The underlying mechanism of disease severity, molecular evolution, and epidemiology lie in genomic sequences that are not fully understood yet. Identification of the evolutionary history, parsimony-informative sites and others genetic variations of this deadly virus will facilitate the development of new strategies to control the local transmission and provide deep insight in the identification of potential therapeutic targets for controlling COVID-19.