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

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.

Unravelling the link between SARS-CoV-2 mutation frequencies, patient comorbidities, and structural dynamics

Genomic surveillance is crucial for tracking emergence and spread of novel variants of pathogens, such as SARS-CoV-2, to inform public health interventions and to enforce control measures. However, in some settings especially in low- and middle- income counties, where sequencing platforms are limited, only certain patients get to be selected for sequencing surveillance. Here, we show that patients with multiple comorbidities potentially harbour SARS-CoV-2 with higher mutation rates and thus deserve more attention for genomic surveillance. The relationship between the patient comorbidities, and type of amino acid mutations was assessed. Correlation analysis showed that there was a significant tendency for mutations to occur within the ORF1a region for patients with higher number of comorbidities. Frequency analysis of the amino acid substitution within ORF1a showed that nsp3 P822L of the PLpro protease was one of the highest occurring mutations. Using molecular dynamics, we simulated that the P822L mutation in PLpro represents a system with lower Root Mean Square Deviation (RMSD) fluctuations, and consistent Radius of gyration (Rg), Solvent Accessible Surface Area (SASA) values-indicate a much stabler protein than the wildtype. The outcome of this study will help determine the relationship between the clinical status of a patient and the mutations of the infecting SARS-CoV-2 virus.

Global landscape of SARS-CoV-2 mutations and conserved regions

BACKGROUND

At the end of December 2019, a novel strain of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) disease (COVID-19) has been identified in Wuhan, a central city in China, and then spread to every corner of the globe. As of October 8, 2022, the total number of COVID-19 cases had reached over 621 million worldwide, with more than 6.56 million confirmed deaths. Since SARS-CoV-2 genome sequences change due to mutation and recombination, it is pivotal to surveil emerging variants and monitor changes for improving pandemic management.

METHODS

10,287,271 SARS-CoV-2 genome sequence samples were downloaded in FASTA format from the GISAID databases from February 24, 2020, to April 2022. Python programming language (version 3.8.0) software was utilized to process FASTA files to identify variants and sequence conservation. The NCBI RefSeq SARS-CoV-2 genome (accession no. NC_045512.2) was considered as the reference sequence.

RESULTS

Six mutations had more than 50% frequency in global SARS-CoV-2. These mutations include the P323L (99.3%) in NSP12, D614G (97.6) in S, the T492I (70.4) in NSP4, R203M (62.8%) in N, T60A (61.4%) in Orf9b, and P1228L (50.0%) in NSP3. In the SARS-CoV-2 genome, no mutation was observed in more than 90% of nsp11, nsp7, nsp10, nsp9, nsp8, and nsp16 regions. On the other hand, N, nsp3, S, nsp4, nsp12, and M had the maximum rate of mutations. In the S protein, the highest mutation frequency was observed in aa 508-635(0.77%) and aa 381-508 (0.43%). The highest frequency of mutation was observed in aa 66-88 (2.19%), aa 7-14, and aa 164-246 (2.92%) in M, E, and N proteins, respectively.

CONCLUSION

Therefore, monitoring SARS-CoV-2 proteomic changes and detecting hot spots mutations and conserved regions could be applied to improve the SARS-CoV-2 diagnostic efficiency and design safe and effective vaccines against emerging variants.

Mutational analysis of the spike protein of SARS-COV-2 isolates revealed atomistic features responsible for higher binding and infectivity

Introduction: The perpetual appearance of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-COV-2), and its new variants devastated the public health and social fabric around the world. Understanding the genomic patterns and connecting them to phenotypic attributes is of great interest to devise a treatment strategy to control this pandemic. Materials and Methods: In this regard, computational methods to understand the evolution, dynamics and mutational spectrum of SARS-CoV-2 and its new variants are significantly important. Thus, herein, we used computational methods to screen the genomes of SARS-CoV-2 isolated from Pakistan and connect them to the phenotypic attributes of spike protein; we used stability-function correlation methods, protein-protein docking, and molecular dynamics simulation. Results: Using the Global initiative on sharing all influenza data (GISAID) a total of 21 unique mutations were identified, among which five were reported as stabilizing while 16 were destabilizing revealed through mCSM, DynaMut 2.0, and I-Mutant servers. Protein-protein docking with Angiotensin-converting enzyme 2 (ACE2) and monoclonal antibody (4A8) revealed that mutation G446V in the receptor-binding domain; R102S and G181V in the N-terminal domain (NTD) significantly affected the binding and thus increased the infectivity. The interaction pattern also revealed significant variations in the hydrogen bonding, salt bridges and non-bonded contact networks. The structural-dynamic features of these mutations revealed the global dynamic trend and the finding energy calculation further established that the G446V mutation increases the binding affinity towards ACE2 while R102S and G181V help in evading the host immune response. The other mutations reported supplement these processes indirectly. The binding free energy results revealed that wild type-RBD has a TBE of -60.55 kcal/mol while G446V-RBD reported a TBE of -73.49 kcal/mol. On the other hand, wild type-NTD reported -67.77 kcal/mol of TBE, R102S-NTD reported -51.25 kcal/mol of TBE while G181V-NTD reported a TBE of -63.68 kcal/mol. Conclusions: In conclusion, the current findings revealed basis for higher infectivity and immune evasion associated with the aforementioned mutations and structure-based drug discovery against such variants.

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.

