Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 16 January 2026

Sec. Antimicrobials, Resistance and Chemotherapy

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1718340

This article is part of the Research TopicIntegrative Microbial and Chemical Genomics to Decipher Antibiotic Resistance Mechanisms and Developing Innovative Antimicrobial ApproachesView all 9 articles

Whole-genome sequencing-based source tracing and infection control of Serratia marcescens blood culture events in pediatric patients


Wang Zhang*Wang Zhang1*Caihua MaCaihua Ma1Falin XuFalin Xu2Chenjing ZhaoChenjing Zhao2Ling Wang*Ling Wang1*
  • 1Department of Infection Prevention and Control, The Third Affiliated Hospital of Zhengzhou University, Zhengzhou, China
  • 2Department of Pediatrics, The Third Affiliated Hospital of Zhengzhou University, Zhengzhou, China

Objective: Serratia marcescens (S. marcescens) is an opportunistic pathogen increasingly associated with nosocomial infections in immunocompromised pediatric patients. This study aimed to investigate the epidemiological and molecular epidemiological characteristics of S. marcescens-positive blood cultures and to provide evidence for targeted infection prevention strategies.

Methods: Eleven cases of S. marcescens-positive blood cultures were identified across eight pediatric wards of a tertiary hospital in China in October 2024. Clinical and microbiological assessments were conducted to distinguish true infections from contamination. Antimicrobial susceptibility testing was performed, and whole-genome sequencing (WGS) was used to characterize resistance, virulence, and plasmid replicons. Core genome single-nucleotide polymorphism (SNP)-based phylogenies, Bayesian temporal inference, and transmission tree reconstruction were used to explore genetic relatedness, transmission dynamics, and cryptic cases.

Results: Two cases were diagnosed as hospital-acquired sepsis, while nine were classified as contamination. The isolates exhibited intrinsic resistance to first- and second-generation cephalosporins but remained susceptible to carbapenems. Seven resistance genes—including aac(6’)-Ic, oqxB, and tet(41)—and four virulence genes, including cheY and fliM, were identified. Nine isolates carried IncFII-type plasmids. Core genome SNP analysis revealed minimal genetic divergence, with the most recent common ancestor traceable to late 2023, suggesting approximately 10 months of silent transmission. Transmission tree inference further indicated the presence of undetected cryptic cases. Following the implementation of bundled interventions—including sink replacements and disinfectant use, reinforcement of hand hygiene, and environmental decontamination—no new cases occurred during 1 month of follow-up.

Conclusion: This study identified a highly clonal S. marcescens cluster with evidence of prolonged, unnoticed circulation and cross-ward transmission in the pediatric units. These findings underscore the hidden persistence of this pathogen in the hospital environment and the need for strengthened infection control measures.

Introduction

Serratia marcescens is a Gram-negative, facultatively anaerobic bacillus commonly found in natural environments and in moist areas of healthcare facilities, such as faucets, wet towels, and various medical devices (Bourdin et al., 2023). Although it exhibits limited pathogenicity in immunocompetent individuals, it can cause severe healthcare-associated infections, including bloodstream infections, pneumonia, and urinary tract infections, in neonates, critically ill patients, and individuals receiving prolonged invasive procedures (Crivaro et al., 2007). In recent years, S. marcescens has gained increasing clinical attention due to its multidrug resistance and diverse virulence factors, with particularly elevated transmission risk in pediatric and neonatal wards (Reddy et al., 2021; Arslan et al., 2010). The organism frequently shows resistance to cephalosporins, aminoglycosides, and fluoroquinolones (Martineau et al., 2018), and utilizes flagellar proteins and secretion systems to enhance adhesion and immune evasion (Ghaith et al., 2018). Moreover, its strong environmental persistence enables prolonged survival and potential formation of silent transmission chains within hospitals, complicating outbreak control (Shaik Ismail et al., 2025).

Pediatric wards represent high-risk settings for healthcare-associated infections because children, especially neonates, have immature immune systems and frequently undergo invasive procedures (Shane et al., 2017; Patel and Kaufman, 2015). During blood collection, agitation and movement in infants often result in inadequate skin disinfection, contributing to contamination of blood culture specimens (Zackeroff et al., 2024; Pimentel de Araujo et al., 2017). Although S. marcescens is one of the Serratia species with the most abundant genomic data available (Williams et al., 2022), its potential silent circulation in pediatric units is often underestimated (Tsagris et al., 2012; Cristina et al., 2019; Bechmann et al., 2022). In cases lacking typical clinical manifestations, undetected carriers or unnoticed links in transmission chains may further obscure epidemiological investigations (Parker et al., 2023). Therefore, elucidating the dissemination characteristics of S. marcescens in pediatric settings is essential for strengthening early warning and infection-control strategies (Mosalli et al., 2022).

Traditional infection-control investigations rely on epidemiological tracing and environmental surveillance but often fail to resolve fine-scale genetic relationships between isolates or identify cryptic transmission events (De Silva et al., 2016; Fitzpatrick et al., 2016). Whole-Genome Sequencing (WGS) enables precise strain classification (Harch et al., 2025) and has become widely used in healthcare-associated infection source-tracking (Sundermann et al., 2022). High-resolution core-genome single-nucleotide polymorphism (SNP) analysis can reveal phylogenetic relatedness among isolates and, together with virulence gene, antimicrobial resistance gene, and plasmid profiling, allows accurate reconstruction of transmission pathways (Alam et al., 2015; Mitchell et al., 2024). Bayesian modeling and transmission-dynamic analysis can further estimate the duration of transmission and infer undetected cases (Didelot et al., 2014). The application of WGS has shifted infection-control practices from experience-based assessment to molecularly informed, high-precision identification (Brlek et al., 2024).

This study investigates a cluster of 11 S. marcescens bloodstream culture-positive events that occurred in the pediatric department of a tertiary hospital in 2024. By applying WGS, we aimed to differentiate true infections from contamination, clarify strain-relatedness and transmission characteristics, identify potential silent transmitters, and provide evidence to improve the interpretation of pediatric blood culture results and optimize infection-prevention strategies.

Materials and methods

Patient source

In October 2024, a total of 11 cases of S. marcescens-positive blood cultures were reported across eight pediatric wards in a tertiary hospital. Among these, three cases were identified in the Department of Pediatric Neurology, two in Pediatric Hematology, and one case each in Pediatric Cardionephrology and Endocrinology, Pediatric Respiratory Medicine Ward 1, Pediatric Respiratory Medicine Ward 2, the Department of Infectious Diseases, Pediatric Cardiothoracic Surgery, and Pediatric Orthopedics. This retrospective observational study was reviewed and approved by the hospital’s institutional ethics committee. The overall study design and analytical workflow are illustrated in Supplementary Figure 1, encompassing case collection, epidemiological investigation, laboratory testing, WGS and molecular analyses, transmission dynamics inference, as well as evaluation of intervention measures and clinical outcomes.

