Enabling genomic island prediction and comparison in multiple genomes to investigate bacterial evolution and outbreaks

Association of pathogenic determinants of Fusobacterium necrophorum with bacteremia, and Lemierre's syndrome

Fusobacterium necrophorum is a Gram-negative anaerobic bacterium responsible for localized infections of the oropharynx that can evolve into bacteremia and/or septic thrombophlebitis of the jugular vein or peritonsillar vein, called Lemierre's syndrome. To identify microbial genetic determinants associated with the severity of this life-threatening disease, 70 F. necrophorum strains were collected and grouped into two categories according to the clinical presentation: (i) localized infection, (ii) bacteremia with/without Lemierre's syndrome. Comparative genomic analyses revealed two clades with distinct genetic content, one clade being significantly enriched with isolates from subjects with bacteremia. To identify genetic determinants contributing to F. necrophorum pathogenicity, genomic islands and virulence factor orthogroups (OVFs) were predicted. The presence/absence profiles of OVFs did not group isolates according to their clinical category, but rather according to their phylogeny. However, a variant of lktA, a key virulence factor, with a frameshift deletion that results in two open reading frames, was associated with bacteremia. Moreover, a genome-wide association study identified three orthogroups associated with bacteremic strains: (i) cas8a1, (ii) a sodium/solute symporter, and (iii) a POP1 domain-containing protein. Further studies must be performed to assess the functional impact of lktA mutation and of these orthogroups on the physiopathological mechanisms of F. necrophorum infections.

Genetic diversity, determinants, and dissemination of Burkholderia pseudomallei lineages implicated in melioidosis in Northeast Thailand

Melioidosis is an often-fatal neglected tropical disease caused by an environmental bacterium Burkholderia pseudomallei. However, our understanding of the disease-causing bacterial lineages, their dissemination, and adaptive mechanisms remains limited. To address this, we conduct a comprehensive genomic analysis of 1,391 B. pseudomallei isolates collected from nine hospitals in northeast Thailand between 2015 and 2018, and contemporaneous isolates from neighbouring countries, representing the most densely sampled collection to date. Our study identifies three dominant lineages, each with unique gene sets potentially enhancing bacterial fitness in the environment. We find that recombination drives lineage-specific gene flow. Transcriptome analyses of representative clinical isolates from each dominant lineage reveal increased expression of lineage-specific genes under environmental conditions in two out of three lineages. This underscores the potential importance of environmental persistence for these dominant lineages. The study also highlights the influence of environmental factors such as terrain slope, altitude, and river direction on the geographical dispersal of B. pseudomallei. Collectively, our findings suggest that environmental persistence may play a role in facilitating the spread of B. pseudomallei, and as a prerequisite for exposure and infection, thereby providing useful insights for informing melioidosis prevention and control strategies.

Comprehensive genome analysis of five novel flavobacteria: Flavobacterium piscisymbiosum sp. nov., Flavobacterium pisciphilum sp. nov., Flavobacterium flavipigmentatum sp. nov., Flavobacterium lipolyticum sp. nov. and Flavobacterium cupriresistens sp. nov

