Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Applied and Environmental Microbiology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Food Microbiology

Salmonella enterica subsp. enterica Serovar Heidelberg Food Isolates Associated with a Salmonellosis Outbreak Have Enhanced Stress Tolerance Capabilities

Andrea J. Etter, Alyssa M. West, John L. Burnett, Sophie Tongyu Wu, Deklin R. Veenhuizen, Raeya A. Ogas, Haley F. Oliver
Johanna Björkroth, Editor
Andrea J. Etter
aDepartment of Food Science, Purdue University, West Lafayette, Indiana, USA
bInterdisciplinary Life Sciences Program, Purdue University, West Lafayette, Indiana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Andrea J. Etter
Alyssa M. West
aDepartment of Food Science, Purdue University, West Lafayette, Indiana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John L. Burnett
aDepartment of Food Science, Purdue University, West Lafayette, Indiana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sophie Tongyu Wu
aDepartment of Food Science, Purdue University, West Lafayette, Indiana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Deklin R. Veenhuizen
aDepartment of Food Science, Purdue University, West Lafayette, Indiana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Raeya A. Ogas
aDepartment of Food Science, Purdue University, West Lafayette, Indiana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Haley F. Oliver
aDepartment of Food Science, Purdue University, West Lafayette, Indiana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Haley F. Oliver
Johanna Björkroth
University of Helsinki
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/AEM.01065-19
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Salmonella enterica serovar Heidelberg is currently the 12th most common serovar of Salmonella enterica causing salmonellosis in the United States and results in twice the average incidence of blood infections caused by nontyphoidal salmonellae. Multiple outbreaks of salmonellosis caused by Salmonella Heidelberg resulted from the same poultry processor, which infected 634 people during 2013 and 2014. The hospitalization and invasive illness rates were 38% and 15%, respectively. We hypothesized that the outbreak strains of Salmonella Heidelberg had enhanced stress tolerance and virulence capabilities. We sourced nine food isolates collected during the outbreak investigation and three reference isolates to assess their tolerance to heat and sanitizers, ability to attach to abiotic surfaces, and invasiveness in vitro. We performed RNA sequencing on three isolates (two outbreak-associated isolates and a reference Salmonella Heidelberg strain) with various levels of heat tolerance to gain insight into the mechanism behind the isolates’ enhanced heat tolerance. We also performed genomic analyses to determine the genetic relationships among the outbreak isolates. Ultimately, we determined that (i) six Salmonella Heidelberg isolates associated with the foodborne outbreak had enhanced heat tolerance, (ii) one outbreak isolate with enhanced heat tolerance also had an enhanced biofilm-forming ability under stressful conditions, (iii) exposure to heat stress increased the expression of Salmonella Heidelberg multidrug efflux and virulence genes, and (iv) outbreak-associated isolates were likely transcriptionally primed to better survive processing stresses and, potentially, to cause illness.

IMPORTANCE This study provides a deep analysis of the intrinsic stress tolerance and virulence capabilities of Salmonella Heidelberg that may have contributed to the length and severity of a recent salmonellosis outbreak. Additionally, this study provides a comprehensive analysis of the transcriptomic response of S. enterica strains to heat stress conditions and compares baseline stationary-phase gene expression among outbreak- and non-outbreak-associated Salmonella Heidelberg isolates. These data can be used in assay development to screen isolates for stress tolerance and subsequent survival. This study adds to our understanding of the strains associated with the outbreak and informs ongoing regulatory discussions on Salmonella in poultry.

INTRODUCTION

Salmonellae are frequently associated with poultry and eggs (1–3), though contact with live beef calves (4) and ingestion of contaminated pork (5), fruits and vegetables (6, 7), spices (8), and peanut butter (9) have all been implicated in salmonellosis outbreaks. However, Salmonella enterica serovars Heidelberg, Typhimurium, and Enteritidis are all commonly associated with poultry and eggs, especially chicken and turkey (10, 11). No serovar of S. enterica is currently considered an adulterant in raw poultry. However, to meet USDA performance standards, no more than 9.8% of young broilers (5/51 sampled), 25% of ground chicken samples, and 15.4% of chicken parts sampled can be positive for S. enterica (12). To achieve this benchmark, poultry processors employ a number of pathogen reduction steps during slaughter and processing, including the use of high temperatures and the use of antimicrobials in scald, chill, and wash water (13, 14). Nonetheless, these measures do not entirely eliminate S. enterica from the finished product.

Some isolates may be more likely to survive processing control strategies due to intrinsic characteristics. For instance, Salmonella Senftenberg 775W is extremely heat resistant (15), which would result in better survival of improper cooking or processing measures. Additionally, acid tolerance varies considerably among S. enterica strains, which could affect the efficacy of some antimicrobials (13, 15). Finally, S. enterica strains have widely different abilities to form biofilms, which are far more resistant to stress than planktonic cells, providing protection to the bacteria from mechanical stress, heat, sanitizers, and desiccation (16–24). Proper handling kills Salmonella. However, research indicates that few U.S. consumers use meat thermometers when cooking poultry, and many wash carcasses in the kitchen sink and share cutting boards between poultry and other foods (25–27). Consequently, bacteria which survive pathogen reduction measures during processing can cross-contaminate food, non-food-contact surfaces, and ready-to-eat foods, resulting in potential illness. In view of this, understanding the intrinsic factors that may support the survival of the pathogen from processing control strategies is crucial to improving their efficacy, especially for serovars causing invasive illnesses.

Recently, an outbreak of salmonellosis resulting from Salmonella Heidelberg in chicken underscored the importance of studying differences in stress tolerance among salmonellae. This outbreak was traced back to a single processor in the western United States and lasted for more than a year, despite the processor meeting and eventually surpassing performance standards (2, 28, 29). In fact, as the company increased efforts to reduce the S. enterica prevalence in its products, the percentage of S. enterica isolates detected overall dropped, but the percentage of S. enterica isolates belonging to the outbreak strains increased over almost 5 months before declining (29). Seven Salmonella Heidelberg strains were implicated in the outbreak, and several of these were multidrug resistant (MDR) (2). The outbreak had an above-average hospitalization rate and invasive illness rate for Salmonella Heidelberg (a 38% versus 26% average hospitalization rate and a 15% versus 13% average invasive illness rate), suggesting that the strains might have had above-average virulence (2, 30). Given the length and scope of the outbreak, we hypothesized that the outbreak strains of Salmonella Heidelberg may have had enhanced stress tolerance capacities.

The purpose of this study was to identify intrinsic Salmonella Heidelberg characteristics that may have contributed to the length, scope, and severity of the 2013-2014 Salmonella Heidelberg outbreak from chicken (2, 28, 29). To accomplish this, we assessed the heat tolerance, sanitizer tolerance, and attachment ability of nine food isolates associated with the outbreak and compared these characteristics to those of reference strains. We also evaluated the invasion capacity of four outbreak-associated food isolates and defined the transcriptomic profile of three Salmonella Heidelberg isolates (two outbreak-associated food isolates and strain SL476) during stationary phase and after heat shock to identify potential mechanisms of enhanced heat and stress tolerance. Finally, we compared the phenotypic and transcriptomic characteristics with epidemiological and genomic data to draw broader inferences about the isolates.

RESULTS

Genomic analysis of Salmonella Heidelberg outbreak-associated isolates.(i) Isolates contain multiple plasmids and resistance genes. All six isolates (R1-0001, R1-0002, R1-0003, R1-0006, R1-0007, and SL476) had evidence of at least one plasmid (Table 1). Reads mapping to ColpVC, a plasmid found in Salmonella Enteritidis (GenBank accession no. JX133088.1), were the most common, followed by reads mapping to an IncHI2/IncHI2A resistance plasmid from Serratia marcescens (GenBank accession no. BX664015.1). Four isolates had reads mapping to an IncI1 plasmid from Salmonella Typhimurium (GenBank accession no. AP005147.1), while two had reads mapping to an IncX1 plasmid from Escherichia coli (GenBank accession no. EU370913.1). Finally, R1-0002 had reads mapping to an IncA/C2 blaCYM-2 resistance plasmid from Klebsiella pneumoniae (GenBank accession no. JN157804.1). PlasmidFinder failed to recognize the second recognized SL476 plasmid; a BLAST analysis of the nucleotide sequences in the National Center for Biotechnology Information (NCBI) database identified it as matching Escherichia coli strain Sanji plasmid pSJ_3 (GenBank accession no. CP011066.1). All Salmonella Heidelberg isolates possessed the fosfomycin resistance gene (fosA7). Aminoglycoside resistance gene aac(6′)-Iaa was also present across the isolates (Table 1). However, its nucleotide sequence had only 97.3% similarity with the nucleotide sequence of the reference gene. Alignment of the isolate sequences with the reference nucleotide sequence revealed that their sequences all contained seven point mutations leading to the following amino acid substitutions: R8K, H10N, D12E, H42Y, V76A, T105A, and S117T. All of these are within the predicted acetyltransferase domain. Furthermore, all isolates contained a missense mutation in parC (T57S) associated with resistance to nalidixic acid and ciprofloxacin.

View this table:
  • View inline
  • View popup
TABLE 1

Genomic characteristics for seven Salmonella Heidelberg outbreak-associated isolates and SL476

(ii) Genomic analyses reveal isolate relationships. With a seven-gene multilocus sequence typing (MLST) scheme, all Salmonella Heidelberg isolates typed as MLST sequence type (ST) 15 (Table 1). However, using core genome MLST (cgMLST), the isolates all had different STs. A phylogenetic analysis of the seven outbreak isolates and SL476 demonstrated that isolate pairs from the same pulsed-field gel electrophoresis (PFGE) type generally clustered together (Fig. 1). When the strains were compared, isolates from PFGE type JF6X01.0041 were most closely related to JF6X01.0258 isolate R1-0006, while JF6X01.0045 isolates clustered most closely with JF6X01.0326 isolate R1-0007. SL476 (PFGE type unknown) was not closely related to any of the isolates. Analyzing the nucleotide difference (ND) matrix, there were substantial differences between isolates of the same PFGE strains (see Table S1 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507). There were 134 total NDs among the isolates, with 58 NDs appearing between isolates of JF6X01.0045 (R1-0002, R1-0004), while isolates from JF6X01.0041 (R1-0001, R1-0005) had 24 NDs and R1-0003 and R1-0006 (JF6X01.0258) had 22 NDs. The isolates differed from SL476 by 109 to 134 NDs. A larger analysis including food isolates from two other PFGE types involved in the outbreak (see Fig. S1 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507) showed that the isolates with PFGE patterns matching those of isolates involved in the 2013-2014 outbreak clustered more closely with each other than with SL476, regardless of their collection year.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Phylogenetic tree showing the relationships for the seven Salmonella Heidelberg outbreak-related isolates for which genomic sequences were publicly available as well as SL476. Relationships were assessed by use of the NDtree (v.1.2) phylogenetic analysis pipeline using a neighbor-joining algorithm. The tree was visualized using the FigTree (v.1.4.4) program.

