An evaluation of the species and subspecies of the genus Salmonella with whole genome sequence data: Proposal of type strains and epithets for novel S. enterica subspecies VII, VIII, IX, X and XI

Comprehensive blueprint of Salmonella genomic plasticity identifies hotspots for pathogenicity genes

Understanding the dynamic evolution of Salmonella is vital for effective bacterial infection management. This study explores the role of the flexible genome, organised in regions of genomic plasticity (RGP), in shaping the pathogenicity of Salmonella lineages. Through comprehensive genomic analysis of 12,244 Salmonella spp. genomes covering 2 species, 6 subspecies, and 46 serovars, we uncover distinct integration patterns of pathogenicity-related gene clusters into RGP, challenging traditional views of gene distribution. These RGP exhibit distinct preferences for specific genomic spots, and the presence or absence of such spots across Salmonella lineages profoundly shapes strain pathogenicity. RGP preferences are guided by conserved flanking genes surrounding integration spots, implicating their involvement in regulatory networks and functional synergies with integrated gene clusters. Additionally, we emphasise the multifaceted contributions of plasmids and prophages to the pathogenicity of diverse Salmonella lineages. Overall, this study provides a comprehensive blueprint of the pathogenicity potential of Salmonella. This unique insight identifies genomic spots in nonpathogenic lineages that hold the potential for harbouring pathogenicity genes, providing a foundation for predicting future adaptations and developing targeted strategies against emerging human pathogenic strains.

Detection and Phenotypic Antimicrobial Susceptibility of Salmonella enterica Serotypes in Dairy Cattle Farms in the Po Valley, Northern Italy

The presence of Salmonella spp. in dairy cattle farms poses a major risk to animal health and welfare. This study focused on Salmonella detection in dairy farms located in the Cremona and Mantua provinces (northern Italy) in samples collected and submitted to laboratories in 2021-2022. A total of 2710 samples from different sources, including calf carcasses/organs (n = 128), rectal swabs (n = 1937), feces (n = 390), bulk milk (n = 93), and overshoes/swabs (n = 127) for environmental sampling, were analyzed for the presence of Salmonella spp. and were included in the present study. Our results indicate that Salmonella was most commonly firstly identified from calf carcasses and organs (61.67%) and that the serotypes most frequently detected in dairies were S. Dublin (38.33%), S. Typhimurium (23.33%), and S. Typhimurium monophasic variant (14.17%). The most common pathological findings in calf carcasses were enteritis, hepatosplenomegaly, and pneumonia. The antimicrobial resistance pattern analyzed using the MIC assay of 51 Salmonella isolates revealed the presence of multi-resistant strains, which pose a major risk to public and animal health.

Notification of bacterial strains made available by the United Kingdom National Collection of Type Cultures in 2022

Here, we report on the one hundred and twenty-five bacterial strains made available by the National Collection of Type Cultures in 2022 alongside a commentary on the strains, their provenance and significance.

What's on a prophage: analysis of Salmonella spp. prophages identifies a diverse range of cargo with multiple virulence- and metabolism-associated functions

The gain of mobile elements, such as prophages, can introduce cargo to the recipient bacterium that could facilitate its persistence in or expansion to a new environment, such as a host. While previous studies have focused on identifying and characterizing the genetic diversity of prophages, analyses characterizing the cargo that prophages carry have not been extensively explored. We characterized prophage regions from 303 Salmonella spp. genomes (representing 254 unique serovars) to assess the distribution of prophages in diverse Salmonella. On average, prophages accounted for 3.7% (0.1%-8.8%) of the total genomic content of each isolate. Prophage regions annotated as Gifsy 1 and Salmon Fels 1 were the most commonly identified intact prophages, suggesting that they are common throughout the Salmonella genus. Among 21,687 total coding sequences (CDSs) from intact prophage regions in subsp. enterica genomes, 7.5% (median; range: 1.1%-47.6%) were categorized as having a function not related to prophage integration or phage structure, some of which could potentially provide a functional attribute to the host Salmonella cell. These predicted functions could be broadly categorized into CDSs involved in: (i) modification of cell surface structures (i.e., glycosyltransferases); (ii) modulation of host responses (e.g., SodC/SodA, SopE, ArtAB, and typhoid toxin); (iii) conferring resistance to heavy metals and antimicrobials; (iv) metabolism of carbohydrates, amino acids, and nucleotides; and (v) DNA replication, repair, and regulation. Overall, our systematic analysis of prophage cargo highlights a broader role for prophage cargo in influencing the metabolic, virulence, and resistance characteristics of Salmonella.

