Genomic insights into the pathogenicity and environmental adaptability of Enterococcus hirae R17 isolated from pork offered for retail sale

Abstract Genetic information about Enterococcus hirae is limited, a feature that has compromised our understanding of these clinically challenging bacteria. In this study, comparative analysis was performed of E. hirae R17, a daptomycin‐resistant strain isolated from pork purchased from a retail market in Beijing, China, and three other enterococcal genomes (Enterococcus faecium DO, Enterococcus faecalis V583, and E. hirae ATCC ™9790). Some 1,412 genes were identified that represented the core genome together with an additional 139 genes that were specific to E. hirae R17. The functions of these R17 strain‐specific coding sequences relate to the COGs categories of carbohydrate transport and metabolism and transcription, a finding that suggests the carbohydrate utilization capacity of E. hirae R17 may be more extensive when compared with the other three bacterial species (spp.). Analysis of genomic islands and virulence genes highlighted the potential that horizontal gene transfer played as a contributor of variations in pathogenicity in this isolate. Drug‐resistance gene prediction and antibiotic susceptibility testing indicated E. hirae R17 was resistant to several antimicrobial compounds, including bacitracin, ciprofloxacin, daptomycin, erythromycin, and tetracycline, thereby limiting chemotherapeutic treatment options. Further, tolerance to biocides and metals may confer a phenotype that facilitates the survival and adaptation of this isolate against food preservatives, disinfectants, and antibacterial coatings. The genomic plasticity, mediated by IS elements, transposases, and tandem repeats, identified in the E. hirae R17 genome may support adaptation to new environmental niches, such as those that are found in hospitalized patients. A predicted transmissible plasmid, pRZ1, was found to carry several antimicrobial determinants, along with some predicted pathogenic genes. These data supported the previously determined phenotype confirming that the foodborne E. hirae R17 is a multidrug‐resistant pathogenic bacterium with evident genome plasticity and environmental adaptability.

Enterococci, due to their fermentative potential, were thought to have the capacity when present in food matrices to extend shelflife and to contribute to the improvement of both flavor and texture (Dziewit et al., 2015). However, production of toxic compounds by some of these bacteria raises concerns about the safety of these isolates in food production (Camargo et al., 2014;Choi & Woo, 2013;Franz et al., 2001). Additionally, several studies reported that enterococci can transfer resistance to even more virulent bacteria, such as Staphylococcus aureus (Bellanger et al., 2014;Durand, Brueckner, Sampadian, Willett, & Belliveau, 2014). Hence, the ability of enterococci to act as reservoirs of antibiotic resistance genes within the food chain poses an important food safety risk (Delpech et al., 2012;Gousia, Economou, Bozidis, & Papadopoulou, 2015;Jahan, Krause, & Holley, 2013;Pesavento, Calonico, Ducci, Magnanini, & Lo Nostro, 2014;Stensland et al., 2015). Some Enterococcus spp. are now regarded as opportunistic zoonotic pathogens, arising from their environmental adaptability, contributing toward resistance to unfavorable factors, such as high temperature, high salt concentration, high acid and alkali environment, along with its occurrence and detection in food and foodproducing animals (Anderson et al., 2015;Hammerum, 2012;Kelesidis, 2015;Larsen et al., 2010;Manson, Hancock, & Gilmore, 2010).
The aim of this study was to describe the genome of a daptomycinresistant E. hirae R17 recovered from a food sample, and to investigate the genetic basis underpinning its antimicrobial resistance phenotype, virulence, and environmental adaptability.

