# Phenotypic and genotypic within-host diversity of Pseudomonas aeruginosa urinary isolates

Agnès Cottalorda, Sandrine Dahyot, Anaïs Soares, Kevin Alexandre, Isabelle Zorgniotti, Manuel Etienne, Estelle Jumas-Bilak, Martine Pestel-Caron

► To cite this version:

Agnès Cottalorda, Sandrine Dahyot, Anaïs Soares, Kevin Alexandre, Isabelle Zorgniotti, et al. Phenotypic and genotypic within-host diversity of Pseudomonas aeruginosa urinary isolates. Scientific Reports, 2022, 12, 10.1038/s41598-022-09234-5. insu-03665281

HAL Id: insu-03665281

https://insu.hal.science/insu-03665281

Submitted on 11 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished 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_1.png" /></figure>

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

Check for updates

OPEN

# Phenotypic and genotypic
within-host diversity
of *Pseudomonas aeruginosa* urinary
isolates

Agnès Cottalorda<fnref n="1" />✉, Sandrine Dahyot², Anaïs Soares², Kevin Alexandre³,
Isabelle Zorgniotti⁴, Manuel Etienne³, Estelle Jumas-Bilak⁴ & Martine Pestel-Caron²

This study aimed to assess phenotypic and molecular inter-patient and within-host diversity of *Pseudomonas aeruginosa* isolates responsible for urinary tract infection (UTI) or asymptomatic bacteriuria (AB). Clinical data of 120 consecutive *P. aeruginosa* UTI (n = 40) and AB (n = 80) were prospectively analyzed. Up to five *P. aeruginosa* isolates per sample were collected. Antimicrobial susceptibility testing (AST) was determined for all isolates (n = 591); a subset of 358 was characterized by multilocus sequence typing. 444 isolates (75%) were non-multidrug resistant (MDR), 113 (19%) were MDR, and 34 (6%) were extensively drug resistant. A genetically highly diverse population was observed (64 sequence types [STs]), without strict correlation between genotypes and clinical settings. 35 patients (28%; 12 UTIs and 23 ABs) presented distinct antimicrobial resistance (AMR) profiles within a given urine sample, significantly associated with previous carbapenem and fluroquinolones exposure; five of them also exhibited polyclonal UTI or AB (with isolates belonging to two STs). *P. aeruginosa* urinary isolates of these 120 patients were highly diverse, in terms of AMR as well as genetic background. Both within-host AMR and molecular diversity can complicate AST, treatment and control of *P. aeruginosa* UTI.

*Pseudomonas aeruginosa* is an opportunistic pathogen responsible for 7–10% of healthcare-associated urinary tract infections (UTI)<fnref n="1" />⁻³. *P. aeruginosa* UTI commonly affects patients with underlying conditions, such as urinary tract abnormalities²,⁴, and are particularly promoted by urinary indwelling catheters².

Because of its ability to develop resistance to multiple classes of antibiotics, multidrug resistance (MDR) is frequently reported in *P. aeruginosa* urinary strains⁴,⁵. This phenomenon, associated with *P. aeruginosa* ability to form biofilm, leads to infections that are difficult to treat⁶.

To explore the diffusion of resistant strains and clonal relatedness of *P. aeruginosa* isolates, a large variety of genotyping methods has been developed such as Multilocus sequence typing (MLST). MLST has successfully identified epidemic high-risk clones resistant to antibiotics and responsible for healthcare-associated infections⁷. On the other hand, within-host diversity of *P. aeruginosa* isolates has been identified in cystic fibrosis (CF) infections, which may complicate patient's care⁸⁻¹⁰. In contrast, few studies have explored the genetic diversity of *P. aeruginosa* urinary strains; they often included a limited number of strains and mainly focused on multiresistant ones¹¹,¹². Thus, little is known about both molecular epidemiology and within-sample diversity of *P. aeruginosa* urinary isolates.

In this context, the aims of this study were (1) to investigate antimicrobial resistance (AMR) and genetic background of *P. aeruginosa* isolates collected from UTI or asymptomatic bacteriuria (AB), and (2) to explore within-host diversity of these isolates.