IMPORTANCE

Lysogenic bacteriophages (phages) can integrate their genome into a bacterial host's genome, potentially introducing genetic elements that can affect the fitness of the host bacterium. The functions of prophage-encoded genes are important to understand as these genes could be mobilized and transferred to a new host. Using a large genomic dataset representing >300 isolates from all known subspecies and species of Salmonella, our study contributes important new findings on the distribution of prophages and the types of cargo that diverse Salmonella prophages carry. We identified a number of coding sequences (CDSs) annotated as having cell surface-modifying attributes, suggesting that prophages may have played an important role in shaping Salmonella's diverse surface antigen repertoire. Furthermore, our characterization of prophages suggests that they play a broader role in facilitating the acquisition and transfer of CDSs associated with metabolism, DNA replication and repair, virulence factors, and to a lesser extent, antimicrobial resistance.

Genomic investigation of Salmonella enterica Serovar Welikade from a pediatric diarrhea case first time in Shanghai, China

BACKGROUND

Salmonella, an important foodborne pathogen, was estimated to be responsible for 95.1 million cases and 50,771 deaths worldwide. Sixteen serovars were responsible for approximately 80% of Salmonella infections in humans in China, and infections caused by a few uncommon serovars have been reported in recent years, though not with S. Welikade. This study reports the first clinical case caused by S. Welikade in China and places Chinese S. Welikade isolates in the context of global isolates via genomic analysis. For comparison, S. Welikade isolates were also screened in the Chinese Local Surveillance System for Salmonella (CLSSS). The minimum inhibitory concentrations (MICs) of 28 antimicrobial agents were determined using the broth microdilution method. The isolates were sequenced on an Illumina platform to identify antimicrobial resistance genes, virulence genes, and phylogenetic relationships.

RESULTS

The S. Welikade isolate (Sal097) was isolated from a two-year-old boy with acute gastroenteritis in 2021. Along with the other two isolates found in CLSSS, the three Chinese isolates were susceptible to all the examined antimicrobial agents, and their sequence types (STs) were ST5123 (n = 2) and ST3774 (n = 1). Single nucleotide polymorphism (SNP)-based phylogenetic analysis revealed that global S. Welikade strains can be divided into four groups, and these three Chinese isolates were assigned to B (n = 2; Sal097 and XXB1016) and C (n = 1; XXB700). In Group B, the two Chinese ST5123 isolates were closely clustered with three UK ST5123 isolates. In Group C, the Chinese isolate was closely related to the other 12 ST3774 isolates. The number of virulence genes in the S. Welikade isolates ranged from 59 to 152. The galF gene was only present in Group A, the pipB2 gene was only absent from Group A, the avrA gene was only absent from Group B, and the allB, sseK1, sspH2, STM0287, and tlde1 were found only within Group C and D isolates. There were 15 loci unique to the Sal097 isolate.

CONCLUSION

This study is the first to characterize and investigate clinical S. Welikade isolates in China. Responsible for a pediatric case of gastroenteritis in 2021, the clinical isolate harbored no antimicrobial resistance and belonged to phylogenetic Group B of global S. Welikade genomes.

Pangenomic analysis identifies correlations between Akkermansia species and subspecies and human health outcomes

Aim: Akkermansia are common members of the human gastrointestinal microbiota. The prevalence of these mucophilic bacteria, especially Akkermansia muciniphila (A. muciniphila), correlates with immunological and metabolic health. The genus Akkermansia in humans includes species with significantly larger genomes than A. muciniphila, leading us to postulate that this added genetic content may influence how they impact human metabolic and immunological health. Methods: We conducted a pangenomic analysis of 234 Akkermansia complete or near-complete genomes. We also used high-resolution species and subspecies assignments to reanalyze publicly available metagenomic datasets to determine if there are relationships between Akkermansia species and A. muciniphila clades with various disease outcomes. Results: Analysis of genome-wide average nucleotide identity, 16S rRNA gene identity, conservation of core Akkermansia genes, and analysis of the fatty acid composition of representative isolates support the partitioning of the genus Akkermansia into several species. In addition, A. muciniphila sensu stricto, the most prevalent Akkermansia species in humans, should be subdivided into two subspecies. For a pediatric cohort, we observed species-specific correlations between Akkermansia abundance with baseline obesity or after various interventions. For inflammatory bowel disease cohorts, we identified a decreased abundance of Akkermansia in patients with ulcerative colitis or Crohn's disease, which was species and subspecies-dependent. In patients undergoing immune checkpoint inhibitor therapies for non-small cell lung carcinoma, we observed a significant association between one A. muciniphila subspecies and survival outcomes. Conclusion: Our findings suggest that the prevalence of specific Akkermansia species and/or subspecies can be crucial in evaluating their association with human health, particularly in different disease contexts, and is an important consideration for their use as probiotics.

