Coding-complete genomes of 18 SARS-CoV-2 Omicron JN.1, JN.1.4, and JN.1.11 sub-lineages in Bangladesh.

We report 18 coding-complete genome sequences of emerging SARS-CoV-2 Omicron sub-lineages JN.1, JN.1.4, and JN.1.11 from Bangladesh. Nasopharyngeal swab samples were obtained from individuals with COVID-19 symptoms between December 2023 and January 2024. Whole genome sequencing was performed following the ARTIC Network-based protocol using Oxford Nanopore Technology.

Coding-complete genomes of XBB.1.16, XBB.2.3, FL.4 (alias of XBB.1.9.1.4), and XBB.3 of SARS-CoV-2 Omicron isolated from Bangladesh.

We announce the coding-complete genomes of four different strains of SARS-CoV-2 Omicron lineages, XBB.1.16, XBB.2.3, FL.4 (alias of XBB.1.9.1.4), and XBB.3. These strains were obtained between October 2022 and May 2023 from nasopharyngeal swabs of four Bangladeshi individuals, while one of them had a travel history. Genomic data were produced by implementing ARTIC Network-based amplicon sequencing using the Oxford Nanopore Technology.

Natural selection shapes the evolution of SARS-CoV-2 Omicron in Bangladesh.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved to give rise to a highly transmissive and immune-escaping variant of concern, known as Omicron. Many aspects of the evolution of SARS-CoV-2 and the driving forces behind the ongoing Omicron outbreaks remain unclear. Substitution at the receptor-binding domain (RBD) in the spike protein is one of the primary strategies of SARS-CoV-2 Omicron to hinder recognition by the host angiotensin-converting enzyme 2 (ACE2) receptor and avoid antibody-dependent defense activation. Here, we scanned for adaptive evolution within the SARS-CoV-2 Omicron genomes reported from Bangladesh in the public database GISAID (www.gisaid.org; dated 2 April 2023). The ratio of the non-synonymous (Ka) to synonymous (Ks) nucleotide substitution rate, denoted as ω, is an indicator of the selection pressure acting on protein-coding genes. A higher proportion of non-synonymous to synonymous substitutions (Ka/Ks or ω > 1) indicates positive selection, while Ka/Ks or ω near zero indicates purifying selection. An equal amount of non-synonymous and synonymous substitutions (Ka/Ks or ω = 1) refers to neutrally evolving sites. We found evidence of adaptive evolution within the spike (S) gene of SARS-CoV-2 Omicron isolated from Bangladesh. In total, 22 codon sites of the S gene displayed a signature of positive selection. The data also highlighted that the receptor-binding motif within the RBD of the spike glycoprotein is a hotspot of adaptive evolution, where many of the codons had ω > 1. Some of these adaptive sites at the RBD of the spike protein are known to be associated with increased viral fitness. The M gene and ORF6 have also experienced positive selection. These results suggest that although purifying selection is the dominant evolutionary force, positive Darwinian selection also plays a vital role in shaping the evolution of SARS-CoV-2 Omicron in Bangladesh.

Genetic diversity of Salmonella Paratyphi A isolated from enteric fever patients in Bangladesh from 2008 to 2018.