<footnote n="1">¹GRAM 2.0, Normandie Univ, UNIROUEN, UNICAEN, 76000 Rouen, France. ²GRAM 2.0, CHU Rouen, Department of Microbiology, Normandie Univ, UNIROUEN, UNICAEN, 76000 Rouen, France. ³GRAM 2.0, CHU Rouen, Department of Infectious Diseases, Normandie Univ, UNIROUEN, UNICAEN, 76000 Rouen, France. ⁴Team Pathogènes Hydriques Santé Environnement, UMR 5569 HydroSciences Montpellier, University of Montpellier, Montpellier, France. ⁵email: agnes.cottalorda@chu-rouen.fr</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Overall resistance rates of the 591 *P. aeruginosa* urinary isolates to antimicrobial groups. Overall resistance rates are presented for each antimicrobial group in bold, above the bar chart. Percentage of non-MDR (MultiDrug Resistant), MDR, and XDR (eXtensively Drug Resistant) profiles are indicated for each antimicrobial group inside the bar chart.</figcaption></figure>

## Results

**Demographic and clinical characteristics of patients.** The demographic and clinical characteristics of the 120 consecutive patients suffering from *P. aeruginosa* bacteriuria considered as UTI or AB are presented in Table S1. The mean age of the patients was 64 years, ranging from 1 month to 101 years. 64% (77/120) were male. Mean Charlson index was 5.4. Patients were mostly hospitalized in medical (46/120, 38%) or in surgery (32/120, 27%) wards. Within the 6 months before inclusion, 86% were hospitalized, 68% underwent urinary tract manipulation, and 78% received at least one antibiotic treatment (Table S1). The antibiotic classes previously received by patients were mostly penicillins (72/120, 60%), and cephalosporins (43/120, 36%) (Table S2).

Compared to AB, demographic or clinical characteristics significantly associated with UTI were male sex (78% versus 58%), bacteraemia (30% versus 3%), and uncomplicated diabetes (32% versus 13%) (*P* < 0.05) (Table S1).

**Phenotypic and molecular characteristics of isolates.** Antimicrobial susceptibility testing (AST). Among the 591 isolates tested (up to five per urine sample), resistance rates ranged from 13 to 34% depending on the antimicrobial groups (Fig. 1). Low rates (≤ 20%) were observed for cephalosporins, monobactam, and aminoglycosides while rates were higher (≥ 25%) for fluoroquinolones (34%), carbapenems (27%) and penicillins (26%). 444 isolates (75%) were non-MDR (including 318 wild-type isolates [54%]), 113 (19%) were MDR and 34 (6%) were extensively drug resistant (XDR) (Fig. 1).

MDR and/or XDR profiles were significantly more observed for patients with previous urine culture positive for *P. aeruginosa* or antibiotic treatment (especially cephalosporins and fluoroquinolones) (*P* < 0.05) (Table S2). In addition, patients with resistant isolates received a higher mean number of different antimicrobial groups in the past 6 months (i.e. 1.5, 2.6, 3.3 for patients with non-MDR, MDR, and XDR isolates, respectively).

*Genotyping by MLST.* 64 STs were identified for the 358 isolates tested (two to three per urine sample). 46 STs were singletons, as identified in only one patient (Table S3). The most predominant STs were ST395 and ST308 (11 patients each), followed by ST253 (7 patients), ST309 and ST235 (6 patients each).

29 STs were identified for the 120 UTI isolates, and 46 STs for the 238 AB isolates. No strict correlation was found between STs and clinical setting, as 11 of the 18 STs identified in at least two patients were observed in both UTI and AB contexts (Fig. 2A). Despite this, ST309 (6 patients) and ST298 (4 patients) tended to be associated with AB and UTI, respectively. Furthermore, 4 STs (ST308, ST235, ST111 and ST175) were mainly associated with MDR and XDR profiles (Fig. 2B).

