Front. Public Health Frontiers in Public Health Front. Public Health 2296-2565 Frontiers Media S.A. 10.3389/fpubh.2026.1736106 Original Research Circulation of ESBL-producing Klebsiella pneumoniae at the human–animal–environment interface in farms on the central coast of Peru Benavides Julio A. 1 2 * Methodology Conceptualization Validation Writing – review & editing Supervision Resources Formal analysis Visualization Writing – original draft Funding acquisition Project administration Salgado-Caxito Marília 1 Investigation Visualization Data curation Formal analysis Writing – review & editing Methodology Writing – original draft Escobar-Páramo Patricia 1 Data curation Formal analysis Visualization Writing – review & editing Writing – original draft Léon Daphne 1 3 Data curation Writing – review & editing Investigation Writing – original draft Jara Luis M. 3 Investigation Data curation Formal analysis Writing – review & editing Bardales-Mendoza Olga 4 Writing – review & editing Investigation Murga Clara 3 Investigation Writing – review & editing Medrano Patricia 3 Investigation Writing – review & editing Pérez Veronica 3 5 Writing – review & editing Investigation Aylas-Jurado Brenda 3 Writing – review & editing Investigation Su-Tello Roberto 3 Investigation Writing – review & editing Najarro Juana 3 Investigation Writing – review & editing Osorio-Linares Raiza 6 Writing – review & editing Investigation Salvador-Tasayco Elías 6 Writing – review & editing Investigation Shiva Carlos 3 * Formal analysis Writing – original draft Methodology Supervision Validation Funding acquisition Conceptualization Writing – review & editing Resources Project administration UMR MIVEGEC, IRD, CNRS, University of Montpellier, Montpellier, France One Health Institute and Doctorado en Medicina de la Conservación, Faculty of Life Sciences, Universidad Andrés Bello, Santiago, Chile Facultad de Medicina Veterinaria y Zootecnia, Universidad Peruana Cayetano Heredia, Lima, Peru Facultad de Educación, Universidad Peruana Cayetano Heredia, Lima, Peru Emerge, Unidad de Investigación en Enfermedades Emergentes y Cambio Climático, Facultad de Salud Pública y Administración, Universidad Peruana Cayetano Heredia, Lima, Peru Facultad de Medicina Veterinaria y Zootecnia, Universidad Nacional San Luis Gonzaga, Ica, Peru Correspondence: Julio A. Benavides, julio.benavides@ird.fr; Carlos Shiva, carlos.shiva@upch.pe 05 02 2026 2026 14 1736106 30 10 2025 26 12 2025 02 01 2026 Copyright © 2026 Benavides, Salgado-Caxito, Escobar-Páramo, Léon, Jara, Bardales-Mendoza, Murga, Medrano, Pérez, Aylas-Jurado, Su-Tello, Najarro, Osorio-Linares, Salvador-Tasayco and Shiva. 2026 Benavides, Salgado-Caxito, Escobar-Páramo, Léon, Jara, Bardales-Mendoza, Murga, Medrano, Pérez, Aylas-Jurado, Su-Tello, Najarro, Osorio-Linares, Salvador-Tasayco and Shiva https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Antibiotic-resistant bacteria, such as Extended-spectrum β-lactamase-producing Klebsiella pneumoniae (ESBL-Kp), represent a major threat to public health. Despite increasing reports of community-acquired ESBL-Kp infections, the dynamics of ESBL-Kp in low- and middle-income countries communities remains poorly understood. This study aimed to detect the fecal carriage of ESBL-Kp at the human–animal–environment interface in Peru and to characterize underlying molecular mechanisms involved. We detected 15 ESBL-Kp isolates among 652 (2.3%) fecal and water samples (i.e., 244 dairy cows from 25 farms, 261 pigs from 25 farms, 58 dogs, 39 farmers, and 50 water sources). ESBL-Kp was isolated from two humans, two dogs, four cows, three pigs, and four water sources. Genomic analyses identified 10 different ESBL-Kp sequence types (STs), including high-risk clones ST11, ST307, and ST37, as well as ST45, ST111, ST1, ST348, ST789, ST881, and ST983, and three CTX-M-encoding genes (blaCTX-M-15, blaCTX-M-27, blaCTX-M-14). Four ESBL-Kp STs (ST348, ST983, ST789, and ST11) were found in more than one source, both within and between farms. In particular, two ESBL-Kp ST983 isolates, one from a human and one from a cow on the same farm, differed by 37 SNPs and had almost identical genetic profiles, suggesting bacterial clonal exchange across host species or contamination from a common source. Likewise, two nearly identical ST348 isolates were recovered from a cow and a water source on the same farm, suggesting ESBL-Kp exchange between livestock and water. Our results highlight the circulation of ESBL-Kp across humans, animals, and water in rural environments in Peru, representing one of the first genomic studies exploring potential cross-species ESBL-Kp relatedness in Latin America. Our study supports the need to implement sanitary measures to limit the cross-species spread of antimicrobial-resistant bacteria and to reduce water source contamination in Peru.

 antimicrobial resistance domestic animal extended-spectrum β-lactamase farmer Latin America livestock One Health water Institut de Recherche pour le Développement 10.13039/100012947 PE501078124-2022 The author(s) declared that financial support was received for this work and/or its publication. This research was funded by PROCIENCIA-CONCYTEC (grant number PE501078124-2022) and by the Jeune Équipe Associée à l’Institut de Recherche pour le Développement-IRD (JEAI) STOPAMR: ‘Stratégies One Health pour la prévention de l’antibiorésistance en milieu rural’. MS-C is funded by a postdoctoral scholarship awarded by the IRD. DL is funded by a PhD scholarship from the ‘Escuela Doctoral Franco Peruana en Ciencias de la Vida’ of the Universidad Peruana Cayetano Heredia and the IRD. section-at-acceptance Infectious Diseases: Epidemiology and Prevention Introduction

Extended-spectrum β-lactamase-producing Klebsiella pneumoniae (ESBL-Kp) represents an important threat to public health, and has been traditionally regarded as a hospital-associated pathogen (1). However, growing evidence suggests that ESBL-Kp can also be acquired in community settings, highlighting the need to understand its circulation within the community (2, 3). Klebsiella pneumoniae (Kp) is naturally present in the microbiome of humans and animals, where it can also act as an opportunistic nosocomial pathogen, and is estimated to cause approximately one-third of all Gram-negative bacterial infections (4). Kp can cause a wide range of diseases including pneumonia, urinary tract infections (UTIs), cystitis, bloodstream infections, and septicemia (5). The rise of multidrug-resistant Kp, particularly ESBL-Kp, poses a serious public health threat due to limited therapeutic options available, as well as its associated increases in mortality, hospital length of stay, and healthcare costs (5, 6).

ESBL-Kp circulates in hospital settings at variable levels across regions. The prevalence of hospital-acquired infections in high-income countries (HICs) ranges from 8% in Australia and New Zealand, 12% in the USA, 13% in Canada, and up to 21% in Western Europe (7). In contrast, low- and middle-income countries (LMICs) report higher prevalence, ranging from 23% in Asia excluding China, 32% in Eastern Europe, and up to 40% in the Middle East and Africa (7). Latin America has also reported high hospital ESBL-Kp prevalence ranging from 40 to 60% (8–10). ESBL-Kp also causes community-acquired infections, including pneumonia, UTIs, and bacteremia (11–13). Understanding the dynamics of ESBL-Kp in community settings is essential to mitigate its burden. For example, a high prevalence of ESBL among community-acquired Kp infections has been estimated in LMICs including Asian countries (e.g., 17% in bloodstream infections in China (14); 10–15% in UTIs in India and Iran (15, 16)). The estimated proportion of ESBL-Kp in community-acquired Kp infections in Latin America is high compared to other regions, ranging from 18% in community-acquired UTIs caused by Kp in Cuba (17), to 37% in community-acquired bloodstream infections caused by Kp in Brazil (18). Thus, understanding the dynamics of ESBL-Kp in Latin American communities can contribute to reducing its burden.

Several studies have reported the fecal carriage of ESBL-Kp in healthy human populations. For example, a prevalence of <1% of ESBL-Kp has been detected in healthy humans of HICs [e.g., Spain, France, Sweden, and Norway (19)]. In contrast, the limited data available from LMICs suggest higher level of community intestinal carriage of ESBL-Kp in African and Asian countries [e.g., India: 6% (20); Ethiopia: 6% (21); Nepal: 10% (22); Chad: 13% (23)], while there is a lack of epidemiological data on the human fecal carriage of ESBL-Kp in Latin America. Kp can also colonize animals (24–26) and persist or disseminate through water sources (27), underscoring the need for integrated One Health approaches to study the spread of ESBL-Kp in the community. However, the transmission dynamics of ESBL-Kp at the human–animal–environment interface are poorly understood (28). In rural areas of Latin America, inadequate sanitation and inappropriate antibiotic use are widespread (29, 30), potentially favoring the selection and dissemination of ESBL-Kp across the human–animal–environment interface.