Criteria for determining infection or contamination

The diagnostic criteria for hospital-acquired infections were based on the Trial Guidelines for the Diagnosis of Hospital Infections issued by the Ministry of Health of the People’s Republic of China in 2001. The classification of blood culture results as infection or contamination was performed in accordance with the 2022 Expert Consensus on the Clinical Practice of Blood Culture Techniques for the Diagnosis of Bloodstream Infections. The comprehensive assessment included: (1) the patient’s clinical manifestations and laboratory results; (2) whether the case met the diagnostic criteria for bloodstream infection; (3) whether features indicative of contamination were present, such as single-bottle positivity, lack of corresponding clinical evidence of infection, or absence of infection during follow-up. Final determinations were made jointly by at least two senior infectious disease specialists.

Epidemiological investigation

An epidemiological investigation team was promptly established following the incident. The team comprised personnel from the Department of Hospital Infection Control, the Medical Affairs Office, the Nursing Department, the Laboratory Department, and clinical staff from the affected wards. Using the Xinglin Real-Time Hospital Infection Surveillance System (Hangzhou Xinglin Information Technology Co., Ltd., hereafter referred to as the Xinglin system), the team retrieved clinical data and relevant information for all patients with S. marcescens-positive blood cultures.

The investigation protocol included: (1) verification of the classification of each case as either true infection or contamination; (2) collection of demographic information and clinical data for each patient; (3) analysis of temporal, spatial, and population distribution patterns. In addition, interviews were conducted with attending physicians, nurses, specimen transport personnel, and laboratory microbiologists from the involved departments to review procedures for specimen collection, transportation, and identification, and to identify potential risk factors associated with infection or contamination.

Environmental hygiene surveillance

Environmental sampling was conducted for all 11 S. marcescens-positive patients and their corresponding wards. Samples were obtained from the patients themselves, healthcare workers who had direct contact with the patients, general service staff, specimen couriers, as well as the surrounding patient environment and shared ward areas. Sampling sites included: hands of patients and other personnel; frequently touched surfaces such as bedside tables and headboard-mounted equipment; ward sinks; and in-use disinfectants. Following collection, samples were inoculated onto Columbia blood agar plates and chocolate blood agar plates (Zhengzhou Antu Bioengineering Co., Ltd.) and incubated at 35–37°C for 24–48 h. Bacterial identification was subsequently performed based on culture results.

Isolate identification and antimicrobial susceptibility testing

Bacterial identification was performed using the Auto MS 1,000 automated matrix-assisted laser desorption/ionization time-of-flight mass spectrometry system (Zhengzhou Antu Bioengineering Co., Ltd.). Antimicrobial susceptibility testing was conducted using the Kirby-Bauer disk diffusion method, and inhibition zone diameters were interpreted according to the 2023 guidelines of the Clinical and Laboratory Standards Institute (CLSI). Results were classified as susceptible (S), intermediate (I), or resistant (R). Escherichia coli ATCC 25922, provided by the National Center for Clinical Laboratories under the National Health Commission, was used as the quality control strain.

WGS and assembly of isolates

A total of 11 S. marcescens isolates were collected and subjected to WGS. Following genomic DNA extraction, sequencing libraries were prepared using the Illumina TruSeq DNA PCR-Free Library Prep Kit and sequenced on the Illumina NovaSeq 6000 platform (Beijing Novogene Bioinformatics Technology Co., Ltd.) with paired-end 150 bp reads (PE150). The average sequencing depth per isolate was approximately 100×. Raw reads were quality-controlled using fastp v0.23, with a quality threshold of Q30 ≥ 85% to ensure data reliability.

High-quality sequencing reads were de novo assembled using a complementary multi-assembler strategy. SOAPdenovo v2.04 and ABySS v2.3 were primarily used to construct draft chromosomal assemblies and optimize long-range scaffold structure from high-coverage short-read data, while SPAdes v3.15 was used to improve the assembly of repetitive regions and low-abundance plasmid contigs. The three short-read assemblies were then integrated using CISA, and the resulting scaffolds were polished with GapCloser v1.12 to fill gaps and correct local errors. For selected representative isolates, long-read sequencing data generated with the Oxford Nanopore MinION platform were incorporated, and hybrid assemblies were performed using Unicycler v0.4.9 to increase the continuity of complex genomic regions and plasmid structures.

Analysis of antibiotic resistance genes, virulence genes, and plasmid replicons

Virulence genes in S. marcescens isolates were identified using the Virulence Factor Database (VFDB),1 while resistance genes were detected through the Comprehensive Antibiotic Resistance Database (CARD).2 BLAST-based screening was conducted using ABRicate v0.8.10, with gene presence confirmed when both sequence identity and coverage exceeded 80%.

Plasmid detection was performed using PlasmidFinder v2.1, with the Enterobacteriaceae database selected for scanning. A minimum identity threshold of 95% and a minimum coverage threshold of 60% were applied; replicons meeting these criteria were considered indicative of the presence of plasmids belonging to corresponding incompatibility groups.

To further characterize plasmid features, plasmid typing was conducted using MOB-suite v3.0, providing multidimensional classification data based on Replicon-MOB-Mobility profiles. These results were visualized using a heatmap to illustrate plasmid variability among isolates. Plasmid homology analysis was conducted using BRIG v0.95, which generated circular alignment maps using representative plasmids as references, thereby revealing conserved core regions and divergent variable regions. Subsequently, the igraph package was employed to construct a plasmid network map, illustrating the cross-distribution of plasmids among patients and wards and indicating potential transmission chains.

Homology and phylogenetic analysis of isolates

Homology and phylogenetic relationships among the isolates in this study were analyzed based on SNPs. A core gene SNP-based phylogenetic tree was constructed using the maximum likelihood method to determine precise SNP distances among isolates. SNP calling was performed using Snippy v4.3, which aligned the sequences of interest to a reference dataset composed of 18 S. marcescens isolates from clinical blood cultures and catheters in China, obtained from GenBank. The core genome, including conserved sites, was extracted and filtered to remove recombination regions using Gubbins v2.3.4. SNP distances were calculated with snp-dists v0.7.0. Phylogenetic tree visualization and annotation were completed using the online tool ChiPlot.

Inference of isolate transmission events

Bayesian timed phylogenetic analysis was employed to infer the transmission dynamics of S. marcescens isolates within the healthcare setting. Time-calibrated phylogenetic trees were constructed using BEAST v2.4.4, with core SNP multi-locus alignments from 11 isolates serving as input and sampling dates used to define temporal nodes. The optimal substitution model HKY+G (Gamma categories = 4) was selected via jModelTest. A relaxed log-normal molecular clock model and a Bayesian skyline coalescent prior were applied. The Markov Chain Monte Carlo (MCMC) was run for 100,000,000 generations with sampling every 10,000 generations. Parameter convergence was assessed using Tracer v1.7. The maximum clade credibility tree was generated with TreeAnnotator v2.4.4, with node heights indicating the median divergence time.

