Global population structure of the Serratia marcescens complex and identification of hospital-adapted lineages in the complex

The ICU environment contributes to the endemicity of the "Serratia marcescens complex" in the hospital setting

Serratia marcescens is an opportunistic pathogen historically associated with sudden outbreaks in intensive care units (ICUs) and the spread of carbapenem-resistant genes. However, the ecology of S. marcescens populations in the hospital ecosystem remains largely unknown. We combined epidemiological information of 1,432 Serratia spp. isolates collected from sinks of a large ICU that underwent demographic and operational changes (2019-2021) and 99 non-redundant outbreak/non-outbreak isolates from the same hospital (2003-2019) with 165 genomic data. These genomes were grouped into clades (1-4) and subclades (A and B) associated with distinct species: Serratia nematodiphila (1A), S. marcescens (1B), Serratia bockelmannii (2A), Serratia ureilytica (2B), S. marcescens/Serratia nevei (3), and S. nevei (4A and 4B). They may be classified into an S. marcescens complex (SMC) due to the similarity between/within subclades (average nucleotide identity >95%-98%), with clades 3 and 4 predominating in our study and publicly available databases. Chromosomal AmpC β-lactamase with unusual basal-like expression and prodigiosin-lacking species contrasted classical features of Serratia. We found persistent and coexisting clones in sinks of subclades 4A (ST92 and ST490) and 4B (ST424), clonally related to outbreak isolates carrying blaVIM-1 or blaOXA-48 on prevalent IncL/pB77-CPsm plasmids from our hospital since 2017. The distribution of SMC populations in ICU sinks and patients reflects how Serratia species acquire, maintain, and enable plasmid evolution in both "source" (permanent, sinks) and "sink" (transient, patients) hospital patches. The results contribute to understanding how water sinks serve as reservoirs of Enterobacterales clones and plasmids that enable the persistence of carbapenemase genes in healthcare settings, potentially leading to outbreaks and/or hospital-acquired infections.IMPORTANCEThe "hospital environment," including sinks and surfaces, is increasingly recognized as a reservoir for bacterial species, clones, and plasmids of high epidemiological concern. Available studies on Serratia epidemiology have focused mainly on outbreaks of multidrug-resistant species, overlooking local longitudinal analyses necessary for understanding the dynamics of opportunistic pathogens and antibiotic-resistant genes within the hospital setting. This long-term genomic comparative analysis of Serratia isolated from the ICU environment with isolates causing nosocomial infections and/or outbreaks within the same hospital revealed the coexistence and persistence of Serratia populations in water reservoirs. Moreover, predominant sink strains may acquire highly conserved and widely distributed plasmids carrying carbapenemase genes, such as the prevalent IncL-pB77-CPsm (pOXA48), persisting in ICU sinks for years. The work highlights the relevance of ICU environmental reservoirs in the endemicity of certain opportunistic pathogens and resistance mechanisms mainly confined to hospitals.

Genetic barriers more than environmental associations explain Serratia marcescens population structure

Bacterial species often comprise well-separated lineages, likely emerged and maintained by genetic isolation and/or ecological divergence. How these two evolutionary actors interact in the shaping of bacterial population structure is currently not fully understood. In this study, we investigate the genetic and ecological drivers underlying the evolution of Serratia marcescens, an opportunistic pathogen with high genomic flexibility and able to colonise diverse environments. Comparative genomic analyses reveal a population structure composed of five deeply-demarcated genetic clusters with open pan-genome but limited inter-cluster gene flow, partially explained by Restriction-Modification (R-M) systems incompatibility. Furthermore, a large-scale research on hundred-thousands metagenomic datasets reveals only a partial habitat separation of the clusters. Globally, two clusters only show a separate gene composition coherent with ecological adaptations. These results suggest that genetic isolation has preceded ecological adaptations in the shaping of the species diversity, an evolutionary scenario coherent with the Evolutionary Extended Synthesis.

Systematic analysis of plasmids of the Serratia marcescens complex using 142 closed genomes

Plasmids play important roles in bacterial genome diversification. In the Serratia marcescens complex (SMC), a notable contribution of plasmids to genome diversification was also suggested by our recent analysis of >600 draft genomes. As accurate analyses of plasmids in draft genomes are difficult, in this study we analysed 142 closed genomes covering the entire complex, 67 of which were obtained in this study, and identified 132 plasmids (1.9-244.4 kb in length) in 77 strains. While the average numbers of plasmids in clinical and non-clinical strains showed no significant difference, strains belonging to clade 2 (one of the two hospital-adapted lineages) contained more plasmids than the others. Pangenome analysis revealed that of the 28 954 genes identified, 12.8 % were plasmid-specific, and 1.4 % were present in plasmids or chromosomes depending on the strain. In the latter group, while transposon-related genes were most prevalent (31.4 % of the function-predicted genes), genes related to antimicrobial resistance and heavy metal resistance accounted for a notable proportion (22.7 %). Mash distance-based clustering separated the 132 plasmids into 23 clusters and 50 singletons. Most clusters/singletons showed notably different GC contents compared to those of host chromosomes, suggesting their recent or relatively recent appearance in the SMC. Among the 23 clusters, 17 were found in only clinical or only non-clinical strains, suggesting the possible preference of their distribution on the environmental niches of host strains. Regarding the host strain phylogeny, 16 clusters were distributed in two or more clades, suggesting their interclade transmission. Moreover, for many plasmids, highly homologous plasmids were found in other species, indicating the broadness of their potential host ranges, beyond the genus, family, order, class or even phylum level. Importantly, highly homologous plasmids were most frequently found in Klebsiella pneumoniae and other species in the family Enterobacteriaceae, suggesting that this family, particularly K. pneumoniae, is the main source for plasmid exchanges with the SMC. These results highlight the power of closed genome-based analysis in the investigation of plasmids and provide important insights into the nature of plasmids distributed in the SMC.

