Skip to content
1887

Journal of Medical Microbiology logo Journal of Medical Microbiology

Volume 74, Issue 5

A high-resolution genomic and phenotypic analysis of resistance evolution of an strain from a critically unwell patient treated with piperacillin/tazobactam Open Access

Abstract

Resistance to the -lactam/-lactamase inhibitor (BL/BLI) combination antibiotic piperacillin/tazobactam (TZP) predominantly occurs via -lactamase enzymes, also leading to resistance to third-generation cephalosporins (3GCs). However, if -lactamases inactive against 3GCs and inhibited by tazobactam are expressed at high levels, leading to enzyme hyperproduction, the surplus enzyme escapes inhibition by tazobactam and inactivates the antibiotic piperacillin.

Understanding this mechanism is clinically relevant, as enzyme hyperproduction can emerge upon antibiotic administration, resulting in treatment failure despite initial resistance profiles supporting TZP use.

Our aim was to determine whether this was a case of within-patient evolution and by what mechanism or an acquisition of a second unrelated, more resistant, strain.

Whole-genome sequencing was performed on the isolate to determine the genetic basis of resistance. We also assessed the impact of TZP exposure on the amplification of the gene and monitored the stability of gene copy number over 5 days in the absence of antibiotic pressure. In addition, we determined the MICs of ceftriaxone and TZP, with TZP MIC contextualized in relation to gene copy number and resistance levels.

We report the identification of an isolate that developed resistance to TZP during patient treatment but maintained sensitivity to ceftriaxone. We show that TZP resistance evolved via IS-mediated duplication of a containing transposable unit on a plasmid, resulting in hyperproduction of TEM-1 -lactamase, and that ten copies of induce resistance greater than 32 times the MIC. Furthermore, under experimental conditions, exposure to TZP further increases amplification of , whereas, in the absence of TZP, gene copy number of IS and remains stable over 5 days, despite a 48,205 bp genome size increase compared to the pre-amplification isolate. We additionally detect phenotypic changes that might indicate host adaptation potentially linked to the additional genes that are amplified.

Our analysis advances the understanding of infections caused by isolates evolving -lactamase hyperproduction, which represents a complex problem in both detection and treatment. As 40% of antibiotics active against WHO priority pathogens in the pre-clinical pipeline are BL/BLI combinations, further investigations are of urgent concern.

  • Received:
  • Accepted:
  • Published Online:
Keyword(s): antimicrobial resistance , enzyme hyperproduction , Escherichia coli , gene amplification , IS26 , TEM-1 , β-lactam and β-lactamase
Funding
This study was supported by the:
  • Royal Academy of Engineering (Award RCSRF2021\11\15)
    • Principal Award Recipient: PaulA. Hoskisson
  • Biotechnology and Biological Sciences Research Council (Award BB/V011278/2)
    • Principal Award Recipient: EvaHeinz
  • Biotechnology and Biological Sciences Research Council (Award BB/V011278/1)
    • Principal Award Recipient: EvaHeinz
  • Wellcome Trust (Award 217303/Z/19/Z)
    • Principal Award Recipient: EvaHeinz
  • Medical Research Foundation (Award MR/R015678/1)
    • Principal Award Recipient: AliceJ. Fraser
© 2025 The Authors
  • This is an open-access article distributed under the terms of the Creative Commons Attribution License. This article was made open access via a Publish and Read agreement between the Microbiology Society and the corresponding author’s institution.
Loading

Article metrics loading...

/content/journal/jmm/10.1099/jmm.0.002018
2025-05-19
2026-04-13

Metrics

Loading full text...

Full text loading...

/deliver/fulltext/jmm/74/5/jmm002018.html?itemId=/content/journal/jmm/10.1099/jmm.0.002018&mimeType=html&fmt=ahah

