Dynamics of SARS-CoV-2 variants of concern (VOC) in Bangladesh during the first half of 2021

SARS-CoV-2 genomic surveillance and reliability of PCR single point mutation assay (SNPsig® SARS-CoV-2 EscapePLEX CE) for the rapid detection of variants of concern in Cameroon

Background

Surveillance of SARS-CoV-2 variants of concern (VOCs) and lineages is crucial for decision-making. Our objective was to study the SARS-CoV-2 clade dynamics across epidemiological waves and evaluate the reliability of SNPsig® SARS-CoV-2 EscapePLEX CE in detecting VOCs in Cameroon.

Material and methods

A laboratory-based study was conducted on SARS-CoV-2 positive nasopharyngeal specimens cycle threshold (Ct)≤30 at the Chantal BIYA International Reference Centre in Yaoundé-Cameroon, between April-2020 to August-2022. Samples were analyzed in parallel with Sanger sequencing and (SNPsig® SARS-CoV-2 EscapePLEX CE), and performance characteristics were evaluated by Cohen's coefficient and McNemar test.

Results

Of the 130 sequences generated, SARS-CoV-2 clades during wave-1 (April-November 2020) showed 97 % (30/31) wild-type lineages and 3 % (1/31) Gamma-variant; wave-2 (December-2020 to May-2021), 25 % (4/16) Alpha-variant, 25 % (4/16) Beta-variant, 44 % (7/16) wild-type and 6 % (1/16) mu; wave-3 (June-October 2021), 94 % (27/29) Delta-variant, 3 % (1/29) Alpha-variant, 3 % (1/29) wild-type; wave-4 (November-2021 to August-2022), 98 % (53/54) Omicron-variant and 2 % (1/54) Delta-variant. Omicron sub-variants were BA.1 (47 %), BA.5 (34 %), BA.2 (13 %) and BA.4 (6 %). Globally, the two genotyping methods accurately identified the SARS-CoV-2 VOCs (P = 0.17, McNemar test; Ka = 0.67).

Conclusion

Genomic surveillance reveals a rapid dynamic in SARS-CoV-2 strains between epidemiological waves in Cameroon. For wide-spread variant surveillance in resource-limited settings, SNPsig® SARS-CoV-2 EscapePLEX CEkit represents a suitable tool, pending upgrading for distinguishing Omicron sub-lineages.

Molecular insights into the SARS-CoV-2 Omicron variant from Bangladesh suggest diverse and continuous evolution

The study analyzed the molecular dynamics of the circulating SARS-CoV-2 Omicron variant from its identification in November 2021 to January 2023. The SARS-CoV-2 sequences from Bangladesh revealed three distinct waves of the Omicron variant. More than 50 sub-lineages of Omicron variant were introduced into the country, with the majority belonging to the major lineages of BA.1-like (24.91%), BA.2-like (43.35%), BA.5-like (5.76%), XBB (10.47%), and "Others and Unassigned" (18.64%). Furthermore, the relative frequencies over time revealed that Omicron lineages existed for a short period of time before being replaced by other sub-lineages. Many potential mutations were found in the receptor binding domain of the Spike protein including G339D/H, S371 L/F, K417 N, T478K, E484A, Q493R, Q498R, and N501Y. In conclusion, the SARS-CoV-2 Omicron variant from Bangladesh showed diverse genetic features and continuous evolution. Therefore, the choice of vaccine and monitoring of hospitalized patients is important alongside genetic characterization of the circulating SARS-CoV-2.

Temporal dynamics and fatality of SARS-CoV-2 variants in Bangladesh

Background and Aims

Since the beginning of the SARS-CoV-2 pandemic, multiple new variants have emerged posing an increased risk to global public health. This study aimed to investigate SARS-CoV-2 variants, their temporal dynamics, infection rate (IFR) and case fatality rate (CFR) in Bangladesh by analyzing the published genomes.

Methods