At the level of transmission dynamics, TransPhylo v1.4 was applied to reconstruct transmission trees based on the time-scaled phylogeny, allowing inference of the number of unsampled cases and the structure of the transmission chains. The model was configured with gamma-distributed priors for generation time and sampling delay, and posterior sampling was used to infer transmission links between detected isolates and potential cryptic cases.

Results

General characteristics of S. marcescens-positive blood cultures

In October 2024, a total of 205 blood culture samples were submitted from the pediatric department of a tertiary hospital, among which 11 isolates of S. marcescens were identified, yielding a detection rate of 5.37%. This rate was significantly higher than the annual detection rate in 2023 (0.35%). Among the 11 cases, 2 were classified as hospital-acquired infections, while the remaining 9 were determined to be blood culture contaminations. The cases were distributed across eight pediatric wards, including Pediatric Neurology, Pediatric Hematology, Pediatric Cardionephrology and Endocrinology, Pediatric Respiratory Medicine Wards 1 and 2, the Department of Infectious Diseases, Pediatric Cardiothoracic Surgery, and Pediatric Orthopedics, indicating widespread distribution across multiple units (Table 1).

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics of 11 patients with Serratia marcescens positive blood cultures.

Population distribution

Among the 11 affected patients, 7 were male and 4 were female, with ages ranging from 2 months to 15 years. For the 9 cases identified as blood culture contaminations, 8 samples were collected on the day of admission and 1 prior to elective surgery, and none of these patients exhibited clinical signs of sepsis. In contrast, the two hospital-acquired infections both developed sepsis more than 48 h after admission, and one patient required transfer to the intensive care unit due to clinical deterioration. All patients received antimicrobial therapy according to institutional protocols, with definitive treatment primarily consisting of meropenem or imipenem based on susceptibility results.

Temporal and spatial distribution

Eleven cases of S. marcescens-positive blood cultures were primarily clustered into two outbreak periods during October 2024. The first cluster occurred within the 7-day period starting on the 7th, during which seven positive cases were reported. Significant overlap in hospitalization periods was observed, with three consecutive cases identified in the Department of Pediatric Neurology. Notably, Patient 1 and Patient 4 shared the same room with adjacent beds, indicating a potential risk of intra-ward transmission. In the following 2 days, additional cases were identified in the Pediatric Hematology Department and the second unit of the Pediatric Respiratory Department. The second cluster occurred between the 27th and 29th, with four cases reported during this interval. On the 27th alone, three cases were identified, including one confirmed hospital-acquired infection in the Department of Pediatric Cardiothoracic Surgery. On the 29th, one additional case was detected during a routine preoperative screening in the Pediatric Orthopedics Department. Overall, two distinct incidence peaks were observed in early and late October (Figure 1A).

FIGURE 1
Chart A illustrates infection and contamination cases over time with the investigation starting on October 10, 2024. Red bars represent infection cases, and blue bars represent contamination cases. Chart B displays isolation dates by hospital ward, using colored shapes to indicate departments and infection types, with red for hospital-acquired and blue for contamination. Department shapes include circles for Internal Medicine, triangles for Infectious Diseases, and squares for Surgery.

Figure 1. Temporal and spatial distribution of S. marcescens hospital-acquired sepsis and blood culture contamination. (A) Timeline of S. marcescens-positive blood cultures. Red bars indicate hospital-acquired sepsis cases; blue bars indicate contamination cases. (B) Spatiotemporal distribution by ward and date. Each symbol represents one isolate; different shapes correspond to department type (circle = internal medicine, triangle = infectious diseases, square = surgery), and colors denote infection classification (red = hospital-acquired, blue = contamination). The y-axis indicates ward location (building-floor).

In terms of spatial distribution, the 11 cases were predominantly located in two adjacent hospital buildings. Of these, 10 cases occurred on different floors of Building A, spanning multiple departments, including Pediatric Neurology, Pediatric Hematology, both units of Pediatric Respiratory Medicine, Pediatric Cardionephrology and Endocrinology, Infectious Diseases, and Pediatric Orthopedics. This indicates a scattered distribution across multiple wards and floors. The remaining case, a hospital-acquired infection, occurred in the Department of Pediatric Cardiothoracic Surgery located in Building B (Figure 1B). In summary, the cases exhibited two temporally distinct outbreak peaks and a spatial concentration across multiple floors of Building A, accompanied by inter-ward distribution.

Antimicrobial susceptibility results

The 11 isolated strains of S. marcescens exhibited highly consistent resistance profiles (Figure 2A). As expected for this species, all isolates exhibited intrinsic resistance to first- and second-generation β-lactam antibiotics (e.g., cefazolin, cefotetan, cefoxitin, cefuroxime), as well as to nitrofurantoin. However, they remained susceptible to 15 other antimicrobial agents, including five additional β-lactams, three carbapenems (imipenem, meropenem, doripenem), three aminoglycosides (amikacin, gentamicin, tobramycin), two quinolones (ciprofloxacin, levofloxacin), and two tetracyclines (minocycline, doxycycline).

FIGURE 2
Panel A shows a heatmap of antibiotic resistance for different strains (SM01 to SM12) against various antibiotics. Red indicates resistance, while blue indicates susceptibility. Panel B depicts a bar chart of minimal inhibitory concentrations (MIC) for six antibiotic classes, with nitrofurans showing the highest MIC around 500 micrograms per milliliter.

Figure 2. Antimicrobial resistance characteristics of clinical S. marcescens isolates. (A) Heatmap of antimicrobial susceptibility for 11 S. marcescens isolates. Blue indicates susceptibility; red indicates resistance. (B) Distribution of MICs across antimicrobial classes.

Further analysis of the minimum inhibitory concentration (MIC) distribution revealed uniformly elevated MIC values for first- and second-generation cephalosporins, consistent with the known intrinsic resistance of S. marcescens to these agents, whereas MIC values for third-generation cephalosporins and carbapenems remained low (Figure 2B). A subset of strains showed intermediate resistance tendencies to aminoglycosides, although their MIC values did not meet the threshold for formal resistance. According to the continuous antimicrobial resistance surveillance conducted in our hospital from 2019 to 2023, S. marcescens isolates have consistently exhibited intrinsic resistance to first- and second-generation cephalosporins, with no detection of resistance to third-generation cephalosporins or carbapenems and no evidence of ESBL- or carbapenemase-mediated mechanisms. The 11 isolates analyzed in this study showed the same susceptibility pattern, indicating that no new resistance mechanisms or increases in resistance levels have emerged in this event.

Furthermore, none of the 11 isolates displayed phenotypic features suggestive of AmpC β-lactamase derepression, and no characteristic patterns associated with AmpC induction or derepression were observed, supporting the absence of AmpC-related resistance mechanisms.

