Preplanned Studies: Hidden Spread of RareSalmonellaSerovars Isolated from Healthy Individuals — Yulin City, Guangxi Zhuang Autonomous Region, China, 2013–2022

Gastrointestinal infections, predominantly foodborne, are among the most common illnesses worldwide, affecting an estimated 550 million people annually, including 220 million children younger than five years of age (1).Salmonella entericais a major pathogen implicated in these infections.Salmonellacomprises more than 2,500 serovars.SalmonellaTyphimurium andSalmonellaEnteritidis are frequently associated with diarrheal illnesses (2). Recent reports indicate that outbreaks can also be attributed to rareSalmonellaserovars. For example,S.Ealing,S.Welikade,S.Isangi, andS.Coeln have been linked to outbreaks in the UK, France, China, and Norway, respectively (3–6).

One Health is an integrated approach to balancing and optimizing the health of people, animals, and ecosystems. It leverages the close, interdependent links among these fields to create new surveillance and disease control methods (7). Healthy individuals, particularly food handlers, could be potential sources of hiddenSalmonelladissemination. Therefore, monitoringSalmonellaprevalence among healthy individuals is essential for comprehensively understanding and managing the associated risks.

Monitoring the prevalence of multidrug-resistant (MDR)Salmonellaisolates is a common method for assessing the risk of salmonellosis. This approach facilitates timely detection and intervention in cases of salmonellosis, thereby limiting further pathogen dissemination. Whole-genome sequencing (WGS) is widely used to predict and identify serovars, virulence factors, and antimicrobial resistance genes (ARGs) inSalmonellasurveillance.

Using resistance profiling and genomic analysis ofSalmonellaisolates collected from healthy individuals during a decade-long, large-scale, laboratory-based investigation, we identified five rareSalmonellaserovars reported for the first time in China and revealed the prevalence of MDR isolates and the hidden dissemination of these serovars.

Between 2013 and 2020, fecal samples were systematically collected from individuals undergoing occupational health examinations in Yulin City, Guangxi Zhuang Autonomous Region, China.Salmonellaspp. were cultured from these samples according to a previously established protocol (8). Serovars were determined using slide agglutination and by predicting genome sequences usingSalmonellaIn SilicoTyping Resource (SISTR). A systematic literature review using keyword searches for “Salmonella” and “serotype” was performed. Corresponding genomes were retrieved from the EnteroBase and GenBank databases to verify the novelty of the identified serovars within China.

Two commercial antimicrobial susceptibility testing panels (BD Phoenix™ NMIC-413 and BD Phoenix™ RUONMIC-801) were used with the microbroth dilution method, following the manufacturer’s instructions. For clarity, only the common panel of 11 antibiotics is presented, namely: Amikacin (AMK), Ceftazidime (CAZ), Chloramphenicol (CHL), Ciprofloxacin (CIP), Colistin (COL), Ertapenem (ETP), Meropenem (MEM), Piperacillin-Tazobacta (TZP), Tetracycline (TET), Tigecycline (TIG), Trimethoprim-Sulfamethoxazole (SXT). Results were interpreted in accordance with the Clinical and Laboratory Standards Institute (CLSI) 2022 guidelines.Escherichia coliATCC 25922 served as the quality control strain.

Genomic DNA was extracted from 15 bacterial isolates using the Wizard Genomic DNA Purification Kit (Promega) after incubation. Paired-end libraries with insert sizes of 300-500 bp were constructed and sequenced on the MGI200 platform according to the manufacturer’s protocol. The sequencing reads for each isolate were assembled using SPAdes v3.5.0. ARGs, plasmid profiles, and virulence factors were identified using Abricate software and databases, including CARD, PlasmidFinder, and VFDB.

To trace the origins of theS.Ouakam andS.Ealing serotype epidemics, 79S.Ouakam and 114S.Ealing isolates were sourced from the EnteroBase database. These isolates were analyzed using single-nucleotide polymorphisms (SNPs) and hierarchical clustering of core genome multilocus sequence typing (HierCC). The analyses were performed using the online analytical pipeline within EnteroBase. The resulting data were used to generate a phylogenetic tree and visualized using the iTOL v6 platform.