Bacterial genomic analyses can contribute to identifying the genetic mechanisms driving ESBL-Kp circulation in the community (e.g., Kp clones, mobile genetic elements), as well as elucidating transmission pathways. A diverse pool of ESBL-Kp sequence types (STs) has been reported to carry different ESBL-encoding genes (e.g., CTX-M, SHV, TEM) worldwide (31). However, hospital-acquired ESBL-Kp infections are still dominated by high-risk clones, including ST307, ST15, ST147, ST405, and ST11, usually encoding blaCTX-M-15 on plasmids (e.g., IncFIB/IncFII plasmid replicons) (32, 33). Latin America also reflects a dominance of Kp high-risk clones, including ST11, ST147, and ST307 in human infections (34), particularly carrying the blaCTX-M-15 gene (35). However, less is known about the ESBL-Kp STs circulating in animals and environmental sources (36). In Peru, previous genomic surveillance in the southern region revealed that all multidrug-resistant Kp isolates recovered from three hospitals between 2022 and 2023 were ESBL producers, predominantly belonging to the high-risk pandemic clones ST147, ST37, and ST629, and carrying blaCTX-M, blaSHV, and blaTEM alleles (37). Furthermore, the emergence of ST307 carrying the blaCTX-M-15 gene was detected in a tertiary hospital (38), and ESBL-Kp cases were associated with neonatal sepsis (39, 40).

In 2017, we isolated the ESBL-Kp high-risk Kpclone ST307 from the feces of a pig on a farm located in coastal Peru, representing the first report of ESBL-Kp circulation among livestock in the country (41). Moreover, farms in this region are characterized by a high prevalence of fecal carriage of ESBL-producing Escherichia coli in both cows and pigs (42). However, to our knowledge, no study has simultaneously investigated the circulation of ESBL-Kp among humans, animals, and environmental sources in Peru, particularly in rural settings.

The aim of this study was therefore to assess the occurrence of ESBL-Kp across humans, domestic animals, and water samples from small-scale farms in rural Peru. We further performed whole-genome sequencing to characterize the genetic determinants of resistance and to evaluate evidence of cross-species transmission.

 Methodology Study regions and sampling

This study was part of a larger research project focusing on understanding the circulation of ESBL-producing Enterobacterales across the human–animal–environment interface in coastal Peru, which previously reported a high prevalence of ESBL-E. coli in pigs (San Bartolo district) and cows (Huaura district) on small-scale farms in Lima (42). In this study, we focused on assessing the occurrence of fecal carriage of ESBL-Kp in farmers and animals, as well as the ESBL-Kp contamination in water collected from these farms. Based on our previous studies on ESBL-Enterobacteriales [48% prevalence of ESBL-E. coli (43)], we assumed a prevalence of ESBL-Kp fecal carriage of 40%, a confidence level of 90%, and a population ranging between 100,000 and 1 million animals, resulting in a total sample size of 260 livestock animals. As we aimed to collect approximately 10 livestock samples per farm (i.e., 90% probability of detecting at least one positive animal assuming a conservative 10% prevalence), we sampled 50 farms. In those farms, we aimed to collect samples from all dogs, all humans working directly with livestock, and one water sample. Dairy farms (n = 25) were selected as follows: 14 farms were randomly selected from a list of 100 member farms belonging to a local dairy association (the ‘Asociación de Ganaderos de la Irrigación San Felipe’) in the Vegueta district of Huaura in Lima; one farm was selected through an ongoing collaboration with the Universidad Peruana Cayetano Heredia (UPCH) in the Lurin district; and 10 farms were randomly selected from a list provided by a veterinarian working with UPCH in the Cañete district. Pig farms (n = 25) were randomly selected from a list of 50 farm owners in the La Chutana district of San Bartolo in Lima. Sampling was conducted between May and July 2023. A total of 652 samples were collected, including 505 fecal samples from livestock (244 cows and 261 pigs), 58 fecal samples from dogs, 39 fecal swabs from farmers who voluntarily participated, and 50 water samples (500 mL) from livestock water sources. Farm locations and details of husbandry practices (e.g., herd size, biosecurity level, antibiotic use) are provided in Salgado-Caxito et al. (42).

 ESBL-<italic>Kp</italic> isolation and antimicrobial susceptibility test

Microbiological procedures to isolate ESBL-producing Enterobacterales from fecal samples were performed as described in Salgado-Caxito et al. (42). Briefly, each fecal sample was stored in Stuart transport medium, plated on MacConkey agar (Himedia®, Maharashtra, India) supplemented with 4 mg/L of cefotaxime, and incubated at 37 °C for 24–48 h. For water samples, 500 mL were filtered through a 0.45 μm pore-size membrane filter (Millipore, USA) (44). The filters were then processed using the same procedure as that used for fecal samples. From plates showing bacterial growth, we selected all morphotypes compatible with Klebsiella spp., which resulted in the selection of a single colony per plate.

Bacterial species identification was performed using the VITEK®2 Advanced Expert System™ (AES), according to the manufacturer’s instructions. Phenotypic ESBL production was assessed using the combination disk test with cefotaxime (30 μg) and ceftazidime (30 μg) disks, tested alone and in combination with clavulanic acid (30/10 μg), following Clinical and Laboratory Standards Institute (CLSI) guidelines. Klebsiella pneumoniae subsp. pneumoniae ATCC® 700,603™ was used as a positive control for ESBL production, and Escherichia coli ATCC® 25,922™ (Microbiologics, St. Cloud, MN, USA) was used as the negative control.

Susceptibility to antimicrobials was also evaluated using the VITEK®2 AES. The following β-lactam antibiotics were included in the study: ceftazidime, cefepime, piperacillin/tazobactam, ertapenem, and meropenem. We also tested the fluoroquinolone ciprofloxacin and the aminoglycoside amikacin, two antibiotic classes widely used in swine and cattle farms in Peru and other low-income countries (30, 45, 46). Minimum inhibitory concentrations (MICs) were determined using the VITEK®2 (software version 9.02) in Natural Resistance (NATR) mode, based on Clinical and Laboratory Standards Institute (CLSI) guidelines (47).

 Whole genome sequencing (WGS)

Isolates confirmed as ESBL-Kp were submitted for short-read (Illumina®) genome sequencing at MicrobesNG (Birmingham, UK). DNA extraction, library preparation, genome sequencing, read trimming, assembly, and quality control were carried out by MicrobesNG (Birmingham, UK). Briefly, 5–40 μL of bacterial culture preserved in DNA/RNA Shield (Zymo Research, USA) was mixed with an extraction buffer containing lysozyme and RNase A and incubated at 37 °C for 25 min. Then, Proteinase K and a second aliquot of RNase A were added, and the incubation was repeated at 65 °C for 5 min. A DNA library was prepared using the Nextera XT Library Preparation Kit (Illumina Inc.), and short-read sequencing was performed on the Illumina NovaSeq 6,000 (Illumina, San Diego, USA) using 2 × 250 bp paired-end reads, with a minimum coverage of 30×, following the manufacturer’s protocol. Short-reads were adapter-trimmed using Trimmomatic (version 0.30) (48). Trimmed short-reads were then de novo assembled using SPAdes (version 3.14.1) (49), and contigs were annotated using Prokka (version 1.11) (50). Genome quality assessment of assemblies was performed with QUAST (51) and BUSCO (version 5.3.2, Enterobacterales_odb10 database) (52). Species identification and contamination screening were performed using Kraken2 (version 2.1.1) (53). Complete information about the sequencing and assembly processes, including the MicrobesNG pipeline, procedures, and software used, is publicly available at https://microbesng.com.

In silico analyses of the ESBL-Kp genomes were conducted using the Bacterial and Viral Bioinformatics Resource Center (BV-BRC1) (54) and the Center for Genomic Epidemiology (CGE2). We used BV-BRC to obtain a comprehensive genome report, including assembly statistics (completeness and contamination), genome annotations (i.e., genes related to antimicrobial resistance (AMR), drug targets, transporters, and virulence factors), and data visualization (55). All ESBL-Kp genomes were considered of good quality, with 97.1–98.7% consistency and 100% completeness based on BV-BRC assembly statistics (56). Mean coverage depth was 62.8 × (range: 45.3x-127x). Genome sizes averaged 5.52 Mb (range: 5.3–5.74 Mb), with contig N50 values averaging 3.57 Mb (range: 1.61–7.66 Mb), and an average GC content of 57.1%.

The BV-BRC tool was also used to predict phenotypic antibiotic resistance based on genomic features. The CGE was used to detect plasmid replicon types [PlasmidFinder v2.0.1 (57, 58)] and to identify Kp sequence types (STs) [MLST v2.0.9; (56, 58–63)]. The CGE tool MobileElementFinder (MGE) v1.0.3 was used to identify whether ESBL-encoding genes were associated with plasmids or chromosomes (64). AMR genes were predicted using a 95% identity threshold and a minimum coverage of 80%, while the same identity threshold with a minimum coverage of 60% was applied for plasmid detection (65, 66). For virulence genes, predictions were performed using a 95% identity threshold, the BLOSUM80 matrix, E-value < 0.01, and >70% coverage (67). All databases were accessed in March 2024.