References

  1. Bush K, Bradford PA. β-Lactams and β-lactamase inhibitors: an overview. Cold Spring Harb Perspect Med 2016; 6:a025247 [View Article] [PubMed]
     [Google Scholar]
  2. National institute for health care excellence (NICE) British National Formulary Treatment Summaries: cephalosporins; 2023 https://bnf.nice.org.uk/treatment-summaries/cephalosporins/
  3. National institute for health care excellence (NICE) British National Formulary Treatment Summaries: carbapenems; 2023 https://bnf.nice.org.uk/treatment-summaries/carbapenems
  4. Murray CJ, Ikuta KS, Sharara F, Swetschinski L, Aguilar GR et al. Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet 2022; 399:629–655 [View Article] [PubMed]
     [Google Scholar]
  5. Barton GJ, Morecroft CW, Henney NC. A survey of antibiotic administration practices involving patients with sepsis in UK critical care units. Int J Clin Pharm 2020; 42:65–71 [View Article] [PubMed]
     [Google Scholar]
  6. Wilson APR. Sparing carbapenem usage. J Antimicrob Chemother 2017; 72:2410–2417 [View Article] [PubMed]
     [Google Scholar]
  7. Suzuki Y, Sato T, Fukushima Y, Nakajima C, Suzuki Y et al. Contribution of β-lactamase and efflux pump overproduction to tazobactam-piperacillin resistance in clinical isolates of Escherichia coli. Int J Antimicrob Agents 2020; 55:105919 [View Article] [PubMed]
     [Google Scholar]
  8. Crowley B, Benedí VJ, Doménech-Sánchez A. Expression of SHV-2 β-lactamase and of reduced amounts of OmpK36 porin in Klebsiella pneumoniae results in increased resistance to cephalosporins and carbapenems. Antimicrob Agents Chemother 2002; 46:3679–3682 [View Article] [PubMed]
     [Google Scholar]
  9. Sanbongi Y, Suzuki T, Osaki Y, Senju N, Ida T et al. Molecular evolution of β-lactam-resistant Haemophilus influenzae: 9-year surveillance of penicillin-binding protein 3 mutations in isolates from Japan. Antimicrob Agents Chemother 2006; 50:2487–2492 [View Article] [PubMed]
     [Google Scholar]
  10. Lee J, Oh CE, Choi EH, Lee HJ. The impact of the increased use of piperacillin/tazobactam on the selection of antibiotic resistance among invasive Escherichia coli and Klebsiella pneumoniae isolates. Int J Infect Dis 2013; 17:e638–43 [View Article] [PubMed]
     [Google Scholar]
  11. van Duin D, Bonomo RA. Ceftazidime/avibactam and ceftolozane/tazobactam: second-generation β-lactam/β-lactamase inhibitor combinations. Clin Infect Dis 2016; 63:234–241 [View Article] [PubMed]
     [Google Scholar]
  12. Turner MS, Andersson P, Bell JM, Turnidge JD, Harris T et al. Plasmid-borne blaSHV genes in Klebsiella pneumoniae are associated with strong promoters. J Antimicrob Chemother 2009; 64:960–964 [View Article] [PubMed]
     [Google Scholar]
  13. Hubbard A, Mason J, Roberts P, Parry CM, Corless C et al. Piperacillin/tazobactam resistance in a clinical isolate of Escherichia coli due to IS26-mediated amplification of blaTEM-1B. Nat Commun 2020; 11:4915 [View Article] [PubMed]
     [Google Scholar]
  14. Hansen KH, Andreasen MR, Pedersen MS, Westh H, Jelsbak L et al. Resistance to piperacillin/tazobactam in Escherichia coli resulting from extensive IS26-associated gene amplification of blaTEM-1. J Antimicrob Chemother 2019; 74:3179–3183 [View Article] [PubMed]
     [Google Scholar]
  15. Schechter LM, Creely DP, Garner CD, Shortridge D, Nguyen H et al. Extensive gene amplification as a mechanism for piperacillin-tazobactam resistance in Escherichia coli. mBio 2018; 9:e00583-18 [View Article] [PubMed]
     [Google Scholar]
  16. San Millan A, Escudero JA, Gifford DR, Mazel D, MacLean RC. Multicopy plasmids potentiate the evolution of antibiotic resistance in bacteria. Nat Ecol Evol 2016; 1:1–8 [View Article]
     [Google Scholar]
  17. Rosenkilde CE, Munck C, Porse A, Linkevicius M, Andersson DI et al. Collateral sensitivity constrains resistance evolution of the CTX-M-15 β-lactamase. Nat Commun 2019; 10:618 [View Article] [PubMed]
     [Google Scholar]
  18. Edwards T, Heinz E, van Aartsen J, Howard A, Roberts P et al. Piperacillin/tazobactam resistant, cephalosporin susceptible Escherichia coli bloodstream infections are driven by multiple acquisition of resistance across diverse sequence types. Microbiology 2020 [View Article]
     [Google Scholar]
  19. Maillard A, Delory T, Bernier J, Villa A, Chaibi K et al. Effectiveness of third-generation cephalosporins or piperacillin compared with cefepime or carbapenems for severe infections caused by wild-type AmpC β-lactamase-producing Enterobacterales: a multi-centre retrospective propensity-weighted study. Int J Antimicrob Agents 2023; 62:106809 [View Article] [PubMed]
     [Google Scholar]
  20. Harris PN, Tambyah PA, Lye DC, Mo Y, Lee TH et al. Effect of piperacillin-tazobactam vs meropenem on 30-day mortality for patients with E. coli or Klebsiella pneumoniae bloodstream infection and ceftriaxone resistance: a randomized clinical trial. JAMA 2018; 320:984–994 [View Article] [PubMed]
     [Google Scholar]
  21. Kahlmeter G, Brown D, Goldstein F, MacGowan A, Mouton J et al. European committee on antimicrobial susceptibility testing (EUCAST) technical notes on antimicrobial susceptibility testing. Clin Microbiol Infect 2006; 12:501–503 [View Article] [PubMed]
     [Google Scholar]
  22. The European Committee on Antimicrobial Susceptibility Testing Breakpoint tables for the interpretation of MICs and zone diameters. Version 10.0; 2020 www.eucast.org
  23. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 2014; 30:2114–2120 [View Article] [PubMed]
     [Google Scholar]
  24. Andrews S. FastQC: a quality control tool for high throughput sequence data. Cambridge, United Kingdom: Babraham Bioinformatics, Babraham Institute; 2010
  25. Oxford Nanopore Technologies Guppy. n.d https://communitynanoporetechcom/downloads
  26. Wick RR, Menzel P. Filtlong; 2017 https://github.com/rrwick/Filtlong
  27. Wick R. Porechop: adapter trimmer for Oxford Nanopore reads; 2018 https://github.com/rrwick/Porechop
  28. Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol 2017; 13:e1005595 [View Article] [PubMed]
     [Google Scholar]
  29. Wick RR, Schultz MB, Zobel J, Holt KE. Bandage: interactive visualization of de novo genome assemblies. Bioinformatics 2015; 31:3350–3352 [View Article] [PubMed]
     [Google Scholar]
  30. Oxford Nanopore Technologies Medaka. n.d https://github.com/nanoporetech/medaka
  31. Wick RR, Holt KE. Polypolish: short-read polishing of long-read bacterial genome assemblies. PLoS Comput Biol 2022; 18:e1009802 [View Article] [PubMed]
     [Google Scholar]
  32. Hu J, Fan J, Sun Z, Liu S. NextPolish: a fast and efficient genome polishing tool for long-read assembly. Bioinformatics 2020; 36:2253–2255 [View Article]
     [Google Scholar]
  33. Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. 13033997 arXiv preprint arXiv 2013
     [Google Scholar]
  34. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J et al. The sequence alignment/map format and SAMtools. bioinformatics 2009; 25:2078–2079 [View Article]
     [Google Scholar]
  35. Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P et al. Artemis: sequence visualization and annotation. Bioinformatics 2000; 16:944–945 [View Article] [PubMed]
     [Google Scholar]
  36. Carver TJ, Rutherford KM, Berriman M, Rajandream M-A, Barrell BG et al. ACT: the artemis comparison tool. Bioinformatics 2005; 21:3422–3423 [View Article] [PubMed]
     [Google Scholar]
  37. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 2014; 30:2068–2069 [View Article] [PubMed]
     [Google Scholar]
  38. Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE et al. Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci USA 1998; 95:3140–3145 [View Article] [PubMed]
     [Google Scholar]
  39. Joensen KG, Tetzschner AM, Iguchi A, Aarestrup FM, Scheutz F. Rapid and easy in silico serotyping of Escherichia coli isolates by use of whole-genome sequencing data. J Clin Microbiol 2015; 53:2410–2426 [View Article] [PubMed]
     [Google Scholar]
  40. Lee I, Ouk Kim Y, Park S-C, Chun J. OrthoANI: an improved algorithm and software for calculating average nucleotide identity. Int J Syst Evol Microbiol 2016; 66:1100–1103 [View Article] [PubMed]
     [Google Scholar]
  41. Carattoli A, Zankari E, García-Fernández A, Voldby Larsen M, Lund O et al. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother 2014; 58:3895–3903 [View Article] [PubMed]
     [Google Scholar]
  42. Bortolaia V, Kaas RS, Ruppe E, Roberts MC, Schwarz S et al. ResFinder 4.0 for predictions of phenotypes from genotypes. J Antimicrob Chemother 2020; 75:3491–3500 [View Article] [PubMed]
     [Google Scholar]
  43. Alcock BP, Huynh W, Chalil R, Smith KW, Raphenya AR et al. CARD 2023: expanded curation, support for machine learning, and resistome prediction at the comprehensive antibiotic resistance database. Nucleic Acids Res 2023; 51:D690–D699 [View Article] [PubMed]
     [Google Scholar]
  44. Deatherage DE, Barrick JE. Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseq. Methods Mol Biol 2014; 1151:165–188 [View Article] [PubMed]
     [Google Scholar]
  45. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001; 25:402–408 [View Article] [PubMed]
     [Google Scholar]
  46. Edwards T, Williams C, Teethaisong Y, Sealey J, Sasaki S et al. A highly multiplexed melt-curve assay for detecting the most prevalent carbapenemase, ESBL, and AmpC genes. Diagn Microbiol Infect Dis 2020; 97:115076 [View Article] [PubMed]
     [Google Scholar]
  47. Edwards T, Sasaki S, Williams C, Hobbs G, Feasey NA et al. Speciation of common gram-negative pathogens using a highly multiplexed high resolution melt curve assay. Sci Rep 2018; 8:1114 [View Article] [PubMed]
     [Google Scholar]
  48. R Core Team R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; n.d https://www.R-project.org
  49. Villanueva RAM, Chen ZJ. ggplot2: elegant graphics for data analysis. Taylor & Francis; 2019
  50. Wilkins D. Gggenes: draw gene arrow maps in “ggplot2”. R package version 0.5.0; 2023
  51. Team RC. R: a language and environment for statistical computing; 2013
  52. Zhou K, Tao Y, Han L, Ni Y, Sun J. Piperacillin-tazobactam (TZP) resistance in Escherichia coli due to hyperproduction of TEM-1 β-lactamase mediated by the promoter Pa/Pb. Front Microbiol 2019; 10:833 [View Article] [PubMed]
     [Google Scholar]
  53. Römling U, Galperin MY, Gomelsky M. Cyclic di-GMP: the first 25 years of a universal bacterial second messenger. Microbiol Mol Biol Rev 2013; 77:1–52 [View Article] [PubMed]
     [Google Scholar]
  54. Tuckerman JR, Gonzalez G, Sousa EHS, Wan X, Saito JA et al. An oxygen-sensing diguanylate cyclase and phosphodiesterase couple for c-di-GMP control. Biochemistry 2009; 48:9764–9774 [View Article]
     [Google Scholar]
  55. Cimdins A, Simm R, Li F, Lüthje P, Thorell K et al. Alterations of c-di-GMP turnover proteins modulate semi-constitutive rdar biofilm formation in commensal and uropathogenic Escherichia coli. Microbiologyopen 2017; 6:e00508 [View Article] [PubMed]
     [Google Scholar]
  56. Bokranz W, Wang X, Tschäpe H, Römling U. Expression of cellulose and curli fimbriae by Escherichia coli isolated from the gastrointestinal tract. J Med Microbiol 2005; 54:1171–1182 [View Article] [PubMed]
     [Google Scholar]
  57. Lajhar SA, Brownlie J, Barlow R. Characterization of biofilm-forming capacity and resistance to sanitizers of a range of E. coli O26 pathotypes from clinical cases and cattle in Australia. BMC Microbiol 2018; 18:41 [View Article] [PubMed]
     [Google Scholar]
  58. Ross P, Weinhouse H, Aloni Y, Michaeli D, Weinberger-Ohana P et al. Regulation of cellulose synthesis in Acetobacter xylinum by cyclic diguanylic acid. Nature 1987; 325:279–281 [View Article] [PubMed]
     [Google Scholar]
  59. Richter AM, Possling A, Malysheva N, Yousef KP, Herbst S et al. Local c-di-GMP signaling in the control of synthesis of the E. coli biofilm exopolysaccharide pEtN cellulose. J Mol Biol 2020; 432:4576–4595 [View Article] [PubMed]
     [Google Scholar]
  60. Quail M, Haydon D, Guest J. The pdhR-aceEF-lpd operon of Escherichia coli expresses the pyruvate dehydrogenase complex. Mol Microbiol 1994; 12:95–104 [View Article] [PubMed]
     [Google Scholar]
  61. Steiert P, Stauffer L, Stauffer G. The lpd gene product functions as the L protein in the Escherichia coli glycine cleavage enzyme system. J Bacteriol 1990; 172:6142–6144 [View Article] [PubMed]
     [Google Scholar]
  62. Kim Y, Ingram L, Shanmugam K. Dihydrolipoamide dehydrogenase mutation alters the NADH sensitivity of pyruvate dehydrogenase complex of Escherichia coli K-12. J Bacteriol 2008; 190:3851–3858 [View Article] [PubMed]
     [Google Scholar]
  63. Harmer CJ, Moran RA, Hall RM. Movement of IS26-associated antibiotic resistance genes occurs via a translocatable unit that includes a single IS26 and preferentially inserts adjacent to another IS26. mBio 2014; 5:e01801-14 [View Article] [PubMed]
     [Google Scholar]
  64. Pedró L, Baños RC, Aznar S, Madrid C, Balsalobre C et al. Antibiotics shaping bacterial genome: deletion of an IS91 flanked virulence determinant upon exposure to subinhibitory antibiotic concentrations. PLoS One 2011; 6:e27606 [View Article] [PubMed]
     [Google Scholar]
  65. Zienkiewicz M, Kern-Zdanowicz I, Carattoli A, Gniadkowski M, Cegłowski P. Tandem multiplication of the IS26-flanked amplicon with the bla(SHV-5) gene within plasmid p1658/97. FEMS Microbiol Lett 2013; 341:27–36 [View Article] [PubMed]
     [Google Scholar]
  66. Moran RA, Baomo L, Doughty EL, Guo Y, Ba X et al. Extended-spectrum β-lactamase genes traverse the Escherichia coli populations of intensive care unit patients, staff, and environment. Microbiol Spectr 2023; 11:e0507422 [View Article] [PubMed]
     [Google Scholar]
  67. Shropshire WC, Aitken SL, Pifer R, Kim J, Bhatti MM et al. IS26-mediated amplification of blaOXA-1 and blaCTX-M-15 with concurrent outer membrane porin disruption associated with de novo carbapenem resistance in a recurrent bacteraemia cohort. J Antimicrob Chemother 2021; 76:385–395 [View Article] [PubMed]
     [Google Scholar]
  68. Chua SL, Sivakumar K, Rybtke M, Yuan M, Andersen JB et al. C-di-GMP regulates Pseudomonas aeruginosa stress response to tellurite during both planktonic and biofilm modes of growth. Sci Rep 2015; 5:10052 [View Article] [PubMed]
     [Google Scholar]
  69. Kwan BW, Osbourne DO, Hu Y, Benedik MJ, Wood TK. Phosphodiesterase DosP increases persistence by reducing cAMP which reduces the signal indole. Biotechnol Bioeng 2015; 112:588–600 [View Article] [PubMed]
     [Google Scholar]
  70. White A, Surette M. Comparative genetics of the rdar morphotype in Salmonella. J Bacteriol 2006; 188:8395–8406 [View Article] [PubMed]
     [Google Scholar]
  71. Liu C, Dong N, Huang X, Huang Z, Cai C et al. Emergence of the clinical rdar morphotype carbapenem-resistant and hypervirulent Klebsiella pneumoniae with enhanced adaption to hospital environment. Sci Total Environ 2023; 889:164302 [View Article]
     [Google Scholar]
  72. White A, Gibson D, Kim W, Kay W, Surette M. Thin aggregative fimbriae and cellulose enhance long-term survival and persistence of Salmonella. J Bacteriol 2006; 188:3219–3227 [View Article] [PubMed]
     [Google Scholar]
  73. Shropshire WC, Amiji H, Bremer J, Anand SS, Strope B et al. Genetic determinants underlying the progressive phenotype of β-lactam/β-lactamase inhibitor resistance in Escherichia coli. bioRxiv 2023 [View Article]
     [Google Scholar]
  74. Jiang J, Chen L, Chen X, Li P, Xu X et al. Carbapenemase-encoding gene copy number estimator (CCNE): a tool for carbapenemase gene copy number estimation. Microbiol Spectr 2022; 10:e0100022 [View Article] [PubMed]
     [Google Scholar]
  75. Rodríguez Villodres A, Gálvez Benítez L, Arroyo MJ, Méndez G, Mancera L et al. Ultrasensitive and rapid identification of ESRI developer- and piperacillin/tazobactam-resistant Escherichia coli by the MALDIpiptaz test. Emerg Microbes Infect 2022; 11:2034–2044 [View Article] [PubMed]
     [Google Scholar]
/content/journal/jmm/10.1099/jmm.0.002018
Loading
A high-resolution genomic and phenotypic analysis of resistance evolution of an Escherichia coli strain from a critically unwell patient treated with piperacillin/tazobactam
J. Med. Microbiol. 74, 002018 (2025); https://doi.org/10.1099/jmm.0.002018
/content/journal/jmm/10.1099/jmm.0.002018
/content/journal/jmm/10.1099/jmm.0.002018
Loading

Data & Media loading...

Supplements

Supplementary material 1

PDF

Most cited Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An error occurred
Approval was partially successful, following selected items could not be processed due to error