| Enterococcus detection and identification
Pork meat, sewage samples, and surface swabs of the ground, walls, and chopping board were collected from a stallholder at a free-trade market in Beijing in 2015 and analyzed microbiologically for the presence of Enterococcus spp., Salmonella spp., Campylobacter spp., and Staphylococcus spp. A quantity of 25 g samples of pork meat was aseptically weighed and separately placed in a stomacher bag (Filtra-Bag, VWR, Luqiao Inc., Beijing, China). Then, 225 ml of brain heart infusion broth (BHI) broth was added and samples were homogenized by stomaching (BagMixer 400, Intersciences Inc., Markham, ON, Canada) for 1 min. The bags were statically incubated for 16 hr at 37°C (Jahan et al., 2013). Broth cultures for presumptive Enterococcus spp. were streaked onto selective isolation indoxylβd-glucoside (mEI) agar (Luqiao Inc., Beijing, China). Plates were then incubated at 42°C for 40-48 hr and then typical Enterococcus spp. colonies were screened and subsequently confirmed by VITEK 2 (bioMérieux, Marcy, l'Etoile, France), Bruker MALDI Biotyper (Germany), and 16S rDNA genebased sequencing as reported previously (Peng, Wang, Hu, & Li, phenotype confirming that the foodborne E. hirae R17 is a multidrug-resistant pathogenic bacterium with evident genome plasticity and environmental adaptability.

| Antimicrobial susceptibility testing
Susceptibility to a panel of antimicrobial agents was determined by broth microdilution, and interpreted according to the criteria based on the Clinical & Laboratory Standards Institute (CLSI) interpretive standards (CLSI, 2014). The minimum inhibitory concentrations (MIC) of 11 antimicrobial compounds were measured for 86 Enterococcus spp. isolates, and these included ampicillin, bacitracin, chloramphenicol, ciprofloxacin, daptomycin, erythromycin, high-level streptomycin, high-level gentamicin, penicillin, tetracycline, and vancomycin. The breakpoint for bacitracin susceptibility was defined as an MIC ≤ 32 μg/ml (Tran, Munita, & Arias, 2015). Enterococcus faecalis ATCC™29212 was used as a control microorganism for these experiments.

| Complete genome sequencing, assembly, and annotation
The complete genome sequence was determined for E. hirae R17, and published previously (Peng et al., 2016a). The genome of E. hirae R17 was sequenced using the Pacific Biosciences RS II sequencing platform (Pacific Biosciences, Menlo Park, CA, USA). Single-molecule real-time (SMRT ® ) sequencing was conducted using the C4 sequencing chemistry and P6 polymerase with one SMRT ® cell. De novo assembly of the PacBio reads were carried out using continuous long reads (CLR) following the hierarchical genome assembly process (HGAP) workflow (PacBioDevNet; Pacific Biosciences) as available in SMRT ® Analysis v2.3 (Chin et al., 2013). The functions of predicted proteins were annotated based on homologs (using SwissProt, http://www.uniprot. org/uniprot/), and clusters of ortholog groups were determined (using COG, http://www.ncbi.nlm.nih.gov/COG/). The NCBI Prokaryotic Genome Annotation Pipeline (PGAP) was employed to identify coding sequences (CDS) based on the best-placed reference protein set and GeneMarkS+. The complete genome sequence of E. hirae R17 was deposited in GenBank under the accession number CP015516 (chromosome) and CP015517 (plasmid pRZ1).

| Comparative genomic analysis
Three complete and annotated genomes of E. faecium DO (accession number from CP003583 to CP003586), E. faecalis V583 (accession number from AE016830 to AE016833), and E. hirae ATCC™9790 (accession number CP003504.1 and HQ724512.1) were available from GenBank and used as reference sequences in this study. The GC content of each sequence was calculated using the GC-Profile web server (Gao & Zhang, 2006). Orthologs between all three reference genomes and E. hirae R17 were identified using OrthoMCL software by reference to their protein sets (Fischer et al., 2011). Significant similarity was defined as E-value ≤1e-05, the length of the protein alignment ≥60%, and amino acid sequence identity ≥50%, and the best hits were selected and further analyzed.

| Features of the E. hirae R17 genome
An E. hirae strain with DAP resistance was screened, designated R17, and its complete genome sequence was determined. As shown in pseudogenes, and one CRISPR-cas sequence. The chromosome had a GC content of 36.96%, and it showed a clear GC skew at the origin of replication ( Figure 1). The size of the plasmid pRZ1 was found to be 73,574 bp with a GC content of 35.57%, encoding 66 ORFs, which covered 85.9% of the plasmid, and also encoding seven pseudogenes.

| Comparative genome analysis of E. hirae R17 with three other Enterococcus genomes
A 16S rRNA and an ANI phylogenetic tree describing the species belonging to genus Enterococcus were constructed based on 16S rRNA sequences and whole genomic sequences, respectively, and in which species within the same taxonomic groups were tightly clustered ( Figure 2). For the 16S rRNA phylogenetic tree, E. hirae R17 and reference strain E. hirae ATCC™9790 were observed to cluster together. Although E. hirae R17 is considered distantly related to isolates E. faecium DO and E. faecalis V583 as shown in the phylogenetic tree, its 16S rRNA sequence was 99% identical to the sequences in those two strains. The 16S rRNA phylogenetic tree also revealed that the E. hirae cluster was closely related to E. durans, but distant from E. faecium, E. mundtii, and E. faecalis ( Figure 2a). The ANI tree provides an in-depth phylogenetic relationship at the species and strain levels ( Figure 2b) loci. Orthologs analysis indicated that a total of 2,153 genes were shared between E. hirae R17 and ATCC™9790 (Figure 1), and some 278 genes and nine genomic regions (containing more than five consecutive strain-specific genes) were unique to E. hirae R17 (Table S1 and Figure 1). In total, 1,412 genes were shared across all four genomes, and 139 genes and five regions were unique in E. hirae R17 (Table S2 and Figure 1). The majority of these strain-specific genes are hypothetical proteins. The inner-and intraspecies genome comparison indicated that gene exchanges appeared to be a common feature among the Enterococcus spp. Strain-specific proteins of known function mostly belonged to the COG categories of carbohydrate transport and metabolism, transcription, and cell wall/membrane/envelope biogenesis (Table S1 and Table S2), suggesting that carbohydrate utilization and cell wall/membrane/envelope functions of E. hirae R17 may be different from the other three bacteria.
Two E. hirae R17-specific genomic regions (genomic region #4 [located from 2,304,674 to 2,319,265 bp] and #5 [located from 2,336,013 to 2,343,132 bp] shown in Table S2) were identified as putative GI, which may be the result of HGT because both regions contain transposases that mapped closely together. In particular, region #4 consisted of 10 hypothetical proteins, one for cell wall/membrane/ envelope biogenesis, one for signal transduction mechanisms, and one for replication, recombination and repair. Region #5 contained glycosyl transfer-related genes, and may be involved in cell wall/membrane/ envelope biogenesis.

| Genomic islands identified in E. hirae R17
GIs differ in GC content, codon usage, and k-mer frequencies from the rest of the genome due to their origins (Dhillon et al., 2015). A total of 13 areas of the sequence including 122,195 bp were predicted as GIs in E. hirae R17, of which 12 GIs were loci on the chromosome and one was found on the plasmid (Table S3). Genomic islands five and six may be a single GI as they are located very close together (separated by only 4,263 bp). Similarly, GI-9, -10, -11, and -12 likely form a single GI because their encoding ORFs are close together.
As shown in Table S3, many of the hypothetical proteins, pseudogenes, and transposon-related proteins were found in the 12 predicted GIs on the chromosome. Genomic island seven contained 10 hypothetical proteins, containing four possible pseudo genes, two predicted transposases, and one putative transposon, Tn552 DNAinvertase bin3. In addition to seven hypothetical proteins, GI-6 contained the HTH-type transcriptional repressor RghR, a toxin-encoding gene with ambiguous function, and one tyrosine recombinase XerClike gene. On GI-2 and -4, four and six genes were identified, respectively, encoding hypothetical proteins. Seven genes including six hypothetical proteins and one uncharacterized phage-related protein Lin1259/Lin1739 were identified on GI-5. Genomic islands 9, 10, 11, and 12 contained 29 ORFs, including 15 that were hypothetical proteins, one transposase, and one viral-enhancing factor. The hyaluronan synthase gene (AND73241.1), was considered to be a virulence factor in clinical isolates of E. faecium (Arias et al., 2009), and this gene was also located on GI-10. The dispersion of transposases identified on these GIs may relate to their acquisition by HGT.
Genomic islands three and eight consisted of six and five known genes, respectively. Genomic island three contained three galactosidaserelated proteins, two regulatory proteins, and one histidine kinase.
Three potassium-transporting related proteins, one sensor protein, and a tagatose aldolase were located in GI-8. There were 25 ribosomal proteins in GI-1, a common feature among the genomic islands of foodborne enterococcal genomes previously sequenced (Bonacina et al., 2016).  The sole GI identified on plasmid pRZ1 contained erythromycinand bacitracin-resistance-related genes, and these may function to

| Resistance to antimicrobial agents
Using the genome to query databases containing antibiotic resistanceencoding genes showed that the E. hirae R17 sequence possessed 33 resistance-encoding genes. Of these, 27 of these genes were located to the chromosome and six to the plasmid ( Interestingly, E. hirae R17 has two copies of the bacitracinresistance related gene bcrD, one on the chromosome and one on the plasmid. Alignments showed that the identity of the two predicted BcrD was only 54%. The most similar sequences to the BcrD from the chromosome (AND72321.1) and from the plasmid (AND73778.1) were separately used to make a phylogenetic tree, and the result showed that the BcrD on chromosome was far in phylogenetic distance from the BcrD on the plasmid (Figure 3). The chromosomal BcrD of E. hirae R17 was close to that of E. hirae and E. durans. However, the bcrD encoded on the plasmid was more similar to the genes of

| Antibacterial biocide and metal resistance genes
The predicted antibacterial biocide and metal resistance genes are listed in Table 3. Potential resistance genes related to hydrochloric acid, hydrogen peroxide, copper, silver, iron, and selenium were encoded on the chromosome, and the resistance genes related to ethidium bromide, sodium dodecyl sulfate, and tetraphenylphosphonium, as well as the cetrimide resistance-related gene lmrS were found on the plasmid. Interestingly, the predicted lmrS gene overlapped the

| Mobile genetic elements and genomic plasticity
The genome of E. hirae R17 contained 14 transposase genes, of which nine genes were located on the chromosome and five were located on the plasmid (Table 4) The numbers of tandem repeats (TR) identified hinted at genomic plasticity and bacterial environmental adaptation (Munita, Bayer, & Arias, 2015). Some 127 TR were identified with period sizes ranging from 6 to 441 bp on the chromosome and two tandem repeats with period sizes of 8 bp on the plasmid (Table S5). Many tandem repeats of E. hirae R17 were located in protein-coding regions, and thus could be part of the mechanism for the generation of targeted gene variation.

| Plasmid pRZ1 in E. hirae R17
Alignment of the 73.6 kb length sequence of plasmid pRZ1 in NCBI using BLASTN showed the highest query cover rate was 57%, which indicated that pRZ1 was a new plasmid. T A B L E 3 Biocide and metal resistance genes in E. hirae R17 predicted ORFs in the plasmid pRZ1 sequence, most were hypothetical proteins (Figure 4). Three conjugal transfer proteins indicate that plasmid pRZ1 may be capable of conjugal transfer to a new host.
A 16.3 kb-region (AND73768.1-AND73783.1) containing four IS6 family elements and six genes associated with the resistance of bacitracin, macrolide antibiotics, and tetracycline (Figure 4), that showed similarity to the chromosomal sequence of Streptococcus gallolyticus UCN34, S. aureus subsp. aureus JS395, and E. faecium E506 with 44% match length, indicating that the architecture of the region was unique. Plasmid pRZ1 harbored a type III secretion system protein, PrgN (AND73763.1), and a type IV secretory protein (AND73805.1), which may contribute to the pathogenicity of the plasmid.

| DISCUSSION
Enterococci are bacteria known to be involved in the dissemination of resistance to antimicrobial compounds . In earlier reports, E. hirae were found more frequently in municipal waste water, river, soil, and animals, and less frequently identified in sources such as fruits, vegetable, and food products of animal origin (Chubiz, Lee, Delaney, & Marx, 2012;Munita et al., 2014). Therefore, the isolation and identification of E. hirae R17 from a retail pork sample was of interest because of its MDR phenotype and potential pathogenicity. In this study, a genomic comparison between E. hirae R17 and other selected enterococcal reference genomes was made. GIs, virulence-encoding genes, antimicrobial resistance genes, among others were investigated.
Genome comparison of the foodborne E. hirae R17 to inter-and intraspecies isolates allowed us to determine the phylogeny of E. hirae and to accurately quantify the genomic diversity between strains of the genus Enterococcus. The ANI phylogenetic tree provides stronger discriminatory power than the 16S rRNA tree and enables enterococcal strain differentiation. The phylogenetic relationships among Enterococcus strains constructed by ANI are more consistent with the genome-and protein-based phylogenetic trees published previously than the 16S rRNA tree (Bonacina et al., 2016). The function of the E. hirae R17-specific proteins suggested that this isolate may have unique carbohydrate utilization processes and a cell wall/membrane/ envelope constitution. In a recent report, the overall genomes of six pig fecal E. hirae isolates were found to be highly similar to each other and to the type strain E. hirae ATCC™9790, but the gene content and structural differences are likely related to ecological niche specialization and may influence their competitiveness among strains (Katyal et al., 2016).
Genomic islands are known to contribute to the development of pathogenicity and the diversity that can be contained within a single bacterial species. Several studies have reported that the enterococci that cause infection in hospitalized patients are different from those isolates that colonize the gastrointestinal tract of the healthy human host, suggesting that the gain and loss of mobile genetic elements, rather than evolutionary descent, may be the most important driving force in determining virulence-associated properties in isolates (Qin et al., 2012;Waters et al., 2011). In this study, two possible virulence factors, a toxin-encoding gene and a hyaluronidase gene, located on GIs, strongly indicated that these pathogenesis-related genes originated by HGT.
Capsular and other cell envelope polysaccharides of several Grampositive bacteria are known to play important roles in virulence and protective immunity (Thurlow, Thomas, Fleming, & Hancock, 2009a;Thurlow, Thomas, & Hancock, 2009b). Among the predicted virulence factors identified in the genome of E. hirae R17, some may contribute to pathogenesis and evasion of the host innate immune response Lee et al., 2015). Similar to other Gram-positive bacteria, MSCRAMM (Microbial Surface Components Recognizing Adhesive Matrix Molecules)-like cell wall-anchored proteins were identified in E. hirae R17, including the biofilm-associated Ebp pili, known to be important for endocarditis and UTI in animal models (Sillanpaa et al., 2010). The hyaluronidase function encoded by the hyl gene is primarily found in clinical E. faecium isolates and is much less common in isolates recovered from other sources (Rice et al., 2003).
The conjugative transfer of a plasmid containing the hyl gene makes the transconjugant more virulent in a mouse peritonitis models, suggesting a role of hyaluronidase in virulence (Arias et al., 2009). The endocarditis enterococcal antigen (EfaA) is commonly found both among E. faecalis and E. faecium, and is presumed to be involved in the adhesion of enterococci to biotic and abiotic surfaces or evasion of the immune response (Abriouel et al., 2008). In addition, the binding to extracellular matrix components of a number of other surface proteins may also contribute to the virulence of E. hirae R17 by acting as adhesions.
The clinical importance of the genus Enterococcus is directly related to its antibiotic resistance, which significantly contributes to the risk of colonization and infection. Relative to the streptococci, enterococci are intrinsically resistant to many commonly used antimicrobial agents, including most cephalosporins, all semisynthetic penicillins, and clindamycin. The intrinsic resistance to clindamycin is mediated by the gene product of lsa, which is also found in the E. hirae R17 genome, although its mechanism of action remains poorly defined (Bonacina et al., 2016;Hollenbeck & Rice, 2012). Enterococcal isolates are often tolerant to β-lactam antibiotics, likely due to the expression of lowaffinity penicillin-binding proteins (PBPs) (Arias & Murray, 2012).
Although four genes encoding PBPs (denoted as PBP2b, PBP1b, PBP1a, and PBP2x) were identified in the E. hirae R17 chromosome, this bacterium is susceptible to penicillin and ampicillin.
In addition to intrinsic resistance, enterococci have been extraordinarily successful at exploiting mechanisms of gene mobility and exchange to rapidly acquire resistance to virtually any antimicrobial agent in current use . The regular use of antimicrobials agents in modern animal production potentially increases selection pressure on enterococci to become resistant for niche adaptation and also increases the risk of drug-resistant genes spreading from those bacteria in food-producing animals to humans (Bonacina et al., 2016;Koluman & Dikici, 2013;Peng et al., 2014;Van Boeckel et al., 2015). Analysis of the drug-resistance genes identified in E. hirae R17 suggests that this bacterium acquires these features by horizontal means involving plasmids, transposons, and insertion sequences. As F I G U R E 4 A schematic representation of the genetic map of plasmid pRZ1. Genes are denoted by arrows and are colored based on gene function: green arrows represent hypothetical proteins; blue arrows represent transposases; yellow arrows represent predicted antibiotic resistance genes; red arrows represent other encoding genes of known function; blue stripped represent genomic island (GIP1). The innermost circle presents the GC skew [(G-C)/(G+C)], values >0 (red), values <0 (green). The outer circle represents the GC content, above median GC content (red), less than the median (green) an example, a 15 kb resistance gene island in plasmid pRZ1 was found to be composed of four IS6 elements, which may function to disseminate genes to neighboring bacteria. This cluster contained genes encoding resistance to bacitracin (bcrD, bcrB, bcrA, and bcrR), macrolides (ermB), and tetracycline (tetK and tetT).
Bacitracin is used both in clinical practice and as an animal growth promoter in some countries. There were recent reports of high-level bacitracin resistance in the genus Enterococcus (Kelesidis, 2015). The bcr operon is often clustered in one region on the enterococcal genome; however, it remains unclear why there is one orphan bcrD gene on the chromosome of E. hirae R17 and a complete bcr operon on the plasmid.
DAP is a lipopeptide antibiotic frequently used as a last-resort antibiotic against vancomycin-resistant enterococci. There are limited treatment options for DAP-resistant enterococcal isolates (Tang et al., 2015). Other than E. faecalis and E. faecium, this is the first report of DAP-resistant bacterium in E. hirae. The genetic basis and supporting mechanism for this DAP resistance of enterococci is complex and poorly understood. The changes in regulatory systems that control the cell envelope response to stress may be critical to the development of DAP resistance in enterococci clinical isolates . In the chromosome of E. hirae R17, the DAP resistance-associated genes mprF and cls were found, and these may contribute to the observed DAP resistance of E. hirae R17 (Arias et al., 2012). Further studies on the mechanisms of DAP resistance in E. hirae R17 are warranted, because DAP is not typically used in animal production in China due to its high price.
Antibacterial biocides and metals can contribute to the development and maintenance of antibiotic resistance in bacterial communities through co-selection (Dortet, Anguel, Fortineau, Richard, & Nordmann, 2013;Pal et al., 2014). Based on the analysis of the genome of E. hirae R17, genes encoding resistance to several heavy metals, including copper, iron, silver, and also selenium, were identified. can exploit the instability of TRs to reversibly shut down or modulate the function of specific genes, allowing adaptation to changing environments on short evolutionary time scales without an increased overall mutation rate . Therefore, the presence of IS sequences and several TRs indicate that E. hirae R17 has strong genomic plasticity with potential for environmental adaption. The results of this study also suggest that the R17 plasmid may be responsible for much of the horizontal gene transfer that has allowed antibiotic, virulence, and antibacterial biocide traits to converge in clinical adapted lineages.
Mobile elements may provide an accessory pool of adaptive genes that, if acquired by a bacterium, may have the potential to enhance survival in select environments. Paradoxically, the CRISPR-cas system in bacteria confers resistance against foreign genetic elements, providing adaptive immunity against foreign DNA including plasmids and bacteriophage (Chuang et al., 2014). One hypothesis suggests that an inverse relationship exists between the acquisition of one or more antimicrobial resistance-encoding genes and the possession of complete CRISPR loci. This appears to hold for members of the recently emerged high-risk E. faecium and E. faecalis lineages, as such bacteria lack these loci. Surprisingly, this hypothesis does not appear to hold in the case of E. hirae R17, as a CRISPR1-cas loci was identified ( Table 1).
The incorporation of these components may have been stimulated by modern antibiotic therapy, promoting the acquisition of foreign elements conferring antibiotic resistance, decreasing genome stability/ increasing plasticity, and enabling the colonization of new habitats, including the antibiotic-laden hospital environment (O'Driscoll & Crank, 2015).
It has been suggested previously that the pan-genome of E. faecium is essentially unlimited in size. This may indicate that E. faecium has the potential to acquire and incorporate exogenous DNA into its genome (van Schaik et al., 2010). In contrast, our previously limited genomic information of E. hirae has hindered the prediction of core and pan genomic information of this species. Additionally, there is a need for a multistrain genomic analysis of E. hirae to describe its fundamental biology and to determine the causes of its emergence as a possible opportunistic pathogen.

| CONCLUSION
The E. hirae R17 genome sequence provides a starting point to study the evolutionary history and the potential pathobiology of this species, particularly those isolates recovered from food-related samples. We found that the isolate originating from retailed pork acquired multiple antibiotic resistance genes and also has genes that may play roles in colonization and the infection of animal food consumers. The genomic flexibility of E. hirae R17, with its IS elements, transposase, and TRs, allows rapid adaptation to new environmental niches such as those that are found in food sterilization and storage sites. Genome-wide studies, which will be facilitated by the sequence data presented here, are therefore needed to increase our understanding of the basic biology of E. hirae and to identify genes that may present risks for public health.