Resistance genes, virulence genes, and plasmid profiles

The resistance gene and virulence gene profiles of all 11 S. marcescens isolates were completely identical. A total of seven resistance genes were identified, including the aminoglycoside-modifying enzyme gene aac(6’)-Ic, the quinolone efflux pump gene oqxB, the streptomycin resistance gene srt-2, the tetracycline resistance gene tet(41), and three genes associated with multidrug efflux and regulation—mexI, CRP, and H-NS. Regarding virulence genes, all isolates harbored flgH, fliG, fliM, and cheY, genes involved in flagellar function and chemotaxis, suggesting a potential role in adhesion, motility, and host invasion. The resistance and virulence gene profiles of 28 S. marcescens isolates—comprising strains from blood cultures, catheters, and this study—are shown in Figure 3A.

FIGURE 3
Panel A shows a heatmap indicating sequence identity across different plasmid positions for samples SM01 to SM10, with colors ranging from purple to yellow. Panel B is a bar chart depicting the presence (1) or absence (0) of specific plasmid features, including replicon IncFII, MOBP, and mobility types for each sample. Panel C illustrates a network graph connecting different samples with nodes colored by hospital department, such as Pediatric Neurology and Infectious Diseases.

Figure 3. Plasmid alignment and typing analysis of S. marcescens isolates. (A) Homology heatmap based on reference plasmid sequences. Rows represent plasmid-carrying isolates; columns correspond to reference plasmid positions. The color gradient indicates sequence identity (0–1) across sliding windows. (B) Replicon-MOB-Mobility typing results. All plasmids carried IncFII replicons and were classified as MOBP-type. Among them, three isolates (SM01, SM03, SM07) harbored conjugative plasmids, while the remaining carried mobilizable plasmids. (C) Plasmid similarity network based on average nucleotide identity (ANI). Node colors indicate the wards/departments from which the isolates were obtained (Neurology, Hematology ward, Respiratory, Pediatric Nephrology, Infectious Diseases, and Pediatric Cardiothoracic Surgery), and node size is proportional to plasmid length. Edges represent plasmid similarity, with thicker lines indicating higher similarity; thick black edges denote conjugative plasmid-associated connections, whereas thin edges indicate general similarity links.

Among the nine S. marcescens isolates carrying plasmids, sequence alignment based on BRIG revealed high homology in core regions across plasmids, while variable regions exhibited deletions or divergence, indicating structural variation (Figure 3A). All plasmids carried the IncFII-type replicon and were classified as MOBP-type plasmids. Further plasmid mobility typing revealed that three isolates (SM01, SM03, SM07) harbored conjugative plasmids, while the remaining six possessed mobilizable plasmids (Figure 3B). These results suggest that although the plasmids exhibited overall structural conservation, certain variable regions differed, and the plasmids retained a limited but notable potential for horizontal transmission.

Plasmid similarity network analysis further revealed the relationships among plasmids from isolates obtained in different wards (Figure 3C). Isolates from the Department of Neurology (SM01, SM03, SM04) formed a highly homologous cluster, suggesting the presence of potential clonal transmission within this unit. SM02 originated from the Pediatric Nephrology Department, SM09 originated from the Pediatric Cardiothoracic Surgery Department, and SM08 from the Infectious Diseases Department. Isolates from the Hematology ward (SM05, SM07) also exhibited strong homology, forming a distinct subgroup. Isolates from the Respiratory Department (SM10) displayed plasmid-level similarity links to strains from both Neurology and Hematology, indicating the possible cross-ward distribution of homologous plasmids. The Infectious Diseases isolate (SM08) was also connected to multiple isolates in the network, further supporting potential inter-ward dissemination. Node size reflects plasmid length, and certain conjugative plasmids, such as those in SM01, SM03, and SM07, served as “hub nodes” connecting multiple isolates, highlighting their elevated potential for dissemination.

Homology analysis

The 11 isolates in this study exhibited significant genetic divergence from 17 reference S. marcescens strains derived from patient blood or catheter blood in Chinese healthcare institutions, forming an independent clonal cluster (Figure 4A). Based on core genome SNP alignment, two phylogenetic trees were constructed. As shown in Figure 4B, the 11 S. marcescens isolates demonstrated minimal genetic variation, with a maximum SNP distance of only 16. These isolates were divided into two branches. The hospital-acquired strains SM01 and SM09 showed no SNP differences and clustered within the same branch. The contamination-related isolates SM03 and SM05 formed a distinct clonal group, with a pairwise core-genome SNP distance of only 1 SNP between them (Figure 4B). Notably, the contamination strains SM06 and SM07 shared identical SNP profiles with the hospital-acquired isolate SM09. These findings reveal a high degree of homology and a closely linked transmission chain among S. marcescens strains involved in this nosocomial outbreak.

FIGURE 4
Two phylogenetic trees labeled A and B with heatmaps. Tree A shows the relationship among bacterial strains, resistance genes, and geographical origins. Strains in the study and reference strains are highlighted. The heatmap indicates the presence of resistance and virulence genes. Tree B shows the genetic relationships between strains with a scale for single nucleotide polymorphisms (SNPs). The accompanying heatmap shows SNP differences between strains, with varying shades of blue. Both sections categorize strains as contamination or hospital-acquired.

Figure 4. Core genome phylogenetic trees, resistance/virulence gene profiles, and SNP distance matrix of S. marcescens isolates. (A) Maximum-likelihood phylogenetic tree based on core genome alignment, including 11 S. marcescens isolates from this study (orange) and 17 reference strains (green). The presence of resistance genes (red) and virulence genes (blue) is shown as a heatmap to the right of the phylogenetic tree. The far-right column indicates isolate classification (hospital-acquired vs. contamination) and the geographic origin of reference strains. (B) SNP-based phylogenetic tree and pairwise SNP distance heatmap for the 11 isolates in this study. Color intensity reflects the number of SNP differences, as indicated by the scale bar. Hospital-acquired isolates are labeled; the others represent contamination cases. Numbers within the heatmap cells denote pairwise core-genome SNP distances between isolates, for example, SM03 and SM05 differ by only 1 SNP.

Analysis of isolate transmission events

Bayesian-timed phylogenetic trees revealed that the most recent common ancestor of the S. marcescens isolates in this study emerged in late December 2023 (2023.96). This finding indicates that the source of infection had been circulating widely in the hospital environment for at least 10 months before detection in October 2024, without effective elimination. This prolonged presence likely facilitated widespread intra-hospital transmission, ultimately leading to the observed cluster of cases. The 11 isolates diverged into three clades across three time points in 2024. The hospital-acquired strains SM01 and SM09 were located on the same divergence branch at a specific time point, despite originating from different hospital buildings. This suggests the possibility of a shared transmission chain. In contrast, blood culture contamination isolates were distributed across all three divergence points (Figure 5A).