We retrieved 6610 complete whole genome sequences of the SARS-CoV-2 from the GISAID (Global Initiative on Sharing all Influenza Data) platform from March 2020 to October 2022, and performed different in-silico bioinformatics analyses. The clade and Pango lineages were assigned by using Nextclade v2.8.1. SARS-CoV-2 infections and fatality data were collected from the Institute of Epidemiology Disease Control and Research (IEDCR), Bangladesh. The average IFR was calculated from the monthly COVID-19 cases and population size while average CFR was calculated from the number of monthly deaths and number of confirmed COVID-19 cases.

Results

SARS-CoV-2 first emerged in Bangladesh on March 3, 2020 and created three pandemic waves so far. The phylogenetic analysis revealed multiple introductions of SARS-CoV-2 variant(s) into Bangladesh with at least 22 Nextstrain clades and 107 Pangolin lineages with respect to the SARS-CoV-2 reference genome of Wuhan/Hu-1/2019. The Delta variant was detected as the most predominant (48.06%) variant followed by Omicron (27.88%), Beta (7.65%), Alpha (1.56%), Eta (0.33%) and Gamma (0.03%) variant. The overall IFR and CFR from circulating variants were 13.59% and 1.45%, respectively. A time-dependent monthly analysis showed significant variations in the IFR (p = 0.012, Kruskal-Wallis test) and CFR (p = 0.032, Kruskal-Wallis test) throughout the study period. We found the highest IFR (14.35%) in 2020 while Delta (20A) and Beta (20H) variants were circulating in Bangladesh. Remarkably, the highest CFR (1.91%) from SARS-CoV-2 variants was recorded in 2021.

Conclusion

Our findings highlight the importance of genomic surveillance for careful monitoring of variants of concern emergence to interpret correctly their relative IFR and CFR, and thus, for implementation of strengthened public health and social measures to control the spread of the virus. Furthermore, the results of the present study may provide important context for sequence-based inference in SARS-CoV-2 variant(s) evolution and clinical epidemiology beyond Bangladesh.

A Statistical Synopsis of COVID-19 Components and Descriptive Analysis of Their Socio-Economic and Healthcare Aspects in Bangladesh Perspective

The aim of the work is to analyze the socio-economic and healthcare aspects that arise in the contemporary COVID-19 situation from Bangladesh perspective. We elaborately discuss the successive COVID-19 occurrences in Bangladesh with consequential information. The components associated with the COVID-19 commencement and treatment policy with corresponding features and their consequences are patently delineated. The effect of troublesome issues related to the treatment is detailed with supporting real-time data. We elucidate the applications of modern technologies advancement in epidemiological aspects and their existent compatibility in Bangladesh. We statistically analyze the real-time data through figurative and tabular approaches. Some relevant measures of central tendency and dispersion are utilized to explore the data structure and its observable specifications. For a clear manifestation, Z- scores of the COVID-19 components are analyzed through the Box-Whisker plot. We have discovered that the gathered data exhibit features that are unsatisfactory for the normal distribution, are highly positively skewed, and are predominated by the earliest occurrences. Infections and deaths were initially lower than the global average, but they drastically rose in the first quarter of 2021 and persisted for the remainder of the year. Substantial preventive results were produced by the region-wisetime-worthy moves. In the fourth quarter of 2021, the infections and deaths noticeably decreased, and the number of recoveries was highly significant. In the middle of 2022, a lethal rise in infections was observed in Bangladesh and that was quickly stabilized, and the pandemic ingredients were under control. According to our assessment, some concluding remarks are made at the end of this work.

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.

Characterization and Structural Prediction of Proteins in SARS-CoV-2 Bangladeshi Variant Through Bioinformatics

