VFDB 2019: a comparative pathogenomic platform with an interactive web interface

From sequence to significance: A thorough investigation of the distinctive genome features uncovered in C. Werkmanii strain NIB003

Citrobacter werkmanii (C. werkmanii), an opportunistic urinary bacterium that causes diarrhea, is poorly understood. Our research focuses on genetic features that are crucial to disease development, such as pathogenic interactions, antibiotic resistance, virulence genes and genetic variation. Following its morphological, biochemical, and molecular identification, the whole genome of C. werkmanii strain NIB003 was sequenced in Bangladesh for the first time. Despite having around 80% whole genome conservation, the research shows that the Bangladeshi strain forms a separate phylogenetic cluster. This emphasises the genetic variability within C. werkmanii, resulting in particular modifications at the strain level and changes in its ability to cause disease. The results of the genetic diversity analysis indicate that the Bangladeshi sequenced genome is more diverse than the other strains due to the existence of unique features, such as the presence of t-RNA binding domain and N-6 adenine-specific DNA methylases.

Comprehensive physiological and genomic characterization of a potential probiotic strain, Lactiplantibacillus plantarum ILSF15, isolated from the gut of tribes of Odisha, India

Characterizing probiotic features of organisms isolated from diverse environments can lead to the discovery of novel strains with promising functional features and health attributes. The present study attempts to characterize a novel probiotic strain isolated from the gut of the tribal population of Odisha, India. Based on 16S rRNA-based phylogeny, the strain was identified as a species of the Lactiplantibacillus genus and was named Lactiplantibacillus plantarum strain ILSF15. The current investigation focuses on elucidating this strain's genetic and physiological properties associated with probiotic attributes such as biosafety risk, host adaptation/survival traits, and beneficial functional features. The novel strain was observed, in vitro, exhibiting features such as acid/bile tolerance, adhesion to the host enteric epithelial cells, cholesterol assimilation, and pathogen exclusion, indicating its ability to survive the harsh environment of the human GIT and resist the growth of harmful microorganisms. Additionally, the L. plantarum ILSF15 strain was found to harbor genes associated with the metabolism and synthesis of various bioactive molecules, including amino acids, carbohydrates, lipids, and vitamins, highlighting the organism's ability to efficiently utilize diverse resources and contribute to the host's nutrition and health. Several genes involved in host adaptation/survival strategies and host-microbe interactions were also identified from the ILSF15 genome. Moreover, L. plantarum strains, in general, were found to have an open pangenome characterized by high genetic diversity and the absence of specific lineages associated with particular habitats, signifying its versatile nature and potential applications in probiotic and functional food industries.

psk1 virulence gene-induced pulmonary and systemic tuberculosis in a young woman with normal immune function: A case report

BACKGROUND

Tuberculosis is a chronic infectious disease and an important public health problem. Despite progress in controlling tuberculosis, the incidence of tuberculosis in China is still very high, with 895000 new cases annually. This case report describes the investigation of a case of severe disseminated tuberculosis in a young adult with normal immune function, conducted to ascertain why a Mycobacterium tuberculosis (M. tuberculosis) strain caused such severe disease.

CASE SUMMARY