Eight isolates were obtained through a study on culture-dependent bacteria from fish farms and identified as members of the genus Flavobacterium based on pairwise analysis of the 16S rRNA gene sequences. The highest pairwise identity values were calculated as 98.8 % for strain F-30 T and Flavobacterium bizetiae, 99.0 % for strain F-65 T and Flavobacterium branchiarum, 98.7 % for strain F-126 T and Flavobacterium tructae, 98.2 % for strain F-323 T and Flavobacterium cupreum while 99.7 % identity level was detected for strain F-70 T and Flavobacterium geliluteum. In addition, strains F-33, Fl-77, and F-70 T shared 100 % identical 16S rRNA genes, while strains F-323 T and Fl-318 showed 99.9 % identity. A polyphasic approach including comparative analysis of whole-genome data was employed to ascertain the taxonomic provenance of the strains. In addition to the morphological, physiological, biochemical and chemotaxonomic characteristics of the strains, the overall genome-relatedness indices of dDDH and ANI below the established thresholds confirmed the classification of the strains as five novel species within the genus Flavobacterium. The comprehensive genome analyses of the strains were also conducted to determine the biosynthetic gene clusters, virulence features and ecological distribution patterns. Based on the polyphasic characterisations, including comparative genome analyses, it is concluded that strains F-30 T, F-65 T, F-70 T, F126T and F-323 T represent five novel species within the genus Flavobacterium for which Flavobacterium piscisymbiosum sp. nov. F-30 T (=JCM 34194 T = KCTC 82254 T), Flavobacterium pisciphilum sp. nov. F-65 T (=JCM 34197 T = KCTC 82257 T), Flavobacterium flavipigmentatum sp. nov. F-70 T (Fl-33 = Fl-77 = JCM 34198 T = KCTC 82258 T), Flavobacterium lipolyticum sp. nov. F-126 T (JCM 34199 T = KCTC 82259 T) and Flavobacterium cupriresistens sp. nov. F-323 T (Fl-318 = JCM 34200 T = KCTC 82260 T), are proposed.

In Silico Analysis of the Phylogenetic and Physiological Characteristics of Sphingobium indicum B90A: A Hexachlorocyclohexane-Degrading Bacterium

The study focuses on the in silico genomic characterization of Sphingobium indicum B90A, revealing a wealth of genes involved in stress response, carbon monoxide oxidation, β-carotene biosynthesis, heavy metal resistance, and aromatic compound degradation, suggesting its potential as a bioremediation agent. Furthermore, genomic adaptations among nine Sphingomonad strains were explored, highlighting shared core genes via pangenome analysis, including those related to the shikimate pathway and heavy metal resistance. The majority of genes associated with aromatic compound degradation, heavy metal resistance, and stress response were found within genomic islands across all strains. Sphingobium indicum UT26S exhibited the highest number of genomic islands, while Sphingopyxis alaskensis RB2256 had the maximum fraction of its genome covered by genomic islands. The distribution of lin genes varied among the strains, indicating diverse genetic responses to environmental pressures. Additionally, in silico evidence of horizontal gene transfer (HGT) between plasmids pSRL3 and pISP3 of the Sphingobium and Sphingomonas genera, respectively, has been provided. The manuscript offers novel insights into strain B90A, highlighting its role in horizontal gene transfer and refining evolutionary relationships among Sphingomonad strains. The discovery of stress response genes and the czcABCD operon emphasizes the potential of Sphingomonads in consortia development, supported by genomic island analysis.

The mobilome of Staphylococcus aureus from wild ungulates reveals epidemiological links at the animal-human interface

Staphylococcus aureus thrives at animal-human-environment interfaces. A large-scale work from our group indicated that antimicrobial resistance (AMR) in commensal S. aureus strains from wild ungulates is associated with agricultural land cover and livestock farming, raising the hypothesis that AMR genes in wildlife strains may originate from different hosts, namely via exchange of mobile genetic elements (MGE). In this work, we generate the largest available dataset of S. aureus draft genomes from wild ungulates in Portugal and explore their mobilome, which can determine important traits such as AMR, virulence, and host specificity, to understand MGE exchange. Core genome multi-locus sequence typing based on 98 newly generated draft genomes and 101 publicly available genomes from Portugal demonstrated that the genomic relatedness of S. aureus from wild ungulates assigned to livestock-associated sequence types (ST) is greater compared to wild ungulate isolates assigned to human-associated STs. Screening of host specificity determinants disclosed the unexpected presence in wildlife of the immune evasion cluster encoded in φSa3 prophage, described as a human-specific virulence determinant. Additionally, two plasmids, pAVX and pETB, previously associated with avian species and humans, respectively, and the Tn553 transposon were detected. Both pETB and Tn553 encode penicillin resistance through blaZ. Pangenome analysis of wild ungulate isolates shows a core genome fraction of 2133 genes, with isolates assigned to ST72 and ST3224 being distinguished from the remaining by MGEs, although there is no reported role of these in adaptation to wildlife. AMR related gene clusters found in the shell genome are directly linked to resistance against penicillin, macrolides, fosfomycin, and aminoglycosides, and they represent mobile ARGs. Altogether, our findings support epidemiological interactions of human and non-human hosts at interfaces, with MGE exchange, including AMR determinants, associated with putative indirect movements of S. aureus among human and wildlife hosts that might be bridged by livestock.

