Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 26, Number 12—December 2020

Transmission of Multidrug-Resistant Salmonella enterica Subspecies enterica 4,[5],12:i:- Sequence Type 34 between Europe and the United States

Ehud ElnekaveComments to Author , Samuel L. Hong, Seunghyun Lim, Dave Boxrud, Albert Rovira, Alison E. Mather, Andres Perez, and Julio Alvarez
Author affiliations: Hebrew University of Jerusalem, Jerusalem, Israel (E. Elnekave); University of Minnesota, St. Paul, Minnesota, USA (E. Elnekave, S. Lim, A. Rovira, A. Perez, J. Alvarez); University of Leuven, Leuven, Belgium (S.L. Hong); Minnesota Department of Health, St. Paul (D. Boxrud); Quadram Institute Bioscience, Norwich, UK (A.E. Mather); University of East Anglia, Norwich (A.E. Mather); Universidad Complutense, Madrid, Spain (J. Alvarez)

Cite This Article


Multidrug-resistant Salmonella enterica subspecies enterica 4,[5],12:i:- sequence type 34 represents a worldwide public health risk. To determine its origin in the United States, we reconstructed a time-scaled phylogeny with a discrete trait geospatial model. The clone in the United States was introduced from Europe on multiple occasions in the early 2000s.

Since the late 1990s, reports of an emerging multidrug-resistant Salmonella enterica subspecies enterica serotype 4,[5],12:i:- strain have been published in Europe (1). This strain is a monophasic variant of Salmonella Typhimurium, predominantly resistant to ampicillin, streptomycin, sulfonamides, and tetracycline (ASSuT). Its rapid increase in North America after 1998 has also been described (2). However, precise knowledge of the time of introduction and the initial influx of clinical cases caused by this serotype in the United States is not available because of inconsistent reporting before 2004 (3).

Previously, on the basis of high genetic similarity between Salmonella 4,[5],12:i:- sequence type (ST) 34 isolates from the United States and Europe and Salmonella Typhimurium strains from Europe, we suggested a European origin for the Salmonella 4,[5],12:i:- ST34 clade (4). With this study, we aimed to reconstruct a time-scaled phylogeny of the emerging ST34 clade by using a Bayesian modeling approach to determine its origin in United States.

The Study

Figure 1

Demographic reconstruction and phylogenetic analysis of Salmonella enterica subspecies enterica 4,[5],12:i:- sequence type 34 isolates. A) Demographic reconstruction (subset 2) shows the population exponential growth over time. The red line indicates the median effective population size with 95% highest posterior density credible interval (gray). B) Time-scaled phylogenetic analysis of isolates in subset 2 (n = 110 sequences after duplicates removal). Isolates were collected from multiple sources in the United States (blue) and Europe (red) during 2008–2017.An asymmetric discrete trait analysis model was used to reconstruct the locations on the nodes. The nodes, branches, and tree tips were annotated according to the collection location. The 95% highest posterior density credible intervals of node heights are indicted with transparent purple bars. The posterior probability for all inferred ancestral locations was >70%. The isolate source (food product, gray; human, purple; livestock, yellow; and other, green) is depicted in the heatmap appended to the tree tips.

Figure 1. Demographic reconstruction and phylogenetic analysis of Salmonella enterica subspecies enterica4,[5],12:i:- sequence type 34 isolates. A) Demographic reconstruction (subset 2) shows the population exponential growth...

Figure 2

Summary of posterior estimates of all 10 subsets of sequences of Salmonella enterica subspecies enterica 4,[5],12:i:- sequence type (ST) 34 collected from multiple sources in the United States and Europe during 2008–2017. A) Inferred time (year) of the most recent common ancestor of the emerging Salmonella 4,[5],12:i:- ST34 clade. B) Estimated mutation rate (uncorrelated log-normally distributed mean parameter). C–F) Number of unobserved transitions from Europe to the United States (C) and United States to Europe (E) along each branch (Markov jumps) and total phylogenetic tree length spent (Markov rewards) in Europe (D) and the United States (F). Violin plots illustrate the posterior distribution and probability density of each subset. Dashed red vertical lines indicate average posterior value; red shaded areas indicate average 95% highest posterior density credible interval of all subsets.

Figure 2. Summary of posterior estimates of all 10 subsets of sequences of Salmonella enterica subspecies enterica4,[5],12:i:- sequence type (ST) 34 collected from multiple sources in...

