Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genomic comparison of diverse Salmonella serovars isolated from swine

  • Sushim K. Gupta,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Writing – original draft, Writing – review & editing

    Affiliation Bacterial Epidemiology and Antimicrobial Resistance Unit, USDA-ARS, Athens, GA, United States of America

  • Poonam Sharma,

    Roles Formal analysis, Validation, Writing – review & editing

    Affiliation Bacterial Epidemiology and Antimicrobial Resistance Unit, USDA-ARS, Athens, GA, United States of America

  • Elizabeth A. McMillan,

    Roles Investigation, Validation, Writing – review & editing

    Affiliations Bacterial Epidemiology and Antimicrobial Resistance Unit, USDA-ARS, Athens, GA, United States of America, Department of Microbiology, University of Georgia, Athens, GA, United States of America

  • Charlene R. Jackson,

    Roles Conceptualization, Funding acquisition, Validation, Writing – review & editing

    Affiliation Bacterial Epidemiology and Antimicrobial Resistance Unit, USDA-ARS, Athens, GA, United States of America

  • Lari M. Hiott,

    Roles Methodology, Validation, Writing – review & editing

    Affiliation Bacterial Epidemiology and Antimicrobial Resistance Unit, USDA-ARS, Athens, GA, United States of America

  • Tiffanie Woodley,

    Roles Methodology, Validation, Writing – review & editing

    Affiliation Bacterial Epidemiology and Antimicrobial Resistance Unit, USDA-ARS, Athens, GA, United States of America

  • Shaheen B. Humayoun,

    Roles Writing – review & editing

    Affiliation Bacterial Epidemiology and Antimicrobial Resistance Unit, USDA-ARS, Athens, GA, United States of America

  • John B. Barrett,

    Roles Methodology

    Affiliation Bacterial Epidemiology and Antimicrobial Resistance Unit, USDA-ARS, Athens, GA, United States of America

  • Jonathan G. Frye ,

    Contributed equally to this work with: Jonathan G. Frye, Michael McClelland

    Roles Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing

    Jonathan.Frye@ars.usda.gov

    ‡ These authors are joint senior authors on this work.

    Affiliation Bacterial Epidemiology and Antimicrobial Resistance Unit, USDA-ARS, Athens, GA, United States of America

  • Michael McClelland

    Contributed equally to this work with: Jonathan G. Frye, Michael McClelland

    Roles Conceptualization, Funding acquisition, Writing – review & editing

    ‡ These authors are joint senior authors on this work.

    Affiliation Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, CA, United States of America

Abstract

Food animals act as a reservoir for many foodborne pathogens. Salmonella enterica is one of the leading pathogens that cause food borne illness in a broad host range including animals and humans. They can also be associated with a single host species or a subset of hosts, due to genetic factors associated with colonization and infection. Adult swine are often asymptomatic carriers of a broad range of Salmonella servoars and can act as an important reservoir of infections for humans. In order to understand the genetic variations among different Salmonella serovars, Whole Genome Sequences (WGS) of fourteen Salmonella serovars from swine products were analyzed. More than 75% of the genes were part of the core genome in each isolate and the higher fraction of gene assign to different functional categories in dispensable genes indicated that these genes acquired for better adaptability and diversity. High concordance (97%) was detected between phenotypically confirmed antibiotic resistances and identified antibiotic resistance genes from WGS. The resistance determinants were mainly located on mobile genetic elements (MGE) on plasmids or integrated into the chromosome. Most of known and putative virulence genes were part of the core genome, but a small fraction were detected on MGE. Predicted integrated phage were highly diverse and many harbored virulence, metal resistance, or antibiotic resistance genes. CRISPR (Clustered regularly interspaced short palindromic repeats) patterns revealed the common ancestry or infection history among Salmonella serovars. Overall genomic analysis revealed a great deal of diversity among Salmonella serovars due to acquired genes that enable them to thrive and survive during infection.

Introduction

