
==== Front
mSystems
mSystems
msystems
mSystems
2379-5077
American Society for Microbiology 1752 N St., N.W., Washington, DC

39158311
msystems00586-24
10.1128/msystems.00586-24
msystems.00586-24
Research Article
clinical-microbiologyClinical MicrobiologyGenomic analysis of antibiotic resistance genes and mobile genetic elements in eight strains of nontyphoid Salmonella
Liu Haibing 1 Conceptualization Funding acquisition Investigation Methodology Software Visualization Writing – original draft Writing – review and editing
Zheng Lijie 1 Investigation Methodology
Fan Huimin 1 Formal analysis Investigation
https://orcid.org/0000-0001-5036-5647
Pang Ji 1 Conceptualization Project administration Resources Supervision pangji910313@163.com

1 Department of Clinical Laboratory, The Affiliated People’s Hospital of Jiangsu University , Zhenjiang, Jiangsu, China
Editor Gilbert Jack A. University of California San Diego , La Jolla, California, USA

Address correspondence to Ji Pang, pangji910313@163.com
The authors declare no conflict of interest.

9 2024
19 8 2024
19 8 2024
9 9 e00586-2422 4 2024
12 7 2024
Copyright © 2024 Liu et al.
2024
Liu et al.
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

ABSTRACT