Pathology, microbiology, and genetic diversity associated with Erysipelothrix rhusiopathiae and novel Erysipelothrix spp. infections in southern sea otters (Enhydra lutris nereis)

Erysipelothrix spp., including E. rhusiopathiae, are zoonotic bacterial pathogens that can cause morbidity and mortality in mammals, fish, reptiles, birds, and humans. The southern sea otter (SSO; Enhydra lutris nereis) is a federally-listed threatened species for which infectious disease is a major cause of mortality. We estimated the frequency of detection of these opportunistic pathogens in dead SSOs, described pathology associated with Erysipelothrix infections in SSOs, characterized the genetic diversity and antimicrobial susceptibility of SSO isolates, and evaluated the virulence of two novel Erysipelothrix isolates from SSOs using an in vivo fish model. From 1998 to 2021 Erysipelothrix spp. were isolated from six of >500 necropsied SSOs. Erysipelothrix spp. were isolated in pure culture from three cases, while the other three were mixed cultures. Bacterial septicemia was a primary or contributing cause of death in five of the six cases. Other pathology observed included suppurative lymphadenopathy, fibrinosuppurative arteritis with thrombosis and infarction, bilateral uveitis and endophthalmitis, hypopyon, petechia and ecchymoses, mucosal infarction, and suppurative meningoencephalitis and ventriculitis. Short to long slender Gram-positive or Gram-variable bacterial rods were identified within lesions, alone or with other opportunistic bacteria. All six SSO isolates had the spaA genotype-four isolates clustered with spaA E. rhusiopathiae strains from various terrestrial and marine animal hosts. Two isolates did not cluster with any known Erysipelothrix spp.; whole genome sequencing revealed a novel Erysipelothrix species and a novel E. rhusiopathiae subspecies. We propose the names Erysipelothrix enhydrae sp. nov. and Erysipelothrix rhusiopathiae ohloneorum ssp. nov. respectively. The type strains are E. enhydrae UCD-4322-04 and E. rhusiopathiae ohloneorum UCD-4724-06, respectively. Experimental injection of tiger barbs (Puntigrus tetrazona) resulted in infection and mortality from the two novel Erysipelothrix spp. Antimicrobial susceptibility testing of Erysipelothrix isolates from SSOs shows similar susceptibility profiles to isolates from other terrestrial and aquatic animals. This is the first description of the pathology, microbial characteristics, and genetic diversity of Erysipelothrix isolates recovered from diseased SSOs. Methods presented here can facilitate case recognition, aid characterization of Erysipelothrix isolates, and illustrate assessment of virulence using fish models.

Bifidobacterium longum subsp. iuvenis subsp. nov., a novel subspecies isolated from the faeces of weaning infants

The species Bifidobacterium longum currently comprises four subspecies: B. longum subsp. longum , B. longum subsp. infantis , B. longum subsp. suis and B. longum subsp. suillum . Recently, several studies on B. longum suggested the presence of a separate clade containing four strains isolated from infants and one from rhesus macaque. These strains shared a phylogenetic similarity to B. longum subsp. suis DSM 20210T and B. longum subsp. suillum JCM1995T [average nucleotide identity (ANI) of 98.1 %) while showed an ANI of 96.5 % with both B. longum subsp. infantis and B. longum subsp. longum . The current work describes five novel additional B. longum strains isolated from Bangladeshi weaning infants and demonstrates their common phylogenetic origin with those of the previously proposed separated clade. Based on polyphasic taxonomic approach comprising loci multilocus sequence analysis and whole genome multilocus sequence typing, all ten examined strains have been confirmed as a distinct lineage within the species B. longum with B. longum subsp. suis and B. longum subsp. suillum as closest subspecies. Interestingly, these strains are present in weaning infants and primates as opposed to their closest relatives which have been typically isolated from pig and calves. These strains, similarly to B. longum subsp. infantis , show a common capacity to metabolize the human milk oligosaccharide 3-fucosyllactose. Moreover, they harbour a riboflavin synthesis operon, which differentiate them from their closest subspecies, B. longum subsp. suis and B. longum subsp. suillum . Based on the consistent results from genotypical, ecological and phenotypical analyses, a novel subspecies with the name Bifidobacterium longum subsp. iuvenis, with type strain NCC 5000T (=LMG 32752T=CCOS 2034T), is proposed.