Salmonella enterica subsp. enterica has the ability to infect a wide range of hosts, including both animals and humans. In the latter group, the bacterium causes foodborne illnesses ranging from mild diarrhea and gastroenteritis to severe systemic infections such as enteric fever [1], second only to norovirus as the causative agent of foodborne illness; Salmonella infections in humans are the leading cause of hospitalization and deaths from foodborne illness in the United States [2]. The Centers for Disease Control and Prevention has estimated 1.2 million foodborne illnesses with 19,000 hospitalizations and approximately 380 deaths in the United States annually (https://www.cdc.gov/salmonella/index.html). S. enterica ranks as the leading cause of foodborne disease as measured by the combined cost of illness and quality-adjusted life-year [3]. The direct economic losses owing to salmonellosis in the US exceed an estimated $3.5 billion per year. More than 2,600 Salmonella serovars has been confirmed by the agglutination properties of the somatic O, flagellar H, and capsular Vi antigens [4].

Inappropriate use of antimicrobials in food animal production during treatment and prevention of diseases and for growth promotion contribute to resistance, including acquisition of antibiotic resistance (AR) genes through horizontal gene transfer (HGT). Multi drug resistant (MDR) bacteria are recognized as a major threat to public health [5] and require a comprehensive approach to combat them [6]. Resistance determinants are often present on mobile genetic elements (MGE), such as plasmids, integrons etc. and can be transferred among multiple bacterial genera [7, 8]. A wide range of plasmids that carry AR and virulence genes have been reported in Salmonella [911]. Here we catalog all the antibiotic resistance genes and their organization in a diverse set of Salmonella from swine.

Phage play a profound role in bacterial evolution as they assist in transfer of antibiotic resistance (AR) and, virulence genes, including inserting their genome into the host’s DNA [12]. These lysogens can also become lytic by replicating and then killing their temporary host bacteria. A large number of phage have been reported from Salmonella, such as Fels-1, Gifsy-2, P22, FelixO1, etc., and some of them carry several virulence and resistance gene cassettes [5, 1316]. Here, we investigate phage distribution among Salmonella serovars to determine the resistance and virulence genes associated with phage. Here we also characterize the CRISPR (clustered regularly interspaced short palindromic repeats) elements, which act as an adaptive immune system against exogenous DNA [17], including phage DNA. In addition, analysis of CRISPR sequences improve the discriminatory power of molecular characterization of Salmonella [18].

Few studies have been conducted to identify genetic factors associated with Salmonella serovars isolated from swine using WGS [19, 20]. In order to understand the genetic variations among different Salmonella serovars, the WGS of fourteen Salmonella serovars isolated from swine and swine swab were analyzed. The WGS data was assessed for resistance and virulence determinants, and their association with MGE was predicted. WGS was further analyzed to identify CRISPR elements and phage associated resistance and virulence genes.

Material and methods

Isolate selection and antimicrobial susceptibility testing

Fourteen Salmonella isolates were collected by the National Antimicrobial Resistance Monitoring System (NARMS) in year 2004–2005. All isolates were streaked onto Mueller-Hinton (MH) agar (Oxoid, Cambridge, UK). A single colony was selected and subsequently inoculated in MH broth (Oxoid, Cambridge, UK), and incubated for 16–18 h at 37°C with shaking (250 rpm). All isolates were subjected to susceptibility testing via the SensititreTM semi-automated antimicrobial susceptibility system (TREK Diagnostic Systems, Inc.) using a custom-made panel including amikacin, gentamicin, kanamycin, streptomycin, ampicillin, amoxicillin-clavulanic acid, ceftiofur, ceftriaxone, cefoxitin, sulfamethoxazole/sulfisoxazole, trimethoprim-sulfamethoxazole, chloramphenicol, ciprofloxacin, nalidixic acid, and tetracycline [21]. The isolates were subjected to preliminary biochemical screening to distinguish the different serogroups using serogroup-specific antisera (Difco Laboratories, Detroit, MI) and serotyping was used to identify serovars at the National Veterinary Services Laboratories, APHIS, USDA (Ames, IA).

Genome sequencing and analysis

Genomic DNA was isolated using GenElute bacterial genomic DNA kit (Sigma-Aldrich, St. Louis, MO) according to manufacturer's protocols. The quality of DNA was tested using Nanodrop and quantified using an Invitrogen Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA). The extracted DNA was stored in 1X TE buffer (pH 8.0) at -20°C until further use. Paired-end sequencing was performed on an Illumina HiSeq2500 (Illumina Inc., San Diego, CA) as described elsewhere [22] and the reads were deposited in the SRA database under the bioproject PRJNA254816. The SRA accession numbers is available in Table 1. The reads were retrieved from the SRA dataset and reads with a phred score ≥30 were de novo assembled using the A5-pipeline_version_20141120 [23]. The genome and plasmid contigs from the assembly were sorted on the basis of remote blast using an in-house python script. The genomic contigs were annotated using PROKKA [24]. In order to reconstruct the pan-genome the assembled genomes were annotated with Prokka and used as input for Roary with the identity cut-off of 95% [25]. Roary generated clusters of homologous gene groups from which core, accessory and unique genes were predicted. The Clusters of Orthologous Groups of proteins (COGs) database was used for the functional annotation [26]. The amino acid sequences generated from the Prokka was used as input for functional annotation based on orthologous group using WebMGA online server (http://weizhong-lab.ucsd.edu/metagenomic-analysis/server/cog/). The assembled contigs and deduced amino acid sequences from Prokka were used to predict the acquired antibiotic and metal resistance genes using the ARG-ANNOT [27] and BacMet [28] respectively, with e-value (1e -10), coverage ≥90% and identity ≥90%. The adherence-associated gene cluster were identified using the Virulence Factor Database (VFDB) [29] using a threshold of ≥95% identity and ≥95% coverage. Lysogenic phage were predicted using the web-based tool PHASTER [30]. CRISPR regions were predicted using CRISPR finder [31] and blastn was done to identify similar phage among Salmonella serovars and linear comparison of similar phage was done using Easyfig [32]. The Integrall integron database (http://integrall.bio.ua.pt) was used to analyze and assign integron sequences [33]. Genes associated with plasmid replicon were identified using PlasmidFinder to identify the target sequence in the genomes of each isolate [34].

thumbnail
Table 1. Assembly statistics of the Salmonella enterica serovars isolated from swine and their products.

https://doi.org/10.1371/journal.pone.0224518.t001

Results

This study was the part of retrospective study to maximize the understanding of the AR gene distribution, MGE and genome diversity of the Salmonella [22] from food animals. Fourteen Salmonella isolates from swine were selected based on differences in phenotypic AR profile, serotype from the collection of NARMS isolates between 2004 and 2005.

General features and sequenced genomes

The genome statistics are presented in Table 1. In brief, the average number of contigs per genome was 95 (range: 61 to 201 contigs). The median assembly coverage ranged between 53 to 95 fold which was adequate to produce bacterial draft genomes [35]. The average GC content of each genome was 52.1% which was consistent with that of the complete S. enterica chromosome [36]. The average number of coding sequences (CDSs) per isolate was 4,480 and the highest and lowest number of CDS was obtained for S. Agona (4,680) and S. Give (4,308), respectively. Additionally, average nucleotide identity matrix showed >98% identity among the serovars and >99.9% identity to their respective reference genome sequences from NCBI Gen-Bank. The high average nucleotide identity (ANI) [37] values of pair wise genome comparisons again confirmed that these serovars were nearly identical to their corresponding reference serovars. The reads were deposited in the GenBank Sequence Read Archive (SRA) database and details are provided in Table 1.

Antimicrobial resistance genes, integrons, and plasmids

Antibiotic susceptibility assays were performed according to Clinical and Laboratory Standards Institute (CLSI) standards. Eleven isolates were resistant to at least one of the tested antibiotics and nine of them were multidrug resistant (MDR; resistant to two or more antimicrobials); the remaining three Salmonella serovars Give, Infantis, and Manhattan were susceptible to all tested antimicrobials. None of the isolates were resistant to azithromycin, ciprofloxacin or nalidixic acid. The susceptibility testing results and corresponding predicted AR genes results are summarized in Table 2. High frequency of resistance to tetracycline (9/14, 64.3%) was observed among the isolates, followed by resistance to beta-lactams (7/14, 50.0%), aminoglycosides (7/14, 50.0%) and sulfonamides (7/14, 50.0%). Three isolates 3/14, 21%). (S. Agona, S. Ohio and S. Typhimurium) were resistant to two antibiotics (chloramphenicol and sulfamethoxazole) while two isolates (2/14, 14.3%); S. Anatum and S. Worthington) were only resistant to tetracycline. The correlation between phenotypically confirmed AR and in silico predicted AR genes of the serovars is presented in Table 2. Resistance determinants were predicted with ARG-ANNOT. Hits with > 90% identity and 100% coverage were considered as positive resistance determinants. A total of 73 AR genes (59 in Table 1 and 14 aac6-I variants) were predicted from the genome sequences and at least one AR gene was predicted in each isolate (Table 1). The aminoglycoside resistance gene variants Y and aa of aac6-I was detected in all the isolates and these cryptic variants were the only aminoglycoside resistance gene variants predicted in Salmonella serovars Give, Infantis, and Mannhattan exhibited no phenotypic resistance. Among eight streptomycin resistant isolates strA and strB were detected in five isolates while four isolates carry aadA gene variants and all these three genes (strA, strB, and aadA) were detected in S. Heidelberg. Seven isolates were resistant to beta-lactams and three genes encoding beta-lactamases were identified in these seven isolates, with the most common being blaTEM-1 ((3/14), 21.4%) and blaCMY-2 ((3/14), 21.4%), followed by blaCARB-2 ((2/14), 14.3%). The Amber class A β-lactamase gene (blaCARB-2) conferring resistance to ampicillin were detected in serovar Tennessee and Typhimurium. Amber classA potential ESBL (extended spectrum of beta-lactamase) gene (blaTEM-1) conferring resistance to ampicillin was detected in serovars Bovismorbificans, Heidelberg and Ohio. Amber class C genes (blaCMY-2) that conferred resistance to ampicillin, amoxicillin, and ceftiofur were detected in serovars Agona, Minnesota and Ohio. Four different tetracycline resistance gene allele’s tetA, tetB, tetC and tetG were identified in the analyzed Salmonella serovars. With tetA and tetB the most frequently detected tetracycline resistances, occurring in 4/14, 28.6% and 3/14, 21.4% of isolates, respectively, followed by tetC and tetG each in 1/14, 7.1% of isolates. The macrolide resistance gene ereA was detected in Tennessee, while trimethoprim gene variants (dfrA1, A7, A12, and A16) were detected in serovars Ohio, Agona, Heidelberg, and Tennessee respectively. The sulfonamide resistance gene alleles were identified in seven Salmonella serovars; only sul1 was detected in serovar Bovismorbificans, Choleraesuis, Cubana and Typhimurium ((4/14); 28.6%), while sul1 and sul2 were detected in serovar Agona and Ohio ((2/14); 14.3%) and sul3 was detected in serovar Heidelberg ((1/14); 7.1%). The floR gene was detected in three of fourteen chloramphenicol resistant isolates and only these three serovars Agona, Ohio, and Typhimurium conferred resistance to chloramphenicol. No mutations in chromosomal genes gyrA, gyrB, parE: conferring resistance to ciprofloxacin and nalidixic acid and rRNA genes 23S rRNA, rplD, rplVb conferring resistance to macrolide were observed; a point mutation in parC at position T57S was observed in most of the studied serovars except Bovismorbificans and Typhimurium; however, none of these isolates were resistant to nalidixic acid (S2 Table).

thumbnail
Table 2. Antibiotic resistance phenotype and predicted antibiotic resistance and virulence genes in different Salmonella serovars.

https://doi.org/10.1371/journal.pone.0224518.t002