Exploring the comparative genome of rice pathogen Burkholderia plantarii: unveiling virulence, fitness traits, and a potential type III secretion system effector

This study presents a comprehensive genomic analysis of Burkholderia plantarii, a rice pathogen that causes blight and grain rot in seedlings. The entire genome of B. plantarii KACC 18964 was sequenced, followed by a comparative genomic analysis with other available genomes to gain insights into its virulence, fitness, and interactions with rice. Multiple secondary metabolite gene clusters were identified. Among these, 12 demonstrated varying similarity levels to known clusters linked to bioactive compounds, whereas eight exhibited no similarity, indicating B. plantarii as a source of potentially novel secondary metabolites. Notably, the genes responsible for tropolone and quorum sensing were conserved across the examined genomes. Additionally, B. plantarii was observed to possess three complete CRISPR systems and a range of secretion systems, exhibiting minor variations among the analyzed genomes. Genomic islands were analyzed across the four genomes, and a detailed study of the B. plantarii KACC 18964 genome revealed 59 unique islands. These islands were thoroughly investigated for their gene contents and potential roles in virulence. Particular attention has been devoted to the Type III secretion system (T3SS), a crucial virulence factor. An in silico analysis of potential T3SS effectors identified a conserved gene, aroA. Further mutational studies, in planta and in vitro analyses validated the association between aroA and virulence in rice. Overall, this study enriches our understanding of the genomic basis of B. plantarii pathogenicity and emphasizes the potential role of aroA in virulence. This understanding may guide the development of effective disease management strategies.

Increased environmental microbial diversity reduces the disease risk of a mosquitocidal pathogen

IMPORTANCE

The host-specific microbiotas of animals can both reduce and increase disease risks from pathogens. In contrast, how environmental microbial communities affect pathogens is largely unexplored. Aquatic habitats are of interest because water enables environmental microbes to readily interact with animal pathogens. Here, we focused on mosquitoes, which are important disease vectors as terrestrial adults but are strictly aquatic as larvae. We identified a pathogen of mosquito larvae from the field as a strain of Chromobacterium haemolyticum. Comparative genomic analyses and functional assays indicate this strain and other Chromobacterium are mosquitocidal but are also opportunistic pathogens of other animals. We also identify a critical role for diversity of the environmental microbiota in disease risk. Our study characterizes both the virulence mechanisms of a pathogen and the role of the environmental microbiota in disease risk to an aquatic animal of significant importance to human health.

Genetic diversity, determinants, and dissemination of Burkholderia pseudomallei lineages implicated in melioidosis in northeast Thailand

Melioidosis is an often-fatal neglected tropical disease caused by an environmental bacterium Burkholderia pseudomallei. However, our understanding of the disease-causing bacterial lineages, their dissemination, and adaptive mechanisms remains limited. To address this, we conducted a comprehensive genomic analysis of 1,391 B. pseudomallei isolates collected from nine hospitals in northeast Thailand between 2015 and 2018, and contemporaneous isolates from neighbouring countries, representing the most densely sampled collection to date. Our study identified three dominant lineages with unique gene sets enhancing bacterial fitness, indicating lineage-specific adaptation strategies. Crucially, recombination was found to drive lineage-specific gene flow. Transcriptome analyses of representative clinical isolates from each dominant lineage revealed heightened expression of lineage-specific genes in environmental versus infection conditions, notably under nutrient depletion, highlighting environmental persistence as a key factor in the success of dominant lineages. The study also revealed the role of environmental factors - slope of terrain, altitude, direction of rivers, and the northeast monsoons - in shaping B. pseudomallei geographical dispersal. Collectively, our findings highlight persistence in the environment as a pivotal element facilitating B. pseudomallei spread, and as a prelude to exposure and infection, thereby providing useful insights for informing melioidosis prevention and control strategies.