BACKGROUND: The proportion of enteric fever cases caused by Salmonella Paratyphi A is increasing and may increase further as we begin to introduce typhoid conjugate vaccines (TCVs). While numerous epidemiological and genomic studies have been conducted for S. Typhi, there are limited data describing the genomic epidemiology of S. Paratyphi A in especially in endemic settings, such as Bangladesh. PRINCIPAL FINDINGS: We conducted whole genome sequencing (WGS) of 67 S. Paratyphi A isolated between 2008 and 2018 from eight enteric disease surveillance sites across Bangladesh. We performed a detailed phylogenetic analysis of these sequence data incorporating sequences from 242 previously sequenced S. Paratyphi A isolates from a global collection and provided evidence of lineage migration from neighboring countries in South Asia. The data revealed that the majority of the Bangladeshi S. Paratyphi A isolates belonged to the dominant global lineage A (67.2%), while the remainder were either lineage C (19.4%) or F (13.4%). The population structure was relatively homogenous across the country as we did not find any significant lineage distributions between study sites inside or outside Dhaka. Our genomic data showed presence of single point mutations in gyrA gene either at codon 83 or 87 associated with decreased fluoroquinolone susceptibility in all Bangladeshi S. Paratyphi A isolates. Notably, we identified the pHCM2- like cryptic plasmid which was highly similar to S. Typhi plasmids circulating in Bangladesh and has not been previously identified in S. Paratyphi A organisms. SIGNIFICANCE: This study demonstrates the utility of WGS to monitor the ongoing evolution of this emerging enteric pathogen. Novel insights into the genetic structure of S. Paratyphi A will aid the understanding of both regional and global circulation patterns of this emerging pathogen and provide a framework for future genomic surveillance studies.

Genome-wide analysis provides a deeper understanding of the population structure of the Salmonella enterica serotype Paratyphi B complex in Bangladesh.

The Salmonella enterica serotype Paratyphi B complex causes a wide range of diseases, from gastroenteritis to paratyphoid fever, depending on the biotypes Java and sensu stricto. The burden of Paratyphi B biotypes in Bangladesh is still unknown, as these are indistinguishable by Salmonella serotyping. Here, we conducted the first whole-genome sequencing (WGS) study on 79 Salmonella isolates serotyped as Paratyphi B that were collected from 10 nationwide enteric disease surveillance sites in Bangladesh. Placing these in a global genetic context revealed that these are biotype Java, and the addition of these genomes expanded the previously described PG4 clade containing Bangladeshi and UK isolates. Importantly, antimicrobial resistance (AMR) genes were scarce amongst Bangladeshi S. Java isolates, somewhat surprisingly given the widespread availability of antibiotics without prescription. This genomic information provides important insights into the significance of S. Paratyphi B biotypes in enteric disease and their implications for public health.

Genotypic and phenotypic profiles of antibiotic-resistant bacteria isolated from hospitalised patients in Bangladesh.