FIGURE 5
Panel A shows a phylogenetic tree with various samples labeled SM01 to SM12, including dates and evolutionary distances. Panel B displays a minimum spanning tree with samples connected by lines, showing relationship clusters and relevant dates.

Figure 5. Phylogenetic and transmission inference of a S. marcescens bloodstream culture cluster in pediatric wards. (A) Time-scaled phylogenetic tree constructed using BEAST. The x-axis represents nucleotide differences (SNP distance), and the y-axis denotes the inferred timeline. The most recent common ancestor (TMRCA) was estimated to have emerged in late December 2023 (approximately 2023.96). Hospital-acquired isolates SM01 and SM09 are highlighted in red and clustered on the same divergence branch. (B) Transmission tree inferred by TransPhylo based on the BEAST-derived phylogeny. Solid circles represent sampled isolates; hollow circles indicate unsampled individuals inferred by the model (U1-U4).

Using TransPhylo to infer transmission events from the BEAST-derived timed phylogeny, approximately four unsampled individuals (represented as hollow nodes; 95% CrI = 1–10) were predicted in addition to the 11 sampled isolates (solid nodes) (Figure 5B). These unsampled nodes, dated between January and October 2024, served as intermediaries connecting various sampled strains. For instance, node U4 was located between SM01 and SM09, forming the path “SM01 → U4 → SM09,” which explains how these genetically related isolates could belong to the same transmission chain despite differing in collection time and hospital ward. Other unsampled nodes (e.g., U1-U3) appeared earlier, from early 2024 to the summer months, suggesting the presence of cryptic transmission before the observed cluster.

Environmental sampling and the effectiveness of bundled infection control measures

Environmental sampling and comprehensive infection control measures were implemented concurrently with the epidemiological investigation. A total of 125 specimens were collected, including 72 hand swabs (11 from patients, 18 from physicians, 27 from nurses, 11 from ancillary staff, and 5 from specimen transport personnel), 30 samples from patients’ surroundings, 15 from ward sinks, and 8 from disinfectants in use. S. marcescens was not identified in any of these samples. Following a multidisciplinary meeting, bundled infection control measures were introduced, alongside interventions aimed at reducing blood culture contamination. On-site investigation revealed that pediatric ward nurses were using povidone-iodine for skin disinfection before blood collection, but did not allow sufficient drying time before puncture. Although S. marcescens was not detected in sink samples, a high burden of environmental microorganisms was observed. The implemented measures included: (1) replacing povidone-iodine with a faster-drying ≥ 2 g/L chlorhexidine-ethanol (70%) solution for skin disinfection; (2) replacing sinks, faucets, and pipelines in all eight wards with S. marcescens-positive blood cultures; (3) ensuring that healthcare workers, patients, and caregivers used alcohol-based hand sanitizers, with reinforced education, training, and supervision on hand hygiene; (4) strengthening hand hygiene monitoring for patients, caregivers, and medical staff; (5) enhancing surveillance of S. marcescens bloodstream infections; (6) isolating S. marcescens blood culture-positive patients and their caregivers in single rooms; (7) reviewing and revising the clinical specimen submission process; (8) increasing scrutiny of blood culture collection procedures by nursing staff; (9) enforcing stricter management of ward consultation physicians; and (10) upgrading terminal disinfection equipment in the wards. During the 1-month follow-up period after these interventions, no additional cases of S. marcescens hospital-acquired bloodstream infections or blood culture contamination were reported in the pediatric wards.

Discussion

S. .marcescens has become an important pathogen in recent healthcare-associated infections, with multiple outbreaks reported in neonatal and intensive care settings (Aracil-Gisbert et al., 2024; Zhu et al., 2025). However, small-scale and cryptic transmission events occurring in general pediatric wards remain insufficiently characterized. In this study, we documented 11 blood culture–positive S. marcescens events across eight pediatric units in a tertiary hospital. Two patients developed confirmed bloodstream infections, whereas the remaining cases were initially judged as contamination; nevertheless, the isolates showed highly homogeneous genomic characteristics, indicating an underestimated clonal transmission chain. The detection rate in October 2024 (5.37%) was markedly higher than that of 2023 (0.35%), demonstrating a clear clustering pattern. Given the immature immunity of pediatric patients and the high risk of contamination during blood collection, general pediatric wards represent susceptible and easily overlooked environments. These findings underscore the need to strengthen infection surveillance in non-critical care settings.

The high number of S. marcescens–positive blood cultures labeled as contaminants, with only two true bloodstream infections, aligns with NICU outbreak reports where infant gut/mucosal colonization is common but invasive disease is uncommon (Montagnani et al., 2015; Mathew et al., 2015; Bourdin et al., 2023). In our cluster, the near-identical genomic and resistance profiles of all 11 isolates indicate a single hospital-adapted lineage colonizing multiple children; under such colonization pressure, difficult pediatric sampling and suboptimal povidone-iodine antisepsis likely caused inoculation into culture bottles, and the disappearance of S. marcescens–positive cultures after introducing chlorhexidine–alcohol skin preparation supports contamination related to colonization and technique rather than unrecognized bacteremia (Min et al., 2014; Marlowe et al., 2010).

WGS has been repeatedly shown to rapidly confirm clonality, reconstruct transmission chains, and guide infection-control interventions in NICU and ICU outbreaks (Henares et al., 2025; Merla et al., 2024). Large ICU and nationwide genomic investigations also demonstrate that S. marcescens can persist in hospital environments for extended periods, often without a definable point source (Zhu et al., 2025; Aracil-Gisbert et al., 2024; Taxt et al., 2025). These findings align with our BEAST/TransPhylo inference that the transmission chain likely originated in late 2023. A key distinction of this study is that the event occurred in general pediatric wards, where most cases appeared as blood culture contamination but were genomically indistinguishable. This suggests that true transmission in low-acuity settings may be more easily missed. Furthermore, compared with recent genomic studies highlighting the diversity of AmpC activity and resistance evolution (Chen et al., 2025), our isolates carried multiple resistance genes but lacked AmpC derepression mutations and remained susceptible to third-generation cephalosporins and carbapenems. This indicates an early-stage lineage with strong transmissibility but not yet advanced resistance, contributing new insights into S. marcescens transmission in non-critical pediatric units.