Gene flow and species boundaries of the genus Salmonella

The genus Salmonella comprises two species, Salmonella bongori and Salmonella enterica, which are infectious to a wide variety of animal hosts. The diversity within S. enterica has been further partitioned into 6-10 subspecies based on such features as host range, geography, and most recently, genetic relatedness and phylogenetic affiliation. Although Salmonella pathogenicity is attributable to large numbers of acquired virulence factors, the extent of homologous exchange in the species at large is apparently constrained such that the species and subspecies form distinct clusters of strains. To explore the extent of gene flow within and among subspecies, and to ultimately define true biological species, we evaluated patterns of recombination in over 1,000 genomes currently assigned to the genus. Those Salmonella subspecies containing sufficient numbers of sequenced genomes to allow meaningful analysis-i.e., subsp. enterica and diarizonae-were found to be reproductively isolated from one another and from all other subspecies. Based on the configuration of genomic sequence divergence among subspecies, it is expected that each of the other Salmonella subspecies will also represent a biological species. Our findings argue against the application of prescribed nucleotide-identity thresholds to delineate bacterial species and contend that the Biological Species Concept should not be disregarded for bacteria, even those, like Salmonella, that demonstrate complex patterns of species and subspecies divergence. IMPORTANCE The Biological Species Concept (BSC), which defines species boundaries based on the capacity for gene exchange, is widely used to classify sexually reproducing eukaryotes but is generally thought to be inapplicable to bacteria due to their completely asexual mode of reproduction. We show that the genus Salmonella, whose thousands of described serovars were formerly considered to be strictly clonal, undergoes sufficient levels of homologous recombination to be assigned to species according to the BSC. Aside from the two recognized species, Salmonella enterica and Salmonella bongori, several (and likely all) of the subspecies within S. enterica are reproductively isolated from one another and should each be considered a separate biological species. These findings demonstrate that species barriers in bacteria can form despite high levels of nucleotide identity and that commonly applied thresholds of genomic sequence identity are not reliable indicators of bacterial species status.

Uncharacterized and lineage-specific accessory genes within the Proteus mirabilis pan-genome landscape

Proteus mirabilis is a Gram-negative bacterium recognized for its unique swarming motility and urease activity. A previous proteomic report on four strains hypothesized that, unlike other Gram-negative bacteria, P. mirabilis may not exhibit significant intraspecies variation in gene content. However, there has not been a comprehensive analysis of large numbers of P. mirabilis genomes from various sources to support or refute this hypothesis. We performed comparative genomic analysis on 2,060 Proteus genomes. We sequenced the genomes of 893 isolates recovered from clinical specimens from three large US academic medical centers, combined with 1,006 genomes from NCBI Assembly and 161 genomes assembled from Illumina reads in the public domain. We used average nucleotide identity (ANI) to delineate species and subspecies, core genome phylogenetic analysis to identify clusters of highly related P. mirabilis genomes, and pan-genome annotation to identify genes of interest not present in the model P. mirabilis strain HI4320. Within our cohort, Proteus is composed of 10 named species and 5 uncharacterized genomospecies. P. mirabilis can be subdivided into three subspecies; subspecies 1 represented 96.7% (1,822/1,883) of all genomes. The P. mirabilis pan-genome includes 15,399 genes outside of HI4320, and 34.3% (5,282/15,399) of these genes have no putative assigned function. Subspecies 1 is composed of several highly related clonal groups. Prophages and gene clusters encoding putatively extracellular-facing proteins are associated with clonal groups. Uncharacterized genes not present in the model strain P. mirabilis HI4320 but with homology to known virulence-associated operons can be identified within the pan-genome. IMPORTANCE Gram-negative bacteria use a variety of extracellular facing factors to interact with eukaryotic hosts. Due to intraspecies genetic variability, these factors may not be present in the model strain for a given organism, potentially providing incomplete understanding of host-microbial interactions. In contrast to previous reports on P. mirabilis, but similar to other Gram-negative bacteria, P. mirabilis has a mosaic genome with a linkage between phylogenetic position and accessory genome content. P. mirabilis encodes a variety of genes that may impact host-microbe dynamics beyond what is represented in the model strain HI4320. The diverse, whole-genome characterized strain bank from this work can be used in conjunction with reverse genetic and infection models to better understand the impact of accessory genome content on bacterial physiology and pathogenesis of infection.