The renowned respiratory disease induced by the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) has become a global epidemic in just less than a year by the first half of 2020. The subsequent efficient human-to-human transmission of this virus eventually affected millions of people worldwide. The most devastating thing is that the infection rate is continuously uprising and resulting in significant mortality especially among the older age population and those with health co-morbidities. This enveloped, positive-sense RNA virus is chiefly responsible for the infection of the upper respiratory system. The virulence of the SARS-CoV-2 is mostly regulated by its proteins such as entry to the host cell through fusion mechanism, fusion of infected cells with neighboring uninfected cells to spread virus, inhibition of host gene expression, cellular differentiation, apoptosis, mitochondrial biogenesis, etc. But very little is known about the protein structures and functionalities. Therefore, the main purpose of this study is to learn more about these proteins through bioinformatics approaches. In this study, ORF10, ORF7b, ORF7a, ORF6, membrane glycoprotein, and envelope protein have been selected from a Bangladeshi Corona-virus strain G039392 and a number of bioinformatics tools (MEGA-X-V10.1.7, PONDR, ProtScale, ProtParam, SCRIBER, NetSurfP v2.0, IntFOLD, UCSF Chimera, and PyMol) and strategies were implemented for multiple sequence alignment and phylogeny analysis with 9 different variants, predicting hydropathicity, amino acid compositions, protein-binding propensity, protein disorders, and 2D and 3D protein modeling. Selected proteins were characterized as highly flexible, structurally and electrostatically extremely stable, ordered, biologically active, hydrophobic, and closely related to proteins of different variants. This detailed information regarding the characterization and structure of proteins of SARS-CoV-2 Bangladeshi variant was performed for the first time ever to unveil the deep mechanism behind the virulence features. And this robust appraisal also paves the future way for molecular docking, vaccine development targeting these characterized proteins.

Asymptomatic RT-PCR positive COVID-19 patients in orthopaedic pre-operative evaluation during the peak of the second wave

Introduction

Asymptomatic COVID-19 patients are the most challenging and feared obstacles in resuming these surgical procedures. The purpose of this study was to evaluate the proportion of asymptomatic carriers detected by RT-PCR in pre-operative orthopaedic evaluation during the peak of the second wave.

Methods

514 asymptomtomatic COVID-19 patients, negative for TOCC (Travel, Profession, Cluster, Contact) risk factors were observed retrospectively. A nasopharyFngeal RT-PCR test was obtained 48 to 72 h before the surgery in all cases. Possible risk factors for a positive test was identified.

Results

The detected asymptomatic COVID-19 infection rate during the peak of the second wave among the pre-operative orthopaedic patients was 12.3%. Younger age, female gender, longer duration of admission to RT-PCR test interval were found to be significant (p= < 0.05) risk factors for asymptomatic RT-PCR to be positive. The hazard ratio (HR) for being asymptomatic RT-PCR positive was 4.3 (p = 0. 025), while the RT-PCR was performed at 14 days, but the HR increased to 9.2 (p = 0.049) when the test was performed after 45 days.

Conclusion

According to our findings, pre-operative testing to rule out COVID-19 should be regarded as a critical step in preventing the disease clusters in hospitals.

Codon usage divergence in Delta variants (B.1.617.2) of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spreads all over the world and brings great harm to humans in many countries. Many new SARS-CoV-2 variants appeared during its transmission. In the present study, the Delta variants (B.1.617.2) of SARS-CoV-2, which have appeared in many countries, were considered for analysis. In order to evaluate the evolutionary divergence of the Delta variants(B.1.617.2), the codon usage divergence in Delta variants (B.1.617.2) of SARS-CoV-2 was compared to that of the SARS-CoV-2 genomes emerged before June 2020. All Delta variants (B.1.617.2) and 350 early genomes of SARS-CoV-2 in the NCBI database were downloaded. Codon usage pattern including the basic composition, the GC ratio of the third position (GC3) and the first two positions (GC12) in codons, overall GC contents, the effective number of codons (ENC), the codon bias index (CBI), the relative synonymous codon usage (RSCU) values, etc., of all concerned important gene sequences were all calculated. Codon usage divergence of them was calculated via summing their standard deviations. The results suggested that base compositions in both Delta variants (B.1.617.2) of SARS-CoV-2 and the early SARS-CoV-2 genomes were similar to each other. However, the internal codon usage divergence for most genes in Delta variants (B.1.617.2) was significantly wider than that of SARS-CoV-2. The RSCU values were further used to explore the synonymous and non-synonymous mutations in the sequences of the Delta variants (B.1.617.2), and the results showed the synonymous mutations are more obvious than the non-synonymous in the concerned sequences. The related codon usage divergence analysis is helpful for further study on the adaptability and disease prognosis of the SARS-CoV-2 variants.