**Within-host diversity.** Within-host diversity was assessed by comparing AST profiles and STs between isolates from a given sample. Among the 120 patients, 85 (71%) presented no within-host AMR diversity (71 with non-MDR, 10 with MDR, and 4 with XDR isolates), whereas 35 presented within-host AMR diversity (Fig. 3A). Minor discrepancies were observed for 21 of them. Interestingly, 14 patients (12%) harbored isolates with two
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Minimum spanning tree of the 358 *P. aeruginosa* urinary isolates typed by multilocus sequence typing (MLST). The sequences were concatenated and analyzed with BioNumerics. Clustering of MLST profiles was done using a categorical coefficient. The colors used are based (A) on clinical contexts (AB: Asymptomatic Bacteriuria; UTI: Urinary Tract Infection) and (B) on antimicrobial resistance profiles (MDR: MultiDrug Resistant; XDR: eXtensively Drug Resistant). Each circle represents a sequence type (ST) and its size is proportional to the number of isolates. Length of the lines represent the genetic distance between isolates. Thick, short lines connecting two types denote types differing in a single locus; thin, longer lines connect double-locus variants; and dashed lines indicate the most likely connection between two types differing in more than two loci. STs identified in at least two patients are annotated, and the three most prevalent worldwide epidemic high-risk clones<sup>7</sup> (B) are in bold.</figcaption></figure>

distinct AMR profiles, corresponding to 13 major discrepancies (9 non-MDR/MDR, and 4 MDR/XDR) and one very major (non-MDR/XDR) (Fig. 3A). AMR diversity concerned all antibiotic classes (Fig. 3B). Comparison of clinical characteristics between patients with or without AMR diversity revealed that previous exposure to carbapenems or fluoroquinolones was significantly more frequently observed in patients with AMR diversity (*P* < 0.05) (Table 1).

Remarkably, urine sample of five patients (4%; 2 with UTI and 3 with AB) harbored isolates with two distinct STs. These STs differed for 4 to 7 alleles for 4 patients and for a single *locus* (*mutL*) for one patient (Table 2). Within-host genetic diversity was always associated with within-host AMR diversity, with either minor (two patients), major (two patients), or very major discrepancies (one patient). With the limit of the small number of patients, this genetic diversity seemed to occur preferentially in patients with prior hospitalization (all the five patients) and previous antibiotic treatment (four patients) (Table 1).

Discussion

*P. aeruginosa* is frequently responsible for healthcare-associated UTI, but little is known about the genetic back-
ground of the urinary strains or within-host diversity. Thus, we investigated phenotypic and molecular diversity
of *P. aeruginosa* in the setting of UTI or AB, by exploring multiple isolates per urine sample.
<figure><img src="image_5.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Within-host AMR diversity. (A) Number and percent of patients with within-host antimicrobial resistance (AMR) diversity, including minor, major, or very major discrepancies. (B) Percent of patients with within-host AMR diversity according to antibiotic classes. Overall diversity rates are presented for each antimicrobial group in bold, above the bar chart. Percentage of AMR diversity due to distinct categorization (susceptible, intermediate or resistant) or relative standard deviation ≥ 20% of the inhibition zone diameters are indicated for each antimicrobial group inside the bar chart.</figcaption></figure>

Patients enrolled in this study presented characteristics previously described as risk factors for *P. aeruginosa* UTI (i.e. male sex, urological disorder, urinary tract catheter and previous UTI/antibiotic therapy/hospitalization)¹,²,⁴. The rate of MDR/XDR (25%) isolates was close to that reported worldwide (15–30%)⁷ and to that observed in studies based on *P. aeruginosa* urinary isolates (10.2–41.2%)⁴,¹³.

Genetic diversity of our isolates was assessed by MLST, commonly used to describe the global epidemiology of *P. aeruginosa*⁷,¹⁴,¹⁵. A genetically highly diverse population (64 STs) was observed, as previously described for other infections¹⁶, with no strict correlation between genotype and clinical setting (AB or UTI). Although most of the STs were identified in only one patient, 13 of the 18 spread clones reviewed by Oliver et al.⁷ were identified, including the two predominant clones in our population (ST395, and ST308). Some of them were mainly associated with multiresistance, including the three most predominant worldwide epidemic high-risk clones (ST111, ST175, ST235)⁷. Of note, ST234 and ST235, previously described for *P. aeruginosa* urinary isolates¹¹,¹², were identified for some of our isolates.

Within-host diversity was then assessed by comparing phenotypes and genotypes of isolates collected from a given sample.