The majority of Salmonella Heidelberg outbreak-associated isolates have enhanced heat tolerance in stationary phase.We assessed the heat tolerance of nine Salmonella Heidelberg food isolates associated with the outbreak as well as SL476 and Salmonella Typhimurium 14028 at 56°C, the temperature of a hard (i.e., high-heat) scald in poultry processing (Table 2) (31). There were no significant differences in heat tolerance among isolates in logarithmic phase (data not shown). However, six of the nine outbreak-associated Salmonella Heidelberg isolates showed significantly increased heat tolerance (as determined by a repeated-measures analysis of variance [ANOVA] for the survivor curves) compared to that of Salmonella Typhimurium 14028 (P < 0.05). R1-0006 and R1-0007 showed significantly higher heat tolerance than Salmonella Heidelberg SL476 and Salmonella Typhimurium 14028 (P < 0.05) in stationary phase (Fig. 2).

View this table:
  • View inline
  • View popup
TABLE 2

Characteristics of isolates studied

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Survival of Salmonella Heidelberg outbreak-related isolates and S. enterica reference strains at 56°C during stationary phase. Error bars represent the standard deviation from at least three biological replicates.

There was no association between antibiotic resistance and enhanced heat tolerance (see Table S2 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507). R1-0007 is a pansusceptible isolate, while R1-0006 is resistant to kanamycin, streptomycin, and tetracycline (Table 2). Antibiotic-resistant isolates R1-0001, R1-0002, and R1-0005 did not have enhanced heat tolerance, while antibiotic-resistant isolates R1-0003, R1-0006, and R1-0009 and pansusceptible isolates R1-0004, R1-0007 and R1-0008 had significantly higher heat tolerance than Salmonella Typhimurium ATCC 14028 (P < 0.05), and R1-0006 and R1-0007 had significantly higher heat tolerance than SL476 (P < 0.05).