Characterization and genome analysis of Neobacillus mesonae NS-6, a ureolysis-driven strain inducing calcium carbonate precipitation

In this study, a highly promising bacterium was isolated from sandstone oil in the Ordos Basin, named strain NS-6 which exhibited exceptional urease production ability and demonstrated superior efficiency in inducing the deposition of calcium carbonate (CaCO3). Through morphological and physiochemical characteristics analysis, as well as 16S rRNA sequencing, strain NS-6 was identified as Neobacillus mesonae. The activity of urease and the formation of CaCO3 increased over time, reaching a maximum of 7.9 mmol/L/min and 184 mg (4.60 mg/mL) respectively at 32 h of incubation. Scanning Electron Microscopy (SEM) revealed CaCO3 crystals ranging in size from 5 to 6 μm, and Energy Dispersive X-ray (EDX) analysis verified the presence of calcium, carbon, and oxygen within the crystals. X-ray Diffraction (XRD) analysis further confirmed the composition of these CaCO3 crystals as calcite and vaterite. Furthermore, the maximum deposition of CaCO3 by strain NS-6 was achieved using response surface methodology (RSM), amounting to 193.8 mg (4.845 mg/mL) when the concentration of calcium ions was 0.5 mmol/L supplemented with 0.9 mmol/L of urea at pH 8.0. Genome-wide analysis revealed that strain NS-6 possesses a chromosome of 5,736,360 base pairs, containing 5,442 predicted genes, including 3,966 predicted functional genes and 1,476 functionally unknown genes. Genes like ureA, ureB, and ureC related to urea catabolism were identified by gene annotation, indicating that strain NS-6 is a typical urease-producing bacterium and possesses a serial of genes involved in metabolic pathways that mediated the deposition of CaCO3 at genetic level.

Staphylococcus brunensis sp. nov. isolated from human clinical specimens with a staphylococcal cassette chromosome-related genomic island outside of the rlmH gene bearing the ccrDE recombinase gene complex

Novel species of coagulase-negative staphylococci, which could serve as reservoirs of virulence and antimicrobial resistance factors for opportunistic pathogens from the genus Staphylococcus, are recognized in human and animal specimens due to advances in diagnostic techniques. Here, we used whole-genome sequencing, extensive biotyping, MALDI-TOF mass spectrometry, and chemotaxonomy to characterize five coagulase-negative strains from the Staphylococcus haemolyticus phylogenetic clade obtained from human ear swabs, wounds, and bile. Based on the results of polyphasic taxonomy, we propose the species Staphylococcus brunensis sp. nov. (type strain NRL/St 16/872T = CCM 9024T = LMG 31872T = DSM 111349T). The genomic analysis revealed numerous variable genomic elements, including staphylococcal cassette chromosome (SCC), prophages, plasmids, and a unique 18.8 kb-long genomic island SbCIccrDE integrated into the ribosomal protein L7 serine acetyltransferase gene rimL. SbCIccrDE has a cassette chromosome recombinase (ccr) gene complex with a typical structure found in SCCs. Based on nucleotide and amino acid identity to other known ccr genes and the distinct integration site that differs from the canonical methyltransferase gene rlmH exploited by SCCs, we classified the ccr genes as novel variants, ccrDE. The comparative genomic analysis of SbCIccrDE with related islands shows that they can accumulate virulence and antimicrobial resistance factors creating novel resistance elements, which reflects the evolution of SCC. The spread of these resistance islands into established pathogens such as Staphylococcus aureus would pose a great threat to the healthcare system. IMPORTANCE The coagulase-negative staphylococci are important opportunistic human pathogens, which cause bloodstream and foreign body infections, mainly in immunocompromised patients. The mobile elements, primarily the staphylococcal cassette chromosome mec, which confers resistance to methicillin, are the key to the successful dissemination of staphylococci into healthcare and community settings. Here, we present a novel species of the Staphylococcus genus isolated from human clinical material. The detailed analysis of its genome revealed a previously undescribed genomic island, which is closely related to the staphylococcal cassette chromosome and has the potential to accumulate and spread virulence and resistance determinants. The island harbors a set of conserved genes required for its mobilization, which we recognized as novel cassette chromosome recombinase genes ccrDE. Similar islands were revealed not only in the genomes of coagulase-negative staphylococci but also in S. aureus. The comparative genomic study contributes substantially to the understanding of the evolution and pathogenesis of staphylococci.

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 or homologous instances of a query gene cluster of interest amongst a database of target genomes. Subsequently, zol enables reliable, context-specific inference of protein-encoding ortholog groups for individual genes across gene cluster instances. In addition, zol performs functional annotation and computes a variety of statistics for each inferred ortholog group. These programs are showcased through application to: (i) longitudinal tracking of a virus in metagenomes, (ii) discovering novel population-genetic insights of two common BGCs in a fungal species, and (iii) uncovering large-scale evolutionary trends of a virulence-associated gene cluster across thousands of genomes from a diverse bacterial genus.