OBJECTIVES: Characterisation of resistance phenotype and genotype is crucial to understanding the burden and transmission of antimicrobial resistance (AMR). This study aims to determine the spectrum of AMR and associated genes encoding aminoglycoside, macrolide and ß-lactam classes of antimicrobials in bacteria isolated from hospitalised patients in Bangladesh. METHODS: 430 bacterial isolates from patients with respiratory, intestinal, wound infections and typhoid fever, presenting to clinical care from 2015 to 2019, were examined. They included Escherichia coli (n = 85); Staphylococcus aureus (n = 84); Salmonella typhi (n = 82); Klebsiella pneumoniae (n = 42); Streptococcus pneumoniae (n = 36); coagulase-negative staphylococci (n = 28); Enterococcus faecalis (n = 27); Pseudomonas aeruginosa (n = 26); and Acinetobacter baumannii (n = 20). Reconfirmation of these clinical isolates and antimicrobial susceptibility tests was performed. PCR amplification using resistance gene-specific primers was done, and the amplified products were confirmed by Sanger sequencing. RESULTS: 53% of isolates were multidrug-resistant (MDR), including 97% of Escherichia coli. There was a year-wise gradual increase in MDR isolates from 2015 to 2018, and there was an almost twofold increase in the number of MDR strains isolated in 2019 (P = 0.00058). Among the 5 extended-spectrum ß-lactamases investigated, CTX-M-1 was the most prevalent (63%) followed by NDM-1 (22%); Escherichia coli was the major reservoir of these genes. The ermB (55%) and aac(6')-Ib (35%) genes were the most frequently detected macrolide and aminoglycoside resistance genes, respectively. CONCLUSION: MDR pathogens are highly prevalent in hospital settings of Bangladesh.

Population structure and antimicrobial resistance patterns of Salmonella Typhi isolates in urban Dhaka, Bangladesh from 2004 to 2016.

BACKGROUND: Multi-drug resistant typhoid fever remains an enormous public health threat in low and middle-income countries. However, we still lack a detailed understanding of the epidemiology and genomics of S. Typhi in many regions. Here we have undertaken a detailed genomic analysis of typhoid in urban Dhaka, Bangladesh to unravel the population structure and antimicrobial resistance patterns in S. Typhi isolated between 2004-2016. PRINCIPAL FINDINGS: Whole genome sequencing of 202 S. Typhi isolates obtained from three study locations in urban Dhaka revealed a diverse range of S. Typhi genotypes and AMR profiles. The bacterial population within Dhaka were relatively homogenous with little stratification between different healthcare facilities or age groups. We also observed evidence of exchange of Bangladeshi genotypes with neighboring South Asian countries (India, Pakistan and Nepal) suggesting these are circulating throughout the region. This analysis revealed a decline in H58 (genotype 4.3.1) isolates from 2011 onwards, coinciding with a rise in a diverse range of non-H58 genotypes and a simultaneous rise in isolates with reduced susceptibility to fluoroquinolones, potentially reflecting a change in treatment practices. We identified a novel S. Typhi genotype, subclade 3.3.2 (previously defined only to clade level, 3.3), which formed two localized clusters (3.3.2.Bd1 and 3.3.2.Bd2) associated with different mutations in the Quinolone Resistance Determining Region (QRDR) of gene gyrA. SIGNIFICANCE: Our analysis of S. Typhi isolates from urban Dhaka, Bangladesh isolated over a twelve year period identified a diverse range of AMR profiles and genotypes. The observed increase in non-H58 genotypes associated with reduced fluoroquinolone susceptibility may reflect a change in treatment practice in this region and highlights the importance of continued molecular surveillance to monitor the ongoing evolution of AMR in Dhaka. We have defined new genotypes and lineages of Bangladeshi S. Typhi which will facilitate the identification of these emerging AMR clones in future surveillance efforts.

Mutation Spectrum in TPO Gene of Bangladeshi Patients with Thyroid Dyshormonogenesis and Analysis of the Effects of Different Mutations on the Structural Features and Functions of TPO Protein through In Silico Approach.

Although thyroid dyshormonogenesis (TDH) accounts for 10-20% of congenital hypothyroidism (CH), the molecular etiology of TDH is unknown in Bangladesh. Thyroid peroxidase (TPO) is most frequently associated with TDH and the present study investigated the spectrum of TPO mutations in Bangladeshi patients and analyzed the effects of mutations on TPO protein structure through in silico approach. Sequencing-based analysis of TPO gene revealed four mutations in 36 diagnosed patients with TDH including three nonsynonymous mutations, namely, p.Ala373Ser, p.Ser398Thr, and p.Thr725Pro, and one synonymous mutation p.Pro715Pro. Homology modelling-based analysis of predicted structures of MPO-like domain (TPO142-738) and the full-length TPO protein (TPO1-933) revealed differences between mutant and wild type structures. Molecular docking studies were performed between predicted structures and heme. TPO1-933 predicted structure showed more reliable results in terms of interactions with the heme prosthetic group as the binding energies were -11.5 kcal/mol, -3.2 kcal/mol, -11.5 kcal/mol, and -7.9 kcal/mol for WT, p.Ala373Ser, p.Ser398Thr, and p.Thr725Pro, respectively, implying that p.Ala373Ser and p.Thr725Pro mutations were more damaging than p.Ser398Thr. However, for the TPO142-738 predicted structures, the binding energies were -11.9 kcal/mol, -10.8 kcal/mol, -2.5 kcal/mol, and -5.3 kcal/mol for the wild type protein, mutant proteins with p.Ala373Ser, p.Ser398Thr, and p.Thr725Pro substitutions, respectively. However, when the interactions between the crucial residues including residues His239, Arg396, Glu399, and His494 of TPO protein and heme were taken into consideration using both TPO1-933 and TPO142-738 predicted structures, it appeared that p.Ala373Ser and p.Thr725Pro could affect the interactions more severely than the p.Ser398Thr. Validation of the molecular docking results was performed by computer simulation in terms of quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) simulation. In conclusion, the substitutions mutations, namely, p.Ala373Ser, p.Ser398Thr, and p.Thr725Pro, had been involved in Bangladeshi patients with TDH and molecular docking-based study revealed that these mutations had damaging effect on the TPO protein activity.

