Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), a novel evolutionary divergent RNA virus, is responsible for the present devastating COVID-19 pandemic. To explore the genomic ...signatures, we comprehensively analyzed 2,492 complete and/or near-complete genome sequences of SARS-CoV-2 strains reported from across the globe to the GISAID database up to 30 March 2020. Genome-wide annotations revealed 1,516 nucleotide-level variations at different positions throughout the entire genome of SARS-CoV-2. Moreover, nucleotide (nt) deletion analysis found twelve deletion sites throughout the genome other than previously reported deletions at coding sequence of the ORF8 (open reading frame), spike, and ORF7a proteins, specifically in polyprotein ORF1ab (n = 9), ORF10 (n = 1), and 3´-UTR (n = 2). Evidence from the systematic gene-level mutational and protein profile analyses revealed a large number of amino acid (aa) substitutions (n = 744), demonstrating the viral proteins heterogeneous. Notably, residues of receptor-binding domain (RBD) showing crucial interactions with angiotensin-converting enzyme 2 (ACE2) and cross-reacting neutralizing antibody were found to be conserved among the analyzed virus strains, except for replacement of lysine with arginine at 378th position of the cryptic epitope of a Shanghai isolate, hCoV-19/Shanghai/SH0007/2020 (EPI_ISL_416320). Furthermore, our results of the preliminary epidemiological data on SARS-CoV-2 infections revealed that frequency of aa mutations were relatively higher in the SARS-CoV-2 genome sequences of Europe (43.07%) followed by Asia (38.09%), and North America (29.64%) while case fatality rates remained higher in the European temperate countries, such as Italy, Spain, Netherlands, France, England and Belgium. Thus, the present method of genome annotation employed at this early pandemic stage could be a promising tool for monitoring and tracking the continuously evolving pandemic situation, the associated genetic variants, and their implications for the development of effective control and prophylaxis strategies.
The emerged novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) has created a global health crisis that warrants an accurate and detailed characterization of the rapidly ...evolving viral genome for understanding its epidemiology, pathogenesis, and containment. Here, we explored 61,485 sequences of the nucleocapsid (N) protein, a potent diagnostic and prophylactic target, for identifying the mutations to review their roles in real‐time polymerase chain reaction based diagnosis and observe consequent impacts. Compared to the Wuhan reference strain, a total of 1034 unique nucleotide mutations were identified in the mutant strains (49.15%, n = 30,221) globally. Of these mutations, 367 occupy primer binding sites including the 3′‐end mismatch to the primer‐pair of 11 well‐characterized primer sets. Noteworthily, CDC (USA) recommended the N2 primer set contained a lower mismatch than the other primer sets. Moreover, 684 amino acid (aa) substitutions were located across 317 (75.66% of total aa) unique positions including 82, 21, and 83 of those in the RNA binding N‐terminal domain (NTD), SR‐rich region, and C‐terminal dimerization domain, respectively. Moreover, 11 in‐frame deletions, mostly (n = 10) within the highly flexible linker region, were revealed, and the rest was within the NTD region. Furthermore, we predicted the possible consequence of high‐frequency mutations (≥20) and deletions on the tertiary structure of the N protein. Remarkably, we observed that a high frequency (67.94% of mutated sequences) co‐occuring mutations (R203K and G204R) destabilized and decreased overall structural flexibility. The N protein of SARS‐CoV‐2 comprises an average of 1.2 mutations per strain compared to 4.4 and 0.4 in Middle East respiratory syndrome‐related coronavirus and SARS‐CoV, respectively. Despite being proposed as the alternative target to spike protein for vaccine and therapeutics, the ongoing evolution of the N protein may challenge these endeavors, thus needing further immunoinformatics analyses. Therefore, continuous monitoring is required for tracing the ongoing evolution of the SARS‐CoV‐2 N protein in prophylactic and diagnostic interventions.
Highlights
Analyzing 61,485 sequences of the nucleocapsid (N) protein, 1034 unique nucleotide mutations were identified in the mutant strains (49.15%, n=30,221) globally.
Of these mutations, 367 occupy primer binding sites including the 3'‐end mismatch to the primer‐pair of 11 well‐characterized primer sets.
In total, 684 amino acid (aa) substitutions in 317 unique positions including 82, 21, and 83 in RNA binding N‐terminal domain (NTD), SR‐rich region, and C‐terminal dimerization domain (CTD), respectively were detected.
Eleven in‐frame deletions found, mostly (n = 10) within the highly flexible linker region of the SARS‐CoV‐2 N protein.
The N protein of SARS‐CoV‐2 comprises an average of 1.2 mutations per strain compared to 4.4 and 0.4 in MERS‐CoV and SARS‐CoV, respectively.
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiologic agent of the ongoing pandemic of coronavirus disease 2019 (COVID-19), a public health emergency of international concerns ...declared by the World Health Organization (WHO). An immuno-informatics approach along with comparative genomics was applied to design a multi-epitope-based peptide vaccine against SARS-CoV-2 combining the antigenic epitopes of the S, M, and E proteins. The tertiary structure was predicted, refined and validated using advanced bioinformatics tools. The candidate vaccine showed an average of ≥90.0% world population coverage for different ethnic groups. Molecular docking and dynamics simulation of the chimeric vaccine with the immune receptors (TLR3 and TLR4) predicted efficient binding. Immune simulation predicted significant primary immune response with increased IgM and secondary immune response with high levels of both IgG1 and IgG2. It also increased the proliferation of T-helper cells and cytotoxic T-cells along with the increased IFN-γ and IL-2 cytokines. The codon optimization and mRNA secondary structure prediction revealed that the chimera is suitable for high-level expression and cloning. Overall, the constructed recombinant chimeric vaccine candidate demonstrated significant potential and can be considered for clinical validation to fight against this global threat, COVID-19.
The microbiome of the anaerobic digester (AD) regulates the level of energy production. To assess the microbiome diversity and composition in different stages of anaerobic digestion, we collected 16 ...samples from the AD of cow dung (CD) origin. The samples were categorized into four groups (Group-I, Group-II, Group-III and Group-IV) based on the level of energy production (CH
%), and sequenced through whole metagenome sequencing (WMS). Group-I (n = 2) belonged to initial time of energy production whereas Group-II (n = 5), Group-III (n = 5), and Group-IV (n = 4) had 21-34%, 47-58% and 71-74% of CH
, respectively. The physicochemical analysis revealed that level of energy production (CH
%) had significant positive correlation with digester pH (r = 0.92, p < 0.001), O
level (%) (r = 0.54, p < 0.05), and environmental temperature (°C) (r = 0.57, p < 0.05). The WMS data mapped to 2800 distinct bacterial, archaeal and viral genomes through PathoScope (PS) and MG-RAST (MR) analyses. We detected 768, 1421, 1819 and 1774 bacterial strains in Group-I, Group-II, Group-III and Group-IV, respectively through PS analysis which were represented by Firmicutes, Bacteroidetes, Proteobacteria, Actinobacteria, Spirochaetes and Fibrobacteres phyla (> 93.0% of the total abundances). Simultaneously, 343 archaeal strains were detected, of which 95.90% strains shared across four metagenomes. We identified 43 dominant species including 31 bacterial and 12 archaeal species in AD microbiomes, of which only archaea showed positive correlation with digester pH, CH
concentration, pressure and temperature (Spearman correlation; r > 0.6, p < 0.01). The indicator species analysis showed that the species Methanosarcina vacuolate, Dehalococcoides mccartyi, Methanosarcina sp. Kolksee and Methanosarcina barkeri were highly specific for energy production. The correlation network analysis showed that different strains of Euryarcheota and Firmicutes phyla exhibited significant correlation (p = 0.021, Kruskal-Wallis test; with a cutoff of 1.0) with the highest level (74.1%) of energy production (Group-IV). In addition, top CH
producing microbiomes showed increased genomic functional activities related to one carbon and biotin metabolism, oxidative stress, proteolytic pathways, membrane-type-1-matrix-metalloproteinase (MT1-MMP) pericellular network, acetyl-CoA production, motility and chemotaxis. Importantly, the physicochemical properties of the AD including pH, CH
concentration (%), pressure, temperature and environmental temperature were found to be positively correlated with these genomic functional potentials and distribution of ARGs and metal resistance pathways (Spearman correlation; r > 0.5, p < 0.01). This study reveals distinct changes in composition and diversity of the AD microbiomes including different indicator species, and their genomic features that are highly specific for energy production.
Milk microbiomes significantly influence the pathophysiology of bovine mastitis. To assess the association between microbiome diversity and bovine mastitis, we compared the microbiome of clinical ...mastitis (CM, n = 14) and healthy (H, n = 7) milk samples through deep whole metagenome sequencing (WMS). A total of 483.38 million reads generated from both metagenomes were analyzed through PathoScope (PS) and MG-RAST (MR), and mapped to 380 bacterial, 56 archaeal, and 39 viral genomes. We observed distinct shifts and differences in abundance between the microbiome of CM and H milk in phyla Proteobacteria, Bacteroidetes, Firmicutes and Actinobacteria with an inclusion of 68.04% previously unreported and/or opportunistic strains in CM milk. PS identified 363 and 146 bacterial strains in CM and H milk samples respectively, and MR detected 356 and 251 bacterial genera respectively. Of the identified taxa, 29.51% of strains and 63.80% of genera were shared between both metagenomes. Additionally, 14 archaeal and 14 viral genera were found to be solely associated with CM. Functional annotation of metagenomic sequences identified several metabolic pathways related to bacterial colonization, proliferation, chemotaxis and invasion, immune-diseases, oxidative stress, regulation and cell signaling, phage and prophases, antibiotic and heavy metal resistance that might be associated with CM. Our WMS study provides conclusive data on milk microbiome diversity associated with bovine CM and its role in udder health.
Wearing masks or personal protective equipment (PPE) has become an integral part of the occupational life of physicians due to the coronavirus disease 2019 (COVID-19) pandemic. Most physicians have ...been developing various health hazards related to the use of different protective gears. This study aimed to determine the burden and spectrum of various health hazards associated with using masks or PPE and their associated risk factors. This cross-sectional survey was conducted in Dhaka Medical College from March 01-May 30, 2021, among physicians from different public hospitals in Dhaka, Bangladesh. We analyzed the responses of 506 physicians who completed case record forms through Google forms or hard copies. The mean (SD) age of the respondents was 35.4 7.7, and 69.4% were men. Approximately 40% were using full PPE, and 55% were using N-95 masks. A total of 489 (96.6%) patients experienced at least one health hazard. The reported severe health hazards were syncope, severe dyspnea, severe chest pain, and anaphylaxis. Headache, dizziness, mood irritation, chest pain, excessive sweating, panic attack, and permanent facial disfigurement were the minor health hazards reported. Extended periods of work in the COVID-19-unit, reuse of masks, diabetes, obesity, and mental stress were risk factors for dyspnea. The risk factors for headaches were female sex, diabetes, and previous primary headaches. Furthermore, female sex and reusing masks for an extended period (> 6 h) were risk factors for facial disfigurement. The risk factors for excessive sweating were female sex and additional evening office practice for an extended period. Healthcare workers experienced several occupational hazards after using masks and PPE. Therefore, an appropriate policy is required to reduce such risks.
The field of fish microbiome research has rapidly been advancing, primarily focusing on farmed or laboratory fish species rather than natural or marine fish populations. This study sought to reveal ...the distinctive gut bacteriome composition and diversity within the anadromous fish species Tenualosa ilisha (hilsa), which holds the status of being the national fish of Bangladesh. We conducted an analysis on 15 gut samples obtained from 15 individual hilsa fishes collected from three primary habitats (e.g., freshwater = 5, brackish water = 5 and marine water = 5) in Bangladesh. The analysis utilized metagenomics based on 16S rRNA gene sequencing targeting the V3-V4 regions. Our comprehensive identification revealed a total of 258 operational taxonomic units (OTUs). The observed OTUs were represented by six phyla, nine classes, 19 orders, 26 families and 40 genera of bacteria. Our analysis unveiled considerable taxonomic differences among the habitats (freshwater, brackish water, and marine water) of hilsa fishes, as denoted by a higher level of shared microbiota (p = 0.007, Kruskal-Wallis test). Among the identified genera in the gut of hilsa fishes, including Vagococcus, Morganella, Enterobacter, Plesiomonas, Shigella, Clostridium, Klebsiella, Serratia, Aeromonas, Macrococcus, Staphylococcus, Proteus, and Hafnia, several are recognized as fish probiotics. Importantly, some bacterial genera such as Sinobaca, Synechococcus, Gemmata, Serinicoccus, Saccharopolyspora, and Paulinella identified in the gut of hilsa identified in this study have not been reported in any aquatic or marine fish species. Significantly, we observed that 67.50% (27/40) of bacterial genera were found to be common among hilsa fishes across all three habitats. Our findings offer compelling evidence for the presence of both exclusive and communal bacteriomes within the gut of hilsa fishes, exhibiting potential probiotic properties. These observations could be crucial for guiding future microbiome investigations in this economically significant fish species.
The emergence and spread of multidrug-resistant pathogens like Pseudomonas aeruginosa are major concerns for public health worldwide. This study aimed to assess the prevalence of P. aeruginosa in ...clinical, environmental, and poultry sources in Bangladesh, along with their antibiotic susceptibility and the profiling of β-lactamase and virulence genes using standard molecular and microbiology techniques. We collected 110 samples from five different locations, viz., BAU residential area (BAURA; n = 15), BAU Healthcare Center (BAUHCC; n = 20), BAU Veterinary Teaching Hospital (BAUVTH; n = 22), Poultry Market (PM; n = 30) and Mymensingh Medical College Hospital (MCCH; n = 23). After overnight enrichment in nutrient broth, 89 probable Pseudomonas isolates (80.90%) were screened through selective culture, gram-staining and biochemical tests. Using genus- and species-specific PCR, we confirmed 22 isolates (20.0%) as P. aeruginosa from these samples. Antibiogram profiling revealed that 100.0% P. aeruginosa isolates (n = 22) were multidrug-resistant isolates, showing resistance against Doripenem, Penicillin, Ceftazidime, Cefepime, and Imipenem. Furthermore, resistance to aztreonam was observed in 95.45% isolates. However, P. aeruginosa isolates showed a varying degree of sensitivity against Amikacin, Gentamicin, and Ciprofloxacin. The blaTEM gene was detected in 86.0% isolates, while blaCMY, blaSHV and blaOXA, were detected in 27.0%, 18.0% and 5.0% of the P. aeruginosa isolates, respectively. The algD gene was detected in 32.0% isolates, whereas lasB and exoA genes were identified in 9.0% and 5.0% P. aeruginosa isolates. However, none of the P. aeruginosa isolates harbored exoS gene. Hence, this study provides valuable and novel insights on the resistance and virulence of circulating P. aeruginosa within the clinical, environmental, and poultry environments of Bangladesh. These findings are crucial for understanding the emergence of β-lactamase resistance in P. aeruginosa, highlighting its usefulness in the treatment and control of P. aeruginosa infections in both human and animal populations.
We previously reported that SARS-CoV-2 infection reduces human nasopharyngeal commensal microbiomes (bacteria, archaea and commensal respiratory viruses) with inclusion of pathobionts. This study ...aimed to assess the possible changes in the abundance and diversity of resident mycobiome in the nasopharyngeal tract (NT) of humans due to SARS-CoV-2 infections. Twenty-two (n = 22) nasopharyngeal swab samples (including COVID-19 = 8, Recovered = 7, and Healthy = 7) were collected for RNA-sequencing followed by taxonomic profiling of mycobiome. Our analyses indicate that SARS-CoV-2 infection significantly increased (p < 0.05, Wilcoxon test) the population and diversity of fungi in the NT with inclusion of a high proportion of opportunistic pathogens. We detected 863 fungal species including 533, 445, and 188 species in COVID-19, Recovered, and Healthy individuals, respectively that indicate a distinct mycobiome dysbiosis due to the SARS-CoV-2 infection. Remarkably, 37% of the fungal species were exclusively associated with SARS-CoV-2 infection, where S. cerevisiae (88.62%) and Phaffia rhodozyma (10.30%) were two top abundant species. Likewise, Recovered humans NT samples were predominated by Aspergillus penicillioides (36.64%), A. keveii (23.36%), A. oryzae (10.05%) and A. pseudoglaucus (4.42%). Conversely, Nannochloropsis oceanica (47.93%), Saccharomyces pastorianus (34.42%), and S. cerevisiae (2.80%) were the top abundant fungal species in Healthy controls nasal swabs. Importantly, 16% commensal fungal species found in the Healthy controls were not detected in either COVID-19 patients or when they were cured from COVID-19 (Recovered). We also detected several altered metabolic pathways correlated with the dysbiosis of fungal mycobiota in COVID-19 patients. Our results suggest that SARS-CoV-2 infection causes significant dysbiosis of mycobiome and related metabolic functions possibly play a determining role in the progression of SARS-CoV-2 pathogenesis. These findings might be helpful for developing mycobiome-based diagnostics, and also devising appropriate therapeutic regimens including antifungal drugs for prevention and control of concurrent fungal coinfections in COVID-19 patients.
Abstract
The identification of deleterious mutations in different variants of SARS-CoV-2 and their roles in the morbidity of COVID-19 patients has yet to be thoroughly investigated. To unravel the ...spectrum of mutations and their effects within SARS-CoV-2 genomes, we analyzed 5,724 complete genomes from deceased COVID-19 patients sourced from the GISAID database. This analysis was conducted using the Nextstrain platform, applying a generalized time-reversible model for evolutionary phylogeny. These genomes were compared to the reference strain (hCoV-19/Wuhan/WIV04/2019) using MAFFT v7.470. Our findings revealed that SARS-CoV-2 genomes from deceased individuals belonged to 21 Nextstrain clades, with clade 20I (Alpha variant) being the most predominant, followed by clade 20H (Beta variant) and clade 20J (Gamma variant). The majority of SARS-CoV-2 genomes from deceased patients (33.4%) were sequenced in North America, while the lowest percentage (0.98%) came from Africa. The ‘G’ clade was dominant in the SARS-CoV-2 genomes of Asian, African, and North American regions, while the ‘GRY’ clade prevailed in Europe. In our analysis, we identified 35,799 nucleotide (NT) mutations throughout the genome, with the highest frequency (11,402 occurrences) found in the spike protein. Notably, we observed 4150 point-specific amino acid (AA) mutations in SARS-CoV-2 genomes, with D614G (20%) and N501Y (14%) identified as the top two deleterious mutations in the spike protein on a global scale. Furthermore, we detected five common deleterious AA mutations, including G18V, W45S, I33T, P30L, and Q418H, which play a key role in defining each clade of SARS-CoV-2. Our novel findings hold potential value for genomic surveillance, enabling the monitoring of the evolving pattern of SARS-CoV-2 infection, its emerging variants, and their impact on the development of effective vaccination and control strategies.