Genomic insight into Myroides oncorhynchi sp. nov., a new member of the Myroides genus, isolated from the internal organ of rainbow trout (Oncorhynchus mykiss)

The strain M-43^T was isolated from the Oncorhynchus mykiss from a fish farm in Mugla, Turkey. Pairwise 16S rRNA gene sequence analysis was used to identify strain M-43^T. The strain was a member of the genus Myroides sharing the highest 16S rRNA gene sequence identity levels of 98.7%, 98.3%, and 98.3% with the type strains of M. profundi D25^T, M. odoratimimus subsp. odoratimimus CCUG 39352^T and M. odoratimimus subsp. xuanwuensis DSM27251^T, respectively. A polyphasic taxonomic approach including whole genome-based analyses was employed to confirm the taxonomic provenance of strain M-43^T within the genus Myroides. The overall genome relatedness indices (OGRI) for strain M-43^T compared with its most closely related type strains M. odoratimimus subsp. xuanwuensis DSM 27251^T, M. profundi D25^T, and M. odoratimimus subsp. odoratimimus ATCC BAA-634^T, were calculated as 25.3%, 25.1%, and 25% for digital DNA-DNA hybridization (dDDH), 83.3%, 83.6%, and 83.4% for average nucleotide identity (ANI) analyses, respectively. The OGRI values between strain M-43^T and its close neighbors confirmed that the strain represents a novel species in the genus Myroides. The DNA G + C content of the strain is 33.7%. The major fatty acids are iso-C_15:0 and summed feature 9 (iso-C_17:1 ω9c and/or 10-methyl C_16:0). The predominant polar lipids are phosphatidylethanolamine, an amino-lipid and five unidentified lipids. The major respiratory quinone is MK-6. Chemotaxonomic and phylogenomic analyses of this isolate confirmed that the strain represents a novel species for which the name Myroides oncorhynchi sp. nov. is proposed, with M-43^T as the type strain (JCM 34205^T = KCTC 82265^T).

Genomic insight into Chryseobacterium turcicum sp. nov. and Chryseobacterium muglaense sp. nov. isolated from farmed rainbow trout in Turkey

Four strains, designated as C-2, C-17^T, C-39^T and Ch-15, were isolated from farmed rainbow trout samples showing clinical signs during an investigation for a fish-health screening study. The pairwise 16S rRNA gene sequence analysis showed that strain C-17^T shared the highest identity level of 98.1 % with the type strain of Chryseobacterium piscium LMG 23089^T while strains C-2, C-39^T and Ch-15 were closely related to Chryseobacterium balustinum DSM 16775^T with an identity level of 99.3 %. A polyphasic approach involving phenotypic, chemotaxonomic and genome-based analyses was employed to determine the taxonomic provenance of the strains. The overall genome relatedness indices including dDDH and ANI analyses confirmed that strains C-2, C-17^T, C-39^T and Ch-15 formed two novel species within the genus Chryseobacterium. Chemotaxonomic analyses showed that strains C-17^T and C-39^T have typical characteristics of the genus Chryseobacterium by having phosphatidylethanolamine in their polar lipid profile, MK-6 as only isoprenoid quinone and the presence of iso-C_15:0 as major fatty acid. The genome size and G + C content of the strains ranged between 4.4 and 5.0 Mb and 33.5 - 33.6 %, respectively. Comprehensive genome analyses revealed that the strains have antimicrobial resistance genes, prophages and horizontally acquired genes in addition to secondary metabolite-coding gene clusters. In conclusion, based on the polyphasic analyses conducted on the present study, strains C-17^T and C-39^T are representatives of two novel species within the genus Chryseobacterium, for which the names Chryseobacterium turcicum sp. nov. and Chryseobacterium muglaense sp. nov. with the type strains C-17^T (=JCM 34190^T = KCTC 82250^T) and C-39^T (=JCM 34191^T = KCTC 822251^T), respectively, are proposed.

