Serratia marcescens is an important nosocomial pathogen causing various opportunistic infections, such as urinary tract infections, bacteremia and sometimes even hospital outbreaks. The recent emergence and spread of multidrug-resistant (MDR) strains further pose serious threats to global public health. This bacterium is also ubiquitously found in natural environments, but the genomic differences between clinical and environmental isolates are not clear, including those between S. marcescens and its close relatives. In this study, we performed a large-scale genome analysis of S. marcescens and closely related species (referred to as the ‘ S. marcescens complex’), including more than 200 clinical and environmental strains newly sequenced here. Our analysis revealed their phylogenetic relationships and complex global population structure, comprising 14 clades, which were defined based on whole-genome average nucleotide identity. Clades 10, 11, 12 and 13 corresponded to S. nematodiphila , S. marcescens sensu stricto, S. ureilytica and S. surfactantfaciens, respectively. Several clades exhibited distinct genome sizes and GC contents and a negative correlation of these genomic parameters was observed in each clade, which was associated with the acquisition of mobile genetic elements (MGEs), but different types of MGEs, plasmids or prophages (and other integrative elements), were found to contribute to the generation of these genomic variations. Importantly, clades 1 and 2 mostly comprised clinical or hospital environment isolates and accumulated a wide range of antimicrobial resistance genes, including various extended-spectrum β-lactamase and carbapenemase genes, and fluoroquinolone target site mutations, leading to a high proportion of MDR strains. This finding suggests that clades 1 and 2 represent hospital-adapted lineages in the S. marcescens complex although their potential virulence is currently unknown. These data provide an important genomic basis for reconsidering the classification of this group of bacteria and reveal novel insights into their evolution, biology and differential importance in clinical settings.

In the past two decades, Shiga toxin-producing Escherichia coli (STEC) has been responsible for multiple large-scale outbreaks worldwide, affecting thousands of individuals. While surveillance systems in developed countries such as the United States, the United Kingdom, Europe, Australia, Japan, and Canada are well-established, data on STEC prevalence in developing nations remain sparse, partly due to the absence of well-structured molecular diagnostic networks or surveillance systems. This review analyzed 250 studies published between 2014 and 2024 across 39 developing countries in Africa, Asia, Latin America, and the Caribbean, yielding 8986 STEC isolates. Detailed serogroup and serotype data were available for 55.9% of these, with O111, O157, and O26 being most common in humans. In animals, O157:H7 was most frequent, while food isolates mirrored global trends with O157 and O111 dominance. Notably, O145, a serogroup frequently reported in the U.S. and Europe, was absent from the ‘’Top Seven’’ serogroups. Shiga toxin subtypes stx1a and stx2a were most prevalent in human cases. In animal isolates, stx2e was the most prevalent subtype, while stx2c was most commonly found in food samples. We recommend establishing reference laboratories in these regions to improve data quality, strengthen monitoring efforts, and reduce the burden of STEC infections globally.

Bacteriophages (or phages) play major roles in the evolution of bacterial pathogens via horizontal gene transfer. Multiple phages are often integrated in a host chromosome as prophages, not only carrying various novel virulence-related genetic determinants into host bacteria but also providing various possibilities for prophage-prophage interactions in bacterial cells. In particular, Escherichia coli strains such as Shiga toxin (Stx)-producing E. coli (STEC) and enteropathogenic E. coli (EPEC) strains have acquired more than 10 prophages (up to 21 prophages), many of which encode type III secretion system (T3SS) effector gene clusters. In these strains, some prophages are present at a single locus in tandem, which is usually interpreted as the integration of phages that use the same attachment (att) sequence. Here, we present phages integrating into T3SS effector gene cluster-associated loci in prophages, which are widely distributed in STEC and EPEC. Some of the phages integrated into prophages are Stx-encoding phages (Stx phages) and have induced the duplication of Stx phages in a single cell. The identified attB sequences in prophage genomes are apparently derived from host chromosomes. In addition, two or three different attB sequences are present in some prophages, which results in the generation of prophage clusters in various complex configurations. These phages integrating into prophages represent a medically and biologically important type of inter-phage interaction that promotes the accumulation of T3SS effector genes in STEC and EPEC, the duplication of Stx phages in STEC, and the conversion of EPEC to STEC and that may be distributed in other types of E. coli strains as well as other prophage-rich bacterial species.