From a mechanistic standpoint, S. marcescens can persist in wet hospital environments such as sinks and drains, which serve as important environmental reservoirs (Aracil-Gisbert et al., 2024). Outbreaks are often associated with AmpC derepression or plasmid-mediated resistance (Chen et al., 2025). Although our isolates carried aac(6’)-Ic, oqxB, tet(41) and other resistance genes, no AmpC derepression mutations were detected, and carbapenem susceptibility was preserved—consistent with a “high transmission but low resistance” early-phase lineage. Although S. marcescens possesses an inducible AmpC system, WGS of our isolates revealed no ampC, ampD, or ampR mutations associated with derepression. Previous work (Pérez-Gallego et al., 2016) shows that extreme AmpC overexpression, more than 1,000-fold above wild-type levels, occurs only in amidase-deficient backgrounds and can affect resistance and fitness. No such genetic features were identified in our isolates, all of which remained susceptible to third-generation cephalosporins, thereby excluding the possibility of AmpC derepression. Additionally, the widespread presence of motility and adhesion–related virulence genes such as cheY and fliM supports enhanced environmental survival and intra-ward dissemination (Prado et al., 2022; Valiatti et al., 2023). Together, these features explain the “low virulence but high colonization and spread” phenotype observed in this event.

Although nine cases were ultimately classified as blood culture contamination, their genomic similarity to the two bloodstream infection isolates indicates that traditional contamination assessments based solely on clinical presentation and procedural evaluation may be insufficient. Previous studies have similarly reported that early, low-symptom infections during cluster events can be misclassified as contamination (Klucher et al., 2021). Integrating WGS data into contamination adjudication may therefore improve diagnostic accuracy, avoid missed infections, and prevent overlooking ongoing transmission chains.

This study has several limitations. First, the genomic analysis included only 11 S. marcescens isolates, which, while sufficient to define the clonal structure of this outbreak, is smaller than sample sizes typically used in large-scale WGS epidemiological investigations. Thus, the reconstructed transmission network may not capture all routes or unsampled carriers. Second, all isolates were obtained from a single pediatric center during a limited outbreak period, which may restrict the generalizability of the findings. Third, environmental sampling did not yield positive isolates, preventing assessment of whether environmental persistence contributed to the prolonged circulation suggested by the Bayesian analysis. Fourth, we did not systematically assess gastrointestinal colonization (for example, through stool or rectal screening cultures) in the nine patients with blood-culture contamination, and therefore could not evaluate the gut as a potential reservoir in this cluster. Finally, integrating WGS with broader longitudinal sampling, metagenomic surveillance, and real-time infection control assessments would further improve outbreak interpretation. Future studies with larger, multi-center datasets and continuous genomic monitoring will be crucial to refine transmission dynamics and support stronger infection control strategies in pediatric wards.

Conclusion

This study reveals that a genomically uniform S. marcescens lineage can establish sustained transmission even within general pediatric wards, environments not traditionally viewed as conducive to long-term nosocomial persistence. Genomic reconstruction, supported by Bayesian temporal and transmission analyses, indicates that the cluster represented a prolonged and only partially detected transmission process that routine surveillance did not capture (Figure 6). The conserved resistance, virulence, and plasmid architecture suggest early plasmid-stabilized dissemination preceding the emergence of more advanced resistance traits. These findings show that low-acuity pediatric settings can maintain silent propagation of hospital-adapted organisms and illustrate the limitations of relying solely on clinical presentation to distinguish contamination from infection. More broadly, the work underscores the value of integrating real-time genomic epidemiology into pediatric infection-control frameworks to detect emerging lineages at an early stage and to intervene before transmission becomes entrenched.

FIGURE 6
Flowchart illustrating infection control in a hospital setting. It shows case distribution across eight wards, with eleven cases. Antimicrobial susceptibility and distribution are shown, highlighting concerns with cephalosporins and resistance genes in IncFII plasmid. Homology and transmission chain analysis indicate low SNP variation and hidden cases. The intervention section illustrates hygiene measures, leading to zero new cases.

Figure 6. Genomic epidemiology and infection control strategies for a S. marcescens outbreak in a pediatric hospital.

Data availability statement

The original contributions presented in this study are included in this article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were reviewed and approved by the Ethics Committee of the Third Affiliated Hospital of Zhengzhou University (Approval No. 2025-187-01). The requirement for written informed consent was waived because this study was a retrospective observational study using existing clinical and microbiological data, in accordance with institutional and local ethical regulations.

Author contributions

WZ: Conceptualization, Data curation, Formal analysis, Funding acquisition, Writing – original draft. CM: Investigation, Methodology, Project administration, Resources, Writing – original draft. FX: Software, Supervision, Validation, Visualization, Writing – review & editing. CZ: Conceptualization, Methodology, Project administration, Visualization, Writing – review & editing. LW: Formal analysis, Methodology, Software, Visualization, Writing – original draft.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the Science and Technology Research Program of Henan Province (No. 252102311086).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1718340/full#supplementary-material

Supplementary Figure 1 | Study flowchart.

Abbreviations

ANI, Average Nucleotide Identity; CARD, Comprehensive Antibiotic Resistance Database; MIC, Minimum Inhibitory Concentration; MCMC, Markov Chain Monte Carlo; SNP, Single Nucleotide Polymorphism; VFDB, Virulence Factor Database; WGS, Whole Genome Sequencing.

Footnotes

References

Alam, M. T., Read, T. D., Boyle-Vavra, S., Miller, L. G., Eells, S. J., et al. (2015). Transmission and microevolution of USA300 MRSA in U.S. households: Evidence from whole-genome sequencing. mBio 6:e00054. doi: 10.1128/mbio.00054-15

PubMed Abstract | Crossref Full Text | Google Scholar

Aracil-Gisbert, S., Fernández-De-Bobadilla, M. D., Guerra-Pinto, N., Serrano-Calleja, S., Pérez-Cobas, A. E., Soriano, C., et al. (2024). The ICU environment contributes to the endemicity of the “Serratia marcescens complex” in the hospital setting. mBio 15:e0305423. doi: 10.1128/mbio.03054-23

PubMed Abstract | Crossref Full Text | Google Scholar

Arslan, U., Erayman, I., Kirdar, S., Yuksekkaya, S., Cimen, O., Tuncer, I., et al. (2010). Serratia marcescens sepsis outbreak in a neonatal intensive care unit. Pediatrics Int. 52, 208–212. doi: 10.1111/j.1442-200x.2009.02934.x

PubMed Abstract | Crossref Full Text | Google Scholar

Bechmann, L., Böttger, R., Baier, C., Tersteegen, A., Bauer, K., Kaasch, A. J., et al. (2022). Serratia marcescens outbreak in a neonatal intensive care unit associated with contaminated donor milk. Infect. Control Hosp. Epidemiol. 44, 891–897. doi: 10.1017/ice.2022.187

PubMed Abstract | Crossref Full Text | Google Scholar

Bourdin, T., Benoit, M. -È,Monnier, A., Bédard, E., Prévost, M., Charron, D., et al. (2023). Serratia marcescens colonization in a neonatal intensive care unit has multiple sources, with sink drains as a major reservoir. Appl. Environ. Microbiol. 89:e0010523. doi: 10.1128/aem.00105-23

PubMed Abstract | Crossref Full Text | Google Scholar