The man, the plant, and the insect: shooting host specificity determinants in Serratia marcescens pangenome

Introduction

Serratia marcescens is most commonly known as an opportunistic pathogen causing nosocomial infections. It, however, was shown to infect a wide range of hosts apart from vertebrates such as insects or plants as well, being either pathogenic or growth-promoting for the latter. Despite being extensively studied in terms of virulence mechanisms during human infections, there has been little evidence of which factors determine S. marcescens host specificity. On that account, we analyzed S. marcescens pangenome to reveal possible specificity factors.

Methods

We selected 73 high-quality genome assemblies of complete level and reconstructed the respective pangenome and reference phylogeny based on core genes alignment. To find an optimal pipeline, we tested current pangenomic tools and obtained several phylogenetic inferences. The pangenome was rich in its accessory component and was considered open according to the Heaps' law. We then applied the pangenome-wide associating method (pan-GWAS) and predicted positively associated gene clusters attributed to three host groups, namely, humans, insects, and plants.

Results

According to the results, significant factors relating to human infections included transcriptional regulators, lipoproteins, ABC transporters, and membrane proteins. Host preference toward insects, in its turn, was associated with diverse enzymes, such as hydrolases, isochorismatase, and N-acetyltransferase with the latter possibly exerting a neurotoxic effect. Finally, plant infection may be conducted through type VI secretion systems and modulation of plant cell wall synthesis. Interestingly, factors associated with plants also included putative growth-promoting proteins like enzymes performing xenobiotic degradation and releasing ammonium irons. We also identified overrepresented functional annotations within the sets of specificity factors and found that their functional characteristics fell into separate clusters, thus, implying that host adaptation is represented by diverse functional pathways. Finally, we found that mobile genetic elements bore specificity determinants. In particular, prophages were mainly associated with factors related to humans, while genetic islands-with insects and plants, respectively.

Discussion

In summary, functional enrichments coupled with pangenomic inferences allowed us to hypothesize that the respective host preference is carried out through distinct molecular mechanisms of virulence. To the best of our knowledge, the presented research is the first to identify specific genomic features of S. marcescens assemblies isolated from different hosts at the pangenomic level.

The accessory protein TagV is required for full Type VI secretion system activity in Serratia marcescens

The bacterial Type VI secretion system (T6SS) is a dynamic macromolecular structure that promotes inter- and intra-species competition through the delivery of toxic effector proteins into neighbouring cells. The T6SS contains 14 well-characterised core proteins necessary for effector delivery (TssA-M, PAAR). In this study, we have identified a novel accessory component required for optimal T6SS activity in the opportunistic pathogen Serratia marcescens, which we name TagV. Deletion of tagV, which encodes an outer membrane lipoprotein, caused a reduction in the T6SS-dependent antibacterial activity of S. marcescens Db10. Mutants of S. marcescens lacking the core component TssJ, a distinct outer membrane lipoprotein previously considered essential for T6SS firing, retained a modest T6SS activity that could be abolished through deletion of tagV. TagV did not interact with the T6SS membrane complex proteins TssL or TssM, but is proposed to bind to peptidoglycan, indicating that the mechanism by which TagV promotes T6SS firing differs from that of TssJ. Homologues of tagV were identified in several other bacterial genera, suggesting that the accessory function of TagV is not restricted to S. marcescens. Together, our findings support the existence of a second, TssJ-independent mechanism for T6SS firing that is dependent upon the activity of TagV proteins.

[Genome analysis-based studies on bacterial genetic diversity]

There are a huge number of bacterial species on earth, and a huge intra-species genomic diversity are also observed in many bacteria. The high ability of bacteria to acquire foreign DNA and the presence of various mobile genetic elements contribute the generation of such genomic diversity. During the biochemical and genetic analysis of a Pseudomonas aeruginosa toxin, called cytotoxin, and its converting phage, which I first engaged in my research carrier, I became very interested in the genetic diversity of bacteria and mobile genetic elements such as bacteriophages, and realized the usefulness and power of genome analysis. Since then, I have been involved in genome analyses of various pathogenic bacteria such as enterohemorrhagic Escherichia coli (EHEC), commensal bacteria of human and other animals, and bacteria or bacterial communities in natural environments. I was so lucky that I jumped in this research field at the very begging of genome analyses and experienced a very exciting time of surprisingly rapid advance in genome sequencing technologies which revolutionized a wide range of biology. In this article, I first review the main findings which our group obtained from the genome analyses on the P. aeruginosa cytotoxin converting phage and those on the evolution and genomic diversity of EHEC and related bacteria. The results of our analyses of Rickettsiaceae family genomes, which show surprisingly very low genomic diversity, and genome sequence-based analyses of an intrahospital bacterial outbreak and within-host genomic diversity are also summarized.