Nontyphoidal Salmonella (NTS) is the main etiological agent of human nontyphoidal salmonellosis. The aim of this study was to analyze the epidemiological characteristics and horizontal transfer mechanisms of antimicrobial resistance (AMR) genes from eight strains of NTS detected in Zhenjiang City, Jiangsu Province, China. Fecal samples from outpatients with food-borne diarrhea were collected in 2022. The NTS isolates were identified, and their susceptibility was tested with the Vitek 2 Compact system. The genomes of the NTS isolates were sequenced with the Illumina NovaSeq platform and Oxford Nanopore Technologies platform. The AMR genes and mobile genetic elements (MGEs) were predicted with the relevant open access resources. Eight strains of NTS were isolated from 153 specimens, and Salmonella Typhimurium ST19 was the most prevalent serotype. The AMR gene with the highest detection rate was AAC(6')-Iaa (10.5%) followed by TEM-1 (7.9%), sul2 (6.6%), and tet(A) (5.3%). Eleven MGEs carrying 34 AMR genes were identified on the chromosomes of 3 of the 8 NTS, including 3 resistance islands, 6 composite transposons (Tns), and 2 integrons. Eighteen plasmids carrying 40 AMR genes were detected in the 8 NTS strains, including 6 mobilizable plasmids, 3 conjugative plasmids, and 9 nontransferable plasmids, 7 of which carried 10 composite Tns and 3 integrons. This study provided a theoretical basis, from a genetic perspective, for the prevention and control of NTS resistance in Zhenjiang City.

IMPORTANCE

Human nontyphoidal salmonellosis is one of the common causes of bacterial food-borne illnesses, with significant social and economic impacts, especially those caused by invasive multidrug-resistant nontyphoidal Salmonella, which entails high morbidity and mortality. Antimicrobial resistance is mainly mediated by drug resistance genes, and mobile genetic elements play key roles in the capture, accumulation, and dissemination of antimicrobial resistance genes. Therefore, it is necessary to study the epidemiological characteristics and horizontal transfer mechanisms of antimicrobial resistance genes of nontyphoidal Salmonella to prevent the spread of multidrug-resistant nontyphoidal Salmonella.

KEYWORDS

nontyphoidal Salmonella
whole-genome sequencing
antimicrobial resistance gene
mobile genetic element
The Second Phase of Jinshan Medical Personnel Training Program 2022-JSYZ-13 Liu Haibing Rehabilitation Medicine Talent Leading Program KFQM-27 Liu Haibing cover-dateSeptember 2024
==== Body
pmcINTRODUCTION

Nontyphoidal Salmonella (NTS) is a gram-negative, facultatively anaerobic, nonspore-forming bacillus, belonging to the family Enterobacteriaceae, and is one of the commonest etiological agents of acute gastroenteritis disease globally. Most cases of gastroenteritis caused by NTS in humans are self-limiting diarrheal diseases and do not require treatment with antibiotics (1). However, invasive bacterial infections caused by NTS pose a particular threat to immunocompromised individuals, including those infected with human immunodeficiency virus 1 or malaria, infants, the elderly, and those with underlying comorbidities (2). If left untreated, infections with invasive NTS (iNTS) often result in meningitis, septicemia, and even death. NTS is reportedly responsible for 3.4 million enteric infections and over 680,000 associated deaths annually throughout the world (3). The areas with the highest incidence of iNTS disease are sub-Saharan Africa and Europe, with 227 and 102 cases of iNTS per 100,000 population, respectively (4). In China, the incidence of infectious NTS diseases was estimated to be 626.5 cases per 100,000 persons (5).

The emergence of bacteria with antimicrobial resistance (AMR) is a natural phenomenon that can be accelerated by the overuse of antibiotics and the action of selection pressures as well as the injudicious antimicrobial use by humans. The AMR of NTS is mainly caused by the misuse of antibiotics in animal agriculture, which is poorly regulated and defined (6). Many classes of antibiotics used to treat NTS infections in food-producing animals are the same as those used in humans, and drug-resistant NTS can be transmitted to humans via the farm-to-fork route. Between 2015 and 2018, the U.S. Centers for Disease Control and Prevention estimated that 16% of NTS were resistant to one or more antibiotics (7). The resistance rates of NTS to ampicillin (AMP), trimethoprim-sulfamethoxazole (SXT), ciprofloxacin (CIP), and ceftriaxone in China are 57.3%, 24.8%, 11.2%, and 9.9%, respectively (8). AMR is mainly mediated by resistance genes associated with enzymatic modification, energy-dependent efflux, target protection, and changes in the bacterial cell wall. Therefore, it is necessary to study the composition and distribution of the drug resistance genes of NTS to slow the development of AMR.

Mobile genetic elements (MGEs), including plasmids, integrons, phages, transposons (Tns), insertion sequences (ISs), integrative and conjugative elements (ICEs), and genomic islands, play key roles in the capture, accumulation, and dissemination of AMR genes among bacterial populations and in the evolution of bacterial genomes via conjugation, transduction, or transformation (9). The MGE-mediated horizontal transfer of foreign genes also contributes significantly to the evolution of NTS, including the acquisition of AMR and virulence genes. Wiesner et al. (10) described the IncA/C plasmids of a multidrug-resistant S. Typhimurium ST213 strain, which carried a plasmid-borne blaCMY-2 gene encoding an extended-spectrum β-lactamase, which could be transferred within bacterial communities by the conjugal machinery. The AcrAB-TolC efflux pump in Salmonella can be upregulated by the insertion of IS1 or IS10, leading to an increase in drug resistance (11).

Although NTSs that cause human infections are becoming increasingly resistant to antibiotics, there is currently no in-depth research into the composition and transmission of AMR genes in NTS in local areas of China. In recent years, whole-genome sequencing (WGS) techniques have been exploited to provide information about the complete genome of microorganisms, which can be analyzed by using bioinformatics tools to identify and provide deeper insights into molecular epidemiology typing, bacterial resistomes and virulence factors, and pan-genome plasticity mediated by MGEs (12). In this study, we sequenced the whole genomes of eight strains of NTS from outpatients in Zhenjiang City, Jiangsu Province, China, which caused acute gastroenteritis in humans, to analyze the compositions and distributions of their AMR genes and the MGEs associated with the horizontal transfer of these genes. The main objectives of this work were to provide a theoretical basis, from a genetic perspective, for the prevention and control of NTS resistance in Zhenjiang City.

RESULTS

Characterization of NTS isolates

In total, 153 specimens from patients with suspected acute food-borne gastroenteritis were collected, and 8 strains of NTS were isolated from 8 patients, including 6 children (<6 years old) and 2 elderly people (>65 years old), who were particularly at risk of acquiring iNTS disease. Serotype analyses indicated that Typhimurium was the most prevalent serotype among the eight strains (37.5%), and the serotypes of the remaining strains were Goldcoast, London, Enteritidis, Potsdam, and 1,4,[5],12:i:− (a monophasic variant of S. Typhimurium). The other related characteristics of the NTS isolates are shown in Table 1.

TABLE 1 The related characteristics of NTS isolates

Strain name	Sequence type	Serotype	Serovar antigen	
Serum agglutination	SISTRa	SeqSero	
XSK	ST19	Typhimurium	O4,12:Hi:H2	O1,4,[5],12:Hi:H1,2	O4:Hi:H1,2	
CHC	ST19	Typhimurium	O4,12:Hi:H2	O1,4,[5],12:Hi:H1,2	O4:Hi:H1,2	
ZCX	ST19	Typhimurium	O4,12:Hi:H2	O1,4,[5],12:Hi:H1,2	O4:Hi:H1,2	
ZLQ	ST34	1,4,[5],12:i:−	O4,12:Hi	O1,4,[5],12:Hi	O4:Hi	
YZY	ST2529	Goldcoast	O8:Hr:Hl,w	O6,8:Hr:Hl,w	O8:Hr:Hl,w	
FFL	ST155	London	O3,10:Hl,v:H6	O3,10,15:Hl,v:H1,6	O3,10:Hl,v:H1,6	
CYX	ST11	Enteritidis	O9,12:Hg,m	O1,9,12:Hg,m	O9:Hg,m	
ZYX	ST2039	Potsdam	O7:Hl,v:He,n,z15	O6,7,14:Hl,v:He,n,z15	O7:Hl,v:He,n,z15	
a SISTR, Salmonella In Silico Typing Resource.

Epidemiological surveillance of NTS isolates

To investigate the epidemiology of these 8 NTS strains, a minimum spanning tree of 13 strains of Salmonella was constructed based on core-genome multilocus sequence typing (cgMLST). As shown in Fig. 1, three strains of S. Typhimurium ST19 (XSK, CHC, and ZCX) shared the closest evolutionary relationship with S. Typhimurium ST34 YZU0175 (isolated from a diarrhea patient in Yangzhou City) and S. 1,4,[5],12:i:− ST34 ZLQ. No close evolutionary relationships were found for the remaining four NTS.

Fig 1 Minimum spanning tree of NTS isolates was constructed based on cgMLST.

A network diagram of different nodes labeled with various identifiers and colored based on a legend indicating ST groups. Edges between nodes have red numerical values, representing some form of similarity or distance metric.

Composition and distribution of AMR genes

A total of 76 acquired resistance genes were predicted, 40 of which (52.6%) were located on plasmids and the rest on chromosomes. The most prevalent AMR gene was AAC(6')-Iaa (10.5%), encoding resistance to aminoglycosides, followed by TEM-1 (7.9%), encoding resistance to β-lactams; sul2 (6.6%), encoding resistance to sulfonamides; and tet(A) (5.3%), encoding resistance to tetracyclines. The antimicrobial category targeted by the most resistance genes was aminoglycosides (30.3%), followed by sulfonamides (13.2%), and β-lactams and tetracyclines (9.2%) (Fig. 2a). AMR genes were detected in eight strains of NTS, and S. Typhimurium CHC had the most AMR genes (17), followed by S. London FFL (15), S. Typhimurium ZCX (13), and S. Typhimurium XSK (13). Forty-three (56.6%) AMR genes were derived from 3 strains of S. Typhimurium, and 13 AMR genes were shared by S. Typhimurium XSK and CHC. The AMR gene profiles of S. 1,4,[5],12:i:− ZLQ and S. Enteritidis CYX were also similar (Fig. 2b).

Fig 2 Composition and distribution of AMR genes. (a) Distribution of AMR genes located on chromosomes and plasmids of eight strains of NTS. (b) Heatmap of AMR gene distributions in eight strains of NTS.

A bar chart depicts the distribution of AMR genes across chromosomes and plasmids for various antibiotics, categorized by gene type and antibiotic class. A heatmap displays the presence and absence of specific AMR genes across different samples.

The correlation between the resistance phenotypes (Table S1) and the AMR genes targeting the three classes of antibiotics recommended by Clinical and Laboratory Standards Institute (CLSI) was used to evaluate the ability of WGS to predict phenotypic resistance. Overall, phenotypic resistance correlated strongly with the presence of known AMR determinants. The overall sensitivity of the AMR genes detected in predicting resistance across the three classes of antibiotics was 100% (78.5%–100%), the specificity was 80% (49.0%–94.3%), and Kappa value was 0.8 (0.4–1.0). Among these genes, ampicillin-resistance-related genes had the strongest ability to predict the resistance phenotype (Table 2).

TABLE 2 Evaluation of the ability of WGS to predict resistance phenotypesa

Antibiotic	No. of test results	Sensitivity (%), (95% CI)	Specificity (%), (95% CI)	Kappa value, (95% CI)	
PT: R	PT: S	
GT: R	GT: S	GT: R	GT: S	
AMP	6	0	0	2	100 (61.0–100)	100 (34.2–100)	1.0 (0.3–1.0)	
CIP	4	0	1	3	100 (51.0–100)	75 (30.1–95.4)	0.75 (0.1–1.0)	
SXT	4	0	1	3	100 (51.0–100)	75 (30.1–95.4)	0.75 (0.1–1.0)	
Overall	14	0	2	8	100 (78.5–100)	80 (49.0–94.3)	0.80 (0.4–1.0)	
a R, resistant; S, susceptible; PT, phenotype; GT, genotype.

Characterization of plasmids

A total of 18 plasmids were detected in 8 strains of NTS. Among these strains, S. Potsdam ZYX had the greatest number of plasmids (six), followed by S. Typhimurium CHC (four) and S. Enteritidis CYX, S. London FFL, and S. Typhimurium ZCX (two each). Of the 18 plasmids detected, 6 were identified as mobilizable plasmids, 3 as conjugative plasmids, and the remaining as nontransferable plasmids. A replicon typing analysis revealed that the mobilizable plasmids consisted of three Col(pHAD28) plasmids, two Col440I plasmids, and one ColRNAI plasmid. The conjugative plasmids included one IncHI2 plasmid, IncHI2A; one IncFIB(S) plasmid, IncFII(S); and one IncI1-I(Alpha) plasmid. A total of 40 AMR genes were identified across 8 plasmids, with 12 genes present in 2 conjugative plasmids and the remaining genes found in 6 nontransferable plasmids. The conjugative plasmids carrying AMR genes originated from S. Typhimurium ZCX, with pZCX-1 containing eight AMR genes and pZCX-2 containing four AMR genes. The pertinent details of these 18 plasmids are given in Table 3.

TABLE 3 The characterization of NTS plasmids

Plasmid name	Replicon type	Transferability	Antibiotic resistance genes	GenBank accession numbers	
pCHC-1	p0111	No	tet(A)/ AAC (3)-IIe/ OqxA/
OqxB/APH(3')-Ia	OP927715	
pCHC-2	ColRNAI	Mobilizable		OP927716	
pCHC-3	Col440I	Mobilizable		OP927717	
pCHC-4	Col(pHAD28)	Mobilizable		OP927718	
pCYX-1	IncFIB(S)
IncFII(S)	No	TEM-1	OP927719	
pCYX-2	IncX1	No	TEM-1/sul2/APH(3'')-Ib /APH (6)-Id	OP927720	
pFFL-1	IncI1-I(Alpha)	Conjugative		OP927721	
pFFL-2	IncFIB(K)	No	sul1/qacE/aadA16/dfrA27/arr-3/AAC(6')-Ib-cr6/tet(A)/APH (6)-Id/TEM-1/AAC (3)-Iid/catA2/sul2/mphA/floR	OP927722	
pXSK	Col440I	Mobilizable		OP927723	
pZCX-1	IncHI2
IncHI2A	Conjugative	CTX-M-64/arr-3/dfrA27/aadA16/sul1/tet(A)/qacE/ mphA	OP927724	
pZCX-2	IncFIB(S)
IncFII(S)	Conjugative	QnrS2/QnrS2/OqxA/OqxB	OP927725	
pZLQ	Col(pHAD28)	Mobilizable		OP927726	
pZYX-1	IncU, pKPC-CAV1321	No	sul1/dfrB4/qacE	OP927727	
pZYX-2	IncFIB(pECLA)
IncFII(pECLA)	No		OP927728	
pZYX-3	IncQ2	No	QnrS2	OP927729	
pZYX-4	Col(pHAD28)	Mobilizable		OP927730	
pZYX-5	Col440I	No		OP927731	
pZYX-6	Col(pHAD28)	No		OP927732	

Analysis of MGEs carrying AMR genes located on chromosomes

A total of 11 MGEs carrying AMR genes were identified on the chromosomes of 3 of the 8 NTS: S. Typhimurium CHC, S. Typhimurium XSK, and S. 1,4,[5],12:i:− ZLQ. These MGEs included three resistance islands, six composite Tns, and two integrons. Specifically, the chromosome of S. Typhimurium CHC carried 1 resistance island, 2 composite Tns, and 1 integron, with a total of 10 AMR genes located within these MGEs. As shown in Fig. 3a, the primary structure of the resistance island was a composite Tn containing 15 ISs, 9 AMR genes, and 1 integron carrying 4 AMR genes. The chromosome of S. Typhimurium CHC contained another composite Tn flanked by IS15DI, which comprised six ISs and one QnrS1 gene. The MGEs associated with AMR genes carried on the chromosome of S. Typhimurium XSK showed similarities to those of S. Typhimurium CHC, except for an additional tet(A) AMR gene, one IS (TnpA_TnAs1), and three other genes presented in the composite Tn located in the resistance island (Fig. 3b). The chromosome of S. 1,4,[5],12:i:− ZLQ carried five AMR genes on two composite Tns: one containing four AMR genes and five ISs, and the other carrying a tet(B) gene and four ISs. These AMR genes and related MGEs constituted the primary structure of the resistance island (Fig. 3c).

Fig 3 MGEs carrying AMR genes located on (a) chromosome of S. Typhimurium CHC; (b) chromosome of S. Typhimurium XSK; and (c) chromosome of S. 1,4,[5],12:i:− ZLQ.

Three chromosome maps illustrate resistance islands and composite transposons for CHC, XSK, and ZLQ chromosomes. Each map highlights integrons, various resistance genes, and transposon elements.

Analysis of MGEs carrying AMR genes located on plasmids

A total of 13 MGEs carrying AMR genes were identified on 7 plasmids and included 10 Tns and 3 integrons. The plasmid with the greatest number of MGEs was pFFL-2 (four), followed by pZCX-1 (three) and pZYX-1 (two). Nonmobile plasmids, pCYX-1 and pCYX-2, were isolated from S. Enteritidis CYX, and pCYX-1 carried a noncomposite Tn (Tn3) containing a TEM-1 gene (Fig. 4a), and pCYX-2 carried a composite Tn flanked by Tn3 and IS26, which harbored four AMR genes (Fig. 4b). The 14 AMR genes carried by pFFL-2 were distributed in 3 composite Tns, and a single integron containing 6 AMR genes was located in one of the composite Tns. A total of 10 ISs were involved in the formation of the 3 composite Tns, with IS26 accounting for 60% (Fig. 4c). The two conjugative plasmids isolated from S. Typhimurium ZCX, pZCX-1 and pZCX-2, contained three and one MGEs, respectively. The MGEs of pZCX-1 were two composite Tns and one integron, and eight AMR genes distributed across the two composite Tns (Fig. 4d). The MGE located in pZCX-2 was a composite Tn carrying four AMR genes, including two QnrS2 genes flanked by IS26, and the OqxA and OqxB gene flanked by IS15DI (Fig. 4e). The composite Tn located on pCHC-1 carried five AMR genes, which were distinct from those on the chromosome of S. Typhimurium CHC (Fig. 4f). The composite Tn present on pZYX-1 carried one integron, which harbored three AMR genes: dfrB4, qacE, and sul1 (Fig. 4g).

Fig 4 MGEs carrying AMR genes located on (a) pCYX-1, (b) pCYX-2, (c) pFFL-2, (d) pZCX-1, (e) pZCX-2, (f) pCHC-1, and (g) pZYX-1.

Seven plasmid maps highlight resistance islands and composite transposons for pCYX-1, pCYX-2, pZCX-1, pZCX-2, pFFL-2, pCHC-1, and pZYX-1. Each map illustrates integrons, various resistance genes, and transposon elements.

Composition and distribution of ISs

A total of 277 ISs were identified in the 8 NTS strains, 211 of which were located on chromosomes and 66 on plasmids. The ISs with the highest detection rate were ISEch5 (15.9%), followed by IS26 (11.2%), ISEhe3 (8.7%), and ISSen1 (8.7%). ISEch5 and ISSen1 were only detected on chromosomes, and the plasmid-located IS with the highest detection rate was IS26 (Fig. 5a). Of the 277 ISs, 25% were involved in the formation of composite Tns, including 38 ISs located in composite Tns on chromosomes and 32 in composite Tns on plasmids. Among these, IS26 was the most numerous IS (28) that constituted composite Tns, followed by IS15DI (13) and ISKpn26 (4) (Fig. 5b).

Fig 5 Composition and distribution of ISs identified in the eight NTS strains. (a) ISs located on chromosomes and plasmids. (b) Distribution of ISs involved in the formation of composite Tns.

Two bar charts depict the number of insertion sequences in chromosomes and plasmids for various IS types and the breakdown of the most frequent IS types, exhibiting their distribution across chromosomes and plasmids.

Identification of ICEs and integrative or mobilizable elements

An ICE with a length of 80,793 bp was detected on the chromosome of S. 1,4,[5],12:i:− ZLQ, which was composed of 87 genes, including 24 genes related to the type IV secretion system, 2 genes encoding relaxase, 1 gene encoding integrase, 1 gene for the coupling protein to a type IV conjugative protein complex, and 59 genes encoding other proteins (Fig. 6a). A total of five integrative or mobilizable elements (IMEs) were identified in five NTS strains: S. Typhimurium CHC, S. Typhimurium XSK, S. Typhimurium ZCX, S. Goldcoast YZY, and S. London FFL. Each of these IMEs contained two important genes encoding integrase and relaxase, and the relaxase predominantly belonged to the MobH family (Fig. 6b).

Fig 6 Identification of ICEs and IMEs. (a) ICE located on the chromosome of S. 1,4,[5],12:i:− ZLQ. (b) IMEs identified in five strains of NTS.

Gene maps of mobile genetic elements depict ZLQ-ICE, indicating various gene components along the sequence (a). IME elements for XSK, YZY, ZCX, CHC, FFL, and other components with positions marked (b).

Analysis of prophages in NTS

A total of 21 intact prophage regions were detected in the 8 NTS strains, and the isolates with the greatest numbers of prophages detected were S. 1,4,[5],12:i:− ZLQ and S. Typhimurium ZCX (four each), followed by S. Potsdam ZYX, S. London FFL, and S. Enteritidis CYX (three each), S. Typhimurium CHC (two), and S. Goldcoast YZY and S. Typhimurium XSK (one each). Among these intact prophages, the most prevalent profiles were Phage_Gifsy_2 (six), Salmon_118970_sal3 (four), and Phage_Gifsy_1 (three), which were present in six, three, and two genomes, respectively (Fig. 7).

Fig 7 Prophage regions identified in eight strains of NTS.

A circular plot of prophages in various samples. Different gene components are color-coded, including lysis, protease, terminase, coat protein, attachment site, integrase, and more.

DISCUSSION

Food-borne diseases are a serious public health problem in China, entailing considerable morbidity and mortality. In recent years, NTS has been recognized as the dominant pathogen of food-borne diseases in China, and S. Typhimurium and S. Enteritidis are the predominant serotypes, causing approximately 45.5% of NTS infections (13). Similar conclusions were drawn in our study, in that S. Typhimurium was the most prevalent serotype, accounting for 37.5% of the isolates identified. S. 1,4,[5],12:i:− is a monophasic variant of S. Typhimurium lacking the fljB-encoded second-phase H antigen, so S. 1,4,[5],12:i:− is closely related to S. Typhimurium. In this study, the close homology between S. 1,4,[5],12:i:− ZLQ and three strains of S. Typhimurium was also reflected in the epidemiology. We also found that a strain of NTS (S. Typhimurium YZU0175) with the same ST as S. 1,4,[5],12:i:− ZLQ caused diarrhea in a human from Yangzhou City, which is very close to Zhenjiang City.

Theoretically, any drug resistance phenotype of an organism can be derived from the AMR genotype. Therefore, it is necessary to study the composition and distribution of AMR genes to monitor and prevent multidrug-resistant NTS. In this study, a total of 76 AMR genes were predicted, including 71 determinants encoding resistance to 9 different antimicrobial categories and 5 determinants of disinfectant resistance. Most of the AMR genes targeted the antimicrobial categories of aminoglycosides, sulfonamides, β-lactams, and tetracyclines, and the most prevalent determinants in these categories were AAC(6')-Iaa, sul2, TEM-1, and tet(A), respectively. Chen et al. reported that similar AMR gene patterns existed in NTS from Shaoxing City, in which 100% of isolates contained the AAC(6′)-Iaa_1 gene, 65.52% were positive for blaTEM−1B, and 52.87% contained the tet(A) gene (14). Predicting AMR genes with WGS technologies is an emerging technique for identifying drug-resistant bacteria to guide clinical treatment, although there are still many issues. For example, it is difficult to predict AMR accurately with WGS when it is caused by new genetic mutations, new resistance mechanisms, or the increased expression of intrinsic resistance genes (15). In this study, we found that identifying NTS strains clinically resistant to fluoroquinolones, AMP, or SXT with WGS was as reliable as conventional drug sensitivity testing, especially in predicting AMP resistance, for which 100% (61.0%–100%) sensitivity , 100% (34.2%–100%) specificity, and 1.0 (0.3–1.0) Kappa value were achieved. Our overall finding of strong correlations between phenotypic and genotypic resistance in Salmonella is consistent with the study of 640 strains of Salmonella from humans and retail meats, in which the resistance genotypes for AMP, CIP, and SXT correlated with the phenotypes with 99.7%, 100%, and 100% specificity, respectively, and 99.6%, 100%, and 86.4% sensitivity, respectively (16).

As important vectors in the transmission of AMR genes, MGEs play a very significant role in the evolution of NTS genomes, allowing them to adapt rapidly to selective pressures, especially in response to antimicrobial exposure. Transmissible plasmids, which have featured prominently as agents of the horizontal transfer of AMR genes in Salmonella, can be divided into two classes based on self-transmission: (i) conjugative plasmids, which contain the essential components for conjugation; and (ii) mobilizable plasmids, which must coexist in the donor bacterium with a conjugative plasmid to be mobilized by conjugation (17). In this study, the detection rate of transmissible plasmids was high, and they accounted for 50% of the total number of plasmids. The conjugative plasmids pZCX-1 and pZCX-2 are multiresistance plasmids carrying eight and three AMR genes, respectively, which make S. Typhimurium ZCX a high-risk strain for the transmission of AMR genes. The conjugative plasmid pFFL-1, which carries no AMR genes, coexisted with plasmid pFFL-2 in S. London FFL, and the 14 AMR genes carried by pFFL-2 could be transferred to pFFL-1 by composite Tns or integron to transmit by conjugation. No NTS with coexisting mobilizable plasmids and conjugative plasmids was detected in our study. Bacterial transposable elements are divided into four categories: ISs, composite Tns, Tn3 family, and prophages. In general, IS elements are the simplest autonomous mobile elements and can affect bacterial antibiotic resistance by gene inactivation or the transmission of AMR genes as part of a composite Tn. It has been reported that specific ISs composed of composite Tns correlate strongly with AMR genes. For example, ISKpn40 was present near all mcr-3 genes, and a strong association between IS30 and mcr-1.1 was observed (18). Similarly, there was an association between ISs and AMR genes in the NTS we identified. For example, the IS15DI-TnpR_Tn3-TEM-1-IS26 unit and IS15DI-QnrS1-ISAs17-ISEae2-ISEc36-TnpA_Tn3-IS15DI composite Tn were present in the chromosomes of S. Typhimurium XSK and S. Typhimurium CHC, and the IS15DI-TnpR_Tn3-TEM-1-IS26 unit was also detected on the chromosome of S. 1,4,[5],12:i:− ZLQ. Modules composed of specific ISs and AMR genes were identified in multiple plasmids, including the APH (6)-Id-APH(3'')-Ib-sul2-IS26 unit and IS15DI-OqxB-OqxA-IS15DI unit. Transposons, which commonly act as the primary reservoirs of AMR genes, provide functionality for the transfer of bacterial DNA segments to another part of the genome. Composite Tns flanked by ISs are the main form of bacterial transposons, and the distribution of ISs and AMR genes located on chromosomal composite Tns differs from that of plasmid-borne composite Tns. This unbalanced distribution of ISs and AMR genes located in composite Tns was also confirmed in our study. For example, AMR genes, including AAC(6')-Iaa, aadA1, aadA2, sul3, tet(M), and QnrS1, and ISs, including IS5, ISVsa3, ISEc36, and ISEae2, were detected only on chromosomes, whereas AMR genes, such as QnrS2, aadA16, and dfrA27, and ISs, such as IS903B, ISCfr1, and IS5075, were only present on plasmids. Prophage sequences are a major source of genomic variability, providing a large number of AMR genes, virulence genes, and toxin genes for recipient bacteria. Salmonella enterica genomes contain large numbers of prophage profiles characterized by a high degree of variability (19). In the present study, 11 types of prophages were identified, and Salmon_118,970_sal3 and Phage_Gifsy_2 were highly prevalent. These results are consistent with those of previous studies, which reported that Salmon_118,970_sal3 was the most prevalent prophage profile in Salmonella (20). Furthermore, no AMR genes were detected in these prophage profiles. It is well known that integrons play a key role in the transmission of AMR genes and the evolution of S. enterica, especially class I integrons, which are widely distributed in S. enterica isolates (21). The five integrons identified in the present study were also class 1 integrons, which are transferred across bacterial genomes as components of Tns. Our results also showed that the AMR genes carried by two chromosomal integrons were quite different from those carried by three plasmid integrons, and the aadA1, aadA2, cmlA1, and dfrA12 genes were mainly located in chromosomal integrons, whereas sul1, qacE, aadA16, dfrA27, and arr-3 genes were mainly present in plasmid integrons. Resistance islands, defined as genomic islands that contain multiple resistance determinants, have a significant impact on the dissemination of AMR by horizontal gene transfer (22). Salmonella genomic island 1 (SGI1) is the first resistance island found in S. enterica Typhimurium DT104 isolates and contributes significantly to antibiotic resistance by carrying a class 1 integron that contains five AMR genes (23). In the present study, three resistance islands were identified, and almost all chromosomal AMR genes were located in them. They also carried multiple MGEs, including class 1 integrons, composite Tns, and ISs. The three resistance islands carried a greater number of AMR genes than SGI1, so they are potentially important vectors for the horizontal transfer of chromosomal AMR genes in Salmonella.

Conclusions

To date, few studies have investigated the composition characteristics and horizontal transfer mechanisms of AMR genes from NTS in Zhenjiang City. In the present study, the AMR genes and MGEs of eight NTS strains from Zhenjiang City were identified with WGS combined with a bioinformatic analysis. Salmonella Typhimurium ST19 was the leading cause of nontyphoidal salmonellosis in Zhenjiang City. The most prevalent AMR genes included AAC(6')-Iaa, TEM-1, sul2, and tet(A), which undergo horizontal gene transfer mainly on plasmids, composite Tns, and integrons. This study was limited by the low detection rate of NTS. Despite this, our study provides comprehensive information to better manage and prevent the spread of NTS AMR in Zhenjiang City.

MATERIALS AND METHODS

Sample collection and NTS identification

During 2022, fecal samples were collected from outpatients with food-borne diarrhea in Zhenjiang City, whose clinical presentations included diarrhea, vomiting, fever, nausea, and abdominal cramps. Patients taking any type of antibiotic were excluded from the study. Bacterial cultures were established by inoculating Salmonella-Shigella agar (CHROMagar) with the fecal specimens within 1 h of receiving them and incubating them at 37°C for 18–24 h. Presumptive colonies were further screened by testing in triple sugar iron agar, and the isolates were identified with the Vitek 2 Compact system (bioMérieux). The flagellar and lipopolysaccharide antigens were determined by slide agglutination with commercial antiserum (Ningbo Tianrun Bio-Pharmaceutical Co. Ltd), and the serotypes were assigned according to the Kauffmann-White scheme (24).

Antimicrobial susceptibility testing

The antibiotics recommended by the Clinical and Laboratory Standards Institute for the treatment of enteric Salmonella infections include fluoroquinolones, AMP, and SXT. Therefore, the antimicrobial susceptibility of the NTS isolates for these three classes of antibiotics was tested with the Vitek 2 Compact system (bioMérieux), and the tested antimicrobial agents were AMP, CIP, and SXT. The results were interpreted based on the CLSI criteria (2022). Intermediate phenotypes were counted as resistant in this analysis. Escherichia coli ATCC 25922 was used as the control strain.

Genomic analysis of AMR genes and MGEs

The NTS isolates were cultured in nutrient broth (Hangzhou Microbial Reagent Co. Ltd) for 24 h at 36°C, collected, and immediately transported to a commercial laboratory (Personalbio) for sequencing with the Illumina NovaSeq platform and the Oxford Nanopore Technologies (ONT) platform. Genomic analysis of AMR genes and MGEs mediated by WGS technology is divided into two parts as dry lab and wet lab. The term dry lab denotes the bioinformatic analysis, whereas the wet lab is related to molecular biology practical experiments (25).

Wet lab section

The genomic DNA was extracted with the CTAB method. The library for next-generation sequencing was constructed according to the standard preparative procedure of the Illumina TruSeq DNA LT Sample Preparation Kit, and the quality of the library was assessed with an Agilent Bioanalyzer 2100 according to the instructions of the Agilent High Sensitivity DNA Kit (Agilent Technologies Inc.). The QuantiFluor dsDNA System (Promega) was used for the quantitative analysis of the library according to the standard protocol of the Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen). Paired-end sequencing was performed when the library met the quality requirements. The concentration of the genomic DNA used to construct a three-generation sequencing library was assessed with the Qubit dsDNA HS Assay Kit (Invitrogen), after which its purity was evaluated with a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific). Finally, the integrity of the DNA was confirmed with 0.35% agarose gel electrophoresis. The construction of the three-generation sequencing library, the assessment of its quality, and its sequencing were performed according to the standard protocols provided by ONT.

Dry lab section

The quality control of the raw reads for next-generation sequencing was performed with FastQC (26), adapter contamination was removed with Adapter Removal (version 2) (27), and all the reads were quality corrected with SOAPec (version 2.0) (28). The three-generation sequencing data were assembled with Unicycler (29) and Flye (30), and the assembled contigs were corrected with high-quality second-generation sequencing data using Pilon v1.18 (31). The serotypes of the NTS strains were predicted with the web-based tool SeqSero2 (v1.2.1) (32) and the Salmonella In Silico Typing Resource (v1.1.1) (33). The sequence type (ST) was identified with multilocus sequence typing (v2.22.0) (34). The publicly available resistance gene database ResFinder 4.4.2 was used to predict acquired resistance genes, with a threshold of 90% (35). The replicon type and the transferability of plasmids were determined with PlasmidFinder (v2.1.6) (36) and oriTfinder 1.1 (37), respectively. VRprofile 2 (38) was used to identify the MGEs carrying AMR genes, including integrons, composite Tns, and resistance islands. The MGEs that did not carry AMR genes were predicted with other open access resources, including IS elements with Isfinder (39), ICEs and IMEs with ICEberg (40) (with the BLAST default parameters), and prophages with PHASTER (41) (only “intact” prophages were used in this analysis). The minimum spanning tree was constructed with PHYLOViZ online (42) based on core-genome multilocus sequence typing, and five strains of Salmonella isolated from humans in Jiangsu Province were also used in the epidemiological analysis.

Evaluation of the ability of WGS to predict phenotypic resistance

The ability of WGS to predict phenotypic resistance was evaluated by evaluating the correlation between the resistance phenotypes and AMR genes. When the phenotypic results were used as the reference standard, sensitivity was equal to the number of genotypically resistant strains divided by the total number of isolates showing phenotypic resistance, and specificity was equal to the number of isolates that were genotypically susceptible divided by the total number of strains with susceptible phenotypes. Correlation between phenotypic and genotypic antibiotic resistance was determined using Cohen’s Kappa coefficient of agreement and 95% confidence interval (CI) values. The kappa value was interpreted as follows: values of 0.81–1.00 indicated almost perfect agreement, values of 0.61CI0.80 indicated substantial agreement, values of 0.41CI0.60 indicated moderate agreement, values of 0.21CI0.40 indicated fair agreement, and values of ≤0.20 indicated slight to no agreement (43). Sensitivity, specificity, Kappa value, and 95% CI were calculated using OpenEpi version 3.01 (44).

ACKNOWLEDGMENTS

We thank Personalbio (Shanghai, China) for technical assistance and International Science Editing for editing the manuscript.

SUPPLEMENTAL MATERIAL

The following material is available online at https://doi.org/10.1128/msystems.00586-24.

10.1128/msystems.00586-24.SuF1 Table S1 msystems.00586-24-s0001.pdf

Antibiotic resistance patterns of eight strains of nontyphoidal Salmonella.

10.1128/msystems.00586-24.SuF2 Table S2 msystems.00586-24-s0002.xlsx

Gene profiles for Salmonella chromosome.

ASM does not own the copyrights to Supplemental Material that may be linked to, or accessed through, an article. The authors have granted ASM a non-exclusive, world-wide license to publish the Supplemental Material files. Please contact the corresponding author directly for reuse.
==== Refs
REFERENCES

1 McDermott PF, Zhao S, Tate H. 2018. Antimicrobial resistance in nontyphoidal Salmonella. Microbiol Spectr 6 . doi:10.1128/microbiolspec.ARBA-0014-2017
2 Balasubramanian R, Im J, Lee JS, Jeon HJ, Mogeni OD, Kim JH, Rakotozandrindrainy R, Baker S, Marks F. 2019. The global burden and epidemiology of invasive non-typhoidal Salmonella infections. Hum Vaccin Immunother 15 :1421–1426. doi:10.1080/21645515.2018.1504717 30081708
3 Ao TT, Feasey NA, Gordon MA, Keddy KH, Angulo FJ, Crump JA. 2015. Global burden of invasive nontyphoidal Salmonella disease, 2010(1). Emerg Infect Dis 21 :941–949. doi:10.3201/eid2106.140999 25860298
4 Kirk MD, Pires SM, Black RE, Caipo M, Crump JA, Devleesschauwer B, Döpfer D, Fazil A, Fischer-Walker CL, Hald T, Hall AJ, Keddy KH, Lake RJ, Lanata CF, Torgerson PR, Havelaar AH, Angulo FJ. 2015. World Health Organization estimates of the global and regional disease burden of 22 foodborne bacterial, protozoal, and viral diseases, 2010: a data synthesis. PLoS Med 12 :e1001921. doi:10.1371/journal.pmed.1001921 26633831
5 Xu Y, Zhou X, Jiang Z, Qi Y, Ed-Dra A, Yue M. 2020. Epidemiological investigation and antimicrobial resistance profiles of Salmonella isolated from breeder chicken hatcheries in Henan, China. Front Cell Infect Microbiol 10 :497. doi:10.3389/fcimb.2020.00497 33042870
6 Brunelle BW, Bearson SMD, Bearson BL. 2013. Tetracycline accelerates the temporally-regulated invasion response in specific isolates of multidrug-resistant Salmonella enterica serovar Typhimurium. BMC Microbiol 13 :202. doi:10.1186/1471-2180-13-202 24020473
7 McMillan EA, Jackson CR, Frye JG. 2020. Transferable Plasmids of Salmonella enterica associated with antibiotic resistance genes. Front Microbiol 11 :562181. doi:10.3389/fmicb.2020.562181 33133037
8 Song W, Shan Q, Qiu Y, Lin X, Zhu C, Zhuo Z, Wang C, Tong J, Li R, Wan C, Zhu Y, Chen M, Xu Y, Lin D, Wu S, Jia C, Gao H, Yang J, Zhao S, Zeng M, Collaborative Working Group of the Pediatric Subgroup of the China Society of Infectious D. 2022. Clinical profiles and antimicrobial resistance patterns of invasive Salmonella infections in children in China. Eur J Clin Microbiol Infect Dis 41 :1215–1225. doi:10.1007/s10096-022-04476-7 36040531
9 Ghaly TM, Gillings MR. 2022. New perspectives on mobile genetic elements: a paradigm shift for managing the antibiotic resistance crisis. Philos Trans R Soc Lond B Biol Sci 377 :20200462. doi:10.1098/rstb.2020.0462 34839710
10 Wiesner M, Calva E, Fernández-Mora M, Cevallos MA, Campos F, Zaidi MB, Silva C. 2011. Salmonella Typhimurium ST213 is associated with two types of IncA/C plasmids carrying multiple resistance determinants. BMC Microbiol 11 :9. doi:10.1186/1471-2180-11-9 21223599
11 Siguier P, Gourbeyre E, Chandler M. 2014. Bacterial insertion sequences: their genomic impact and diversity. FEMS Microbiol Rev 38 :865–891. doi:10.1111/1574-6976.12067 24499397
12 Karampatakis T, Tsergouli K, Behzadi P. 2024. Pan-genome plasticity and virulence factors: a natural treasure trove for Acinetobacter baumannii. Antibiotics (Basel) 13 :257. doi:10.3390/antibiotics13030257 38534692
13 Zhan Z, Xu X, Gu Z, Meng J, Wufuer X, Wang M, Huang M, Chen J, Jing C, Xiong Z, Zeng M, Liao M, Zhang J. 2019. Molecular epidemiology and antimicrobial resistance of invasive non-typhoidal Salmonella in China, 2007-2016. Infect Drug Resist 12 :2885–2897. doi:10.2147/IDR.S210961 31571942
14 Chen J, Ed-Dra A, Zhou H, Wu B, Zhang Y, Yue M. 2022. Antimicrobial resistance and genomic investigation of non-typhoidal Salmonella isolated from outpatients in Shaoxing city, China. Front Public Health 10 :988317. doi:10.3389/fpubh.2022.988317 36176509
15 Ellington MJ, Ekelund O, Aarestrup FM, Canton R, Doumith M, Giske C, Grundman H, Hasman H, Holden MTG, Hopkins KL, Iredell J, Kahlmeter G, Köser CU, MacGowan A, Mevius D, Mulvey M, Naas T, Peto T, Rolain J-M, Samuelsen Ø, Woodford N. 2017. The role of whole genome sequencing in antimicrobial susceptibility testing of bacteria: report from the EUCAST subcommittee. Clin Microbiol Infect 23 :2–22. doi:10.1016/j.cmi.2016.11.012 27890457
16 McDermott PF, Tyson GH, Kabera C, Chen Y, Li C, Folster JP, Ayers SL, Lam C, Tate HP, Zhao S. 2016. Whole-genome sequencing for detecting antimicrobial resistance in nontyphoidal Salmonella. Antimicrob Agents Chemother 60 :5515–5520. doi:10.1128/AAC.01030-16 27381390
17 Garcillán-Barcia MP, Alvarado A, de la Cruz F. 2011. Identification of bacterial plasmids based on mobility and plasmid population biology. FEMS Microbiol Rev 35 :936–956. doi:10.1111/j.1574-6976.2011.00291.x 21711366
18 Jia C, Wang Z, Huang C, Teng L, Zhou H, An H, Liao S, Liu Y, Huang L, Tang B, Yue M. 2023. Mobilome-driven partitions of the resistome in Salmonella. mSystems 8 :e00883-23. doi:10.1128/msystems.00883-23 37855620
19 Gymoese P, Kiil K, Torpdahl M, Østerlund MT, Sørensen G, Olsen JE, Nielsen EM, Litrup E. 2019. WGS based study of the population structure of Salmonella enterica serovar Infantis. BMC Genomics 20 :870. doi:10.1186/s12864-019-6260-6 31730461
20 Ksibi B, Ktari S, Ghedira K, Othman H, Maalej S, Mnif B, Fabre L, Rhimi F, Hello SL, Hammami A. 2022. Antimicrobial resistance genes, virulence markers and prophage sequences in Salmonella enterica serovar Enteritidis isolated in Tunisia using whole genome sequencing. Curr Res Microb Sci 3 :100151. doi:10.1016/j.crmicr.2022.100151 35909609
21 Miriagou V, Carattoli A, Fanning S. 2006. Antimicrobial resistance islands: resistance gene clusters in Salmonella chromosome and plasmids. Microbes Infect 8 :1923–1930. doi:10.1016/j.micinf.2005.12.027 16716633
22 Juhas M, van der Meer JR, Gaillard M, Harding RM, Hood DW, Crook DW. 2009. Genomic islands: tools of bacterial horizontal gene transfer and evolution. FEMS Microbiol Rev 33 :376–393. doi:10.1111/j.1574-6976.2008.00136.x 19178566
23 Hall RM. 2010. Salmonella genomic islands and antibiotic resistance in Salmonella enterica. Future Microbiol 5 :1525–1538. doi:10.2217/fmb.10.122 21073312
24 Sedrakyan AM, Ktsoyan ZA, Arakelova KA, Zakharyan MK, Hovhannisyan AI, Gevorgyan ZU, Mnatsakanyan AA, Kakabadze EG, Makalatia KB, Chanishvili NA, Pirnay JP, Arakelyan AA, Aminov RI. 2020. Extended-spectrum β-lactamases in human isolates of multidrug-resistant non-typhoidal Salmonella enterica. Front Microbiol 11 :592223. doi:10.3389/fmicb.2020.592223 33414769
25 Behzadi P, Ranjbar R. 2019. DNA microarray technology and bioinformatic web services. Acta Microbiol Immunol Hung 66 :19–30. doi:10.1556/030.65.2018.028 30010394
26 Patel RK, Jain M. 2012. NGS QC Toolkit: a toolkit for quality control of next generation sequencing data. PLoS One 7 :e30619. doi:10.1371/journal.pone.0030619 22312429
27 Schubert M, Lindgreen S, Orlando L. 2016. AdapterRemoval v2: rapid adapter trimming, identification, and read merging. BMC Res Notes 9 :88. doi:10.1186/s13104-016-1900-2 26868221
28 Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, He G, Chen Y, Pan Q, Liu Y, et al. . 2012. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience 1 :18. doi:10.1186/2047-217X-1-18 23587118
29 Wick RR, Judd LM, Gorrie CL, Holt KE. 2017. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol 13 :e1005595. doi:10.1371/journal.pcbi.1005595 28594827
30 Freire B, Ladra S, Parama JR. 2022. Memory-efficient assembly using Flye. IEEE/ACM Trans Comput Biol Bioinform 19 :3564–3577. doi:10.1109/TCBB.2021.3108843 34469305
31 Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, Cuomo CA, Zeng Q, Wortman J, Young SK, Earl AM. 2014. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One 9 :e112963. doi:10.1371/journal.pone.0112963 25409509
32 Zhang S, den Bakker HC, Li S, Chen J, Dinsmore BA, Lane C, Lauer AC, Fields PI, Deng X. 2019. SeqSero2: rapid and improved Salmonella serotype determination using whole-genome sequencing data. Appl Environ Microbiol 85 :e01746-19. doi:10.1128/AEM.01746-19 31540993
33 Yoshida CE, Kruczkiewicz P, Laing CR, Lingohr EJ, Gannon VPJ, Nash JHE, Taboada EN. 2016. The Salmonella In Silico Typing Resource (SISTR): an open web-accessible tool for rapidly typing and subtyping draft Salmonella genome assemblies. PLoS One 11 :e0147101. doi:10.1371/journal.pone.0147101 26800248
34 Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, Marvig RL, Jelsbak L, Sicheritz-Pontén T, Ussery DW, Aarestrup FM, Lund O. 2012. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol 50 :1355–1361. doi:10.1128/JCM.06094-11 22238442
35 Florensa AF, Kaas RS, Clausen P, Aytan-Aktug D, Aarestrup FM. 2022. ResFinder - an open online resource for identification of antimicrobial resistance genes in next-generation sequencing data and prediction of phenotypes from genotypes. Microb Genom 8 :000748. doi:10.1099/mgen.0.000748 35072601
36 Carattoli A, Hasman H. 2020. PlasmidFinder and In Silico pMLST: identification and typing of plasmid replicons in whole-genome sequencing (WGS). Methods Mol Biol 2075 :285–294. doi:10.1007/978-1-4939-9877-7_20 31584170
37 Li X, Xie Y, Liu M, Tai C, Sun J, Deng Z, Ou HY. 2018. oriTfinder: a web-based tool for the identification of origin of transfers in DNA sequences of bacterial mobile genetic elements. Nucleic Acids Res 46 :W229–W234. doi:10.1093/nar/gky352 29733379
38 Wang M, Goh YX, Tai C, Wang H, Deng Z, Ou HY. 2022. VRprofile2: detection of antibiotic resistance-associated mobilome in bacterial pathogens. Nucleic Acids Res 50 :W768–W773. doi:10.1093/nar/gkac321 35524563
39 Siguier P, Varani A, Perochon J, Chandler M. 2012. Exploring bacterial insertion sequences with ISfinder: objectives, uses, and future developments. Methods Mol Biol 859 :91–103. doi:10.1007/978-1-61779-603-6_5 22367867
40 Bi D, Xu Z, Harrison EM, Tai C, Wei Y, He X, Jia S, Deng Z, Rajakumar K, Ou HY. 2012. ICEberg: a web-based resource for integrative and conjugative elements found in bacteria. Nucleic Acids Res 40 :D621–6. doi:10.1093/nar/gkr846 22009673
41 Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, Wishart DS. 2016. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res 44 :W16–W21. doi:10.1093/nar/gkw387 27141966
42 Ribeiro-Gonçalves B, Francisco AP, Vaz C, Ramirez M, Carriço JA. 2016. PHYLOViZ online: web-based tool for visualization, phylogenetic inference, analysis and sharing of minimum spanning trees. Nucleic Acids Res 44 :W246–51. doi:10.1093/nar/gkw359 27131357
43 Wu SH, Xiao YX, Hsiao HC, Jou R. 2022. Development and assessment of a novel whole-gene-based targeted next-generation sequencing assay for detecting the susceptibility of Mycobacterium tuberculosis to 14 drugs. Microbiol Spectr 10 :e0260522. doi:10.1128/spectrum.02605-22 36255328
44 Sullivan KM, Dean A, Soe MM. 2009. OpenEpi: a web-based epidemiologic and statistical calculator for public health. Public Health Rep 124 :471–474. doi:10.1177/003335490912400320 19445426