Brlek, P., Bulić, L., Bračić, M., Projić, P., Škaro, V., Shah, N., et al. (2024). Implementing Whole genome sequencing (WGS) in clinical practice: Advantages, challenges, and future perspectives. Cells 13:504. doi: 10.3390/cells13060504

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, W., Li, Z., Zheng, J., Wang, D., Gao, S., Zhou, W., et al. (2025). A comprehensive analysis of the epidemiological and genomic characteristics of global Serratia Marcescens. J. Glob. Antimicrobial Resist. 40, 81–89. doi: 10.1016/j.jgar.2024.11.013

PubMed Abstract | Crossref Full Text | Google Scholar

Cristina, M. L., Sartini, M., and Spagnolo, A. M. (2019). Serratia marcescens infections in Neonatal intensive care units (NICUs). IJERPH 16:610. doi: 10.3390/ijerph16040610

PubMed Abstract | Crossref Full Text | Google Scholar

Crivaro, V., Bagattini, M., Salza, M. F., Raimondi, F., Rossano, F., Triassi, M., et al. (2007). Risk factors for extended-spectrum β-lactamase-producing Serratia marcescens and Klebsiella pneumoniae acquisition in a neonatal intensive care unit. J. Hosp. Infect. 67, 135–141. doi: 10.1016/j.jhin.2007.07.026

PubMed Abstract | Crossref Full Text | Google Scholar

De Silva, D., Peters, J., Cole, K., Cole, M. J., Cresswell, F., Dean, G., et al. (2016). Whole-genome sequencing to determine transmission of Neisseria gonorrhoeae: An observational study. Lancet Infect. Dis. 16, 1295–1303. doi: 10.1016/s1473-3099(16)30157-8

PubMed Abstract | Crossref Full Text | Google Scholar

Didelot, X., Gardy, J., and Colijn, C. (2014). Bayesian inference of infectious disease transmission from whole-genome sequence data. Mol. Biol. Evol. 31, 1869–1879. doi: 10.1093/molbev/msu121

PubMed Abstract | Crossref Full Text | Google Scholar

Fitzpatrick, M. A., Ozer, E. A., and Hauser, A. R. (2016). Utility of whole-genome sequencing in characterizing acinetobacter epidemiology and analyzing hospital outbreaks. J. Clin. Microbiol. 54, 593–612. doi: 10.1128/jcm.01818-15

PubMed Abstract | Crossref Full Text | Google Scholar

Ghaith, D. M., Zafer, M. M., Ismail, D. K., Al-Agamy, M. H., Bohol, M. F. F., Al-Qahtani, A., et al. (2018). First reported nosocomial outbreak of Serratia marcescens harboring bla IMP–4 and bla VIM–2 in a neonatal intensive care unit in Cairo, Egypt. Infect Drug Resist. 11, 2211–2217. doi: 10.2147/idr.s174869

PubMed Abstract | Crossref Full Text | Google Scholar

Harch, S. A. J., Jenkins, F., Farhat, R., and van Hal, S. J. (2025). Complexities in species identification for Serratia marcescens complex for the modern microbiology laboratory. Microbiol. Spectr. 13:e01361-24. doi: 10.1128/spectrum.01361-24

PubMed Abstract | Crossref Full Text | Google Scholar

Henares, D., Cubero, M., Martinez-de-Albeniz, I., Arranz, A., Rocafort, M., Brotons, P., et al. (2025). Rapid identification of a Serratia marcescens outbreak in a neonatal intensive care unit by third-generation long-read nanopore sequencing. Antimicrob. Resist. Infect. Control 14:63. doi: 10.1186/s13756-025-01582-x

PubMed Abstract | Crossref Full Text | Google Scholar

Klucher, J. M., Davis, K., Lakkad, M., Painter, J. T., and Dare, R. K. (2021). Risk factors and clinical outcomes associated with blood culture contamination. Infect. Control Hosp. Epidemiol. 43, 291–297. doi: 10.1017/ice.2021.111

PubMed Abstract | Crossref Full Text | Google Scholar

Marlowe, L., Mistry, R. D., Coffin, S., Leckerman, K. H., McGowan, K. L., Dai, D., et al. (2010). Blood culture contamination rates after skin antisepsis with chlorhexidine gluconate versus povidone-iodine in a pediatric emergency department. Infect. Control Hosp. Epidemiol. 31, 171–176. doi: 10.1086/650201

PubMed Abstract | Crossref Full Text | Google Scholar

Martineau, C., Li, X., Lalancette, C., Perreault, T., Fournier, E., Tremblay, J., et al. (2018). Serratia marcescens Outbreak in a neonatal intensive care unit: New insights from next-generation sequencing applications. J. Clin. Microbiol. 56:e00235-18. doi: 10.1128/jcm.00235-18

PubMed Abstract | Crossref Full Text | Google Scholar

Mathew, N., Adhisivam, B., and Bhat, B. V. (2015). Clinical profile and outcome of serratia infection among neonates. Indian J. Pediatr. 83, 726–727. doi: 10.1007/s12098-015-1928-7

PubMed Abstract | Crossref Full Text | Google Scholar

Merla, C., Ramus, M., Kuka, A., Mileto, I., Gaiarsa, S., Di Comite, A., et al. (2024). Impact of whole genome sequencing to investigate transmission of Serratia marcescens in neonatal intensive care unit. J. Infect. Public Health 17, 918–921. doi: 10.1016/j.jiph.2024.03.025

PubMed Abstract | Crossref Full Text | Google Scholar

Min, H., Park, C. S., Kim, D. S., and Kim, K. H. (2014). Blood culture contamination in hospitalized pediatric patients: A single institution experience. Korean J. Pediatr. 57:178. doi: 10.3345/kjp.2014.57.4.178

PubMed Abstract | Crossref Full Text | Google Scholar

Mitchell, B. I., Kling, K., Bolon, M. K., Rathod, S. N., Malczynski, M., Ruiz, J., et al. (2024). Identifying Candida auris transmission in a hospital outbreak investigation using whole-genome sequencing and SNP phylogenetic analysis. J. Clin. Microbiol. 62:e0068024. doi: 10.1128/jcm.00680-24

PubMed Abstract | Crossref Full Text | Google Scholar

Montagnani, C., Cocchi, P., Lega, L., Campana, S., Biermann, K. P., Braggion, C., et al. (2015). Serratia marcescens outbreak in a neonatal intensive care unit: Crucial role of implementing hand hygiene among external consultants. BMC Infect. Dis. 15:11. doi: 10.1186/s12879-014-0734-6

PubMed Abstract | Crossref Full Text | Google Scholar

Mosalli, R., Alqarni, S. A., Khayyat, W. W., Alsaidi, S. T., Almatrafi, A. S., Bawakid, A. S., et al. (2022). Respiratory syncytial virus nosocomial outbreak in neonatal intensive care: A review of the incidence, management, and outcomes. Am. J. Infect. Control 50, 801–808. doi: 10.1016/j.ajic.2021.10.027