Antarctic Sphingomonas sp. So64.6b showed evolutive divergence within its genus, including new biosynthetic gene clusters

Introduction

The antibiotic crisis is a major human health problem. Bioprospecting screenings suggest that proteobacteria and other extremophile microorganisms have biosynthetic potential for the production novel antimicrobial compounds. An Antarctic Sphingomonas strain (So64.6b) previously showed interesting antibiotic activity and elicitation response, then a relationship between environmental adaptations and its biosynthetic potential was hypothesized. We aimed to determine the genomic characteristics in So64.6b strain related to evolutive traits for the adaptation to the Antarctic environment that could lead to its diversity of potentially novel antibiotic metabolites.

Methods

The complete genome sequence of the Antarctic strain was obtained and mined for Biosynthetic Gene Clusters (BGCs) and other unique genes related to adaptation to extreme environments. Comparative genome analysis based on multi-locus phylogenomics, BGC phylogeny, and pangenomics were conducted within the closest genus, aiming to determine the taxonomic affiliation and differential characteristics of the Antarctic strain.

Results and discussion

The Antarctic strain So64.6b showed a closest identity with Sphingomonas alpina, however containing a significant genomic difference of ortholog cluster related to degradation multiple pollutants. Strain So64.6b had a total of six BGC, which were predicted with low to no similarity with other reported clusters; three were associated with potential novel antibiotic compounds using ARTS tool. Phylogenetic and synteny analysis of a common BGC showed great diversity between Sphingomonas genus but grouping in clades according to similar isolation environments, suggesting an evolution of BGCs that could be linked to the specific ecosystems. Comparative genomic analysis also showed that Sphingomonas species isolated from extreme environments had the greatest number of predicted BGCs and a higher percentage of genetic content devoted to BGCs than the isolates from mesophilic environments. In addition, some extreme-exclusive clusters were found related to oxidative and thermal stress adaptations, while pangenome analysis showed unique resistance genes on the Antarctic strain included in genetic islands. Altogether, our results showed the unique genetic content on Antarctic strain Sphingomonas sp. So64.6, −a probable new species of this genetically divergent genus–, which could have potentially novel antibiotic compounds acquired to cope with Antarctic poly-extreme conditions.

Exploring the mobilome and resistome of Enterococcus faecium in a One Health context across two continents

Enterococcus faecium is a ubiquitous opportunistic pathogen that is exhibiting increasing levels of antimicrobial resistance (AMR). Many of the genes that confer resistance and pathogenic functions are localized on mobile genetic elements (MGEs), which facilitate their transfer between lineages. Here, features including resistance determinants, virulence factors and MGEs were profiled in a set of 1273 E. faecium genomes from two disparate geographic locations (in the UK and Canada) from a range of agricultural, clinical and associated habitats. Neither lineages of E. faecium, type A and B, nor MGEs are constrained by geographic proximity, but our results show evidence of a strong association of many profiled genes and MGEs with habitat. Many features were associated with a group of clinical and municipal wastewater genomes that are likely forming a new human-associated ecotype within type A. The evolutionary dynamics of E. faecium make it a highly versatile emerging pathogen, and its ability to acquire, transmit and lose features presents a high risk for the emergence of new pathogenic variants and novel resistance combinations. This study provides a workflow for MGE-centric surveillance of AMR in Enterococcus that can be adapted to other pathogens.