High resolution melting curve analysis targeting the HBB gene mutational hot-spot offers a reliable screening approach for all common as well as most of the rare beta-globin gene mutations in Bangladesh.

BACKGROUND: Bangladesh lies in the global thalassemia belt, which has a defined mutational hot-spot in the beta-globin gene. The high carrier frequencies of beta-thalassemia trait and hemoglobin E-trait in Bangladesh necessitate a reliable DNA-based carrier screening approach that could supplement the use of hematological and electrophoretic indices to overcome the barriers of carrier screening. With this view in mind, the study aimed to establish a high resolution melting (HRM) curve-based rapid and reliable mutation screening method targeting the mutational hot-spot of South Asian and Southeast Asian countries that encompasses exon-1 (c.1 - c.92), intron-1 (c.92 + 1 - c.92 + 130) and a portion of exon-2 (c.93 - c.217) of the HBB gene which harbors more than 95% of mutant alleles responsible for beta-thalassemia in Bangladesh. RESULTS: Our HRM approach could successfully differentiate ten beta-globin gene mutations, namely c.79G > A, c.92 + 5G > C, c.126_129delCTTT, c.27_28insG, c.46delT, c.47G > A, c.92G > C, c.92 + 130G > C, c.126delC and c.135delC in heterozygous states from the wild type alleles, implying the significance of the approach for carrier screening as the first three of these mutations account for ~85% of total mutant alleles in Bangladesh. Moreover, different combinations of compound heterozygous mutations were found to generate melt curves that were distinct from the wild type alleles and from one another. Based on the findings, sixteen reference samples were run in parallel to 41 unknown specimens to perform direct genotyping of the beta-thalassemia specimens using HRM. The HRM-based genotyping of the unknown specimens showed 100% consistency with the sequencing result. CONCLUSIONS: Targeting the mutational hot-spot, the HRM approach could be successfully applied for screening of beta-thalassemia carriers in Bangladesh as well as in other countries of South Asia and Southeast Asia. The approach could be a useful supplement of hematological and electrophortic indices in order to avoid false positive and false negative results.

Cloning and overexpression of the 3-hydroxyisobutyrate dehydrogenase gene from pseudomonas putida E23.

The structural gene for NAD+-dependent 3-hydroxyisobutyrate dehydrogenase (EC 1.1.1.31) from Pseudomonas putida E23 was cloned in Escherichia coli cells to obtain a large amount of the enzyme and its nucleotides were sequenced to study its structural relationship with other proteins. The gene encoded a polypeptide containing 295 amino acid residues and was in a cluster with the gene for methylmalonate semialdehyde dehydrogenase. Transformed E. coli cells overproduced 3-hydroxyisobutyrate dehydrogenase, and the recombinant enzyme was purified to homogeneity with a high yield. Lysine and asparagine residues, which are important in catalysis of the 3-hydroxyacid dehydrogenase family, are conserved in this enzyme.

Cloning and sequencing of the serine dehydrogenase gene from Agrobacterium tumefaciens.

The structural gene for NADP+-dependent serine dehydrogenase [EC 1.1.1.-] from Agrobacterium tumefaciens ICR 1600 was cloned into Escherichia coli cells and its complete DNA sequence was analyzed. The gene encodes a polypeptide containing 249 amino acid residues. The enzyme had high sequence similarity to short-chain alcohol dehydrogenases from bacteria and unknown proteins of Haemophilus influenzae, Escherichia coli, and Saccharomyces cerevisiae.