PubMed Abstract | Crossref Full Text | Google Scholar

Parker, L. A., Pruitt, J., Monk, A., Lambert, M. T., Lorca, G. L., and Neu, J. (2023). Oral care in critically ill infants and the potential effect on infant health: An integrative review. Crit. Care Nurse 43, 39–50. doi: 10.4037/ccn2023902

PubMed Abstract | Crossref Full Text | Google Scholar

Patel, M., and Kaufman, D. A. (2015). Anti-lipoteichoic acid monoclonal antibody (pagibaximab) studies for the prevention of staphylococcal bloodstream infections in preterm infants. Expert Opin. Biol. Therapy 15, 595–600. doi: 10.1517/14712598.2015.1019857

PubMed Abstract | Crossref Full Text | Google Scholar

Pérez-Gallego, M., Torrens, G., Castillo-Vera, J., Moya, B., Zamorano, L., Cabot, G., et al. (2016). Impact of AmpC derepression on fitness and virulence: The mechanism or the pathway? mBio 7:e01783-16. doi: 10.1128/mbio.01783-16

PubMed Abstract | Crossref Full Text | Google Scholar

Pimentel, de Araujo, F., Tinelli, M., Battisti, A., Ercoli, A., Anesi, A., et al. (2017). An outbreak of skin infections in neonates due to a Staphylococcus aureus strain producing the exfoliative toxin A. Infection 46, 49–54. doi: 10.1007/s15010-017-1084-2

PubMed Abstract | Crossref Full Text | Google Scholar

Prado, G., Mendes, E. T., Martins, R. C. R., Perdigão-Neto, L. V., Freire, M. P., Marchi, A. P., et al. (2022). Phenotypic and genotypic characteristics of a carbapenem-resistant Serratia marcescens cohort and outbreak: Describing an opportunistic pathogen. Int. J. Antimicrob. Agents 59:106463. doi: 10.1016/j.ijantimicag.2021.106463

PubMed Abstract | Crossref Full Text | Google Scholar

Reddy, K., Bekker, A., Whitelaw, A. C., Esterhuizen, T. M., and Dramowski, A. (2021). A retrospective analysis of pathogen profile, antimicrobial resistance and mortality in neonatal hospital-acquired bloodstream infections from 2009–2018 at Tygerberg Hospital, South Africa. PLoS One 16:e0245089. doi: 10.1371/journal.pone.0245089

PubMed Abstract | Crossref Full Text | Google Scholar

Shaik Ismail, B., Toh, H. X., Seah, J. H., Tan, K. Y., Lee, L. C., Tay, Y. Y., et al. (2025). Serratia marcescens outbreak at a neonatal intensive care unit in an acute care tertiary hospital in Singapore. J. Hosp. Infect. 156, 21–25. doi: 10.1016/j.jhin.2024.10.002

PubMed Abstract | Crossref Full Text | Google Scholar

Shane, A. L., Sánchez, P. J., and Stoll, B. J. (2017). Neonatal sepsis. Lancet 390, 1770–1780. doi: 10.1016/s0140-6736(17)31002-4

PubMed Abstract | Crossref Full Text | Google Scholar

Sundermann, A. J., Chen, J., Miller, J. K., Martin, E. M., Snyder, G. M., Van Tyne, D., et al. (2022). Whole-genome sequencing surveillance and machine learning for healthcare outbreak detection and investigation: A systematic review and summary. Antimicrob. Steward Healthc. Epidemiol. 2:e91. doi: 10.1017/ash.2021.241

PubMed Abstract | Crossref Full Text | Google Scholar

Taxt, A. M., Eldholm, V., Kols, N. I., Haugan, M. S., Raffelsberger, N., Asfeldt, A. M., et al. (2025). A national outbreak of Serratia marcescens complex: Investigation reveals genomic population structure but no source, Norway, June 2021 to February 2023. Eurosurveillance 30:2400291. doi: 10.2807/1560-7917.es.2025.30.5.2400291

PubMed Abstract | Crossref Full Text | Google Scholar

Tsagris, V., Nika, A., Kyriakou, D., Kapetanakis, I., Harahousou, E., Stripeli, F., et al. (2012). Influenza A/H1N1/2009 outbreak in a neonatal intensive care unit. J. Hosp. Infect. 81, 36–40. doi: 10.1016/j.jhin.2012.02.009

PubMed Abstract | Crossref Full Text | Google Scholar

Valiatti, T. B., Bessa-Neto, F. O., Santos, F. F., Silva, R. G. B., Veiga, R., Cassu-Corsi, D., et al. (2023). Clonal dissemination of highly virulent Serratia marcescens strains producing KPC-2 in food-producing animals. One Health 17:100591. doi: 10.1016/j.onehlt.2023.100591

PubMed Abstract | Crossref Full Text | Google Scholar

Williams, D. J., Grimont, P. A. D., Cazares, A., Grimont, F., Ageron, E., Pettigrew, K. A., et al. (2022). The genus Serratia revisited by genomics. Nat. Commun. 13:5195. doi: 10.1038/s41467-022-32929-2

PubMed Abstract | Crossref Full Text | Google Scholar

Zackeroff, S., Nash, D., McDermott, K., Miller, R. R., and Pasquini, G. (2024). Central lines, aseptic batching services, and infection rates: A pharmacy-led initiative of intravenous tube priming within a NICU. Am. J. Health-Syst. Pharm. 81, e777–e782. doi: 10.1093/ajhp/zxae197

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, W., Chen, X., Shen, H., Wei, M., Yang, C., and Gu, L. (2025). Genomic diversity, antimicrobial resistance and dissemination of Serratia marcescens complex in patients admitted to ICUs. Front. Cell. Infect. Microbiol. 15:1672468. doi: 10.3389/fcimb.2025.1672468

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: infection control, infection control (ICP) practitioners, sepsis, Serratia marcescens, whole genome sequencing

Citation: Zhang W, Ma C, Xu F, Zhao C and Wang L (2026) Whole-genome sequencing-based source tracing and infection control of Serratia marcescens blood culture events in pediatric patients. Front. Microbiol. 16:1718340. doi: 10.3389/fmicb.2025.1718340

Received: 06 October 2025; Revised: 23 December 2025; Accepted: 23 December 2025;
Published: 16 January 2026.

Edited by:

João Pedro Rueda Furlan, Federal University of Paraíba, Brazil

Reviewed by:

Cristina Merla, San Matteo Hospital Foundation (IRCCS), Italy
Sonia Aracil-Gisbert, University of Alcalá, Spain
Blanca Pérez Viso, Ramón y Cajal University Hospital, Spain

Copyright © 2026 Zhang, Ma, Xu, Zhao and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ling Wang, d2FuZ2xpbmcwMTIyQHNpbmEuY29t; Wang Zhang, end6ZHNmeUAxNjMuY29t

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.