We used publicly available whole-genome sequences of 1,431 Salmonella 4,[5],12:i:- ST34 isolates from the United States and Europe from 2008–2017, including sequences from 690 isolates from Europe (mainly from the United Kingdom and Denmark) and 741 isolates from multiple US states (Appendix 1). We used BEAST version 1.8.4 (5) to estimate divergence times, mutation rates, and location trait transitions. We applied the modeling approach to 10 subsets of 112 sequences selected from the study population. These sequences represented 33% (474/1,431) of the study population and included 242 sequences from Europe (76% from humans, 8% from food products, 8% from livestock, and 8% from other sources) and 232 from the United States (62% from humans, 13% from food products, 21% from livestock, and 3% from other sources). Time-scaled phylogenies of each subset were reconstructed by using a general time-reversible nucleotide substitution model, an uncorrelated lognormal relaxed molecular clock, and an exponential growth coalescent model with asymmetric trait transitions (Figure 1; Appendix 2 Figures 1–10). All time-scaled phylogenies presented similar topology to a maximum-likelihood phylogeny constructed by using all 1,431 study isolates (based on visual inspection; Appendix 2 Figure 11). Overall, averaged estimates from all subsets were in agreement as follows (Figure 2): the evolutionary rate was 3.64 × 10−7 substitutions/site/year (95% highest posterior density [HPD] 2.65–4.64 × 10−7), which corresponds to an accumulation of »1–2 single-nucleotide polymorphisms/genome/year; the time to most recent common ancestor was 1994 (95% HPD 1988–2000); the number of collection location state transitions (Markov jumps) from Europe to the United States was 7.7 (95% HPD 5.9–9.3) and from the United States to Europe was 0.8 (95% HPD 0–2.2); and the waiting times (in years; Markov rewards) were 519.9 (95% HPD 393.1–667.8) in Europe and 318.6 (95% HPD 234.0–417.6) in the United States. The exponential growth rate of the population was estimated at 0.316/year (averaged across all subset means; Figure 1). In addition, the estimated (averaged) main introduction into the United States was 2004 (95% HPD 2000–2006; Appendix 2 Table 1). The occurrence of several additional smaller introductions was suggested by 48 sequences (6 from Europe and 42 from the United States). The 6 sequences from Europe were recovered from human sources; among sequences from the United States, 28 were from humans, 5 from food products, 6 from livestock, and 3 from other sources. Information on recent international travel was obtained for 22/28 of US isolates from humans, 2 of whom had traveled (1 to the Philippines and the other to France [S. Meyer et al., Minnesota Department of Health, pers. comm., 2019 Sep 23]).

Among the 1,431 Salmonella 4,[5],12:i:- ST34 sequences, 978 (68.34%) had genetic determinants contributing to the ASSuT profile, 108 (7.55%) conferred resistance to quinolones, and 82 (5.73%) conferred resistance to extended-spectrum cephalosporins. The probability of harboring most predominant acquired antimicrobial resistance genes (AARGs) conferring the resistance phenotypes described above was significantly higher for sequences of US isolates (odds ratio 2.37–26.05; Table). Yet associations between the collection location and the presence of blaCTX-M or qnrS1 genes were not significant (Table). In addition, AARGs conferring resistance to colistin (mcr-1/mcr-3/mcr-5; Appendix 1) were detected in isolates from Europe only (n = 5).


Salmonella 4,[5],12:i:- ST34 was introduced into the United States from Europe on multiple occasions since the beginning of the 21st century. The main introduction occurred in 2004; additional independent introductions resulted in small clades for which the predominant sources were human travelers and imported food products. Human travelers (6) and imported food products (7) have been described as potential vehicles for introduction of salmonellae.