Shiga toxin (Stx) is the key virulence factor of Stx-producingEscherichia coli(STEC). All known Stxs (Stx1 and Stx2) are encoded by bacteriophages (Stx phages). Although the genetic diversity of Stx phages has frequently been described, systematic analyses of Stx phages in a single STEC lineage are limited. In this study, focusing on the O26:H11 STEC sequence type 21 (ST21) lineage, where thestx1agene is highly conserved, we analysed the Stx1a phages in 39 strains representative of the entire ST21 lineage and found a high level of variation in Stx1a phage genomes caused by various mechanisms, including replacement by a different Stx1a phage at the same or different locus. The evolutionary timescale of events changing Stx1a phages in ST21 was also determined. Furthermore, by using an Stx1 quantification system developed in this study, we found notable variations in the efficiency of Stx1 production upon prophage induction, which sharply contrasted with the conserved iron regulated Stx1 production. These variations were associated with the Stx1a phage alteration in some cases but not in other cases; thus, Stx1 production in this STEC lineage was determined by differences not only in Stx1 phages but also in host-encoded factors.

Shiga toxin (Stx)-producing Escherichia coli (STEC) are foodborne pathogens causing serious diseases, such as haemorrhagic colitis and haemolytic uraemic syndrome. Although O157:H7 STEC strains have been the most prevalent, incidences of STEC infections by several other serotypes have recently increased. O121:H19 STEC is one of these major non-O157 STECs, but systematic whole genome sequence (WGS) analyses have not yet been conducted on this STEC. Here, we performed a global WGS analysis of 638 O121:H19 strains, including 143 sequenced in this study, and a detailed comparison of 11 complete genomes, including four obtained in this study. By serotype-wide WGS analysis, we found that O121:H19 strains were divided into four lineages, including major and second major lineages (named L1 and L3, respectively), and that the locus of enterocyte effacement (LEE) encoding a type III secretion system (T3SS) was acquired by the common ancestor of O121:H19. Analyses of 11 complete genomes belonging to L1 or L3 revealed remarkable interlineage differences in the prophage pool and prophage-encoded T3SS effector repertoire, independent acquisition of virulence plasmids by the two lineages, and high conservation in the prophage repertoire, including that for Stx2a phages in lineage L1. Further sequence determination of complete Stx2a phage genomes of 49 strains confirmed that Stx2a phages in lineage L1 are highly conserved short-tailed phages, while those in lineage L3 are long-tailed lambda-like phages with notable genomic diversity, suggesting that an Stx2a phage was acquired by the common ancestor of L1 and has been stably maintained. Consistent with these genomic features of Stx2a phages, most lineage L1 strains produced much higher levels of Stx2a than lineage L3 strains. Altogether, this study provides a global phylogenetic overview of O121:H19 STEC and shows the interlineage genomic differences and the highly conserved genomic features of the major lineage within this serotype of STEC.

Among Shiga toxin (Stx)-producing Escherichia coli (STEC) strains of various serotypes, O157:H7 and five major non-O157 STEC (O26:H11, O111:H8, O103:H2, O121:H19 and O145:H28) can be selectively isolated by using tellurite-containing media. While human infections by O165:H25 STEC strains have been reported worldwide, their detection and isolation are not easy, as they are not resistant to tellurite. Systematic whole-genome sequencing (WGS) analyses have not yet been conducted. Here, we defined O165:H25 strains and their close relatives, including O172:H25 strains, as clonal complex 119 (CC119) and performed a global WGS analysis of the major lineage of CC119, called CC119 sensu stricto (CC119ss), by using 202 CC119ss strains, including 90 strains sequenced in this study. Detailed comparisons of 13 closed genomes, including 7 obtained in this study, and systematic analyses of Stx phage genomes in 50 strains covering the entire CC119ss lineage, were also conducted. These analyses revealed that the Stx2a phage, the locus of enterocyte effacement (LEE) encoding a type III secretion system (T3SS), many prophages encoding T3SS effectors, and the virulence plasmid were acquired by the common ancestor of CC119ss and have been stably maintained in this lineage, while unusual exchanges of Stx1a and Stx2c phages were found at a single integration site. Although the genome sequences of Stx2a phages were highly conserved, CC119ss strains exhibited notable variation in Stx2 production levels. Further analyses revealed the lack of SpLE1-like elements carrying the tellurite resistance genes in CC119ss and defects in rhamnose, sucrose, salicin and dulcitol fermentation. The genetic backgrounds underlying these defects were also clarified.