A total of 7,044Salmonellaisolates were recovered from 372,708 individual samples in a retrospective investigation conducted from 2013 to 2022, with an isolation rate of 1.89% among healthy individuals. Among these, 119 distinctSalmonella entericaserovars were identified in healthy individuals from Yulin. Utilizing an extensive literature review, we report the first identification of five rareSalmonellaserovars in China. A total of 14 isolates of rare serovars were sourced from individuals ranging in age from 17 to 60 years, comprising 86% (12/14) males and 14% (2/14) females. Notably, 93% (13/14) of the individuals were engaged in food handling, while the remaining 7% (1/14) were not. The isolates comprised five rare serovars:S.Welikade (n=1),S.Mountpleasant (n=1),S.Coeln (n=1),S.Ealing (n=2), andS.Ouakam (n=9).

S.Welikade,S.Mountpleasant, andS.Ealing isolates were susceptible to all 11 antibiotics tested. In contrast, theS.Coeln isolate and nineS.Ouakam isolates were MDR. Both exhibited resistance to CHL, SXT, and TET. Notably, theS.Coeln isolate was resistant to TIG, an antibiotic typically reserved for severe Gram-negative infections. Meanwhile, all nineS.Ouakam isolates were resistant to CAZ and CIP (Table 1).

Analysis of the 15 isolates’ genomes revealed five distinct plasmid types. The twoS.Ealing isolates lacked plasmids, whereas theS.Mountpleasant isolate harbored an IncQ1 plasmid. TheS.Welikade isolate contained both Col440I and p0111 plasmids. AllS.Ouakam andS.Coeln isolates carried the IncHI2/IncHI2A plasmid (Figure 1).

Figure 1.

Distribution of plasmids and antimicrobial resistance genes (ARGs) among the studied isolates.

TheS.Ouakam isolates possessed diverse ARGs, conferring resistance to aminoglycosides, sulfonamides, lincosamides, tetracyclines, beta-lactams, chloramphenicol, quinolones, and rifamycins. TheS.Coeln isolates exhibited resistance mediated by a broader spectrum of 14 ARGs, attributed to the presence of IncHI2/IncHI2A plasmids. Without these plasmids, theS.Welikade,S.Mountpleasant, andS.Ealing isolates displayed a more restricted resistance profile, predominantly featuring a single ARG,aac(6')-Iaa(Figure 1).

Our investigation revealed 147 distinct virulence genes within theSalmonellaisolates, highlighting multiple mechanisms underlying virulence and pathogenicity. ThecdtBgene for typhoid toxin production was found inS.Welikade, whileastA(EAST1 toxin) was detected in the Sa17631 isolate ofS.Ouakam. Furthermore, all rare serotype isolates harbored gene clusters involved in iron uptake (entA,entB,entC,entD,entE,entF,entS) and themig-14gene, associated with resistance to antimicrobial peptides (Figure 1).

We analyzed a dataset of 34,408 SNPs to elucidate the phylogenetic relationships among 88S.Ouakam genomes. The resulting phylogenetic tree revealed two distinct groups (Figure 2A). To facilitate further analysis, we reconstructed the phylogenetic tree specifically for the isolates of group 2 (Figure 2B). NineS.Ouakam isolates, derived from healthy individuals, formed a cluster within clade α, demonstrating SNP distances ranging from 0 to 58. To define outbreak clusters accurately, we employed the internationally recognized HierCC HC5 standard, which is widely used in infectious disease outbreaks (9). SevenS.Ouakam isolates classified under the same HC5 cluster (type 428586) were obtained from individuals employed by the same company, indicating the potential for hidden transmission ofSalmonellaspp.among this population.

Figure 2.

Phylogenetic tree ofSalmonellaOuakamisolates. (A) Phylogenetic structure of 88SalmonellaOuakam isolates, including 79 isolates from EnteroBase and 9 isolates from Yulin. (B) Phylogenetic tree of isolates belonging to group 2 in Figure 2A.

Note: The isolates in red font represent those collected in this study.

Analysis of 116 isolates allowed for phylogenetic tree reconstruction forS.Ealing, revealing two distinct groups. Within group B, two clades, designated α and β, were identified. Isolates recovered in this study, characterized by HC5 type 435357, clustered within clade α. Two additional isolates within clade α, originating from plants in Paraguay and the Netherlands, corresponded to HC5 types 268261 and 214748, respectively. The remaining four isolates within clade α were all characterized by HC5 type 353819 (Figure 3).

Figure 3.

The phylogenetic tree of 116S. Ealing isolates, including 114 isolates from EnteroBase and 2 isolates from Yulin City, Guangxi Zhuang Autonomous Region, China.

Note: The isolates collected from Yulin are highlighted in red.