To investigate the phylogenetic relationships of the ESBL-Kp isolates, we constructed a core genome MLST (cgMLST) phylogeny using the NGPhylogeny ‘A la carte Workflows’ platform3 (68). Phylogenetic trees were inferred using a maximum-likelihood approach with Smart Model Selection (PhyML + SMS), applying default parameters and 100 bootstrap replicates (Figure 1). In addition, clonal circulation was accessed by building single-nucleotide polymorphism (SNP)-based phylogenetic trees for ESBL-Kp isolates of the same ST using CSI Phylogeny 1.4 from the CGE (69–73). The obtained phylogenetic tree (Figure 2) was built using the altered version of FastTree with default parameters and visualized using FigTree v1.4.4. Because CSI Phylogeny requires at least three genomes to generate a SNP-based phylogeny, publicly available genomes of the same ST from South America were included when only two isolates of a given ST were identified in our dataset. These genomes were retrived from the Klebsiella Pasteur MLST database4. Given potential biases related to assembly differences between our sequences and publicly available assemblies, we re-assembled all genomes included in the SNP-based phylogenetic analyses on the BV-BRC platform using publicly available reads. Metadata for each genome in the phylogenetic analyses (source, year, country) are provided in Figure 2.

Core genome MLST (cgMLST) phylogeny and AMR gene profile of the 15 ESBL-Klebsiella pneumoniae isolates in rural Peru. Left: CgMLST phylogeny. Isolates sharing the same sequence type (ST) and originating from the same farm are highlighted in bold and with identical colors. The cgMLST tree was built using the PhyML+SMS tool implemented in the NGPhylogeny.fr online platform. Node labels indicate bootstrap support values, and the scale bar represents the number of substitutions per site. Figure created using FigTree v1.4.4. Right: Heatmap indicating the presence (blue) or absence (gray) of antimicrobial resistant genes (x-axis) for each isolate in the phylogenetic tree. Figure created using R 4.2.1.

Phylogenetic tree showing relationships among ESBL-Kp from differen farms, animal type, and sequence type (ST). Icons depict animals or water sources. Heatmap on the right represents gene presence, with blue squares indicating gene occurrence across various samples labeled at the bottom. Scale bar signifies genetic distance.