Genomic and phenotypic features of Acinetobacter baumannii isolated from oil reservoirs reveal a novel subspecies specialized in degrading hazardous hydrocarbons

The genus Acinetobacter encompasses biotechnologically relevant species and nosocomial pathogens. In this study, nine isolates recovered from different oil reservoir samples showed the ability to grow with petroleum as the only carbon source and possessed the ability to emulsify kerosene. The whole genomes of the nine strains were sequenced and analyzed. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values of all strains were compared to the reference strains, and the results were below the reference values (<97.88 and 82, respectively), suggesting that the isolates belong to a new subspecies of Acinetobacter baumannii. The name Acinetobacter baumannii oleum ficedula is proposed. A comparison of the whole genome repertoire of 290 Acinetobacter species indicated that the strains in this study resemble non-pathogenic Acinetobacter strains. However, the new isolates resemble A. baumannii when comparing virulence factors. The isolates in this study carry many genes involved in hydrocarbon degradation, indicating the potential to degrade most toxic compounds listed by environmental regulatory agencies such as ATSDR, EPA, and CONAMA. In addition, despite the absence of known biosurfactant or bioemulsifier genes, the strains showed emulsifying activity, suggesting the presence of new pathways or genes related to this process. This study investigated the genomic, phenotypic, and biochemical features of the novel environmental subspecies A. baumannii oleum ficedula, revealing their potential to degrade hydrocarbons and to produce biosurfactants or bioemulsifiers. Applying these environmental subspecies in bioaugmentation strategies sheds light on future approaches to bioremediation. The study shows the importance of genomic analysis of environmental strains and their inclusion in metabolic pathways databases, highlighting unique enzymes/alternative pathways for consuming hazardous hydrocarbons.

The temporal dynamics of antimicrobial-resistant Salmonella enterica and predominant serovars in China

Salmonella enterica is one of the most common bacterial pathogens in humans and animals. Systematic studies on the trends and geographical distribution of antimicrobial-resistant Salmonella and dominant serovars have been well studied in European and American countries while not in China. Here, taking the One-Health strategy, we used >35 000 Salmonella enterica isolates to explore the temporal and spatial dynamics of dominant serovars in China. We found that Salmonella Typhimurium was the dominant serovar causing human infection in China, which was consistent with Australia but inconsistent with North American and European countries. The proportion of Salmonella serovars Typhimurium, London, Rissen, Corvallis, Meleagridis, Kentucky, and Goldcoast showed an increasing trend during 2006-2019. We randomly selected 1962 isolates for comparative genomics and antimicrobial resistance studies and found that the number of antibiotic resistance genes (ARGs) per isolate increased 1.84 and 2.69 times of human and non-human origins, respectively, spanning 14 years. The proportion of antimicrobial-resistant Salmonella isolates had an increasing trend during 2006-2019, especially beta-lactam, quinolone, tetracycline, and rifampicin resistance. Moreover, we found that higher diversity of sequence types (STs) in S. Typhimurium than in other serovars, ST34 from pig and ST19 from chicken origin, were mainly associated with isolates causing child and adult gastro-infection, respectively. Our results fill in the data gap on the trends of dominant serovars and antimicrobial resistance of Salmonella enterica in China. These data provide useful information for public health decision-makers prioritizing interventions for foodborne diseases and food safety.

Evaluation of Genomic Typing Methods in the Salmonella Reference Laboratory in Public Health, England, 2012-2020

We aim to provide an evidence-based evaluation of whole genome sequence (WGS) methods, employed at the Salmonella reference laboratory in England, in terms of its impact on public health and whether these methods remain a fit for purpose test under UKAS ISO 15189. The evaluation of the genomic methods were mapped against the value of detecting microbiological clusters to support the investigation of food-borne outbreaks of Salmonella in England between 2012-2020. The analysis of WGS with both SNP- and allelic-based methods provided an unprecedented level of strain discrimination and detection of additional clusters when comparing to all of the previous typing methods. The robustness of the routine genomic sequencing at the reference laboratory ensured confidence in the microbiological identifications, even in large outbreaks with complex international food distribution networks. There was evidence that the phylogeny derived from the WGS data can be used to inform the provenance of strains and support discrimination between domestic and non-domestic transmission events. Further insight on the evolutionary context of the emerging pathogenic strains was enabled with a deep dive of the phylogenetic data, including the detection of nested clusters. The public availability of the WGS data linked to the clinical, epidemiological and environmental context of the sequenced strains has improved the trace-back investigations during outbreaks. The global expansion in the use of WGS-based typing in reference laboratories has shown that the WGS methods are a fit for purpose test in public health as it has ensured the rapid implementation of interventions to protect public health, informed risk assessment and has facilitated the management of national and international food-borne outbreaks of Salmonella.