Salmonella Heidelberg outbreak-associated isolates attach to an abiotic surface better under stress conditions.We assessed the attachment capability of 10 Salmonella Heidelberg isolates, 4 Salmonella Typhimurium isolates, 4 Salmonella Newport food isolates, Salmonella Typhimurium ATCC 14028, and Salmonella Enteritidis at 21°C (in 1× tryptic soy broth [TSB]; see Fig. S2 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507] and 1/20× TSB) and at 4°C. At 21°C in 1/20× TSB (Fig. 3A), serovar and time were significant (P < 0.05). Salmonella Heidelberg attached significantly less than Salmonella Typhimurium, Enteritidis, and Newport (adjusted P value [Padj] < 0.05). Attachment efficacy increased significantly as time increased from 24 h to 120 h (Padj < 0.05). At the isolate level, Salmonella Enteritidis isolate R1-0010 had the highest attachment, which was significantly higher than that of R1-0001, R1-0002, R1-0003, R1-0005, and SL476 (Padj < 0.05). Salmonella Heidelberg R1-0002 had the lowest attachment, which was significantly lower than that of all the isolates except SL476 and R1-0005 (Padj < 0.05). However, heat-tolerant Salmonella Heidelberg isolates R1-0006 and R1-0007 did not attach significantly less than Salmonella Enteritidis R1-0010. When this experiment was carried out to 9 days, we found similar results (see Fig. S3 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Attachment levels of the S. enterica isolates on a polystyrene microtiter plate at 21°C in 1/20× TSB (A) and 4°C in 1/20× TSB (B), as determined by a crystal violet microtiter assay. Letters above the isolates indicate statistically significant differences among the serovars. Attachment levels of the S. enterica isolates were determined at days 1, 3, and 5. The attachment levels at 21°C in 1× TSB at days 1, 3, and 5, and 21°C in 1/20× TSB at days 5, 7, and 9 are available in Fig. S2 and S3 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507.

At 4°C in 1/20× TSB (Fig. 3B), serovar, time, and serovar × time were significant (P < 0.05). Salmonella Typhimurium had the highest attachment among the serovars, while Newport had the lowest (Padj < 0.05). Salmonella Heidelberg had significantly lower attachment than Salmonella Typhimurium, but its attachment was not significantly different from that of Salmonella Enteritidis or Newport (Padj < 0.05). Salmonella Enteritidis attached significantly better than Salmonella Newport (Padj < 0.05). Attachment at 120 h was significantly lower than at both 24 h and 72 h (Padj < 0.05). At the isolate level at 4°C, heat-tolerant isolate Salmonella Heidelberg R1-0006 had the highest attachment among all the isolates, while SL476 had the lowest attachment, and Salmonella Typhimurium 14028 had attachment comparable to that of most Salmonella Heidelberg and Newport isolates (Padj < 0.05).

Salmonella Heidelberg food isolates do not have enhanced tolerance to a PAA-based sanitizer or increased in vitro invasiveness.Isolates rarely survived or grew in the presence of peroxyacetic acid (PAA) at concentrations greater than 50 ppm, indicating that they were sensitive to the sanitizer at label use concentrations (see Table S3 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507).

We intended to investigate the in vitro virulence of all the outbreak-associated isolates, in addition to the reference isolates (12 total isolates), in a Caco-2 cell model; however, 5 Salmonella Heidelberg outbreak-associated strains grew in up to 250 μg/ml gentamicin when incubated in Dulbecco modified Eagle medium with 10% fetal bovine serum. Consequently, we tested resistance to amikacin disulfate, an aminoglycoside to which the isolates were anticipated to be sensitive. Seven of the 12 isolates were resistant to 250 μg/ml amikacin disulfate (data not shown). Consequently, we proceeded with standard gentamicin protection assays for the six gentamicin-susceptible isolates and omitted the other six from our in vitro studies. We were able to assess the invasiveness of four Salmonella Heidelberg outbreak-associated isolates; two isolates were resistant to both gentamicin and amikacin and therefore were not assessed. We also assessed reference strains SL476 and 14028. While strain was a significant factor in the model (P = 0.049), there were no significant differences in mean invasiveness among the six strains after adjusting for multiple comparisons. Invasiveness ranged from approximately 3% for R1-0009 to just over 1.2% for R1-0006 (see Fig. S4 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507).

Heat stress increases expression of some HSPs, virulence, and stress responses genes in Salmonella Heidelberg R1-0006, R1-0007, and SL476.(i) Increased expression of HSPs. The heat shock protein (HSP) genes hscBA were upregulated in all three isolates, but more so in SL476 than in R1-0006 and R1-0007 (Fig. 4). Conversely, stress response regulator rpoH was expressed at a higher level in R1-0006 and R1-0007 than in SL476. Only R1-0006 showed significantly increased expression of the heat shock protein gene dnaJ, while R1-0007 was the only isolate with significant increases in the expression of ibpB. SL476 showed lower differential expression of the heat shock protein gene ibpA; R1-0006 and R1-0007 showed less differential expression of the genes encoding heat shock proteins GroS, HslR, and SigE. All three isolates showed significantly increased expression of the negative regulators of SigE, RseBC and stress-associated chaperone ClpP. R1-0006 and R1-0007 also had significantly decreased expression of sigE and increased expression of negative regulator rseA.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Changes in expression of genes encoding heat shock proteins (HSPs) and stress tolerance proteins after 5.1 min at 56°C (7.5 min, including the time to temperature). Results are the log2 fold change (FC) of gene expression at 56°C compared to that at 37°C for each isolate. Significant differential gene expression was defined as a 1-log2 change and a Padj value of <0.05.

Among the genes encoding stress tolerance systems, most genes were significantly less expressed during heat stress (Fig. 4), including genes for osmotic stress, cold shock, oxidative stress, and acid stress. Expression of cold shock protein genes cspAE was decreased among all three isolates, and cspC expression was decreased in SL476 and R1-0006. However, R1-0006 had increased expression of phage shock protein genes pspABDE, fusaric acid resistance gene SeHa_C1611, and stress-associated gene deaD; R1-0007 also had increased expression of pspE. Both R1-0006 and R1-0007 had increased expression of general stress protein gene ycfR (bhsA). R1-0007 also had increased expression of copper resistance protein gene cutP, the gene for a manganese-containing catalase (SeHa_C1921), osmotic stress gene putP, and cobalt resistance gene apaG.

(ii) Heat stress induces increased expression of virulence genes. During heat shock, multiple genes associated with S. enterica virulence showed significantly increased expression among the three isolates (Fig. 5). Salmonella pathogenicity island 1 (SPI-1) genes invI, epaR, and yopJ (avrA) had significantly increased expression across all three isolates. Additionally, SPI-1 genes invEGIJ, epaOQR, spaS, and orgB had significantly increased expression in SL476. In R1-0006, SPI-1 genes invGI, epaR, and orgB showed significantly increased expression; however, invH, sicA, prgHIJ, sopB, and sopE2 all had significantly decreased expression compared to that during stationary phase. In R1-0007, SPI-1 genes epaQR and invI also showed significantly increased expression, but invFH, sicA, sipAB, iagB, prgHJ, sopB, sopE, and sopE2 all had significantly decreased expression under heat stress. Multiple SPI-2 genes had increased expression among the three isolates (Fig. 5), but the patterns of gene expression varied substantially. SL476 showed the greatest number of SPI-2 genes with significantly increased expression; expression of sseDEFGJ, and ssaCDEJKLOV significantly increased. In R1-0006, the genes sifB, ssaT, and sseDFG showed significantly increased expression. R1-0007 had the fewest SPI-2 genes with significantly increased expression; sseF had significantly higher expression under heat stress, while ssaG and pipB2 had significantly lower expression. Other secretion systems were also affected (Fig. 5). SPI-5 gene pipA had increased expression across all three isolates. SL476 also showed significantly increased expression of SPI-3 gene mgtC and SPI-5 genes pipAB. Additionally, four type VI secretion system (T6SS) genes had significantly increased expression across all isolates.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Changes in virulence gene expression after 5.1 min at 56°C (7.5 min, including the time to temperature). Results are the log2 fold change of gene expression at 56°C compared to that at 37°C for each isolate. Significant differential gene expression was defined as a 1-log2 change and a Padj value of <0.05.

Three other putative virulence genes—sinH, nleB, and the gene for a putative attachment/invasion locus (SeHa_C0439)—had significantly increased expression during heat stress. Conversely, invasion response regulator uvrY and a virK homologue (SeHa_C1038) showed significantly decreased expression across all three isolates. SL476 and R1-0006 also showed significantly reduced expression of cold shock protein genes cspCE, which have been shown to affect virulence (33), and in SL476 and R1-0007 the expression of another copy of virK (SeHa_C2956) was also significantly decreased. R1-0007 showed significantly reduced expression of virK but significantly increased expression of putative virulence gene mviM.

(iii) Heat stress upregulates multidrug efflux pumps. During heat stress, the gens for multiple antibiotic resistance or multidrug efflux pumps were differentially expressed in all three isolates (Fig. 6). The patterns of expression for antibiotic resistance genes varied across the isolates. Streptomycin resistance gene strB was expressed at a significantly higher level in SL476, while mig-14, which may be involved in polymyxin resistance (34), was expressed at significantly lower levels. Conversely, the expression of streptomycin resistance gene strA was significantly decreased in R1-0006; however, SeHa_C2520, potentially encoding antibiotic resistance gene norA, had significantly increased expression. Multidrug efflux pump genes macAB and mdfA had significantly increased expression across all isolates. The genes mdtB, sbmA, acrD, dinF, and envR had significantly increased expression in SL476, while the expression of mdtHIJK and acrR was significantly decreased under heat shock conditions. In R1-0006, multidrug transporter genes mdtG, dinF, emrB, and envR and acriflavine/quaternary ammonium efflux pump gene smvA showed significantly increased expression under heat stress, but expression of the master efflux pump gene tolC, as well as mdtHJ, acrR, and marR, was significantly decreased. MDR efflux pump genes mdtG, mdtM, mdtQ, dinF, emrB, envR, SeHa_C0564, SeHa_C2503, and smvA also had significantly increased expression in R1-0007, but marAR, acrR, mdtJK, mprA, and ramA had significantly lower expression than they did under control conditions.

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Changes in expression of antibiotic resistance genes and multidrug efflux pumps after 5.1 min at 56°C (7.5 min, including the time to temperature). Results are the log2 fold change of gene expression at 56°C compared to that at 37°C for each isolate. Significant differential gene expression was defined as a 1-log2 change and a Padj value of <0.05.

Heat-tolerant isolates R1-0006 and R1-0007 have higher expression of stress response systems, virulence genes, and multidrug efflux pumps in stationary phase than SL476.(i) Increased expression of stress response genes. HSP gene hslR expression was increased in both R1-0006 and R1-0007 during stationary phase (Fig. 7) compared to that in SL476 during stationary phase. In R1-0006, HSP genes hspQ, hslR, dnaJ, clpB, and clpX were also expressed at significantly higher levels. Acid stress chaperone gene hdeB was also expressed at significantly higher levels. R1-0007 showed significantly higher expression of HSP genes hslR and hslU, while the sigE negative regulator gene rseA and HSP gene cpxP were expressed at significantly lower levels than in SL476 during stationary phase. Putative disulfide bond chaperone gene dsbI and heat-induced spheroplast protein Y gene spy were also expressed at significantly lower levels in R1-0007.

FIG 7
  • Open in new tab
  • Download powerpoint
FIG 7

Differences in expression of genes encoding HSPs and stress proteins during stationary phase at 37°C for each isolate compared to SL476 gene expression levels at the same temperature. Significant differential gene expression was defined as a 1-log2 change and a Padj value of <0.05.

A number of stress-related genes were expressed at significantly higher levels in R1-0006 during stationary phase (Fig. 7), including general stress response protein gene ycfR (bhsA), phosphate starvation protein gene psiF, two copies of phage shock protein gene pspD, and oxidative stress protein genes soxR and dps. Acid stress response genes aaeAB were also expressed at higher levels. Multiple genes involved in osmotic stress (e.g., osmBCYW) were expressed at significantly higher levels. Other genes expressed at significantly higher levels included the gene encoding paraquat-inducible protein B (pqiB_1), comA (ydeJ, SeHa_C1685), and the general stress protein 39 gene (SeHa_C3405). Cold shock protein gene cspD, universal stress protein gene uspF, phosphate starvation protein gene psiB, and copper resistance gene copD were all expressed at significantly lower levels in R1-0006 than in SL576.

Fewer stress-related genes were differentially expressed in R1-0007 in stationary phase compared to SL476 gene expression levels (Fig. 7). oxyR regulon gene dps, phosphate starvation stress protein gene psiF, phage shock protein gene pspD, a putative secondary copy of the pspD gene, and comA (ydeJ, SeHa1685) were expressed at significantly higher levels in R1-0007. Conversely, cold shock protein gene cspD, telluride resistance gene tehB, and phosphate starvation protein gene psiB were downregulated.

(ii) Outbreak-associated isolates have higher expression of virulence genes than SL476 in stationary phase. Multiple SPI-1 and SPI-5 genes were expressed at significantly higher levels in R1-0006 and R1-0007 than in SL476 in stationary phase (Fig. 8). Several fimbrial genes were also expressed at significantly higher levels. In contrast, most identified type IV secretion system (T4SS) genes were expressed at significantly lower levels. SPI-1 virulence genes expressed at significantly higher levels in R1-0006 than in SL476 included spaOPR, sicP, sicA, prgIH (eprIH), prgJ, invAEH, iacP, and effector genes sipAB and sopAB; SPI-1 effector genes sopE and invR (SeHa_C3091) were expressed at significantly lower levels (Fig. 8). SPI-2 gene ssaEJ, SeHa_C1543, and effector genes pipB and pipB2 were expressed at significantly higher levels in R1-0006, while ssaT, ttrB, and pipA were expressed at significantly lower levels. Four putative T6SS (SPI-6) genes were differentially expressed. ompA, sciS (icmF, SeHa_C0323), and a gene matching virG (SeHa_C0348) were expressed at significantly higher levels, while an hcp gene (SeHa_C3375) was expressed at significantly lower levels. Additionally, virulence-associated genes zinc transcriptional regulator gene zur, surface-exposed virulence protein gene bigA, virulence genes nleB, phoP, and virK, and virulence-associated flagellar genes srfABC were expressed at significantly lower levels, while psiF was expressed at significantly higher levels.

FIG 8
  • Open in new tab
  • Download powerpoint
FIG 8

Differences in expression of virulence genes during stationary phase at 37°C for each isolate compared to SL476 gene expression levels at the same temperature. Significant differential gene expression was defined as a 1-log2 change and a Padj value of <0.05.

Fewer SPI-1 virulence genes were expressed at significantly higher levels in R1-0007 than in R1-0006 compared to their levels of expression in SL476 (Fig. 8). SPI-1 genes expressed at significantly higher levels included spaP (epaP), sitAB, sicP, sicA, invAGHIJR, spaS, and effector genes sipAB and sopAB. SPI-2 genes ssaCEKQ and effector genes pipB and pipB2 were all expressed at significantly higher levels. ompA and sciS (icmF, SeHa_C0323), involved in expression of T6SS genes, and mgtAC and slsA of SPI-3 were also expressed at significantly higher levels in R1-0007. Finally, R1-0007 had significantly higher expression of virulence-associated genes nleB and psiF and significantly lower expression of virulence regulator gene phoP and virulence genes srfABC and mviM.

(iii) R1-0006 has higher stationary-phase expression of multidrug efflux pumps than SL476. pmrD, involved in resistance to polymyxin (35), was expressed at significantly lower levels during stationary phase in R1-0006 than in SL476, as was aminoglycoside resistance gene aac(6′)1aa (Fig. 9). R1-0006 had significantly higher levels of expression of multidrug efflux pump genes mdtAB, mdtG, mdtK, baeS, and acrB than SL476, while acrR was expressed at significantly lower levels. Only the efflux pump gene emrD was expressed at significantly lower levels in R1-0007 than in SL476. marAR was expressed at significantly higher levels in both isolates R1-0006 and R1-0007 than in SL476, and multidrug-efflux pump gene mdtIJ and the putative multidrug efflux gene SeHa_C1561 were expressed at significantly lower levels in the two isolates than in SL476.

FIG 9
  • Open in new tab
  • Download powerpoint
FIG 9

Differences in expression of genes encoding antibiotic resistance genes and multidrug efflux pumps during stationary phase at 37°C for each isolate compared to SL476 gene expression levels at the same temperature. Significant differential gene expression was defined as a 1-log2 change and a Padj value of <0.05.

DISCUSSION

Genotypic analyses.Determining genetic differences among the isolates was tricky. A nucleotide difference (ND)-based phylogenetic analysis showed that isolate pairs from the same PFGE type typically clustered together, while cgMLST (36) divided all the isolates into unique STs. A search of the EnteroBase database (36) for the STs in this study found that SL476, R1-0001, R1-0002, R1-0004, R1-0005, R1-0006, and R1-0007 were the only isolates in their STs, while there was just one other isolate in R1-0003’s ST. This suggests that the cgMLST classification system employed by EnteroBase may be overly stringent.

Phylogenetic analyses better illuminated the differences between the isolates. The ND matrix showed the largest differences (58 NDs) between isolates with PFGE type JF6X01.0045; the data indicate that R1-0002 is an MDR isolate, while R1-0004 is pansusceptible. The isolates with PFGE pattern JF6X01.0041 (R1-0001, R1-0005) had identical antibiotic resistances but differed by 24 single nucleotide polymorphisms (SNPs), while R1-0003 and R1-0006 (JF6X01.0258) differed by 22 NDs, with R1-0006 being resistant to only three of the five antimicrobials that R1-0003 had resistance to. Overall, the data indicated that isolate pairs from the same strain were not necessarily highly related. This is not surprising, as the outbreak involved multiple processing establishments, which were in turn sourced from multiple farms (28, 29). Hoffmann et al. found as few as 17 SNPs between isolates from the 2011 outbreak associated with ground turkey (37, 38), while Keefer et al. found that retail and clinical Salmonella Heidelberg isolates with the same PFGE profile had 35 SNPs different among them (39). This makes JF6X01.0045 isolates unusually dissimilar for isolates from the same strain within an outbreak, while isolate pairs from strains JF6X01.0258 and JF6X01.0041 differ within normal parameters. While SL476 was quite dissimilar from the outbreak-associated isolates, the NDtree phylogenetic analysis pipeline determined it to be the most appropriate reference genome for the isolates, affirming it as a good reference for this study.

All isolates, regardless of their phenotypic antibiotic resistance profile, contained aminoglycoside resistance gene aac(6′)-Iaa, which codes for a deactivator of kanamycin and amikacin (40). However, this gene does not provide clinical-level resistance to kanamycin or amikacin (40); in fact, this class of aminoglycoside resistance genes is often not expressed in S. enterica (41, 42). Consequently, its presence alone does not adequately explain the high levels of amikacin resistance that we uncovered when attempting to use amikacin for cell invasion assays. However, alignment of the isolate sequences with the reference sequence revealed that the isolates shared seven point mutations leading to amino acid substitutions. These mutations may have increased amikacin inactivation, especially as several mutations changed the charge type or size of the amino acid R group. Additionally, this gene was expressed in SL476 and R1-0007 and at low levels (mean, 25 copies) in R1-0006 as well. This may explain the resistance to amikacin that we uncovered when attempting to use it for invasion assays.

All the isolates also had a T57S mutation in topoisomerase, encoded by parC. This mutation is associated with resistance to ciprofloxacin and nalidixic acid, which none of the isolates were phenotypically resistant to. However, the literature indicates that gyrA mutations are primarily responsible for resistance to these antimicrobials (43–45). The T57S mutation in the protein encoded by parC does not seem to yield clinical levels of fluoroquinolone resistance, and it is rarely found in strains without gyrA mutations (42, 45, 46). One study remarked that certain parC mutations were associated with specific serovars of S. enterica (45). However, Keefer et al. did not find it among the 12 Salmonella Heidelberg isolates that they studied (39), suggesting that it is not specifically associated with Salmonella Heidelberg as a serovar. Consequently, the significance of this mutation in our isolate set is unclear. However, this underscores that antibiotic resistance data from genomic analyses must be interpreted with caution in the absence of phenotypic data.

Enhanced heat tolerance may contribute to processing intervention survival and therefore to outbreak scope and duration.The heat tolerance of SL476 was similar to that found in published data for Salmonella Enteritidis, but with more extreme tailing of the survival curve (15, 47, 48), while Salmonella Typhimurium ATCC 14028 had a somewhat lower heat tolerance than SL476 (not statistically significant). R1-0006 and R1-0003 were from strain JF6X01.0258, a strain that caused 186/634 outbreak illnesses, yet represented less than 1% of the Salmonella Heidelberg isolates in PulseNet prior to the start of the outbreak (28). This strain was also associated with illnesses from rotisserie chicken (28). Consequently, R1-0006’s enhanced heat tolerance may be epidemiologically important.

Several recent studies have investigated the transcriptional response of S. enterica to stress. Pin et al. exposed Salmonella Typhimurium strain 4/74 stationary-phase cultures to 45°C or pH 5 for 30 min and compared the microarray data from RNA collected (i) during exposure to these conditions, (ii) immediately after reversal to normal conditions, and (iii) at 30 min postreversal with those for control cultures exposed to 25°C and pH 7 (49). They found that acidification changed the transcriptional levels of only a few genes, while heat stress caused widespread transcriptional changes (49). They also found that 85% of the genes upregulated by heat shock remained upregulated 30 min after the reversal to 25°C and that this was associated with increased tolerance to further heat stress or acidification (49). This agrees with our results, which also found global transcriptional changes after heat shock. While R1-0006 and R1-0007 had varied upregulation of HSPs and stress response systems during heat shock, they also had significantly higher (≥2-fold increase, P < 0.05) expression of HSPs and stress proteins in stationary phase at 37°C, potentially priming for adaptation to and survival during heat stress. Sirsat et al. investigated the transcriptome of mid-log-phase Salmonella Typhimurium ATCC 14028 after 10 min at 48°C using a partial microarray (1,152 genes assessed) and also found that several heat shock proteins and stress proteins were upregulated; however, due to the different growth phase and temperature, these results are not directly comparable to our data (50).

Milillo et al. studied the transcriptional response of stationary-phase Salmonella Typhimurium LT2 to heated sodium propionate or sodium acetate (2.5% antimicrobial, pH 4, 55°C) for 1 min (51). They found that multiple genes for heat shock proteins or molecular chaperones were downregulated (51). The authors attributed this to the presence of membrane damage, citing a study which found inhibition of heat shock proteins following membrane damage (51). These results are somewhat similar to the heat shock results that we found for our isolates and suggest that we may have been observing some membrane damage in our isolates. This is not unexpected, given the temperature (56°C). While we did not find data from previous work to which to compare our stationary-phase transcriptomic data, stationary phase is known to induce stress tolerance within S. enterica (52, 53). Therefore, it may be that our more stress-tolerant strains activate this general stress response to stationary phase at higher levels and that this provides cross protection to other stresses.

Salmonella Heidelberg is not considered exceptionally heat tolerant among the salmonellae (15), which is further supported by the fact that SL476 had heat tolerance similar to that of Salmonella Enteritidis (47). Lianou and Koutsoumanis found no significant differences in heat tolerance among S. enterica serovars in a study that included Salmonella Heidelberg, Enteritidis, Typhimurium, and Senftenberg (15). Humpheson et al. noted that within heat-stressed monocultures of Salmonella Enteritidis, a subpopulation of bacteria possessed enhanced heat tolerance, resulting in a survivor tail (47). This suggests that the heat tolerance is often a phenotypic phenomenon and that the isolates that we studied simply had a higher subpopulation of heat-tolerant survivors. In contrast, Nguyen et al. recently identified a heat resistance island in Salmonella Senftenberg, indicating a genotypic basis for heat tolerance in some salmonellae (54). We did not find this island in our heat-tolerant isolates, which is unsurprising, as to date it has been identified only in Salmonella Senftenberg (54). We also did not investigate promoter sequences for mutations which might have led to enhanced gene expression, but we recognize that this could contribute to differential expression. Further, some genes in S. enterica are known to have multiple transcriptional start sites.

Heat stress does not consistently enhance virulence gene expression.Heat stress has been observed to increase virulence gene expression under some circumstances in S. enterica. When Sirsat et al. investigated the transcriptome of mid-log-phase Salmonella Typhimurium ATCC 14028 after 10 min at 48°C using a partial microarray, they found that a number of genes involved in virulence were also upregulated (50). Conversely, when Milillo et al. studied the transcriptional response of stationary-phase Salmonella Typhimurium to heat and acid stress by exposing Salmonella Typhimurium LT2 to heated 2.5% sodium propionate or sodium acetate, pH 4, at 55°C for 1 min (51), they found that transcription of virulence genes was not increased, a result that they attributed to membrane damage (51). Our results had some similarities with those of both studies, with upregulation of some virulence genes (including genes from SPI-1 and SPI-2) under heat stress but with a pattern which varied across isolates. Consequently, the effect of heat stress on virulence gene expression may be isolate dependent.

Stress tolerance may provide an adaptive advantage for attachment under adverse conditions but is not associated with sanitizer tolerance.Attachment ability is known to be highly variable among and within S. enterica serovars. Previous investigators detected serovar-specific differences in attachment; however, they also found high variation within isolates of the same serovars (16, 24). Keelara et al. compared the attachment of Salmonella Typhimurium, Infantis, and Derby strains in a standard crystal violet microtiter assay and found that fewer Salmonella Typhimurium strains than strains of the other two serovars exhibited moderate or strong attachment (55). Conversely, our Salmonella Typhimurium isolates had the highest overall attachment capacity in this study. Studies of Salmonella Heidelberg have generally found it to be an average biofilm former (56, 57). In contrast, our data indicated that all Salmonella Heidelberg isolates except R1-0001 and R1-0005 were weak or nonattaching biofilm formers, using the categories and conditions defined by Stepanović et al. (24 h of attachment at 35°C; we tested the isolates at 37°C; data not shown) (16). However, the fact that R1-0006 demonstrated greater attachment than the other isolates under stressful conditions (e.g., cold, nutrient-poor conditions) may indicate that at least some of the outbreak isolates are better able to form biofilms as a response to certain stressful conditions, such as a chiller in a processing plant or a grocery store cold case.

As enhanced sanitizer tolerance could increase S. enterica survival in the poultry processing environment, we tested Salmonella Heidelberg, Typhimurium, and Newport food isolates, as well as two reference strains, for increased tolerance to a PAA-based sanitizer (see Table 2 for the complete isolate list). None of the isolates that we studied showed tolerance to PAA-based sanitizers (see Table S3 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507). Some Salmonella Typhimurium isolates can develop tolerance to peracetic acid compounds (58), and studies in other pathogens indicate that multidrug-resistant pathogens may be more likely to be sanitizer tolerant (59, 60). However, it is crucial to note that the conditions under which we studied sanitizer tolerance did not contain organic soils common in the poultry processing environment, which can reduce sanitizer efficacy (61). Furthermore, studying the sanitizer tolerance of biofilms from these isolates would also have been useful; however, it was beyond the scope of this study.

Isolates with enhanced tolerance to one stress are more likely to tolerate other stresses.In our study, heat-tolerant isolate R1-0006 showed the highest attachment levels of any isolate in vitro under conditions of cold and nutrient stress. R1-0006 also demonstrated higher expression of heat shock proteins, general stress systems, osmotic stress systems, and some MDR efflux systems. Ultimately, these data suggest that the phenotypic characteristics of this isolate contributed to its fitness and higher overall prevalence than the other S. enterica strains and serovars during the first three phases of the Food Safety and Inspection Service (FSIS) monitoring period (29). Multiple studies have shown that enhanced tolerance to one stress may cross protect the organism against other stresses. For instance, Humphrey et al. found that Salmonella Enteritidis isolates that were more tolerant of heat, acid, or hydrogen peroxide also better survived desiccation (62). Others have demonstrated that the preexposure of Salmonella Enteritidis (63) and Typhimurium (49) to sublethal temperatures increased their tolerance to both heat and acid in subsequent exposures, indicating an overlap in stress protection systems. Hamilton et al. studied the transcriptome of Salmonella Typhimurium in a nonstatic biofilm system and found that, in addition to upregulation of attachment genes, Salmonella Typhimurium upregulated genes involved in the heat and oxidative stress response (64).

Heat stress may provide cross protection against antimicrobial stress.In addition to showing significantly increased expression of heat shock proteins and stress genes, all Salmonella Heidelberg isolates had higher expression of efflux pumps under heat stress. This may indicate that these isolates are more likely to develop antimicrobial tolerance under high-stress conditions, especially since overexpression of mdfA can produce resistance to tetracycline, chloramphenicol, doxorubicin, and norfloxacin, macAB overexpression provides resistance to erythromycin, and the emrAB system provides resistance to nalidixic acid (65). In agreement with this, the AcrAB-TolC antibiotic efflux system has been shown to be upregulated in response to bile and oxidative stress (66, 67). Additionally, macAB have been shown to play a role in virulence in Salmonella Typhimurium and are induced by conditions inside macrophages (65). Consequently, upregulation of these pumps during exposure to stress could increase virulence and antibiotic tolerance or even resistance.

Epidemiological and in vitro data do not support enhanced virulence among outbreak strains.Given the length and scope of the outbreak (28, 29), we hypothesized that the outbreak-associated Salmonella Heidelberg isolates might also have increased virulence; however, we found no evidence of increased invasiveness among the isolates tested. This is similar to the findings in other studies of increased invasiveness in Salmonella Typhimurium DT104, a multidrug-resistant strain associated with increased hospitalization risks (68, 69). Studies with other model systems may be warranted to more effectively evaluate virulence potential.

Whether antibiotic-resistant isolates are more virulent than pansusceptible isolates is disputed (68, 70). In this outbreak, the strains that caused the most hospitalizations and highest blood infection percentages displayed MDR phenotypes (28, 29). Conversely, apart from strain JF6X01.0672, which caused two cases of disease (neither of which resulted in blood infections), the lowest percentage of blood infections was found for pansusceptible strains (or strains for which most isolates were pansusceptible) (28, 29). Due to the high gentamicin resistance of our MDR isolates, we were unable to obtain in vitro results for comparison; however, epidemiological data from this outbreak (28) indicate that MDR isolates of Salmonella Heidelberg may be more invasive.

Data from PulseNet listing the numbers of isolates of each strain (i.e., PFGE PulseNet pattern) derived from blood, stool, and other sites in patients during the outbreak showed that there were substantial differences in blood infection percentages among the seven strains, with the percentage of blood infections varying from 0% to 17% (72). Of the six strains represented in our study, four had blood infection percentages of between 14 and 17% (see Table S2 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507) (72).

The strain responsible for the highest number of outbreak illnesses was JF6X01.0258 (R1-0003, R1-0006) (28); this strain was also associated with increased heat tolerance in our in vitro assays. However, it did not have a high hospitalization rate (31% versus an average for this outbreak of 38%) and had an expected invasive illness rate (14% versus an average of 15% for this outbreak) (72). Heat-tolerant strain JF6X01.0326 (R1-0007) had a 36% hospitalization rate but only a 6% invasive illness rate (28, 72). Finally, heat-susceptible strain JF6X01.0041 (R1-0001, R1-0005) caused 103 illnesses and had a 49% hospitalization rate and a 16% invasive illness rate (28). Invasive illnesses rates (for which the percentage of blood isolates is a proxy measurement [30]) are likely a more accurate epidemiological measure of isolate virulence than hospitalization rates, as hospitalization may be affected by factors unrelated to isolate virulence, such as patient age, climate, or poor hydration management. However, regardless of the measure used, there does not appear to be a strong relationship between heat tolerance and virulence potential.

Summary.Combined, these data suggest that some outbreak-associated Salmonella Heidelberg isolates may be uniquely primed to survive processing-type stresses, but considerable variation existed among the isolates. Data from this work will inform future risk assessments and industry best practices for Salmonella Heidelberg and other S. enterica serovars and may contribute to the development of gene expression-based assays to identify isolates more likely to survive processing stresses.

MATERIALS AND METHODS

Strains used in this study.The strains sourced and used in this study are described in Table 2. All strains were stored in 20% glycerol at −80°C. Strains were recovered in tryptic soy broth (TSB) or tryptic soy agar (TSA) at 37°C (Becton, Dickinson, Franklin Lake, NJ). We obtained nine food-derived isolates representing six of the seven outbreak strains from the USDA Food Safety and Inspection Service through the USDA Agricultural Research Service (28). These isolates were compared to Salmonella Heidelberg reference strain SL476, Salmonella Typhimurium reference strain ATCC 14028, and a Salmonella Enteritidis strain.

Genomic analyses.The U.S. Food and Drug Administration sequenced seven of the nine outbreak isolates that we studied and placed their raw read genomic sequences in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) collection (Table 1). Reads were downloaded and transferred to Purdue’s Research Computing clusters. Reads were trimmed using the Trimmomatic trimmer tool to a quality score (Q) of 20 (73) and assembled using the SPAdes (v.3.10.1) algorithm (74), and the assemblies were assessed using the Quast (v.3.2) program (75). Assemblies were submitted to Center for Genomic Epidemiology (CGE) servers to search for plasmids (PlasmidFinder [v.2.0] [76]) and antibiotic resistance genes and mutations (ResFinder [v.3.1] [77]) and to determine the MLST of the isolates (MLST [v.2.0] [32]). All CGE programs were run on the default settings. As an additional typing mechanism, raw reads were submitted to analysis with the CGE cgMLSTFInder tool (v.1.0) (36). To further assess isolate relationships, the (i) original 8 isolates and (ii) a fuller set of 12 isolates (the same 8 isolates plus 4 isolates from 2013 with PFGE profiles matching those of the two other 2013-2014 outbreak strains [JF6X01.0022, JF6X01.0122]) were downloaded from the NCBI SRA and submitted to NDtree (v.1.2), a raw reads-based phylogenetic tree tool (78–80). NDtree was allowed to choose the best reference genome and chose that with GenBank accession no. NC_011083 (SL476) as the template. The phylogenetic tree was visualized using the FigTree (v.1.4.4) program (available from http://tree.bio.ed.ac.uk/software/figtree/).

Heat shock assays.S. enterica cultures were isolated on TSA and grown overnight for 18 to 24 h. A single colony was inoculated into 5 ml of TSB and grown for 16 h at 37°C at 200 rpm; overnight cultures were serially diluted to 103 CFU/ml. Cells were grown to log phase (6 h) or stationary phase (12 h), as determined by the use of standard growth curves. Logarithmic- and stationary-phase cultures were aliquoted into 2-ml microcentrifuge tubes (VWR, Radnor, PA) and heat shocked in a 56°C water bath (VWR, Radnor, PA). Aliquots were removed, serially diluted, and plated on tryptic soy agar (Becton, Dickinson, Franklin Lake, NJ) for 0, 2.5, 5, and 7.5 min (log phase) and for 0, 15, 30, 45, and 60 min (stationary phase) of total exposure to 56°C. The time to temperature was determined to be 2.4 min using a thermocouple (Omega, Norwalk, CT) and was accounted for in the total exposure time. The plates were incubated at 37°C for 36 ± 1 h. Experiments were conducted in biological triplicate. Statistically significant differences between isolates and reference strains were determined by repeated-measures analysis of variance (ANOVA) followed by Dunnett’s test, with α equal to 0.05. Analysis was performed in RStudio (v.0.98.1091) software.

Salmonella attachment assays.Attachment to abiotic surfaces as an indicator of biofilm-forming potential was measured using a standard crystal violet assay as previously described (81, 82). Isolated colonies of 10 Salmonella Heidelberg isolates, 4 Salmonella Typhimurium isolates, 4 Salmonella Newport food isolates, Salmonella Typhimurium ATCC 14028, and Salmonella Enteritidis were transferred to TSB and incubated at 37°C for 16 to 24 h to an optical density at 600 nm (to OD600) of 0.6 to 0.8. Polystyrene 96-well plates (Corning, Corning, NY) were prepared with 1× (nutrient-rich conditions) or 1/20× (nutrient-limited conditions) TSB, inoculated, and incubated statically at 4°C or 21°C for 1, 2, or 5 days. Attachment in 1/20× TSB at 21°C at 5, 7, and 9 days was also assessed (see Fig. S2 and S3 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507). Attached cells were stained with 0.4% crystal violet and indirectly quantified by determination of the OD600. Experiments were conducted in biological triplicate. Statistical analyses were performed in SAS (v.9.4) software (SAS Institute, Cary, NC). The OD600 was log10 transformed to fit a generalized linear mixed model with the Proc Glimmix procedure. Analysis of attachment capacity was conducted at 4°C (n = 538), and 21°C (n = 536), respectively, at both the serovar and isolate levels (α = 0.05). The main effects “serovar,” “isolate,” “time,” “serovar × time,” and “isolate × time” were tested to determine a significant correlation with attachment capacity. Tukey pairwise comparisons and least-squares means with the Tukey adjustment were used to identify significant correlations.

PAA tolerance assay.Tolerance to a peracetic acid-based sanitizer (BioSide HS, 15%; Madison Chemical Co., Madison, IN) was determined for 18 isolates (Table 2) using a microtiter-based assay as described by Wang et al. (81). Briefly, isolated colonies were grown overnight on TSA at 37°C and then suspended in TSB to an OD600 of 0.6 to 0.7. A 1:100 dilution of the suspension was transferred to 2× TSB for a final concentration of 1× TSB; 50 μl of the diluted suspension was transferred to assay plates containing 50 μl of each sanitizer solution for final concentrations of 25, 50, 100, and 250 ppm sanitizer. The plates were incubated statically at 37°C. The plates were scanned at OD600 at 0 and 24 h to determine the MIC. Experiments were conducted in biological triplicate. Isolates that consistently grew in ≥50 ppm PAA after 24 h (OD600 > 0.1) were considered tolerant (see Table S3 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507). Analysis is limited to descriptive statistics.

Caco-2 cell invasion assays.Caco-2 cell (gentamicin protection) invasions assays were conducted as previously described (83). Salmonella strains were grown in TSB to late log phase (OD600 = 0.4); each inoculum was quantified by serial dilution and plating on TSA. Experiments were conducted in biological triplicate. Data were analyzed by ANOVA with Tukey’s test for strain mean comparisons in SAS (v.9.4) software; significance was defined as a Padj value of <0.05.

RNA sequencing and data processing.RNA was extracted from stationary-phase (37°C; control) and heat-shocked (5.1 min of total exposure at 56°C) R1-0006, R1-0007, and SL476 cells using the TRIzol reagent (Ambion, Foster City, CA) in triplicate for each isolate and condition. A Ribo-Zero kit (Illumina, San Diego, CA) was used to remove rRNA; TURBO DNase (Thermo Fisher, Waltham, MA) was used to remove DNA. DNA removal was assessed using quantitative reverse transcription-PCR for rpoD and gyrB (71, 84–86). Primers, probes, and 2× master mix were purchased from Integrated DNA Technologies [IDT], Coralville, IA) (see Table S4 posted at https://www.researchgate.net/publication/333609668_AJE2_Supplemental_data_20190507). RNA samples with no reverse transcriptase were used to assess residual DNA contamination. cDNA libraries were constructed using a ScriptSeq (v.2) library kit for bacteria using the low-input protocol with the RiboZero kit (Illumina, San Diego, CA), and samples were sequenced on an Illumina HiSeq 2500 platform (San Diego, CA) by the Purdue University Genomics Core (West Lafayette, IN).

Transcriptomic analyses were performed on Purdue University Research Computing clusters. Reads were trimmed to Q20, and adaptors were removed using the Trimmomatic (v.0.36) trimmer tool (73). The trimmed reads were aligned with the SL476 index sequence (created from NCBI RefSeq genome accession no. NC_011083.1 using the Bowtie2 [v.2.1.0] program) using TopHat (v.2.1.1) software (87); read counts were obtained using the HTSeq (v.0.6.1) package (88). Differential expression between control and heat shock samples was determined using the DESeq2 (v.1.12.4) program (89). For all comparisons, genes with a base mean read count of <30 were removed due to low read counts. Genes with a log2 fold change in expression of >1 and a Padj value of <0.05 were considered significantly differentially expressed (DE). Transcriptomes were annotated using the KEGG online database for SL476 as well as a publicly available general feature format (GFF) of SL476. Unannotated gene calls were made using the National Center for Biotechnology Information (NCBI) online protein BLAST (blastp suite) server (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastp&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome) to search the nonredundant protein sequences database. For proteins with ambiguous results, protein sequences (per KEGG SL476 data) of top gene hits were aligned with sequences in the UniProt database (90) using the EMBOSS Needle gene-wise aligner (protein setting, http://www.ebi.ac.uk/Tools/psa/emboss_needle/) (91). Where genes had multiple names, both were included with the current or most common name first.

Data availability.RNA sequencing data (raw read files) for this project are available for download from NCBI GEO under accession number GSE103418 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE103418).

ACKNOWLEDGMENTS

We thank and acknowledge John Johnston, USDA-FSIS scientific liaison, and Todd Ward and James Sweezey of USDA ARS NCAUR for facilitating and providing the Salmonella Heidelberg outbreak-associated isolates. We thank Alice Green and Stephanie Defibaugh-Chavez of the USDA-FSIS and PulseNet participating laboratories for generously sharing epidemiological data. We also thank Peter Pascuzzi and Michael Gribskov for RNA sequencing analysis discussions and the staff of the Purdue Rosen Center for Advanced Computing for training in command-line next-generation sequencing analyses in Unix. Clara Assisi and Susan Hammons assisted with heat shock assays and advised on the RNA extraction process, and Caitlinn Lineback assisted with invasion assays.

The heat shock assays and RNA sequencing were funded by the Purdue Center for Food Safety Engineering through the USDA Agricultural Research Service (award number 59-1935-2-279). A. J. Etter and all other experiments were supported by the USDA National Institute of Food and Agriculture Predoctoral Fellowships program (award number 2017-67011-26041). H. F. Oliver is supported by USDA National Institute of Food and Agriculture Hatch project 2016-67017-24459. This material is based upon work supported by the U.S. Department of Agriculture, Agricultural Research Service, under agreement no. 59-8072-6-001.

Any opinions, findings, conclusion, or recommendations expressed in this publication are those of the authors and do not necessarily reflect the view of the U.S. Department of Agriculture.

FOOTNOTES

    • Received 8 May 2019.
    • Accepted 3 June 2019.
    • Accepted manuscript posted online 7 June 2019.
  • Copyright © 2019 American Society for Microbiology.

All Rights Reserved.

REFERENCES

  1. 1.↵
    CDC. 2012. Multistate outbreak of human Salmonella Heidelberg infections linked to “kosher broiled chicken livers” from Schreiber Processing Corporation. http://www.cdc.gov/salmonella/heidelberg-chickenlivers/011112/index.html. Accessed 27 November 2013.
  2. 2.↵
    CDC. 2014. Multistate outbreak of multidrug-resistant Salmonella Heidelberg infections linked to Foster Farms brand chicken. http://www.cdc.gov/salmonella/heidelberg-10-13/index.html. Accessed 11 November 2018.
  3. 3.↵
    CDC. 2011. Multistate outbreak of human Salmonella Heidelberg infections linked to ground turkey. http://www.cdc.gov/salmonella/heidelberg/111011/index.html. Accessed 27 November 2013.
  4. 4.↵
    CDC. 2018. Multistate outbreak of multidrug-resistant Salmonella Heidelberg infections linked to contact with dairy calves (final update). https://www.cdc.gov/salmonella/heidelberg-11-16/index.html. Accessed 21 September 2018.
  5. 5.↵
    CDC. 2015. Multistate outbreak of multidrug-resistant Salmonella I 4,[5],12:i:− and Salmonella Infantis infections linked to pork (final update). http://www.cdc.gov/salmonella/pork-08-15/index.html. Accessed 20 April 2019.
  6. 6.↵
    CDC. 2015. Multistate outbreak of Salmonella Poona infections linked to imported cucumbers (final update). http://www.cdc.gov/salmonella/poona-09-15/index.html. Accessed 20 April 2019.
  7. 7.↵
    CDC. 2012. Multistate outbreak of Salmonella Typhimurium and Salmonella Newport infections linked to cantaloupe (final update). https://www.cdc.gov/salmonella/typhimurium-cantaloupe-08-12/. Accessed 29 December 2018.
  8. 8.↵
    1. Den Bakker HC,
    2. Switt AIM,
    3. Cummings CA,
    4. Hoelzer K,
    5. Degoricija L,
    6. Rodriguez-Rivera LD,
    7. Wright EM,
    8. Fang R,
    9. Davis M,
    10. Root T,
    11. Schoonmaker-Bopp D,
    12. Musser KA,
    13. Villamil E,
    14. Waechter HN,
    15. Kornstein L,
    16. Furtado MR,
    17. Wiedmann M
    . 2011. 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 77:8648–8655. doi:10.1128/AEM.06538-11.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    CDC. 2009. Multistate outbreak of Salmonella Typhimurium infections linked to peanut butter, 2008-2009 (final update). https://www.cdc.gov/salmonella/2009/peanut-butter-2008-2009.html.
  10. 10.↵
    1. Jackson BR,
    2. Griffin PM,
    3. Cole D,
    4. Walsh KA,
    5. Chai SJ
    . 2013. Outbreak-associated Salmonella enterica serotypes and food commodities, United States, 1998-2008. Emerg Infect Dis 19:1239–1244. doi:10.3201/eid1908.121511.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Zhao S,
    2. White DG,
    3. Friedman SL,
    4. Glenn A,
    5. Blickenstaff K,
    6. Ayers SL,
    7. Abbott JW,
    8. Hall-Robinson E,
    9. McDermott PF
    . 2008. Antimicrobial resistance in Salmonella enterica serovar Heidelberg isolates from retail meats, including poultry, from 2002 to 2006. Appl Environ Microbiol 74:6656–6662. doi:10.1128/AEM.01249-08.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    Food Safety and Inspection Service, USDA. 2016. New performance standards for Salmonella and Campylobacter in not-ready-to-eat comminuted chicken and turkey products and raw chicken parts and changes to related agency verification procedures: response to comments and announcement of implementation schedule. Docket no. FSIS-2014-0023. Fed Regist 81:7285–7300.
    OpenUrl
  13. 13.↵
    1. Bauermeister LJ,
    2. Bowers JW,
    3. Townsend JC,
    4. McKee SR
    . 2008. The microbial and quality properties of poultry carcasses treated with peracetic acid as an antimicrobial treatment. Poult Sci 87:2390–2398. doi:10.3382/ps.2008-00087.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Rothrock MJ, Jr,
    2. Ingram KD,
    3. Gamble J,
    4. Guard J,
    5. Cicconi-Hogan KM,
    6. Hinton A, Jr,
    7. Hiett KL
    . 2015. The characterization of Salmonella enterica serotypes isolated from the scalder tank water of a commercial poultry processing plant: recovery of a multidrug-resistant Heidelberg strain. Poult Sci 94:467–472. doi:10.3382/ps/peu060.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Lianou A,
    2. Koutsoumanis KP
    . 2013. Evaluation of the strain variability of Salmonella enterica acid and heat resistance. Food Microbiol 34:259–267. doi:10.1016/j.fm.2012.10.009.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Stepanović S,
    2. Cirković I,
    3. Ranin L,
    4. Svabić-Vlahović M
    . 2004. Biofilm formation by Salmonella spp. and Listeria monocytogenes on plastic surface. Lett Appl Microbiol 38:428–432. doi:10.1111/j.1472-765X.2004.01513.x.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    1. Solomon EB,
    2. Niemira BA,
    3. Sapers GM,
    4. Annous BA
    . 2005. Biofilm formation, cellulose production, and curli biosynthesis by Salmonella originating from produce, animal, and clinical sources. J Food Prot 68:906–912. doi:10.4315/0362-028X-68.5.906.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    1. Vestby LK,
    2. Moretro T,
    3. Langsrud S,
    4. Heir E,
    5. Nesse LL
    . 2009. Biofilm forming abilities of Salmonella are correlated with persistence in fish meal- and feed factories. BMC Vet Res 5:20. doi:10.1186/1746-6148-5-20.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Rasimus S,
    2. Kolari M,
    3. Rita H,
    4. Hoornstra D,
    5. Salkinoja-Salonen M
    . 2011. Biofilm-forming bacteria with varying tolerance to peracetic acid from a paper machine. J Ind Microbiol Biotechnol 38:1379–1390. doi:10.1007/s10295-010-0921-4.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Kim SH,
    2. Wei CI
    . 2009. Molecular characterization of biofilm formation and attachment of Salmonella enterica serovar Typhimurium DT104 on food contact surfaces. J Food Prot 72:1841–1847. doi:10.4315/0362-028X-72.9.1841.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Romling U,
    2. Bokranz W,
    3. Rabsch W,
    4. Zogaj X,
    5. Nimtz M,
    6. Tschape H
    . 2003. Occurrence and regulation of the multicellular morphotype in Salmonella serovars important in human disease. Int J Med Microbiol 293:273–285. doi:10.1078/1438-4221-00268.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    1. Steenackers H,
    2. Hermans K,
    3. Vanderleyden J,
    4. De Keersmaecker S
    . 2012. Salmonella biofilms: an overview on occurrence, structure, regulation and eradication. Food Res Int 45:502–531. doi:10.1016/j.foodres.2011.01.038.
    OpenUrlCrossRef
  23. 23.↵
    1. Seixas R,
    2. Gabriel M,
    3. Machado J,
    4. Tavares L,
    5. Bernardo F,
    6. Oliveira M
    . 2014. Effect of simulated gastrointestinal conditions on biofilm formation by Salmonella 1,4,[5],12:i:-. Scientific World Journal 2014:153956. doi:10.1155/2014/153956.
    OpenUrlCrossRef
  24. 24.↵
    1. Argarwal RK,
    2. Singh S,
    3. Bhilegaonkar KN,
    4. Singh VP
    . 2011. Optimization of microtiter plate assay for the testing of biofilm formation ability in different Salmonella serotypes. Int Food Res J 18:1493–1498.
    OpenUrl
  25. 25.↵
    1. Abbot JM,
    2. Byrd-Bredbenner C,
    3. Schaffner D,
    4. Bruhn CM,
    5. Blalock L
    . 2009. Comparison of food safety cognitions and self-reported food-handling behaviors with observed food safety behaviors of young adults. Eur J Clin Nutr 63:572–579. doi:10.1038/sj.ejcn.1602961.
    OpenUrlCrossRefPubMed
  26. 26.↵
    FDA. 2010. 2010 food safety survey. https://www.fda.gov/downloads/Food/FoodScienceResearch/ConsumerBehaviorResearch/UCM407008.pdf.
  27. 27.↵
    FDA. 2016. 2016 FDA food safety survey. https://www.fda.gov/downloads/Food/FoodScienceResearch/ConsumerBehaviorResearch/UCM529453.pdf.
  28. 28.↵
    1. Gieraltowski L,
    2. Higa J,
    3. Peralta V,
    4. Green A,
    5. Schwensohn C,
    6. Rosen H,
    7. Libby T,
    8. Kissler B,
    9. Marsden-Haug N,
    10. Booth H,
    11. Kimura A,
    12. Grass J,
    13. Bicknese A,
    14. Tolar B,
    15. Defibaugh-Chávez S,
    16. Williams I,
    17. Wise M
    . 2016. National outbreak of multidrug resistant Salmonella Heidelberg infections linked to a single poultry company. PLoS One 11:e0162369. doi:10.1371/journal.pone.0162369.
    OpenUrlCrossRef
  29. 29.↵
    1. Green A,
    2. Defibaugh-Chavez S,
    3. Douris A,
    4. Vetter D,
    5. Atkinson R,
    6. Kissler B,
    7. Khroustalev A,
    8. Robertson K,
    9. Sharma Y,
    10. Becker K,
    11. Dessai U,
    12. Antoine N,
    13. Allen L,
    14. Holt K,
    15. Gieraltowski L,
    16. Wise M,
    17. Schwensohn C
    . 2018. Intensified sampling in response to a Salmonella Heidelberg outbreak associated with multiple establishments within a single poultry corporation. Foodborne Pathog Dis 15:153–160. doi:10.1089/fpd.2017.2340.
    OpenUrlCrossRef
  30. 30.↵
    1. Jones TF,
    2. Ingram LA,
    3. Cieslak PR,
    4. Vugia DJ,
    5. Tobin-D'Angelo M,
    6. Hurd S,
    7. Medus C,
    8. Cronquist A,
    9. Angulo FJ
    . 2008. Salmonellosis outcomes differ substantially by serotype. J Infect Dis 198:109–114. doi:10.1086/588823.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    1. McKee SR,
    2. Townsend JC,
    3. Bilgili SF
    . 2008. Use of a scald additive to reduce levels of Salmonella Typhimurium during poultry processing. Poult Sci 87:1672–1677. doi:10.3382/ps.2008-00061.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Larsen MV,
    2. Cosentino S,
    3. Rasmussen S,
    4. Friis C,
    5. Hasman H,
    6. Marvig RL,
    7. Jelsbak L,
    8. Sicheritz-Ponten T,
    9. Ussery DW,
    10. Aarestrup FM,
    11. Lund O
    . 2012. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol 50:1355–1361. doi:10.1128/JCM.06094-11.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Michaux C,
    2. Holmqvist E,
    3. Vasicek E,
    4. Sharan M,
    5. Barquist L,
    6. Westermann AJ,
    7. Gunn JS,
    8. Vogel J
    . 2017. RNA target profiles direct the discovery of virulence functions for the cold-shock proteins CspC and CspE. Proc Natl Acad Sci U S A 114:6824–6829. doi:10.1073/pnas.1620772114.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Brodsky IE,
    2. Ernst RK,
    3. Miller SI,
    4. Falkow S
    . 2002. mig-14 is a Salmonella gene that plays a role in bacterial resistance to antimicrobial peptides. J Bacteriol 184:3203–3213. doi:10.1128/jb.184.12.3203-3213.2002.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Roland KL,
    2. Esther CR,
    3. Spitznagel JK
    . 1994. Isolation and characterization of a gene, pmrD, from Salmonella typhimurium that confers resistance to polymyxin when expressed in multiple copies. J Bacteriol 176:3589–3597. doi:10.1128/jb.176.12.3589-3597.1994.
    OpenUrlAbstract/FREE Full Text
  36. 36.↵
    1. Alikhan NF,
    2. Zhou Z,
    3. Sergeant MJ,
    4. Achtman M
    . 2018. A genomic overview of the population structure of Salmonella. PLoS Genet 14:e1007261. doi:10.1371/journal.pgen.1007261.
    OpenUrlCrossRef
  37. 37.↵
    1. Hoffmann M,
    2. Zhao S,
    3. Pettengill J,
    4. Luo Y,
    5. Monday SR,
    6. Abbott J,
    7. Ayers SL,
    8. Cinar HN,
    9. Muruvanda T,
    10. Li C,
    11. Allard MW,
    12. Whichard J,
    13. Meng J,
    14. Brown EW,
    15. McDermott PF
    . 2014. Comparative genomic analysis and virulence differences in closely related Salmonella enterica serotype Heidelberg isolates from humans, retail meats, and animals. Genome Biol Evol 6:1046–1068. doi:10.1093/gbe/evu079.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Hoffmann M,
    2. Luo Y,
    3. Lafon PC,
    4. Timme R,
    5. Allard MW,
    6. McDermott PF,
    7. Brown EW,
    8. Zhao S
    . 2013. Genome sequences of Salmonella enterica serovar Heidelberg isolates isolated in the United States from a multistate outbreak of human Salmonella infections. Genome Announc 1:e00004-12. doi:10.1128/genomeA.00004-12.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    1. Keefer AB,
    2. Xiaoli L,
    3. M’ikanatha NM,
    4. Yao K,
    5. Hoffmann M,
    6. Dudley EG
    . 2019. Retrospective whole-genome sequencing analysis distinguished PFGE and drug-resistance-matched retail meat and clinical Salmonella isolates. Microbiology 165:270–286. doi:10.1099/mic.0.000768.
    OpenUrlCrossRef
  40. 40.↵
    1. Salipante SJ,
    2. Hall BG
    . 2003. Determining the limits of the evolutionary potential of an antibiotic resistance gene. Mol Biol Evol 20:653–659. doi:10.1093/molbev/msg074.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    1. Magnet S,
    2. Courvalin P,
    3. Lambert T
    . 1999. Activation of the cryptic aac(6′)-Iy aminoglycoside resistance gene of Salmonella by a chromosomal deletion generating a transcriptional fusion. J Bacteriol 181:6650–6655.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Neuert S,
    2. Nair S,
    3. Day MR,
    4. Doumith M,
    5. Ashton PM,
    6. Mellor KC,
    7. Jenkins C,
    8. Hopkins KL,
    9. Woodford N,
    10. de Pinna E,
    11. Godbole G,
    12. Dallman TJ
    . 2018. Prediction of phenotypic antimicrobial resistance profiles from whole genome sequences of non-typhoidal Salmonella enterica. Front Microbiol 9:592. doi:10.3389/fmicb.2018.00592.
    OpenUrlCrossRef
  43. 43.↵
    1. Giraud E,
    2. Brisabois A,
    3. Martel JL,
    4. Chaslus-Dancla E
    . 1999. Comparative studies of mutations in animal isolates and experimental in vitro- and in vivo-selected mutants of Salmonella spp. suggest a counterselection of highly fluoroquinolone-resistant strains in the field. Antimicrob Agents Chemother 43:2131–2137. doi:10.1128/AAC.43.9.2131.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Walker RA,
    2. Saunders N,
    3. Lawson AJ,
    4. Lindsay EA,
    5. Dassama M,
    6. Ward LR,
    7. Woodward MJ,
    8. Davies RH,
    9. Liebana E,
    10. Threlfall EJ
    . 2001. Use of a LightCycler gyrA mutation assay for rapid identification of mutations conferring decreased susceptibility to ciprofloxacin in multiresistant Salmonella enterica serotype Typhimurium DT104 isolates. J Clin Microbiol 39:1443–1448. doi:10.1128/JCM.39.4.1443-1448.2001.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    1. Hopkins KL,
    2. Arnold C,
    3. Threlfall EJ
    . 2007. Rapid detection of gyrA and parC mutations in quinolone-resistant Salmonella enterica using pyrosequencing technology. J Microbiol Methods 68:163–171. doi:10.1016/j.mimet.2006.07.006.
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    1. Ling JM,
    2. Chan EW,
    3. Lam AW,
    4. Cheng AF
    . 2003. Mutations in topoisomerase genes of fluoroquinolone-resistant salmonellae in Hong Kong. Antimicrob Agents Chemother 47:3567–3573. doi:10.1128/AAC.47.11.3567-3573.2003.
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    1. Humpheson L,
    2. Adams MR,
    3. Anderson WA,
    4. Cole MB
    . 1998. Biphasic thermal inactivation kinetics in Salmonella Enteritidis PT4. Appl Environ Microbiol 64:459–464.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Doyle ME,
    2. Mazzotta AS
    . 2000. Review of studies on the thermal resistance of Salmonellae. J Food Prot 63:779–795. doi:10.4315/0362-028X-63.6.779.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Pin C,
    2. Hansen T,
    3. Muñoz-Cuevas M,
    4. de Jonge R,
    5. Rosenkrantz JT,
    6. Löfström C,
    7. Aarts H,
    8. Olsen JE
    . 2012. The transcriptional heat shock response of Salmonella Typhimurium shows hysteresis and heated cells show increased resistance to heat and acid stress. PLoS One 7:e51196. doi:10.1371/journal.pone.0051196.
    OpenUrlCrossRef
  50. 50.↵
    1. Sirsat SA,
    2. Baker CA,
    3. Park SH,
    4. Muthaiyan A,
    5. Dowd SE,
    6. Ricke SC
    . 2015. Transcriptomic response of Salmonella Typhimurium heat shock gene expression under thermal stress at 48°C. J Food Res 4:51–56. doi:10.5539/jfr.v4n5p51.
    OpenUrlCrossRef
  51. 51.↵
    1. Milillo SR,
    2. Martin E,
    3. Muthaiyan A,
    4. Ricke SC
    . 2011. Immediate reduction of Salmonella enterica serotype Typhimurium viability via membrane destabilization following exposure to multiple-hurdle treatments with heated, acidified organic acid salt solutions. Appl Environ Microbiol 77:3765–3772. doi:10.1128/AEM.02839-10.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. Bang IS,
    2. Kim BH,
    3. Foster JW,
    4. Park YK
    . 2000. OmpR regulates the stationary-phase acid tolerance response of Salmonella enterica serovar Typhimurium. J Bacteriol 182:2245–2252. doi:10.1128/jb.182.8.2245-2252.2000.
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    1. Rychlik I,
    2. Barrow PA
    . 2005. Salmonella stress management and its relevance to behaviour during intestinal colonisation and infection. FEMS Microbiol Rev 29:1021–1040. doi:10.1016/j.femsre.2005.03.005.
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    1. Nguyen SV,
    2. Harhay GP,
    3. Bono JL,
    4. Smith TP,
    5. Harhay DM
    . 2017. Genome sequence of the thermotolerant foodborne pathogen Salmonella enterica serovar Senftenberg ATCC 43845 and phylogenetic analysis of loci encoding increased protein quality control mechanisms. mSystems 2:e00190-16. doi:10.1128/mSystems.00190-16.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Keelara S,
    2. Thakur S,
    3. Patel J
    . 2016. Biofilm formation by environmental isolates of Salmonella and their sensitivity to natural antimicrobials. Foodborne Pathog Dis 13:509–516. doi:10.1089/fpd.2016.2145.
    OpenUrlCrossRef
  56. 56.↵
    1. de Ornellas Dutka Garcia KC,
    2. de Oliveira Corrêa IM,
    3. Pereira LP,
    4. Silva TM,
    5. de Souza Ribeiro Mioni M,
    6. de Moraes Izidoro AC,
    7. Bastos IHV,
    8. Gonçalves GAM,
    9. Okamoto AS,
    10. Andreatti Filho RL
    . 2017. Bacteriophage use to control Salmonella biofilm on surfaces present in chicken slaughterhouses. Poult Sci 96:3392–3398. doi:10.3382/ps/pex124.
    OpenUrlCrossRef
  57. 57.↵
    1. Kim SH,
    2. Wei CI
    . 2007. Biofilm formation by multidrug-resistant Salmonella enterica serotype Typhimurium phage type DT104 and other pathogens. J Food Prot 70:22–29. doi:10.4315/0362-028X-70.1.22.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Dubois-Brissonnet F,
    2. Naitali M,
    3. Mafu AA,
    4. Briandet R
    . 2011. Induction of fatty acid composition modifications and tolerance to biocides in Salmonella enterica serovar Typhimurium by plant-derived terpenes. Appl Environ Microbiol 77:906–910. doi:10.1128/AEM.01480-10.
    OpenUrlAbstract/FREE Full Text
  59. 59.↵
    1. Moken MC,
    2. McMurry LM,
    3. Levy SB
    . 1997. Selection of multiple-antibiotic-resistant (mar) mutants of Escherichia coli by using the disinfectant pine oil: roles of the mar and acrAB loci. Antimicrob Agents Chemother 41:2770–2772. doi:10.1128/AAC.41.12.2770.
    OpenUrlAbstract/FREE Full Text
  60. 60.↵
    1. Lavilla Lerma L,
    2. Benomar N,
    3. Casado Muñoz MDC,
    4. Gálvez A,
    5. Abriouel H
    . 2015. Correlation between antibiotic and biocide resistance in mesophilic and psychrotrophic Pseudomonas spp. isolated from slaughterhouse surfaces throughout meat chain production. Food Microbiol 51:33–44. doi:10.1016/j.fm.2015.04.010.
    OpenUrlCrossRef
  61. 61.↵
    1. Nagel GM,
    2. Bauermeister LJ,
    3. Bratcher CL,
    4. Singh M,
    5. McKee SR
    . 2013. Salmonella and Campylobacter reduction and quality characteristics of poultry carcasses treated with various antimicrobials in a post-chill immersion tank. Int J Food Microbiol 165:281–286. doi:10.1016/j.ijfoodmicro.2013.05.016.
    OpenUrlCrossRefPubMed
  62. 62.↵
    1. Humphrey TJ,
    2. Slater E,
    3. McAlpine K,
    4. Rowbury RJ,
    5. Gilbert RJ
    . 1995. Salmonella Enteritidis phage type 4 isolates more tolerant of heat, acid, or hydrogen peroxide also survive longer on surfaces. Appl Environ Microbiol 61:3161–3164.
    OpenUrlAbstract/FREE Full Text
  63. 63.↵
    1. Humphrey TJ,
    2. Richardson NP,
    3. Statton KM,
    4. Rowbury RJ
    . 1993. Effects of temperature shift on acid and heat tolerance in Salmonella Enteritidis phage type 4. Appl Environ Microbiol 59:3120–3122.
    OpenUrlAbstract/FREE Full Text
  64. 64.↵
    1. Hamilton S,
    2. Bongaerts RJ,
    3. Mulholland F,
    4. Cochrane B,
    5. Porter J,
    6. Lucchini S,
    7. Lappin-Scott HM,
    8. Hinton JC
    . 2009. The transcriptional programme of Salmonella enterica serovar Typhimurium reveals a key role for tryptophan metabolism in biofilms. BMC Genomics 10:599. doi:10.1186/1471-2164-10-599.
    OpenUrlCrossRefPubMed
  65. 65.↵
    1. Nishino K,
    2. Latifi T,
    3. Groisman EA
    . 2006. Virulence and drug resistance roles of multidrug efflux systems of Salmonella enterica serovar Typhimurium. Mol Microbiol 59:126–141. doi:10.1111/j.1365-2958.2005.04940.x.
    OpenUrlCrossRefPubMedWeb of Science
  66. 66.↵
    1. Baucheron S,
    2. Nishino K,
    3. Monchaux I,
    4. Canepa S,
    5. Maurel MC,
    6. Coste F,
    7. Roussel A,
    8. Cloeckaert A,
    9. Giraud E
    . 2014. Bile-mediated activation of the acrAB and tolC multidrug efflux genes occurs mainly through transcriptional derepression of ramA in Salmonella enterica serovar Typhimurium. J Antimicrob Chemother 69:2400–2406. doi:10.1093/jac/dku140.
    OpenUrlCrossRefPubMed
  67. 67.↵
    1. Nikaido E,
    2. Shirosaka I,
    3. Yamaguchi A,
    4. Nishino K
    . 2011. Regulation of the acrAB multidrug efflux pump in Salmonella enterica serovar Typhimurium in response to indole and paraquat. Microbiology 157:648–655. doi:10.1099/mic.0.045757-0.
    OpenUrlCrossRefPubMedWeb of Science
  68. 68.↵
    1. Carlson SA,
    2. Browning M,
    3. Ferris KE,
    4. Jones BD
    . 2000. Identification of diminished tissue culture invasiveness among multiple antibiotic resistant Salmonella Typhimurium DT104. Microb Pathog 28:37–44. doi:10.1006/mpat.1999.0322.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    1. Carlson SA,
    2. Willson RM,
    3. Crane AJ,
    4. Ferris KE
    . 2000. Evaluation of invasion-conferring genotypes and antibiotic-induced hyperinvasive phenotypes in multiple antibiotic resistant Salmonella Typhimurium DT104. Microb Pathog 28:373–378. doi:10.1006/mpat.2000.0355.
    OpenUrlCrossRefPubMedWeb of Science
  70. 70.↵
    1. Varma JK,
    2. Greene KD,
    3. Ovitt J,
    4. Barrett TJ,
    5. Medalla F,
    6. Angulo FJ
    . 2005. Hospitalization and antimicrobial resistance in Salmonella outbreaks, 1984–2002. Emerg Infect Dis 11:943–946. doi:10.3201/eid1106.041231.
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    1. Fink RC,
    2. Evans MR,
    3. Porwollik S,
    4. Vazquez-Torres A,
    5. Jones-Carson J,
    6. Troxell B,
    7. Libby SJ,
    8. McClelland M,
    9. Hassan HM
    . 2007. FNR is a global regulator of virulence and anaerobic metabolism in Salmonella enterica serovar Typhimurium (ATCC 14028s). J Bacteriol 189:2262–2273. doi:10.1128/JB.00726-06.
    OpenUrlAbstract/FREE Full Text
  72. 72.↵
    PulseNet. 2017. Blood and stool isolates by strain for the 2013-2014 Foster Farms outbreak investigation. Centers for Disease Control and Prevention, Atlanta, GA.
  73. 73.↵
    1. Bolger AM,
    2. Lohse M,
    3. Usadel B
    . 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120. doi:10.1093/bioinformatics/btu170.
    OpenUrlCrossRefPubMedWeb of Science
  74. 74.↵
    1. Bankevich A,
    2. Nurk S,
    3. Antipov D,
    4. Gurevich AA,
    5. Dvorkin M,
    6. Kulikov AS,
    7. Lesin VM,
    8. Nikolenko SI,
    9. Pham S,
    10. Prjibelski AD,
    11. Pyshkin AV,
    12. Sirotkin AV,
    13. Vyahhi N,
    14. Tesler G,
    15. Alekseyev MA,
    16. Pevzner PA
    . 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19:455–477. doi:10.1089/cmb.2012.0021.
    OpenUrlCrossRefPubMed
  75. 75.↵
    1. Gurevich A,
    2. Saveliev V,
    3. Vyahhi N,
    4. Tesler G
    . 2013. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29:1072–1075. doi:10.1093/bioinformatics/btt086.
    OpenUrlCrossRefPubMedWeb of Science
  76. 76.↵
    1. Carattoli A,
    2. Zankari E,
    3. García-Fernández A,
    4. Voldby Larsen M,
    5. Lund O,
    6. Villa L,
    7. Møller Aarestrup F,
    8. Hasman H
    . 2014. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother 58:3895–3903. doi:10.1128/AAC.02412-14.
    OpenUrlAbstract/FREE Full Text
  77. 77.↵
    1. Zankari E,
    2. Hasman H,
    3. Cosentino S,
    4. Vestergaard M,
    5. Rasmussen S,
    6. Lund O,
    7. Aarestrup FM,
    8. Larsen MV
    . 2012. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother 67:2640–2644. doi:10.1093/jac/dks261.
    OpenUrlCrossRefPubMedWeb of Science
  78. 78.↵
    1. Kaas RS,
    2. Leekitcharoenphon P,
    3. Aarestrup FM,
    4. Lund O
    . 2014. Solving the problem of comparing whole bacterial genomes across different sequencing platforms. PLoS One 9:e104984. doi:10.1371/journal.pone.0104984.
    OpenUrlCrossRef
  79. 79.↵
    1. Leekitcharoenphon P,
    2. Nielsen EM,
    3. Kaas RS,
    4. Lund O,
    5. Aarestrup FM
    . 2014. Evaluation of whole genome sequencing for outbreak detection of Salmonella enterica. PLoS One 9:e87991. doi:10.1371/journal.pone.0087991.
    OpenUrlCrossRef
  80. 80.↵
    1. Joensen KG,
    2. Scheutz F,
    3. Lund O,
    4. Hasman H,
    5. Kaas RS,
    6. Nielsen EM,
    7. Aarestrup FM
    . 2014. Real-time whole-genome sequencing for routine typing, surveillance, and outbreak detection of verotoxigenic Escherichia coli. J Clin Microbiol 52:1501–1510. doi:10.1128/JCM.03617-13.
    OpenUrlAbstract/FREE Full Text
  81. 81.↵
    1. Wang J,
    2. Ray AJ,
    3. Hammons SR,
    4. Oliver HF
    . 2015. Persistent and transient Listeria monocytogenes strains from retail deli environments vary in their ability to adhere and form biofilms and rarely have inlA premature stop codons. Foodborne Pathog Dis 12:151–158. doi:10.1089/fpd.2014.1837.
    OpenUrlCrossRefPubMed
  82. 82.↵
    1. Karsi A,
    2. Howe K,
    3. Kirkpatrick TB,
    4. Wills R,
    5. Bailey RH,
    6. Lawrence ML
    . 2008. Development of bioluminescent Salmonella strains for use in food safety. BMC Microbiol 8:10. doi:10.1186/1471-2180-8-10.
    OpenUrlCrossRefPubMed
  83. 83.↵
    1. Nightingale KK,
    2. Windham K,
    3. Martin KE,
    4. Yeung M,
    5. Wiedmann M
    . 2005. Select Listeria monocytogenes subtypes commonly found in foods carry distinct nonsense mutations in inlA, leading to expression of truncated and secreted internalin A, and are associated with a reduced invasion phenotype for human intestinal epithelial cells. Appl Environ Microbiol 71:8764–8772. doi:10.1128/AEM.71.12.8764-8772.2005.
    OpenUrlAbstract/FREE Full Text
  84. 84.↵
    1. O'Leary D,
    2. McCabe EM,
    3. McCusker MP,
    4. Martins M,
    5. Fanning S,
    6. Duffy G
    . 2015. Acid environments affect biofilm formation and gene expression in isolates of Salmonella enterica Typhimurium DT104. Int J Food Microbiol 206:7–16. doi:10.1016/j.ijfoodmicro.2015.03.030.
    OpenUrlCrossRef
  85. 85.↵
    1. Li J,
    2. Overall CC,
    3. Johnson RC,
    4. Jones MB,
    5. McDermott JE,
    6. Heffron F,
    7. Adkins JN,
    8. Cambronne ED
    . 2015. ChIP-Seq analysis of the sigmaE regulon of Salmonella enterica serovar Typhimurium reveals new genes implicated in heat shock and oxidative stress response. PLoS One 10:e0138466. doi:10.1371/journal.pone.0138466.
    OpenUrlCrossRef
  86. 86.↵
    1. Cheng Y,
    2. Pedroso AA,
    3. Porwollik S,
    4. McClelland M,
    5. Lee MD,
    6. Kwan T,
    7. Zamperini K,
    8. Soni V,
    9. Sellers HS,
    10. Russell SM,
    11. Maurer JJ
    . 2015. rpoS-regulated core genes involved in the competitive fitness of Salmonella enterica serovar Kentucky in the intestines of chickens. Appl Environ Microbiol 81:502–514. doi:10.1128/AEM.03219-14.
    OpenUrlAbstract/FREE Full Text
  87. 87.↵
    1. Kim D,
    2. Pertea G,
    3. Trapnell C,
    4. Pimentel H,
    5. Kelley R,
    6. Salzberg SL
    . 2013. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 14:R36. doi:10.1186/gb-2013-14-4-r36.
    OpenUrlCrossRefPubMed
  88. 88.↵
    1. Anders S,
    2. Pyl PT,
    3. Huber W
    . 2015. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics 31:166–169. doi:10.1093/bioinformatics/btu638.
    OpenUrlCrossRefPubMedWeb of Science
  89. 89.↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    . 2014. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15:550. doi:10.1186/s13059-014-0550-8.
    OpenUrlCrossRefPubMed
  90. 90.↵
    The UniProt Consortium. 2017. UniProt: the universal protein knowledgebase. Nucleic Acids Res 45:D158–D169. doi:10.1093/nar/gkw1099.
    OpenUrlCrossRefPubMed
  91. 91.↵
    1. Rice P,
    2. Longden I,
    3. Bleasby A
    . 2000. EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet 16:276–277. doi:10.1016/S0168-9525(00)02024-2.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Salmonella enterica subsp. enterica Serovar Heidelberg Food Isolates Associated with a Salmonellosis Outbreak Have Enhanced Stress Tolerance Capabilities
Andrea J. Etter, Alyssa M. West, John L. Burnett, Sophie Tongyu Wu, Deklin R. Veenhuizen, Raeya A. Ogas, Haley F. Oliver
Applied and Environmental Microbiology Aug 2019, 85 (16) e01065-19; DOI: 10.1128/AEM.01065-19

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Applied and Environmental Microbiology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Salmonella enterica subsp. enterica Serovar Heidelberg Food Isolates Associated with a Salmonellosis Outbreak Have Enhanced Stress Tolerance Capabilities
(Your Name) has forwarded a page to you from Applied and Environmental Microbiology
(Your Name) thought you would be interested in this article in Applied and Environmental Microbiology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Salmonella enterica subsp. enterica Serovar Heidelberg Food Isolates Associated with a Salmonellosis Outbreak Have Enhanced Stress Tolerance Capabilities
Andrea J. Etter, Alyssa M. West, John L. Burnett, Sophie Tongyu Wu, Deklin R. Veenhuizen, Raeya A. Ogas, Haley F. Oliver
Applied and Environmental Microbiology Aug 2019, 85 (16) e01065-19; DOI: 10.1128/AEM.01065-19
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

RNA-Seq
Salmonella enterica serovar Heidelberg
genomics
heat tolerance
outbreak
poultry
stress tolerance

Related Articles

Cited By...

About

  • About AEM
  • Editor in Chief
  • Editorial Board
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Ethics
  • Contact Us

Follow #AppEnvMicro

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

 

Print ISSN: 0099-2240; Online ISSN: 1098-5336