Clinical epidemiology
Clinical epidemiology analysis was performed using 695 Kp infections reported from 2017 to 2022 at Sapporo Medical University Hospital, including 393 “Colonization” and 302 “Infection” cases. “Infection” cases were analysed by classifying them into the presence or absence of immunosuppression, site of infection, and presence or absence of bacteraemia. Contingency tables were analysed using Fisher’s exact test. A p value < 0.05 was considered to indicate significance.
Bacterial isolation, antimicrobial susceptibility testing, and string test
A total of 277 Kp strains were isolated from clinical specimens from hospitalized patients at Sapporo Medical University between 2017 and 2021. These clinical specimens comprised 100 urine samples, 113 respiratory samples, 12 blood samples, and 52 other samples (drainage, tongue coating, skin, vaginal lubricant, pus, and bile). Identification of Kp (K. pneumoniae subsp.) was performed by MALDI Biotyper (Bruker Corporation, Billerica, USA). BIDMC_1, a carbapenem-resistant Kp strain isolated at the Beth Israel Deaconess Medical Center (BIDMC), was provided by BEI Resources (NIAID, NIH, USA).
The antimicrobial susceptibility of Kp strains was tested by the broth microdilution method, and the results were interpreted according to Clinical and Laboratory Standards Institute (CLSI) recommendations. In this study, the following antimicrobial agents were used: ciprofloxacin (Wako Pure Chemical Industry, Tokyo, Japan), ciprofloxacin hydrochloride monohydrate (Tokyo Chemical Industry, Tokyo, Japan), amikacin (Wako Pure Chemical Industry, Tokyo, Japan), kanamycin (Wako Pure Chemical Industry, Tokyo, Japan), and meropenem (Wako Pure Chemical Industry, Tokyo, Japan).
HMV strains were defined by a positive string test as previously described1. A single colony grown overnight on Mueller-Hinton II agar was fished, and the formation of a string > 5 mm in length was defined as a positive result. For detection of hypervirulence factors (serotypes K1 and K2, rmpA, rmpA2, iutA, iroN, and an IncHIB plasmid), multiplex PCR was performed as previously described2.
Serum susceptibility
In this study, we used human serum from individual healthy donors (Cedarlane Laboratories Ltd, Burlington, Canada). The serum MIC was defined as the minimum % serum concentration that prevented the visible growth of microorganisms. We set the resistance breakpoint at 32%, and isolates that exhibited more than 48% of serum MIC were defined as serum-resistant isolates because this concentration is the composition of the total blood in humans. Kp strains were grown in 0.5 ml of tryptic soy broth from an overnight culture. The strains were diluted 10-4-fold (105 CFU/mL) and incubated in plates with different serum concentrations for each well. After 20 h, colonies were visually confirmed because wells with high serum concentrations had high optical density (OD600 nm).
Measurement of mutation frequency
Mutation frequency was measured by rifampicin assay3. The Kp isolates were cultured overnight in tryptic soy broth. The solution was concentrated 10-fold and plated onto plain or 100 mg/L rifampicin-containing Mueller Hinton II agar plates, and the plates were cultured at 37°C for 24 h. After cultivation, the number of colony-forming units (CFU) that grew on the agar plates was counted. The gene mutation frequency was calculated as [CFU on the rifampicin-containing MH agar plate]/[CFU on the plain MH agar plate]. We defined mutator types as hyper (> 10-7), high (from 10-8 to 10-7), moderate (from 10-9 to 10-8), and low (<10-9). Student’s t test was used for the statistical analysis. A p value < 0.05 was considered to indicate significance.
Serial passaging experiments
Serial passaging experiments were performed by incubating Kp isolates (SMKP838, SMKP590, and BIDMC) in 96-well plates with MHBII containing certain concentrations (serial dilutions from original concentrations) of human serum or antimicrobial agents (ciprofloxacin, amikacin, and meropenem), as previously described4. For the experiments using other Kp clinical isolates, we selected 19 serum-susceptible Kp isolates (serum MICs were from 8 to 16%) from among the hyper- (n = 1), high- (n = 9; contains one K. quasipneumoniae), and low-mutators (n= 9) in the serial passaging experiment in the presence of human serum. In the ciprofloxacin assay, we selected 22 ciprofloxacin-susceptible Kp clinical isolates (ciprofloxacin MICs were from 0.03 to 0.25 mg/L) from among the hyper- (n = 1), high- (n = 11; contains one K. quasipneumoniae), and low-mutators (n = 10). We picked the well with the highest concentration (sub-MIC) of human serum or antimicrobial agent in which the bacteria grew and diluted the bacterial culture 100-fold with 0.85% NaCl. Then, 1 µL of the dilution was inoculated in 96-well plates containing 100 µL of MHBII with various concentrations of human serum or antimicrobial agents and cultivated at 37°C for 24 h. This serial passaging was repeated for 20 days in triplicate.
Time-killing assay
Single colonies of SMKP838 and the mutS mutant strains were grown overnight in TSB medium. The culture solution was adjusted to a final concentration of 1 × 105 CFU/mL and incubated for 0-24 h with each serum-containing solution (1/4 × MIC, 1 × MIC, or 2 × MIC) or in solution without serum (Ctl) at 37°C without shaking. The assay result was determined at 0 min, 30 min, 1 h, 3 h, 6 h, and 24 h.
WGS
Genomic DNA was isolated by a DNeasy Blood & Tissue Kit (Qiagen, Hulsterweg, The Netherlands). The DNA library was prepared by a Nextera XT DNA Library Preparation Kit (Illumina, San Diego, CA) for sequencing 300 bp paired-end reads according to the manufacturer’s protocol. An Illumina MiSeq was used for WGS. CTX-M genes were identified using assembled genome data by Resfinder (https://www.genomicepidemiology.org). MLST was performed using the Institute Pasteur MLST database and software (https://bigsdb.pasteur.fr/klebsiella/). Fast average nucleotide identification (FastANI) against the type strain genome was utilized for species identification. Core-genome single-nucleotide polymorphism (SNP)-based phylogenetic analysis was conducted: the Kp ATCC 35657 genome (accession number: CP015134.1) was used as a mapping reference. Mapping and core-genome extraction were performed using BWA version 0.7.17 with the “bwasw” option, SAMtools version 1.6 with the “mpileup” option, and VerScan version 2.3.9 with the “mpileupcns” option. The exclusion of estimated homologous recombination regions was performed using ClonalFrameML version v1.11-2. Snp-dists was used to count the pairwise SNP distance. A phylogenetic tree was generated using FastTree version 2.1.11 and FigTree version 1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/). The number of accumulated gene mutations during serial passaging experiments was analysed by mapping the genome reads to the reference genome (wild-type strain on day 0) obtained from WGS, followed by basic variant detection using CLC Genomics Workbench 21 (QIAGEN).
Bacterial growth determination
Bacterial growth was monitored by measuring the turbidity (that is, the optical density at 600 nm [OD600]) using an Infinite M200 PRO multimode microplate reader (Tecan, Kawasaki, Japan). Strains were grown in 0.5 ml of TSB (Becton Dickinson) overnight at 37°C, and 1 × 105 CFU/ml bacteria were cultured in 0.1 ml of MHBII broth (Becton Dickinson) in a 96-well plate at 37°C with shaking at 140 rpm for 16 h. Bacterial growth curves were created based on measurements every 10 min for 16 h.
Transposon-directed insertion site sequencing (TraDIS)
The SMKP838 transposon library was constructed using the EZ-Tn5™ <KAN-2> Tnp Transposome™ Kit (Epicentre, Wisconsin, USA). The bacteria with transposase introduced by electroporation (2.5 kV/cm, 200 Ω, and 25 μF) were selected by the formation of colonies on MHII agar containing 50 mg/l kanamycin. Over 100,000 colonies were collected, pooled, and frozen at -80°C in TSB with 10% glycerol as stock solutions until use. The transposon mutant library (106 cfu/mL) was inoculated into 1 mL of plain MHBII, MHBII containing 4% or 8% serum, or MHBII containing 40 mg/L surfactant protein A (SPA) and cultured at 37°C for 20 h. Total DNA was isolated using a Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA). Total DNA (500 ng) was used to prepare the DNA library for TraDIS using an NEBNext Ultra II FS DNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, USA). After fragmentation, end repair, 5' phosphorylation, dA-tailing, and adaptor ligation, and size selection (275-475 bp) according to the manufacturer’s protocol, the transposon-inserted genes were amplified by PCR using NEBNext Ultra II Q5 Master Mix (New England Biolabs), 20 nM NEBTnF2fas (5'-TCGACCTGCAGGCATGCAAGCTTCAGGGTTGAGATGTG-3') and NEBTn5-700 (5'-GTGACTGGAGTTCAGACGTGTGCTCTTCCGATC-3') primers, and 20 ng of fragmented DNA as the template with following condition: initial denature at 98°C for 30 sec, 22 cycles of 98°C for 10 sec and 72°C for 1 min 15 sec, and final extension at 72°C for 2 min. After the purification of the PCR product using AMPure XP beads (Beckman Coulter, Brea, CA, USA), enrichment PCR was performed by using a KAPA HiFi HotStart Library Amplification Kit (Roche, Basel, Switzerland), 20 nM NEBNext i700 primers including NEBNext Multiplex Oligos for Illumina (New England Biolabs) and NEBTn5-501-3 (5'- AATGATACGGCGACCACCGAGATCTACACTATAGCCTACACTCTTTCCCTACACGACGCTCTTCCGATCTTCGACCTGCAGGCATGCAAGCTTC-3'), and 20 ng of the purified DNA as the template with following condition: initial denature at 98°C for 45 sec, 10 cycles of 98°C for 15 sec, 60°C for 30 sec, and 72°C for 10 sec, and final extension at 72°C for 30 sec. The PCR products were purified and size selected (average: 650 bp) using AMPure XP beads. These products were pooled, and a NovaSeq600 was used for TraDIS. TraDIS analysis was performed according to a previous study5, and a false discovery rate-adjusted p value (FDRp) < 0.05 (vs. plain medium) was defined as significant.
Genes with significantly lower detected levels in serum or SPA samples (> 2-fold vs. plain medium) were considered putative serum or SPA resistance genes.
Construction of gene deletion mutants
ThemutS-deletion SMKP838 mutant and each serum resistance gene were generated by the λ-Red recombinase system, as previously described, using pKD46-hyg6,7. Each gene was replaced with Mini genes containing kanamycin resistance cassettes (Gene Bridges, Heidelberg, Germany) and 50 bases corresponding to the upstream and downstream regions of the target genes. The gene deletions were confirmed by PCR using specific primers.
Mouse models of lung and bloodstream infection
Ten- to 12-week-old female BALB/c mice were anaesthetized and infected transbronchially with a microsprayer (TORAY PRECISION, Tokyo, Japan) with 50 μl of a 1×108 CFU/ml solution. The mice were immunosuppressed by intraperitoneal injection of 250 mg/kg five days prior to infection and 125 mg/kg one day prior to infection with cyclophosphamide monohydrate (lot No. SKE6784; FUJIFILM Wako Pure Chemical Corporation, Osaka, Japan). In the treatment group, mice were injected subcutaneously with 100 mg/kg ciprofloxacin monohydrochloride (TOKYO CHEMICAL INDUSTRY, Tokyo, Japan) 1 and 24 hours after infection. After 32 hours in the nontreated group and 48 hours in the treated group or when a ‘prelethal critical endpoint’ had been reached, mice were euthanized by cervical dislocation. Lungs were washed with sterile PBS and homogenized with a gentleMACS Dissociator (Miltenyi Biotec). Homogenates were plated for determination of the number of CFU per lung. Blood was collected by puncturing the jugular vein. A drop of blood was added to 1 ml of PBS, vortexed, and then cultured on an appropriate plate. In the experiments, wild type-derived strains were selected by seeding on MHII agar containing 100 mg/l ampicillin sodium and ΔmutS-derived strains on MHII agar containing 50 mg/l kanamycin to eliminate the effects of other indigenous and environmental bacteria. Fifty colonies from each specimen were randomly selected, and their serum MIC was measured.
Construction of mutS-mutated BIDMC_1
Since the λ-Red recombination system was unable to generate mutS-deficient strains of BIDMC1, we used pORTMAGE and constructed mutS-mutated BIDMC1 (BIDMC1 MutS_Tyr37STOP)8. Hygromycin-integrated pORTMAGE was generated as previously described9. Transformants were produced by electroporation. Oligonucleotides (90 bp) for mutS containing the C111T mutation were designed using the MEGA Oligo Design Tool: MegamutS, CGCAACATCCTGACATTCTGCTGTTTTACCGGATGGGGGATTTTTAaGAGCTATTTTATGACGATGCGAAACGCGCCTCGCAGCTGCTCG; bold “a” indicates an introduced base. Gene replacement was confirmed by direct DNA sequencing.
Ethics statement
This study was approved by the Sapporo Medical University Hospital Institutional Review Board (IRB No. 272-70) and Sapporo Medical University Animal Care and Use Committee (Nos. 17-137, 18-083, and 20-006).
Statistical analysis
We used Prism 9 to calculate the significance of differences. Unpaired, two-tailed Student’s t test or a two-tailed Mann‒Whitney U test was used to compare two groups, and Dunn’s comparison test followed by the Kruskal–Wallis test was used to compare three or more groups. In addition, the log-rank test was used for survival data analysis. Statistical methods and P values are described in each figure. A P value < 0.05 was considered to indicate significance. In addition, Fisher’s exact test and Student’s t test were used in the clinical analysis to compare two groups.
Methods references
1. Vila, A. et al. Appearance of Klebsiella pneumoniae liver abscess syndrome in Argentina: case report and review of molecular mechanisms of pathogenesis. Open Microbiol. J. 5, 107-113 (2011).
2 Yu, F. et al. Multiplex PCR Analysis for Rapid Detection of Klebsiella pneumoniae Carbapenem-Resistant (Sequence Type 258 [ST258] and ST11) and Hypervirulent (ST23, ST65, ST86, and ST375) Strains. J Clin Microbiol 56, e00731-18 (2018).
3. Zhou, H. et al. The mismatch repair system (mutS and mutL) in Acinetobacter baylyi ADP1. BMC Microbiol. 20, 40 (2020).
4. Khil, P. P. et al. Dynamic emergence of mismatch repair deficiency facilitates rapid evolution of ceftazidime-avibactam resistance in Pseudomonas aeruginosa acute infection. mBio 10, e01822-19 (2019).
5. Dorman, M. J., Feltwell, T., Goulding, D. A., Parkhill, J. & Short, F. L. The Capsule Regulatory Network of Klebsiella pneumoniae Defined by density-TraDISort. mBio 9, e01863-18, (2018).
6. Datsenko, K. A. & Wanner, B. L. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc. Natl. Acad. Sci. U. S. A. 97, 6640-6645 (2000).
7. Sato, T. et al. Tigecycline nonsusceptibility occurs exclusively in fluoroquinolone-resistant Escherichia coli clinical isolates, including the major multidrug-resistant lineages O25b:H4-ST131-H30R and O1-ST648. Antimicrob. Agents Chemother. 61, e01654-16 (2017).
8. Nyerges, Á. et al. A highly precise and portable genome engineering method allows comparison of mutational effects across bacterial species. Proc. Natl. Acad. Sci. U. S. A. 113, 2502-2507 (2016).
9. Sato, T. et al. Emergence of the Novel Aminoglycoside Acetyltransferase Variant aac(6')-Ib-D179Y and Acquisition of Colistin Heteroresistance in Carbapenem-Resistant Klebsiella pneumoniae Due to a Disrupting Mutation in the DNA Repair Enzyme MutS. mBio 11, e01954-20 (2020).