The date of introduction of the main clade into the United States is later than the first peer-reviewed report of a Salmonella 4,[5],12:i:- infection in the country in 1998 (2). However, given the antimicrobial susceptibility profile of isolates from that report (mostly not ASSuT) (2), they most likely belonged to the nonemerging ST19 clade, which was described elsewhere (4). In addition, the incidence of Salmonella 4,[5],12:i:- in humans increased only modestly (9.5%) during 2006–2011 but increased dramatically (78.3%) during 2011–2016 (8). A similar increase in detection after 2011 was described for clinical cases in swine from the midwestern United States (9). The difference between the date of main introduction into the United States found in this study and the later sharp increase in its prevalence in animals and humans may in part result from changes in reporting practices and increasing awareness (8). However, the increase since 2011 can be the result of rapid propagation of the ST34 population, possibly associated with swine (4). Moreover, White et al. (10) recently reported that according to the National Antibiotic Resistance Monitoring System, the percentage of ASSuT-resistant Salmonella 4,[5],12:i:- from humans increased from 17% in 2009 to 59.1% in 2015 (out of all Salmonella 4,[5],12:i:- clinical isolates from humans). This increase probably resulted to a large extent from ST34 strains, in which this phenotype is predominant. The estimated exponential yearly growth rate determined in our model (0.316/year), which corresponds to a population doubling time of 2.2 years, is in agreement with this dramatic increase of the ST34 population.

The presence of AARGs conferring resistance to quinolones and extended-spectrum cephalosporins has mainly been observed since 2014 and may be biased by the lack of sequences before 2013 (Appendix 2 Figure 12). Yet AARGs conferring resistance to quinolones were not found in Salmonella 4,[5],12:i:- ST34 strains from Europe collected before 2010 (1), and therefore our findings may reflect an increasing prevalence of these resistance determinants. Given time and overall unidirectionality of Salmonella 4,[5],12:i:- ST34 transmission from Europe to the United States, it is likely that the acquisition of AARGs to quinolones occurred independently in the United States and in Europe. Yet introduction of resistant strains from the United States to Europe is also possible. Contributors to the acquisition of resistance in the United States might be the approval for enrofloxacin use in swine in the United States since 2008 (11) and the potential dissemination of plasmids harboring AARGs to quinolones between Salmonella serotypes (12). Independent acquisition of resistance to quinolones by Salmonella in Asia has also been suggested (13). The presence of mcr resistance genes conferring resistance to colistin in sequences from Europe (n = 5) is alarming, given their recent worldwide spread (14). However, further investigation of the travel history associated with these cases may be required because the acquisition of mcr genes may be travel associated (15). The spread of Salmonella 4,[5],12:i:- ST34 from Europe to the United States and the presence of plasmid-mediated resistance genes to key antimicrobial classes such as quinolones, extended-spectrum cephalosporins, and colistin in this clade further highlights its potential risk to public health and emphasizes the need for robust surveillance and mitigation programs for such transboundary pathogens.

Dr. Elnekave is a veterinarian and epidemiologist who has been working as a postdoctoral fellow in the Veterinary Population Medicine Department, University of Minnesota, and is a faculty member in the Koret School of Veterinary Medicine at the Robert H. Smith Faculty of Agricultural, Food and Environmental Sciences at the Hebrew University of Jerusalem, Israel. His primary research interest is the epidemiology of bacterial foodborne pathogens. He uses genetic tools and Bayesian models to study pathogen emergence, evolution, and antimicrobial resistance.



We thank staff members from the Minnesota Department of Health Foodborne, Waterborne, Vectorborne, and Zoonotic Diseases Section for providing the travel history records and for their advice on this work.

