# Limited Transmission of Klebsiella pneumoniae among Humans, Animals, and the Environment in a Caribbean Island, Guadeloupe (French West Indies)

Alexis Dereeper, Gaëlle Gruel, Matthieu Pot, David Couvin, Elodie Barbier, Sylvaine Bastian, Jean-Christophe Bambou, Moana Gelu-Simeon, Séverine Ferdinand, Stéphanie Guyomard-Rabenirina, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Alexis Dereeper, Gaëlle Gruel, Matthieu Pot, David Couvin, Elodie Barbier, et al. Limited Transmission of Klebsiella pneumoniae among Humans, Animals, and the Environment in a Caribbean Island, Guadeloupe (French West Indies). Microbiology Spectrum, 2022, 10 (5), 10.1128/spectrum.01242-22. hal-03831169

HAL Id: hal-03831169

https://hal.inrae.fr/hal-03831169

Submitted on 3 Nov 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

# Limited Transmission of *Klebsiella pneumoniae* among Humans, Animals, and the Environment in a Caribbean Island, Guadeloupe (French West Indies)

Alexis Dereeper,ª Gaëlle Gruel,ª Matthieu Pot,ª David Couvin,ª Elodie Barbier,ᵇ Sylvaine Bastian,ᶜ Jean-Christophe Bambou,ᵈ
Moana Gelu-Simeon,ᵉ Séverine Ferdinand,ª Stéphanie Guyomard-Rabenirina,ª Virginie Passet,ᵇ Frederic Martino,ᵉ Pascal Piveteau,ʰ
Yann Reynaud,ª Carla Rodrigues,ᵇ Pierre-Marie Roger,ij Xavier Roy, k Antoine Talarmin,ª Benoit Tressieres,¹ Marc Valette, g
i>Sylvain Brisse, f>Sébastien Breurecª,ç,j,l

**ABSTRACT** Guadeloupe (French West Indies), a Caribbean island, is an ideal place to study the reservoirs of the *Klebsiella pneumoniae* species complex (KpSC) and identify the routes of transmission between human and nonhuman sources due to its insularity, small population size, and small area. Here, we report an analysis of 590 biological samples, 546 KpSC isolates, and 331 genome sequences collected between January 2018 and May 2019. The KpSC appears to be common whatever the source. Extended-specrum-β-lactamase (ESBL)-producing isolates (21.4%) belonged to *K. pneumoniae sensu stricto* (phylogroup Kp1), and all but one were recovered from the hospital setting. The distribution of species and phylogroups across the different niches was clearly nonrandom, with a distinct separation of Kp1 and *Klebsiella variicola* (Kp3). The most frequent sequence types (STs) (≥5 isolates) were previously recognized as high-risk multidrug-resistant (MDR) clones, namely, ST17, ST307, ST11, ST147, ST152, and ST45. Only 8 out of the 63 STs (12.7%) associated with human isolates were also found in nonhuman sources. A total of 22 KpSC isolates were defined as hypervirulent: 15 associated with human infections (9.8% of all human isolates), 4 (8.9%) associated with dogs, and 3 (15%) associated with pigs. Most of the human isolates (33.3%) belonged to the globally successful sublineage CG23-1. ST86 was the only clone shared by a human and a nonhuman (dog) source. Our work shows the limited transmission of KpSC isolates between human and nonhuman sources and points to the hospital setting as a cornerstone of the spread of MDR clones and antibiotic resistance genes.

**IMPORTANCE** In this study, we characterized the presence and genomic features of isolates of the *Klebsiella pneumoniae* species complex (KpSC) from human and non-human sources in Guadeloupe (French West Indies) in order to identify the reservoirs and routes of transmission. This is the first study in an island environment, an ideal setting that limits the contribution of external imports. Our data showed the limited transmission of KpSC isolates between the different compartments. In contrast, we

**Editor** Nilton Lincopan, Institute of Biomedical Sciences, Universidade de São Paulo

**Copyright** © 2022 Dereeper et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Address correspondence to Sébastien Breurec, sbreurec@gmail.com.

The authors declare no conflict of interest.

**Received** 17 April 2022

**Accepted** 8 August 2022

**Published** 12 September 2022

<footnote>ªTransmission, Reservoir and Diversity of Pathogens Unit, Pasteur Institute of Guadeloupe, Pointe-à-Pitre, France</footnote>

<footnote>ᵇUMR AgroÉcologie, INRAE, Bourgogne Franche-Comté University, Dijon, France</footnote>

<footnote>ᶜLaboratory of Clinical Microbiology, University Hospital Center of Guadeloupe, Pointe-à-Pitre/Les Abymes, France</footnote>

<footnote>ᵈINRAE, ASSET, Petit-Bourg, France</footnote>

<footnote>ᵉHepato-Gastroenterology Department, University Hospital Center of Guadeloupe, Pointe-à-Pitre/Les Abymes, France</footnote>

<footnote>fInstitut Pasteur, University Paris Cité, Biodiversity and Epidemiology of Bacterial Pathogens, Paris, France</footnote>

<footnote>gIntensive Care Department, University Hospital Center of Guadeloupe, Pointe-à-Pitre/Les Abymes, France</footnote>

<footnote>hUR OPAALE, INRAE, Rennes, France</footnote>

<footnote>iInfectious Disease Department, University Hospital Center of Guadeloupe, Pointe-à-Pitre/Les Abymes, France</footnote>

<footnote>jFaculty of Medecine Hyacinthe Bastaraud, University of the Antilles, Pointe-à-Pitre, France</footnote>

<footnote>kVeterinary Clinic, Baie-Mahault, Guadeloupe</footnote>

<footnote>INSERM Center for Clinical Investigation 1424, Pointe-à-Pitre/Les Abymes, France</footnote>
identified the hospital setting as the epicenter of antibiotic resistance due to the
nosocomial spread of successful multidrug-resistant (MDR) *K. pneumoniae* clones and
antibiotic resistance genes. Ecological barriers and/or limited exposure may restrict
spread from the hospital setting to other reservoirs and vice versa. These results
highlight the need for control strategies focused on health care centers, using
genomic surveillance to limit the spread, particularly of high-risk clones, of this im-
portant group of MDR pathogens.

KEYWORDS Klebsiella pneumoniae, One Health, ESBL, genomic, Caribbean

Klebsiella pneumoniae is currently recognized as an increasing threat to public health due to the emergence and spread of multidrug-resistant (MDR) isolates associated with hospital outbreaks. Recently, the WHO released a global priority list of antibiotic-resistant bacteria requiring new control strategies, including carbapenem-resistant *K. pneumoniae* (CRKP) and extended-spectrum-β-lactamase (ESBL)-producing *K. pneumoniae* (Kp-ESBL) isolates. ESBL and carbapenemase genes are located on mobile genetic elements and are frequently associated with genes encoding resistance to many other classes of antimicrobial agents, leading to bacteria that are difficult to treat. In Europe, it was estimated that 84,535 cases of infections with CRKP and Kp-ESBL bacteria occurred in 2015 and that these infections accounted for 5,805 attributable deaths (1). Estimations of the burden of MDR *K. pneumoniae* infections are lacking worldwide, but MDR rates are increasing globally. For example, in Senegal, *K. pneumoniae* isolates were the most common bacteria associated with neonatal bloodstream infections, of which 85% were ESBL producers (2). *K. pneumoniae* is also responsible for more severe invasive community-acquired infections, often in healthy individuals, including pyogenic liver abscess, pneumonia, and meningitis (3). Infections are caused by hypervirulent *K. pneumoniae* (HvKp) isolates, which belong to particular clonal groups (CGs). In the past, antimicrobial resistance (AMR) genes and virulence genes were present in specific nonoverlapping genomic lineages, but the frontiers are now blurring, resulting in the emergence of MDR and hypervirulent phenotypes in single *K. pneumoniae* isolates (4).

Phylogenetic studies have revealed that the former *K. pneumoniae* species is a geneti-
cally heterogeneous group. It has been redefined as the *Klebsiella pneumoniae* species
complex (KpSC), a group of five closely related species distributed into seven phylogroups
(Kp1 to Kp7): *K. pneumoniae sensu stricto* (phylogroup Kp1), *K. quasipneumoniae* subsp.
*quasipneumoniae* (Kp2), *K. quasipneumoniae* subsp. *similipneumoniae* (Kp4), *K. variicola*
subsp. *variicola* (Kp3), *K. variicola* subsp. *tropica* (Kp5), *K. quasivariicola* (Kp6), and *K. africana*
(Kp7) (5). In this report, for simplicity, KpSC refers to the *K. pneumoniae* species complex,
including the five species/seven phylogroups, and *K. pneumoniae* refers to *K. pneumoniae
sensu stricto* (phylogroup Kp1).

In addition to their importance as human pathogens, members of the KpSC can be
found in a wide range of ecological niches, such as soil, water, plants, insects, birds, reptiles,
and the gut of many mammals (6). However, the prevalence and characteristics of KpSC
isolates in these different niches are poorly known due to the lack of large dedicated
research efforts. Despite the urgent public health threat now represented by the KpSC,
knowledge of the dynamics of the transmission of this bacterial complex from environ-
mental and animal reservoirs to humans using a broad ecological approach with whole-ge-
nome sequencing (WGS) is also scarce (7–11).