Molecular insights into novel environmental strains of Klebsiella quasipneumoniae harboring different antimicrobial-resistance genes

Introduction

The emergence of bacterial pathogens in environmental hosts represents a major risk to public health. This study aimed at characterizing seven novel environmental strains of K. quasipneumoniae using a genomic approach which was misidentified by phenotypic methods in a previous batch of 27 species thought to be K. pneumoniae.

Methods

Whole-genome sequencing was performed using the Illumina platform, and the generated raw reads were de novo assembled. Comparative genomic, resistome, virulome, mobilome, and phylogeny were then investigated using dierent bioinformatics tools.

Results

Six strains were identified as K. quasipneumoniae subsp similipneumoniae and one as K. quasipneumoniae subsp. quasipneumoniae. All isolates were resistant to ampicillin, cephalexin, and amoxicillin-clavulanic acid and harbored the fosA, bla _OKP types, oqxB, and oqxA genes. One isolate additionally harbored a gene cassettes consisting of bla _SHV-1, bla _OXA-1, aac(6')-Ib-cr, catB genes. The aminoglycoside-modifying enzyme gene aph(3")-Ia was bracketed by two insertion elements. Plasmid analyses showed that IncFIB_K was the most prevalent plasmid, circulating in six isolates, while one isolate exhibited seven different plasmids. The isolates have virulence genes responsible for capsule formation, lipopolysaccharide, iron uptake aerobactin (iutA), salmochelins (iroE, iroN), enterobactin siderophore, adherence, and biofilm formation (mrkA, mrkB, mrkC, mrkD, mrkF, and mrkH).

Conclusion

Our study highlights the ecology and transmission of K. quasipneumoniae (which have the ability to disseminate to other environmental sources including animals) outside the clinical setting and the contribution of water, vegetables, and table surfaces as potential reservoirs of farm-to-fork transmission of disease via local markets in Khartoum, Sudan.

Sugar-Phosphate Toxicities Attenuate Salmonella Fitness in the Gut

Pathogens are becoming resistant to antimicrobials at an increasing rate, and novel therapeutic strategies are needed. Using Salmonella as a model, we have investigated the induction of sugar-phosphate toxicity as a potential therapeutic modality. The approach entails providing a nutrient while blocking the catabolism of that nutrient, resulting in the accumulation of a toxic intermediate. We hypothesize that this build-up will decrease the fitness of the organism during infection given nutrient availability. We tested this hypothesis using mutants lacking one of seven genes whose mutation is expected to cause the accumulation of a toxic metabolic intermediate. The araD, galE, rhaD, glpD, mtlD, manA, and galT mutants were then provided the appropriate sugars, either in vitro or during gastrointestinal infection of mice. All but the glpD mutant had nutrient-dependent growth defects in vitro, suggestive of sugar-phosphate toxicity. During gastrointestinal infection of mice, five mutants had decreased fitness. Providing the appropriate nutrient in the animal's drinking water was required to cause fitness defects with the rhaD and manA mutants and to enhance the fitness defect of the araD mutant. The galE and mtlD mutants were severely attenuated regardless of the nutrient being provided in the drinking water. Homologs of galE are widespread among bacteria and in humans, rendering the specific targeting of bacterial pathogens difficult. However, the araD, mtlD, and rhaD genes are not present in humans, appear to be rare in most phyla of bacteria, and are common in several genera of Enterobacteriaceae, making the encoded enzymes potential narrow-spectrum therapeutic targets. IMPORTANCE Bacterial pathogens are becoming increasingly resistant to antibiotics. There is an urgent need to identify novel drug targets and therapeutic strategies. In this work we have assembled and characterized a collection of mutations in our model pathogen, Salmonella enterica, that block a variety of sugar utilization pathways in such a way as to cause the accumulation of a toxic sugar-phosphate. Mutations in three genes, rhaD, araD, and mtlD, dramatically decrease the fitness of Salmonella in a mouse model of gastroenteritis, suggesting that RhaD, AraD, and MtlD may be good narrow-spectrum drug targets. The induction of sugar-phosphate toxicities may be a therapeutic strategy that is broadly relevant to other bacterial and fungal pathogens.

EnteroBase: hierarchical clustering of 100 000s of bacterial genomes into species/subspecies and populations

The definition of bacterial species is traditionally a taxonomic issue while bacterial populations are identified by population genetics. These assignments are species specific, and depend on the practitioner. Legacy multilocus sequence typing is commonly used to identify sequence types (STs) and clusters (ST Complexes). However, these approaches are not adequate for the millions of genomic sequences from bacterial pathogens that have been generated since 2012. EnteroBase (http://enterobase.warwick.ac.uk) automatically clusters core genome MLST allelic profiles into hierarchical clusters (HierCC) after assembling annotated draft genomes from short-read sequences. HierCC clusters span core sequence diversity from the species level down to individual transmission chains. Here we evaluate HierCC's ability to correctly assign 100 000s of genomes to the species/subspecies and population levels for Salmonella, Escherichia, Clostridoides, Yersinia, Vibrio and Streptococcus. HierCC assignments were more consistent with maximum-likelihood super-trees of core SNPs or presence/absence of accessory genes than classical taxonomic assignments or 95% ANI. However, neither HierCC nor ANI were uniformly consistent with classical taxonomy of Streptococcus. HierCC was also consistent with legacy eBGs/ST Complexes in Salmonella or Escherichia and with O serogroups in Salmonella. Thus, EnteroBase HierCC supports the automated identification of and assignment to species/subspecies and populations for multiple genera. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Molecular Epidemiology and Virulence of Non-Typhoidal Salmonella in Armenia

In this work, we analysed human isolates of nontyphoidal Salmonella enterica subsp. enterica (NTS), which were collected from salmonellosis cases in Armenia from 1996 to 2019. This disease became a leading food-borne bacterial infection in the region, with the younger age groups especially affected. The isolates were characterised by serotyping, Enterobacterial Repetitive Intergenic Consensus (ERIC-PCR) typing, and whole genome sequencing (WGS). The main serotypes were S. Typhimurium, S. Enteritidis, and S. Arizonae. ERIC-PCR indicated a high degree of clonality among S. Typhimurium strains, which were also multidrug-resistant and produced extended spectrum beta-lactamases. During the study period, the frequency of S. Typhimurium and S. Arizonae isolations decreased, but with the increase in S. Enteritidis and other NTS. A total of 42 NTS isolates were subjected to WGS and explored for virulence-related traits and the corresponding genetic elements. Some virulence and genetic factors were shared by all NTS serotypes, while the main differences were attributed to the serotype-specific diversity of virulence genes, SPIs, virulence plasmids, and phages. The results indicated the variability and dynamics in the epidemiology of salmonellosis and a high virulence potential of human NTS isolates circulating in the region.

Phylogenomic analysis of the genus Leuconostoc

Leuconostoc is a genus of saccharolytic heterofermentative lactic acid bacteria that inhabit plant-derived matrices and a variety of fermented foods (dairy products, dough, milk, vegetables, and meats), contributing to desired fermentation processes or playing a role in food spoilage. At present, the genus encompasses 17 recognized species. In total, 216 deposited genome sequences of Leuconostoc were analyzed, to check the delineation of species and to infer their evolutive genealogy utilizing a minimum evolution tree of Average Nucleotide Identity (ANI) and the core genome alignment. Phylogenomic relationships were compared to those obtained from the analysis of 16S rRNA, pheS, and rpoA genes. All the phylograms were subjected to split decomposition analysis and their topologies were compared to check the ambiguities in the inferred phylogenesis. The minimum evolution ANI tree exhibited the most similar topology with the core genome tree, while single gene trees were less adherent and provided a weaker phylogenetic signal. In particular, the 16S rRNA gene failed to resolve several bifurcations and Leuconostoc species. Based on an ANI threshold of 95%, the organization of the genus Leuconostoc could be amended, redefining the boundaries of the species L. inhae, L. falkenbergense, L. gelidum, L. lactis, L. mesenteroides, and L. pseudomesenteroides. Two strains currently recognized as L. mesenteroides were split into a separate lineage representing a putative species (G16), phylogenetically related to both L. mesenteroides (G18) and L. suionicum (G17). Differences among the four subspecies of L. mesenteroides were not pinpointed by ANI or by the conserved genes. The strains of L. pseudomesenteroides were ascribed to two putative species, G13 and G14, the former including also all the strains presently belonging to L. falkenbergense. L. lactis was split into two phylogenetically related lineages, G9 and G10, putatively corresponding to separate species and both including subgroups that may correspond to subspecies. The species L. gelidum and L. gasicomitatum were closely related but separated into different species, the latter including also L. inhae strains. These results, integrating information of ANI, core genome, and housekeeping genes, complemented the taxonomic delineation with solid information on the phylogenetic lineages evolved within the genus Leuconostoc.

The Number and Type of Chaperone-Usher Fimbriae Reflect Phylogenetic Clade Rather than Host Range in Salmonella

Salmonella is one of the most successful foodborne pathogens worldwide, owing in part to its ability to colonize or infect a wide range of hosts. Salmonella serovars are known to encode a variety of different fimbriae (hairlike organelles that facilitate binding to surfaces); however, the distribution, number, and sequence diversity of fimbriae encoded across different lineages of Salmonella were unknown. We queried whole-genome sequence (WGS) data for 242 Salmonella enterica subsp. enterica (subspecies enterica) isolates from the top 217 serovars associated with isolation from humans and agricultural animals; this effort identified 2,894 chaperone-usher (CU)-type fimbrial usher sequences, representing the most conserved component of CU fimbriae. On average, isolates encoded 12 different CU fimbrial ushers (6 to 18 per genome), although the distribution varied significantly (P = 1.328E-08) by phylogenetic clade, with isolates in section Typhi having significantly fewer fimbrial ushers than isolates in clade A2 (medians = 10 and 12 ushers, respectively). Characterization of fimbriae in additional non-enterica subspecies genomes suggested that 8 fimbrial ushers were classified as being unique to subspecies enterica isolates, suggesting that the majority of fimbriae were most likely acquired prior to the divergence of subspecies enterica. Characterization of mobile elements suggested that plasmids represent an important vehicle facilitating the acquisition of a wide range of fimbrial ushers, particularly for the acquisition of fimbriae from other Gram-negative genera. Overall, our results suggest that differences in the number and type of fimbriae encoded most likely reflect differences in phylogenetic clade rather than differences in host range. IMPORTANCE Fimbriae of the CU assembly pathway represent important organelles that mediate Salmonella's interactions with host tissues and abiotic surfaces. Our analyses provide a comprehensive overview of the diversity of CU fimbriae in Salmonella spp., highlighting that the majority of CU fimbriae are distributed broadly across multiple subspecies and suggesting that acquisition most likely occurred prior to the divergence of subspecies enterica. Our data also suggest that plasmids represent the primary vehicles facilitating the horizontal transfer of diverse CU fimbriae in Salmonella. Finally, the observed high sequence similarity between some ushers suggests that different names may have been assigned to closely related fimbrial ushers that likely should be represented by a single designation. This highlights the need to establish standard criteria for fimbria classification and nomenclature, which will also facilitate future studies seeking to associate virulence factors with adaptation to or differences in the likelihood of causing disease in a given host.

Comprehensive genome analysis of Lentzea reveals repertoire of polymer-degrading enzymes and bioactive compounds with clinical relevance

The genus Lentzea is a rare group of actinobacteria having potential for the exploration of bioactive compounds. Despite its proven ability to produce compounds with medical relevance, Lentzea genome analysis remains unexplored. Here we show a detailed understanding of the genetic features, biosynthetic gene clusters (BGCs), and genetic clusters for carbohydrate-active enzymes present in the Lentzea genome. Our analysis determines the genes for core proteins, non-ribosomal peptide synthetase condensation domain, and polyketide synthases-ketide synthase domain. The antiSMASH-based sequence analysis identifies 692 BGCs among which 8% are identical to the BGCs that produce geosmin, citrulassin, achromosin (lassopeptide), vancosamine, anabaenopeptin NZ857/nostamide A, alkylresorcinol, BE-54017, and bezastatin. The remaining BGCs code for advanced category antimicrobials like calcium-dependent, glycosylated, terpenoids, lipopeptides, thiopeptide, lanthipeptide, lassopeptide, lingual antimicrobial peptide and lantibiotics together with antiviral, antibacterial, antifungal, antiparasitic, anticancer agents. About 28% of the BGCs, that codes for bioactive secondary metabolites, are exclusive in Lentzea and could lead to new compound discoveries. We also find 7121 genes that code for carbohydrate-degrading enzymes which could essentially convert a wide range of polymeric carbohydrates. Genome mining of such genus is very much useful to give scientific leads for experimental validation in the discovery of new-generation bioactive molecules of biotechnological importance.