This work was supported by the Global Food Venture-MnDRIVE Initiative, the National Institute of Food and Agriculture (Animal Health Formula Fund project MIN-62-091) of the US Department of Agriculture, the Rapid Agricultural Response Fund, and the Swine Disease Eradication Center at the University of Minnesota. In addition, the United States-Israel Binational Agricultural Research and Development Fund awarded a Vaadia-BARD postdoctoral fellowship (no. FI-565-17) to E.E. A.E.M. is a Food Standards Agency research Fellow and is supported by the Biotechnology and Biological Sciences Research Council Institute Strategic Programme Microbes in the Food Chain BB/R012504/1 and its constituent project BBS/E/F/000PR10348 (Theme 1, Epidemiology and Evolution of Pathogens in the Food Chain). The Ramón y Cajal postdoctoral contract from the Spanish Ministry of Economy, Industry and Competitiveness (MINECO) (RYC-2016-20422) was awarded to J.A.



  1. Petrovska  L, Mather  AE, AbuOun  M, Branchu  P, Harris  SR, Connor  T, et al. Microevolution of monophasic Salmonella Typhimurium during epidemic, United Kingdom, 2005–2010. Emerg Infect Dis. 2016;22:61724. DOIPubMedGoogle Scholar
  2. Agasan  A, Kornblum  J, Williams  G, Pratt  CC, Fleckenstein  P, Wong  M, et al. Profile of Salmonella enterica subsp. enterica (subspecies I) serotype 4,5,12:i:- strains causing food-borne infections in New York City. J Clin Microbiol. 2002;40:19249. DOIPubMedGoogle Scholar
  3. Centers for Disease Control and Prevention. National Salmonella surveillance overview. Atlanta: US Department of Health and Human Services; 2011.
  4. Elnekave  E, Hong  S, Mather  AE, Boxrud  D, Taylor  AJ, Lappi  V, et al. Salmonella enterica serotype 4,[5],12:i:- in swine in the United States midwest: an emerging multidrug-resistant clade. Clin Infect Dis. 2018;66:87785. DOIPubMedGoogle Scholar
  5. Drummond  AJ, Rambaut  A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214. DOIPubMedGoogle Scholar
  6. Williamson  DA, Lane  CR, Easton  M, Valcanis  M, Strachan  J, Veitch  MG, et al. Increasing antimicrobial resistance in nontyphoidal Salmonella isolates in Australia from 1979 to 2015. Antimicrob Agents Chemother. 2018;62:e02012-17. DOIPubMedGoogle Scholar
  7. Gould  LH, Kline  J, Monahan  C, Vierk  K. Outbreaks of disease associated with food imported into the United States, 1996–2014. Emerg Infect Dis. 2017;23:5258. DOIPubMedGoogle Scholar
  8. Centers for Disease Control and Prevention. National Salmonella surveillance annual report, 2016: Atlanta: US Department of Health and Human Services; 2018.
  9. Hong  S, Rovira  A, Davies  P, Ahlstrom  C, Muellner  P, Rendahl  A, et al. Serotypes and antimicrobial resistance in Salmonella enterica recovered from clinical samples from cattle and swine in Minnesota, 2006 to 2015. PLoS One. 2016;11:e0168016. DOIPubMedGoogle Scholar
  10. White  PL, Green  AL, Holt  KG, Hale  KR. Multidrug-resistant Salmonella enterica subspecies I serovar 4,[5],12:i:- isolates recovered from Food Safety and Inspection Service-regulated products and food animal ceca, 2007-2016. Foodborne Pathog Dis. 2019;16:67986. DOIPubMedGoogle Scholar
  11. Food and Drug Administration. NADA 141–068 in the approved animal drug products (Green Book) [cited 2020 Aug 31].
  12. Elnekave  E, Hong  SL, Lim  S, Hayer  SS, Boxrud  D, Taylor  AJ, et al. Circulation of plasmids harboring resistance genes to quinolones and/or extended-spectrum cephalosporins in multiple Salmonella enterica serotypes from swine in the United States. Antimicrob Agents Chemother. 2019;63:e02602-18. DOIPubMedGoogle Scholar
  13. Mather  AE, Phuong  TLT, Gao  Y, Clare  S, Mukhopadhyay  S, Goulding  DA, et al. New variant of multidrug-resistant Salmonella enterica serovar Typhimurium associated with invasive disease in immunocompromised patients in Vietnam. MBio. 2018;9:e01056-18. DOIPubMedGoogle Scholar
  14. McEwen  SA, Collignon  PJ. Antimicrobial resistance: a One Health perspective. Microbiol Spectr. 2018;6.PubMedGoogle Scholar
  15. Arnott  A, Wang  Q, Bachmann  N, Sadsad  R, Biswas  C, Sotomayor  C, et al. Multidrug-resistant Salmonella enterica 4,[5],12:i:- sequence type 34, New South Wales, Australia, 2016–2017. Emerg Infect Dis. 2018;24:7513. DOIPubMedGoogle Scholar




Cite This Article

DOI: 10.3201/eid2612.200336

Original Publication Date: November 09, 2020

Table of Contents – Volume 26, Number 12—December 2020

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.



Please use the form below to submit correspondence to the authors or contact them at the following address:

Ehud Elnekave, Koret School of Veterinary Medicine, Robert H. Smith Faculty of Agricultural, Food and Environmental Sciences, The Hebrew University, PO Box 12, Rehovot 76100, Israel

Send To

10000 character(s) remaining.


Page created: September 30, 2020
Page updated: November 19, 2020
Page reviewed: November 19, 2020
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.