A previously healthy 28-year-old woman presented to our hospital with a 1-month history of fever and fatigue. She was diagnosed with severe disseminated pulmonary tuberculosis, spinal tuberculosis with paravertebral abscesses, and tuberculous meningitis. M. tuberculosis was isolated from bronchoalveolar lavage fluid. She was treated with standard antituberculous therapy and underwent debridement, bone graft, and internal fixation surgery for spinal tuberculosis. She responded to therapy and regained her ability to walk following the surgery. We analysed the whole-genome sequence of the strain and designated it BLM-A21. Additional M. tuberculosis genomes were selected from the Virulence Factor Database (http://www.mgc.ac.cn/cgi-bin/VFs/genus.cgi?Genus=Mycobacterium) for comparison. An evolutionary tree of the BLM-A21 strain was built using PhyML maximum likelihood software. Further gene analysis revealed that, except for the pks1 gene, BLM-A21 had similar virulence genes to the CDC 1551 and H37Rv strains, which have lower dissemination.

CONCLUSION

We speculate that the pks1 virulence gene in BLM-A21 may be the key virulence gene responsible for the widespread dissemination of M. tuberculosis infection in this previously healthy adult with normal immune function.

Genetic characteristics and potential pathogenic agents in Campylobacter upsaliensis based on genomic analysis

Campylobacter upsaliensis was the most common Campylobacter species in pets' gastrointestinal tracts and has been isolated from patients with bacteremia, hemolytic-uremic syndrome, spontaneous abortion, and Guillain-Barré syndrome. However, the genetic characteristics and the full extent of its significance as a human pathogen remain to be fully understood. This study involved an investigation for genomic analysis of 154 strains from different sources and additional antimicrobial resistance profiles of 26 strains for this species. The genomes contained 1,558-1,971 CDS and the genome sizes were estimated to vary from 1.53 Mb to 1.86 Mb, with an average GC content of  34.71%. The entire analyzed genomes could be divided into three clades (A, B, and C) based on ANI and phylogenomic analysis. Significantly, nearly all strains in Clade B were isolated from patient samples, and the virulence-related sequences FlgD, GmhA, and CdtC might serve as determining factors for the classification of Clade B. Half of the tested isolates had MIC values over 64 μg mL-1 for nalidixic acid, gentamicin, and streptomycin. Isolates from pets in China carried more resistant elements in the genomes. This study both provided a comprehensive profile of C. upsaliensis for its genomic features and suggested some pathogenic agents for human infection with this species.

Dam construction alters the traits of health-related microbes along the Yangtze River

Dams, constructed globally for energy production and water conservation, fragment rivers, and modify flow regimes, thereby altering the composition of biological communities and ecosystem functions. Despite the extensive use of dams, few studies have explored their potential health impacts, particularly concerning changes in health-related genes, such as antibiotic resistance genes (ARGs) and virulence factor genes (VFGs), and their hosts (i.e., ARB and potential pathogens). Understanding these health-related effects is crucial because they can impact human health through water quality and pathogen prevalence. In this study, we investigated the planktonic microbial community in the Three Gorges Reservoir (TGR) and adjacent upstream and downstream areas of the Yangtze River during both the dry and wet season. Our metagenomic analysis showed that dam construction significantly decreased the abundance of ARGs, but it had an insignificant effect on VFGs. The observed reduction in ARGs abundance could be mainly attributed to the decrease abundance of the major ARGs carrier - Limnohabitansin the TGR and downstream areas due to high grazing pressure and fitness cost. Conversely, the abundance of microbes carrying VFGs (potential pathogens) remained stable from upstream to the dam reservoir, which may explain the negligible impact on VFG abundance. Overall, our results provide a detailed understanding of the ecological health implications of dam construction in large river ecosystems.

Epidemiology and antimicrobial resistance profiles of pathogenic Escherichia coli from commercial swine and poultry abattoirs and farms in South Africa: A One Health approach

Pathogenic Escherichia coli (PEC) are important foodborne bacteria that can cause severe illness in humans. The PECs thrive within the intestines of humans as well as animals and may contaminate multiple ecosystems, including food and water, via faecal transmission. Abattoir and farm employees are at high risk of PEC exposure, which could translate to community risk through person-to-person contact. To determine the epidemiology and resistome of PECs in Gauteng and Limpopo provinces of South Africa, 198 swine faecal samples, 220 poultry cloacal swabs, 108 human hand swabs, 11 run-off water samples from abattoirs and farms were collected from four swine and five poultry commercial abattoirs and two swine farms. One effluent sample each was collected from four wastewater treatment plants (WWTP) and a tertiary hospital setting. Phenotypic and genotypic techniques were used including polymerase chain reaction, pulsed-field gel electrophoresis (PFGE) and whole genome sequencing (WGS). Results showed EHEC and EPEC prevalence was 4.1 % (22/542) and 20.8 % (113/542), respectively, with the O26 serogroup detected the most in PEC isolates. According to the PFGE dendrogram, isolates from poultry, human hand swabs and run-off water clustered together. Diverse virulence factors such as the novel stx2k subtype and eae genes were detected among the 36 representative PEC isolates according to WGS. The results showed that 66.7 % (24/36) of sequenced PECs presented with multi-drug resistance (MDR) to β-lactamase 13.9 % (5/36), aminoglycoside 61.1 % (22/36), tetracycline 41.7 % (15/36) and quinolones 38.9 % (14/36). No colistin nor carbapenem resistance was detected. Sequence types (STs) associated with MDR in this study were: ST752, ST189, ST206, ST10, ST48 and ST38. The findings highlight the threat of zoonotic pathogens to close human contacts and the need for enhanced surveillance to mitigate the spread of MDR foodborne PECs.

Characterization of two lytic bacteriophages infecting carbapenem-resistant clinical Klebsiella pneumoniae in Dhaka, Bangladesh

Bacteriophages or bacteria infecting viruses are genetically diverse. Due to the emergence of antimicrobial-resistant bacteria, lytic bacteriophages are gaining enormous attention for treating superbug infections. Klebsiella pneumoniae is one of the eight most significant nosocomial pathogens and is addressed as a critical priority pathogen by WHO, requiring alternative treatment options. We reported two highly lytic bacteriophages, Klebsiella phage Kpn BM7 and the novel Klebsiella phage Kpn BU9, isolated from hospital wastewater and exhibiting lytic activity against different clinical isolates. Whole-genome analysis revealed that phages BM7 and BU9 belong to class Caudoviricetes. Phage BM7, with a genome length of 170,558 bp, is a member of the genus Marfavirus and the species Marfavirus F48. While phage BU9, with a genome length of 60,450 bp, remains unclassified. Neither phage harbors any lysogenic, toxin, or antimicrobial resistance genes. Both phages can steadily survive up to 40 °C and at pH 5-7. The optimal MOI was 0.1 for BM7 and 1 for BU9, with short latent periods of 10 and 25 min and burst sizes of 85 PFU/cell and 12 PFU/cell, respectively. This is the first carbapenem-resistant K. pneumoniae targeting lytic phages to be reported from Bangladesh. This study suggests that BM7 and BU9 are potential candidates for targeting carbapenem-resistant K. pneumoniae.

Extensive genomic characterization, pre-clinical probiotic evaluation, and safety analysis of Bifidobacterium longum subsp. longum BL21 isolated from infant feces

OBJECTIVE

This study aimed to evaluate the safety and probiotic properties of Bifidobacterium longum subsp. longum BL21 isolated from infant feces for use as a commercial probiotic strain.

METHODS

Whole-genome sequencing; physiological and biochemical assessments; enzymatic assays; metabolite, antibiotic sensitivity, cell adhesion and cytotoxicity, and tolerance tests; and a 14-day oral toxicity study were conducted.

RESULTS

BL21 exhibited genetic integrity, and its genome lacked genes related to antibiotic resistance or virulence. It was found to be non-pathogenic, had efficient carbohydrate metabolism and mucin degradation ability, and was free from biogenic amines. It also showed susceptibility to antibiotics, strong cell adhesion, and resilience to adverse conditions. The aforementioned results confirm that BL21 is a functional probiotic strain with genetic stability, enzymatic capabilities, and non-pathogenic properties that mean it is safe for oral consumption, demonstrating that it is a promising candidate for probiotic applications.

CONCLUSION

The study demonstrates that BL21 is a genetically stable, non-pathogenic probiotic strain with metabolic potential. The strain lacks virulence and antibiotic resistance genes, and its resilience to gastrointestinal conditions, as well as the results of the 14-day oral toxicity study, suggest that BL21 is safe for oral consumption. However, further long-term studies and clinical trials are needed to confirm its safety and efficacy for therapeutic use.

Comprehensive genomic analysis and evaluation of in vivo and in vitro safety of Heyndrickxia coagulans BC99

This study aimed to evaluate the safety profile and beneficial effects of Heyndrickxia coagulans strain BC99 (BC99) for potential use in functional foods and pharmaceuticals. We began with whole genome sequencing of BC99, followed by a comprehensive safety assessment comprising genome analysis, hemolysis, cytotoxicity, antibiotic susceptibility tests, and cell adhesion and tolerance studies, along with acute and subacute oral toxicity studies in animal models. BC99 was isolated from a well-characterized collection originating from the feces of a healthy infant. Our results indicated no hemolytic activity on Columbia blood agar plates and broad antibiotic sensitivity, including to gentamicin, ampicillin, chloramphenicol, ciprofloxacin, and others. Cytotoxicity testing confirmed no adverse effects on HT-29 cells and significant adhesive properties to intestinal epithelial cells. Tolerance tests demonstrated over 90% viability of BC99 under simulated gastrointestinal conditions. In vivo studies in mice and rats confirmed the absence of adverse effects following oral administration. Collectively, these findings support BC99's robust tolerance to gastrointestinal environments, strong adhesion capabilities, and a broad spectrum of antibiotic resistance, underlining its potential as a safe and effective agent for gut microbiota modulation and host health enhancement.

Streptococcus canis transcriptomic modifications in host cell entry environments of human keratinocytes

BACKGROUND

Streptococcus canis is a commensal bacterium in companion animals. This microorganism can infect humans who have been in deep contact with or bitten by pet dogs, suggesting that the skin/soft tissue is one of infection entry sites. To understand pathological process in human cells, we aimed to determine S. canis transcriptomic changes in invasive environments of human keratinocytes.

METHODS

We selected one isolate from candidates with whole-genome sequences, based on re-obtained cell invasion ability (CIA) data into human keratinocytes along with bacterial cytotoxicity. RNA-sequencing was conducted for the samples at baselines and 2 h/5 hr post-inoculation using NovaSeq 6000. Global/differential gene expression analyses [principal component analysis (PCA)/k-means clustering analysis/differentially expressed gene (DEG) analyses] were performed. We classified DEGs into their functional categories. To validate transcriptomic results, we did quantitative reverse transcription-polymerase chain reaction (qRT-PCR) assays.

RESULTS

FU1 isolate was selected from seven candidates, based on re-obtained CIA data with less cytotoxicity. Total read bases of 6.17-9.02 Gbp were obtained by RNA-sequencing. PCA and k-means clustering analysis indicated clustering according to their inoculation times. Volcano plots and Venn diagrams revealed that S. canis invasion into keratinocytes produced altered distributions of many genes. Gene ontology enrichment analysis showed most of the gene expressions were downregulated. DEG functional analysis showed the downregulated DEGs belonging to energy production and conversion/carbohydrate transport and metabolism/amino acid transport and metabolism/nucleotide transport and metabolism, with the upregulated DEGs belonging to transcription. qRT-PCR assays for downregulated/upregulated expressions of four genes (pgk-slo/opuAA-kdpB) validated transcriptomic results.

CONCLUSION

Our observations suggest that S. canis can downregulate its metabolism-associated gene expressions in human keratinocyte environments. The observed gene expression changes can imply the latent infection in human cells. Further investigation is needed to elucidate the underlying mechanisms for the latent infection.

Unveiling the overlooked small-sized microbiome in river ecosystems

Enriching microorganisms using a 0.22-μm pore size is a general pretreatment procedure in river microbiome research. However, it remains unclear the extent to which this method loses microbiome information. Here, we conducted a comparative metagenomics-based study on microbiomes with sizes over 0.22 μm (large-sized) and between 0.22 μm and 0.1 μm (small-sized) in a subtropical river. Although the absolute concentration of small-sized microbiome was about two orders of magnitude lower than that of large-sized microbiome, sequencing only large-sized microbiome resulted in a significant loss of microbiome diversity. Specifically, the microbial community was different between two sizes, and 347 genera were only detected in small-sized microbiome. Small-sized microbiome had much more diverse viral community than large-sized fraction. The viruses had abundant ecological functions and were hosted by 825 species of 169 families, including pathogen-related families. Small-sized microbiome had distinct antimicrobial resistance risks from large-sized microbiome, showing an enrichment of eight antibiotic resistance gene (ARG) types as well as the detection of 140 unique ARG subtypes and five enriched risk rank I ARGs. Draft genomes of five major resistant pathogens having diverse ecological and pollutant-degrading functions were only assembled in small-sized microbiome. These findings provide novel insights into river ecosystems, and highlight the overlooked small-sized microbiome in the environment.

Comparative genomic analysis of ESBL-selected and non-selected Escherichia coli in Australian wastewater: Elucidating differences in diversity, antimicrobial resistance, and virulence profiles

Extended-spectrum β-lactamases (ESBLs)-producing E. coli have been proposed as an indicator bacterium for antimicrobial resistance (AMR) surveillance within a OneHealth framework. However, it is important to understand the effects and potential biases ESBL-selection has on E. coli populations. Utilising whole genome sequencing, this study compared 80 ESBL-selected E. coli isolates with 201 non-selected isolates from Australian wastewater. The findings revealed significant variations between these cohorts in genetic diversity, AMR profiles, and carriage of virulence-associated genes (VAGs), plasmids, and the transmissible Locus of Stress Tolerance (tLST), a genomic island that imparts resistance to extreme heat and chlorination. The study highlights the predominance of certain sequence types (STs), particularly ST131 (75 % clade A), in ESBL-selected isolates (40 % vs 2 %) and overall the ESBL-selected isolates were largely multidrug-resistant (MDR), predominantly carrying genes for resistance to aminoglycosides, extended-spectrum β-lactams, fluoroquinolone, macrolides, sulphonamides/trimethoprim, and tetracyclines. The ESBLs identified were almost exclusively blaCTX-M genes, most commonly blaCTX-M-15 > blaCTX-M-27 > blaCTX-M-14. These were predominately carried on IncF plasmids or chromosomally (always ISEcp1 associated), in equal numbers. In contrast, 80 % of non-selected isolates carried no acquired ARGs, and none carried blaCTX-M genes. In both cohorts, extraintestinal pathogenic E. coli (ExPEC) was the dominate pathotype (35 % total) with few (4 % total) intestinal pathogenic E. coli pathotypes identified (aEPEC > ETEC > EAEC). Nevertheless, some clinically important genes were only identified in the non-selected group, namely tigecycline-resistance gene tet(X4) and AmpC ESBL blaCMY-2. Additionally, the presence of tLST, associated with higher metal resistance gene carriage (Ag, As, Cu, Hg, Ni), in a substantial portion of non-selected isolates (20 % vs 0 %), underscores environmental pressures shaping bacterial populations in wastewater ecosystems. These insights are important for developing comprehensive, less biased genomic surveillance strategies to understand and manage public health threats posed by pathogenic E. coli and AMR.

Immunoinformatic approach for multi-epitope vaccine design against Staphylococcus aureus based on hemolysin proteins

Staphylococcus aureus is a common bacterium that causes a variety of infections in humans. This microorganism produces several virulence factors, including hemolysins, which contribute to its disease-causing ability. The treatment of S. aureus infections typically involves the use of antibiotics. However, the emergence of antibiotic-resistant strains has become a major concern. Therefore, vaccination against S. aureus has gained attention as an alternative approach. Vaccination has the advantage of stimulating the immune system to produce specific antibodies that can neutralize bacteria and prevent infection. However, developing an effective vaccine against S. aureus has proven to be challenging. This study aimed to use in silico methods to design a multi-epitope vaccine against S. aureus infection based on hemolysin proteins. The designed vaccine contained four B-cell epitopes, four CTL epitopes, and four HTL epitopes, as well as the ribosomal protein L7/L12 and pan-HLA DR-binding epitope, included as adjuvants. Furthermore, the vaccine was non-allergenic and non-toxic with the potential to stimulate the TLR2-, TLR-4, and TLR-6 receptors. The predicted vaccine exhibited a high degree of antigenicity and stability, suggesting potential for further development as a viable vaccine candidate. The population coverage of the vaccine was 94.4 %, indicating potential widespread protection against S. aureus. Overall, these findings provide valuable insights into the design of an effective multi-epitope vaccine against S. aureus infection and pave the way for future experimental validations.

Whole-genome sequencing-based analysis of Brucella species isolated from ruminants in various regions of Türki̇ye

BACKGROUND

Brucellosis, a zoonotic disease in Türkiye, which has significant direct and indirect impacts on the healthcare system and livestock. This study, which aimed to investigate the differences among Brucella spp. isolates originating from different regions of Türkiye, for implications for public health and veterinary medicine.

METHOD

Twenty-one isolates from ruminants and two isolates from humans obtained from various regions of Türkiye were utilized in the study. The isolates were identified and biotyped using traditional microbiological procedures, and whole-genome sequencing (WGS) was performed. This was followed by single nucleotide polymorphism (SNP)--based phylogenetic analysis and WGS-based analysis of virulence and resistance genes. Additionally, phenotypic antimicrobial resistance and phage susceptibilities were determined. The obtained data were then compared for concordance, ensuring the validity and reliability of the results.

RESULTS

Our study, employing culture methods, polymerase chain reaction (PCR), and WGS analyses, identified 11 Brucella melitensis (bv 3 (n = 9), one each bv 1 and bv 2) and 12 B. abortus (bv 3 (n = 11), bv 9 (n = 1)) isolates All B. abortus isolates were of bovine origin, while the B. melitensis isolates were from sheep (n = 7), goat (n = 1), ram (n = 1), and humans (n = 2). In the whole-genome SNP-based phylogenetic tree, all B. melitensis strains were found to be of the IIb subtype of genotype II associated with the Eastern Mediterranean lineage. Ten different genotypes were identified in the SNP analysis of the isolates, with a maximum SNP difference of 278 and a minimum SNP difference of 4 among these genotypes. According to the WGS-SNP-based phylogenetic tree of B. abortus isolates, they were grouped in clade C1. In the SNP analysis, where ten different genotypes were identified, the SNP difference among these genotypes was a maximum of 316 and a minimum of 6. In the in silico MLST analysis performed with WGS data, B. melitensis isolates were identified as ST8 and ST102 genotypes, while B. abortus isolates were identified as ST2 and ST3 genotypes. The dominant genotypes were ST8 for B. melitensis and ST2 for B. abortus, respectively. Virulence gene analysis conducted based on WGS data of the 23 B. abortus and B. melitensis isolates revealed 43 virulence gene-associated regions in all strains, irrespective of species, host, or isolation year. Although classical resistance-related genes were not detected by WGS-based antimicrobial resistance gene analysis, phenotypic resistance analysis revealed resistance to azithromycin, rifampin, and trimethoprim/sulfamethoxazole in B. abortus and B. melitensis isolates.

CONCLUSION

Both B. melitensis and B. abortus were circulating species in animals and human. The dominant genotypes were ST8 for B. melitensis and ST2 for B. abortus, respectively. All B. melitensis strains were found to be of the IIb subtype of genotype II associated with the Eastern Mediterranean lineage, while B. abortus isolates, they were grouped in clade C1. Further, a comprehensive study with a sufficient number of isolates covering all regions of Türkiye would provide more accurate information about the current epidemiological situation in the country.

Metagenome mining divulges virulent and multidrug resistant Pseudomonas aeruginosa ST242 and Klebsiella michiganensis ST∗1b23 coinfecting an 8-month-old meningitis infant under ICU in Kampala, Uganda, East Africa

Pediatric meningitis is a global health problem, with insufficiently known pathogens and antibiotic resistance (AMR) especially in low-resource settings. Here, we sought to uncover the virulence and AMR of pathogens associated with infant meningitis, treated with ceftriaxone, in Kampala, Uganda. In a bid to isolate Klebsiella oxytoca, we coincidentally recovered a co-culture and challenged it with antibiotic susceptibility testing (AST) on a panel of 14 antibiotics. We then combined metagenome binning with antiSMASH/PRISM genome mining to unveil the pathogens, their virulence and molecular targets in relation to meningitis. From AST, the co-culture exhibited resistance to multiple aminoglycosides, fluroquinolones, and β-lactams, including ceftriaxone, the inherently used drug. From metagenome annotation, the first bin was identified as Pseudomonas aeruginosa ST242 and the second as Klebsiella michiganensis ST∗1b23. Among others, P. aeruginosa ST242 virulence factors include type 3 and type 6 secretion systems, biofilm, and nonribosomal peptides (NRPs) of the pyoverdine synthase operon, targeting claudin-5, a component of the tight junctions of the blood-brain barrier (BBB). The P. aeruginosa ST242 genome portrays intrinsic resistance to beta-lactamases (blaOXA-50 and blaPAO), aminoglycosides [aph(3')-IIb)], fluoroquinolones (crpP), tetracycline (tmexD2) and fosfomycin (fosA), among others. From K. michiganensis ST∗1b23 genome mining we elucidated a yersiniabactin-related metabolite, targeting the ligand-binding domain of the human polymeric immunoglobulin receptor (pIgR) and other components of the BBB. The K. michiganensis ST∗1b23 chromosome encodes the genes blaOXY-1 and OqxA/B, conferring resistance to β-lactams, fluoroquinolones, and trimethoprim respectively. Notably, we found one frameshift and nine substitution mutations conferring carbapenem and cephalosporin resistance mechanisms. Overall, our findings strongly suggest coinfection and a mechanistic crosstalk between P. aeruginosa ST242 and K. michiganensis ST∗1b23 in the pathogenesis of meningitis in this case. Importantly, ceftriaxone could be an inappropriate treatment choice for these pathogens. Hence, serious surveillance and experimental studies will improve the management of pediatric meningitis.

Large-scale genomic analysis of Elizabethkingia anophelis

The recent emergence of Elizabethkingia anophelis as a human pathogen is a major concern for global public health. This organism has the potential to cause severe infections and has inherent antimicrobial resistance. The potential for widespread outbreaks and rapid global spread highlights the critical importance of understanding the biology and transmission dynamics of this infectious agent. We performed a large-scale analysis of available 540 E. anophelis, including one novel strain isolated from raw milk and sequenced in this study. Pan-genome analysis revealed an open and diverse pan-genome in this species, characterized by the presence of many accessory genes. This suggests that the species has a high level of adaptability and can thrive in a variety of environments. Phylogenetic analysis has also revealed a complex population structure, with limited source-lineage correlation. We identified diverse antimicrobial resistance factors, including core-genome and accessory ones often associated with mobile genetic elements within specific lineages. Mobilome analysis revealed a dynamic landscape primarily composed of genetic islands, integrative and conjugative elements, prophage elements, and small portion of plasmids emphasizing a complex mechanism of horizontal gene transfer. Our study underscores the adaptability of E. anophelis, characterized by a diverse range of antimicrobial resistance genes, putative virulence factors, and genes enhancing fitness. This adaptability is also supported by the organism's ability to acquire genetic material through horizontal gene transfer, primarily facilitated by mobile genetic elements such as integrative and conjugative elements (ICEs). The potential for rapid evolution of this emerging pathogen poses a significant challenge to public health efforts.

VFDB 2025: an integrated resource for exploring anti-virulence compounds

With the escalating crisis of bacterial multidrug resistance, anti-virulence therapeutic strategies have emerged as a highly promising alternative to conventional antibiotic treatments. Anti-virulence compounds are specifically designed to target virulence factors (VFs), disarming pathogens without affecting bacterial growth and thus reduce the selective pressure for resistance development. However, due to the complexity of bacterial pathogenesis, no anti-virulence small molecules have been approved for clinical use thus far, despite the documentation of hundreds of potential candidates. To provide valuable reference resources for drug design, repurposing, and target selection, the virulence factor database (VFDB, http://www.mgc.ac.cn/VFs/) has systematically collected public data on anti-virulence compounds through extensive literature mining, and further integrated this information with its existing knowledge of bacterial VFs. To date, the VFDB has curated a comprehensive dataset of 902 anti-virulence compounds across 17 superclasses reported by 262 studies worldwide. By cross-linking the current knowledge of bacterial VFs with information on relevant compounds (e.g. classification, chemical structure, molecular targets and mechanisms of action), the VFDB aims to bridge the gap between chemists and microbiologists, providing crucial insights for the development of innovative and effective antibacterial therapies to combat bacterial infections and address antibiotic resistance.

Deep learning revealed the distribution and evolution patterns for invertible promoters across bacterial lineages

Invertible promoters (invertons) are crucial regulatory elements in bacteria, facilitating gene expression changes under stress. Despite their importance, their prevalence and the range of regulated gene functions are largely unknown. We introduced DeepInverton, a deep learning model that identifies invertons across a broad phylogenetic spectrum without using sequencing reads. By analyzing 68 733 bacterial genomes and 9382 metagenomes, we have uncovered over 200 000 nonredundant invertons and have also highlighted their abundance in pathogens. Additionally, we identified a post-Cambrian Explosion increase of invertons, paralleling species diversification. Furthermore, we revealed that invertons regulate diverse functions, including antimicrobial resistance and biofilm formation, underscoring their role in environmental adaptation. Notably, the majority of inverton identifications by DeepInverton have been confirmed by the in vitro experiments. The comprehensive inverton profiles have deepened our understanding of invertons at pan-genome and pan-metagenome scales, enabling a broad spectrum of applications in microbial ecology and synthetic biology.

First comparative genomics analysis of Corynebacterium auriscanis

BACKGROUND

Corynebacterium auriscanis is a bacterial species frequently isolated from dogs with external otitis or dermatitis and a zoonotic pathogen transmitted by dog bite. It is considered an opportunistic pathogen, but its pathogenicity mechanisms are poorly studied. Comparative genomics can identify virulence and niche factors that could contribute to understanding its lifestyle.

OBJECTIVES

The objectives of this project was to compare genomes of C. auriscanis to identify genes related to its virulence and lifestyle.

METHODS

The genome of strain 32 was sequenced using Illumina HiSeq 2500 (Illumina, CA, USA) and assembled using Unicycler. The two other non-redundant genomes from the same species available in GenBank were included in the analysis. All genomes were annotated and checked for taxonomy, assembly quality, mobile elements, CRISPR-Cas systems, and virulence and antimicrobial resistance genes. The virulence genes in the three genomes were compared to the ones from other pathogens commonly isolated with C. auriscanis.

FINDINGS

The species has 42 virulence factors that can be classified as niche factors, due to the absence of true virulence factors found in primary pathogens. The gene rbpA could confer basal levels of resistance to rifampin.

MAIN CONCLUSIONS

The absence of true virulence factors in the three genomes suggests C. auriscanis has an opportunistic pathogen lifestyle.

Virulome and phylogenomic profiling of a novel Burkholderia pseudomallei strain from an Indian clinical isolate

Highly pathogenic Burkholderia pseudomallei is the causative agent of melioidosis, a neglected tropical disease endemic in Southeast Asian tropical region. This bacterium encompasses diverse virulence factors which further undergo dynamic gene-expression flux as it transits through distinct environmental niches within the host which may lead to manifestation of differential clinical symptoms. B. pseudomallei, is classified as a Tier 1 select agent in the United States and regarded as a risk group 3 organism in India with the potential to be used as bioweapon. Considering these facts, it is vital to uncover both physiological and genetic heterogeneity of B. pseudomallei, particularly to identify any novel virulence factors that may contribute to pathogenicity. B. pseudomallei strain CM000113 was isolated from a clinical case in India, characterized it for its physiological, biochemical, and prominently genetic traits through WGS. It has a type 2 morphotype with faster doubling time and high biofilm producing capacity as compared to Pseudomonas aeruginosa. The genome size is 7.3 Mbp and it is phylogenetically close to B. pseudomallei strain Mahidol 1106a and Burkholderia mallei Turkey 2. We observed genetic heterogeneity, as key virulence factors that were identified shows sequence dissimilarity with reference strains. Additionally, presence of genomic islands, harbouring two virulence factors, GmhA and GmhB2, associated with pathogenesis indicates possibility of horizontal gene transfer. These results emphasize the need for an extensive study focusing the genome of B. pseudomallei and its associated heterogeneity, to identify molecular biomarkers aiding to develop point-of-care diagnostic kits for early diagnosis of melioidosis.

Whole genome sequencing reveals antibiotic resistance pattern and virulence factors in Klebsiella quasipneumoniae subsp. Similipneumoniae from Hospital wastewater in South-West, Nigeria

Klebsiella quasipneumoniae is a distinct species from K. pneumoniae, even though it is sometimes mistaken phenotypically for the latter in clinical situations. K. quasipneumoniae is a pathogen and this study aims at understanding the genomic antibiotic resistance and virulence characteristics of Klebsiella quasipneumoniae subsp. similipneumoniae (B105 strain) isolated from tertiary hospital wastewater and the potential risks associated with its environmental spread. The Illumina platform was used for whole-genome sequencing (WGS), the generated raw reads (de novo) was assembled using RAPT NCBI, while other standardized bioinformatics tools were utilized to validate and examine the landscape of the genome's antibiotic resistance and virulence factors. The K. quasipneumoniae subsp. similipneumoniae (B105 strain), belonged to sequence type 1422 and was resistant to ampicillin, amoxicillin-clavulanic acid, ceftazidime, cefepime, meropenem, tetracycline, but susceptible to gentamicin. The annotated genome acknowledged the presence of blaOKP-B-2, ompK 36, fosA5, oqxAB, virulence genes responsible for capsule formation, lipopolysaccharide, iron uptake aerobactin (iutA), salmochelins (iroE, iroN), enterobactin siderophore, efllux pump (acrA, acrB) adherence, (mrkC, mrkD, and fimD) and two plasmids replicon IncFIB(K) and IncR. The study resonates the inadequacy of conventional microbiological identification methods to distinguish K. pneumoniae and K. quasipneumoniae and at the same time heightens the importance of using a genomic platform to extol the identity of K. quasipneumoniae subsp. similipneumoniae strain. Furthermore, the peculiarities of the acquired antimicrobial resistance and virulence genes, in this strain, are a potential risk to the environment.

Deciphering the Genetic Architecture of Staphylococcus warneri Prophage vB_G30_01: A Comprehensive Molecular Analysis

The current knowledge of Staphylococcus warneri phages is limited, with few genomes sequenced and characterized. In this study, a prophage, vB_G30_01, isolated from Staphylococcus warneri G30 was characterized and evaluated for its lysogenic host range. The phage was studied using transmission electron microscopy and a host range. The phage genome was sequenced and characterized in depth, including phylogenetic and taxonomic analyses. The linear dsDNA genome of vB_G30_01 contains 67 predicted open reading frames (ORFs), classifying it within Bronfenbrennervirinae. With a total of 10 ORFs involved in DNA replication-related and transcriptional regulator functions, vB_G30_01 may play a role in the genetics and transcription of a host. Additionally, vB_G30_01 possesses a complete set of genes related to host lysogeny and lysis, implying that vB_G30_01 may influence the survival and adaptation of its host. Furthermore, a comparative genomic analysis reveals that vB_G30_01 shares high genomic similarity with other Staphylococcus phages and is relatively closely related to those of Exiguobacterium and Bacillus, which, in combination with the cross-infection assay, suggests possible cross-species infection capabilities. This study enhances the understanding of Staphylococcus warneri prophages, providing insights into phage-host interactions and potential horizontal gene transfer.

Unveiling MurM inhibitors in Enterococcus faecalis V583: a promising approach to tackle antibiotic resistance

Enterococcus faecalis is commonly found in the GI tract of humans and animals. It causes various infections, especially in hospital environments, and shows growing antibiotic resistance. This study utilized a subtractive proteomics approach to find out the potential drug targets in E. faecalis. Unique metabolic pathways were analysed and compared to the host to minimize adverse effects. Among twenty nine pathogenic specific and seventy three host-pathogen common pathways identified using the KEGG database, sixty seven essential proteins were found through the DEG BLAST search. PSORTB predicted that forty cytoplasmic proteins could be suitable as druggable targets. Further analysis identified fourteen proteins with virulence properties using the VFDB BLAST. Among these, seven proteins with more than ten antigenic sites were subjected to DrugBank BLAST, identifying three novel and four existing drug targets. One of the crucial drug targets, MurM, was selected due to its critical role in peptidoglycan biosynthesis. The reason for selecting MurM is crucial for addressing antibiotic resistance, disrupting bacterial cell wall synthesis, and attaining selective antimicrobial activity. MurM belongs to the mixed αβ class with two functional domains. The possible binding site residues of MurM are Trp31, Lys35, Trp38, Arg215, and Tyr219. Virtual screening identified potential lead candidates for MurM, and four were selected based on their physiochemical, pharmacokinetic, and structural properties. This study provides valuable insights into identifying and analysing a potential drug target, the MurM protein, and its inhibitors in E. faecalis V583.

Genetic characterization of Neisseria meningitidis isolates recovered from patients with invasive meningococcal disease in Lithuania

Introduction

Neisseria meningitidis is a gram-negative bacterium responsible for life-threatening invasive infections known as invasive meningococcal disease and is associated with high fatality rates and serious lifelong disabilities among survivors.

Methods

This study aimed to characterize N. meningitidis isolates cultured from blood and cerebrospinal fluid collected between 2009 and 2021 in Lithuania, assess their genomic relationships with European strains, and evaluate the possibility of using a cost-effective method for strain characterization, thus improving the national molecular surveillance of invasive meningococcal disease. In total, 321 N. meningitidis isolates were collected and analyzed using multilocus restriction typing (MLRT). Amplification of the penA gene and restriction fragment length polymorphism analysis were performed to identify the modified penA genes. Based on the MLRT genotyping results, we selected 10 strains for additional analysis using whole-genome sequencing. The sequenced genomes were incorporated into a dataset of publicly available N. meningitidis genomes to evaluate genomic diversity and establish phylogenetic relationships within the Lithuanian and European circulating strains.

Results

We identified 83 different strains using MLRT genotyping. Genomic diversity of N. meningitidis genomes analysed revealed 21 different sequence types (STs) circulating in Lithuania. Among these, ST34 was the most prevalent. Notably, three isolates displayed unique combinations of seven housekeeping genes and were identified as novel STs: ST16969, ST16901, and ST16959. The analyzed strains were found to possess virulence factors not commonly found in N. meningitidis. Six distinct penA profiles were identified, each with different frequencies. In the present study, we also identified N. meningitidis strains with new penA, NEIS0123, NEIS1320, NEIS1525, NEIS1600, and NEIS1753 loci variants. In our study, using the cgMLST scheme, Minimum Spanning Tree (MST) analysis did not identify significant geographic relationships between Lithuanian N. meningitidis isolates and strains from Europe.

Discussion

Discussion: To our knowledge, this is the first study to employ whole genome sequencing (WGS) method for a comprehensive genetic characterization of invasive N. meningitidis isolates from Lithuania. This approach provides a more detailed and precise analysis of genomic relationships and diversity compared to prior studies relying on traditional molecular typing methods and antigen analysis.

Isolation, characterization, and application of a novel Pseudomonas fluorescens phage vB_PF_Y1-MI in contaminated milk

The food industry has incurred substantial losses from contamination by Pseudomonas fluorescens, emphasizing the critical importance of implementing effective control strategies. Phages are potential sterilizers due to their specific killing abilities and the difficulty bacteria face in developing resistance. However, a significant barrier to their development is the lack of diversity among phage types. In this study, we characterized a novel lytic P. fluorescens phage, named vB_PF_Y1-MI. Phage vB_PF_Y1-MI displayed a latent period of nearly 10 min and a high burst size of 1493 PFU/cell. This phage showed good activity over a wide range of temperature (up to 70 °C) and pH (3-12). The genome of phage vB_PF_Y1-MI spans 93,233 bp with a GC content of 45%. It encompasses 174 open-reading frames and 19 tRNA genes, while no lysogeny or virulence-associated genes were detected. Phylogenetic analysis positions it as a novel unassigned evolutionary lineage within the Caudoviricetes class among related dsDNA phages. Our study provides foundational insights into vB_PF_Y1-MI and emphasizes its potential as an effective biological control agent against P. fluorescens. This research offers crucial theoretical groundwork and technical support for subsequent efforts in preventing and controlling P. fluorescens contamination.

Complete genome sequences of 12 lytic phages against multidrug-resistant Enterococcus faecalis

We report the genome sequences of 12 Enterococcus faecalis phages isolated in Kenya, belonging to the genus Copernicusvirus, Efquatrovirus, Saphexavirus, and Kochikohdavirus. They have double-stranded DNA with lengths varying from 17,979 to 147,374 bp and G+C content from 33.14% to 40.05%. The genomes contain 28-250 coding sequences.

Clinical sequelae of gut microbiome development and disruption in hospitalized preterm infants

Aberrant preterm infant gut microbiota assembly predisposes to early-life disorders and persistent health problems. Here, we characterize gut microbiome dynamics over the first 3 months of life in 236 preterm infants hospitalized in three neonatal intensive care units using shotgun metagenomics of 2,512 stools and metatranscriptomics of 1,381 stools. Strain tracking, taxonomic and functional profiling, and comprehensive clinical metadata identify Enterobacteriaceae, enterococci, and staphylococci as primarily exploiting available niches to populate the gut microbiome. Clostridioides difficile lineages persist between individuals in single centers, and Staphylococcus epidermidis lineages persist within and, unexpectedly, between centers. Collectively, antibiotic and non-antibiotic medications influence gut microbiome composition to greater extents than maternal or baseline variables. Finally, we identify a persistent low-diversity gut microbiome in neonates who develop necrotizing enterocolitis after day of life 40. Overall, we comprehensively describe gut microbiome dynamics in response to medical interventions in preterm, hospitalized neonates.

A core microbiome signature as an indicator of health

The gut microbiota is crucial for human health, functioning as a complex adaptive system akin to a vital organ. To identify core health-relevant gut microbes, we followed the systems biology tenet that stable relationships signify core components. By analyzing metagenomic datasets from a high-fiber dietary intervention in type 2 diabetes and 26 case-control studies across 15 diseases, we identified a set of stably correlated genome pairs within co-abundance networks perturbed by dietary interventions and diseases. These genomes formed a "two competing guilds" (TCGs) model, with one guild specialized in fiber fermentation and butyrate production and the other characterized by virulence and antibiotic resistance. Our random forest models successfully distinguished cases from controls across multiple diseases and predicted immunotherapy outcomes through the use of these genomes. Our guild-based approach, which is genome specific, database independent, and interaction focused, identifies a core microbiome signature that serves as a holistic health indicator and a potential common target for health enhancement.

Genomic analysis of Enterobacteriaceae from colorectal cancer patients at a tertiary hospital in Ghana: a case-control study

Colorectal cancer (CRC) is a severe gastrointestinal cancer and a leading cause of cancer-related deaths in Ghana. The potential role of gut Enterobacteriaceae in the increasing incidence of CRC in Ghana is yet to be thoroughly investigated. In this study, Enterobacteriaceae from CRC patients and healthy control participants were analyzed by whole genome sequencing to identify genomic features that are associated with CRC. Socio-demographic data showed a significant association between age and alcohol consumption and CRC. Escherichia coli was the most abundant Enterobacteriaceae isolated from the study participants and they were predominantly intestinal commensals. Escherichia coli isolates belonging to phylogroup D encoded the highest number of virulence genes. The agn43 and int genes were widespread in Escherichia coli isolates from the CRC patients. Multilocus sequence types of potentially pathogenic Escherichia coli from the CRC patients also encoded genes involved in aggregation, adherence and biofilm formation. The ampC2 and ampH antimicrobial resistance genes were also widespread in the genome of the Escherichia coli isolates. This study highlights the virulence tendencies of Escherichia coli from CRC patients and their ability to transfer virulence determinants to other Enterobacteriaceae residing in the gut.

Expression of the locus of enterocyte effacement genes during the invasion process of the atypical enteropathogenic Escherichia coli 1711-4 strain of serotype O51:H40

Atypical enteropathogenic Escherichia coli (aEPEC) is a significant cause of diarrhea in low- and middle-income countries. Certain aEPEC strains, including the Brazilian representative strain of serotype O51:H40 called aEPEC 1711-4, can use flagella to attach to, invade, and persist in T84 and Caco-2 intestinal cells. It can also translocate from the gut to extraintestinal sites in a rat model. Although various aspects of the virulence of this strain were studied and the requirement of a type III secretion system for the efficiency of the invasion process was demonstrated, the expression of the locus of enterocyte effacement (LEE) genes during the invasion and intracellular persistence remains unclear. To address this question, the expression of flagella and the different LEE operons was evaluated during kinetic experiments of the interaction of aEPEC 1711-4 with enterocytes in vitro. The genome of the strain was also sequenced. The results showed that flagella expression remained unchanged, but the expression of eae and escJ increased during the early interaction and invasion of aEPEC 1711-4 into Caco-2 cells, and there was no change 24 h post-infection during the persistence period. The number of actin accumulation foci formed on HeLa cells also increased during the 6-h analysis. No known gene related to the invasion process was identified in the genome of aEPEC 1711-4, which was shown to belong to the global EPEC lineage 10. These findings suggest that the LEE components and the intimate adherence promoted by intimin are necessary for the invasion and persistence of aEPEC 1711-4, but the detailed mechanism needs further study.IMPORTANCEAtypical enteropathogenic Escherichia coli (aEPEC) is a major cause of diarrhea, especially in low- and middle-income countries, like Brazil. However, due to the genome heterogeneity of each clonal group, it is difficult to comprehend the pathogenicity of this strain fully. Among aEPEC strains, 1711-4 can invade eukaryotic cells in vitro, cross the gut barrier, and reach extraintestinal sites in animal models. By studying how different known aEPEC virulence factors are expressed during the invasion process, we can gain insight into the commonalities of this phenotype among other aEPEC strains. This will help in developing preventive measures to control infections caused by invasive strains. No known virulence-encoding genes linked to the invasion process were found. Nevertheless, additional studies are still necessary to evaluate the role of other factors in this phenotype.

Antibiotic resistomes and ecological risk elimination in field-scale constructed wetland revealed by integrated metagenomics and metatranscriptomics

Constructed wetlands (CWs) are identified as significant ecological systems for the potential control of antibiotic resistance genes (ARGs) in the environment. However, the precise mechanisms governing removal, persistence, expression, and associated risks of ARGs during wetland treatment remain poorly understood. In this study, the distribution, mobility, expression, and hosts of ARGs in water, sediments, and plants of a field-scale CW and its parallel natural river were systematically investigated through metagenomic and metatranscriptomic approaches. Results showed that both the abundance and diversity of ARGs in water gradually decreased along the way of CW, reaching a final abundance removal rate of 72.28 % in the effluent. Source tracking analysis indicted that the reduction of ARGs in water was mainly achieved by the dynamic accumulation of ARGs in sediments and plants of the CW. Proteobacteria were identified as primary hosts for ARGs, particularly in sediments and plants during CW treatment. Moreover, although ESKAPE pathogens carrying multiple ARGs persisted in all media throughout the CW treatment, ARG expression levels and risk of water were also significantly decreased after CW treatment. Collectively, our comprehensive multi-omics study would enhance the understanding of ARG removal by CWs, offering insights for controlling antimicrobial resistance in wastewater treatment system.

Three prolonged outbreaks of metallo-β-lactamase-producing Pseudomonas aeruginosa in an Upper Austrian hospital, 2017-2023

In spring 2022, an increase in metallo-β-lactamase-producing Pseudomonas aeruginosa (MBL-Pa) infections was detected in a hospital in Upper Austria. To identify the source of infection and to stop further transmissions, an epidemiological outbreak investigation including whole-genome sequencing (WGS)-based typing was conducted. The final case definition included cases admitted to the hospital between 2020 and 2023 with an MBL-Pa in one of the three genomic clusters identified. In addition, the investigation was extended to include historical cases from 2017. Core genome multilocus sequence typing was performed to assess the genetic relatedness between the isolates. Fifty-four clinical P. aeruginosa isolates and eight P. aeruginosa isolates from the hospital environment were obtained. All but nine isolates grouped into one of three genomic clusters (ST235/blaVIM-1, ST111/blaVIM-2, or ST621/blaIMP-13), which were considered to be distinct, prolonged outbreaks involving 47 out of 52 cases. The most likely source of infection for cluster 1 (ST111/blaVIM-2) and cluster 2 (ST235/blaVIM-1) was sinks in the intensive care unit (ICU) washroom. Cluster 3 clone (ST621/blaIMP-13) could have originated in the urology ward in 2020 and then spread to the ICU years later. However, the nosocomial origin of this clone could not be proven. In March 2023, following the implementation of control measures (gowning, patient isolation, screening, and daily disinfection), no further MLB-Pa was detected, and the outbreaks were considered to be over. As ICUs play an important role in the transmission of P. aeruginosa, emphasis should be placed on genomic surveillance, infection prevention, and control in such wards.

IMPORTANCE

The significance of our work lies in the successful resolution of three prolonged outbreaks of MBL-Pa infections in a hospital in Upper Austria. Through a comprehensive epidemiological investigation coupled with WGS-based typing of P. aeruginosa isolates, the study identified three distinct genomic clusters responsible for prolonged outbreaks involving 47 cases. The investigation pinpointed sinks in the ICU washroom as the likely source of infection for two of the clusters. The study demonstrates the effectiveness of control measures such as hand hygiene, gowning, patient isolation, screening, and disinfection in stopping further transmission and bringing the outbreaks to a close. This underscores the critical role of genomic surveillance and control measures, particularly in high-risk settings like ICUs, in reducing nosocomial transmission of MBL-Pa infections.

Genome characterization of Shewanella algae in Hainan Province, China

Shewanella algae is an emerging marine zoonotic pathogen. In this study, we first reported the Shewanella algae infections in patients and animals in Hainan Province, China. Currently, there is still relatively little known about the whole-genome characteristics of Shewanella algae in most tropical regions, including in southern China. Here, we sequenced the 62 Shewanella algae strains isolated from Hainan Province and combined with the whole genomes sequences of 144 Shewanella algae genomes from public databases to analyze genomic features. Phylogenetic analysis revealed that Shewanella algae is widely distributed in the marine environments of both temperate and tropical countries, exhibiting close phylogenetic relationships with genomes isolated from patients, animals, and plants. Thereby confirming that exposure to marine environments is a risk factor for Shewanella algae infections. Average nucleotide identity analysis indicated that the clonally identical genomes could be isolated from patients with different sample types at different times. Pan-genome analysis identified a total of 21,909 genes, including 1,563 core genes, 8,292 strain-specific genes, and 12,054 accessory genes. Multiple putative virulence-associated genes were identified, encompassing 14 categories and 16 subcategories, with 171 distinct virulence factors. Three different plasmid replicon types were detected in 33 genomes. Eleven classes of antibiotic resistance genes and 352 integrons were identified. Antimicrobial susceptibility testing revealed a high resistance rate to imipenem and colistin among the strains studied, with 5 strains exhibiting multidrug resistance. However, they were all sensitive to amikacin, minocycline, and tigecycline. Our findings clarify the genomic characteristics and population structure of Shewanella algae in Hainan Province. The results offer insights into the genetic basis of pathogenicity in Shewanella algae and enhance our understanding of its global phylogeography.

Multi-omics analysis of excessive nitrogen fertilizer application: Assessing environmental damage and solutions in potato farming

Potatoes (Solanum tuberosum L.) are the third largest food crop globally and are pivotal for global food security. Widespread N fertilizer waste in potato cultivation has caused diverse environmental issues. This study employed microbial metagenomic sequencing to analyze the causes behind the declining N use efficiency (NUE) and escalating greenhouse gas emissions resulting from excessive N fertilizer application. Addressing N fertilizer inefficiency through breeding has emerged as a viable solution for mitigating overuse in potato cultivation. In this study, transcriptome and metabolome analyses were applied to identify N fertilizer-responsive genes. Metagenomic sequencing revealed that excessive N fertilizer application triggered alterations in the population dynamics of 11 major bacterial phyla, consequently affecting soil microbial functions, particularly N metabolism pathways and bacterial secretion systems. Notably, the enzyme levels associated with NO3- increased, and those associated with NO and N2O increased. Furthermore, excessive N fertilizer application enhanced soil virulence factors and increased potato susceptibility to diseases. Transcriptome and metabolome sequencing revealed significant impacts of excessive N fertilizer use on lipid and amino acid metabolism pathways. Weighted gene co‑expression network analysis (WGCNA) was adopted to identify two genes associated with N fertilizer response: PGSC0003DMG400021157 and PGSC0003DMG400009544.

Sanitizer Resistance and Persistence of Listeria monocytogenes Isolates in Tree Fruit Packing Facilities

The foodborne pathogen Listeria monocytogenes can persist in produce processing environments, which increases the risk for food contamination. Increased resistance to antimicrobials commonly used in cleaning and sanitizing procedures may contribute to L. monocytogenes' persistence in these environments. This study aimed to evaluate sanitizer resistance in L. monocytogenes isolates collected from three tree fruit packing facilities (F1, F2, and F3) during packing seasons 2020-2021 (Y1) and 2021-2022 (Y2), and to assess evidence of persistence based on the genomic similarity of isolates to historical isolates collected in previous years. L. monocytogenes isolates collected in 2020-2022 (n = 44) were tested for resistance to peroxyacetic acid (PAA) and a proprietary biofilm-removing agent using a broth microdilution assay. Further, L. monocytogenes isolates were whole genome sequenced and screened for the presence of antimicrobial resistance and virulence genes, as well as to assess the genomic similarity of isolates using the CFSAN SNP bioinformatic pipeline. Over half (57%) of the tested isolates had a PAA minimum inhibitory concentration (MIC) of 250 ppm, which was similar to the applied concentration of the PAA sanitizer in the three facilities (230 ppm). In contrast, 80% of tested isolates had a biofilm remover MIC of 0.13 ppm, which was substantially below the concentration applied in the facilities (137 ppm). Genomes of all tested isolates carried antimicrobial resistance (fosX, lin, mdrL, mprF, and norB) and virulence (inlA, inlB, plcA, plcB, prfA, hly, mpl, and iap) genes. L. monocytogenes isolates collected between 2020 and 2022 belonged to three distinct lineages, with 22 multilocus sequence types (MLSTs) belonging to 22 different clonal complexes. Genomic similarity analysis with historical isolates collected from the same facilities in 2016-2017 demonstrated a 5-year persistence of the genotypes ST 1003 and ST 554 in F2, which were no longer detected in 2022. Overall, our results highlight the need to re-evaluate sanitizer concentrations to effectively control persistent L. monocytogenes strains in tree fruit packing facilities.

Mechanisms of horizontal gene transfer and viral contribution to the fate of intracellular and extracellular antibiotic resistance genes in anaerobic digestion supplemented with conductive materials under ammonia stress

The addition of conductive materials (CMs) is an effective strategy for mitigating ammonia inhibition during anaerobic digestion (AD). However, the introduction of CMs can result in increased antibiotic resistance genes (ARGs) pollution, potentially facilitated by enhanced horizontal gene transfer (HGT). The complex dynamics of intracellular and extracellular ARGs (iARGs/eARGs) and the mechanisms underlying their transfer, mediated by CMs, in ammonia-stressed AD systems remain unclear. In this study, we investigated the effects of three commonly used CMs-nano magnetite (Mag), nano zero-valent iron (nZVI), and granular activated carbon (GAC)-on the fate of iARGs and eARGs during the AD of waste activated sludge under ammonia stress. The results revealed an unexpected enrichment of iARGs by 1.5 %-10.9 % and a reduction of eARGs by 14.1 %-25.2 % in CM-supplemented AD. This discrepancy in the dynamics of iARGs and eARGs may be attributed to changes in microbial hosts and the horizontal transfer of ARGs. Notably, CMs activated prophages within antibiotic-resistant bacteria (ARB) and their symbiotic partners involved in vitamin B12 provision, leading to the lysis of ARB and the subsequent release of eARGs for transformation. Additionally, the abundance of potentially mobile ARGs, which co-occurred with mobile genetic elements, increased by 56.6 %-134.5 % with CM addition, highlighting an enhanced potential for the HGT of ARGs. Specifically, Mag appeared to promote both transformation and conjugation processes, while nZVI only promoted conjugation. Moreover, none of the three CMs had any discernible impact on transduction. GAC proved superior to both nano Mag and nZVI in controlling the enrichment of iARGs, reducing eARGs, and limiting HGTs simultaneously. Overall, these findings provide novel insights into the role of viruses and the mechanisms of ARG spread in CM-assisted AD, offering valuable information for developing strategies to mitigate ARG pollution in practical applications.

Host genetics and microbiota data analysis in colorectal cancer research

Colorectal cancer (CRC) is a heterogeneous disease with a complex aetiology influenced by a myriad of genetic and environmental factors. Despite advances in CRC research, it is a major burden of disease, with the second highest incidence and third leading cause of cancer deaths worldwide. To individualise diagnosis, prognosis, and treatment of CRC, developing new strategies combining precision medicine and bioinformatic procedures is promising. Precision medicine is based on omics technologies and aims to individualise the management of CRC based on patient host genetic characteristics and microbiota. Bioinformatics is central to the application of personalised medicine because it enables the analysis of large datasets generated by these technologies. At the level of host genetics, bioinformatics allows the identification of mutations, genes, molecular pathways, biomarkers and drugs relevant to colorectal carcinogenesis. At the microbiota level, bioinformatics is fundamental to analysing microbial communities' composition and functionality and developing biomarkers and personalised microbiota-based therapies. This paper explores the host and microbiota genetic data analysis in CRC research.

Deep longitudinal lower respiratory tract microbiome profiling reveals genome-resolved functional and evolutionary dynamics in critical illness

The lower respiratory tract (LRT) microbiome impacts human health, especially among critically ill patients. However, comprehensive characterizations of the LRT microbiome remain challenging due to low microbial mass and host contamination. We develop a chelex100-based low-biomass microbial-enrichment method (CMEM) that enables deep metagenomic profiling of LRT samples to recover near-complete microbial genomes. We apply the method to 453 longitudinal LRT samples from 157 intensive care unit (ICU) patients in three geographically distant hospitals. We recover 120 high-quality metagenome-assembled genomes (MAGs) and associated plasmids without culturing. We detect divergent longitudinal microbiome dynamics and hospital-specific dominant opportunistic pathogens and resistomes in pneumonia patients. Diagnosed pneumonia and the ICU stay duration were associated with the abundance of specific antibiotic-resistance genes (ARGs). Moreover, CMEM can serve as a robust tool for genome-resolved analyses. MAG-based analyses reveal strain-specific resistome and virulome among opportunistic pathogen strains. Evolutionary analyses discover increased mobilome in prevailing opportunistic pathogens, highly conserved plasmids, and new recombination hotspots associated with conjugative elements and prophages. Integrative analysis with epidemiological data reveals frequent putative inter-patient strain transmissions in ICUs. In summary, we present a genome-resolved functional, transmission, and evolutionary landscape of the LRT microbiota in critically ill patients.

Colonisation by multidrug-resistant organisms in health workers in primary care: narrow spectrum oral antimicrobials are a risk factor

BACKGROUND

Limited information exists on carriage of multidrug-resistant organisms (MDRO) by health workers (HWs) in primary care settings. This study aims to determine the prevalence of MDRO carriage among HWs in primary care and to identify associated risk factors.

METHODS

A cross-sectional study was conducted across all 12 primary care units in São Caetano do Sul-SP, Brazil, from October to December 2023. Self-collected samples (nasal, oropharyngeal, and inguinal) were obtained. Environment cultures (potable water, sewage and stream water) were evaluated. Stenotrophomonas maltophilia isolates (human and environmental) were typed.

RESULTS

The study included 265/288 (92%) of HWs in primary care teams, mostly women with a median age of 47 years (IQR 38-57); 78% had no comorbidities. MDRO colonisation was found in 8.7% (23 HWs). The following bacteria were found: S. maltophilia (n = 9; 3.4%) in inguinal swabs; methicillin-resistant Staphylococcus aureus (n = 8; 3%) from all sites; extended-spectrum ß-lactamase-producing bacteria (n = 5; 2%) in inguinal swabs; and vancomycin-resistant enterococci in an inguinal swab (n = 1; 0.4%). Previous antibiotic use was significantly associated with MDRO colonisation (OR 2.91, 95% CI 1.19-7.09, p = 0.018), mainly narrow spectrum oral beta-lactams and macrolides. S. malthophilia was polyclonal and human and environmental isolates differed.

CONCLUSION

Colonisation by MRSA, VRE, and ESBL-producing bacteria was low; however, 4% were surprisingly colonized by polyclonal S. maltophilia. This pathogen may also suggest using narrow-spectrum rather than the expected broad-spectrum antimicrobials. Antibiotic use was the only risk factor found, mainly with oral narrow-spectrum drugs.

Core and accessory genomic traits of Vibrio cholerae O1 drive lineage transmission and disease severity

In Bangladesh, Vibrio cholerae lineages are undergoing genomic evolution, with increased virulence and spreading ability. However, our understanding of the genomic determinants influencing lineage transmission and disease severity remains incomplete. Here, we developed a computational framework using machine-learning, genome scale metabolic modelling (GSSM) and 3D structural analysis, to identify V. cholerae genomic traits linked to lineage transmission and disease severity. We analysed in-patients isolates from six Bangladeshi regions (2015-2021), and uncovered accessory genes and core SNPs unique to the most recent dominant lineage, with virulence, motility and bacteriophage resistance functions. We also found a strong correlation between V. cholerae genomic traits and disease severity, with some traits overlapping those driving lineage transmission. GSMM and 3D structure analysis unveiled a complex interplay between transcription regulation, protein interaction and stability, and metabolic networks, associated to lifestyle adaptation, intestinal colonization, acid tolerance and symptom severity. Our findings support advancing therapeutics and targeted interventions to mitigate cholera spread.

Sequence Types and Antimicrobial Resistance Profiles of Salmonella Typhimurium in the Food Chain in Singapore

Salmonella remains a significant foodborne pathogen globally with S. Typhimurium presenting as a frequently occurring serovar. This study aimed to characterize 67 S. Typhimurium isolates from humans, food, farms, and slaughterhouses collected in Singapore from 2016 to 2017. Using whole-genome sequencing analysis, the isolates were found to belong to either ST19 (n = 33) or ST36 (n = 34). ST36 predominated in human intestinal and chicken isolates, while human extra-intestinal and non-chicken food isolates belonged to ST19. Plasmids were predicted in 88.1% (n = 59) of the isolates with the most common incompatibility group profiles being IncFIB(S), IncFII(S) and IncQ1. IncFIB(S) (adjusted p-value < 0.05) and IncFII(S) (adjusted p-value < 0.05) were significantly more prevalent in ST19 isolates, while Col156 (adjusted p-value < 0.05) was more significantly found in ST36 isolates. ST36 isolates exhibited higher resistance to multiple antibiotic classes such as penicillins, phenicols, folate pathway inhibitors, aminoglycosides, β-lactam/β-lactamase inhibitor combinations, tetracyclines, and fluoroquinolones. Phylogenetics analysis suggested potential shared routes of transmission among human, chicken, farm and slaughterhouse environments. Taken together, this study offers a cross-sectional epidemiological insight into the genomic epidemiology and antimicrobial landscape of S. Typhimurium isolates in Singapore, informing strategies for future public health and food safety surveillance.

Phenotypic and genomic analysis of the hypervirulent methicillin-resistant Staphylococcus aureus ST630 clone in China

Methicillin-resistant Staphylococcus aureus (MRSA) sequence type 630 (ST630) is a rarely reported lineage worldwide. This study aimed to trace the dissemination of the emerging MRSA ST630 clones in China and investigate their virulence potential. We collected 22 ST630-MRSA isolates from across China and performed whole-genome sequencing analysis and virulence characterization on these isolates. Epidemiological results showed that MRSA ST630 isolates were primarily isolated from pus/wound secretions, mainly originating from Jiangxi province, and carried diverse virulence and drug resistance genes. Staphylococcal cassette chromosome mec type V (SCCmec V) predominated (11/22, 50.0%) among the MRSA ST630 isolates. Interestingly, nearly half (45.5%) of the 22 ST630-MRSA isolates tested lacked intact SCCmec elements. Phylogenetic analysis demonstrated that ST630-MRSA could be divided into two distinct clades, with widespread dissemination mainly in Chinese regions. Five representative isolates were selected for phenotypic assays, including hemolysin activity, real-time fluorescence quantitative PCR, western blot analysis, hydrogen peroxide killing assay, blood killing assay, cell adhesion and invasion assay, and mouse skin abscess model. The results showed that, compared to the USA300-LAC strain, ST630 isolates exhibited particularly strong invasiveness and virulence in the aforementioned phenotypic assays. This study described the emergence of a highly virulent ST630-MRSA lineage and improved our insight into the molecular epidemiology of ST630 clones in China.IMPORTANCEMethicillin-resistant Staphylococcus aureus (MRSA) sequence type 630 (ST630) is an emerging clone with an increasing isolation rate in China. This study raises awareness of the hypervirulent MRSA ST630 clones in China and alerts people to their widespread dissemination. ST630-staphylococcal cassette chromosome mec V is a noteworthy clone in China, and we present the first comprehensive genetic and phenotypic analysis of this lineage. Our findings provide valuable insights for the prevention and control of infections caused by this emerging MRSA clone.

Gut microbiota-bile acid crosstalk contributes to intestinal damage after nitrate exposure in Bufo gargarizans tadpoles

Bile acids (BAs) are amphipathic steroid acids whose production and diversity depend on both host and microbial metabolism. Nitrate (NO3-) is a widespread pollutant in aquatic ecosystems, which can cause rapid changes in microbial community structure and function. However, the effect of gut microbiota reshaped by nitrate‑nitrogen (NO3-N) on BAs profiles remains unclarified. To test this, intestinal targeted BAs metabolomics and fecal metagenomic sequencing were performed on Bufo gargarizans tadpoles treated with different concentrations of NO3-N. NO3-N exposure induced a reduction in the abundance of microbiota with bile acid-inducible enzymes (BAIs) and/or hydroxysteroid dehydrogenases (HSDHs), thus inhibiting the conversion of primary BAs to secondary BAs. Inhibition of BAs biotransformation decreased protective hydrophilic BAs (UDCA) and increased toxic hydrophobic BAs (CA and CDCA), which may contribute to intestinal histopathological damage. Moreover, we found that NO3-N treatment increased microbial virulence factors and decreased Glycoside hydrolases, further highlighting the deleterious risk of NO3-N. Overall, this study shed light on the complex interactions of NO3-N, gut microbiota, and BAs, and emphasized the hazardous effects of NO3-N pollution on the health of amphibians.

Faecal microbiota and cytokine profiles of rural Cambodian infants linked to diet and diarrhoeal episodes

The gut microbiota of infants in low- to middle-income countries is underrepresented in microbiome research. This study explored the faecal microbiota composition and faecal cytokine profiles in a cohort of infants in a rural province of Cambodia and investigated the impact of sample storage conditions and infant environment on microbiota composition. Faecal samples collected at three time points from 32 infants were analysed for microbiota composition using 16S rRNA amplicon sequencing and concentrations of faecal cytokines. Faecal bacterial isolates were subjected to whole genome sequencing and genomic analysis. We compared the effects of two sample collection methods due to the challenges of faecal sample collection in a rural location. Storage of faecal samples in a DNA preservation solution preserved Bacteroides abundance. Microbiota analysis of preserved samples showed that Bifidobacterium was the most abundant genus with Bifidobacterium longum the most abundant species, with higher abundance in breast-fed infants. Most infants had detectable pathogenic taxa, with Shigella and Klebsiella more abundant in infants with recent diarrhoeal illness. Neither antibiotics nor infant growth were associated with gut microbiota composition. Genomic analysis of isolates showed gene clusters encoding the ability to digest human milk oligosaccharides in B. longum and B. breve isolates. Antibiotic-resistant genes were present in both potentially pathogenic species and in Bifidobacterium. Faecal concentrations of Interlukin-1alpha and vascular endothelial growth factor were higher in breast-fed infants. This study provides insights into an underrepresented population of rural Cambodian infants, showing pathogen exposure and breastfeeding impact gut microbiota composition and faecal immune profiles.

In Vitro and In Vivo Studies on a Mononuclear Ruthenium Complex Reveals It is a Highly Effective, Fast-Acting, Broad-Spectrum Antimicrobial in Physiologically Relevant Conditions

The crystal structure of a previously reported antimicrobial RuII complex that targets bacterial DNA is presented. Studies utilizing clinical isolates of Gram-negative bacteria that cause catheter-associated urinary tract infection, (CA)UTI, in media that model urine and plasma reveal that good antimicrobial activity is maintained in all conditions tested. Experiments with a series of Staphylococcus aureus clinical isolates show that, unlike the majority of previously reported RuII-based antimicrobial leads, the compound retains its potent activity even in MRSA strains. Furthermore, experiments using bacteria in early exponential growth and at different pHs reveal that the compound also retains its activity across a range of conditions that are relevant to those encountered in clinical settings. Combinatorial studies involving cotreatment with conventional antibiotics or a previously reported analogous dinuclear RuII complex showed no antagonistic effects. In fact, although all combinations show distinct additive antibacterial activity, in one case, this effect approaches synergy. It was found that the Galleria Mellonella model organism infected with a multidrug resistant strain of the ESKAPE pathogen Acinetobacter baumannii could be successfully treated and totally cleared within 48 h after a single dose of the lead complex with no detectable deleterious effect to the host.

zol & fai: large-scale targeted detection and evolutionary investigation of gene clusters

Many universally and conditionally important genes are genomically aggregated within clusters. Here, we introduce fai and zol, which together enable large-scale comparative analysis of different types of gene clusters and mobile-genetic elements (MGEs), such as biosynthetic gene clusters (BGCs) or viruses. Fundamentally, they overcome a current bottleneck to reliably perform comprehensive orthology inference at large scale across broad taxonomic contexts and thousands of genomes. First, fai allows the identification of orthologous instances of a query gene cluster of interest amongst a database of target genomes. Subsequently, zol enables reliable, context-specific inference of ortholog groups for individual protein-encoding genes across gene cluster instances. In addition, zol performs functional annotation and computes a variety of evolutionary statistics for each inferred ortholog group. Importantly, in comparison to tools for visual exploration of homologous relationships between gene clusters, zol can scale to thousands of gene cluster instances and produce detailed reports that are easy to digest. To showcase fai and zol, we apply them for: (i) longitudinal tracking of a virus in metagenomes, (ii) discovering novel population-level genetic insights of two common BGCs in the fungal species Aspergillus flavus, and (iii) uncovering large-scale evolutionary trends of a virulence-associated gene cluster across thousands of genomes from a diverse bacterial genus.

Preventive and Therapeutic Effects of Lactiplantibacillus plantarum HD02 and MD159 through Mast Cell Degranulation Inhibition in Mouse Models of Atopic Dermatitis

As the relationship between the gut microbiome and allergies becomes better understood, targeted strategies to prevent and treat allergies through gut microbiome modulation are being increasingly developed. In the study presented herein, we screened various probiotics for their ability to inhibit mast cell degranulation and identified Lactiplatibacillus plantarum HD02 and MD159 as effective candidates. The two strains significantly attenuated vascular permeability induced by mast cell degranulation in a passive cutaneous anaphylaxis (PCA) model and, in the MC903-induced murine atopic dermatitis (AD) model, demonstrated comparable preventive effects against allergies, reducing blood levels of MCPT-1 (mast cell protease-1) and total IgE. In the house dust mite (HDM)-induced murine AD model, both L. plantarum HD02 and MD159 showed therapeutic effects, with L. plantarum HD02 demonstrating superior efficacy. Nevertheless, L. plantarum MD159 better suppressed transepidermal water loss (TEWL). Furthermore, L. plantarum HD02 and MD159 significantly increased the number of splenic Foxp3+ regulatory T cells, with L. plantarum MD159 having a more pronounced effect. However, only L. plantarum HD02 achieved a reduction in immune cells in the draining lymph nodes. Our findings highlight L. plantarum HD02 and MD159 as promising candidates for the prevention and treatment of allergies, demonstrating significant efficacy in suppressing mast cell degranulation, reducing the number of allergy biomarkers, and modulating immune responses in experimental models of AD. Their distinct mechanisms of action suggest potential complementary roles in addressing allergic diseases, underscoring their therapeutic promise in clinical applications.

DTVF: A User-Friendly Tool for Virulence Factor Prediction Based on ProtT5 and Deep Transfer Learning Models

Virulencefactors (VFs) are key molecules that enable pathogens to evade the immune systems of the host. These factors are crucial for revealing the pathogenic processes of microbes and drug discovery. Identification of virulence factors in microbes become an important problem in the field of bioinformatics. To address this problem, this study proposes a novel model DTVF (Deep Transfer Learning for Virulence Factor Prediction), which integrates the ProtT5 protein sequence extraction model with a dual-channel deep learning model. In the dual-channel deep learning model, we innovatively integrate long short-term memory (LSTM) with convolutional neural networks (CNNs), creating a novel integrated architecture. Furthermore, by incorporating the attention mechanism, the accuracy of VF detection was significantly enhanced. We evaluated the DTVF model against other excellent-performing models in the field. DTVF demonstrates superior performance, achieving an accuracy rate of 84.55% and an AUROC of 92.08% on the benchmark dataset. DTVF shows state-of-the-art performance in this field, surpassing the existing models in nearly all metrics. To facilitate the use of biologists, we have also developed an interactive web-based user interface version of DTVF based on Gradio.

Unveiling potential roles of earthworms in mitigating the presence of virulence factor genes in terrestrial ecosystems

Earthworms can redistribute soil microbiota, and thus might affect the profile of virulence factor genes (VFGs) which are carried by pathogens in soils. Nevertheless, the knowledge of VFG profile in the earthworm guts and its interaction with earthworm gut microbiome is still lacking. Herein, we characterized earthworm gut and soil microbiome and VFG profiles in natural and agricultural ecosystems at a national scale using metagenomics. VFG profiles in the earthworm guts significantly differed from those in the surrounding soils, which was mainly driven by variations of bacterial communities. Furthermore, the total abundance of different types of VFGs in the earthworm guts was about 20-fold lower than that in the soils due to the dramatic decline (also by approximately 20-fold) of VFG-carrying bacterial pathogens in the earthworm guts. Additionally, five VFGs related to nutritional/metabolic factors and stress survival were identified as keystones merely in the microbe-VFG network in the earthworm guts, implying their pivotal roles in facilitating pathogen colonization in earthworm gut microhabitats. These findings suggest the potential roles of earthworms in reducing risks related to the presence of VFGs in soils, providing novel insights into earthworm-based bioremediation of VFG contamination in terrestrial ecosystems.