Class 1 integrons were detected in six isolates that were grouped into five different integron profiles (ln127, in167, in363, in1581, and ln1582). Most of these integrons carried a quaternary ammonium compound resistance gene, qacE, and a sulfonamide resistance gene, sul1, as additional genes after the 3`end of the conserved segment. The class 1 integron genesin different Salmonella serovars are summarized in Table 3.

thumbnail
Table 3. Class 1 integrons identified in different Salmonella isolates.

https://doi.org/10.1371/journal.pone.0224518.t003

The presence of plasmids in the assembled contigs were confirmed by in silico replicon typing. The replicon type IncA/C2 was identified in MDR serovars Agona and Ohio and both these isolates harbored identical class 1 integrons In363. The replicon type IncHI2 was detected in serovars Bovisformicans and Heidelberg while IncI1 was identified in serovars Infantis, Minnesota and Worthington. The replicon type IncFII was identified in serovars Choleraesuis, Heidelberg and Typhimurium while IncN and IncP6 were detected in serovars Tennessee, and Minnesota, respectively. No replicon types were identified in serovars Cubana, Give, or Manhattan. Additional replicon col440II, and IncQ was identified in S. Choleraesuis. Eight different resistance gene clusters were detected in different Salmonella serovars and are presented in Fig 1.

thumbnail
Fig 1. Antibiotic resistance gene clusters in different Salmonella serotypes.

Bleomycin and EtBr resistance gene (Green) flanked by antibiotic resistance genes marked in green. The arrow shows orientation of the genes and genes are color coded to define different classes of antibiotic resistance, mobile and other genes categories.

https://doi.org/10.1371/journal.pone.0224518.g001

The AR genes clusters were highly diverse, with the exception of a > 35 kb (99.9% identical) region of homologous genes carrying tetA, strA, strB, sul2, and blaCMY-2, which were detected in Salmonella serovar Agona and Ohio. The region of homologous genes carrying strA/strB, and strA/strB/sul2 were detected in Salmonella serovar Bovismorbificans and Choleraesuis respectively, while aadA/ebr/sul1 and ereA/aadA2/blaCARB-2/dfrA16 were detected in Salmonella serovar Cubana and Tennessee respectively. The gene clusters aphA2/sph, and sul1/floR/tetA were detected in Salmonella serovar Heidelberg and Typhimurium respectively.

The pan-genome and functional genes comparison

Pangenome analysis was initiated using 62,730 genes of 14 Salmonella serovars that resulted into 8,174 clusters, of these 3,456 (42.3%) and 2,201 (26.9%) clusters were part of the core and accessory gene respectively. The remaining 2,517 (30.8%) independent genes identified as the unique genes. The core, accessory and unique genes are represented as inner circle, outer circle, and petals in the floral diagram respectively (Fig 2) and the details of core, accessory and unique genes in each serovars is given in S1 Fig.

thumbnail
Fig 2. Floral venn diagram showing the pangenome of 14 Salmonella serovars from swine.

The orthologus genes identified in all serovars presented in the center as core genes, orthologus genes identified among the serovars but not in all serovars presented in the periphery as accessary genes and each petals represents the unique genes in respective serotypes.

https://doi.org/10.1371/journal.pone.0224518.g002

The unique genes among serovars ranged from 93 to 268, with highest and lowest genes in serovars Worthington and Bovismorbificans respectively. The functional distribution of genes among serovars were examined using the COG database [26]. Functions encoded by the genes in these serotypes revealed > 75% of predicted ORFs were assigned to 24 COG functional group (Fig 3).

thumbnail
Fig 3. Distribution of functional classes of predicted genes according to the clusters of orthologous groups in Salmonella serovars.

Different colors define different COG categories of genes.

https://doi.org/10.1371/journal.pone.0224518.g003

The most abundant COG functional group among serovars was the transport and metabolism of carbohydrates (G) and the highest (426) and lowest number of genes (385) in this category was noticed in S. Typhimurium and S. Minnesota respectively. The next abundant COGs were Transcription (K) followed by amino acid metabolism and transport (E). The highest number of genes in K and E category were noticed in S. Worthington (346) and S. Ohio (381) respectively, however the lowest number of genes in K (322) and E (370) COGs were noticed in S. Minnesota. A single COG functional gene was noticed in all the S. serotypes from extracellular structure categories. The distribution of functionally characterized COG genes in core and accessory genome in Salmonella serovars (Fig 4) revealed that the functional genes in the core genome ranged 76% to 82%.

thumbnail
Fig 4. Distribution of clusters of orthologous groups in core, accessory and unique genes in the genomes of different Salmonella serovars.

Common numbers of core genes in each COG categories were observed and presented in Fig 4(A), The bars represents the number of clusters of orthologous groups assigned genes present in accessory (grey bars) and in the unique genes (black bars).

https://doi.org/10.1371/journal.pone.0224518.g004

Uniform distributions of functional genes from different COG category was observed for core and accessory genes among Salmonella serovars. However, the percent abundances of genes in category G (carbohydrate transport and metabolism; 26–28%) and L (Replication, recombination and repair 21–23%) among accessory genomes were greater among Salmonella serotypes compared to other functional categories. The core genomes in each S. serotype were commonly enriched in COG categories D, F, H, I, J, O and Q relative to those seen in the accessory genomes.

Virulence genes, phage and CRISPR in Salmonella serovars

Genes with >95% sequence homology and coverage when compared to the virulence factor database (VFDB) were considered matches in this analysis. The complete virulence gene profiles of each Salmonella isolate is shown in S1 Table and some of the gene variants that were not detected in all the isolates is shown in Table 2. The majority of the virulence genes were detected in all the Salmonella serovars. These virulence genes were very similar, with small variations in nucleotide/amino acid sequences, while few virulence genes were missing in single or multiple Salmonella serovars. The fimbrial adherence gene operons bcfA/B/C/D/E/F/G, fimI/C/D/H, invA/B/C/E/F/G/H/I/J, csgBA and csgD/E/F/A were the most common. Most of the isolates were positive for the long polar fimbrial operon (lpfA/B/C/D/E) except Salmonella serovars Anatum, Give, Minnesota, and Ohio; however, lpfD was highly diverse among all of the Salmonella isolates. The plasmid-encoded fimbrial operon (pefA/B/C/D) and Salmonella plasmid virulence (spvB/C/R) genes were detected in S. Typhimurium and none of these genes were detected in other isolates. The pathogenicity island 1 (SPI-1), encoding type III secretion system (T3SS) secreted effector genes sipA/B/C/D, sopA, sopB, sopE2, and misL, were highly similar in all the isolates with the exception of sipD genes that shared 91% similarity to the reference gene (S. Typhimurium; NP_461804) in S. Minnesota. AvrA was identified in all the isolates except S. Choleraesuis and S. Ohio. The SPI-2 encoded T3SS genes, including spiC, sifA, sifB, sseF, sseG, sseL, pipB, pipB2, sopD2, and slrP, were detected in all the isolates. GogB was only detected in Salmonella serovars Choleraesuis and Typhimurium. Other effector genes, including ratB, sseK1, sseK2, and sspH2, were detected in some of the Salmonella serovars (S1 Table). The sspH2 effector homolog was divergent (85.52% identical) in S. Ohio. The sensory systems genes phoP/Q were identical in all the serovars with an exception in S. Worthington where a secondary mutation was observed at nucleotide position 450 (A-T) with coverage depth >50x.

Genome contents were further compared, and it was observed that several genes were confined to one serovar, or were highly diverse even if they were annotated as the same gene (S3 Table). For example a 27 kb conting in S. Agona harboring subset of type VI secretion system (T6SS). A unique pathogenesis protein, kcpA, was identified on a phage in S. Anatum. Homologs of the safA (saf-pilin pilus formation protein) was identified in 6 of the 14 Salmonella serovars (Anatum, Bovismorbificans, Cubana, Give, Johanesberg and Manhattan), but these homologs were so diverse that they were included in the list of unique genes. Highly diverse toxin-antitoxin (TA) systems type II genes were detected in different Salmonella serovars, including vapB (Agona and Typhimurium), higA (Choleraesuis and Cubana), yfjZ (Minnesota and Manhattan) and cbtA (Bovismorbificans and Tennessee). UDP-L-Ara4N formyltransferase (arnA), a bifunctional genes facilitating polymyxin resistance, was detected in S. Choleraesuis. Highly diverse colonization factor antigen I subunit E (cfaE) gene was detected in Salmonella serovar Choleraesuis, Give, Minnesota and Worthington. Microcin-M immunity protein (cmi) was detected in S. Heidelberg.

The PHASTER analysis identified intact, questionable, and incomplete phage in genomes of the Salmonella isolates. Only intact phage were analyzed in this study and a total of 52 intact lysogenic phage were predicted by PHASTER and they are listed in Table 4. Most phage were highly diverse. The most common families included relatives of Gifsy-1 (NC_010392), Gifsy-2 (NC_010393), and Salmon_SEN34 (NC_028699). Gifsy-1 like phage were detected in Salmonella serovars Bovismorbificans, Infantis, Typhimurium, and Worthington, Gifsy-2 like phage were detected in Salmonella serovars Choleraesuis, Heidelberg and Typhimurium, and the Salmon_SEN34 like phage were detected in Salmonella serovars Choleraesuis, Cubana and Give. Individual phage found in only one of the 14 included in Cubana (Entero_Psp3 [NC_005340]), Heidelberg (Escher_D108 [NC_013594]; Entero_P4 [NC_001609]), Minnesota (Entero_I2_2 [NC_001332]; Entero_SfV [NC_003444]; Salmon_vB_SemP_Emek [NC_018275]), Manhattan (Salmon_epsilon34 [NC_011976]; Entero_186 [NC_001317]) and Typhimurium (Entero_ST104 [NC_005841]). Of particular interest is the fact that some of the phage carried non-phage “cargo” genes. The AR beta-lactamase gene blaTEM-1 was seen in S. Heidelberg (Salmon_118970_sal4_NC030919), S. Ohio (Entero_lato_NC001422). The blaCARB-2 gene was seen in S. Tennessee (Entero_lato_NC001422) and S. Typhimurium (Salmon_119870_sal3_NC031940). The tetracycline gene tetC was seen in S. Anatum (Salmon_118970_sal3_NC031940); and tetA was seen in S. Cubana (Entero_lato_NC001422). The aminoglycoside gene aadA2 was seen in S. Tennessee (Entero_lato_NC001422) and S. Typhimurium (Salmon_118970_sal3_NC031940). Our analysis further confirmed that blaCARB-2 was harbored on the class 1 integrons In1582 and In167 that integrated into the possible phage Entero_lato_NC001422 (S. Tennessee) and Salmon_119870_sal3_NC031940 (S. Typhimurium), respectively. We also observed that metal resistance genes merDACPTR and silA/C/P/P2 were on the possible prophage Entero_lato_NC_001422 and Salmon_118970_sal4_NC030919 of S. Cubana and S. Heidelberg, respectively (Fig 5).

thumbnail
Fig 5. Examples of integrated bacteriophage.

Cargo genes are annotated, such as antibiotic and other resistance genes (Hg resistance gene [merA/C/P/T/ and mer regulators merD/R; S. Cubana], Ag resistance genes [silA/C/P/P2, SRP; S. Heidelberg], quaternary ammonium compound gene [qacC; S. Tennessee and S. Typhimurium]). The arrow shows orientation of the genes in the contigs and genes are color coded to define different gene categories.

https://doi.org/10.1371/journal.pone.0224518.g005

thumbnail
Table 4. Phage harboring virulence and resistance genes in Salmonella serovars.

https://doi.org/10.1371/journal.pone.0224518.t004

These phage were further analyzed for the presence of virulence genes and results are summarized in Table 4. Eight of fourteen isolates harbored at least one virulence gene on phage. The comparisons for the homologous regions in prophage revealed 10 different phage shared similar regions, although these prophage were not 100% identical, but showed similarity ranging between 40 to 98%. The linear comparisons of the phage are presented in Fig 6. The identical regions were mostly part of the phage related genes.

thumbnail
Fig 6. Alignment of closely homologous intact integrated phage among Salmonella serotypes.

A. Phage genes synteny in five serotypes. B. Phage gene arrangement along with mobile genes (transposase, recombinase) in six serotypes. The arrow shows orientation of the genes and genes are color coded to define mobile, phage and other genes.

https://doi.org/10.1371/journal.pone.0224518.g006

CRISPR loci were identified in the majority of isolates and these loci were compared for homology across serovars. The common spacers were observed towards the ancestral end of the CRISPR array e.g. common spacer was identified in eight serovars (Agona, Cubana, Choleraesuis, Heidelberg, Infantis, Manhattan, Typhimurium and Worthington). Spacer numbers varied between serovars and across the CRISPR loci (Table 5).

thumbnail
Table 5. Prophage target identified in spacer sequences in CRISPR arrays in Salmonella serovars.

https://doi.org/10.1371/journal.pone.0224518.t005

The sequential identical spacers towards the ancestral end was identified among these serovars for eg. two spacers detected in Salmonella serovar Cubana and Worthington, three spacers detected in Salmonella serovar Agona, Heidelberg and Typhimurium, four spacers detected in Salmonella serovar Heidelberg and Typhimurium. Unique spacer arrays were observed in serovar Minnesota, Ohio and Tennessee (the last spacer i.e. spacer 57 was identical to spacer 2 of serovar Anatum). The signature protein of the Type I CRISPR systems Cas3 gene was detected in all the S. serovars and variations were noticed in this gene when compared with S. typhimurium LT2 cas3 as summarized in Table 5. Spacer sequences were also analyzed to detect target prophage with the help of Blast using UniProt phage sequences (https://www.uniprot.org/). Prophage targets were observed in most CRISPR loci, with the most common prophage targets being tail, head, and capsid proteins. Besides this, the other targets were integrase and recombinase regions.

Discussions

The implementation of WGS allows broader inference of pathogen characteristics including prediction of antibiotic resistance and virulence profiles from the sequences. WGS has been previously used for the prediction of AR genes in a wide range of microbes including Salmonella [38]. Phenotypic and genotypic correlation analysis confirmed high concordance (97%) between phenotypically confirmed and in silico predicted AR genes. No absolute concordance was observed due to insertional inactivation if tet(B) gene in Choleraesuis and no resistance phenotype was shown by this isolate (Table 2). Resistance to aminoglycosides is either due to inactivation or modification by acetyltransferases, phosphotransferases, and nucleotidyltransferases [39]. The aac6-I variants Y and aa of was detected in all the isolates and located on the chromosome. This cryptic gene has previously been reported in Salmonella that was due to a deletion in the promoter region [40]. The strA/B and aadA genes are frequently associated with MGE that easily disseminate aminoglycoside resistance genes in Salmonella and other gram negative bacteria; streptomycin resistance has also been used as an important epidemiological marker to indicate the possibility of MDR in pathogens [4143]. Many class 1 integrons harbor aadA gene variants [44]; four aadA variants were associated with class 1 integrons in this analysis. Other aminoglycoside resistance genes that confer resistance to kanamycin (aph(3”)-Ia) and hygromycin B (aph(4’)-Ia) were predicted in this analysis and these genes have also been previously identified in WGS data of Salmonella; moreover, aph(3”)-Ia was the most frequently detected gene from food sources [19]. All three β-lactam resistance genes (blaTEM-1, blaCMY-2 and blaCARB-2) are highly prevalent in Salmonella isolated from U.S. animals, and humans. These genes are the mostly horizontally acquired β-lactamases in Salmonella [42, 45]. Tetracycline is widely used in food animals to combat respiratory infections and have also been used as growth enhancer and an additive in feed [46](45,46). The genes conferring resistance to tetracycline are widespread among Salmonella serovars and easily transferred among wide range of microbes through HGT [47, 48]. All isolates except one with tet genes were resistant to tetracycline in this analysis and only tetracycline efflux pump genes were identified. No other tet genes (ribosomal protection protein or inactivated enzyme) were detected which suggested that the efflux pump mediated resistance was the major tetracycline resistance mechanism among these Salmonella serovars. The high prevalence of these efflux pump genes have been reported in Salmonella species [49] which further confirmed our findings (Table 2). The erythromycin resistance gene ere(A) was detected on a class 1 integrons that harbored dfrA, blaCARB-2 and aadA2 resistance determinants (Table 3). The ere(A) gene has previously been reported in class 1 integrons [50]. All dfr variants (A1, A7, A12 and A16), are associated with different profiles of class 1 integrons encoding resistance to trimethoprim. The class 1 integron carrying dfrA gene variants and other resistance genes have been reported from several Salmonella serovars Anatum, Choleraesuis, Corvallis, Eppendorf, Gallinarum, Kentucky, Rissen, Stanley, Schwarzenrund, Typhimurium, and Weltevreden [5154]. Identical class 1 integron gene cassette (In363) were detected in Salmonella serovar Agona and Ohio (Table 3), that indicated that inter or intra species transfer of integrons help in the spread of antimicrobial resistance genes among bacteria. The fosfomycin resistance fosA gene alleles are located on the chromosome.The floR genes has been isolated from a wide range of animals and humans and it has been previously reported from various Salmonella serovars including S. Typhimurium. Chloramphenicol was used to treat MDR infections in human and animals [5557]. The sulfisoxazole resistance genes sul1 (located within 3`-conserved segment [3`-CS] of class 1 integrons) and sul2 (associated with small multicopy plasmids or large transmissible multiresistance plasmids) are the most frequently found genes for sulfonamide resistance among sulfonamide-resistant isolates from food animals and humans, whereas sul3 is detected in various large Salmonella plasmids [58, 59]. The sulfonamide resistance genes sul1 and sul2 were the most abundant AR genes detected in this analysis, which was in accordance with the previous study where high prevalence of these genes has been reported in the MDR Salmonella isolated from animals, humans and retail meats from Northern America [60]. A secondary mutation in parC without primary mutation in gyrA is ineffective to make these isolates resistant to nalidixic acid; however, these isolates may become highly resistant once primary mutation occurs in gyrA gene as reported previously where several fold increase in resistance was observed when a secondary mutation in parC was seen in an isolate that had a primary mutation in gyrA [61, 62].

The AR genes carrying contigs were distinct and none of them were identical among Salmonella serovars except AR contigs of S. Agona and S. Ohio that shared 99.9% sequence similarity and gave best hit with virulence-resistance plasmids (incA/C2) of S. Typhimurium [63]. This >35.0 kb fragment carried tetA, tetR, strA, strB, sul2 and blaCMY-2. The strA, strB, sul2 AR genes was flanked by Tn7 elements in this partial incA/C2 plasmid and this was the most frequently recorded resistance gene combination among Salmonella serovars located on small broad-host-range plasmids, as well as detected on the chromosome in Salmonella [49]. Another 5.0 kb fragment from S. Bovismorbificans carrying strA and strB best matched with the genome of S. Enteritidis (CP020442 and CP022069) and showed 99% coverage and 100% identity to the region and appear to be widely distributed in Salmonella and other gram-negative bacteria [49]. These genes have been described as being part of transposon Tn5393 and have also been identified in bacteria circulating in humans, animals, and plants [64]. Another contig carrying aadA, ebr and sul1 was detected in S. Cubana and gave best hit with E. coli plasmid pBM0133 (KJ170699)/ pDGO101 and this array of genes is harbored by class 1 integron In1581 (Table 3) [65, 66]. The contig from S. Heidelberg that harbored neomycin (aphA2), bleomycin and streptomycin (sph) resistance genes best matched with the E. coli transposons Tn5 (U00004). These genes (aphA2 and sph) are encoded by Tn5 and used as selectable markers in cloning vectors for both eukaryotes and prokaryotes [67, 68]. A 8.69 kb contig from S. Choleraesuis carrying strA, strB and sul2 gave best hit with TY474p3 plasmid (CP002490) and it was a closed plasmid [49]. Another 5.5Kb contig carrying ereA, aadA2, blaCARB-2 and dfrA16 from S. Tennessee gave best hit with E. coli integron intl1 (KX57988). This integron carrying ereA and aadA2 is rarely reported [69], while a 5.0Kb fragment carrying floR and tetA in S. Typhimurium, this fragments was the part of Salmonella genomic island 1 (SGI1) that best matched with the regions of complete genome of Salmonella Typhimurium (CP028318:4845261–4850261, and CP014979:42459–47459). SGI1 has been identified in several Salmonella serovars including Typhimurium, Agona, Paratyphi B, Albany, Meleagridis and Newport [7074]. The detection of SGI1 in different Salmonella serovars from various sources including animals and humans across the world indicates frequent dissemination of the SGI1 through horizontal transfer [75].

The size of the core genome in this analysis was in accordance with the previous studies where the size of the pan-genome expanded slightly while the size of core genome shrunken with addition of new genomes. The inter serovars core genome size (3, 224 core genes) within the 35 Salmonella subsp. enterica was lower while the intra serovar core genome size (3, 836 core genes) for S. Typhimurium was higher compared to our analysis [76, 77]. Low number of unique genes was directly related to the reduced genomes size, where HGT events were less and S. Bovismorbificans well correlated in this context with lowest numbers of unique genes and lower genome size. We could not draw any correlation between high number of unique genes and with large genome size. However, moderate genome length with high unique genes revealed its evolutionary significance. Uniform distribution of genes was noticed for other COG categories among S. serotypes. The abundance of E, G, and K cog categories of genes increases the diversity, uniqueness and their complex transcriptional regulatory networks that support morphological and physiological differentiation. The abundance of G and L functional categories of genes in dispensable genes clearly showed that these functional genes acquired for the better adaptability and diversity [78].

The majority of virulence genes were part of the core genes and essential for pathogenicity and infection. WGS data helped to correctly interpret the variations in genes among Salmonella serovars for better understanding of pathogen and evolution. This analysis clearly identified some of the genes specifically detected in one or several Salmonella serovars. For example, the pef operon comprising pefABCD genes that is needed to form structural fimbria and mediate the binding of bacteria to the microvilli of enterocytes and spvBCR that enhance virulence were solely identified in S. Typhimurium; these operons were not detected in S. Choleraesuis. However, they have been reported in virulence plasmid pSCV50 [79] harbored in S. Choleraesuis. RatB which encodes a secreted protein associated with intestinal colonization and persistence was detected in all the analyzed serovars except Agona, Anatum, Cubana, Tennessee, and Worthington [80] while gogB an anti-inflammatory effector that limits tissue damage during infection was detected in serovars Typhimurium and Choleraesuis which are predominantly associated with swine [81, 82]. Several studies have demonstrated the loss of genes associated with host specificity, thus some microbes only adapted to become a better pathogen to specific hosts while losing the ability to infect other potential hosts [83]. The lpf operon has been reported in Salmonella and is involved in the adhesion to the small intestine [84]. This operon was detected in most Salmonella serovars with high variations in lpfD genes that are highly host specific; variations in these genes could be an important factor that may influence the host range. The conserved nature of phoP/Q genes suggests that its pivotal regulatory role in Salmonella promotes several phenotypes including cationic antimicrobial peptides resistance which enhance outer membrane barrier function, and the ability to cause invasive and systemic disease [85, 86].

The genes which are confined to one serovar or their highly variable variants detected among serovars were mostly acquired genes and many of them were harbored on phage or plasmid. The T6SS is a versatile secretion system widespread among Gram-negative bacteria and directly involved in a variety of cellular processes, including virulence in several bacterial pathogens [87, 88]. Homology search resulted in best hits with Salmonella serovars Weltevreden (LN890520), Agona (CP015024), and Sloterdijk (CP012349). These loci have been reported to encode T6SS harbored on SPI-19 in Agona and Weltevreden [89]. safA plays a critical role in host recognition; highly diverse safA genes in six different Salmonella serovars confirmed that these serovars could be the part of swine asymptomatic carriers. Whereas highly identical (>98%) safA was found in S. Choleraesuis and S. Typhimurium which may recognize similar hosts as these two serovars are predominantly associated with swine [90]. The cfa gene in S. Choleraesuis was highly conserved and were identical with other S. Choleraesuis in Gen-Bank database (CP012344 and CP007639) and it may be essential for the attachment to a specific host [91]. The antimicrobial peptides encoding cmi from S. Heidelberg best matched with Escherichia coli plasmids harboring this gene (CP035337 and CP035355) and these plasmids harbored several AR genes. We believe cmi was transferred to S. Heidelberg as an accessory gene with other AR genes.

The integration of phage genetic elements in the genomes results in genome diversification of closely related bacterial strains [92] as these phage harbors AR and virulence factor that can enhance bacterial fitness [93]. Phage-associated resistance and virulence genes have been reported in Salmonella revealing the role of phage-transferred genes in pathogenicity and resistance [94, 95]. Gifsy-1 and Gifsy-2 like phage were highly prevalent among Salmonella serovars and these phage have previously been reported in Salmonella [96, 97]. In this study, we found a number of phage that encode for resistance and virulence genes. For example, we identified a 75.1 kb phage that harbor class 1 integrons carrying aadA2, blaCARB-2, sul1 AR genes and encodes a glucosyl transferases gtrA in the serovar Typhimurium, this putative virulence gene has previously been reported in P22-like and P2-like prophage in different Salmonella serovars [98].

The analysis and interpretation of the CRISPR region is complex. Spacers can only be added to the 5′ end and deleted anywhere in the locus, leaving the 3′-spacer as ancestral. The conserved spacers at the ancestral ends allowed us to understand the common infection history or common ancestory of these serovars though there is a possibility of the degradation of many internal spacers [17]. A careful analysis of CRISPR locus alignments provided an evolutionary framework to view the phylogenetic patterns of CRISPR diversity as we observed deletion of spacers in serovars Choleraesuis, Heidelberg and Typhimurium. Prophage targets were observed in most CRISPR loci; and about 87% (33/38) of spacer sequences were found to match with sequences from phage suggesting that phage were an important type of HGT in Salmonella. In this study, all eight cas genes pattern was not identified in S. Worthington as draft genome sequences was used and missing sequencing data cannot be ruled out. Cas3 involved in the cleavage of invading DNA and considered as an important component of the CRISPR mechanism and it was the only cas gene detected in all the serovars [99]. Comparison of cas3 among the serovars confirmed two distinct pattern,that was in accordance with earlier report where two distinct cas3 pattern was reported [100].

Vehicles of dissemination (plasmid, phage, integrons etc.) are crucial for microbial evolution and they play a major role in genetic innovation and genome evolution. Resistance determinants were located either on plasmids, integrons or integrons were also associated with phage which indicated that they were acquired horizontally. Highly diverse class 1 integrons containing different AR gene cassettes were detected among Salmonella serovars and these AR genes accounted for resistance to relative drugs. This suggests that these determinants of AR are functional and crucial for the survival of the microbes. A fraction of the virulence determinants were also identified on mobile elements and this suggests that these determinant are crucial for better survival. The common ancestral ends in CRISPR suggest that these isolates shared common ancestry and microbes are continuously acquiring new phage to counter the CRISPR regulation.

Conclusions

The insight into the WGS data enabled us to virtually access the genetic content for the assessment of genetic diversity, pathogenesis, evolution, serotyping, virulence and resistance profiling that enhanced our understanding of the genomic diversity among Salmonella serovars. The pangenome size was in accordance with the earlier report however the core genome was comparatively higher due to less isolates used in this analysis. Most of the virulence factors were part of the core genomes indicating absolutely required factor for virulence; however, a small fraction of virulence factors were highly diverse (safA) or confined to one isolates (cfa, pefABCD, and spvBCF) indicating variation in pathogenicity among Salmonella serovars. Inter and the intra species transferability of MGE enhance their reachability to the specific or the broader host. Variable MGE information from the WGS enables understanding of the dynamics of horizontal gene transfer carrying resistomes that were widely distributed in other Gram-negative bacteria. Detection of class 1 integrons carrying resistance determinants (In1582 and In167), high sequence diversity among related prophage carrying various resistance and virulence genes, suggesting extensive horizontal gene transfer that enhanced bacterial fitness. WGS enabled us to obtain a comprehensive resistome profile that can be useful to develop antibiotic resistance combat strategies. Comprehensive genome analysis with large dataset could profoundly enhance our understanding of the adaptability and survival mechanism in Salmonella serovars.

Supporting information

S1 Table. Virulence profiles of Salmonella serovars.

https://doi.org/10.1371/journal.pone.0224518.s001

(XLSX)

S2 Table. Analysis of mutations in Chromosomal genes conferring resistance to antibiotics in Salmonella serovars.

https://doi.org/10.1371/journal.pone.0224518.s002

(XLSX)

S3 Table. Gene content in Salmonella serovars confined to one serovars.

https://doi.org/10.1371/journal.pone.0224518.s003

(XLSX)

S1 Fig. Bar graph represents the distributions of the genes in core, accessory and unique gene in different Salmonella serovars from swine.

https://doi.org/10.1371/journal.pone.0224518.s004

(TIF)

Acknowledgments

SKG thanks Calvin Williams for IT support at USNPRC facility. MM thanks Weiping Chu, Steffen Porwollik, and Prerak Desai, who were involved in genomic DNA assessment, sequence storage, and confirmation.

References

  1. 1. Chan K, Baker S, Kim CC, Detweiler CS, Dougan G, Falkow S. Genomic comparison of Salmonella enterica serovars and Salmonella bongori by use of an S. enterica serovar typhimurium DNA microarray. J Bacteriol. 2003;185(2):553–63. pmid:12511502; PubMed Central PMCID: PMC145314.
  2. 2. Scallan E, Hoekstra RM, Angulo FJ, Tauxe RV, Widdowson MA, Roy SL, et al. Foodborne illness acquired in the United States—major pathogens. Emerg Infect Dis. 2011;17(1):7–15. pmid:21192848; PubMed Central PMCID: PMC3375761.
  3. 3. Minor T, Lasher A, Klontz K, Brown B, Nardinelli C, Zorn D. The Per Case and Total Annual Costs of Foodborne Illness in the United States. Risk Anal. 2015;35(6):1125–39. pmid:25557397.
  4. 4. Guibourdenche M, Roggentin P, Mikoleit M, Fields PI, Bockemuhl J, Grimont PA, et al. Supplement 2003–2007 (No. 47) to the White-Kauffmann-Le Minor scheme. Res Microbiol. 2010;161(1):26–9. pmid:19840847.
  5. 5. Davies J, Davies D. Origins and evolution of antibiotic resistance. Microbiol Mol Biol Rev. 2010;74(3):417–33. pmid:20805405; PubMed Central PMCID: PMC2937522.
  6. 6. Davies SC, Fowler T, Watson J, Livermore DM, Walker D. Annual Report of the Chief Medical Officer: infection and the rise of antimicrobial resistance. Lancet. 2013;381(9878):1606–9. pmid:23489756.
  7. 7. Fricke WF, Mammel MK, McDermott PF, Tartera C, White DG, Leclerc JE, et al. Comparative genomics of 28 Salmonella enterica isolates: evidence for CRISPR-mediated adaptive sublineage evolution. J Bacteriol. 2011;193(14):3556–68. pmid:21602358; PubMed Central PMCID: PMC3133335.
  8. 8. Carattoli A. Plasmids and the spread of resistance. Int J Med Microbiol. 2013;303(6–7):298–304. pmid:23499304.
  9. 9. Chu C, Hong SF, Tsai C, Lin WS, Liu TP, Ou JT. Comparative physical and genetic maps of the virulence plasmids of Salmonella enterica serovars typhimurium, enteritidis, choleraesuis, and dublin. Infect Immun. 1999;67(5):2611–4. pmid:10225928; PubMed Central PMCID: PMC116011.
  10. 10. Garcia P, Guerra B, Bances M, Mendoza MC, Rodicio MR. IncA/C plasmids mediate antimicrobial resistance linked to virulence genes in the Spanish clone of the emerging Salmonella enterica serotype 4,[5],12:i. J Antimicrob Chemother. 2011;66(3):543–9. pmid:21177672.
  11. 11. Hoffmann M, Pettengill JB, Gonzalez-Escalona N, Miller J, Ayers SL, Zhao S, et al. Comparative Sequence Analysis of Multidrug-Resistant IncA/C Plasmids from Salmonella enterica. Front Microbiol. 2017;8:1459. pmid:28824587; PubMed Central PMCID: PMC5545573.
  12. 12. Suttle CA. Marine viruses—major players in the global ecosystem. Nat Rev Microbiol. 2007;5(10):801–12. Epub 2007/09/15. pmid:17853907.
  13. 13. Boyd EF, Brussow H. Common themes among bacteriophage-encoded virulence factors and diversity among the bacteriophages involved. Trends Microbiol. 2002;10(11):521–9. Epub 2002/11/07. pmid:12419617.
  14. 14. Hooton SP, Timms AR, Rowsell J, Wilson R, Connerton IF. Salmonella Typhimurium-specific bacteriophage PhiSH19 and the origins of species specificity in the Vi01-like phage family. Virol J. 2011;8:498. Epub 2011/11/04. pmid:22047448; PubMed Central PMCID: PMC3220722.
  15. 15. Kropinski AM, Sulakvelidze A, Konczy P, Poppe C. Salmonella phages and prophages—genomics and practical aspects. Methods Mol Biol. 2007;394:133–75. Epub 2008/03/28. pmid:18363236.
  16. 16. Price-Carter M, Roy-Chowdhury P, Pope CE, Paine S, de Lisle GW, Collins DM, et al. The evolution and distribution of phage ST160 within Salmonella enterica serotype Typhimurium. Epidemiol Infect. 2011;139(8):1262–71. Epub 2010/10/19. pmid:20950514.
  17. 17. Barrangou R, Marraffini LA. CRISPR-Cas systems: Prokaryotes upgrade to adaptive immunity. Mol Cell. 2014;54(2):234–44. Epub 2014/04/29. pmid:24766887; PubMed Central PMCID: PMC4025954.
  18. 18. Horvath P, Romero DA, Coute-Monvoisin AC, Richards M, Deveau H, Moineau S, et al. Diversity, activity, and evolution of CRISPR loci in Streptococcus thermophilus. J Bacteriol. 2008;190(4):1401–12. Epub 2007/12/11. pmid:18065539; PubMed Central PMCID: PMC2238196.
  19. 19. McDermott PF, Tyson GH, Kabera C, Chen Y, Li C, Folster JP, et al. Whole-Genome Sequencing for Detecting Antimicrobial Resistance in Nontyphoidal Salmonella. Antimicrob Agents Chemother. 2016;60(9):5515–20. pmid:27381390; PubMed Central PMCID: PMC4997858.
  20. 20. Pornsukarom S, van Vliet AHM, Thakur S. Whole genome sequencing analysis of multiple Salmonella serovars provides insights into phylogenetic relatedness, antimicrobial resistance, and virulence markers across humans, food animals and agriculture environmental sources. BMC Genomics. 2018;19(1):801. pmid:30400810; PubMed Central PMCID: PMC6218967.
  21. 21. CLSI. Performance standards for antimicrobial susceptibility testing: 25th informational supplement (m100-S25). Clinical and Laboratory Standards Institute. 2015.
  22. 22. McMillan EAG, S. K.; Williams L.; Jove T.; Hiott L. M.; Woodley T. A.; Barrett J. B.; Jackson C. R.; Wasilenko J. L.; Simmons M.; Tillman G. E.; McClelland M.; Frye J. G. Antimicrobial Resistance Genes, Cassettes, and Plasmids Present in Salmonella enterica Associated With United States Food Animals. Frontiers in Microbiology. 2019;10:832. pmid:31057528
  23. 23. Coil D, Jospin G, Darling AE. A5-miseq: an updated pipeline to assemble microbial genomes from Illumina MiSeq data. Bioinformatics. 2015;31(4):587–9. pmid:25338718.
  24. 24. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9. Epub 2014/03/20. pmid:24642063.
  25. 25. Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MT, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31(22):3691–3. Epub 2015/07/23. pmid:26198102; PubMed Central PMCID: PMC4817141.
  26. 26. Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003;4:41. pmid:12969510; PubMed Central PMCID: PMC222959.
  27. 27. Gupta SK, Padmanabhan BR, Diene SM, Lopez-Rojas R, Kempf M, Landraud L, et al. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. Antimicrob Agents Chemother. 2014;58(1):212–20. Epub 2013/10/23. pmid:24145532; PubMed Central PMCID: PMC3910750.
  28. 28. Pal C, Bengtsson-Palme J, Rensing C, Kristiansson E, Larsson DG. BacMet: antibacterial biocide and metal resistance genes database. Nucleic Acids Res. 2014;42(Database issue):D737–43. pmid:24304895; PubMed Central PMCID: PMC3965030.
  29. 29. Chen L, Yang J, Yu J, Yao Z, Sun L, Shen Y, et al. VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 2005;33(Database issue):D325–8. Epub 2004/12/21. pmid:15608208; PubMed Central PMCID: PMC539962.
  30. 30. Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, et al. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44(W1):W16–21. pmid:27141966; PubMed Central PMCID: PMC4987931.
  31. 31. Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35(Web Server issue):W52–7. Epub 2007/06/01. pmid:17537822; PubMed Central PMCID: PMC1933234.
  32. 32. Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27(7):1009–10. pmid:21278367; PubMed Central PMCID: PMC3065679.
  33. 33. Moura A, Soares M, Pereira C, Leitao N, Henriques I, Correia A. INTEGRALL: a database and search engine for integrons, integrases and gene cassettes. Bioinformatics. 2009;25(8):1096–8. pmid:19228805.
  34. 34. Carattoli A, Zankari E, Garcia-Fernandez A, Voldby Larsen M, Lund O, Villa L, et al. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58(7):3895–903. pmid:24777092; PubMed Central PMCID: PMC4068535.
  35. 35. Kisand V, Lettieri T. Genome sequencing of bacteria: sequencing, de novo assembly and rapid analysis using open source tools. BMC Genomics. 2013;14:211. pmid:23547799; PubMed Central PMCID: PMC3618134.
  36. 36. McClelland M, Sanderson KE, Spieth J, Clifton SW, Latreille P, Courtney L, et al. Complete genome sequence of Salmonella enterica serovar Typhimurium LT2. Nature. 2001;413(6858):852–6. pmid:11677609.
  37. 37. Konstantinidis KT, Tiedje JM. Genomic insights that advance the species definition for prokaryotes. Proc Natl Acad Sci U S A. 2005;102(7):2567–72. pmid:15701695; PubMed Central PMCID: PMC549018.
  38. 38. Tyson GH, McDermott PF, Li C, Chen Y, Tadesse DA, Mukherjee S, et al. WGS accurately predicts antimicrobial resistance in Escherichia coli. J Antimicrob Chemother. 2015;70(10):2763–9. pmid:26142410.
  39. 39. Ramirez MS, Tolmasky ME. Aminoglycoside modifying enzymes. Drug Resist Updat. 2010;13(6):151–71. pmid:20833577; PubMed Central PMCID: PMC2992599.
  40. 40. Magnet S, Courvalin P, Lambert T. Activation of the cryptic aac(6')-Iy aminoglycoside resistance gene of Salmonella by a chromosomal deletion generating a transcriptional fusion. J Bacteriol. 1999;181(21):6650–5. pmid:10542165; PubMed Central PMCID: PMC94128.
  41. 41. Carattoli A. Animal reservoirs for extended spectrum beta-lactamase producers. Clin Microbiol Infect. 2008;14 Suppl 1:117–23. pmid:18154535.
  42. 42. Frye JG, Jackson CR. Genetic mechanisms of antimicrobial resistance identified in Salmonella enterica, Escherichia coli, and Enteroccocus spp. isolated from U.S. food animals. Front Microbiol. 2013;4:135. pmid:23734150; PubMed Central PMCID: PMC3661942.
  43. 43. Scholz P, Haring V, Wittmann-Liebold B, Ashman K, Bagdasarian M, Scherzinger E. Complete nucleotide sequence and gene organization of the broad-host-range plasmid RSF1010. Gene. 1989;75(2):271–88. pmid:2653965.
  44. 44. Cambray G, Guerout AM, Mazel D. Integrons. Annu Rev Genet. 2010;44:141–66. pmid:20707672.
  45. 45. Benacer D, Thong KL, Watanabe H, Puthucheary SD. Characterization of drug resistant Salmonella enterica serotype Typhimurium by antibiograms, plasmids, integrons, resistance genes and PFGE. J Microbiol Biotechnol. 2010;20(6):1042–52. pmid:20622506.
  46. 46. Frye JG, Lindsey RL, Meinersmann RJ, Berrang ME, Jackson CR, Englen MD, et al. Related antimicrobial resistance genes detected in different bacterial species co-isolated from swine fecal samples. Foodborne Pathog Dis. 2011;8(6):663–79. pmid:21385089.
  47. 47. Lopes GV, Michael GB, Cardoso M, Schwarz S. Antimicrobial resistance and class 1 integron-associated gene cassettes in Salmonella enterica serovar Typhimurium isolated from pigs at slaughter and abattoir environment. Vet Microbiol. 2016;194:84–92. pmid:27142182.
  48. 48. Roberts MC. Update on acquired tetracycline resistance genes. FEMS Microbiol Lett. 2005;245(2):195–203. pmid:15837373.
  49. 49. Daly M, Villa L, Pezzella C, Fanning S, Carattoli A. Comparison of multidrug resistance gene regions between two geographically unrelated Salmonella serotypes. J Antimicrob Chemother. 2005;55(4):558–61. pmid:15722395.
  50. 50. Peters ED, Leverstein-van Hall MA, Box AT, Verhoef J, Fluit AC. Novel gene cassettes and integrons. Antimicrob Agents Chemother. 2001;45(10):2961–4. pmid:11557503; PubMed Central PMCID: PMC90765.
  51. 51. Abatcha MG, Effarizah ME, Rusul G. Antibiotic susceptibility and molecular characterization of Salmonella enterica serovar Paratyphi B isolated from vegetables and processing environment in Malaysia. Int J Food Microbiol. 2019;290:180–3. pmid:30342248.
  52. 52. Hsu SC, Chiu TH, Pang JC, Hsuan-Yuan CH, Chang GN, Tsen HY. Characterisation of antimicrobial resistance patterns and class 1 integrons among Escherichia coli and Salmonella enterica serovar Choleraesuis strains isolated from humans and swine in Taiwan. Int J Antimicrob Agents. 2006;27(5):383–91. pmid:16621462.
  53. 53. Khemtong S, Chuanchuen R. Class 1 integrons and Salmonella genomic island 1 among Salmonella enterica isolated from poultry and swine. Microb Drug Resist. 2008;14(1):65–70. pmid:18328001.
  54. 54. Kwon HJ, Kim TE, Cho SH, Seol JG, Kim BJ, Hyun JW, et al. Distribution and characterization of class 1 integrons in Salmonella enterica serotype Gallinarum biotype Gallinarum. Vet Microbiol. 2002;89(4):303–9. pmid:12383639.
  55. 55. Adesiji YO, Deekshit VK, Karunasagar I. Antimicrobial-resistant genes associated with Salmonella spp. isolated from human, poultry, and seafood sources. Food Sci Nutr. 2014;2(4):436–42. pmid:25473501; PubMed Central PMCID: PMC4221842.
  56. 56. Almeida F, Seribelli AA, Medeiros MIC, Rodrigues DDP, de MelloVarani A, Luo Y, et al. Phylogenetic and antimicrobial resistance gene analysis of Salmonella Typhimurium strains isolated in Brazil by whole genome sequencing. PLoS One. 2018;13(8):e0201882. pmid:30102733; PubMed Central PMCID: PMC6089434.
  57. 57. El-Sharkawy H, Tahoun A, El-Gohary AEA, El-Abasy M, El-Khayat F, Gillespie T, et al. Epidemiological, molecular characterization and antibiotic resistance of Salmonella enterica serovars isolated from chicken farms in Egypt. Gut Pathog. 2017;9:8. pmid:28203289; PubMed Central PMCID: PMC5301364.
  58. 58. Guerra B, Junker E, Helmuth R. Incidence of the recently described sulfonamide resistance gene sul3 among German Salmonella enterica strains isolated from livestock and food. Antimicrob Agents Chemother. 2004;48(7):2712–5. pmid:15215132; PubMed Central PMCID: PMC434208.
  59. 59. Infante B, Grape M, Larsson M, Kristiansson C, Pallecchi L, Rossolini GM, et al. Acquired sulphonamide resistance genes in faecal Escherichia coli from healthy children in Bolivia and Peru. Int J Antimicrob Agents. 2005;25(4):308–12. pmid:15784310.
  60. 60. Glenn LM, Lindsey RL, Folster JP, Pecic G, Boerlin P, Gilmour MW, et al. Antimicrobial resistance genes in multidrug-resistant Salmonella enterica isolated from animals, retail meats, and humans in the United States and Canada. Microb Drug Resist. 2013;19(3):175–84. pmid:23350745; PubMed Central PMCID: PMC4665089.
  61. 61. Komp Lindgren P, Marcusson LL, Sandvang D, Frimodt-Moller N, Hughes D. Biological cost of single and multiple norfloxacin resistance mutations in Escherichia coli implicated in urinary tract infections. Antimicrob Agents Chemother. 2005;49(6):2343–51. pmid:15917531; PubMed Central PMCID: PMC1140522.
  62. 62. Morgan-Linnell SK, Zechiedrich L. Contributions of the combined effects of topoisomerase mutations toward fluoroquinolone resistance in Escherichia coli. Antimicrob Agents Chemother. 2007;51(11):4205–8. pmid:17682104; PubMed Central PMCID: PMC2151436.
  63. 63. Tamamura Y, Tanaka K, Akiba M, Kanno T, Hatama S, Ishihara R, et al. Complete nucleotide sequences of virulence-resistance plasmids carried by emerging multidrug-resistant Salmonella enterica Serovar Typhimurium isolated from cattle in Hokkaido, Japan. PLoS One. 2013;8(10):e77644. pmid:24155970; PubMed Central PMCID: PMC3796477.
  64. 64. Sundin GW, Bender CL. Dissemination of the strA-strB streptomycin-resistance genes among commensal and pathogenic bacteria from humans, animals, and plants. Mol Ecol. 1996;5(1):133–43. pmid:9147689.
  65. 65. Antunes P, Machado J, Sousa JC, Peixe L. Dissemination of sulfonamide resistance genes (sul1, sul2, and sul3) in Portuguese Salmonella enterica strains and relation with integrons. Antimicrob Agents Chemother. 2005;49(2):836–9. pmid:15673783; PubMed Central PMCID: PMC547296.
  66. 66. Parsons Y, Hall RM, Stokes HW. A new trimethoprim resistance gene, dhfrX, in the In7 integron of plasmid pDGO100. Antimicrob Agents Chemother. 1991;35(11):2436–9. pmid:1804022; PubMed Central PMCID: PMC245401.
  67. 67. Steiniger-White M, Rayment I, Reznikoff WS. Structure/function insights into Tn5 transposition. Curr Opin Struct Biol. 2004;14(1):50–7. pmid:15102449.
  68. 68. Wright GD, Thompson PR. Aminoglycoside phosphotransferases: proteins, structure, and mechanism. Front Biosci. 1999;4:D9–21. pmid:9872733.
  69. 69. Krauland MG, Marsh JW, Paterson DL, Harrison LH. Integron-mediated multidrug resistance in a global collection of nontyphoidal Salmonella enterica isolates. Emerg Infect Dis. 2009;15(3):388–96. pmid:19239750; PubMed Central PMCID: PMC2666292.
  70. 70. Boyd D, Cloeckaert A, Chaslus-Dancla E, Mulvey MR. Characterization of variant Salmonella genomic island 1 multidrug resistance regions from serovars Typhimurium DT104 and Agona. Antimicrob Agents Chemother. 2002;46(6):1714–22. pmid:12019080; PubMed Central PMCID: PMC127246.
  71. 71. Doublet B, Lailler R, Meunier D, Brisabois A, Boyd D, Mulvey MR, et al. Variant Salmonella genomic island 1 antibiotic resistance gene cluster in Salmonella enterica serovar Albany. Emerg Infect Dis. 2003;9(5):585–91. pmid:12737743; PubMed Central PMCID: PMC2972765.
  72. 72. Doublet B, Weill FX, Fabre L, Chaslus-Dancla E, Cloeckaert A. Variant Salmonella genomic island 1 antibiotic resistance gene cluster containing a novel 3'-N-aminoglycoside acetyltransferase gene cassette, aac(3)-Id, in Salmonella enterica serovar newport. Antimicrob Agents Chemother. 2004;48(10):3806–12. pmid:15388438; PubMed Central PMCID: PMC521890.
  73. 73. Ebner P, Garner K, Mathew A. Class 1 integrons in various Salmonella enterica serovars isolated from animals and identification of genomic island SGI1 in Salmonella enterica var. Meleagridis. J Antimicrob Chemother. 2004;53(6):1004–9. pmid:15117931.
  74. 74. Meunier D, Boyd D, Mulvey MR, Baucheron S, Mammina C, Nastasi A, et al. Salmonella enterica serotype Typhimurium DT 104 antibiotic resistance genomic island I in serotype paratyphi B. Emerg Infect Dis. 2002;8(4):430–3. pmid:11971780; PubMed Central PMCID: PMC2730239.
  75. 75. Doublet B, Boyd D, Mulvey MR, Cloeckaert A. The Salmonella genomic island 1 is an integrative mobilizable element. Mol Microbiol. 2005;55(6):1911–24. pmid:15752209.
  76. 76. Jacobsen A, Hendriksen RS, Aaresturp FM, Ussery DW, Friis C. The Salmonella enterica pan-genome. Microb Ecol. 2011;62(3):487–504. pmid:21643699; PubMed Central PMCID: PMC3175032.
  77. 77. Fu S, Hiley L, Octavia S, Tanaka MM, Sintchenko V, Lan R. Comparative genomics of Australian and international isolates of Salmonella Typhimurium: correlation of core genome evolution with CRISPR and prophage profiles. Sci Rep. 2017;7(1):9733. pmid:28851865; PubMed Central PMCID: PMC5575072.
  78. 78. Dilucca M, Cimini G, Giansanti A. Essentiality, conservation, evolutionary pressure and codon bias in bacterial genomes. Gene. 2018;663:178–88. pmid:29678658.
  79. 79. Yu H, Wang J, Ye J, Tang P, Chu C, Hu S, et al. Complete nucleotide sequence of pSCV50, the virulence plasmid of Salmonella enterica serovar Choleraesuis SC-B67. Plasmid. 2006;55(2):145–51. pmid:16257053.
  80. 80. Kingsley RA, Humphries AD, Weening EH, De Zoete MR, Winter S, Papaconstantinopoulou A, et al. Molecular and phenotypic analysis of the CS54 island of Salmonella enterica serotype typhimurium: identification of intestinal colonization and persistence determinants. Infect Immun. 2003;71(2):629–40. pmid:12540539; PubMed Central PMCID: PMC145368.
  81. 81. Alsop J. An outbreak of salmonellosis in a swine finishing barn. J Swine Health Prod. 2005;13:4.
  82. 82. Pilar AV, Reid-Yu SA, Cooper CA, Mulder DT, Coombes BK. GogB is an anti-inflammatory effector that limits tissue damage during Salmonella infection through interaction with human FBXO22 and Skp1. PLoS Pathog. 2012;8(6):e1002773. pmid:22761574; PubMed Central PMCID: PMC3386239.
  83. 83. Holt KE, Thomson NR, Wain J, Langridge GC, Hasan R, Bhutta ZA, et al. Pseudogene accumulation in the evolutionary histories of Salmonella enterica serovars Paratyphi A and Typhi. BMC Genomics. 2009;10:36. pmid:19159446; PubMed Central PMCID: PMC2658671.
  84. 84. Ledeboer NA, Frye JG, McClelland M, Jones BD. Salmonella enterica serovar Typhimurium requires the Lpf, Pef, and Tafi fimbriae for biofilm formation on HEp-2 tissue culture cells and chicken intestinal epithelium. Infect Immun. 2006;74(6):3156–69. pmid:16714543; PubMed Central PMCID: PMC1479237.
  85. 85. Bender JK, Wille T, Blank K, Lange A, Gerlach RG. LPS structure and PhoQ activity are important for Salmonella Typhimurium virulence in the Galleria mellonella infection model [corrected]. PLoS One. 2013;8(8):e73287. pmid:23951347; PubMed Central PMCID: PMC3738532.
  86. 86. Murata T, Tseng W, Guina T, Miller SI, Nikaido H. PhoPQ-mediated regulation produces a more robust permeability barrier in the outer membrane of Salmonella enterica serovar typhimurium. J Bacteriol. 2007;189(20):7213–22. pmid:17693506; PubMed Central PMCID: PMC2168427.
  87. 87. Boyer F, Fichant G, Berthod J, Vandenbrouck Y, Attree I. Dissecting the bacterial type VI secretion system by a genome wide in silico analysis: what can be learned from available microbial genomic resources? BMC Genomics. 2009;10:104. pmid:19284603; PubMed Central PMCID: PMC2660368.
  88. 88. Pukatzki S, McAuley SB, Miyata ST. The type VI secretion system: translocation of effectors and effector-domains. Curr Opin Microbiol. 2009;12(1):11–7. pmid:19162533.
  89. 89. Blondel CJ, Jimenez JC, Contreras I, Santiviago CA. Comparative genomic analysis uncovers 3 novel loci encoding type six secretion systems differentially distributed in Salmonella serotypes. BMC Genomics. 2009;10:354. pmid:19653904; PubMed Central PMCID: PMC2907695.
  90. 90. Zeng L, Zhang L, Wang P, Meng G. Structural basis of host recognition and biofilm formation by Salmonella Saf pili. Elife. 2017;6. pmid:29125121; PubMed Central PMCID: PMC5700814.
  91. 91. Li YF, Poole S, Rasulova F, McVeigh AL, Savarino SJ, Xia D. A receptor-binding site as revealed by the crystal structure of CfaE, the colonization factor antigen I fimbrial adhesin of enterotoxigenic Escherichia coli. J Biol Chem. 2007;282(33):23970–80. pmid:17569668.
  92. 92. Lang AS, Zhaxybayeva O, Beatty JT. Gene transfer agents: phage-like elements of genetic exchange. Nat Rev Microbiol. 2012;10(7):472–82. pmid:22683880; PubMed Central PMCID: PMC3626599.
  93. 93. Penades JR, Chen J, Quiles-Puchalt N, Carpena N, Novick RP. Bacteriophage-mediated spread of bacterial virulence genes. Curr Opin Microbiol. 2015;23:171–8. pmid:25528295.
  94. 94. Ehrbar K, Hardt WD. Bacteriophage-encoded type III effectors in Salmonella enterica subspecies 1 serovar Typhimurium. Infect Genet Evol. 2005;5(1):1–9. pmid:15567133.
  95. 95. Schmieger H, Schicklmaier P. Transduction of multiple drug resistance of Salmonella enterica serovar typhimurium DT104. FEMS Microbiol Lett. 1999;170(1):251–6. pmid:9919675.
  96. 96. Figueroa-Bossi N, Bossi L. Inducible prophages contribute to Salmonella virulence in mice. Mol Microbiol. 1999;33(1):167–76. pmid:10411733.
  97. 97. Stanley TL, Ellermeier CD, Slauch JM. Tissue-specific gene expression identifies a gene in the lysogenic phage Gifsy-1 that affects Salmonella enterica serovar typhimurium survival in Peyer's patches. J Bacteriol. 2000;182(16):4406–13. pmid:10913072; PubMed Central PMCID: PMC94610.
  98. 98. Villafane R, Zayas M, Gilcrease EB, Kropinski AM, Casjens SR. Genomic analysis of bacteriophage epsilon 34 of Salmonella enterica serovar Anatum (15+). BMC Microbiol. 2008;8:227. pmid:19091116; PubMed Central PMCID: PMC2629481.
  99. 99. Beloglazova N, Petit P, Flick R, Brown G, Savchenko A, Yakunin AF. Structure and activity of the Cas3 HD nuclease MJ0384, an effector enzyme of the CRISPR interference. EMBO J. 2011;30(22):4616–27. pmid:22009198; PubMed Central PMCID: PMC3243599.
  100. 100. Touchon M, Rocha EP. The small, slow and specialized CRISPR and anti-CRISPR of Escherichia and Salmonella. PLoS One. 2010;5(6):e11126. pmid:20559554; PubMed Central PMCID: PMC2886076.