Evolutionary Dynamics and Epidemiology of Endemic and Emerging Coronaviruses in Humans, Domestic Animals, and Wildlife

Diverse coronavirus (CoV) strains can infect both humans and animals and produce various diseases. CoVs have caused three epidemics and pandemics in the last two decades, and caused a severe impact on public health and the global economy. Therefore, it is of utmost importance to understand the emergence and evolution of endemic and emerging CoV diversity in humans and animals. For diverse bird species, the Infectious Bronchitis Virus is a significant one, whereas feline enteric and canine coronavirus, recombined to produce feline infectious peritonitis virus, infects wild cats. Bovine and canine CoVs have ancestral relationships, while porcine CoVs, especially SADS-CoV, can cross species barriers. Bats are considered as the natural host of diverse strains of alpha and beta coronaviruses. Though MERS-CoV is significant for both camels and humans, humans are nonetheless affected more severely. MERS-CoV cases have been reported mainly in the Arabic peninsula since 2012. To date, seven CoV strains have infected humans, all descended from animals. The severe acute respiratory syndrome coronaviruses (SARS-CoV and SARS-CoV-2) are presumed to be originated in Rhinolopoid bats that severely infect humans with spillover to multiple domestic and wild animals. Emerging alpha and delta variants of SARS-CoV-2 were detected in pets and wild animals. Still, the intermediate hosts and all susceptible animal species remain unknown. SARS-CoV-2 might not be the last CoV to cross the species barrier. Hence, we recommend developing a universal CoV vaccine for humans so that any future outbreak can be prevented effectively. Furthermore, a One Health approach coronavirus surveillance should be implemented at human-animal interfaces to detect novel coronaviruses before emerging to humans and to prevent future epidemics and pandemics.

SARS-CoV-2 variant evolution in the United States: High accumulation of viral mutations over time likely through serial Founder Events and mutational bursts

Since the first case of COVID-19 in December 2019 in Wuhan, China, SARS-CoV-2 has spread worldwide and within a year and a half has caused 3.56 million deaths globally. With dramatically increasing infection numbers, and the arrival of new variants with increased infectivity, tracking the evolution of its genome is crucial for effectively controlling the pandemic and informing vaccine platform development. Our study explores evolution of SARS-CoV-2 in a representative cohort of sequences covering the entire genome in the United States, through all of 2020 and early 2021. Strikingly, we detected many accumulating Single Nucleotide Variations (SNVs) encoding amino acid changes in the SARS-CoV-2 genome, with a pattern indicative of RNA editing enzymes as major mutators of SARS-CoV-2 genomes. We report three major variants through October of 2020. These revealed 14 key mutations that were found in various combinations among 14 distinct predominant signatures. These signatures likely represent evolutionary lineages of SARS-CoV-2 in the U.S. and reveal clues to its evolution such as a mutational burst in the summer of 2020 likely leading to a homegrown new variant, and a trend towards higher mutational load among viral isolates, but with occasional mutation loss. The last quartile of 2020 revealed a concerning accumulation of mostly novel low frequency replacement mutations in the Spike protein, and a hypermutable glutamine residue near the putative furin cleavage site. Finally, end of the year data and 2021 revealed the gradual increase to prevalence of known variants of concern, particularly B.1.1.7, that have acquired additional Spike mutations. Overall, our results suggest that predominant viral genomes are dynamically evolving over time, with periods of mutational bursts and unabated mutation accumulation. This high level of existing variation, even at low frequencies and especially in the Spike-encoding region may become problematic when super-spreader events, akin to serial Founder Events in evolution, drive these rare mutations to prominence.