It first revealed that almost one-third of our patients exhibited within-sample AMR diversity that could be associated with phenotypic plasticity. This observation is a major concern because it can lead to treatment failure. This result supports the recommendation to pool several distinguishable colonies of *P. aeruginosa* when performing AST¹⁷, or even to carry out a direct AST by disk diffusion from urine specimen¹⁸, in order to identify AMR diversity that could impact the choice of antibiotic treatment. This latter method provides a rapid solution for determining AST and facilitates identification of resistant subpopulations¹⁸. The AMR diversity was previously observed by Mowat et al.⁹ who analyzed sets of 40 *P. aeruginosa* cystic fibrosis (CF) isolates per sputum sample. Despite the lower number of isolates per sample explored here, the high AMR diversity rate suggests that the urinary tract could be a complex environment that favor adaptation, even if the selective pressure in the urinary tract is different than during chronic long-term CF infections.
<table><caption>**Table 1.** Clinical characteristics of patients with phenotypic and molecular within-host diversity. MDR multidrug resistant, XDR extensively drug resistant, AMR antimicrobial resistance, min minimum, max maximum, SD standard deviation. *Charlson index was measured only for adult cases (i.e. 114 patients, 37 with UTI, and 77 with AB). **P value < 0.05.</caption><thead><tr><th>Characteristics</th><th>No diversity (n=85)</th><th>AMR diversity (n=35)</th><th>AMR diversity associated with genetic diversity (n=5)</th></tr></thead><tbody><tr><td><strong>Demographics no. (%)</strong></td><td></td><td></td><td></td></tr><tr><td>Age (years) [min; max]</td><td>64 [0; 97]</td><td>64 [1; 101]</td><td>54 [36; 83]</td></tr><tr><td>Male gender</td><td>55 (65%)</td><td>22 (63%)</td><td>4</td></tr><tr><td><strong>Clinical features no. (%)</strong></td><td></td><td></td><td></td></tr><tr><td>Urinary tract infection</td><td>28 (33%)</td><td>12 (34%)</td><td>2</td></tr><tr><td>Urinary tract catheter</td><td>42 (49%)</td><td>18 (51%)</td><td>3</td></tr><tr><td>Bacteraemia caused by <em>P. aeruginosa</em></td><td>8/43 (19%)</td><td>1/17 (6%)</td><td>0/2</td></tr><tr><td><strong>Comorbid diseases no. (%)</strong></td><td></td><td></td><td></td></tr><tr><td>Urinary comorbidity</td><td>40 (47%)</td><td>21 (60%)</td><td>3</td></tr><tr><td>Charlson index (avg ± SD)*</td><td>5.4 ± 3.0</td><td>5.4 ± 2.5</td><td>3.2 ± 2.2</td></tr><tr><td>Other comorbidity</td><td>9 (11%)</td><td>3 (9%)</td><td>0</td></tr><tr><td><strong>Clinical histories within the previous 6 months no. (%)</strong></td><td></td><td></td><td></td></tr><tr><td>Previous urine culture positive to <em>P. aeruginosa</em></td><td>8 (9%)</td><td>7 (20%)</td><td>2</td></tr><tr><td>Previous antibiotic treatment</td><td>65 (77%)</td><td>28 (80%)</td><td>4</td></tr><tr><td>Penicillins</td><td>50 (59%)</td><td>22 (63%)</td><td>1</td></tr><tr><td>Cephalosporins</td><td>31 (36%)</td><td>12 (34%)</td><td>1</td></tr><tr><td>Carbapenems</td><td>5 (6%)</td><td>7 (20%)**</td><td>2</td></tr><tr><td>Aminoglycosides</td><td>17 (20%)</td><td>8 (23%)</td><td>0</td></tr><tr><td>Quinolones</td><td>11 (13%)</td><td>11 (31%)**</td><td>1</td></tr><tr><td>Other antimicrobials</td><td>35 (41%)</td><td>17 (49%)</td><td>2</td></tr><tr><td>Previous hospitalization</td><td>73 (86%)</td><td>30 (86%)</td><td>5</td></tr><tr><td>Previous urinary tract manipulation</td><td>59 (69%)</td><td>23 (66%)</td><td>3</td></tr><tr><td>Including previous urological surgery</td><td>17 (20%)</td><td>4 (11%)</td><td>1</td></tr></tbody></table>

