Genomic characterization of Salmonella Cerro ST367, an emerging Salmonella subtype in cattle in the United States
BMC Genomics volume 15, Article number: 427 (2014)
Within the last decade, Salmonella enterica subsp. enterica serovar Cerro (S. Cerro) has become one of the most common serovars isolated from cattle and dairy farm environments in the northeastern US. The fact that this serovar is commonly isolated from subclinically infected cattle and is rarely associated with human disease, despite its frequent isolation from cattle, has led to the hypothesis that this emerging serovar may be characterized by reduced virulence. We applied comparative and population genomic approaches to (i) characterize the evolution of this recently emerged serovar and to (ii) gain a better understanding of genomic features that could explain some of the unique epidemiological features associated with this serovar.
In addition to generating a de novo draft genome for one Salmonella Cerro strain, we also generated whole genome sequence data for 26 additional S. Cerro isolates, including 16 from cattle operations in New York (NY) state, 2 from human clinical cases from NY in 2008, and 8 from diverse animal sources (7 from Washington state and 1 from Florida). All isolates sequenced in this study represent sequence type ST367. Population genomic analysis showed that isolates from the NY cattle operations form a well-supported clade within S. Cerro ST367 (designated here “NY bovine clade”), distinct from isolates from Washington state, Florida and the human clinical cases. A molecular clock analysis indicates that the most recent common ancestor of the NY bovine clade dates back to 1998, supporting the recent emergence of this clone.
Comparative genomic analyses revealed several relevant genomic features of S. Cerro ST367, that may be responsible for reduced virulence of S. Cerro, including an insertion creating a premature stop codon in sopA. In addition, patterns of gene deletion in S. Cerro ST367 further support adaptation of this clone to a unique ecological or host related niche.
Our results indicate that the increase in prevalence of S. Cerro ST367 is caused by a highly clonal subpopulation and that S. Cerro ST367 is characterized by unique genomic deletions that may indicate adaptation to specific ecological niches and possibly reduced virulence in some hosts.
Genomic characteristics associated with the emergence or reemergence of pathogens in livestock operations can be subdivided into two categories; (i) genomic features that increase the adaptation to a host, or facilitate the jump to a new host species, or (ii) genomic features that provide increased adaptation to environmental factors in the livestock environment, such as antibiotic resistance. Comparative and population genomic studies are particularly suited to determine which features are responsible for the emergence of certain pathogens. For instance, Price et al.  showed that a putative host jump, from humans to livestock, in a clonal complex in Staphylococcus aureus was associated with the loss of phage-carried human virulence genes and with the acquisition of tetracycline and methicillin resistance.
Salmonella enterica is one of the most frequent causes of bacterial foodborne illness and death in the United States . In Salmonella, examples of emergent clones include S. Typhimurium DT 104, a multidrug resistant clone, which has seen a global epidemic spread from 1990 , and S. enterica serovar 4,5,12:i:–, a monophasic variant of S. Typhimurium, which showed a global increase in the mid-1990s . In this study, we present comparative and population genomic research on S. enterica subsp. enterica serovar Cerro (S. Cerro). S. Cerro is rarely associated with human disease, with only one outbreak reported in the US so far that could be solely attributed to this serovar ; an additional outbreak was recently reported and it was linked to multiple serovars, including S. Cerro . However, this Salmonella serovar has emerged over the last decade as one of the most abundant Salmonella serovars in cattle operations in the northeastern US , including one of the most common serovars among subclinical dairy cattle and in the dairy farm environment  in the northeastern United States. Most of the S. Cerro isolated from cattle and farms represent one pulsed field electrophoresis (PFGE) type, indicating that a single clonal lineage is involved in this emergence . It is unknown what causes S. Cerro to be associated with cattle and why it is rarely involved in human disease. Therefore, we hypothesize that S. Cerro has distinct genomic characteristics that explain its association with cattle and limited association with human disease.
Results and discussion
De novo assembly shows that S.Cerro FSL R8-0235 has a genome size of approximately 4.7 Mbp, contains six prophage regions and represents MLST sequence type ST367
After exclusion of contigs fewer than 200 bp, the total length of the S. Cerro FSL R8-0235 draft de novo assembly was 4,675,817 bp. The assembly consisted of 126 contigs, with a contig N50 of 292,947 bp, and a maximum contig length of 691,181 bp. The average coverage depth of the assembly was 96X. One contig, contig 016, contained genes of an IncI1-like plasmid, however it is unclear whether this is an integrated or extrachromosomal plasmid. In addition to genes involved in plasmid transfer, stability and replication, this plasmid also carries genes encoding a resistance nodulation division (RND) efflux pump . However, none of the isolates sequenced in this study showed resistance to single or multiple antimicrobial agents. No evidence for the existence of additional plasmids within the genome was found. This may be at least partially due to the presence of a DNA phosphorothioation-dependent restriction modification (RM) system in all S. Cerro strains examined in this study. While this RM system has been well characterized in S. Cerro , a PSI-BLAST search reveals this type of RM system is very rare among Salmonella, and only found in a limited number of sequenced Salmonella strains of serovars Saintpaul (SARA23, str. 9712, str. JO2008), Namur (str.05-2929) and Panama (ATCC 7378).
Prediction of lysogenic prophages and prophage remnants in the S. Cerro FSL R8-0235 genome was performed using Prophinder . Six putative prophages or remnants of prophages, ranging in length from 5.78 to 31.52 Kb, were predicted to be present in the S. Cerro FSL R8-0235 genome (Table 1). The six prophage regions, which we refer to as prophage 57014, 57017, 57018, 57023, 57024, and 57025, were compared, using RAST , to previously sequenced genomes to identify homologous regions. Prophages 57023 and 57025 (Table 1) are similar in composition to a S. Typhimurium ES18-like bacteriophage, while 57014 shows similarity to an Enterobacteria P22-like prophage. While typical Salmonella-associated prophages, such as Gifsy-1, Gifsy-2, Fels-1, and Fels-2 [13–15] were not predicted to be present in the S. Cerro FSL R8-0235 genome, prophage 57024 shared many genes with a prophage found in Photorhabdus luminescens subsp. laumondii TT01, which has been described as a successful insect pathogen as well as symbiont of soil entomopathogenic nematodes .
Genome assembly based multi locus sequence typing (MLST) was performed using the online tool  of the Center for Genomic Epidemiology (Lyngby, Denmark; http://www.genomicepidemiology.org/) and an additional BLASTN search. This analysis revealed that S. Cerro FSL R8-0235 belongs to sequence type (ST) 367. According to the Salmonella MLST database (http://mlst.warwick.ac.uk/mlst/dbs/Senterica) ST367 is associated with a S. Cerro isolate from a human case in Germany in 1985. The database also contains an accession of the type strain of S. Cerro, isolated from swine in 1936 in Uruguay. This strain belongs to ST1291 and displays a different allelic type at each of the seven MLST loci. S. Cerro therefore is very likely to be polyphyletic, which makes interpretation of historical references without genomic or MLST sequence data difficult. Because all isolates sequenced in this study belong to ST367, we will refer to these isolates as S. Cerro ST367 from here on. Timme et al.  recently published sequence data for another S. Cerro ST367 strain (strain 818; NZ_AOZJ00000000); this group showed that, among all serovars that have been sequenced so far, S. Adelaide FSL A4-669 is most closely related to S. Cerro ST367 which is consistent with our study (see below).
Population genomic analysis of 27 Salmonella serovar Cerro isolates suggests a recent clonal expansion of a bovine-associated S.Cerro lineage
To infer whether the S. Cerro isolates associated with bovine hosts and cattle-associated environments form separate subpopulations from S. Cerro isolated from other sources, we obtained whole genome sequencing data for 26 additional isolates (Table 2). After removal of putative recombinogenic regions, as identified by BratNextGen , and SNPs that were present in fewer than 90% of the isolates, 343 SNPs were left for analysis. To assess the presence of a temporal signal in the dataset, a Path-O-Gen (available from http://tree.bio.ed.ac.uk/software/pathogen/) analysis was performed using a maximum likelihood tree inferred from the SNP data set. This analysis showed a correlation (Pearson’s Correlation Coefficient 0.80, R2 = 0.645) between the time of isolation of the individual isolates and the root-to-tip divergence, indicating a temporal signal for this dataset and justifying a molecular clock based phylogenetic analysis. A Bayesian analysis, assuming a relaxed molecular clock and a constant population size, inferred the mean mutation rate for the core genome of the 27 S. Cerro isolates to be 2.4 × 10−7/site/year (95% Highest Probability Density (HPD) 1.5 × 10−7 – 3.3 × 10−7). This mutation rate is comparable to mutation rates estimated for Buchnera aphidicola and Helicobacter pylori, but about twice as fast as recently inferred for S. Agona . The New York bovine isolates are found in a well-supported (posterior probability 1.0) clade (NY bovine clade; see Figure 1), well separated from the isolates from Washington state, Florida, and the human clinical isolates from New York state. This may indicate that, although isolates of S. Cerro of the bovine-associated clade were prevalent in farm environments, and thus farm personnel would be frequently exposed to this clone, this clone was not responsible for the human cases in New York state represented by these two isolates. The time of the emergence of the most recent ancestor (MRCA) of the NY bovine clade is estimated to be 1998 (95% HPD 1991–2003). The NY bovine clade is further split up into two clades: (i) a clade with two isolates from northeastern New York (Figure 1: clade 1) and (ii) a clade with 15 bovine associated isolates from western NY state (Figure 1: clade 2). The MRCA of the latter clade dates back to 2002 (95% HPD 1999–2005). Within clade 2, two well supported clusters were identified (marked ‘a’ and ‘b’ in Figure 1). Specifically, ‘cluster a’ contains six isolates that were isolated from Steuben county (NY) and the neighboring Livingston county (NY). This finding suggests a phylogeographic signal in the dataset, which should facilitate more detailed tracing of the emergence of S. Cerro ST367 throughout the northeastern US with a larger sampling and a population genomic analysis.
Genome sequence analysis reveals a stepwise evolution, of S. Cerro ST367 to a bovine-associated clade, characterized by deletion of selected operons and acquisition of a premature stop codon in sopA
Loss or gain of genes within bacterial populations may indicate niche adaption of bacterial subpopulations . To infer patterns of gene loss, we mapped reads of the 27 S. Cerro isolates against well-annotated genomes such as those of S. Typhimurium LT2, S. Typhi CT18, and S. Choleraesuis SC-B67. In addition we mapped the reads of S. Adelaide FSL A4-669  against these genomes, to determine if the patterns of absence were also observed in the most recent common ancestor of this serovar and the S. Cerro population studied here. Reads of the 27 S. Cerro isolates mapped to 86, 88, and 90% of the coding sequences in S. Typhi CT18, S. Typhimurium LT2, and S. Choleraesuis SC-B67, respectively. This is very similar to the percentage of genes shared (89%) between S. Typhimurium LT2 and S. Typhi CT18  and falls in the higher end of the range observed by Jacobsen et al.  for a wide variety of Salmonella serovars. The genome size, and the high number of shared genes thus suggest that the lineage of S. Cerro studied here did not experience notable genome reduction.
Mapping of sequence reads of the isolates of the S. Cerro population further revealed a pattern of gene absence generally conserved within the S. Cerro population sampled here, suggesting that most of the genomic characteristics associated with the emergence of S. Cerro among bovine-associated habitats were present in the MRCA of this S. Cerro clade. Interestingly, loss of some SPIs (Salmonella Pathogenicity Islands) that were found here to be absent or partially absent (gene deletions) from the Cerro population studied (i.e., ST367), but are present in S. Typhimurium LT2 or S. Typhi CT18, have been associated with attenuation of virulence. Specifically, the genomic island at S. Typhi SPI-10 locus is completely absent from the S. Cerro ST367 isolates examined here; this SPI has been associated with virulence in mice . Chaudhuri et al.  also showed that significant reduction of fitness of S. Typhimurium SL1344 is observed during intestinal colonization of cattle when genes in SPI-10 (in particular STM4489) are disrupted by transposon insertion. Genes homologous to (i) STM2230.1c to STM2240 of SPI-12, and (ii) STM3117, STM3123, and STM3119 to STM3121 of SPI-13 were also found to be absent from S. Cerro ST367; these SPIs have been associated with systemic infection of mice in S. Typhimurium , and replication in macrophages (SPI-13: ). Furthermore, disruption of STM2231 in SPI-12 and STM3123 in SPI-13 was previously shown to cause significant reduction in fitness in S. Typhimurium SL1344 during intestinal colonization of cattle . In addition, homologs of STM0293, STM0294 and STM0299 are deleted in S. Cerro ST367. These genes are found in SPI-16, a SPI associated with intestinal persistence in mice . Disruption of STM0293 in S. Typhimurium has been shown to cause reduced fitness with regard to intestinal colonization of cattle . Most of the SPI-related genes found to be absent in S. Cerro ST367 were confirmed to be present in S. Adelaide FSL A4-669, suggesting loss of these genes/SPIs occurred after the divergence of S. Adelaide from the most recent common ancestor of S. Cerro ST367. We found evidence for the presence of four complete toxin-antitoxin (TA) modules (STM 2954.1 N-2955.S; STM4030.S-4031; STM3777-78 and STM4449-50) within the S. Cerro genomes studied here. This is interesting as De la Cruz et al.  suggested that TA modules in Salmonella play a role in virulence, and that the number of genomically encoded TA modules is correlated with pathogenicity of individual strains. By comparison, the number of TA modules in S. enterica subsp. enterica ranges from 5 (S. Paratyphi B SPB7) to 10 (S. Typhimurium LT2), making S. Cerro ST367 one of the subsp. enterica serotypes with the lowest number of TA modules. The number of TA modules in S. Cerro ST367 is similar to that observed in Salmonella enterica subsp. arizonae, a subspecies which is predominantly found in cold blooded hosts and does generally not seem to cause illness in warm blooded hosts . Complete or partial absence of some SPIs in all S. Cerro ST367 and the low number of TA modules in the genome, thus suggests a putative shift of S. Cerro in host and/or tissue tropism before the emergence of the NY bovine-associated clade.
The hypothesis that the S. Cerro population studied here shows unique host and/or tissue tropism characteristics is also supported by the finding that all 27 S. Cerro ST367 isolates sequenced here were found to carry a premature stop codon in sopA, causing a truncation of the gene from 782 aa (in S. Typhimurium LT2) to 433 aa. Previous studies have shown that SopA is involved in virulence during bovine gastrointestinal infections by S. Typhimurium and S. Dublin [34, 35], and that sopA mutations are implicated in reduced polymorphonuclear (PMN) cell migration [34, 36], and fluid secretion in ileal loops in calves . Premature stop codons in sopA have been found in S. Typhi, S. Paratyphi A, and S. Gallinarum and it has been suggested that loss of a functional SopA has been an important factor in the virulence and adaptation of these serovars to a systemic niche in certain hosts [37, 38].
Interestingly, the one base-pair insertion responsible for the premature stop codon occurs within a ~10 bp region of sopA that also contains deletions in S. Typhi and S. Paratyphi A (Figure 2). While S. Typhi and S. Paratyphi A contain additional mutations that may have caused loss of function of SopA , the occurrence of deletions in the same region in S. Cerro sopA suggests this is a replication error prone region in the genome. A conserved domain search (http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) against the Conserved Domain Database  of the aa sequence of the truncated SopA in S. Cerro ST367 revealed the premature stop is situated in the SopA central domain  of the gene. Furthermore, the truncated SopA protein lacks the capsase-3 cleavage sites, which have been demonstrated to be important in induction of PMN transepithelial migration in S. Typhimurium . Although specifically the disruption of the main functional domain in SopA and the loss of the capsase-3 cleavage sites suggest loss of function of SopA in S. Cerro ST367, further molecular genetic experiments have to be conducted to reveal if truncation of SopA in S. Cerro ST367 has lead to loss of function of this gene, and how it affects host cell invasion (as suggested by Raffatelu et al. ) and other SopA associated aspects of Salmonella virulence.
Read mapping also showed one gene cluster to be stepwise deleted in the NY bovine clade (Figure 1). This gene cluster contains homologs of the S. Typhimurium LT2 genes STM1633 to STM1637. This gene cluster encodes a D-alanine transporter and has been recently shown to be required for intracellular survival in murine macrophage-like cells , and disruption of STM1637 has been shown to cause a significant reduction in fitness in intestinal colonization in cattle in S. Typhimurium . This gene cluster is present in all 10 Cerro ST367 isolates that do not belong to the bovine clade. Two isolates (FSL R8-2008, FSL R8-2639) lack two genes (STM1633, STM1634) in this gene cluster. These two isolates represent a clade that split off early from the remaining NY bovine-associated-population. The remaining 15 isolates in this clade lack the entire gene cluster (Figure 1). The (partial) absence of the D-alanine transporter gene cluster is currently the only genomic feature that differentiates the NY bovine clade from the remaining population (including isolates from the NY human cases).
S. Cerro displays reduced invasiveness of human epithelial cells compared to other Salmonellaserovars commonly isolated from bovine sources
The comparative genomic analyses described above suggest S. Cerro lacks several functional genes and genomic elements that are involved in invasion and intracellular survival. To assess if strains of S. Cerro ST367 population (Table 2) studied here are impaired in their ability to invade human intestinal epithelial cells, Caco-2 cells were infected with S. Typhimurium (n = 4), S. Newport (n = 4), S. Kentucky (n = 4), and S. Cerro (n = 4). Each serovar was represented by one isolate each from a bovine clinical case, a subclinically infected bovine host, an environmental sample and a human clinical case. S. Cerro isolates were significantly less invasive than isolates of serovars Typhimurium (P < 0.0001) and Newport (P < 0.0001), but not significantly different from S. Kentucky (P = 0.0734) (Figure 3). However, the overall invasiveness of S. Kentucky seems to be skewed by the presence of one isolate from a human clinical case, which shows very low invasion. When this outlier is excluded from the analysis, the S. Cerro isolates are also significantly less invasive than S. Kentucky (p = 0.004). Thus, consistent with our genomic analyses, S. Cerro ST367 seems to be less invasive in human intestinal epithelial cells than the serovars examined here. Future studies on the ability of S. Cerro to invade bovine intestinal epithelial cells and to cause illness in cattle will be necessary though to determine whether S. Cerro or specific subtypes within S. Cerro truly show attenuated bovine virulence.
Comparative genomic analyses of 27 Salmonella Cerro isolates indicate that this serovar lacks several genes that have previously been shown to be involved in the ability of Salmonella serovars to cause intestinal infection. Reduced invasion of human intestinal epithelial cells, as compared to other serovars, further supports the reduced ability of this serovar to cause intestinal infection, however, further experiments are necessary to determine which genes are responsible for this phenotype. Altogether, these results suggest that the emergence of S. Cerro ST367 among livestock operations in the northeastern United States may not be due to increased adaptation to the bovine host, nor to increased antibiotic resistance. Instead, the frequent isolation of this serovar on cattle farms  may reflect that this serovar was able to disperse rapidly as no efforts were undertaken to control its spread (possibly due to a lack of clinical signs, which left infections undetected). Alternatively, or in addition, S. Cerro (or some subtypes within S. Cerro) may have unique phenotypic characteristics that were not discovered through the comparative genomic analyses conducted here, but that facilitate environmental survival or dispersal.
The 27 S. Cerro isolates for genome sequencing (n = 1) and re-sequencing (n = 26) were isolated from 1986 to 2008 from human cases and domesticated and wild animals in 3 different states (i.e., New York, Washington, and Florida; Table 2).
Genome sequencing, assembly and annotation
The genome of S. Cerro FSL R8-0235 was sequenced using the SOLiD™ system (Applied Biosystems, Foster City). Mate-paired 50 bp reads were obtained and a de novo assembly was performed as detailed in Den Bakker et al. . Contigs longer than 200 bp were submitted to the NCBI Prokaryotic Genomes Automatic Annotation Pipeline (PGAAP)  for automated annotation. Unpaired 50 bp reads for the genomes of the additional 26 S. Cerro ST367 isolates were obtained using the SOLiD™ system (Applied Biosystems, Foster City) as detailed in Den Bakker et al. .
SOLiD™ read mapping, population genetics analysis, and read mapping based gene presence/absence analysis
SOLiD™ reads were mapped against a reference genome (FSL R8_0235) using PerM . ComB  was used to for SNP calling and creation of consensus sequences. Regions with coverage less than 10X were masked in the consensus sequences. Consensus sequences created with ComB were used as input for the BratNextGen  recombination detection software, using 100 replicates of 50 iterations each. SNPs in regions that were predicted to be involved in a recombination event with P < 0.01 were excluded from the analysis.
A maximum likelihood (ML) tree based on the SNP data was created in MEGA 5 , and this ML tree was used to test for the presence of a temporal signal in the dataset using Path-O-Gen 1.4 (available from http://tree.bio.ed.ac.uk/software/pathogen/). BEAST version 1.7.5  was used to create a tip-dated phylogeny of the S. Cerro isolates. Four different models differing in assumptions on mutation rate and effective population size (strict clock, constant population size; strict clock, Gaussian Markov random field (GMRF) model ; relaxed clock, constant population size; relaxed clock, GMRF model) were run for 10 million generations each and compared using the Bayes factor as implemented in Tracer version 1.5 (A. Rambaut available from http://tree.bio.ed.ac.uk/software/tracer/).
Read mapping based gene presence/absence analysis was performed by mapping SOLiD™ reads to selected reference genomes using PerM . Coverage per annotated gene feature in the reference genome was subsequently obtained using the ‘coverage’ tool from the BEDtools suite .
Caco-2 cell invasion assays of S. Cerro, S. Kentucky, S. Typhimurium, and S.Newport
To compare the ability of S. Cerro isolates to invade human intestinal epithelial cells, Caco-2 cells were infected with S. Typhimurium (n = 4), S. Newport (n = 4), S. Kentucky (n = 4), and S. Cerro (n = 4), see Additional file 1. Salmonella Typhimurium ATCC® 14028 was used as a positive control and its sirA isogenic mutant as a negative control. All isolates were susceptible to gentamicin as determined by antimicrobial susceptibility testing (MIC values between 0.25 and 1 μg/ml) by the Cornell University Animal Health Diagnostic Center. Salmonella isolates were grown on Luria Bertani (LB) plates at 37°C for 16 hours. A colony was transferred into 5 mL LB broth and incubated 18 hours at 37°C, without shaking. After 18 hours of incubation, 1 mL of each culture was pelleted by centrifugation and re-suspended in 1 mL of Phosphate Buffered Saline (PBS) pH 7.4. Bacterial cells were diluted and Caco-2 cells were inoculated at an MOI of 10. Each strain was inoculated in triplicate in each of the 3 experiments conducted. Appropriate dilutions were plated on LB for calculation of the initial inoculum.
For all the experiments Caco-2 cells were maintained in Dulbecco’s Modified Eagle Medium (DMEM) 20% FBS 1% non-essential amino acids at 37°C and 5.0% CO2, for no more than 50 passages. The 24-well plates were seeded at a concentration of 5.0 × 104 cells/well and incubated at 37°C and 5% CO2 for 48 hours. Thirty minutes before the cells were inoculated with Salmonella, media in the 24-well plate was replaced with fresh media. Caco-2 cells were inoculated, and incubated at 37°C and 5% CO2 for 1 hour, followed by 3 washes with pre-warmed PBS. Fresh media was distributed into each well followed by a 15 minute incubation at 37°C and 5% CO2. Finally, media with gentamicin (50 μg/mL) was added and the cells were incubated for 1 hour at 37°C and 5% CO2. The cells were then lysed by vigorously pipetting 500 μL of chilled water in each well. The bacterial suspensions recovered were plated on LB and incubated at 37°C overnight. Invasion efficiency was calculated as [CFU recovered/CFU infected] × 100. Statistical analysis was performed using SAS software (SAS Institute Inc., Cary, NC, USA). The invasion efficiencies were analyzed using one-way analysis of variance (ANOVA), Tukey post hoc test, and the data was log-transformed to satisfy ANOVA assumptions of normality.
Availability of supporting data section
All raw read files have been deposited in the Sequence Read Archive of the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/Traces/sra/) under entries PRJNA185435, PRJNA187190–187196, PRJNA187371, PRJNA187373, PRJNA187542, PRJNA187545, PRJNA187919–187921, PRJNA187962, PRJNA187963, PRJNA187965–187974, and PRJNA73959. The de novo assembly of strain FSL R8-0235 has been deposited as a Whole Genome Shotgun project at DDBJ/EMBL/GenBank under the accession JMIJ00000000. The version described in this paper is version JMIJ01000000.
Price LB, Stegger M, Hasman H, Aziz M, Larsen J, Andersen PS, Pearson T, Waters AE, Foster JT, Schupp J, Gillece J, Driebe E, Liu CM, Springer B, Zdovc I, Battisti A, Franco A, Zmudzki J, Schwarz S, Butaye P, Jouy E, Pomba C, Porrero MC, Ruimy R, Smith TC, Robinson DA, Weese JS, Arriola CS, Yu F, Laurent F, et al: Staphylococcus aureus CC398: host adaptation and emergence of methicillin resistance in livestock. mBio. 2012, 3: e00305-11.
Scallan E, Hoekstra RM, Angulo FJ, Tauxe RV, Widdowson MA, Roy SL, Jones JL, Griffin PM: Foodborne illness acquired in the United States-major pathogens. Emerging Infect Dis. 2011, 17: 7-15. 10.3201/eid1701.P11101.
Threlfall EJ: Epidemic Salmonella typhimurium DT 104–a truly international multiresistant clone. J Antimicrob Chemother. 2000, 46: 7-10. 10.1093/jac/46.1.7.
Moreno Switt AI, Soyer Y, Warnick LD, Wiedmann M: Emergence, distribution, and molecular and phenotypic characteristics of Salmonella enterica serotype 4,5,12:i. Foodborne Pathog Dis. 2009, 6: 407-415. 10.1089/fpd.2008.0213.
CDC: Epidemiologic notes and reports Salmonellosis associated with Carne Seca - New Mexico. MMWR Morb Mortal Wkly Rep. 2013, 34: 645-646.
Gicquelais RE, Morris JF, Matthews S, Gladden L, Safi H, Grayson C, Slayton RB, Newton AE, Bordonaro R, Wheeler JG, Smith N, Bosch SA, Haselow DT, Arkansas Department of Health: Multiple-serotype Salmonella outbreaks in two state prisons - Arkansas, August 2012. MMWR Morb Mortal Wkly Rep. 2014, 63: 169-173.
Hoelzer K, Cummings KJ, Wright EM, Rodriguez-Rivera LD, Roof SE, Moreno Switt AI, Dumas N, Root T, Schoonmaker-Bopp DJ, Grohn YT, Siler JD, Warnick LD, Hancock DD, Davis MA, Wiedmann M: Salmonella Cerro isolated over the past twenty years from various sources in the US represent a single predominant pulsed-field gel electrophoresis type. Vet Microbiol. 2011, 150: 389-393. 10.1016/j.vetmic.2011.01.026.
Cummings KJ, Warnick LD, Elton M, Rodriguez-Rivera LD, Siler JD, Wright EM, Grohn YT, Wiedmann M: Salmonella enterica serotype Cerro among dairy cattle in New York: an emerging pathogen?. Foodborne Pathog Dis. 2010, 7: 659-665. 10.1089/fpd.2009.0462.
Piddock LJV: Multidrug-resistance efflux pumps - not just for resistance. Nat Rev Microbiol. 2006, 4: 629-636. 10.1038/nrmicro1464.
Xu T, Yao F, Zhou X, Deng Z, You D: A novel host-specific restriction system associated with DNA backbone S-modification in Salmonella. Nucleic Acids Res. 2010, 38: 7133-7141. 10.1093/nar/gkq610.
Lima-Mendez G, van Helden J, Toussaint A, Leplae R: Prophinder: a computational tool for prophage prediction in prokaryotic genomes. Bioinformatics. 2008, 24: 863-865. 10.1093/bioinformatics/btn043.
Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, Meyer F, Olsen GJ, Olson R, Osterman AL, Overbeek RA, McNeil LK, Paarmann D, Paczian T, Parrello B, Pusch GD, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O: The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008, 9: 75-10.1186/1471-2164-9-75.
McClelland M, Sanderson KE, Spieth J, Clifton SW, Latreille P, Courtney L, Porwollik S, Ali J, Dante M, Du F, Hou S, Layman D, Leonard S, Nguyen C, Scott K, Holmes A, Grewal N, Mulvaney E, Ryan E, Sun H, Florea L, Miller W, Stoneking T, Nhan M, Waterston R, Wilson RK: Complete genome sequence of Salmonella enterica serovar Typhimurium LT2. Nature. 2001, 413: 852-856. 10.1038/35101614.
Figueroa-Bossi N, Uzzau S, Maloriol D, Bossi L: Variable assortment of prophages provides a transferable repertoire of pathogenic determinants in Salmonella. Mol Microbiol. 2001, 39: 260-271. 10.1046/j.1365-2958.2001.02234.x.
Reen FJ, Boyd EF, Porwollik S, Murphy BP, Gilroy D, Fanning S, McClelland M: Genomic comparisons of Salmonella enterica serovar Dublin, Agona, and Typhimurium strains recently isolated from milk filters and bovine samples from Ireland, using a Salmonella microarray. Appl Environ Microbiol. 2005, 71: 1616-1625. 10.1128/AEM.71.3.1616-1625.2005.
Duchaud E, Rusniok C, Frangeul L, Buchrieser C, Givaudan A, Taourit S, Bocs S, Boursaux-Eude C, Chandler M, Charles J-F, Dassa E, Derose R, Derzelle S, Freyssinet G, Gaudriault S, Medigue C, Lanois A, Powell K, Siguier P, Vincent R, Wingate V, Zouine M, Glaser P, Boemare N, Danchin A, Kunst F: The genome sequence of the entomopathogenic bacterium Photorhabdus luminescens. Nat Biotechnol. 2003, 21: 1307-1313. 10.1038/nbt886.
Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, Marvig RL, Jelsbak L, Sicheritz-Ponten T, Ussery DW, Aarestrup FM, Lund O: Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol. 2012, 50: 1355-1361. 10.1128/JCM.06094-11.
Timme RE, Pettengill J, Allard MW, Strain E, Barrangou R, Wehnes C, Kessel JV, Karns J, Musser SM, Brown EW: Phylogenetic diversity of the enteric pathogen Salmonella enterica subsp. enterica inferred from genome-wide reference-free SNP characters. Genome Biol Evol. 2013, 5: 2109-2123. 10.1093/gbe/evt159.
Marttinen P, Hanage WP, Croucher NJ, Connor TR, Harris SR, Bentley SD, Corander J: Detection of recombination events in bacterial genomes from large population samples. Nucleic Acids Res. 2012, 40: e6-10.1093/nar/gkr928.
Moran NA, McLaughlin HJ, Sorek R: The dynamics and time scale of ongoing genomic erosion in symbiotic bacteria. Science. 2009, 323: 379-382. 10.1126/science.1167140.
Morelli G, Didelot X, Kusecek B, Schwarz S, Bahlawane C, Falush D, Suerbaum S, Achtman M: Microevolution of Helicobacter pylori during prolonged infection of single hosts and within families. PLoS Genet. 2010, 6: e1001036-10.1371/journal.pgen.1001036.
Zhou Z, McCann A, Litrup E, Murphy R, Cormican M, Fanning S, Brown D, Guttman DS, Brisse S, Achtman M: Neutral genomic microevolution of a recently emerged pathogen, Salmonella enterica serovar Agona. PLoS Genet. 2013, 9: e1003471-10.1371/journal.pgen.1003471.
Luo C, Walk ST, Gordon DM, Feldgarden M, Tiedje JM, Konstantinidis KT: Genome sequencing of environmental Escherichia coli expands understanding of the ecology and speciation of the model bacterial species. Proc Natl Acad Sci U S A. 2011, 108: 7200-7205. 10.1073/pnas.1015622108.
Den Bakker HC, Moreno Switt AI, Govoni G, Cummings CA, Ranieri ML, Degoricija L, Hoelzer K, Rodriguez-Rivera LD, Brown S, Bolchacova E, Furtado MR, Wiedmann M: Genome sequencing reveals diversification of virulence factor content and possible host adaptation in distinct subpopulations of Salmonella enterica. BMC Genomics. 2011, 12: 425-10.1186/1471-2164-12-425.
Sabbagh SC, Forest CG, Lepage C, Leclerc J-M, Daigle F: So similar, yet so different: uncovering distinctive features in the genomes of Salmonella enterica serovars Typhimurium and Typhi. FEMS Microbiol Lett. 2010, 305: 1-13. 10.1111/j.1574-6968.2010.01904.x.
Jacobsen A, Hendriksen RS, Aaresturp FM, Ussery DW, Friis C: The Salmonella enterica Pan-genome. Microb Ecol. 2011, 62: 487-504. 10.1007/s00248-011-9880-1.
Bishop AL, Baker S, Jenks S, Fookes M, Gaora PO, Pickard D, Anjum M, Farrar J, Hien TT, Ivens A, Dougan G: Analysis of the hypervariable region of the Salmonella enterica genome associated with tRNA(leuX). J Bacteriol. 2005, 187: 2469-2482. 10.1128/JB.187.7.2469-2482.2005.
Chaudhuri RR, Morgan E, Peters SE, Pleasance SJ, Hudson DL, Davies HM, Wang J, van Diemen PM, Buckley AM, Bowen AJ, Pullinger GD, Turner DJ, Langridge GC, Turner AK, Parkhill J, Charles IG, Maskell DJ, Stevens MP: Comprehensive assignment of roles for Salmonella Typhimurium genes in intestinal colonization of food-producing animals. PLoS Genet. 2013, 9: e1003456-10.1371/journal.pgen.1003456.
Haneda T, Ishii Y, Danbara H, Okada N: Genome-wide identification of novel genomic islands that contribute to Salmonella virulence in mouse systemic infection. FEMS Microbiol Lett. 2009, 297: 241-249. 10.1111/j.1574-6968.2009.01686.x.
Shi L, Adkins JN, Coleman JR, Schepmoes AA, Dohnkova A, Mottaz HM, Norbeck AD, Purvine SO, Manes NP, Smallwood HS, Wang H, Forbes J, Gros P, Uzzau S, Rodland KD, Heffron F, Smith RD, Squier TC: Proteomic analysis of Salmonella enterica serovar Typhimurium isolated from RAW 264.7 macrophages: identification of a novel protein that contributes to the replication of serovar Typhimurium inside macrophages. J Biol Chem. 2006, 281: 29131-29140. 10.1074/jbc.M604640200.
Bogomolnaya LM, Santiviago CA, Yang H-J, Bäumler AJ, Andrews-Polymenis HL: “Form variation” of the O12 antigen is critical for persistence of Salmonella Typhimurium in the murine intestine. Mol Microbiol. 2008, 70: 1105-1119. 10.1111/j.1365-2958.2008.06461.x.
De la Cruz MA, Zhao W, Farenc C, Gimenez G, Raoult D, Cambillau C, Gorvel J-P, Méresse S: A toxin-antitoxin module of Salmonella promotes virulence in mice. PLoS Pathog. 2013, 9: e1003827-10.1371/journal.ppat.1003827.
Katribe E, Bogomolnaya LM, Wingert H, Andrews-Polymenis H: Subspecies IIIa and IIIb Salmonellae are defective for colonization of murine models of salmonellosis compared to Salmonella enterica subsp. I serovar Typhimurium. J Bacteriol. 2009, 191: 2843-2850. 10.1128/JB.01223-08.
Wood MW, Jones MA, Watson PR, Siber AM, McCormick BA, Hedges S, Rosqvist R, Wallis TS, Galyov EE: The secreted effector protein of Salmonella Dublin, SopA, is translocated into eukaryotic cells and influences the induction of enteritis. Cell Microbiol. 2000, 2: 293-303. 10.1046/j.1462-5822.2000.00054.x.
Zhang S, Santos RL, Tsolis RM, Stender S, Hardt W-D, Bäumler AJ, Adams LG: The Salmonella enterica serotype Typhimurium effector proteins SipA, SopA, SopB, SopD, and SopE2 act in concert to induce diarrhea in calves. Infect Immun. 2002, 70: 3843-3855. 10.1128/IAI.70.7.3843-3855.2002.
Srikanth CV, Wall DM, Maldonado-Contreras A, Shi HN, Zhou D, Demma Z, Mumy KL, McCormick BA: Salmonella pathogenesis and processing of secreted effectors by caspase-3. Science. 2010, 330: 390-393. 10.1126/science.1194598.
Thomson NR, Clayton DJ, Windhorst D, Vernikos G, Davidson S, Churcher C, Quail MA, Stevens M, Jones MA, Watson M, Barron A, Layton A, Pickard D, Kingsley RA, Bignell A, Clark L, Harris B, Ormond D, Abdellah Z, Brooks K, Cherevach I, Chillingworth T, Woodward J, Norberczak H, Lord A, Arrowsmith C, Jagels K, Moule S, Mungall K, Sanders M, et al: Comparative genome analysis of Salmonella Enteritidis PT4 and Salmonella Gallinarum 287/91 provides insights into evolutionary and host adaptation pathways. Genome Res. 2008, 18: 1624-1637. 10.1101/gr.077404.108.
Holt KE, Thomson NR, Wain J, Langridge GC, Hasan R, Bhutta ZA, Quail MA, Norbertczak H, Walker D, Simmonds M, White B, Bason N, Mungall K, Dougan G, Parkhill J: Pseudogene accumulation in the evolutionary histories of Salmonella enterica serovars Paratyphi A and Typhi. BMC Genomics. 2009, 10: 36-10.1186/1471-2164-10-36.
Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH: CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39 (Database issue): D225-D229.
Diao J, Zhang Y, Huibregtse JM, Zhou D, Chen J: Crystal structure of SopA, a Salmonella effector protein mimicking a eukaryotic ubiquitin ligase. Nat Struct Mol Biol. 2008, 15: 65-70. 10.1038/nsmb1346.
Raffatellu M, Wilson RP, Chessa D, Andrews-Polymenis H, Tran QT, Lawhon S, Khare S, Adams LG, Bäumler AJ: SipA, SopA, SopB, SopD, and SopE2 contribute to Salmonella enterica serotype Typhimurium invasion of epithelial cells. Infect Immun. 2005, 73: 146-154. 10.1128/IAI.73.1.146-154.2005.
Osborne SE, Tuinema BR, Mok MCY, Lau PS, Bui NK, Tomljenovic-Berube AM, Vollmer W, Zhang K, Junop M, Coombes BK: Characterization of DalS, an ATP-binding cassette transporter for D-alanine, and its role in pathogenesis in Salmonella enterica. J Biol Chem. 2012, 287: 15242-15250. 10.1074/jbc.M112.348227.
Angiuoli SV, Gussman A, Klimke W, Cochrane G, Field D, Garrity G, Kodira CD, Kyrpides N, Madupu R, Markowitz V, Tatusova T, Thomson N, White O: Toward an online repository of Standard Operating Procedures (SOPs) for (meta)genomic annotation. OMICS. 2008, 12: 137-141. 10.1089/omi.2008.0017.
Den Bakker HC, Moreno Switt AI, Cummings CA, Hoelzer K, Degoricija L, Rodriguez-Rivera LD, Wright EM, Fang R, Davis M, Root T, Schoonmaker-Bopp D, Musser KA, Villamil E, Waechter H, Kornstein L, Furtado MR, Wiedmann M: A whole genome SNP based approach to trace and identify outbreaks linked to a common Salmonella enterica subsp. enterica serovar Montevideo Pulsed Field Gel Electrophoresis type. Appl Environ Microbiol. 2011, 77: 8648-8655. 10.1128/AEM.06538-11.
Chen Y, Souaiaia T, Chen T: PerM: efficient mapping of short sequencing reads with periodic full sensitive spaced seeds. Bioinformatics. 2009, 25: 2514-2521. 10.1093/bioinformatics/btp486.
Souaiaia T, Frazier Z, Chen T: ComB: SNP calling and mapping analysis for color and nucleotide space platforms. J Comput Biol. 2011, 18: 795-807. 10.1089/cmb.2011.0027.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.
Gill MS, Lemey P, Faria NR, Rambaut A, Shapiro B, Suchard MA: Improving Bayesian population dynamics inference: a coalescent-based model for multiple loci. Mol Biol Evol. 2013, 30: 713-724. 10.1093/molbev/mss265.
Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26: 841-842. 10.1093/bioinformatics/btq033.
We would like to thank Matt Ranieri for his help with DNA extraction. We further would like to thank Dr. Margaret Davis (Washington State University), Dr. Warnick’s laboratory (Cornell University) and the New York State Department of Health for providing isolates used in this study. Partial funding for this work has been provided by a USDA NIFA Special Research Grant to MW (2010-34459-20756).
LDR-R, AMS, MW and HDB declare that they have no competing interests. Life Technologies Corporation partially funded this study by providing sequencing reagents and instruments, and by compensating its employees (CAC, LD, RF and MRF), who participated in study design, data collection and analysis, decision to publish, and preparation of the manuscript.
LDR-R, AMS, CAC, MRF and MW conceived the study. AMS performed DNA isolation. LD and RF performed the genome sequencing. LDR-R, CAC, and HDB performed the genome sequence analysis. LDR-R performed the invasion assays. LDR-R, MW, and HDB wrote the paper. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Rodriguez-Rivera, L.D., Moreno Switt, A.I., Degoricija, L. et al. Genomic characterization of Salmonella Cerro ST367, an emerging Salmonella subtype in cattle in the United States. BMC Genomics 15, 427 (2014). https://doi.org/10.1186/1471-2164-15-427
- Multi Locus Sequence Typing
- Human Intestinal Epithelial Cell
- Gaussian Markov Random Field
- High Probability Density
- Prophage Region