Guadeloupe, a tropical French overseas territory located in the Caribbean, is consid-
ered a very high-resource territory (https://hdr.undp.org). Data on the KpSC are scarce
and recent for this island. An unusually high prevalence of HvKp (24%) was observed
in adults admitted to the intensive care units of two university hospitals in Guadeloupe
and Martinique (another French Caribbean territory) for spontaneous community-
acquired bacterial meningitis (12). Guadeloupe also faced the emergence of hospital-
acquired carbapenemase-producing *K. pneumoniae* infections (13) and an increased
incidence of nosocomial Kp-ESBL infections (14, 15). Guadeloupe is an ideal place to
study the reservoirs of KpSC isolates and identify the routes of transmission to the
human population due to its insularity, small area (1,436 km²), and small population size (395,700 inhabitants in 2019). The primary objectives of our study were to determine the genomic features of a collection of isolated nonhuman *K. pneumoniae* strains regardless of the putative antibiotic resistance phenotype and to compare them with contemporaneous clinical isolates in order to define the reservoirs of clinical isolates and whether recent transmissions of this pathogen can be detected. The secondary objectives were to determine (i) the prevalence of KpSC members from different sources and (ii) their antibiotic susceptibility patterns.

## RESULTS

**KpSC isolates, phylogenetic diversity, and AMR.** For healthy food-producing animals, a total of 199 fecal samples from 124 pigs and 75 beef cattle were collected from 28 farms and the slaughterhouse (34 additional farms). The prevalences of KpSC isolates were 52.0% in bovines (39/75) and 24.2% in pigs (30/124). KpSC isolates were recovered from all nine poultry farms investigated (11 isolates).

For pets, a single rectal swab was taken from 149 dogs and 73 cats from the main animal shelter of Guadeloupe (n = 15) and 7 veterinary clinics (n = 170). For the identification of risk factors for the fecal carriage of ESBL-producing KpSC, 37 pets were eliminated because they displayed clinical signs of diarrhea and/or received antibiotic treatment in the previous month. Most of the pets were seen for vaccination (n = 79; 42.7%), surgery (n = 33; 17.8%), a preventive health visit (n = 28; 15.1%), or skin and soft tissue infection (n = 13; 7.0%). The rates of fecal carriage of KpSC were 27.4% (20/73) among cats and 49.0% (73/149) among dogs.

Dogs (P = 0.013), compared to cats, were significantly associated with KpSC fecal carriage (see Table S1 in the supplemental material).

Considering the environment, totals of 85 locally produced fresh vegetables (n = 45), flowering plants (n = 18), fruits (n = 11 [tomatoes only]), and aromatic herbs (n = 1 [thyme]) were collected. The prevalences of KpSC isolates were 90.1% (40/44) in vegetables, 66.6% (12/18) in flowering plants, and 36.4% (4/11) in fruits (tomatoes). The only aromatic herb tested was positive (thyme). Of the 44 raw water samples tested from 29 catchment points, KpSC isolates were detected in 15 samples (34.1%). Totals of 54 soil samples and 21 water samples from rivers or natural ponds located in proximity were investigated at 21 sites located throughout Guadeloupe. The frequencies of detection were 25.9% (14/54) and 33.3% (7/21), respectively. We did not find any significant association between the presence of KpSC isolates and the level of anthropic pressure (P = 0.795).

A total of 279 contemporaneous human KpSC isolates were recovered. These bacteria were isolated mainly from urine (n = 137; 49.1%), blood (n = 65; 23.3%), and wounds (n = 62; 22.2%). Six were associated with community liver abscess (10.5% of all community isolates), and 2 were associated with meningitis (3.5%). Of the 222 (79.6%) isolates associated with nosocomial infections, most of them were collected from patients hospitalized in medical wards (n = 110; 49.6%), intensive care units (n = 82; 36.9%), and emergency units (n = 62; 27.9%).

Low levels of resistance to antibiotics were observed for KpSC isolates whatever the source, except for those associated with hospital-acquired infections, which displayed significantly higher rates of resistance rates, whatever the antibiotic tested (P < 0.027). Full details of resistance to antibiotics according to the source are available in Table S2. All isolates with decreased phenotypic susceptibility to ertapenem were observed in the hospital setting (n = 10), 3 of which displayed carbapenemase genes (2 bla<sub>KPC-2</sub> and 1 bla<sub>NDM-1</sub>). Kp-ESBL isolates (21.4%; 117/546) were associated with hospital-acquired infections only, except for 1 isolate recovered from the feces of one cat. The latter isolate also displayed resistance to fluoroquinolones. Resistance to this major family of antibiotics was not found in any isolate from other animals or the environment. All isolates collected from the environment were assigned a wild-type resistance phenotype, except for one isolate from a vegetable and two from soil (resistance to amoxicillin-
<figure><img src="image_5.png" /><figcaption>**FIG 1** Species and phylogroup distributions of 433 isolates of the *K. pneumoniae* species complex collected in Guadeloupe (French West Indies), according to source. Species were defined by real-time PCR.</figcaption></figure>

clavulanic acid). Globally, the wild-type resistance phenotype corresponded to more than 80% of the isolates detected, except for the hospital setting, where this wild-type phenotype represented 31.8% (Table S2).

Phylogroup typing was performed by real-time PCR (RT-PCR) for all isolates except those recovered from the hospital setting, where ~50% were tested. The 433 KpSC isolates investigated were assigned to *K. pneumoniae* (Kp1 [*n* = 245; 56.6%], Kp2 [*n* = 62; 14.3%], Kp3 or Kp5 [*n* = 81; 18.7%], and Kp4 [*n* = 45; 10.4%]) (Fig. 1 and Table S3). Although the distribution of species and phylogroups across the various sources was clearly nonrandom, each phylogroup was isolated from all sources. A high prevalence of Kp1 and a low prevalence of Kp3 or Kp5 were observed in isolates associated with humans, domestic animals, pigs, and poultry, whereas the opposite relative frequencies of these groups were observed for vegetables, soil/rivers/natural ponds, catchment water, and bovines (Fig. 1 and Table S3). Kp4 was observed mainly in soil (*n* = 7; 50.0%), and Kp2 was observed mainly in food-producing animals (27.3% in poultry, 33.3% in bovines, and 36.7% in pigs) and on vegetables (45.2%). No significant difference was observed when comparing isolates from human and companion animals (*P* = 0.092), those from hospital and community settings (*P* = 0.09), and those from bovines and the environment (*P* = 0.79).

**Isolate genomic diversity.** WGS was performed on a total of 331 isolates corresponding to a random selection of 66% of those collected from animals and the environment and on 55% of isolates associated with human infections. A full description of the isolates is displayed in Data Set S1 in the supplemental material. Phylogenetic analysis (Fig. 2) was performed to determine the frequencies of Kp1 (*n* = 196; 59.2%), Kp2 (*n* = 44; 13.3%), Kp3 (*n* = 50; 15.1%), and Kp4 (*n* = 41; 12.3%). Assignments were 100% concordant with the phylogrouping results for 433 isolates (see above).

High genetic diversity was found in all sources. A total of 218 sequence types (STs) were identified, with 163 being represented by a single isolate. Nineteen out of 218 STs were high-risk STs (90 isolates; 27.2% of the total), 17 of which were associated with MDR and 2 of which were associated with hypervirulence (ST23 and ST86) (Table 1). All STs represented by 5 or more isolates belonged to previously recognized successful clones. The most frequent STs were high-risk MDR clones, namely, ST17 (*n* = 16; 4.8%), ST307 (*n* = 12; 3.6%), ST11 (*n* = 11; 3.3%), ST147 (*n* = 9; 2.7%), ST152 (*n* = 8; 2.4%), and ST45 (*n* = 7; 2.1%). Most of them were human hospital-acquired Kp-ESBL isolates (47/63; 74.6%).
<figure><img src="image_6.png" /><figcaption>**FIG 2** Phylogenetic tree of 331 isolates of the *K. pneumoniae* species complex from Guadeloupe (French West Indies). A maximum likelihood tree was constructed using RAxML based on the core-genome alignment and drawn with iTOL. Leaves are colored according to the phylogroup (Kp1, Kp2, Kp3, and Kp4), and annotation tracks are displayed as follows: source (human/animal/environment), nosocomial/community origin, ESBL production, and hypervirulence.</figcaption></figure>

Only eight STs (five belonged to *K. pneumoniae*, and three belonged to Kp3) were shared between nonhuman and human isolates, with six being strictly shared with pets: three high-risk MDR clones (ST17, ST37, and ST45), one high-risk hypervirulent clone (ST86), and four minor clones (ST2551, ST3600, ST4174, and ST5685). The frequency of STs from clinical isolates also detected in nonhuman samples was 12.7% (8/63). The two isolates assigned to ST45 were common to humans and dogs. They were almost identical (<10 core-genome multilocus sequence typing [cgMLST] allelic mismatches) and were isolated at 1-month intervals (Fig. 3 and Table S4). It should be noted that an ESBL gene (*bla*<sub>CTX-M-15</sub>) was recovered only in the human isolate. The other shared isolates were genetically more distant. Within the hospital setting, we observed the spread of closely related isolates (<10 single nucleotide polymorphisms [SNPs]) belonging to the four main STs (ST11, ST17, ST45, and ST307) during the 24-month period in different units of the hospital (Fig. 3 and Table S4). For the 4 main STs, isolates from Guadeloupe and those from other Caribbean islands (16) were not directly related based on their genomic sequences (Table S4).

**AMR genes and plasmids.** In agreement with the phenotypic antibiotic susceptibility patterns, a few isolates (7.9%; 14/178) with a gene(s) or mutation(s) conferring resistance to antimicrobials were recovered outside the human sources (see Table S5 in
<table><caption>**TABLE 1** High-risk clones (sequence types) collected from different sources in Guadeloupe (French West Indies)</caption><thead><tr><th rowspan="2">Clone</th><th colspan="5">No. of clones isolated from source</th></tr><tr><th>Humans</th><th>Cats</th><th>Dogs</th><th>Pigs</th><th>Total</th></tr></thead><tbody><tr><td>Multidrug resistant</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ST11</td><td>11</td><td></td><td></td><td></td><td>11</td></tr><tr><td>ST13</td><td>4</td><td></td><td></td><td></td><td>4</td></tr><tr><td>ST14</td><td>2</td><td></td><td></td><td></td><td>2</td></tr><tr><td>ST15</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>ST17</td><td>14</td><td>2</td><td></td><td></td><td>16</td></tr><tr><td>ST20</td><td>2</td><td></td><td></td><td></td><td>2</td></tr><tr><td>ST29</td><td></td><td></td><td></td><td>1</td><td>1</td></tr><tr><td>ST36</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>ST37</td><td>1</td><td></td><td>3</td><td></td><td>4</td></tr><tr><td>ST39</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>ST45</td><td>6</td><td></td><td>1</td><td></td><td>7</td></tr><tr><td>ST101</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>ST147</td><td>9</td><td></td><td></td><td></td><td>9</td></tr><tr><td>ST152</td><td>8</td><td></td><td></td><td></td><td>8</td></tr><tr><td>ST258</td><td>2</td><td></td><td></td><td></td><td>2</td></tr><tr><td>ST307</td><td>12</td><td></td><td></td><td></td><td>12</td></tr><tr><td>ST392</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>Total</td><td>76</td><td>2</td><td>4</td><td>1</td><td>83</td></tr><tr><td>Hypervirulent</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ST23</td><td>5</td><td></td><td></td><td></td><td>5</td></tr><tr><td>ST86</td><td>1</td><td></td><td>1</td><td></td><td>2</td></tr><tr><td>Total</td><td>6</td><td></td><td>1</td><td></td><td>7</td></tr></tbody></table>

the supplemental material). Out of the 331 isolates investigated, 258 (77.9%) were classified into category 0 (low level of resistance). All of the remaining isolates (resistance scores of 1, 2, and 3 [$n = 73$]) belonged to Kp1 and were isolated mostly from humans ($n = 69$; 94.5%) and the hospital setting ($n = 64$; 87.7%). All isolates carrying an ESBL gene ($n = 61$; 18.4%) were recovered from human isolates associated with nosocomial infections, except for one collected from a cat, in agreement with the results of phenotypic antibiotic susceptibility testing. The ESBL genes were *bla*<sub>CTX-M-15</sub> ($n = 59$) and *bla*<sub>SHV-12</sub> ($n = 3$). One nosocomial isolate displayed *bla*<sub>CTX-M-15</sub> and *bla*<sub>SHV-12</sub>. With regard to STs including ≥5 isolates, ST17 (15/16), ST147 (7/9), ST152 (5/8), and ST307 (12/12) had high absolute rates of 3rd-generation cephalosporin resistance (3GCR) of over 78%. MDR high-risk clones were associated with a high number of resistance genes, with 11 out of these 13 STs displaying a mean of 6 or more genes.

Two of the three carbapenemase-producing *K. pneumoniae* isolates were sequenced (Kp018 and Kp020); both harbored a *bla*<sub>KPC-2</sub> gene, in agreement with PCR typing results (17). They were isolated in the hospital setting. These two isolates were closely related (25 cgMLST allelic mismatches, assigned to the common health care-associated clone ST258). They were classified into the highest-resistance category, harboring the carbapenemase gene *bla*<sub>KPC-2</sub> and an alteration of the *mgrB* gene through total (Kp018) or partial (Kp020) deletion, known to confer colistin resistance. They were resistant to colistin (MIC > 4 mg/L). There was no other mechanism of colistin resistance detected in the other isolates, including the acquisition of the *mcr* gene.

Plasmid replicons were found in 67.3% (103/153) of human isolates, 67.8% (38/56) of isolates from pets, 48.3% (29/60) of isolates from food-producing animals, and 45.2% (28/62) of environmental isolates. IncFII (10.6%), IncFIA(HI1) (10.6%), ColRNAI (11.2%), IncR (53.0%), and IncFIB(K) (55.0%) were the most frequently recovered plasmid replicons, with IncFIA (HI1), IncFIB(K), IncFII, and IncR being present in all sources (Table S6).
<figure><img src="image_7.png" /><figcaption>**FIG 3** Unrooted phylogenetic trees of the genomes of the six most represented sequence types (STs) recovered in Guadeloupe (French West Indies). All these STs belong to *K. pneumoniae sensu stricto* (Kp1). Isolates are colored according to their source (i.e., human in blue and animal in orange). Using a threshold of <10 SNPs, single strains are framed with dotted lines.</figcaption></figure>

**Virulence genes.** The prevalences of operons coding for acquired siderophores were 24.2% for yersiniabactin (80/331; 3/80 were incomplete), 6.0% for aerobactin (20/331), and 5.7% for salmochelin (19/331). The genotoxin colibactin loci were present in 4.5% (15/331; 1/15 was incomplete) of the isolates (see Table S7 in the supplemental material). Virulence genes were recovered mainly in human and pet isolates. By comparing human and pet isolates, yersiniabactin (43.8% versus 17.9% [*P* = 0.001]) and colibactin (9.8% versus 0% [*P* = 0.013]) frequencies were significantly higher in humans, but no significant difference was observed for aerobactin (9.2% versus 7.1% [*P* = 0.785]) and salmochelin (5.2% versus 1.9% [*P* = 0.450]). Five isolates from pigs displayed virulence genes, two with yersiniabactin and three with aerobactin, as did one from soil (yersiniabactin). The three aerobactin-positive isolates from pigs carried an *iuc3* aerobactin gene. They came from different farms and belonged to three different genetic backgrounds (ST29, ST432, and ST827). In addition, the *iuc3* gene was recovered on a contig predicted to be plasmidic by combining the MOB-recon and PlasFlow tools. For two out of the three isolates, *iuc3* and an IncFIB(K) replicon were present in the same contig. It should be noted that a presumed plasmid carrying *iuc3* was recovered in one isolate from a dog. The isolate carried an IncFIB(K) replicon but was not recovered on the *iuc3*-positive contig.

In all, a total of 22 isolates (6.6%) were defined as hypervirulent: 15 were associated with community (n = 12; 21.4%) and nosocomial (n = 3; 3.1%) human infections, 4 (8.9%) were associated with dogs, and 3 (15%) were associated with pigs. The genomic features of HvKp isolates are displayed in Table 2. ST23 (subclade I) was the ST most frequently recovered in human isolates (33.3% of human HvKp isolates). Nine isolates (34.6%) harbored the regulator of mucoid phenotype gene *rmpA*, whereas
<table><caption>**TABLE 2** Genomic features of 22 hypervirulent *K. pneumoniae* species complex isolates collected in Guadeloupe (French West Indies) according to sourceª</caption><thead><tr><th rowspan="3">Feature</th><th colspan="5">No. (%) of isolates with feature</th></tr><tr><th colspan="2">Humans</th><th rowspan="2">Dogs<br>(n = 4)</th><th rowspan="2">Pigs<br>(n = 3)</th><th rowspan="2">Total<br>(n = 22)</th></tr><tr><th>Community<br>(n = 12)</th><th>Nosocomial<br>(n = 3)</th></tr></thead><tbody><tr><td>Phylogroup</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>K. pneumoniae (Kp1)</td><td>12 (100)</td><td>3 (100)</td><td>3 (75)</td><td>3 (100)</td><td>21 (95.5)</td></tr><tr><td>K. quasipneumoniae subsp.<br>quasipneumoniae (Kp2)</td><td>0 (0)</td><td>0 (0)</td><td>1 (25)</td><td>0 (0)</td><td>1 (4.5)</td></tr><tr><td>Virulence gene</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Colibactin (clb)</td><td>10 (83.3)</td><td>2 (66.7)</td><td>0 (0)</td><td>0 (0)</td><td>12 (54.5)</td></tr><tr><td>Aerobactin (iuc)</td><td>12 (100)</td><td>3 (100)</td><td>2 (50)</td><td>3 (100)</td><td>20 (90.9)</td></tr><tr><td>Salmochelin (iro)</td><td>12 (100)</td><td>3 (100)</td><td>4 (100)</td><td>0 (0)</td><td>19 (86.4)</td></tr><tr><td>rmpA</td><td>5 (41.7)</td><td>2 (66.7)</td><td>2 (50)</td><td>0 (0)</td><td>9 (40.9)</td></tr><tr><td>rmpA2</td><td>0 (0)</td><td>0 (0)</td><td>1 (25)</td><td>0 (0)</td><td>1 (4.5)</td></tr><tr><td>CG and associated ST</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CG5 ST5</td><td>0 (0)</td><td>0 (0)</td><td>1 (25)</td><td>0 (0)</td><td>1 (4.5)</td></tr><tr><td>CG23 ST23</td><td>3 (25)</td><td>2 (66.7)</td><td>0 (0)</td><td>0 (0)</td><td>5 (22.7)</td></tr><tr><td>CG29 ST29</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>1 (33.3)</td><td>1 (4.5)</td></tr><tr><td>CG35 ST5750</td><td>0 (0)</td><td>1 (33.3)</td><td>0 (0)</td><td>0 (0)</td><td>1 (4.5)</td></tr><tr><td>CG60 ST60</td><td>0 (0)</td><td>0 (0)</td><td>1 (25)</td><td>0 (0)</td><td>1 (4.5)</td></tr><tr><td>CG65 ST3253</td><td>2 (16.7)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>2 (9.1)</td></tr><tr><td>CG66 ST66</td><td>2 (16.7)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>2 (9.1)</td></tr><tr><td>CG66 ST3252</td><td>1 (8.3)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>1 (4.5)</td></tr><tr><td>CG86 ST86</td><td>1 (8.3)</td><td>0 (0)</td><td>1 (25)</td><td>0 (0)</td><td>2 (9.1)</td></tr><tr><td>CG260 ST260</td><td>1 (8.3)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>1 (4.5)</td></tr><tr><td>CG380 ST380</td><td>2 (16.7)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>2 (9.1)</td></tr><tr><td>CG432 ST432</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>1 (33.3)</td><td>1 (4.5)</td></tr><tr><td>CG446 ST446</td><td>0 (0)</td><td>0 (0)</td><td>1 (25)</td><td>0 (0)</td><td>1 (4.5)</td></tr><tr><td>CG827 ST827</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>1 (33.3)</td><td>1 (4.5)</td></tr><tr><td>Capsular type</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>K1</td><td>3 (25)</td><td>2 (66.7)</td><td>1 (25)</td><td>0 (0)</td><td>6 (27.3)</td></tr><tr><td>K2</td><td>7 (58.3)</td><td>1 (33.3)</td><td>1 (25)</td><td>0 (0)</td><td>9 (40.9)</td></tr><tr><td>Other</td><td>2 (16.7)</td><td>0 (0)</td><td>2 (50)</td><td>3 (100)</td><td>7 (31.8)</td></tr></tbody></table>

yersiniabactin, colibactin, and aerobactin operons were present together in 11 (50.0%) isolates. All isolates belonged to *K. pneumoniae*, except for one Kp2 isolate associated with fecal carriage in a dog. These hypervirulence-associated operons were observed mostly in isolates associated with human infections (15/22; 68.2%). The remaining hypervirulent isolates were associated with fecal carriage in dogs (n = 4) and pigs (n = 3). Human infections corresponded mostly to community-acquired infections (12/15; 80%) and were associated mostly with liver abscess (n = 5) and meningitis (n = 2). The human hypervirulent isolates were associated mostly with STs typically associated with capsular serotypes K1 (n = 5; 22.7%) and K2 (n = 9; 40.9%), with the most frequently recovered genetic backgrounds being ST23 (n = 5), ST66 (n = 2), ST86 (n = 2), ST380 (n = 2), and ST3253 (n = 2). ST23 isolates belonged to the CG23-I lineage (data not shown) according to Lam et al. (18) and carried a yersiniabactin-encoding mobile element (ICEKp10), including genes coding for yersiniabactin and colibactin. All HvKp isolates exhibited the wild-type AMR phenotype.

**Virulence genes and MDR convergence.** Virulence and antibiotic resistance elements were always observed in distinct isolates, except for one human *K. pneumoniae* isolate (KP586) where a convergence of virulence and MDR was observed. It was assigned to ST392, belonging to the successful MDR clonal group CG147. It was community acquired and displayed the aerobactin virulence operon, *bla*<sub>CTX-M-15</sub>, and numerous other resistance genes (*aadA2*, *aac3-11a*, *bla*<sub>OXA-1</sub>, *bla*<sub>SHV-67</sub>, *bla*<sub>TEM-30</sub>, *catB4*, *dfrA14*, *dfrA32*, *tetA*, *sull*, and *sull*).

<footnote>ªCG, clonal group; ST, sequence type.</footnote>
**Genome-wide association study.** We performed a genome-wide association study (GWAS) on Kp1 isolates in order to identify genes associated with human sources. A genomic region was significantly associated with virulence and host specificity (human). This region includes a cluster of genes located in an integrative conjugative element that mobilizes the *ybt* locus, which encodes the biosynthesis of the siderophore yersiniabactin and its receptor (see Table S8 in the supplemental material). We also found 5 genes significantly associated with human infections encoding proteins with catalytic activity (acornitate hydratase, cytochrome *bo*₃ ubiquinol oxidase, acetoin:2,6-dichlorophenolindophenol oxidoreductase, and dihydrolipoamide dehydrogenase) and 1 encoding a putative methyltransferase. Two genes were significantly associated with nonhuman isolates, including a gene encoding a core component of a type VI secretion system (T6SS) (T6SS baseplate subunit TssK).

DISCUSSION

The increasing level of AMR is a major health hazard for humans and animals (19). Tackling AMR transmission requires investigations of the nonclinical reservoirs and their relative contribution to human infections through the so-called One Health approach. WGS combined with phylogenetic analysis is a powerful approach to provide detailed insights into bacterial transmission dynamics. Despite the urgent public health threat represented by the KpSC, only a few studies so far have identified their genomic features using WGS from nonhuman sources and compared these features with those obtained from contemporaneous and colocalized clinical isolates (7–11). We observed a limited overlap of STs between clinical isolates and local nonhuman isolates, consistent with the results of previous studies (here, 12.7% of human STs; range in previous studies, 5 to 15%). In a large survey of an Italian contemporary KpSC collection from a well-defined geographical region, direct transmission from animal or environmental reservoirs represents a small fraction (<1%) of human infections (20). This highlights the difficulty in identifying direct transmission events for a pathogen characterized by its high genetic diversity without a dominance of specific successful lineages. However, we identify evidence of sporadic transmissions between animals and humans in our set of isolates, such as the presence of 3 hospital-acquired human MDR high-risk *K. pneumoniae* clones (ST17, ST37, and ST45) in pets and 1 successful human hypervirulent *K. pneumoniae* clone (ST86) in a dog. It should be noted that all but one of the STs (one bovine) were shared only between pets and humans, highlighting the potential risk of companion animals in the transmission of KpSC isolates to humans and vice versa, as previously described (21, 22). The presence of the same isolate (<10 SNPs) assigned to ST45 in a dog and in a human collected 1 month apart highlights the importance of the application of basic hygiene rules for contact with companion animals.

Our results illustrate the efficiency of the genomic approach to distinguish epidemiologically related isolates from unrelated ones within hospitals. Genetically closely related isolates were recovered during the 24-month study across different units of the hospital, suggesting an environmental reservoir and long-term transmission. Most nosocomial isolates belonged to high-risk MDR genetic backgrounds, namely, ST11, ST17, ST45, and ST307. These clones have emerged as important vehicles for the worldwide dissemination of AMR determinants (23, 24), including in the Caribbean islands (16).

Although members of the KpSC can be found in a large variety of ecological niches (25), knowledge of the prevalence and distribution of species and phylogroups belonging to this complex according to the source is limited due to the lack of large-scale systematic sampling efforts. Our findings show that the KpSC is ubiquitous, as shown previously in the environment (25), food (26), and the intestines of mammals (21, 27, 28). To the best of our knowledge, dogs were significantly associated for the first time with a higher risk of KpSC carriage than cats. Consistent with the results of previous studies (8, 20), species and phylogroups were not randomly distributed. In particular, *K. pneumoniae sensu stricto* (Kp1) and Kp3 (*K. variicola* subsp. *variicola*) were clearly separated according to ecological niche:
a high prevalence of Kp1 and a low prevalence of Kp3 in humans and domestic and food-
producing animals (except for bovines) and the contrary in vegetables and the environ-
ment. This may imply ecological barriers that limit the spread of clones and antibiotic
resistance genes. However, the distributions of species/phylogroups within the same
source were not completely consistent across studies (8, 20), highlighting the need for fur-
ther work with the inclusion of isolates from wider environmental and animal sources from
various geographical areas. For example, we observed a high prevalence of Kp3 among
our bovine isolates, probably due to transient flora related to the consumption of raw
plants rather than a specific adaptation to colonize their intestine, consistent with the
results of a previous study (8) but not with the results of another one (20).

As KpSC isolates were largely recovered in human, animal, and environmental reser-
voirs in Guadeloupe, we hypothesize that the KpSC could be a major vector for the
amplification and spread of antibiotic resistance genes due to its abilities to move
between ecological niches, capture and maintain plasmids carrying AMR genes for a
long time, and transfer plasmids within KpSC members but also to other important
Gram-negative bacteria (6). However, a high level of resistance to antibiotics was rarely
found in isolates collected outside the hospital setting but also in isolates of species
other than *K. pneumoniae*, illustrated by the almost exclusive presence of ESBL and car-
bapenemase genes in human *K. pneumoniae* isolates, in agreement with the results of
previous studies in Guadeloupe (29–31). Nevertheless, further studies using a shotgun
metagenomics approach are needed to access the whole resistome in different ecosys-
tems. Despite the risk of the occasional emergence of novel resistance mechanisms in
KpSC isolates from environmental sources, our findings strongly support that the noso-
comial setting is central to KpSC resistance dissemination, as observed in Europe for
carbapenemase-producing *K. pneumoniae* (32), and that KpSC resistance circulates less
frequently between the different compartments.

The *bla*<sub>CTX-M-15</sub> ESBL gene, first detected in 1999 in India, was the most widely dis-
tributed ESBL gene in our set of isolates, in agreement with the worldwide situation.
AMR genes and plasmids are often associated with certain *K. pneumoniae* genetic line-
ages, as highlighted by the success of *K. pneumoniae* ST258 being intricately linked
with *bla*<sub>KPC</sub> (33). Further work will be done to characterize the plasmids carrying ESBL
genes using Oxford Nanopore technologies in order to study the dissemination of plas-
mids and their ESBL genes within the hospital setting.

High levels of virulence also tended to be rare in species other than *K. pneumoniae*. Most of the human HvKp isolates belonged to *K. pneumoniae* ST23 and more precisely to sublineage CG23-1, which emerged in approximately 1928 following the acquisition of ICEKp10 (encoding yersiniabactin and colibactin) and then spread worldwide (18). The uncommon ST66 genetic background (34) was also identified. A reservoir of HvKp was not found outside humans, but dogs could be an important link in the chain of the transmission of this pathogen (9% of hypervirulent isolates were recovered from this species), as highlighted by the presence of ST86 associated with a dog and a human meningitis case. Aerobactin (*iuc operon*) was present in 90% of our HvKp isolates. It is considered a critical siderophore system of HvKp as it contributes predominantly to hypervirulence in laboratory experiments and mouse models of disease, while the inactivation of other siderophore systems has minimal effects (35, 36). Surprisingly, a high frequency of aerobactin (15%) in *K. pneumoniae* isolates from pigs was observed, as previously described in Thailand (37) and Germany (11), probably reflecting an adaptation conferred by this siderophore to porcine hosts. Although our isolates were isolated from pigs from different farms and belonged to different STs, they harbored an *iuc3* gene carried by an IncFIB(K) plasmid (for at least two isolates), as observed in Germany (11). These observations suggest that successful IncFIB(K)/*iuc3*-carrying plasmids have spread across wide geographical distances and occur in different *K. pneumoniae* lineages associated with domestic pigs. The potential risk to animal and human health should also be investigated. Unsurprisingly, a significant association was observed between human infections and the *ybt* locus encoding yersiniabactin,
the most common *K. pneumoniae* high-virulence factor, present in around one-third of clinical isolates (8, 38). The emergence of potentially high-risk MDR and hypervirulent lineages within the hospital setting in Guadeloupe should be monitored, as illustrated by the acquisition of virulence genes in an MDR genetic background (ST392), even if still rarely observed. For the other five genes found to be significantly associated with human infections by GWASs, we did not find any explanation.

Our study represents a large contemporaneous and colocalized sampling and sequencing effort on an island characterized by its small population and area and known to be a hot spot for the spread of health care-associated MDR *K. pneumoniae*. Despite its insular character, which is expected to promote mainly transmissions with local clones and a restricted contribution from the outside, we found limited evidence for direct transmission between human and nonhuman sources (animals and the environment). In contrast, the nosocomial context seems to be a cornerstone of the dissemination of MDR clones and antibiotic resistance genes.

## MATERIALS AND METHODS

**Fecal samples from healthy food-producing animals and pets.** The study design and methods for selecting healthy food-producing animals were described previously (39). Briefly, between January 2018 and May 2019, fecal samples from pigs and beef cattle were collected randomly just after excretion. Fecal material from 17 hen houses (representing 53,000 poultry) was sampled by walking on litter approximately 100 m around a flock in boot socks. In all, the animals originated from 11 pig farms, 8 beef cattle farms, and 9 poultry farms distributed throughout the island and from the only slaughterhouse in Guadeloupe (34 additional farms) for cattle and pigs. Sixty-four percent of farmers declared antibiotic use during the previous 1 year for curative treatment, with the most commonly used antibiotic being tetracycline (69.0%). No ethics committee approval was necessary as no invasive procedure was conducted on live animals according to French national law for the protection of animals (no. 2013-118), which reproduces European directive 2010/63/EU on the protection of animals used for experimental and other scientific purposes.

For pets, from June to September 2019, a single rectal swab was taken from dogs and cats. The animals were included from the main animal shelter of Guadeloupe and seven veterinary clinics located throughout the territory, among animals sent for preventive health services, vaccination, or medical consultation. With regard to the identification of risk factors for the fecal carriage of ESBL-producing KpSC, pets with clinical signs of diarrhea and/or with antibiotic treatment in the previous month were excluded. The information collected for each animal included age, place of residence, general health, and lifestyle (indoors or wandering free outdoors and close contact with other animals or not). The project was approved by the Committee for Ethics in Animal Experiments of the French West Indies and Guyana (reference no. HC_2020_1).

**Fresh fruits, fresh vegetables, flowering plants, aromatic herbs, and water and soil samples.** Locally produced fresh vegetables, flowering plants, fruits, and aromatic herbs were collected aseptically at four local markets during eight campaigns from January to June 2018. The collected samples spanned producers from 14 municipalities throughout Guadeloupe. Data related to the market and the farm of origin were recorded. When the farm of origin was unknown, multiple vendors in different parts of the market were selected to minimize the likelihood that samples came from the same source.

Raw water samples were collected during the same period at 29 catchment points in 11 municipalities, in partnership with the regional health agency and the hygiene laboratory of the Pasteur Institute of Guadeloupe. These samples corresponded to drinking water before treatment.

From October to December 2019, soil was sampled near rivers and natural ponds: the surface layer (0 to 10 cm) was collected after the removal of plants, pebbles, and conspicuous roots. At each site, two or three samples of soil were taken at least 3 m apart. One sample of water from rivers or natural ponds was collected in proximity. Sampling sites were classified by Q-GIS software into two groups according to their degree of anthropogenic pressure: (i) wilderness with no human presence or countryside with limited human activities, (ii) human-perturbed landscapes with a matrix of agriculture and livestock activities, and (iii) urban and suburban areas with high levels of human activity.

All samples were kept at 2°C to 8°C and processed at the laboratory within 24 h.

**Clinical isolates.** Between January 2018 and December 2019, 279 contemporaneous KpSC isolates were collected from patients admitted to the University Hospital of Guadeloupe, a 900-bed teaching hospital. Isolation was performed as part of the routine activity of the hospital bacteriological diagnostic laboratory. All presumptive HvKp isolates, defined as KpSC isolates associated with community-acquired monomicrobial liver abscess or other monomicrobial invasions of normally sterile sites (e.g., meningitis), were included in the same period. The following metadata were anonymously recorded: date of hospital admission, ward of hospitalization, date and site of sampling, and antimicrobial susceptibility testing results. Isolates were considered to be community acquired if they were recovered by culture from a sample obtained within 48 h after admission in a patient with no risk factors for nosocomial acquisition in the previous year, namely, hospitalization or surgery, the use of an indwelling catheter or a percutaneous device, or frequent exposure to health care facilities for an underlying chronic disorder. All other
isolates were considered to be hospital acquired. The study protocol was approved by the ethic commit-
tee of the University Hospital of Guadeloupe (reference no. A5_19_12_05_TRAMID).

**K. pneumoniae** species complex isolation and antimicrobial susceptibility testing. Fruits, vegeta-
bles, flowering plants, and aromatic herbs after mixing and soil and fecal samples (stool and boot sock
samples and endorectal swabs) were inoculated into Luria-Bertani (LB) broth with amoxicillin (10-mg/L
final concentration). For water samples, 100 mL of serially diluted samples was filtered through a 0.45-
µm membrane filter (Millipore, Guyancourt, France), and the membranes were placed into 9 mL of LB
broth with amoxicillin (10-mg/L final concentration). After incubation for 18 h at 37°C, 100 µL of the
enrichment culture was plated onto two selective media, Simmons citrate agar (Becton, Dickinson, USA)
with 1% inositol (SCAI) medium agar plates for KpSC detection and chromogenic agar (CCA) with cef-
triaxone at 4 mg/L (CHROMagar) for the detection of 3GCR KpSC. Presumptive Enterobacteriaceae colo-
nies on selective SCAI medium (large, yellow, glossy colonies) and selective CCA with ceftriaxone (pink
colonies), corresponding to oxidase-negative and Gram-negative bacilli, were isolated randomly and
identified by matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF
MS) on an Axima Performance system (Shimadzu Corp., Japan).

Three colonies were identified randomly for each identical morphology. Susceptibility to amoxicillin (10 µg), amoxicillin-clavulanic acid (20 µg/10 µg), ticarcillin (75 µg), cefotaxime (5 µg), ceftazidime (10 µg), cefepime (30 µg), cefoxitin (30 µg), aztreonam (30 µg), ertapenem (10 µg), gentamicin (15 µg), amikacin (30 µg), trimethoprim-sulfamethoxazole (1.25/23.75 µg), nalidixic acid (30 µg), and ciprofloxacin (5 µg) was tested by the disk diffusion method on Mueller-Hinton agar (Bio-Rad, Marnes-la-Coquette, France), and the production of ESBL was detected by the double-disk synergy test, according to 2020 guidelines of CA-SFM/EUCAST (https://www.sfm-microbiologie.org/2020/10/02/casfm-eucast-v1-2-octobre-2020/). Isolates with a resistant or intermediate phenotype were classified together for analysis. Growth inhibition diameters were measured with the Adagio automated system (Bio-Rad, Marnes-la-Coquette, France). Susceptibility to colistin was determined using a MicroScan Walkaway Plus system (Beckman Coulter, USA). Resistant isolates were defined by an MIC of >4 mg/L. If more than one KpSC isolate with the same antibiotic susceptibility pattern was isolated from the same sample, only the first one was analyzed. An MDR KpSC isolate was defined as an isolate resistant to three or more antimicrobial classes (40).

**DNA extraction, K. pneumoniae species complex identification, and carbapenemase resistance gene screening.** DNA was extracted with a DNA minikit (Qiagen, Germany). A real-time PCR method based on specific sets of primers and probes was used to identify *K. pneumoniae* isolates to the species and phylogroup levels. The protocols are provided in Text S1 in the supplemental material. All isolates were handled using this method, except for those associated with nosocomial infections, of which about half were randomly selected due to the large number of isolates collected. For isolates with decreased susceptibility to ertapenem, carbapenemase genes were searched for using PCR amplification according to methods described in a previous study (17).

**Genome sequencing and data analysis.** WGS was carried out at the Plateforme de Microbiologie Mutualisée de l'Institut Pasteur (Paris, France). Reads were trimmed and filtered with AlienTrimmer software (41), yielding a mean estimated coverage of 86-fold. Genomic assemblies were performed using SPAdes software (42), and the quality of the assembly was evaluated using QUAST software (43). Genomes with a cumulative size of contigs of >6 Mb (expected size of ~5.5 Mb) or a number of contigs of >500 were discarded, as we suspected the presence of multiple clones. The mean N<sub>50</sub> was 197,864 bp.

KpSC phylogrouping, sequence type (ST) assignment, and antibiotic resistance gene detection were
performed using Kleborate (44). This tool classified isolates according to the content of resistance genes
as follows: 0 for no ESBL and no carbapenemase, 1 for ESBL positive, 2 for carbapenemase positive, and
3 for carbapenemase with colistin resistance (44). Plasmid replicons were identified using the
PlasmidFinder database available from ABRicate software, using a minimum coverage and a minimum
identity of 90% (45). Kleborate was also used to search for yersiniabactin, colibactin, aerobactin, and sal-
mochelin operons and the presence of *rmpA* and *rmpA2* and to predict capsular types (44). BIGSdb
(https://bigsdb.pasteur.fr/klebsiella/) was used to check for the presence of intact *iucABCD-iutA* (aerobac-
tin) and *iroBCDN* (salmochelin) operons. To study the genetic support of *iuc* genes in more detail, the
corresponding contigs were predicted to be plasmid or chromosomally associated by combining 2 dif-
ferent software tools (MOB-recon and PlasFlow) with default thresholds (46, 47).

The definitions described previously by Huynh et al. (48) for hypervirulent isolates and successful clones were used. Hypervirulent isolates were defined as isolates harboring at least one of the *rmpA* and *rmpA2* genes and/or at least one complete operon among *iucABCD-iutA* and *iroBCDN*. Successful clones were defined as those belonging to an ST represented at least 10 times in NCBI genomes and mentioned in the title or abstract of at least five publications in NCBI PubMed ("Klebsiella" + "pneumoniae" + "STxxx").

**Pangenome and phylogenetic analyses.** Filtered raw contigs were assigned to the chromosome or plasmid using MOB-suite (47), and chromosomal contigs were then ordered and oriented using RaGOO (49) with a publicly available complete assembled genome available for each *Klebsiella* species, according to the species assigned previously by Kleborate (44). Genome annotation was performed using Prokka (50), and pangenome analysis was performed using Roary software (51). Pangenome matrix representation was done using the Roary plots python utility (https://github.com/sanger-pathogens/Roary/tree/master/contrib/roary Plots).

Based on the core-genome alignment provided by Roary and after the identification of recombinant
regions and the reduction of the alignment using ClonalFrameML (52), a global phylogenetic tree was
constructed using RAxML (53) and visualized using iTOL (54). Estimation of the number of single nucleo-
tide polymorphisms (SNPs) between isolates assigned to an identical ST was performed using the
PathogenWatch platform (https://pathogen.watch), an online global database for genomic surveillance of KpSC isolates (55). Contigs were uploaded, and collections were created separately by ST, in order to generate SNP difference matrix files on the one hand and Newick files to be visualized as unrooted trees using iTOL on the other hand. Two isolates were considered clonal when the number of SNPs between them was <10. Additional *Klebsiella* strains obtained from other Caribbean islands (16) were added to the collection. For this purpose, raw FASTQ files were retrieved from the SRA and preliminarily assembled using Unicycler to then be uploaded to PathogenWatch.

The pangenome matrix from Roary consists of gene presence or absence for each genome. It was used as the input for Scoary V1.6.16 (56) in order to search for genes associated with humans or other sources. Due to the bias that might be introduced by the high predominance of Kp1 in human isolates, we focused our pangenome-wide association studies (pan-GWASs) on Kp1 only. Genes returning a Bonferroni-corrected *P* value of ≤0.05 were considered to be significantly present/absent and were further investigated.

**Statistical analyses.** Results were expressed as numbers and frequencies. In bivariate analyses, the χ² test (or Fisher's exact test when appropriate) was used to compare categorical data between groups. A logistic regression model was performed to identify factors associated with the presence of KpSC isolates and to calculate crude and adjusted odds ratios and their 95% confidence intervals. Factors with a *P* value of <0.20 in the bivariate analysis were retained for the multivariate analysis. For all tests, we considered a *P* value of <0.05 to be significant. Statistical analyses were performed using SPSS (V21; IBM SPSS Statistics, Chicago, IL).

**Data availability.** Reads were deposited in the NCBI SRA public archives under BioProject accession no. PRJNA778230.

## SUPPLEMENTAL MATERIAL

Supplemental material is available online only.

**SUPPLEMENTAL FILE 1**, PDF file, 1.4 MB.

**SUPPLEMENTAL FILE 2**, XLSX file, 0.03 MB.

**SUPPLEMENTAL FILE 3**, XLS file, 0.03 MB.

**SUPPLEMENTAL FILE 4**, XLSX file, 0.2 MB.

## ACKNOWLEDGMENTS

We thank the Plateforme de Microbiologie Mutualisée (P2M) of the Institut Pasteur for Illumina sequencing. We thank all of the students, technicians, and clinicians involved in this work at the Pasteur Institute of Guadeloupe and the University Hospital Center of Guadeloupe.

Sylvain Brisse and Sébastien Breurec conceived and designed the study. Gaëlle Gruel, Matthieu Pot, Stéphanie Guyomard-Rabenirina, Sylvaine Bastian, Moana Gelu-Simeon, Séverine Ferdinand, Marc Valette, and Sébastien Breurec collected biological samples, isolates, and epidemiological data. Elodie Barbier, Carla Rodrigues, Pascal Piveteau, and Sylvain Brisse provided laboratory support and validated the analysis of the RT-PCR assays. Virginie Passet and Sylvain Brisse performed the MLST and cgMLST curation/analyses. Alexis Dereeper, Pascal Piveteau, Yann Reynaud, Sylvain Brisse, and Sébastien Breurec analyzed the data. Benoit Tressieres performed the statistical analyses. Alexis Dereeper and Sébastien Breurec wrote the initial version of the manuscript. All authors provided input to the manuscript and reviewed the final version.

We declare that there are no conflicts of interest.

This work was supported by a FEDER grant financed by the European Union and Guadeloupe Region (Programme Opérationnel FEDER-Guadeloupe-Conseil Régional 2014–2020, grant no. 2015-FED-192). Carla Rodrigues was supported financially by the MedVetKlebs project, a component of the One Health European Joint Programme (EJP), which has received funding from the European Union's Horizon 2020 Research and Innovation Programme under grant agreement no. 773830, and by a Pasteur-Roux fellowship from the Institut Pasteur.

## REFERENCES

1. Cassini A, Högberg LD, Plachouras D, Quattrocchi A, Hoxha A, Simonsen GS, Colomb-Cotinat M, Kretzschmar ME, Devleesschauwer B, Cecchini M, Ouakrim DA, Oliveira TC, Struelens MJ, Suetens C, Monnet DL, Burden of AMR Collaborative Group. 2019. Attributable deaths and disability-

adjusted life-years caused by infections with antibiotic-resistant bacteria in the EU and the European Economic Area in 2015: a population-level modelling analysis. Lancet Infect Dis 19:56–66. https://doi.org/10.1016/S1473-3099(18)30605-4.
2. Breurec S, Bouchiat C, Sire JM, Moquet O, Bercion R, Cisse MF, Glaser P, Ndiaye O, Ka S, Salord H, Seck A, Sy HS, Michel R, Garin B. 2016. High third-generation cephalosporin resistant Enterobacteriaceae prevalence rate among neonatal infections in Dakar, Senegal. BMC Infect Dis 16:587. https://doi.org/10.1186/s12879-016-1935-y.

3. Shon AS, Bajwa RP, Russo TA. 2013. Hypervirulent (hypermucoviscous) Klebsiella pneumoniae: a new and dangerous breed. Virulence 4:107–118. https://doi.org/10.4161/viru.22718.

4. Wyres KL, Wick RR, Judd LM, Froumine R, Tokolyi A, Gorrie CL, Lam MMC, Duchene S, Jenney A, Holt KE. 2019. Distinct evolutionary dynamics of horizontal gene transfer in drug resistant and virulent clones of Klebsiella pneumoniae. PLoS Genet 15:e1008114. https://doi.org/10.1371/journal.pgen.1008114.

5. Rodrigues C, Passet V, Rakotondrasoa A, Diallo TA, Criscuolo A, Brisse S. 2019. Description of Klebsiella africanensis sp. nov., Klebsiella variicola subsp. tropicalensis subsp. nov. and Klebsiella variicola subsp. variicola subsp. nov. Res Microbiol 170:165–170. https://doi.org/10.1016/j.resmic.2019.02.003.

6. Wyres KL, Holt KE. 2018. Klebsiella pneumoniae as a key trafficker of drug resistance genes from environmental to clinically important bacteria. Curr Opin Microbiol 45:131–139. https://doi.org/10.1016/j.mib.2018.04.004.

7. Davis GS, Waits K, Nordstrom L, Weaver B, Aziz M, Gauld L, Grande H, Bigler R, Horwinski J, Porter S, Stegger M, Johnson JR, Liu CM, Price LB. 2015. Intermingled Klebsiella pneumoniae populations between retail meats and human urinary tract infections. Clin Infect Dis 61:892–899. https://doi.org/10.1093/cid/civ428.

8. Holt KE, Wertheim H, Zadoks RN, Baker S, Whitehouse CA, Dance D, Jenney A, Connor TR, Hsu LY, Severin J, Brisse S, Cao H, Wilksch J, Gorrie C, Schultz MB, Edwards DJ, Nguyen KV, Nguyen TV, Dao TT, Mensink M, Minh VL, Nhu NT, Schultsz C, Kuntaman K, Newton PN, Moore CE, Strugnell RA, Thomson NR. 2015. Genomic analysis of diversity, population structure, virulence, and antimicrobial resistance in Klebsiella pneumoniae, an urgent threat to public health. Proc Natl Acad Sci U S A 112: E3574–E3581. https://doi.org/10.1073/pnas.1501049112.

9. Ludden C, Moradigaravand D, Jamrozy D, Gouliouris T, Blane B, Naydenova P, Hernandez-Garcia J, Wood P, Hadjirin N, Radakovic M, Crawley C, Brown NM, Holmes M, Parkhill J, Peacock SJ. 2020. A One Health study of the genetic relatedness of Klebsiella pneumoniae and their mobile elements in the east of England. Clin Infect Dis 70:219–226. https://doi.org/10.1093/cid/ciz174.

10. Runcharoen C, Moradigaravand D, Blane B, Paksanont S, Thammachote J, Anun S, Parkhill J, Chantratita N, Peacock SJ. 2017. Whole genome sequencing reveals high-resolution epidemiological links between clinical and environmental Klebsiella pneumoniae. Genome Med 9:6. https://doi.org/10.1186/s13073-017-0397-1.

11. Klaper K, Hammerl JA, Rau J, Pfeifer Y, Werner G. 2021. Genome-based analysis of Klebsiella spp. isolates from animals and food products in Germany, 2013–2017. Pathogens 10:573. https://doi.org/10.3390/pathogens10050573.

12. Piednoir P, Clarac U, Rolle A, Bastian S, Gruel G, Martino F, Mehdaoui H, Valette M, Breurec S, Carles M. 2020. Spontaneous community-acquired bacterial meningitis in adults admitted to the intensive care units in the Caribbean French West Indies: unusual prevalence of Klebsiella pneumoniae. Int J Infect Dis 100:473–475. https://doi.org/10.1016/j.ijid.2020.09.1420.

13. Bastian S, Nordmann P, Creton E, Malpote E, Thiery G, Martino F, Breurec S, Dortet L. 2015. First case of NDM-1 producing Klebsiella pneumoniae in Caribbean islands. Int J Infect Dis 34:53–54. https://doi.org/10.1016/j.ijid.2015.03.002.

14. Arnaud I, Maugat S, Jarlier V, Astagneau P, National Early Warning, Investigation and Surveillance of Healthcare-Associated Infections Network (RAISIN)/Multidrug Resistance Study Group. 2015. Ongoing increasing temporal and geographical trends of the incidence of extended-spectrum beta-lactamase-producing Enterobacteriaceae infections in France, 2009 to 2013. Euro Surveill 20:30014. https://doi.org/10.2807/1560-7917.ES.2015.20.36.30014.

15. Le Terrier C, Vinetti M, Bonjean P, Richard R, Jarrige B, Pons B, Madeux B, Piednoir P, Ardisson F, Elie E, Martino F, Valette M, Ollier E, Breurec S, Carles M, Thiery G. 2021. Impact of a restrictive antibiotic policy on the acquisition of extended-spectrum beta-lactamase-producing Enterobacteriaceae in an endemic region: a before-and-after, propensity-matched cohort study in a Caribbean intensive care unit. Crit Care 25:261. https://doi.org/10.1186/s13054-021-03660-z.

16. Heinz E, Brindle R, Morgan-McCalla A, Peters K, Thomson NR. 2019. Caribbean multi-centre study of Klebsiella pneumoniae: whole-genome sequencing, antimicrobial resistance and virulence factors. Microb Genom 5:e000266. https://doi.org/10.1099/mgen.0.000266.

17. Poirel L, Walsh TR, Cuvillier V, Nordmann P. 2011. Multiplex PCR for detection of acquired carbapenemase genes. Diagn Microbiol Infect Dis 70: 119–123. https://doi.org/10.1016/j.diagmicrobio.2010.12.002.

18. Lam MMC, Wyres KL, Duchene S, Wick RR, Judd LM, Gan YH, Hoh CH, Archuleta S, Molton JS, Kalimuddin S, Koh TH, Passet V, Brisse S, Holt KE. 2018. Population genomics of hypervirulent Klebsiella pneumoniae clonal-group 23 reveals early emergence and rapid global dissemination. Nat Commun 9:2703. https://doi.org/10.1038/s41467-018-05114-7.

19. Antimicrobial Resistance Collaborators. 2022. Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet 399: 629–655. https://doi.org/10.1016/S0140-6736(21)02724-0.

20. Thorpe T, Booton R, Kallonen T, Gibbon MJ, Couto N, Passet V, Lopez Fernandez JS, Rodrigues C, Matthews L, Mitchell S, Reeve E, David S, Merla C, Corbella M, Ferrari C, Comandatore F, Marone P, Brisse S, Sassera S, Corander J, Feil EJ. 2021. One Health or Three? Transmission modelling of Klebsiella isolates reveals ecological barriers to transmission between humans, animals and the environment. bioRxiv. https://doi.org/10.1101/2021.08.05.455249.

21. Marques C, Belas A, Aboim C, Cavaco-Silva P, Trigueiro G, Gama LT, Pomba C. 2019. Evidence of sharing of Klebsiella pneumoniae strains between healthy companion animals and cohabiting humans. J Clin Microbiol 57:e01537-18. https://doi.org/10.1128/JCM.01537-18.

22. Garcia-Fierro R, Drapeau A, Dazas M, Saras E, Rodrigues C, Brisse S, Madec J-Y, Haenni M. 2022. Comparative phylogenomics of ESBL-, AmpC- and carbapenemase-producing Klebsiella pneumoniae originating from companion animals and humans. J Antimicrob Chemother 77:1263–1271. https://doi.org/10.1093/jac/dkac041.

23. Breurec S, Guessennd N, Timinouni M, Le TAH, Cao V, Ngandjio A, Randrianirina F, Thiberge JM, Kinana A, Dufougeray A, Perrier-Gros-Claude JD, Boisier P, Garin B, Brisse S. 2013. Klebsiella pneumoniae resistant to third-generation cephalosporins in five African and two Vietnamese major towns: monoclonal population structure with two major international clonal groups, CG15 and CG258. Clin Microbiol Infect 19:349–355. https://doi.org/10.1111/j.1469-0691.2012.03805.x.

24. Peirano G, Chen L, Kreiswirth BN, Pitout JDD. 2020. Emerging antimicrobial-resistant high-risk Klebsiella pneumoniae clones ST307 and ST147. Antimicrob Agents Chemother 64:e01148-20. https://doi.org/10.1128/AAC.01148-20.

25. Hu Y, Anes J, Devineau S, Fanning S. 2021. Klebsiella pneumoniae: prevalence, reservoirs, antimicrobial resistance, pathogenicity, and infection. A hitherto unrecognized zoonotic bacterium. Foodborne Pathol Dis 18: 63–84. https://doi.org/10.1089/fpd.2020.2847.

26. Rodrigues C, Hauser K, Cahill N, Ligowska-Marzeta M, Centorotola G, Cornacchia A, Garcia Fierro R, Haenni M, Nielsen EM, Piveoteau P, Barbier E, Morris D, Pomilio F, Brisse S. 2022. High prevalence of Klebsiella pneumoniae in European food products: a multicentric study comparing culture and molecular detection methods. Microbiol Spectr 10:e02376-21. https://doi.org/10.1128/spectrom.02376-21.

27. Zadoks RN, Griffiths HM, Munoz MA, Ahlstrom C, Bennett GJ, Thomas E, Schukken YH. 2011. Sources of Klebsiella and Raoultella species on dairy farms: be careful where you walk. J Dairy Sci 94:1045–1051. https://doi.org/10.3168/jds.2010-3603.

28. Zhong XS, Li YZ, Ge J, Xiao G, Mo Y, Wen YQ, Liu JP, Xiong YQ, Qiu M, Huo ST, Cheng MJ, Chen Q. 2020. Comparisons of microbiological characteristics and antibiotic resistance of Klebsiella pneumoniae isolates from urban rodents, shrews, and healthy people. BMC Microbiol 20:12. https://doi.org/10.1186/s12866-020-1702-5.

29. Guyomard-Rabenirina S, Reynaud Y, Pot M, Albina E, Couvin D, Ducat C, Gruel G, Ferdinand S, Legreneur P, Le Hello S, Malpote E, Sadikalay S, Talarmin A, Breurec S. 2020. Antimicrobial resistance in wildlife in Guadeloupe (French West Indies): distribution of a single bla<sub>CTX-M-1</sub>/Incl1/ST3 plasmid among humans and wild animals. Front Microbiol 11:1524. https://doi.org/10.3389/fmicb.2020.01524.

30. Guyomard-Rabenirina S, Weill FX, Le Hello S, Bastian S, Berger F, Ferdinand S, Legreneur P, Loraux C, Malpote E, Muanza B, Richard V, Talarmin A, Breurec S. 2019. Reptiles in Guadeloupe (French West Indies) are a reservoir of major human Salmonella enterica serovars. PLoS One 14:e0220145. https://doi.org/10.1371/journal.pone.0220145.

31. Pot M, Reynaud Y, Couvin D, Ducat C, Ferdinand S, Gravey F, Gruel G, Guerin F, Malpote E, Breurec S, Talarmin A, Guyomard-Rabenirina S. 2021. Wide distribution and specific resistance pattern to third-generation cephalosporins of Enterobacter cloacae complex members in humans and in the environment in Guadeloupe (French West Indies). Front Microbiol 12:628058. https://doi.org/10.3389/fmicb.2021.628058.

32. David S, Reuter S, Harris SR, Glasner C, Feltwell T, Argimon S, Abudahab K, Goater R, Giani T, Errico G, Asbury M, Sjunnebo S, EuSCAPE Working
Group, ESGEM Study Group, Feil EJ, Rossolini GM, Aanensen DM, Grundmann H. 2019. Epidemic of carbapenem-resistant Klebsiella pneumoniae in Europe is driven by nosocomial spread. Nat Microbiol 4: 1919-1929. https://doi.org/10.1038/s41564-019-0492-8.

33. Buckner MMC, Saw HTH, Osagie RN, McNally A, Ricci V, Wand ME, Woodford N, Ivens A, Webber MA, Piddock LJV. 2018. Clinically relevant plasmid-host interactions indicate that transcriptional and not genomic modifications ameliorate fitness costs of Klebsiella pneumoniae carbapenemase-carrying plasmids. mBio 9:e02303-17. https://doi.org/10.1128/mBio.02303-17.

34. Rodrigues C, d'Humieres C, Papin G, Passet V, Ruppe E, Brisse S. 2020. Community-acquired infection caused by the uncommon hypervirulent Klebsiella pneumoniae ST66-K2 lineage. Microb Genom 6:mgen000419. https://doi.org/10.1099/mgen.0.000419.

35. Russo TA, Olson R, Macdonald U, Metzger D, Maltese LM, Drake EJ, Gulick AM. 2014. Aerobactin mediates virulence and accounts for increased siderophore production under iron-limiting conditions by hypervirulent (hypermucoviscous) Klebsiella pneumoniae. Infect Immun 82:2356-2367. https://doi.org/10.1128/IAI.01667-13.

36. Russo TA, Olson R, MacDonald U, Beanan J, Davidson BA. 2015. Aerobactin, but not yersiniabactin, salmochelin, or enterobactin, enables the growth/survival of hypervirulent (hypermucoviscous) Klebsiella pneumoniae ex vivo and in vivo. Infect Immun 83:3325-3333. https://doi.org/10.1128/IAI.00430-15.

37. Leangapichart T, Lunha K, Jiwakanon J, Angkititrakul S, Järhult JD, Magnusson U, Sunde M. 2021. 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 76:2012-2016. https://doi.org/10.1093/jac/dkab118.

38. Lam MMC, Wick RR, Wyres KL, Gorrie CL, Judd LM, Jenney AWJ, Brisse S, Holt KE. 2018. Genetic diversity, mobilisation and spread of the yersinia-bactin-encoding mobile element ICEKp in Klebsiella pneumoniae populations. Microb Genom 4:e000196. https://doi.org/10.1099/mgen.0.000196.

39. Gruel G, Sellin A, Riveiro H, Pot M, Breurec S, Guyomard-Rabenirina S, Talarmin A, Ferdinand S. 2021. Antimicrobial use and resistance in Escherichia coli from healthy food-producing animals in Guadeloupe. BMC Vet Res 17:116. https://doi.org/10.1186/s12917-021-02810-3.

40. Magiorakos AP, Srinivasan A, Carey RB, Carmeli Y, Falagas ME, Giske CG, Harbarth S, Hindler JF, Kahlmeter G, Olsson-Liljequist B, Paterson DL, Rice LB, Stelling J, Struelens MJ, Vatopoulos A, Weber JT, Monnet DL. 2012. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin Microbiol Infect 18:268-281. https://doi.org/10.1111/j.1469-0691.2011.03570.x.

41. Criscuolo A, Brisse S. 2014. AlienTrimmer removes adapter oligonucleotides with high sensitivity in short-insert paired-end reads. Commentary on Turner (2014) Assessment of insert sizes and adapter content in FASTQ data from NexteraXT libraries. Front Genet 5:130. https://doi.org/10.3389/fgene.2014.00130.

42. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA. 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19:455-477. https://doi.org/10.1089/cmb.2012.0021.

43. Gurevich A, Saveliev V, Vyahhi N, Tesler G. 2013. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29:1072-1075. https://doi.org/10.1093/bioinformatics/btt086.

44. Lam MMC, Wick RR, Watts SC, Cerdeira LT, Wyres KL, Holt KE. 2021. A genomic surveillance framework and genotyping tool for Klebsiella pneumoniae and its related species complex. Nat Commun 12:4188. https://doi.org/10.1038/s41467-021-24448-3.

45. Carattoli A, Zankari E, Garcia-Fernandez A, Voldby Larsen M, Lund O, Villa L, Moller Aarestrup F, Hasman H. 2014. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother 58:3895-3903. https://doi.org/10.1128/AAC.02412-14.

46. Krawczyk PS, Lipinski L, Dziembowski A. 2018. PlasFlow: predicting plasmid sequences in metagenomic data using genome signatures. Nucleic Acids Res 46:e35. https://doi.org/10.1093/nar/gkx1321.

47. Robertson J, Nash JHE. 2018. MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies. Microb Genom 4:e000206. https://doi.org/10.1099/mgen.0.000206.

48. Huynh BT, Passet V, Rakotondrasoa A, Diallo T, Kerleguer A, Hennart M, Lauzanne A, Herindrainy P, Seck A, Bercion R, Borand L, Pardos de la Gandara M, Delarocque-Astagneau E, Guillemot D, Vray M, Garin B, Collard JM, Rodrigues C, Brisse S. 2020. Klebsiella pneumoniae carriage in low-income countries: antimicrobial resistance, genomic diversity and risk factors. Gut Microbes 11:1287-1299. https://doi.org/10.1080/19490976.2020.1748257.

49. Alonge M, Soyk S, Ramakrishnan S, Wang X, Goodwin S, Sedlazeck FJ, Lippman ZB, Schatz MC. 2019. RaGOO: fast and accurate reference-guided scaffolding of draft genomes. Genome Biol 20:224. https://doi.org/10.1186/s13059-019-1829-6.

50. Seemann T. 2014. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30:2068-2069. https://doi.org/10.1093/bioinformatics/btu153.

51. Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MT, Fookes M, Falush D, Keane JA, Parkhill J. 2015. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics 31:3691-3693. https://doi.org/10.1093/bioinformatics/btv421.

52. Didelot X, Wilson DJ. 2015. ClonalFrameML: efficient inference of recombination in whole bacterial genomes. PLoS Comput Biol 11:e1004041. https://doi.org/10.1371/journal.pcbi.1004041.

53. Stamatakis A. 2014. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30:1312-1313. https://doi.org/10.1093/bioinformatics/btu033.

54. Letunic I, Bork P. 2021. Interactive Tree of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res 49:W293-W296. https://doi.org/10.1093/nar/gkab301.

55. Argimon S, David S, Underwood A, Abrudan M, Wheeler NE, Kekre M, Abudahab K, Yeats CA, Goater R, Taylor B, Harste H, Muddyman D, Feil EJ, Brisse S, Holt K, Donado-Godoy P, Ravikumar KL, Okeke IN, Carlos C, Aanensen DM, NIHR Global Health Research Unit on Genomic Surveillance of Antimicrobial Resistance. 2021. Rapid genomic characterization and global surveillance of Klebsiella using Pathogenwatch. Clin Infect Dis 73:S325-S335. https://doi.org/10.1093/cid/ciab784.

56. Brynildsrud O, Bohlin J, Scheffer L, Eldholm V. 2016. Rapid scoring of genes in microbial pan-genome-wide association studies with Scoary. Genome Biol 17:238. https://doi.org/10.1186/s13059-016-1108-8.