<table><caption>**Table 2.** Allelic profiles of *P. aeruginosa* isolates from patients with polyclonal urinary tract infection or asymptomatic bacteriuria. The seven alleles (*acsA*, *aroE*, *guaA*, *mutL*, *nuoD*, *ppsA* and *trpE*) that defined sequence type (ST) are indicated here. Alleles that differ between the 2 STs of the isolates of a given patient are in bold.</caption><thead><tr><th>Patients</th><th>Clinical context</th><th>STs</th><th>Nb of isolates</th><th>acsA</th><th>aroE</th><th>guaA</th><th>mutL</th><th>nuoD</th><th>ppsA</th><th>trpE</th><th>Nb of alleles differentiating STs</th></tr></thead><tbody><tr><td rowspan="2">A</td><td rowspan="2">UTI</td><td>253</td><td>2</td><td>4</td><td>4</td><td>16</td><td>12</td><td>1</td><td>6</td><td>3</td><td rowspan="2">4</td></tr><tr><td>207</td><td>1</td><td><strong>47</strong></td><td>4</td><td><strong>5</strong></td><td><strong>33</strong></td><td>1</td><td>6</td><td><strong>40</strong></td></tr><tr><td rowspan="2">B</td><td rowspan="2">AB</td><td>2406</td><td>2</td><td>29</td><td>8</td><td>36</td><td>67</td><td>1</td><td>6</td><td>3</td><td rowspan="2">1</td></tr><tr><td>3232</td><td>1</td><td>29</td><td>8</td><td>36</td><td><strong>228</strong></td><td>1</td><td>6</td><td>3</td></tr><tr><td rowspan="2">C</td><td rowspan="2">AB</td><td>308</td><td>2</td><td>13</td><td>4</td><td>5</td><td>5</td><td>12</td><td>7</td><td>15</td><td rowspan="2">5</td></tr><tr><td>207</td><td>1</td><td><strong>47</strong></td><td>4</td><td>5</td><td><strong>33</strong></td><td><strong>1</strong></td><td><strong>6</strong></td><td><strong>40</strong></td></tr><tr><td rowspan="2">D</td><td rowspan="2">UTI</td><td>683</td><td>2</td><td>6</td><td>5</td><td>11</td><td>3</td><td>4</td><td>4</td><td>1</td><td rowspan="2">7</td></tr><tr><td>3233</td><td>1</td><td><strong>87</strong></td><td><strong>34</strong></td><td><strong>114</strong></td><td><strong>37</strong></td><td><strong>86</strong></td><td><strong>91</strong></td><td><strong>170</strong></td></tr><tr><td rowspan="2">E</td><td rowspan="2">AB</td><td>483</td><td>2</td><td>116</td><td>5</td><td>6</td><td>5</td><td>3</td><td>12</td><td>68</td><td rowspan="2">7</td></tr><tr><td>446</td><td>1</td><td><strong>18</strong></td><td><strong>4</strong></td><td><strong>5</strong></td><td><strong>3</strong></td><td><strong>1</strong></td><td><strong>17</strong></td><td><strong>13</strong></td></tr></tbody></table>

Then, our study is the first to demonstrate that *P. aeruginosa* UTI or AB can be polyclonal, as two distinct STs were identified within a urinary sample for five patients (4%). For four of them, at least four alleles distinguished the two clones. Analysis of previous urine samples could have been of interest to establish whether these clones already coexisted, especially for the two patients who had a history of *P. aeruginosa* bacteriuria. Moreover, these five patients were previously hospitalized, and four of them received antimicrobial treatments that could have selected resistant strains. This is supported by the fact that genetic diversity was always correlated with AMR diversity.

Nevertheless, the number of patients harboring polyclonal bacteriuria was smaller than that described by Waine et al.¹⁰ in the CF context, where 2 to 3 STs were identified for half of the 60 patients included. However, that study was based on a much larger number of isolates (from 27 to 46 per sputum). Further studies with many
more isolates per urine sample are therefore needed to evaluate within-sample diversity to a greater extent, and to identify the impact of polyclonal UTI on patient outcome.

On the other hand, the resolution of MLST is limited due to the use of a small number of genes. It would be interesting to further characterize isolates with the same ST by whole genome sequencing given the high discriminatory power of this technology¹⁹,²⁰. In this way, genomic comparisons of *P. aeruginosa* strains isolated from CF²¹ and non-CF patients²² have showed that a parental strain could diversify into distinct sublineages. This phenomenon could be identified for patients enrolled here (13% of them exhibiting a previous urine culture positive to *P. aeruginosa*) and would be interesting to analyze in order to better understand mechanisms of *P. aeruginosa* adaptation to the urinary tract. Furthermore, metagenomic studies would be interesting to describe interactions of *P. aeruginosa* with the urinary microbiome²³.