SNP-based phylogeny and pairwise SNP differences of ESBL-Kp isolates sharing the same ST within and across farms. Left: SNP-based phylogenetic trees were built with FastTree on the CGE platform (CSIPhylogeny) with default parameters. Colored labels indicate farm of origin of each isolate. The scale bar represents the number of substitutions per site. Right: Pairwise SNP distance matrices generated from the SNP-based phylogenetic trees shown on the left. Public Latin American genomes were included when fewer than three isolates per ST were available and were retrieved from the BIGSdb Pasteur database (https://bigsdb.pasteur.fr). Figure created using FigTree v1.4.4. and Inkscape.

Phylogenetic trees and corresponding SNP matrices for four ESBL-KP sequence types isolated in this study: ST348, ST983, ST789, and ST11. Each tree shows relationships between samples sources, farms and the same ST previously isolated in Latin America. The matrices display SNP differences between ESBL-Kp.

We also performed a SNP-based phylogenetic analysis comparing the ESBL-Kp ST307 isolate detected in a pig in Peru in our previous study in 2017 (41) with a newly isolated ST307 from this study. We included two other ST307 genomes, one from a human UTI in a hospital in Brazil (2016) (ERR4822787) and one from a bloodstream infection in Colombia (2014) (ERR4795646), representing the geographically closest ST307 genomes available in the database that included raw sequencing reads. In all SNP-based phylogenies, the isolate with the largest assembled genome size was used as the reference. The resulting Newick files were visualized and edited using FigTree v1.4.45.

 Ethics approval statement

This project was approved by the human and animal ethics committees (Comité Institucional de Ética en Humanos and Comité Institucional de Ética para el Uso de Animales) of the Universidad Peruana Cayetano Heredia under protocols 591–49-22 (human ethics committee) and 050–12-22 (animal ethics committee). Written informed consent was obtained from all farmers involved in this study.

 Results

A total of 15 ESBL-Kp isolates were recovered from 652 samples screened (2.3%). ESBL-Kp isolates were obtained from all sampled sources, including two humans, two dogs, four cows, three pigs, three water samples from dairy farms, and one water sample from a swine farm (Table 1). Antimicrobial susceptibility tests showed that all ESBL-Kp isolates exhibited phenotypic resistance to ceftazidime and cefepime, and the majority were resistant to ampicillin/sulbactam (11/15) and ciprofloxacin (12/15). None of the isolates were resistant to ertapenem, meropenem, or amikacin. However, some isolates exhibited intermediate susceptibility to ampicillin/sulbactam (4/15), ertapenem (1/15), and ciprofloxacin (2/15).

Source, farm of origin, and genetic profile of the 15 ESBL-Kp isolates from livestock farms in Peru.

Source Farm ST AMR extrinsic genes Virulence genes Plasmid replicons
Water San Bartolo-G8 ST1 aac(6′)-Ib10, adeF, aph(3′)-Ia, aph(3″)-Ib, aph(6)-Id, blaCTX-M-15, blaOXA-1, blaOXA-900, blaSHV-1, blaTEM-1, catB3, dfrA14, fosA6, oqxA, qnrB1, sul2, tet(C) entD, iroE, iutA, rfbA, rfbB, rfbD, vipB/tssC IncFIB(K)
Cow Huaura-G4 ST11 aac(6′)-Ib10, aadA16, aph(3′)-Ia, arr-3, blaCTX-M-27, blaDHA-1, blaSHV-11, dfrA27, floR, fosA6, oqxA, oqxB, qnrB4, sul1, tet(D) iroE, iutA, rfbK1, vipB/tssC Col440I, IncFIB(K), IncFIB(pKPHS1), IncFII(K), repB(R1701)
Cow Huaura-G5 ST11 aac(6′)-Ib10, aadA11, adeF, aph(3′)-Ia, arr-3, blaCTX-M-27, blaDHA-1, blaSHV-11, dfrA27, floR, fosA6, oqxA, oqxB, qnrB4, sul1, tet(D) iroE, iutA, rfbK1, vipB/tssC Col(pHAD28), IncFIB(K), IncFIB(pKPHS1), IncFII(K), repB(R1701)
Dog Huaura-G5 ST11 aac(6′)-Ib10, aadA16, adeF, aph(3′)-Ia, arr-3, blaCTX-M-27, blaDHA-1, blaSHV-11, dfrA27, floR, fosA6, oqxA, oqxB, qnrB4, sul1, tet(D) iroE, iutA, rfbK1, vipB/tssC Col(pHAD28), IncFIB(K), IncFIB(pKPHS1), IncFII(K), repB(R1701)
Water Huaura-G4 ST37 aac(6′)-Ib10, aadA11, aadA2, adeF, aph(3′)-Ia, arr-3, blaCTX-M-27, blaSHV-11, dfrA12, dfrA27, floR, fosA6, oqxA, oqxB, qnrB2, sul1 entD, impA/tssA, iroE, iutA IncFIA(HI1), IncFIB(K)(pCAV1099-114), IncR, repB(R1701)
Human San Bartolo-G5 ST45 adeF, aph(3″)-Ib, aph(6)-Id, blaCTX-M-15, blaSHV-1, blaTEM-1, fosA6, oqxA, sul2 entD, fyuA, impA/tssA, iroE, irp1, irp2, iutA, rfbA, rfbB, rfbD, rfbK1, vipB/tssC, ybtA, ybtE, ybtP, ybtQ, ybtS, ybtT, ybtU, ybtX IncFIB(K)
Pig San Bartolo-G4 ST111 aadA, adeF, blaCTX-M-15, blaOXA-900, blaSHV-11, blaTEM-1, dfrA17, fosA6, oqxA, oqxB, sul1, tet(C) entD, fyuA, iroE, irp1, irp2, iutA, rfbA, rfbB, rfbD, vipB/tssC, ybtA, ybtE, ybtP, ybtQ, ybtS, ybtT, ybtU, ybtX IncFIB(K)
Pig San Bartolo-G4 ST307 aac(3)-IIc, aac(6′)-Ib10, aadA2, adeF, aph(3′)-Ia, aph(3″)-Ib, aph(6)-Id, blaCTX-M-15, blaOXA-1, blaSHV-28, blaTEM-1, catB3, dfrA12, floR, fosA6, oqxA, qnrB1, sul1, sul2, tet(A) iroE, iutA, rfbA, rfbB, rfbD, vipB/tssC IncFIB(K)
Cow Huaura-G3 ST348 aac(6′)-Ib10, adeF, aph(3″)-Ib, aph(6)-Id, blaCTX-M-15, blaOXA-1, blaSHV-11, blaTEM-1, catB3, dfrA14, fosA6, oqxA, oqxB, qnrB1, sul2, tet(C) fyuA, impA/tssA, iroE, irp1, irp2, iutA, rfbA, rfbB, rfbD, rfbK1, vgrG/tssI, vipB/tssC, ybtA, ybtE, ybtP, ybtQ, ybtS, ybtT, ybtU, ybtX IncFIB(K), IncFII(K)
Water Huaura-G3 ST348 aac(6′)-Ib10, adeF, aph(3″)-Ib, aph(6)-Id, blaCTX-M-15, blaOXA-1, blaSHV-11, blaTEM-1, catB3, dfrA14, fosA6, oqxA, qnrB1, sul2, tet(C) fyuA, impA/tssA, iroE, irp1, irp2, iutA, rfbA, rfbB, rfbD, rfbK1, vgrG/tssI, vipB/tssC, ybtA, ybtE, ybtP, ybtQ, ybtS, ybtT, ybtU, ybtX IncFIB(K), IncFII(K)
Pig San Bartolo-G3 ST789 aac(6′)-Ib10, aadA2, aph(3′)-Ia, blaCTX-M-15, blaOXA-1, blaOXA-900, blaSHV-25, blaTEM-1, catB3, dfrA12, fosA6, oqxA, oqxB, qnrB11, qnrB12, qnrB19, qnrB5, qnrB50, qnrB54, qnrB80, sul1, tet(A) entD, iroE, irp1, rfbA, rfbB, rfbD, vipB/tssC Col(pHAD28), IncFIB(K), IncFII(K)
Dog San Bartolo-G6 ST789 aac(6′)-Ib10, aadA2, adeF, aph(3′)-Ia, blaCTX-M-15, blaOXA-1, blaOXA-900, blaSHV-25, blaTEM-1, catB3, dfrA12, fosA6, oqxA, qnrB11, qnrB12, qnrB19, qnrB21, qnrB25, qnrB50, qnrB54, qnrB71, qnrB73, qnrB80, sul1, tet(A) entD, iroE, iutA, rfbA, rfbB, rfbD, vipB/tssC Col(pHAD28), IncFIB(K), IncFII(K)
Water Huaura-G7 ST881 aac(3)-IId, adeF, aph(3″)-Ib, aph(6)-IdV, blaCTX-M-14, blaSHV-27, fosA6, oqxA, qnrS1, sul1, sul2, tet(C) iroE, iutA, rfbA, rfbB, rfbD, rfbK1, vipB/tssC IncFIB(K)
Cow Huaura-G8 ST983 adeF, aph(3″)-Ib, aph(6)-Id, blaCTX-M-15, blaSHV-168, blaTEM-1, dfrA14, fosA6, oqxA, qnrB1, sul2, tet(C) entD, fyuA, iroE, irp1, irp2, iutA, vipB/tssC, ybtA, ybtE, ybtP, ybtQ, ybtS, ybtT, ybtU, ybtX IncFIB(K)
Human Huaura-G8 ST983 adeF, aph(3″)-Ib, aph(6)-Id, blaCTX-M-15, blaSHV-168, blaTEM-1, dfrA14, fosA6, oqxA, oqxB, qnrB1, sul2, tet(C) entD, fyuA, iroE, irp1, irp2, iutA, vipB/tssC, ybtA, ybtE, ybtP, ybtQ, ybtS, ybtT, ybtU, ybtX IncFIB(K)

Only extrinsic AMR genes are listed, and virulence genes are shown only when they are absent in at least one isolate (i.e., not present in 100% of isolates).

 Antimicrobial resistance genes

Genomic analyses identified 94 unique AMR genes. Of these, 28 (30%) were classified as intrinsic, 54 (57%) as extrinsic (acquired), and 12 (13%) as indeterminate, mostly corresponding to chromosomal mutations or resistance determinants not clearly attributable to horizontal transfer. The distribution of extrinsic AMR genes differed across sources. Dog isolates exhibited the highest number of extrinsic AMR genes per genome (16 and 25 genes), followed by pigs (12, 20, and 22 genes), cows (12 and 16 genes), and water-derived isolates (12 to 17 genes). In contrast, the two human isolates harbored only 9 and 13 genes. For the 7 antibiotics tested, the predicted resistance from the genomes using BV-BRC agreed with the phenotypic resistance obtained through the VITEK®2 AES system in 89.8% of all combinations (Supplementary Table S2).

β-lactam resistance genes detected included variants of blaCTX-M, blaSHV, blaTEM, blaOXA, and blaDHA, which together accounted for resistance to ceftazidime, cefepime, and ampicillin/sulbactam. ESBL-encoding genes were detected in all isolates including blaCTX-M-15 in 10 isolates, blaCTX-M-27 in four isolates, and blaCTX-M-14 in one isolate. Narrow-spectrum β-lactamases were also common, including blaTEM-1 (10 isolates), blaOXA-1 (6 isolates), and blaOXA-900 (4 isolates), as well as six distinct SHV alleles (blaSHV-11 in seven isolates, blaSHV-168 in two isolates, blaSHV-25 in two isolates, blaSHV-1 in two isolates, blaSHV-27 in one isolate, and blaSHV-28 in one isolate). The plasmid-mediated AmpC gene blaDHA-1 was detected in the three ST11 isolates. Common resistance genes to other antibiotic classes included aminoglycoside-modifying enzymes such as aadA, aph(3′), aph(6)-Id, and aac(6′)-Ib; plasmid-mediated quinolone resistance genes (qnrB); sulfonamides (sul1 or sul2); fosfomycin (FosA6); tetracycline (tet(A), tet(C), or tet(D)); trimethoprim (dfrA12, dfrA14, dfrA17, or dfrA27). In contrast, genes conferring resistance to phenicol (floR in 5 isolates) were less common. All AMR and virulence genes detected in these isolates are provided in Supplementary Table S1.

 Virulence-associated genes

A total of 71 non-redundant virulence-associated genes were identified across the 15 ESBL-Kp genomes. However, none of the ESBL-Kp isolates carried core hypervirulence regulators such as iucA, iroB, peg-344, rmpA, and rmpA2. Virulence genes detected were associated with siderophore-mediated iron acquisition systems (26 genes, 37%, e.g., iroE) and adhesion factors (19 genes, 27%, e.g., fimK and mrkI), followed by capsule/LPS-related loci (7 genes, 10%) and Type VI secretion system (T6SS) components (12 genes, 17%, e.g., tssG, vasE/tssK). Water-derived isolates carried 54–70 virulence-associated genes, cow isolates 54–70, pig isolates between 56 and 68, and dog isolates 54–57. Virulence gene profiles were largely similar across sampled sources. The presence of each AMR and virulence gene was used to compare isolates belonging to the same ST (Figure 1).

 Plasmid replicons

The PlasmidFinder tool predicted nine different plasmid replicons in the ESBL-Kp isolates, including Col(pHAD28), Col440I, IncFIA(HI1), IncFIB(K), IncFIB(K)(pCAV1099-114), IncFIB(pKPHS1), IncFII(K), IncR, and repB(R1701). The plasmid replicon IncFIB(K) was present in all genomes. IncFII(K) was also frequently identified (7/15) and occurred exclusively in co-occurrence with IncFIB(K). Smaller plasmids such as Col(pHAD28) and repB(R1701) were detected in four isolates, while IncR and Col440I were detected in only one genome each (Table 1). The MGE tool identified that only the ESBL gene blaCTX-M-27 was located in the same contig as a plasmid, including three associations with plasmid replicon IncFII(K), and one association with the plasmid replicon IncFIA(Hl1), together with the floR gene. The genes qnrB5 and qnrB81 were associated with plasmid replicon Col440I in one isolate, while blaOXA-1 was associated with IncFII(K) in one isolate.

 Sequence types

The CGE MLST tool identified ten different K. pneumoniae sequence types (STs), including ST1, ST11, ST37, ST45, ST111, ST307, ST348, ST789, ST881, and ST983 (Table 1). Four STs (ST348 (n = 2), ST983 (n = 2), ST789 (n = 2), and ST11 (n = 3)) were found in more than one source (Table 1 and Figure 1). No STs were shared between dairy and swine farms. The cgMLST phylogeny showed that isolates with the same ST were closely related, whether they were recovered from the same farm (ST348, ST983, ST11) or from different farms (ST789, ST11) (Figure 1). However, several nodes on the cgMLST phylogenetic tree had a low bootstrap value, limiting robust inferences of relationships across different STs.

SNP analyses and the comparison of genetic profiles (including AMR genes, virulence genes, and plasmid replicons) showed a high degree of genomic similarity among isolates belonging to the same Kp ST (Figure 2). On farm Huaura-G8, the two ST983 isolates (one recovered from a human and one from a cow) differed by 37 SNPs and shared identical virulence gene and plasmid replicon profiles, differing only by a single AMR gene (oqxB). Similarly, on farm Huaura-G3, the two ST348 isolates from a cow and the water source differed by 26 SNPs and had nearly identical genetic profiles, again differing only by the oqxB gene. For ST789, isolates recovered from a pig (Farm San Bartolo-G3) and a dog (Farm San Bartolo-G6) differed by 28 SNPs. These isolates carried the same plasmid replicons but differed by two virulence genes (irp1, iutA) and seven AMR genes, including five qnrB alleles, adeF, and oqxB. All three ST11 isolates differed by 24 SNPs or less, and all harbored the AmpC β-lactamase gene blaDHA-1. The ST11 isolated from a dog on farm Huaura-G5 differed by 21 SNPs from the ST11 isolate recovered from a cow on the same farm, and by 21 SNPs from a ST11 isolate obtained from one cow of another farm (Huaura G4). The two ST11 isolates from cows of different farms differed by 24 SNPs. Across the ST11 isolates, gene content variation was limited, with differences of no more than two AMR genes and two virulence genes. The two ST11 isolates from the same farm (Huaura G5) carried the same five plasmid replicons (Col(pHAD28), IncFIB(K), IncFIB(pKPHS1), IncFII(K), repB(R1701)). In contrast, the remaining ST11 isolate from the other farm carried Col440I instead of Col(pHAD28) (Table 1).

The SNP phylogeny of the four ESBL-Kp ST307 isolates revealed pairwise SNP differences ranging from 50 to 300. In particular, the ST307 isolate identified in this study differed by 251 SNPs from the ST307 isolate previously recovered in 2017 from a pig in the same region.

 Discussion

The dynamics of ESBL-Klebsiella pneumoniae (ESBL-Kp) in the community remain poorly understood, despite increasing evidence of community-acquired infections caused by these bacteria. In particular, the extent of ESBL-Kp circulation across the human–animal–environment interface in LMICs, including those in Latin America, remains unknown. In this study, we detected ESBL-Kp across all sampled source types in a rural setting in coastal Peru, including humans, pigs, cows, dogs, and water. This finding suggests that ESBL-Kp circulates broadly across sources, despite a low overall isolation frequency (2.3%). Genomic analyses revealed a diverse population structure, comprising 10 different ESBL-Kp sequence types, with the majority of isolates (10/15) carrying the ESBL gene blaCTX-M-15. The highly genomic similarity and closely related genetic profiles observed among the four STs detected across different host species and environmental sources suggest clonal exchanges of ESBL-Kp at this farmer–animal–water interface.

ESBL-Kp is of critical importance to public health because of the difficulties in treating infections caused by this bacterium (74). The proportion of hospital- and community-acquired infections caused by ESBL-Kp is relatively low (<5%) among Kp infections (75, 76). Therefore, a low fecal carriage rate of ESBL-Kp in human and animal populations is also expected (19, 26). Our findings support this expectation: only a small proportion of human and animal fecal samples (<3%) carried ESBL-Kp. Despite this low detection rate, ESBL-Kp was found across all sampled source types (humans, livestock, dogs, and water) and across multiple farms and localities. This widespread distribution highlights the need for a better understanding of the selective pressures and circulation pathways that maintain these bacteria in this rural environment. The absence of apparent clinical signs in the sampled animals and humans, together with the absence of the five genes (iucA, iroB, peg-344, rmpA, and rmpA2) associated with hypervirulence in the detected bacteria (77), suggests that these strains are classical, likely commensal K. pneumoniae.

The molecular mechanisms driving the spread of ESBL-Kp were first attributed to SHV- and TEM-type β-lactamases genes in the Americas, followed by the expansion of CTX-M-encoding genes, particularly blaCTX-M-15 in North America (28, 78). Our genomic analyses revealed a high diversity of ESBL-Kp STs including high-risk global clones ST307 and ST11 (32, 33), and three CTX-M-encoding alleles, confirming the predominant circulation of blaCTX-M-15. A high clonal diversity of ESBL-Kp has also been described in hospital settings (79) and in community-acquired infections (31, 80). In our study, only the gene blaCTX-M-27 was detected on the same contig as a plasmid, suggesting that blaCTX-M-15 and blaCTX-M-14 may be chromosomally located. However, long-read sequencing would be needed to establish the exact location of these genes and to confirm this hypothesis, particularly since blaCTX-M-15 is often carried in plasmids (81, 82). The ESBL-encoding genes blaCTX-M-15 and blaCTX-M-14, and the narrow-spectrum β-lactamases blaSHV-1, blaSHV-11, and blaSHV-27, have been previously reported in Kp from hospitals in Peru (37, 39, 40, 83). This indicates that multiple molecular mechanisms responsible for ESBL-Kp are selected in both hospital and community settings, including rural areas. We also found the high-risk clone ESBL-Kp ST307 in a pig carrying the same ESBL-encoding genes as an ST307 previously isolated from a pig in the same area in 2017 (41). Although the genomic differences between these isolates (i.e., 251 SNPs) do not support them as being the same clone, these findings suggest that ESBL-Kp has been circulating within this animal population for several years rather than representing a recent introduction.

Circulation of ESBL-Enterobacterales across the components of the human–animal–environment interface can have important consequences for both human and animal health. However, such transmission is often difficult to detect, since most studies do not sample multiple sources at the same time and location (84). In this study, we identified four ESBL-Kp STs present in more than one source. Phylogenomic analyses of these STs were consistent with potential cross-species clonal circulation events. In particular, we detected a farmer and a cow from the same farm harboring almost identical ESBL-Kp ST983 genomes and genetic profiles Kp, suggesting either cross-species transmission or contamination from a common source. This result is in line with previous evidence of human–dog transmission in household settings in high-income countries (85), but contrasts with findings from Guadeloupe island, where limited cross-species transmission was observed (86). We also detected very similar ESBL-Kp ST11 isolates collected from a dog and a cow from the same farm, as well as from a cow on a different farm. This suggests potential circulation of the same ESBL-Kp clone between farms, possibly mediated by the movement of animals, people, or contaminated materials. Together, these findings highlight the need to identify practices and sanitary conditions that facilitate bacterial transmission between humans and their livestock in rural settings.

The environment of humans and animals, including water sources, can also participate in the circulation of ESBL-Kp (87, 88). However, the role of water in the spread of ESBL-Kp in rural environments in Latin America remains poorly studied, particularly in Peru (87). We detected ESBL-Kp in the water sources of four farms (one swine farm and three dairy farms). These isolates belonged to four different STs (ST1, ST348, ST881, and ST37). ESBL-Kp ST348 was nearly identical to an isolate from a cow collected on the same farm. The results of this study support the need to design One Health-oriented interventions to interrupt human–animal–environment circulation, including measures to reduce water contamination by ESBL-Kp.

Our study represents a first step in understanding ESBL-Kp circulation in Peruvian communities. However, several limitations should be addressed by future research. First, the overall ESBL-Kp genetic diversity circulating at this interface is likely underestimated, particularly since we relied on a single-colony selection per host. Although a low within-host diversity of ESBL-Kp has been estimated in humans (89), studying within-host diversity could allow a better understanding of the full extend of ESBL-Kp diversity in this interface, particularly for animal or water samples. Second, although this study included a large number of sample, its cross-sectional design does not allow the evaluation of temporal trends, which could contribute to inferring directionality of bacterial circulation across hosts and environments. Longitudinal studies however, are difficult to conduct and require sustained funding and long-term farmer engagement. Finally, future studies could further characterize the molecular mechanisms involved in the spread of the detected ESBL-encoding genes. For example, the use of long-read whole-genome sequencing would allow determining whether some of the detected genes (e.g., blaCTX-M-15) are located on plasmids.

 Data availability statement

Sequences are available at the European Bioinformatics Institute (EMBL-EBI) using the project accession number PRJEB107091.

 Ethics statement

The studies involving humans were approved by the Comité Institucional de Ética en Humanos of the Universidad Peruana Cayetano Heredia. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. The animal studies were approved by the Comité Institucional de Ética para el Uso de Animales of the Universidad Peruana Cayetano Heredia. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.

 Author contributions

JB: Methodology, Conceptualization, Validation, Writing – review & editing, Supervision, Resources, Formal analysis, Visualization, Writing – original draft, Funding acquisition, Project administration. MS-C: Investigation, Visualization, Data curation, Formal analysis, Writing – review & editing, Methodology, Writing – original draft. PE-P: Data curation, Formal analysis, Visualization, Writing – review & editing, Writing – original draft. DL: Data curation, Writing – review & editing, Investigation, Writing – original draft. LJ: Investigation, Data curation, Formal analysis, Writing – review & editing. OB-M: Writing – review & editing, Investigation. CM: Investigation, Writing – review & editing. PM: Investigation, Writing – review & editing. VC: Writing – review & editing, Investigation. BA-J: Writing – review & editing, Investigation. RS-T: Investigation, Writing – review & editing. JN: Investigation, Writing – review & editing. RO-L: Writing – review & editing, Investigation. ES-T: Writing – review & editing, Investigation. CS: Formal analysis, Writing – original draft, Methodology, Supervision, Validation, Funding acquisition, Conceptualization, Writing – review & editing, Resources, Project administration.

 Acknowledgments

We thank all livestock associations and farmers involved in this study for their cooperation, and the students and laboratory technicians from the Universidad Peruana Cayetano Heredia (Lima, Peru) and the Hospital Arnaud de Villeneuve (Montpellier, France) for their assistance with laboratory analyses.

 Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Generative AI statement

The author(s) declared that Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2026.1736106/full#supplementary-material

 References PeñaC PujolM ArdanuyC RicartA PallarésR LiñaresJ . An outbreak of hospital-acquired Klebsiella pneumoniae bacteraemia, including strains producing extended-spectrum β-lactamase. J Hosp Infect. (2001) 47:53–9. doi: 10.1053/jhin.2000.0862, 11161899 ChongY ShimodaS ShimonoN. Current epidemiology, genetic evolution and clinical impact of extended-spectrum β-lactamase-producing Escherichia coli and Klebsiella pneumoniae. Infect Genet Evol. (2018) 61:185–8. doi: 10.1016/j.meegid.2018.04.005 RamatlaT MafokwaneT LekotaK MonyamaM KhasapaneG SerageN . “One health” perspective on prevalence of co-existing extended-spectrum β-lactamase (ESBL)-producing Escherichia coli and Klebsiella pneumoniae: a comprehensive systematic review and meta-analysis. Ann Clin Microbiol Antimicrob. (2023) 22:88. doi: 10.1186/s12941-023-00638-3, 37740207 Mohd AsriNA AhmadS MohamudR Mohd HanafiN Mohd ZaidiNF IrekeolaAA . Global prevalence of nosocomial multidrug-resistant Klebsiella pneumoniae: a systematic review and Meta-analysis. Antibiotics. (2021) 10:1508. doi: 10.3390/antibiotics10121508 ChangD SharmaL Dela CruzCS ZhangD. Clinical epidemiology, risk factors, and control strategies of Klebsiella pneumoniae infection. Front Microbiol. (2021) 12:750662. doi: 10.3389/fmicb.2021.750662 RozenkiewiczD Esteve-PalauE Arenas-MirasM GrauS DuranX SorlíL . Clinical and economic impact of community-onset urinary tract infections caused by ESBL-producing Klebsiella pneumoniae requiring hospitalization in Spain: an observational cohort study. Antibiotics. (2021) 10:585. doi: 10.3390/antibiotics10050585 KarlowskyJA LobSH DeRykeCA SiddiquiF YoungK MotylMR . Prevalence of ESBL non-CRE Escherichia coli and Klebsiella pneumoniae among clinical isolates collected by the SMART global surveillance programme from 2015 to 2019. Int J Antimicrob Agents. (2022) 59:106535. doi: 10.1016/j.ijantimicag.2022.106535, 35091052 GalesAC CastanheiraM JonesRN SaderHS. Antimicrobial resistance among gram-negative bacilli isolated from Latin America: results from SENTRY antimicrobial surveillance program (Latin America, 2008–2010). Diagn Microbiol Infect Dis. (2012) 73:354–60. doi: 10.1016/j.diagmicrobio.2012.04.007, 22656912 SarmientoAF MejiaMN RozasAB GiraldoCA MálagaG. Frecuencia y factores de riesgo para bacteriemia por enterobacterias productoras de betalactamasa de espectro extendido en pacientes de un hospital público de Lima, Perú. Rev Peru Med Exp Salud Publica. (2018) 35:62–7. doi: 10.17843/rpmesp.2018.351.3601 Tejada-LlacsaPJ HuarcayaJM MelgarejoGC GonzalesLF CahuanaJ PariRM . Caracterización de infecciones por bacterias productoras de BLEE en un hospital de referencia nacional. Anales Facult Medicina. (2015) 76:161–6. doi: 10.15381/anales.v76i2.11143 Cardenas-AlvarezJ BalaylaG TrianaA Diaz LankenauR Franco-ParedesC Henao-MartínezAF . Clinical Spectrum and outcomes of cryptogenic Klebsiella pneumoniae liver abscess in the Americas: a scoping review. Pathogens. (2023) 12:661. doi: 10.3390/pathogens12050661 RichelsenR SmitJ AnruPL SchønheyderHC NielsenH. Incidence of community-onset extended-spectrum β-lactamase-producing Escherichia coli and Klebsiella pneumoniae infections: an 11-year population-based study in Denmark. Infect Dis. (2020) 52:547–56. doi: 10.1080/23744235.2020.1763452 QiuY LinD XuY ChengY WangF ZhuQ . Invasive Klebsiella pneumoniae infections in community-settings and healthcare settings. Infect Drug Resist. (2021) 14:2647–56. doi: 10.2147/IDR.S315871 34285517 QuanJ ZhaoD LiuL ChenY ZhouJ JiangY . High prevalence of ESBL-producing Escherichia coli and Klebsiella pneumoniae in community-onset bloodstream infections in China. J Antimicrob Chemother. (2017) 72:273–80. doi: 10.1093/jac/dkw372, 27624571 LohiyaA KantS KapilA GuptaSK MisraP RaiSK. Pattern of antibiotic resistance among community derived isolates of Enterobacteriaceae using urine sample: a study from northern India. J Clin Diagn Res. (2015) 9:LC15–9. doi: 10.7860/JCDR/2015/14230.6254 PereraPDVM GamageS De SilvaHSM JayatillekeSK de SilvaN AydinA . Phenotypic and genotypic distribution of ESBL, AmpC β-lactamase and carbapenemase-producing Enterobacteriaceae in community-acquired and hospital-acquired urinary tract infections in Sri Lanka. J Global Antimicrob Resist. (2022) 30:115–22. doi: 10.1016/j.jgar.2022.05.024 Cabrera-RodríguezLE Miralles-SuarezAI Ones-RoqueR Torres-HerreraY Pantaleón-HernándezM. Prevalence of extended-spectrum β-lactamase-producing Escherichia coli and Klebsiella pneumoniae isolates in outpatients with urinary tract infection in a Cuban municipality. Rev Fac Med. (2023) 71:e99767. doi: 10.15446/revfacmed.v71n2.99767 MonteiroA d SS SilvaM d O GalvãoVS BomfimAP de AraújoLG SilvaCMP . High proportions of multidrug-resistant Klebsiella pneumoniae isolates in community-acquired infections, Brazil. Sci Rep. (2025) 15:9698. doi: 10.1038/s41598-025-92549-w RaffelsbergerN BuczekDJ SvendsenK SmåbrekkeL PöntinenAK LöhrIH . Community carriage of ESBL-producing Escherichia coli and Klebsiella pneumoniae: a cross-sectional study of risk factors and comparative genomics of carriage and clinical isolates. mSphere. (2023) 8:e00025–3. doi: 10.1128/msphere.00025-23 BabuR KumarA KarimS WarrierS NairSG SinghSK . Faecal carriage rate of extended-spectrum β-lactamase-producing Enterobacteriaceae in hospitalised patients and healthy asymptomatic individuals coming for health check-up. J Global Antimicrob Resist. (2016) 6:150–3. doi: 10.1016/j.jgar.2016.05.007, 27530858 AmareA EshetieS KasewD MogesF. High prevalence of fecal carriage of extended-spectrum beta-lactamase and carbapenemase-producing Enterobacteriaceae among food handlers at the University of Gondar, Northwest Ethiopia. PLoS One. (2022) 17:e0264818. doi: 10.1371/journal.pone.0264818, 35298493 SapkotaB YadavSK DhunganaG AnsariS MishraSK. Intestinal carriage of extended-Spectrum β-lactamase- (ESBL-) possessing Escherichia coli and Klebsiella species among Nepalese health science and non-health science students. Can J Infect Dis Med Microbiol. (2021) 2021:4767429. doi: 10.1155/2021/4767429, 33897921 Ouchar MahamatO TidjaniA LounnasM HideM BenavidesJ SomasseC . Fecal carriage of extended-spectrum β-lactamase-producing Enterobacteriaceae in hospital and community settings in Chad. Antimicrob Resist Infect Control. (2019) 8:169. doi: 10.1186/s13756-019-0626-z, 31695911 WuX LiuJ FengJ ShabbirMAB FengY GuoR . Epidemiology, environmental risks, virulence, and resistance determinants of Klebsiella pneumoniae from dairy cows in Hubei, China. Front Microbiol. (2022) 13:858799. doi: 10.3389/fmicb.2022.858799 KaspersenH UrdahlAM Franklin-AlmingFV IlagHK HetlandMAK BernhoffE . Population dynamics and characteristics of Klebsiella pneumoniae from healthy poultry in Norway. Front Microbiol. (2023) 14:1193274. doi: 10.3389/fmicb.2023.1193274 LeangapichartT LunhaK JiwakanonJ AngkititrakulS JärhultJD MagnussonU . Characterization of Klebsiella pneumoniae complex isolates from pigs and humans in farms in Thailand: population genomic structure, antibiotic resistance and virulence genes. J Antimicrob Chemother. (2021) 76:2012–6. doi: 10.1093/jac/dkab118, 33829268 WarethG NeubauerH. The animal-foods-environment interface of Klebsiella pneumoniae in Germany: an observational study on pathogenicity, resistance development and the current situation. Vet Res. (2021) 52:16. doi: 10.1186/s13567-020-00875-w, 33557913 Guzmán-BlancoM LabarcaJA VillegasMV GotuzzoE. Extended spectrum β-lactamase producers among nosocomial Enterobacteriaceae in Latin America. Braz J Infect Dis. (2014) 18:421–33. doi: 10.1016/j.bjid.2013.10.005 HartingerSM Medina-PizzaliML Salmon-MulanovichG LarsonAJ Pinedo-BardalesM VerasteguiH . Antimicrobial resistance in humans, animals, water and household environs in rural Andean Peru: exploring dissemination pathways through the one health lens. Int J Environ Res Public Health. (2021) 18:4604. doi: 10.3390/ijerph18094604 BenavidesJA StreickerDG GonzalesMS Rojas-PaniaguaE ShivaC. Knowledge and use of antibiotics among low-income small-scale farmers of Peru. Prev Vet Med. (2021) 189:105287. doi: 10.1016/j.prevetmed.2021.105287 ZhangJ ZhouK ZhengB ZhaoL ShenP JiJ . High prevalence of ESBL-producing Klebsiella pneumoniae causing community-onset infections in China. Front Microbiol. (2016) 7:01830. doi: 10.3389/fmicb.2016.01830 Navon-VeneziaS KondratyevaK CarattoliA. Klebsiella pneumoniae: a major worldwide source and shuttle for antibiotic resistance. FEMS Microbiol Rev. (2017) 41:252–75. doi: 10.1093/femsre/fux013, 28521338 WyresKL LamMMC HoltKE. Population genomics of Klebsiella pneumoniae. Nat Rev Microbiol. (2020) 18:344–59. doi: 10.1038/s41579-019-0315-1, 32055025 MedinaVA García-VegaAS RodríguezF Triviño-GarcíaES RobledoJA LealAL . Epidemiological genomics of Klebsiella pneumoniae isolates from hospitals across Colombia. NPJ Antimicrob Resist. (2025) 3:64. doi: 10.1038/s44259-025-00127-x SennatiS SantellaG Di ConzaJ PallecchiL PinoM GhiglioneB . Changing epidemiology of extended-spectrum β-lactamases in Argentina: emergence of CTX-M-15. Antimicrob Agents Chemother. (2012) 56:6003–5. doi: 10.1128/aac.00745-12 22908156 SartoriL SelleraFP MouraQ CardosoB CerdeiraL LincopanN. Multidrug-resistant CTX-M-15-positive Klebsiella pneumoniae ST307 causing urinary tract infection in a dog in Brazil. J Glob Antimicrob Resist. (2019) 19:96–7. doi: 10.1016/j.jgar.2019.09.003, 31520809 Aguilar-AncoriEG Marin-CarrascoM Campo-PfuyoLI Muñiz-DuranJG Espinoza-CulupúA. Identification of pandemic ST147, ESBL-type β-lactamases, carbapenemases, and virulence factors in Klebsiella pneumoniae isolated from southern Peru. Sci Rep. (2025) 15:14870. doi: 10.1038/s41598-025-97464-8 Palomino-KobayashiLA Egoávil-EspejoR Salvador-LujánG YáñezP GavilánRG PonsMJ . Detection of rare β-lactamase blaSCO-1 from a Klebsiella pneumoniae high-risk clone in Peru. JAC-Antimicrob Resist. (2025) 7:dlaf089. doi: 10.1093/jacamr/dlaf089, 40421127 HerbozoC JulcaI FloresF HernandezR ZegarraJ. Incidence and microbiological characteristics of neonatal late onset sepsis in a neonatal intensive care unit in Peru. Int J Infect Dis. (2021) 108:171–5. doi: 10.1016/j.ijid.2021.05.012, 34004330 GarcíaC AstocondorL Rojo-BezaresB JacobsJ SáenzY. Molecular characterization of extended-Spectrum β-lactamase-producer Klebsiella pneumoniae isolates causing neonatal Sepsis in Peru. Am J Trop Med Hyg. (2016) 94:285–8. doi: 10.4269/ajtmh.15-0373, 26643537 Fuentes-CastilloD ShivaC LincopanN SanoE FontanaH StreickerDG . Global high-risk clone of extended-spectrum β-lactamase (ESBL)-producing Klebsiella pneumoniae ST307 emerging in livestock in Peru. Int J Antimicrob Agents. (2021) 58:106389. doi: 10.1016/j.ijantimicag.2021.106389 Salgado-CaxitoM LéonD BardalesO JaraLM MedranoP MurgaC . Unexplained high prevalence of ESBL-Escherichia coli among cattle and pigs in Peru. Antibiotics. (2025) 14:867. doi: 10.3390/antibiotics14090867 BenavidesJA GodreuilS Opazo-CapurroA MahamatOO FalconN OravcovaK . Long-term maintenance of multidrug-resistant Escherichia coli carried by vampire bats and shared with livestock in Peru. Sci Total Environ. (2022) 810:152045. doi: 10.1016/j.scitotenv.2021.152045, 34883172 MonteroL IrazabalJ CardenasP GrahamJP TruebaG. Extended-spectrum beta-lactamase producing-Escherichia coli isolated from irrigation waters and produce in Ecuador. Front Microbiol. (2021) 12:709418. doi: 10.3389/fmicb.2021.709418/full HosainMZ KabirSML KamalMM. Antimicrobial uses for livestock production in developing countries. Vet World. (2021) 14:210–21. doi: 10.14202/vetworld.2021.210-221 ChowdhuryS HassanMM AlamM SattarS BariMS SaifuddinAKM . Antibiotic residues in milk and eggs of commercial and local farms at Chittagong, Bangladesh. Vet World. (2015) 8:467–71. doi: 10.14202/vetworld.2015.467-471 LewisJS, editor. Performance standards for antimicrobial susceptibility testing. 34th ed. Malvern, PA: Clinical and Laboratory Standards Institute. (2024)., 24642063 BolgerAM LohseM UsadelB. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. (2014) 30:2114–20. doi: 10.1093/bioinformatics/btu170, 24695404 PrjibelskiA AntipovD MeleshkoD LapidusA KorobeynikovA. Using SPAdes De Novo Assembler. Curr Protoc Bioinformatics. (2020) 70:e102. doi: 10.1002/cpbi.102, 32559359 SeemannT. Prokka: rapid prokaryotic genome annotation. Bioinformatics. (2014) 30:2068–9. doi: 10.1093/bioinformatics/btu153, 24642063 GurevichA SavelievV VyahhiN TeslerG. QUAST: quality assessment tool for genome assemblies. Bioinformatics. (2013) 29:1072–5. doi: 10.1093/bioinformatics/btt086, 23422339 ManniM BerkeleyMR SeppeyM SimãoFA ZdobnovEM. BUSCO update: novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. Mol Biol Evol. (2021) 38:4647–54. doi: 10.1093/molbev/msab199, 34320186 WoodDE LuJ LangmeadB. Improved metagenomic analysis with kraken 2. Genome Biol. (2019) 20:257. doi: 10.1186/s13059-019-1891-0, 31779668 OlsonRD AssafR BrettinT ConradN CucinellC DavisJJ . Introducing the bacterial and viral bioinformatics resource center (BV-BRC): a resource combining PATRIC, IRD and ViPR. Nucleic Acids Res. (2023) 51:D678–89. doi: 10.1093/nar/gkac1003 WattamAR DavisJJ AssafR BoisvertS BrettinT BunC . Improvements to PATRIC, the all-bacterial bioinformatics database and analysis resource center. Nucleic Acids Res. (2017) 45:D535–42. doi: 10.1093/nar/gkw1017, 27899627 GriffithsD FawleyW KachrimanidouM BowdenR CrookDW FungR . Multilocus sequence typing of Clostridium difficile. J Clin Microbiol. (2010) 48:770–8. doi: 10.1128/jcm.01796-09, 20042623 CarattoliA ZankariE García-FernándezA Voldby LarsenM LundO VillaL . In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. (2014) 58:3895–903. doi: 10.1128/AAC.02412-14, 24777092 CamachoC CoulourisG AvagyanV MaN PapadopoulosJ BealerK . BLAST+: architecture and applications. BMC Bioinform. (2009) 10:421. doi: 10.1186/1471-2105-10-421, 20003500 LarsenMV CosentinoS RasmussenS FriisC HasmanH MarvigRL . Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol. (2012) 50:1355–61. doi: 10.1128/JCM.06094-11, 22238442 BartualSG SeifertH HipplerC LuzonMAD WisplinghoffH Rodríguez-ValeraF. Development of a multilocus sequence typing scheme for characterization of clinical isolates of Acinetobacter baumannii. J Clin Microbiol. (2005) 43:4382–90. doi: 10.1128/jcm.43.9.4382-4390.2005 LemeeL DhalluinA Pestel-CaronM LemelandJF PonsJL. Multilocus sequence typing analysis of human and animal Clostridium difficile isolates of various toxigenic types. J Clin Microbiol. (2004) 42:2609–17. doi: 10.1128/JCM.42.6.2609-2617.2004, 15184441 WirthT FalushD LanR CollesF MensaP WielerLH . Sex and virulence in Escherichia coli: an evolutionary perspective. Mol Microbiol. (2006) 60:1136–51. doi: 10.1111/j.1365-2958.2006.05172.x, 16689791 JaureguyF LandraudL PassetV DiancourtL FrapyE GuigonG . Phylogenetic and genomic diversity of human bacteremic Escherichia coli strains. BMC Genomics. (2008) 9:560. doi: 10.1186/1471-2164-9-560, 19036134 JohanssonMHK BortolaiaV TansirichaiyaS AarestrupFM RobertsAP PetersenTN. Detection of mobile genetic elements associated with antibiotic resistance in Salmonella enterica using a newly developed web tool: MobileElementFinder. J Antimicrob Chemother. (2021) 76:101–9. doi: 10.1093/jac/dkaa390, 33009809 Suarez-YanaT Salgado-CaxitoM HayerJ Rojas-SerenoZE Pino-HurtadoMS Campaña-BurguetA . ESBL-producing Escherichia coli prevalence and sharing across seabirds of Central Chile. Sci Total Environ. (2024) 951:175475. doi: 10.1016/j.scitotenv.2024.175475, 39142400 HayerJ Salgado-CaxitoM Opazo-CapurroA MuñozPG MillánJ PiñeiroA . Multiple clonal transmissions of clinically relevant extended-spectrum beta-lactamase-producing Escherichia coli among livestock, dogs, and wildlife in Chile. J Glob Antimicrob Resist. (2023) 34:247–52. doi: 10.1016/j.jgar.2023.07.009, 37463613 PatiNB DoijadSP SchultzeT MannalaGK YaoY JaiswalS . Enterobacter bugandensis: a novel enterobacterial species associated with severe clinical infection. Sci Rep. (2018) 8:5392. doi: 10.1038/s41598-018-23069-z LemoineF CorreiaD LefortV Doppelt-AzeroualO MareuilF Cohen-BoulakiaS . NGPhylogeny.Fr: new generation phylogenetic services for non-specialists. Nucleic Acids Res. (2019) 47:W260–5. doi: 10.1093/nar/gkz303 LiH. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. (2014). Available online at: https://figshare.com/articles/poster/Aligning_sequence_reads_clone_sequences_and_assembly_con_gs_with_BWA_MEM/963153/1 (Accessed October 22, 2025) LiH HandsakerB WysokerA FennellT RuanJ HomerN . The sequence alignment/map format and SAMtools. Bioinformatics. (2009) 25:2078–9. doi: 10.1093/bioinformatics/btp352, 19505943 QuinlanAR HallIM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. (2010) 26:841–2. doi: 10.1093/bioinformatics/btq033, 20110278 MarçaisG DelcherAL PhillippyAM CostonR SalzbergSL ZiminA. MUMmer4: A fast and versatile genome alignment system. PLoS Comput Biol. (2018) 14:e1005944. doi: 10.1371/journal.pcbi.1005944, 29373581 PriceMN DehalPS ArkinAP. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS One. (2010) 5:e9490. doi: 10.1371/journal.pone.0009490, 20224823 World Health Organization. WHO Publishes list of bacteria for which new antibiotics are urgently needed. (2017). Available online at: https://www.who.int/news/item/27-02-2017-who-publishes-list-of-bacteria-for-which-new-antibiotics-are-urgently-needed (Accessed October 22, 2025). KoWC PatersonDL SagnimeniAJ HansenDS Von GottbergA MohapatraS . Community-acquired Klebsiella pneumoniae bacteremia: global differences in clinical patterns. Emerg Infect Dis. (2002) 8:160–6. doi: 10.3201/eid0802.010025 ToubianaJ TimsitS FerroniA GrasseauM NassifX LortholaryO . Community-onset extended-Spectrum β-lactamase-producing Enterobacteriaceae invasive infections in children in a University Hospital in France. Medicine (Baltimore). (2016) 95:e3163. doi: 10.1097/MD.0000000000003163, 27015202 TangY DuP DuC YangP ShenN RussoTA . Genomically defined hypervirulent Klebsiella pneumoniae contributed to early-onset increased mortality. Nat Commun. (2025) 16:2096. doi: 10.1038/s41467-025-57379-4, 40025046 WangG HuangT SurendraiahPKM WangK KomalR ZhugeJ . CTX-M β-lactamase–producing Klebsiella pneumoniae in suburban new York City, New York, USA. Emerg Infect Dis. (2013) 19:1803–10. doi: 10.3201/eid1911.121470, 24188126 XercavinsM JiménezE PadillaE RieraM FreixasN Boix-PalopL . High clonal diversity of ESBL-producing Klebsiella pneumoniae isolates from clinical samples in a non-outbreak situation. A cohort study. Antimicrob Resist Infect Control. (2020) 9:5. doi: 10.1186/s13756-019-0661-9 LepuschitzS HauserK SchrieblA SchlagenhaufenC StögerA ChakeriA . Fecal Klebsiella pneumoniae carriage is intermittent and of high clonal diversity. Front Microbiol. (2020) 11:581081. doi: 10.3389/fmicb.2020.581081, 33324367 ZhuoC LiX qiang, ZongZ yong, ZhongNS. Epidemic plasmid carrying blaCTX-M-15 in Klebsiella penumoniae in China. PLoS One 2013;8:e52222 CalboE GarauJ. The changing epidemiology of hospital outbreaks due to ESBL-producing Klebsiella Pneumoniae: the CTX-M-15 type consolidation. Future Microbiol. (2015) 10:1063–75. doi: 10.2217/fmb.15.22, 26059626 KrappF CuicapuzaD SalvatierraG ButeauJP AmaroC AstocondorL . Emerging carbapenem-resistant Klebsiella pneumoniae in a tertiary care hospital in Lima, Peru. Microbiol Spectr. (2025). doi: 10.1128/spectrum.01825-24 MilenkovM ProuxC RasolofoarisonTL RakotomalalaFA RasoanandrasanaS RahajamananaVL . Implementation of the WHO tricycle protocol for surveillance of extended-spectrum β-lactamase producing Escherichia coli in humans, chickens, and the environment in Madagascar: a prospective genomic epidemiology study. Lancet Microbe. (2024) 5:100850. doi: 10.1016/S2666-5247(24)00065-X, 38908389 MarquesC BelasA AboimC Cavaco-SilvaP TrigueiroG GamaLT . Evidence of sharing of Klebsiella pneumoniae strains between healthy companion animals and cohabiting humans. J Clin Microbiol. (2019) 57:e01537–18. doi: 10.1128/JCM.01537-18 DereeperA GruelG PotM CouvinD BarbierE BastianS . Limited transmission of Klebsiella pneumoniae among humans, animals, and the environment in a Caribbean Island, Guadeloupe (French West Indies). Microbiol Spectr. (2022) 10:e01242–22. ChoS JacksonCR FryeJG. Freshwater environment as a reservoir of extended-spectrum β-lactamase-producing Enterobacteriaceae. J Appl Microbiol. (2023) 134:lxad034. doi: 10.1093/jambio/lxad034, 36806844 Davidova-GerzovaL LausovaJ SukkarI NechutnaL KubackovaP KrutovaM . Multidrug-resistant ESBL-producing Klebsiella pneumoniae complex in Czech hospitals, wastewaters and surface waters. Antimicrob Resist Infect Control. (2024) 13:141. doi: 10.1186/s13756-024-01496-0, 39593189 Within-host genetic diversity of extended-spectrum beta-lactamase-producing Enterobacterales in long-term colonized patients | nature communications. Nat Commun. (2023) 14:8495. doi: 10.1038/s41467-023-44285-w

Edited by: Corina-Diana Ceapă, National Autonomous University of Mexico, Mexico

Reviewed by: Tobias Manuel Appel, University Hospital Frankfurt, Germany

Jhonathan Bazalar-Gonzales, National University of San Marcos, Peru

https://www.bv-brc.org

http://www.genomicepidemiology.org/

http://NGPhylogeny.fr

https://bigsdb.pasteur.fr/, data submitted up to 2024-12-31

https://tree.bio.ed.ac.uk/software/figtree/