Our study presents the limitations of a monocentric study. Its results may not be generalized to other geographic areas. Nevertheless, isolates were collected from several clinical wards over a period of 2 years, and the high proportion of singleton STs identified here (72%) was consistent with the *P. aeruginosa* population structure²⁴.

In conclusion, this study provides strong evidence of within-host diversity of *P. aeruginosa* urinary isolates, whether phenotypic or genotypic. These findings can complicate diagnosis, treatment and control of *P. aeruginosa* UTI.

## Methods

**Pseudomonas aeruginosa strain collection and patient characteristics.** From June 2016 to August 2018, we prospectively collected isolates for all monomicrobial urine cultures positive for *P. aeruginosa* (regardless of the level of bacteriuria and leucocyturia) in the Rouen University Hospital. Clinical data were prospectively collected from the hospital's computerized medical records. The diagnosis of UTI or AB was assigned according to the diagnosis retained by the physician in charge of the patient, and confirmed by a committee made up of microbiologists and infectious disease physicians. A total of 120 patients were consecutively included, 40 with UTI and 80 with AB (Table S4).

For each urine culture, two (in case of low-level bacteriuria, i.e. $2.10^2$ CFU/mL) to five single colonies representative of all morphotypes observed were selected. Bacterial isolates were identified by matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry (Bruker Daltonik GmbH, Bremen, Germany) and stored at -80 °C for further studies.

**Antibiotic susceptibility.** We assessed the activity of 14 antibiotics for the 591 isolates by disk diffusion method on Mueller-Hinton agar (Bio-Rad, Marnes-la-Coquette, France) according to the European Committee on Antimicrobial Susceptibility Testing recommendations (https://www.eucast.org/), as previously described²⁵. Antimicrobial disks (Bio-Rad) belonged to distinct antimicrobial groups: penicillins (ticarcillin, ticarcillin-clavulanate, piperacillin, piperacillin-tazobactam), cephalosporins (ceftazidime, cefepime), monobactam (aztreonam), carbapenems (imipenem, meropenem), aminoglycosides (gentamicin, tobramycin, amikacin), and fluoroquinolones (ciprofloxacin, lévofloxacine). For each experiment, a quality control strain (*P. aeruginosa* strain ATCC 27853) was performed. A relative standard deviation of the inhibition zone diameter between intra-sample isolates was calculated for the 15 antibiotics tested.

The proportion of MDR (not susceptible to at least three antimicrobial groups) and XDR (not susceptible to the six antimicrobial groups tested) isolates was determined according to consensus recommendations²⁶. The other isolates—resistant to none, one or two groups—were categorized as non-MDR.

Within-host AMR diversity was assessed between intra-sample isolates by comparing AMR profile (non-MDR, MDR, or XDR) and/or according to the relative standard deviation of the inhibition zone diameters within a given profile. Diversity was defined as follows: a major discrepancy as a change to the next category (from non-MDR to MDR or MDR to XDR), while a change from non-MDR to XDR was considered as a very major discrepancy. A minor discrepancy was defined for isolates exhibiting the same AMR profile but with a distinct categorization (susceptible, intermediate or resistant) to one or more antibiotic groups or with a relative standard deviation ≥ 20% for at least one of the 15 antibiotics tested.

To identify clinical risk factors for multidrug resistance (non-MDR versus MDR/XDR), patients exhibiting isolates with different AMR profiles were classified according to the most resistant one (i.e. patients with non-MDR and MDR isolates were classified as MDR patients).

**MLST.** MLST was performed for two (in case of low-level bacteriuria) to three isolates per urine culture selected to represent the phenotypic diversity (in terms of morphotype and AMR profile).

MLST was carried out as previously described²⁷, except for the use of newly designed primers (Table S5). Briefly, bacterial DNA was extracted by InstaGene Matrix kit (Bio-Rad) according to the manufacturer's recommendations. PCR products were purified and sequenced by GENEWIZ Europe (Leipzig, Germany). Sequencing data were aligned with BioEdit software (http://www.mbio.ncsu.edu/BioEdit/bioedit.html). Allelic profiles and corresponding sequence types (STs) were assigned using the international PubMLST database (https://pubmlst.org/paeruginosa/). New alleles and STs were submitted to this database. Clinical and microbial data were uploaded into BioNumerics software (Version 7.6, Applied Maths, Belgium) in order to construct minimum spanning trees based on concatenated sequences.

**Statistical analyzes.** The Pearson's Chi-squared test was used to compare categorical data while ANOVA test was used to analyze continuous data. All tests were two-tailed, and significance was considered at *P* value < 0.05. Statistical analyzes were performed using R software (v.3.5.1).
**Statements on study approvals.** We confirm that all methods and protocols used in this study were carried out in accordance with relevant guidelines and regulations. Clinical data and isolates were anonymised so that they were irretrievably unrelated to an identifiable subject. According to French regulation on observational database analyzes, no specific informed consent was required from patients. The study was approved and registered by the Clinical Research and Innovation Delegation of the Rouen University Hospital under number 2018/413/OB. Ethical approval for this study was obtained from the Ethics Committee for Research of the Rouen University Hospital, Rouen, France (No. E2021-77).

Received: 29 September 2021; Accepted: 25 February 2022

Published online: 30 March 2022

**References**

1. Zhang, X., Niu, S. & Zhang, L. Antimicrobial susceptibilities and clinical characterization of pseudomonas aeruginosa isolates from urinary tract infections. *Urol. Int.* **93**, 464–469 (2014).

2. Djordjevic, Z., Folic, M. M., Zivic, Z., Markovic, V. & Jankovic, S. M. Nosocomial urinary tract infections caused by *Pseudomonas aeruginosa* and Acinetobacter species: Sensitivity to antibiotics and risk factors. *Am. J. Infect. Control* **41**, 1182–1187 (2013).

3. Lamas Ferreiro, J. L. et al. *Pseudomonas aeruginosa* urinary tract infections in hospitalized patients: Mortality and prognostic factors. *PLoS ONE* **12**, e0178178 (2017).

4. Gomila, A. et al. Risk factors and prognosis of complicated urinary tract infections caused by *Pseudomonas aeruginosa* in hospitalized patients: A retrospective multicenter cohort study. *Infect. Drug Resist.* **11**, 2571–2581 (2018).

5. Jiménez-Guerra, G. et al. Urinary tract infection by Acinetobacter baumannii and Pseudomonas aeruginosa: Evolution of antimicrobial resistance and therapeutic alternatives. *J. Med. Microbiol.* **67**, 790–797 (2018).

6. Morales, E. et al. Hospital costs of nosocomial multi-drug resistant *Pseudomonas aeruginosa* acquisition. BMC Health Serv. Res. **12**, 122 (2012).

7. Oliver, A., Mulet, X., López-Causapé, C. & Juan, C. The increasing threat of *Pseudomonas aeruginosa* high-risk clones. *Drug Resist. Updat. Rev. Comment. Antimicrob. Anticancer Chemother.* **21–22**, 41–59 (2015).

8. Markussen, T. et al. Environmental heterogeneity drives within-host diversification and evolution of *Pseudomonas aeruginosa*. *MBio* **5**, e01592-14 (2014).

9. Mowat, E. et al. *Pseudomonas aeruginosa* population diversity and turnover in cystic fibrosis chronic infections. *Am. J. Respir. Crit. Care Med.* **183**, 1674–1679 (2011).

10. Waine, D. J., Honeybourne, D., Smith, E. G., Whitehouse, J. L. & Dowson, C. G. Cross-sectional and longitudinal multilocus sequence typing of *Pseudomonas aeruginosa* in cystic fibrosis sputum samples. *J. Clin. Microbiol.* **47**, 3444–3448 (2009).

11. Pobiega, M. et al. Molecular characterization of carbapenem-resistant *Pseudomonas aeruginosa* strains isolated from patients with urinary tract infections in Southern Poland. *Diagn. Microbiol. Infect. Dis.* **83**, 295–297 (2015).

12. Osawa, K. et al. Molecular characteristics of carbapenem-resistant *Pseudomonas aeruginosa* isolated from urine in Hyogo, Japan. *Int. J. Urol.* **26**, 127–133 (2019).

13. Gajdács, M., Burián, K. & Terhes, G. Resistance levels and epidemiology of non-fermenting gram-negative bacteria in urinary tract infections of inpatients and outpatients (RENFUTI): A 10-year epidemiological snapshot. *Antibiot. Basel Switz.* **8**, 143 (2019).

14. Park, Y. & Koo, S. H. Epidemiology, molecular characteristics, and virulence factors of carbapenem-resistant *Pseudomonas aeruginosa* isolated from patients with urinary tract infections. *Infect. Drug Resist.* **15**, 141–151 (2022).

15. Pérez, A. et al. High incidence of MDR and XDR Pseudomonas aeruginosa isolates obtained from patients with ventilator-associated pneumonia in Greece, Italy and Spain as part of the MagicBullet clinical trial. *J. Antimicrob. Chemother.* **74**, 1244–1252 (2019).

16. Parkins, M. D., Somayaji, R. & Waters, V. J. Epidemiology, biology, and impact of clonal *Pseudomonas aeruginosa* infections in cystic fibrosis. *Clin. Microbiol. Rev.* **31**, e00019-18 (2018).

17. Winstanley, C., O'Brien, S. & Brockhurst, M. A. *Pseudomonas aeruginosa* evolutionary adaptation and diversification in cystic fibrosis chronic lung infections. *Trends Microbiol.* **24**, 327–337 (2016).

18. Coorevits, L., Boelens, J. & Claeys, G. Direct susceptibility testing by disk diffusion on clinical samples: A rapid and accurate tool for antibiotic stewardship. *Eur. J. Clin. Microbiol. Infect. Dis. Off. Publ. Eur. Soc. Clin. Microbiol.* **34**, 1207–1212 (2015).

19. Sullivan, C. B., Diggle, M. A. & Clarke, S. C. Multilocus sequence typing: Data analysis in clinical microbiology and public health. *Mol. Biotechnol.* **29**, 245–254 (2005).

20. Chen, J.-W., Lau, Y. Y., Krishnan, T., Chan, K.-G. & Chang, C.-Y. Recent advances in molecular diagnosis of *Pseudomonas aeruginosa* infection by state-of-the-art genotyping techniques. *Front. Microbiol.* **9**, 1104 (2018).

21. Bianconi, I. et al. Persistence and microevolution of *Pseudomonas aeruginosa* in the cystic fibrosis lung: A single-patient longitudinal genomic study. *Front. Microbiol.* **9**, 3242 (2018).

22. Hilliam, Y. et al. *Pseudomonas aeruginosa* adaptation and diversification in the non-cystic fibrosis bronchiectasis lung. *Eur. Respir. J.* **49**, 1602108 (2017).

23. Lewis, D. A. et al. The human urinary microbiome; bacterial DNA in voided urine of asymptomatic adults. *Front. Cell. Infect. Microbiol.* **3**, 41 (2013).

24. Horcajada, J. P. et al. Epidemiology and treatment of multidrug-resistant and extensively drug-resistant *Pseudomonas aeruginosa* infections. *Clin. Microbiol. Rev.* **32**, e00031-e119 (2019).

25. Cottalorda, A. et al. Within-host microevolution of *Pseudomonas aeruginosa* urinary isolates: A seven-patient longitudinal genomic and phenotypic study. *Front. Microbiol.* **11**, 611246 (2020).

26. Magiorakos, A.-P. et al. 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 (2012).

27. Curran, B., Jonas, D., Grundmann, H., Pitt, T. & Dowson, C. G. Development of a multilocus sequence typing scheme for the opportunistic pathogen *Pseudomonas aeruginosa*. *J. Clin. Microbiol.* **42**, 5644–5649 (2004).

**Acknowledgements**

This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors. The authors are greatly thankful to Mathilde Réveillon, François Croustillères, and Delphine Feyeux for technical assistance; Fabien Aujoulat for his help with MLST experiments and Maxime Grand for his help with statistical analyzes.
**Author contributions**

Author contributions included conception and design (A.C., E.J.B., M.P.C., S.D.), collection of strains and clinical data (A.C., A.S., K.A., M.E., M.P.C., S.D.), experiments (A.C., I.Z.), analysis and interpretation of data (A.C., E.J.B., K.A., M.E., M.P.C., S.D.), and manuscript preparation (A.C., A.S., K.A., M.P.C., S.D.). All authors amended and approved the final version of the manuscript.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-09234-5.

**Correspondence** and requests for materials should be addressed to A.C.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022