# Within-Host Microevolution of Pseudomonas aeruginosa Urinary Isolates: A Seven-Patient Longitudinal Genomic and Phenotypic Study

Agnès Cottalorda, Marie Leoz, Sandrine Dahyot, François Gravey, Maxime Grand, Thomas Froidure, Fabien Aujoulat, Simon Le Hello, Estelle Jumas-Bilak, Martine Pestel-Caron

► To cite this version:

Agnès Cottalorda, Marie Leoz, Sandrine Dahyot, François Gravey, Maxime Grand, et al. Within-Host Microevolution of Pseudomonas aeruginosa Urinary Isolates: A Seven-Patient Longitudinal Genomic and Phenotypic Study. Frontiers in Microbiology, 2021, 11, pp.611246. 10.3389/fmicb.2020.611246. hal-03163322

HAL Id: hal-03163322

https://hal.umontpellier.fr/hal-03163322

Submitted on 9 Mar 2021

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

# Within-Host Microevolution of *Pseudomonas aeruginosa* Urinary Isolates: A Seven-Patient Longitudinal Genomic and Phenotypic Study

Agnès Cottalorda<fnref n="1" />*, Marie Leoz<fnref n="1" />, Sandrine Dahyot<fnref n="1" />,², François Gravey³, Maxime Grand<fnref n="1" />, Thomas Froidure<fnref n="1" />, Fabien Aujoulat⁴, Simon Le Hello³, Estelle Jumas-Bilak⁴ and Martine Pestel-Caron<fnref n="1" />,²

## OPEN ACCESS

**Edited by:**
Yasir Muhammad,
King Abdulaziz University,
Saudi Arabia

**Reviewed by:**
Jason Sahl,
Northern Arizona University,
United States
Volker Behrends,
University of Roehampton,
United Kingdom
Irene Bianconi,
University of Trento, Italy
Gloria Soberón-Chávez,
National Autonomous University
of Mexico, Mexico

*Correspondence:
Agnès Cottalorda
agnes.cottalorda@hotmail.fr

**Specialty section:**
This article was submitted to
Evolutionary and Genomic
Microbiology,
a section of the journal
Frontiers in Microbiology

Received: 28 September 2020

Accepted: 14 December 2020

Published: 14 January 2021

**Citation:**
Cottalorda A, Leoz M, Dahyot S, Gravey F, Grand M, Froidure T, Aujoulat F, Le Hello S, Jumas-Bilak E and Pestel-Caron M (2021) *Within-Host Microevolution of Pseudomonas aeruginosa Urinary Isolates: A Seven-Patient Longitudinal Genomic and Phenotypic Study*. *Front. Microbiol.* 11:611246.
doi: 10.3389/fmicb.2020.611246

**Background:** *Pseudomonas aeruginosa* is responsible for up to 10% of healthcare associated urinary tract infections (UTI), which can be difficult to treat and can lead to bacterial persistence. While numerous whole genome sequencing (WGS) analyses have explored within-host genomic adaptation and microevolution of *P. aeruginosa* during cystic fibrosis (CF) infections, little is known about *P. aeruginosa* adaptation to the urinary tract.

**Results:** Whole genome sequencing was performed on 108 *P. aeruginosa* urinary isolates, representing up to five isolates collected from 2 to 5 successive urine samples from seven patients hospitalized in a French hospital over 48-488 days. Clone type single nucleotide polymorphisms (ctSNPs) analysis revealed that each patient was colonized by a single clone type (<6000 SNPs between two isolates) at a given time and over time. However, 0-126 SNPs/genome/year were detected over time. Furthermore, large genomic deletions (1-5% of the genome) were identified in late isolates from three patients. For 2 of them, a convergent deletion of 70 genes was observed. Genomic adaptation (SNPs and deletion) occurred preferentially in genes encoding transcriptional regulators, two-component systems, and carbon compound catabolism. This genomic adaptation was significantly associated with a reduced fitness, particularly in artificial urine medium, but no strict correlation was identified between genomic adaptation and biofilm formation.

**Conclusion:** This study provides the first insight into *P. aeruginosa* within-host evolution in the urinary tract. It was driven by mutational mechanisms and genomic deletions and could lead to phenotypic changes in terms of fitness and biofilm production. Further metabolomic and phenotypic analyses are needed to describe in-depth genotype-phenotype associations in this complex and dynamic host-environment.

**Keywords:** *Pseudomonas aeruginosa*, urinary isolates, whole genome sequencing, within-host evolution, single nucleotide polymorphism, large genomic deletion

<footnote n="1">¹ Normandie Université, UNIROUEN, Groupe de Recherche sur l'Adaptation Microbienne (GRAM 2.0 EA2656), Rouen, France, ² Department of Bacteriology, Rouen University Hospital, Normandie Université, UNIROUEN, Groupe de Recherche sur l'Adaptation Microbienne (GRAM 2.0 EA2656), Rouen, France, ³ Normandie Université, UNICAEN, Groupe de Recherche sur l'Adaptation Microbienne (GRAM 2.0 EA2656), Caen, France, ⁴ Team Pathogènes Hydriques Santé Environnement, UMR 5569 HydroSciences Montpellier, Unité de Bactériologie, UFR Pharmacie, University of Montpellier, Montpellier, France</footnote>
# INTRODUCTION

Pseudomonas aeruginosa is an opportunistic Gram negative pathogen, responsible for 7–10% of healthcare associated urinary tract infections (UTIs) (Lamas Ferreiro et al., 2017). These UTIs preferentially occur in patients with underlying conditions, including patients with indwelling urinary catheter (Djordjevic et al., 2013; Gomila et al., 2018). P. aeruginosa can adhere and form biofilm on such surfaces (Olivares et al., 2019), which limits antibiotic diffusion and promotes emergence of antimicrobial resistant strains (Olivares et al., 2019; Soares et al., 2020). Moreover, this allows the microorganism to elude the immune system and persist causing frequent relapses (Gomila et al., 2018). Thus, readmissions are more frequent when UTIs are caused by P. aeruginosa (Gomila et al., 2018).

Advances in whole genome sequencing (WGS) techniques have made it possible to explore within-host genomic adaptation and microevolution of P. aeruginosa, particularly during cystic fibrosis (CF) infections (Marvig et al., 2015a; Bartell et al., 2019; La Rosa et al., 2019). Within-host evolution in CF patients is driven by three evolutionary modes, including general or directed diversity (Bartell et al., 2019), and convergent molecular evolution (Marvig et al., 2015b; Bartell et al., 2019). For example, the whole genome analysis of 474 longitudinal P. aeruginosa isolates collected from 34 CF patients led to the identification of a convergent evolution of 52 genes (Marvig et al., 2015b). These genes were involved in regulatory networks, central metabolism, antimicrobial resistance and secretion of virulence factors. Adaptive trajectories can also result in loss events which were described to be more frequent than gene acquisition (Rau et al., 2012; Bianconi et al., 2018) and could be characterized by large genomic deletions (Rau et al., 2012; Cabot et al., 2016; Hocquet et al., 2016).

Even though identification of genotype-phenotype links can be difficult for complex traits governed by multiple regulatory networks (Bartell et al., 2019), phenotypic changes have been previously described during persistent CF infections. They included a decreased secretion of virulence factors, the emergence of hypermutable strains often associated with an increased antimicrobial resistance, a switch from planktonic to sessile lifestyle, or a phenotypic conversion to mucoid or rugose small colony variants (Cabot et al., 2016; Hocquet et al., 2016; Winstanley et al., 2016; Klockgether and Tümmler, 2017; Pestrak et al., 2018; Colque et al., 2020).

In contrast, no study has explored to date P. aeruginosa within-host genomic and phenotypic evolution in the urinary tract. Indeed, most of the whole genome sequences of P. aeruginosa urinary isolates currently consist of genome announcements (Tielen et al., 2014; Wibberg et al., 2014, 2015; Vorhölter et al., 2015; Hussain et al., 2017).

In this context, we sequenced the whole genome of sequential urinary isolates of seven patients colonized and/or infected by P. aeruginosa to analyze the nature and extent of their genomic changes and provide an overview of within-host P. aeruginosa diversity at a given time and over time. We have also explored the

impact of genomic changes on phenotypic traits such as fitness and biofilm formation.

# MATERIALS AND METHODS

## Bacterial Isolates and Clinical Data

Two to five monomicrobial urine cultures (named I-V) of P. aeruginosa (regardless of the level of leucocyturia) were prospectively and consecutively collected from seven patients (A-G) over a 27-month period (from June 2016 to August 2018), with up to five isolates per sample (1–5), representing a total of 108 isolates (Table 1 and Figure 1). Isolates from the first sample collected were considered as early isolates while those from samples obtained at the subsequent times were considered as late isolates. Clinical data were collected from computerized medical and biological records of the hospital in parallel with the collection of the urinary isolates. According to the French regulation on observational database analyses, this study did not need specific informed consent requirements; it was registered by the Directorate of Clinical Research and Innovation of the Rouen University Hospital Centre under the number 2018/413/OB. The patients had episodes of UTI and/or asymptomatic bacteriuria (AB) averaging 192 days apart (ranging from 48 to 488 days).

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 analyses.

## Antibiotic Susceptibility Testing

Antibiotic susceptibility testing (AST) was performed for all isolates by disk diffusion method on Mueller-Hinton agar (Bio-Rad®, Marnes-la-Coquette, France) according to the criteria of the European Committee on Antimicrobial Susceptibility Testing<fnref n="1" />. Antibiotic disks were purchased from Bio-Rad®. Consensus recommendations (Magiorakos et al., 2012) were used to evaluate the proportion of Multi Drug-Resistant (MDR, not susceptible to at least three antimicrobial groups) and extensively Drug Resistant (XDR, not susceptible to at least six antimicrobial groups) isolates, considering the following six antimicrobial groups: penicillin (ticarcillin, ticarcillin-clavulanate, piperacillin, piperacillin-tazobactam), cephalosporins (ceftazidime, cefepime), monobactams (aztreonam), carbapenems (imipenem, meropenem), aminoglycosides (gentamicin, tobramycin, amikacin), and fluoroquinolones (ciprofloxacin, levofloxacin). The other isolates – resistant to none, one or two classes of antipseudomonal agents – were defined as non-MDR. The P. aeruginosa strain ATCC 27853 was used as a quality control strain in each experiment.

Of note in a parallel study (Cottalorda et al., unpublished), the antimicrobial resistance profiles, identified for isolates collected from the first urine sample, were also explored for these seven patients.

<footnote n="1">¹https://www.eucast.org/</footnote>
<table><caption>**TABLE 1** | Clinical context associated with the 108 *P. aeruginosa* urinary isolates studied.</caption><thead><tr><th>Patient</th><th>A</th><th>B</th><th>C</th><th>D</th><th>E</th><th>F</th><th>G</th></tr></thead><tbody><tr><td>Age range (years)</td><td>70-74</td><td>65-69</td><td>80-84</td><td>65-69</td><td>60-64</td><td>15-19</td><td>75-79</td></tr><tr><td>Urinary comorbidity<sup>a</sup></td><td>N</td><td>Y</td><td>Y</td><td>Y</td><td>Y</td><td>Y</td><td>Y</td></tr><tr><td>Immunocompromised conditions<sup>b</sup></td><td>Y</td><td>Y</td><td>N</td><td>Y</td><td>Y</td><td>N</td><td>N</td></tr><tr><td><strong>Clinical context [time]<sup>c</sup></strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Urine sample I</td><td><strong>UTI</strong> [D<sub>0</sub>]</td><td>AB<sup>*</sup> [D<sub>0</sub>]</td><td><strong>UTI</strong><sup>*</sup> [D<sub>0</sub>]</td><td>AB [D<sub>0</sub>]</td><td><strong>UTI</strong><sup>*</sup> [D<sub>0</sub>]</td><td>AB<sup>*</sup> [D<sub>0</sub>]</td><td><strong>UTI</strong><sup>*</sup> [D<sub>0</sub>]</td></tr><tr><td>Urine sample II</td><td><strong>UTI</strong> [D<sub>14</sub>]</td><td>AB<sup>*</sup> [D<sub>33</sub>]</td><td><strong>UTI</strong><sup>*</sup> [D<sub>48</sub>]</td><td>AB [D<sub>259</sub>]</td><td>AB [D<sub>28</sub>]</td><td><strong>UTI</strong><sup>*</sup> [D<sub>113</sub>]</td><td>AB<sup>*</sup> [D<sub>46</sub>]</td></tr><tr><td>Urine sample III</td><td>AB [D<sub>82</sub>]</td><td>AB<sup>*</sup> [D<sub>296</sub>]</td><td></td><td></td><td>AB [D<sub>70</sub>]</td><td></td><td>AB<sup>*</sup> [D<sub>90</sub>]</td></tr><tr><td>Urine sample IV</td><td></td><td>AB<sup>*</sup> [D<sub>362</sub>]</td><td></td><td></td><td>AB [D<sub>81</sub>]</td><td></td><td><strong>UTI</strong><sup>*</sup> [D<sub>271</sub>]</td></tr><tr><td>Urine sample V</td><td></td><td>AB<sup>*</sup> [D<sub>488</sub>]</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Letters (A–G) identify patients. Roman numerals (I–V) correspond to urine sample numbers in the series collected from a given patient. <sup>a</sup>Urinary comorbidity includes neurogenic bladder, benign prostate hypertrophy, urinary tract diversion, and renal transplant. <sup>b</sup>Immunocompromised conditions include immunodepression (such as AIDS, lymphoma, leukemia) or immunosuppressive treatment. <sup>c</sup>Days (D) between a given urine sample and the first one. *Indicates urinary catheterization. Y, yes; N, no; UTI, urinary tract infection; AB, asymptomatic bacteriuria.

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

# DNA Extraction, Library Construction, Genome Sequencing, and De novo Assembly

## PacBio Library Construction and Sequencing

For each patient, the first sample isolate with the higher rate of antimicrobial resistance (named X-I-1 where X was a capital letter from A to G) was sequenced with the Pacific Biosciences (PacBio) technology. It was considered as the intra-patient reference sequence. Genomic DNA was extracted using the NucleoSpin Microbial DNA kit (Macherey-Nagel, Illkirch-Graffenstaden, France) according to the manufacturer's recommendations. The SMRTbell library was prepared using a SMRTbell Express 2 Template prep kit (Pacific Biosciences®, California Inc., United States), following the “procedure and

checklist preparing Multiplexed Microbial Libraries using SMRTbell Express Template prep kit 2.0” protocol. Single-molecule Real-time long reads sequencing was performed at the Gentyane Sequencing Platform (Clermont-Ferrand, France) with a PacBio Sequel Sequencer (Pacific Biosciences®, Menlo Park, CA, United States). De novo assembly was performed using “Microbial Assembly” pipeline in PacBio SMRT Link v8.0.0.80529, using default parameters, except “Genome Size,” which was set to X.

## Illumina NextSeq 500 Library Construction

The genome of the 101 remaining isolates was sequenced by the “Plateforme de Microbiologie Mutualisée P2M” (Institut Pasteur, Paris, France) using the MagNA Pure 96 Instrument (Roche Diagnostics, Meylan, France) for genomic DNA extraction,
and the Nextera XT library kit (Illumina Inc.), the NextSeq 500 sequencing system (Illumina Inc.). Paired-end reads were submitted to pre-processing using `fqCleaner` and to *de novo* assembly using SPADES, with a 50X minimum average sequencing depth.

## Genome Annotation, Genotyping, Clone Type Definition

Genome annotation was performed for the seven isolates sequenced by PacBio using the Prokaryotic Genome Annotation Pipeline (PGAP) version 2019-08-22 (Tatusova et al., 2016)<fnref n="2" />.

*In silico* Multilocus Sequence Typing (MLST) genotyping was performed according to the scheme established by Curran et al. (2004), using KMA software version 1.2.12 (Clausen et al., 2018) with default settings<fnref n="3" />, and the online database (updated on 2019-03-18) of the Center for Genomic Epidemiology (Larsen et al., 2012).

Phylogenetic analysis of the 108 Illumina and PacBio assemblies produced by SPADES was performed using Parsnp software version 1.2 (Treangen et al., 2014) with default parameters using the A-I-1 strain as reference. The iTOL software (Letunic and Bork, 2019) was used to visualize and edit the phylogenetic tree and identify phylogenetic clusters (bootstrap of 100) inside the *P. aeruginosa* urinary population.

Illumina reads were aligned against every PacBio reference sequences using Snippy software version 4.4.0 (available on https://github.com/tseemann/snippy/). Clone type SNPs (ctSNPs), bounded by at least 50 exact base-pair matches on both sides, were extracted from the alignments using a custom script<fnref n="4" />. Isolates were considered to belong to different clone types when more than 6 000 ctSNPs were found between their genomes (Marvig et al., 2015b).

PlasmidFinder 2.1, available on the Center for Genomic Epidemiology database<fnref n="5" /> was used to identify plasmid in genome sequences.

## SNPs Analysis

An intra-patient SNPs analysis was performed using Snippy software, by mapping Illumina reads to their corresponding PacBio reference sequence. This intra-patient mapping approach was chosen over the fast Parsnp analysis of genome assemblies (well fit for inter-patient scaled phylogenies) because it allowed a fine-scale analysis of the SNPs from raw reads, and maximized the size of the intra-patient core genomes. SNPs analysis was manually curated. When SNPs were found in all the isolates sequenced by Illumina, they were considered specific to their common PacBio reference sequence even though we cannot exclude that some were in fact artifacts due to the use of distinct sequencing approaches (PacBio vs. Illumina). Snippy called the SnpEff software (Cingolani et al., 2012) to annotate and predict the effect of each mutation. To identify the corresponding PAO1 locus tags, gene sequences from each PacBio reference

genome were blasted against all PAO1 gene sequences retrieved from https://www.pseudomonas.com/. Gene sequences that with no blast hit against PAO1 were classified as “hypothetical, unclassified, unknown.” Genes where SNPs occurred were classified into functional categories using PseudoCAP available on the *Pseudomonas* Genome Database<fnref n="6" />. SNPs were classified as low (e.g., synonymous mutations), moderate (e.g., missense mutations), and high effect (e.g., stop gained). SNPs that were not identified in early isolates (collected from the first urine sample for each patient) but present in at least two of late isolates were further described as “late stage SNPs.”

The core alignments produced by Snippy software for each patient were used to infer within-host phylogenies: Neighbor Joining trees were produced using MEGA version X (Kumar et al., 2018), with a Tamura Nei evolutionary model, rate heterogeneity following a Gamma distribution (parameter = 1), and complete deletion of gaps and missing data.

## Characterization of Genomic Deletions

Genomic deletions were visualized using CGView software (Stothard and Wishart, 2005), based on the core alignments produced by Snippy. Deletions were also identified using KMA software version 1.2.12 (see text footnote 3), with Illumina reads as query and the CDS identified by pgap from their corresponding PacBio reference genome as database. The proportion of the core and accessory genome were calculated based on the core alignments produced by Snippy. To identify the corresponding PAO1 locus tags, gene sequences from each PacBio reference genome were blasted against all PAO1 gene sequences retrieved from https://www.pseudomonas.com/. Gene sequences that with no blast hit against PAO1 were classified as “hypothetical, unclassified, unknown.” The deleted genes were classified into functional categories using PseudoCAP available on the *Pseudomonas* Genome Database (see text footnote 6).

Genomic deletions were confirmed by PCR amplification. Briefly, genomic DNA was extracted using the InstaGene™ Matrix kit (Bio-Rad) following the manufacturer's recommendations. PCRs were then performed on a Veriti Thermal Cycler (Applied Biosystems, Foster City, CA, United States) in a final volume of 25 µL containing 12.5 µL GoTaq® Green Master Mix (Promega, Charbonnières-les-Bains, France), 0.50 mM of each primer (Supplementary Table 1) and 5 µL of DNA. In each case, three primer pairs were used: two that bordered each side of the potentially deleted fragment and allowed the amplification for isolates without genomic deletion, and one outside the deleted region, which allowed the amplification for isolates with genomic deletion. Two genomes without deletion [the intra-patient reference (X-I-1) and when possible, an isolate collected from the last urine sample] obtained from the same patient were used as a control.

## Biofilm Formation Assay

Biofilm formation was assessed by a crystal violet staining assay as previously described (O'Toole, 2011) with minor modifications. Briefly, strains were grown overnight in trypticase soya broth

<footnote n="2">²https://github.com/ncbi/pgap/</footnote>

<footnote n="3">³https://bitbucket.org/genomicepidemiology/kma.git</footnote>

<footnote n="4">⁴https://github.com/marieleoz/pairwise-SNPs-using-Snippy</footnote>

<footnote n="5">⁵https://cge.cbs.dtu.dk/services/PlasmidFinder/</footnote>

<footnote n="6">⁶https://www.pseudomonas.com/</footnote>
(Thermo Fisher Scientific, Les Ulis, France) supplemented with 1% glucose (w/v) (VWR, Leuven, Belgium) (Brooks and Keevil, 1997). Overnight cultures were adjusted to an optical density at 600nm (OD<sub>600nm</sub>) of 0.01. In a 96-well plastic microtiter plate (Thermo Fisher Scientific, Illkirch-Graffenstaden, France), 100 µL of the standardized samples were inoculated per well in triplicate. After 24 h of static incubation at 37°C, supernatants were removed, and the cells that remained attached to the microtiter wells were washed twice with water by submersion. Cells were stained with 125 µL of an aqueous solution of 0.1% crystal violet. After 15 min of incubation at room temperature, the plate was rinsed thrice with water by submersion and dried for 2 h. Cells were washed and resuspended in 50% (v/v) ethanol. OD<sub>590nm</sub> of the resuspended biofilm was determined and normalized to the OD<sub>600nm</sub> of the corresponding planktonic cells. All experiments were performed at least three times. The *P. aeruginosa* strain PA14 was used as a positive control and fresh medium as a negative control for each experiment. All data were expressed as percent of biofilm formation relative to the positive control. Isolates forming less than 50% of biofilm compared to PA14 were considered as low biofilm-producers, whereas isolates forming more than 50% were considered as high biofilm-producers. Statistical analyses were performed by Kruskal-Wallis test on R® software (v.3.5.1).

## Fitness Assays

Growth curves were performed in two distinct media: trypticase soya broth, and artificial urine medium (AUM) according to Brooks and Keevil' protocol (Brooks and Keevil, 1997) with minor modifications. Briefly, 2-(N-morpholino)ethanesulfonic acid (MES buffer) at 100 mM was added to the previously described AUM in order to limit precipitates formation. Overnight cultures were adjusted to an OD<sub>600nm</sub> of 0.01. In a 96-well plastic microtiter plate, 200 µL of standardized samples were inoculated per well in triplicate. The plate was placed inside a humidity cassette and incubated at 37°C for 24 h under continuous double-orbital shaking (108 rpm) into the *Spark®* multimode microplate reader (Tecan, Männedorf, Switzerland). OD<sub>600nm</sub> was measured every 15 min. All experiments were performed at least three times. The *P. aeruginosa* strain PA14 was used as an internal control and fresh medium as a negative control for each experiment. Statistical analyses were performed by Kruskal-Wallis test on R® software (v.3.5.1). The generation time ratio was defined for a given patient as the mean generation time of isolates with genomic deletion on the mean generation time of isolates without one.

## RESULTS

### Clinical Data

During a 27-month period, seven hospitalized patients (A-G) with more than one urine sample positive to *P. aeruginosa* were included in this study (Figure 1).

Chronic asymptomatic bacteriuria (AB) was observed for two patients (B and D) while chronic UTI was observed for one patient (C). A switch from UTI to AB (patients A and E), from

AB to UTI (patient F) or successive switches from UTI and AB (patient G) were identified (Table 1). Of note, four patients (B, C, F, and G) required long term urinary catheterization. Six patients had urinary disorders, including benign prostate hypertrophy (C, D, G), urinary tract diversion (B, F), and kidney transplant (E). All patients received at least one antibiotic treatment during the 6 months before urine collection, and five of them (A, B, C, E, and G) also received antibiotics in the interval of their urine samples (Table 1 and Figure 2A and Supplementary Table 2).

## Antimicrobial Resistance Profiles

Isolates from the first urine sample exhibited a non-MDR phenotype for four patients (C, E, F, G), a MDR phenotype for one patient (D), and a XDR phenotype for one patient (B). Within-sample diversity was observed between isolates collected from patient A (four MDR-isolates and one XDR) (Figure 2B).

Over time, the rate of antimicrobial resistance was similar for isolates of two patients (B, E) while it decreased for patient F, and increased for patients C and G (Figure 2B). Within-sample diversity of antimicrobial resistance profiles increased over time for three patients (A, D, G) (Figure 2B and Supplementary Table 3). No plasmid that could explain the acquisition of these resistance profiles have been identified in the genome sequences.

Surprisingly, no strict association was identified between antimicrobial resistance (Figure 2B) and previous antibiotic therapy (Figure 2A). Indeed, resistance to a given group of antibiotics was not systematically associated with previous antimicrobial therapy by one molecule of that group, suggesting acquisition of cross-resistance. Similarly, the increase of antimicrobial resistance of late isolates was not completely explained by antibiotic treatments during the period between urine samples.

## Genotyping and Phylogeny of *P. aeruginosa* Urinary Population

Sequencing data are presented in Supplementary Table 4. The number of predicted genes per genome sequenced by PacBio and annotated with PGAP ranged from 5,991 to 6,886.

In silico MLST analysis performed for the 108 urinary isolates revealed six sequence types (STs): ST111, ST232, ST235, ST308, ST633, and ST3227. All isolates from a given patient, within a given sample but also over time, belonged to the same ST. Two patients (D and E) had isolates with the same ST (ST308). The phylogenetic analysis of the 108 core genomes (draft assemblies) led to the identification of six clusters that were congruent with the 6 STs identified (Supplementary Figure 1A).

Clone type single nucleotide polymorphisms analysis led to the identification of six clone types, in accordance with the STs identified. The average number of ctSNPs between isolates from a given patient was seven (from 0 to 123) while it was 26,242 between patients (from 412 to 33,638) (Supplementary Table 5). Of note, ST308 genomes belonged to a single clone type (number of ctSNPs between patients D and E from 412 to 1,212), even though isolates were collected from distinct patients. However, there was no clinical support for a cross infection between patients D and E, since they were hospitalized in separate units at
distinct times and were not epidemiologically linked. Sequences sampled from each patient formed two distinct and well supported phylogenetic clusters (Supplementary Figure 1B).

## Within-Host Evolution by SNPs

To characterize over-time evolution of within-host single-clone types, we investigated all SNPs identified for each isolate, depending on the time of urine sampling (Supplementary Table 6). For the seven patients, the mean number of SNPs was higher for the late isolates (8.6 ± 25) than for the isolates collected

from the first urine sample (0.8 ± 1). The mutation rate ranged from 0 to 126 SNPs per genome per year (Table 2).

D-isolates had the highest number of SNPs, with an average of 2 SNPs per isolate in the first urine sample, increasing to 90 in the second one (Table 2). Phylogenetic analysis of the core-genome sequences sampled from this patient revealed distinct sublineages in the second urine sample, consistent with the high number of SNPs observed (Figure 3). Indeed, one late isolate (D-II-1) was phylogenetically related to the early isolates, and did not share any SNP with the other four late isolates that clustered separately.

<figure><img src="image_4.png" /><figcaption>TABLE 2 | Average number of SNPs per urine sample for the seven patients studied.</figcaption></figure>

<table><thead><tr><th>Patient</th><th></th><th>A</th><th>B</th><th>C</th><th>D</th><th>E</th><th>F</th><th>G</th></tr></thead><tbody><tr><td>Number of positions investigated in each core genome</td><td></td><td>6,598,293</td><td>6,803,544</td><td>6,266,825</td><td>6,854,231</td><td>6,890,906</td><td>6,472,838</td><td>6,416,922</td></tr><tr><td rowspan="5">Urine culture</td><td>I</td><td>0.8</td><td>1.4</td><td>0.2</td><td>2</td><td>0.6</td><td>0</td><td>0.4</td></tr><tr><td>II</td><td>0.2 [10]</td><td>2.8 [31]</td><td>1.2 [9]</td><td>89.6 [126]</td><td>0 [0]</td><td>7.4 [24]</td><td>0 [0]</td></tr><tr><td>III</td><td>0.6 [3]</td><td>6.4 [8]</td><td></td><td></td><td>0.4 [2]</td><td></td><td>0 [0]</td></tr><tr><td>IV</td><td></td><td>9 [9]</td><td></td><td></td><td>0.2 [1]</td><td></td><td>2.7 [4]</td></tr><tr><td>V</td><td></td><td>5.8 [4]</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Illumina reads were mapped to the corresponding PacBio reference using Snippy software. The total number of positions investigated in each intra-patient core genome is indicated. The mean number of SNPs/year (mutation rate, in square brackets) was calculated for isolates from each urine culture, based on the time elapsed since the previous urine culture.
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Phylogenetic analysis of the isolates collected from patient D. The core-genome alignment produced by snippy was used for phylogenetic inference using MEGA software. The phylogenetic tree was rooted on the PacBio reference sequence D-I-1. Isolates with genomic deletion are identified by a red triangle. Bootstrap values are indicated when >70. The length of the scale bar represents 0.0000020 substitutions per site, corresponding to 13.71 SNPs in this 6,854,231 nucleotide-long core genome alignment.</figcaption></figure>

Interestingly, no mutation was detected in DNA mismatch repair genes (like *mutS* or *mutL*) but 2 SNPs were identified in genes involved in DNA replication, recombination, modification and repair (a missense SNP, Asp559Asn, in *uvrA*, and a nonsense SNP at position 309/530 in the gene homologous of PA1405 at locus tag H7Q90_28210).

To identify genetic mutations that could be involved in phenotypic adaptation, we focused on late stage SNPs (i.e., SNPs that were not identified in early isolates but present in at least two of late isolates) with moderate (non-synonymous mutations) or high (stop gained) effect.

No such SNP was detected in isolates collected from three patients (A, C, and E), while 2–79 late stage SNPs were identified for the other patients (B, D, F, and G). No pathoadaptive genes reflecting convergent evolution by these SNPs were identified for the isolates collected from the four patients. However, SNPs preferentially occurred in genes encoding transcriptional regulators or involved in two component systems (Table 3). However, no pathoadaptive genes reflecting convergent evolution by SNPs was identified between isolates collected from distinct patients.

Single nucleotide polymorphisms were identified in genes involved in virulence for isolates collected from patient D, F, and G. For D-isolates, eight SNPs were identified in genes involved in motility and adherence (*fleS*, *fimT*), in type III and VI secretion systems [*exsA* (nonsense SNP), *icmF1*, *hsiG3*], in ferric uptake regulator (FUR) regulon (PA2094 and PA2402), or in stringent response (antitoxin *higA*). For F-isolates, three SNPs were identified in genes involved in alginate biosynthesis or alginate regulation and biofilm formation [*algG*, *algB*, and *mucA* (nonsense SNP)]. These isolates presented a mucoid morphotype after 48 h of culture (data not shown), which could be explained by *mucA* inactivation (an anti-sigma factor of the *algU* gene involved in alginate biosynthesis). Finally, one SNP was observed in the *fecI* sigma factor of the FUR regulon in G-isolates.

Of note, SNP occurring in genes involved in antimicrobial resistance were only identified for D-isolates (*mexS*).

## Within-Host Evolution by Genomic Deletion

For three patients (A, D, and F), at least two of the late isolates had a smaller genome size than isolates collected from the first urine sample (Supplementary Table 7), as shown by the read mapping on the intra-patient reference sequence (Figures 4A,B). Deletions represented up to 5% of the size of the reference genome for A-isolates, and about 1% for D- and F-isolates. All these genomic deletions were experimentally confirmed by PCR (Supplementary Table 1).

The largest genomic deletion (361 kbp) was observed in two late A-isolates (A-III-3 and -5), and led to the loss of 375 coding protein sequences (CDS). Genomic deletions of 35 kbp (48 CDS) and 101 kbp (96 CDS) and were observed for D-II-2 to D-II-5 isolates and F-II-1 to F-II-5 isolates, respectively. The deleted fragment of D- and F-isolates was immediately bordered by direct repeats (DR). For the F-isolates, the 1677 bp DR sequence corresponded to an insertion sequence (IS) that belongs to the IS1182 transposase family. Interestingly, the same IS was found at seven positions in the F-I-1 genome. Furthermore, the deletions in A- and D-isolates led to the loss of 11 and 19 phage related CDS, respectively.

For A- and F-isolates, 70 deleted CDS were common (homologous to PA1965 to PA2034 genes), suggesting a convergent evolution (Table 4). This region included genes involved in antimicrobial resistance (*mexX/Y/Z* genes), regulatory systems (transcriptional regulators or two component systems), carbon compound catabolism (*liu*, *exa*), and the genes *hmgA* (encoding a homogentisate 1,2-dioxygenase) and *galU* (encoding an UDP-glucose phosphorylase involved in lipopolysaccharide biosynthesis) (Cabot et al., 2016;
<table><caption>**TABLE 3** | PseudoCAP function classification of genes affected by late stage SNPs with moderate or high effect identified in isolates collected from four patients (B, D, F, G).</caption><thead><tr><th>PseudoCAP manually assigned functions</th><th>B- isolates</th><th>D- isolates</th><th>F- isolates</th><th>G- isolates</th></tr></thead><tbody><tr><td>Adaptation, protection</td><td>0%</td><td>2%</td><td>11%</td><td>0%</td></tr><tr><td>Amino acid biosynthesis and metabolism</td><td>0%</td><td>2%</td><td>0%</td><td>0%</td></tr><tr><td>Antibiotic resistance and susceptibility</td><td>0%</td><td>2%</td><td>0%</td><td>0%</td></tr><tr><td>Biosynthesis of cofactors, prosthetic groups and carriers</td><td>0%</td><td>2%</td><td>0%</td><td>0%</td></tr><tr><td>Carbon compound catabolism</td><td>0%</td><td>4%</td><td>0%</td><td>33%</td></tr><tr><td>Cell division</td><td>0%</td><td>2%</td><td>0%</td><td>0%</td></tr><tr><td>Cell wall/LPS/capsule</td><td>0%</td><td>4%</td><td>33%</td><td>0%</td></tr><tr><td>Central intermediary metabolism</td><td>0%</td><td>1%</td><td>0%</td><td>0%</td></tr><tr><td>Chemotaxis</td><td>33%</td><td>1%</td><td>0%</td><td>0%</td></tr><tr><td>DNA replication, recombination, modification and repair</td><td>0%</td><td>2%</td><td>0%</td><td>0%</td></tr><tr><td>Energy metabolism</td><td>0%</td><td>1%</td><td>0%</td><td>0%</td></tr><tr><td>Fatty acid and phospholipid metabolism</td><td>0%</td><td>0%</td><td>0%</td><td>33%</td></tr><tr><td>Membrane proteins</td><td>0%</td><td>13%</td><td>0%</td><td>0%</td></tr><tr><td>Motility and Attachment</td><td>33%</td><td>2%</td><td>0%</td><td>0%</td></tr><tr><td>Nucleotide biosynthesis and metabolism</td><td>0%</td><td>2%</td><td>0%</td><td>0%</td></tr><tr><td>Protein secretion/export apparatus</td><td>0%</td><td>3%</td><td>0%</td><td>0%</td></tr><tr><td>Putative enzymes</td><td>0%</td><td>6%</td><td>0%</td><td>0%</td></tr><tr><td>Secreted Factors (toxins, enzymes, alginate)</td><td>0%</td><td>1%</td><td>11%</td><td>0%</td></tr><tr><td>Transcription, RNA processing and degradation</td><td>0%</td><td>3%</td><td>0%</td><td>0%</td></tr><tr><td>Transcriptional regulators</td><td>0%</td><td>7%</td><td>22%</td><td>33%</td></tr><tr><td>Translation, post-translational modification, degradation</td><td>0%</td><td>4%</td><td>0%</td><td>0%</td></tr><tr><td>Transport of small molecules</td><td>0%</td><td>5%</td><td>0%</td><td>0%</td></tr><tr><td>Two-component regulatory systems</td><td>33%</td><td>3%</td><td>11%</td><td>0%</td></tr><tr><td>Hypothetical, unclassified, unknown</td><td>0%</td><td>27%</td><td>11%</td><td>0%</td></tr></tbody></table>

Percent of SNPs identified per gene function are presented here. PseudoCAP functions are available on the Pseudomonas Genome Database (https://www.pseudomonas.com/). When over than 5% of genes for a given PseudoCAP function were affected (except for "Hypothetical, unclassified, unknown" function), values were indicated in bold.

López-Causapé et al., 2018; Wardell et al., 2019). Interestingly, isolates with the *hmgA* deletion produced a red pigment (pyomelanin) when they were grown 24 h on TSA (data not shown), as previously reported (Cabot et al., 2016).

Moreover, in A-isolates, the deleted region contained genes encoding EAL domain-containing proteins, which are involved in c-di-GMP degradation (a second messenger involved in biofilm formation), and genes involved in virulence (*hcnABC*, encoding an hydrogen cyanide synthase; the effector *exoY* of the type III secretion system; and *cupA1* to *cupA5*, involved in *fimbriae* synthesis).

## Phenotypic Analyses

Fitness and biofilm formation assays were performed for A-, D-, and F-isolates (n = 35 isolates) for which genomic deletions occurred over time.

A- and D-isolates were low biofilm-producers (Figure 5). An increased biofilm formation was identified for A-Δ isolates compared to the others (*P* < 0.05). In contrast, the ability of the D-Δ isolates was reduced (*P* < 0.05). In parallel, early F-isolates were high biofilm-producers while deleted late F-isolates produced significantly less biofilm (*P* < 0.01).

Isolates with a genomic deletion had a significant lower fitness in TSB broth (for D- and F-isolates) and in AUM (for isolates from all patients) (Figure 6). The reverse trend was observed for A-isolates in TSB broth (*P* < 0.05). Interestingly, the generation time ratio (defined as the mean generation time of isolates with genomic deletion on the mean generation time of isolates without one for a given patient) was near to 2 in AUM (1.7 for A-isolates, 2.3 for D-isolates, and 1.9 for F-isolates). Thus, isolates with genomic deletion grew about half as fast as isolates without a genomic deletion in AUM.

## DISCUSSION

*Pseudomonas aeruginosa* is a major cause of healthcare associated UTI, characterized by frequent relapses (Djordjevic et al., 2013; Gomila et al., 2018). Numerous studies explored mechanisms of *P. aeruginosa* persistence during chronic pulmonary infections in CF patients (Caballero et al., 2015; Marvig et al., 2015b; Bianconi et al., 2018; Bartell et al., 2019), but little is known about its pathoadaptive adaptation in the urinary tract. To this end, up to five isolates of *P. aeruginosa* were collected in longitudinal urine samples from seven hospitalized patients.

All patients included in our study were colonized or infected over time by a single clone type, well-predicted by *in silico* MLST. Nevertheless, distinct sublineages have been identified in the second urine sample collected from one patient (D), suggesting a progressive replacement of the naive initial *P. aeruginosa* population over time. This is consistent with longitudinal studies of CF infections, which showed that chronic *P. aeruginosa* infections were often associated with diversification to distinct variants (Marvig et al., 2015b; Bianconi et al., 2018; Bartell et al., 2019). Even though our sampling period was shorter (488 days at most) than these studies (where inclusion could occurred nearly a decade apart), Bartell et al. (2019) pointed that rapid adaptation occur within 2–3 years timeline after colonization in CF patients.

To our knowledge, this is the first comprehensive whole genome analysis of within-host evolution of *P. aeruginosa* in the urinary tract. It has highlighted the molecular mechanisms of its adaptation, which was mostly driven by two distinct mechanisms: SNPs and large genomic deletions.

Single nucleotide polymorphisms analysis revealed that an average of less than 10 SNPs differentiated early from late isolates for a given patient. The number of SNPs identified was much lower than that detected for isolates collected during CF infections (up to 27,249 SNPs per strain in Bianconi et al.'s study) (Bianconi et al., 2018). This is probably because we used here
A

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

B

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

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

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

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

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

**FIGURE 4** | Genomic deletions identified in late isolates collected from three patients (A, D, F). (A) Genome coverage for isolates collected from three patients (A, D, F). Genome coverage is represented according to Illumina read mapping by Snippy on their corresponding PacBio reference sequence. Green, blue, and purple colors indicate that the isolates are from sample I, II, and III, respectively. Samples and genomic regions where large deletions occurred in late isolates are framed by red rectangles. (B) Proportion of the core and accessory genomes for isolates collected from three patients (A, D, F). The proportions were calculated based on the core alignments produced by Snippy. The size of the core genome was 6,598,293 bp, 6,854,231 bp, and 6,472,838 bp, for A-, D-, and F-isolates, respectively.

<table><caption>**TABLE 4** | PseudoCAP function classification of deleted genes identified in isolates collected from three patients (A, D, F).</caption><thead><tr><th>PseudoCAP manually assigned functions</th><th>A-isolates</th><th>D-isolates</th><th>F-isolates</th><th>A- and F- isolates</th></tr></thead><tbody><tr><td>Adaptation, Protection</td><td>1%</td><td>0%</td><td>1%</td><td>0%</td></tr><tr><td>Amino acid biosynthesis and metabolism</td><td>2%</td><td>0%</td><td>3%</td><td>4%</td></tr><tr><td>Antibiotic resistance and susceptibility</td><td>1%</td><td>0%</td><td>3%</td><td>2%</td></tr><tr><td>Biosynthesis of cofactors, prosthetic groups and carriers</td><td>2%</td><td>0%</td><td><strong>6%</strong></td><td><strong>8%</strong></td></tr><tr><td>Carbon compound catabolism</td><td><strong>5%</strong></td><td>0%</td><td><strong>13%</strong></td><td><strong>16%</strong></td></tr><tr><td>Cell wall/LPS/Capsule</td><td>0%</td><td>0%</td><td>1%</td><td>0%</td></tr><tr><td>Central intermediary metabolism</td><td>2%</td><td>0%</td><td>2%</td><td>2%</td></tr><tr><td>Chaperones and heat shock proteins</td><td>1%</td><td>0%</td><td>1%</td><td>1%</td></tr><tr><td>DNA replication, recombination, modification and repair</td><td>1%</td><td>0%</td><td>0%</td><td>0%</td></tr><tr><td>Energy metabolism</td><td>1%</td><td>0%</td><td>1%</td><td>1%</td></tr><tr><td>Fatty acid and phospholipid metabolism</td><td>1%</td><td>0%</td><td>1%</td><td>1%</td></tr><tr><td>Membrane protein</td><td><strong>7%</strong></td><td>0%</td><td><strong>9%</strong></td><td><strong>10%</strong></td></tr><tr><td>Motility and Attachment</td><td>1%</td><td>0%</td><td><strong>5%</strong></td><td>0%</td></tr><tr><td>Nucleotide biosynthesis and metabolism</td><td>0%</td><td>0%</td><td>0%</td><td>0%</td></tr><tr><td>Putative enzymes</td><td><strong>7%</strong></td><td>0%</td><td><strong>7%</strong></td><td><strong>8%</strong></td></tr><tr><td>Secreted Factors (toxins, enzymes, alginate)</td><td>0%</td><td>0%</td><td>0%</td><td>0%</td></tr><tr><td>Transcriptional regulators</td><td><strong>7%</strong></td><td>0%</td><td><strong>11%</strong></td><td><strong>11%</strong></td></tr><tr><td>Translation, post-translational modification, degradation</td><td>0%</td><td>0%</td><td>1%</td><td>1%</td></tr><tr><td>Transport of small molecules</td><td><strong>7%</strong></td><td>0%</td><td><strong>9%</strong></td><td><strong>6%</strong></td></tr><tr><td>Two-component regulatory systems</td><td>1%</td><td>0%</td><td>4%</td><td><strong>5%</strong></td></tr><tr><td>Hypothetical, unclassified, unknown</td><td>52%</td><td>100%</td><td>24%</td><td>23%</td></tr></tbody></table>

Percent of deleted genes per gene function are presented here. PseudoCAP functions are available on the Pseudomonas Genome Database (https://www.pseudomonas.com/). A convergent evolution by deletion was identified in a common region between A- and F- isolates. Percent of deleted genes per gene function are indicated by "A- and F- isolates." When over than 5% of genes for a given PseudoCAP function were affected (except for "Hypothetical, unclassified, unknown" function), values were indicated in bold.
<figure><img src="image_12.png" /><figcaption>**FIGURE 5** | Mean biofilm formation of 35 *P. aeruginosa* urinary isolates collected from three patients (A, D, and F) compared to PA14. Biofilm formation was evaluated by crystal violet in microtiter plates and was compared between isolates with genomic deletions (-Δ) and without genomic deletion for each patients. A- Δ: A-III-3, and A-III-5 isolates; D- Δ: D-II-2 to D-II-5 isolates; F- Δ: F-II-1 to F-II-5 isolates. Significance was identified as *P* < 0.05 (*) or *P* < 0.01 (**) by Kruskal-Wallis test.</figcaption></figure>

the intra-patient reference sequence to limit the phylogenetic
distance between isolates and the reference genome while CF
studies usually use a reference sequence (*P. aeruginosa* DK2 or
PAO1) (Marvig et al., 2015b; Bianconi et al., 2018). Nevertheless,
we cannot exclude that the first urine sample of each patient
included here did not correspond to the first episode of UTI
or AB and therefore we could have underestimated within-
host evolution.

However, in order to quantitatively compare evolution by
SNPs over time, we have evaluated the number of SNPs per
genome per year. For most patients, the mean number of SNPs
per genome per year (ranging from 0 to 126) could have been
higher than that previously described for chronic CF infections
(1-3) (Marvig et al., 2015a). Interestingly, the highest rate
observed for one patient (126, patient D) reached that of 89-
106 identified for hypermutator strains (Marvig et al., 2015a).
For this patient, no SNP was identified in DNA mismatch repair
genes (*mutS* and *mutL*) that often characterize hypermutator
strains (Colque et al., 2020), but two SNPs in genes involved
in DNA replication, recombination, modification and repair
were identified, as previously described (Qin et al., 2018). The
hypermutable phenotype would be interesting to confirm *in vitro*
in further analyses.

Late stage SNPs with moderate or high effect were identified
in four patients, all from urine samples collected more than
100 days apart (B, D, F, G). In contrast to CF infections (Marvig
et al., 2015b), no pathoadaptive genes that could suggest a
convergent evolution were identified. This could be due to
the heterogeneity of the patients included here, particularly in
terms of age, catheterization, and urinary disorders, which could
greatly influence the bacterial environment and *P. aeruginosa*
within-host adaptation. Nevertheless, SNPs in genes encoding

transcriptional regulators and two component systems were
identified in three of the four patients, suggesting that
remodeling of regulatory networks might be important in
*P. aeruginosa* adaptation to the urinary tract. Similarly, 10 of
the 52 pathoadaptive genes identified during a longitudinal
study in CF infections encoded transcriptional regulators
(Marvig et al., 2015b).

In addition to mutations, and similarly to Bianconi et al. (2018), we identified that bacterial persistence could be associated with loss events. Indeed, genomic deletions were identified in isolates collected from three of the seven patients studied (A, D, F) over up to 100 days apart. The size of the genomic deletions observed in our study (35–361 kbp) was slightly smaller than those previously reported during CF infections, which could reach up to 8% (525 kbp) of the genome (Rau et al., 2012).

Analysis of the deleted genes showed a convergent evolution
on 70 CDS for two patients (A, F). Among them, deletion of
the *hmgA* and *galU* genes have been previously reported *in vitro*
(after meropenem exposure) or *in vivo* particularly during CF
infections (Rau et al., 2012; Cabot et al., 2016; Hocquet et al., 2016;
Wardell et al., 2019). Deletion of *hmgA* resulted in pyomelanin
hyperproduction, which was also phenotypically observed for the
isolates of our study that presented the deletion of this gene.
The *galU* gene is involved in polysaccharide biosynthesis, and
its deleted mutant was shown to be devoid of O-antigen, and
promoted persistence of *P. aeruginosa* in chronic lung infection
by immune evasion (Rodríguez-Rojas et al., 2009; Hocquet et al.,
2016). Furthermore, as *galU* has previously been described as
necessary for phage adsorption, its deletion led to a phage
resistance mechanism (Le et al., 2014).

The convergent evolution in patients A and F led to the
deletion of genes PA1965 to PA2034. Interestingly, this genomic
<figure><img src="image_13.png" /><figcaption>**FIGURE 6** | Generation time of 35 *P. aeruginosa* urinary isolates collected from three patients (A, D, and F) compared to PA14. Generation time (G) (min) in (A) TSB broth and (B) artificial urine medium of isolates with (-Δ) and without genomic deletion. A- Δ: A-III-3, and A-III-5 isolates; D- Δ: D-II-2 to D-II-5 isolates; F- Δ: F-II-1 to F-II-5 isolates. Significance was identified as *P* < 0.05 (*) or *P* < 0.01 (**) by Kruskal-Wallis test.</figcaption></figure>

area was previously described as subject to genomic deletion. In a study focused on six pyomelanin-producing mutants of *P. aeruginosa* (Hocquet et al., 2016), 46 genes in the vicinity of *hmgA* (from PA1974 to PA2019) were deleted in all the pyomelanin-producing isolates. Similarly, Wardell et al. (2019) identified large genomic deletions (from 225 kb up to 480 kb) that led to the loss of the PA2022 to PA2208 genes in 13 meropenem-resistant mutants. Thus, the deletion of genes in the vicinity of *hmgA* seems to provide an advantage for *P. aeruginosa* adaptation to human host. Gene function analysis identified that convergent evolution by deletion mostly affected genes involved in carbon compound catabolism, membrane proteins, and transcriptional regulators. Reduction in the catabolic repertoire was previously described in clinical CF isolates over time (La Rosa et al., 2019). In the same way, in an oxygen-limiting environment such as urinary tract, *P. aeruginosa* could divert its energy sources from sugar-derived carbon sources to denitrification or arginine fermentation (Tielen et al., 2013).

We then studied whether the large genomic deletions observed for some isolates had a phenotypic impact. We first studied biofilm formation, since chronic colonization is known to be associated with a switch from planktonic to sessile lifestyle and biofilm formation (Olivares et al., 2019). In two patients (A, D); who didn't undergo urine catheterization, all isolates - whether they had the genomic deletion or not - were low-biofilm producers. These results were correlated with the observation made by Bartell et al. (2019) where biofilm formation didn't significantly correlate with duration of colonization. Isolates with genomic deletion from the third patient (F), who had intermittent self-catheterization, produced significantly less biofilm than isolates without deletion. This phenotypic observation was probably not explained by the function of the deleted genes as none were involved in biofilm formation. However, the deleted isolates also presented SNPs in genes involved in alginate biosynthesis and regulation (*algG*, *algB*, and *mucA*), which could be involved in the reduction of biofilm production.

Next, we analyzed the fitness of our isolates in rich and artificial urine medium as several pathoadaptive genes (involved in central metabolism or antimicrobial resistance) have been shown to optimize pathogen fitness during CF infections (Marvig et al., 2015a). We showed that isolates with genomic deletion presented a lower fitness than isolates without deletion, particularly in AUM. These results are similar to those of Rodríguez-Rojas et al. (2009), who identified that isolates with genomic deletion presented a lower fitness (during a 48 h of co-infection with the wild-type strain in an acute respiratory model) that contributed to a decreased clearance and an increased persistence of *P. aeruginosa* in chronic lung infections. Further analyses based on continuous cultures are needed to get closer to the urinary dynamic state.

Numerous studies have showed that *P. aeruginosa* adaptation to pulmonary airway during chronic infections was associated with mutations in genes involved in antimicrobial, like those of the MexAB-OprM efflux pump (*mexA*, *mexB*, *nalD*) or genes involved in fluoroquinolones resistance (*gyrA*, *gyrB*) (Marvig et al., 2015b; López-Causapé et al., 2018; Bartell et al., 2019; Colque et al., 2020). In contrast, in our study, isolates from only one patient (D) harbored late stage SNP in a gene involved in antimicrobial resistance. This mutation in *mexS*, involved in regulation of the MexEF-OprN efflux pump, was correlated with increased within-sample diversity in terms of antimicrobial resistance profiles. In parallel, genomic deletions were associated for two patients (A, F) with the loss of the *mexX*, *mexY*, *mexZ* genes (involved in the MexXY-OprM efflux pump and aminoglycosides resistance) (López-Causapé et al., 2018), but they were correlated with phenotypic aminoglycoside-resistance for only one patient. Thus, acquisition of antimicrobial resistance doesn't appear to provide a main advantage for *P. aeruginosa* adaptation in the urinary tract, even if a longer follow-up of the patients would have been useful to confirm this hypothesis. However, further genotypic analyses by real-time reverse transcription-PCR would be needed to determine the level of expression of genes potentially involved in antimicrobial resistance.
This study presented some limitations. First, a higher number of isolates per urine sample may be of interest to explore intra-host diversity and limit the risk of identifying late stage SNPs that may have been present in isolates from the first sample not analyzed in this study. Second, using two separate methodologies for genome sequencing may have biased the identification of the SNPs. Nevertheless, this approach allowed the mapping of the genomes on an intra-patient genome reference that was sequenced by a long-read technology. Furthermore, it could be interesting to explore other genomic adaptation mechanisms, such as short insertions-deletions. Third, genomic analyses required the use of numerous successive bioinformatics tools, each with its own sensitivity and specificity.

Finally, all patients included in our study were infected or colonized by a given clone type over time, even though diversification could occur. *P. aeruginosa* adaptation to the urinary tract frequently occurred by SNPs and/or deletion in genes encoding transcriptional regulators, two-component systems, and carbon compound catabolism. Limitation of oxygen and nutrients in the urinary tract could drive *P. aeruginosa* to adapt its metabolism in favor of non-carbon energy sources and then to remodel its regulatory networks. Further metabolomic and phenotypic studies (such as biofilm formation in artificial urine medium) would be of interest to identify the regulatory pathways used by *P. aeruginosa* for adaptation to urinary conditions. On the other hand, metagenomic analyses could be interesting to explore the dynamics and interplay of the urinary microbiome with *P. aeruginosa* persistent bacteriuria.

DATA AVAILABILITY STATEMENT

The datasets presented in this study can be found in
online repositories. The names of the repository/repositories

and accession number(s) can be found below: PRJNA6
56414.

ETHICS STATEMENT

According to the French regulation on observational database analyses, this study did not need specific informed consent requirements; it was registered by the Directorate of Clinical Research and Innovation of the Rouen University Hospital Centre under the number 2018/413/OB.

AUTHOR CONTRIBUTIONS

AC, SD, MP-C, and EJ-B: conception and design. AC, MG, and TF: phenotypic experiments. AC and ML: genomic analyses. AC, ML, MG, FG, FA, and SL: interpretation of data. AC, ML, SD, and MP-C: manuscript preparation. All authors amended and approved the final version of the manuscript.

ACKNOWLEDGMENTS

The authors would like to thank the “Centre Régional Informatique et d’Applications Numériques de Normandie” for high-performance computing resources dedicated to our researches.

SUPPLEMENTARY MATERIAL

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

REFERENCES

Bartell, J. A., Sommer, L. M., Haagensen, J. A. J., Loch, A., Espinosa, R., Molin, S., et al. (2019). Evolutionary highways to persistent bacterial infection. *Nat. Commun.* 10:629. doi: 10.1038/s41467-019-08504-7

Bianconi, I., D'Arcangelo, S., Esposito, A., Benedet, M., Piffer, E., Dinnella, G., et al. (2018). Persistence and microevolution of *Pseudomonas aeruginosa* in the cystic fibrosis lung: a single-patient longitudinal genomic study. *Front. Microbiol.* 9:3242. doi: 10.3389/fmicb.2018.03242

Brooks, T., and Keevil, C. W. (1997). A simple artificial urine for the growth of urinary pathogens. *Lett. Appl. Microbiol.* 24, 203–206. doi: 10.1046/j.1472-765X. 1997.00378.x

Caballero, J. D., Clark, S. T., Coburn, B., Zhang, Y., Wang, P. W., Donaldson, S. L., et al. (2015). Selective sweeps and parallel pathoadaptation drive *Pseudomonas aeruginosa* evolution in the cystic fibrosis lung. *mBio* 6:e00981-15. doi: 10.1128/mBio.00981-15

Cabot, G., Zamorano, L., Moyà, B., Juan, C., Navas, A., Blázquez, J., et al. (2016). Evolution of *Pseudomonas aeruginosa* antimicrobial resistance and fitness under low and high mutation rates. *Antimicrob. Agents Chemother.* 60, 1767–1778. doi: 10.1128/AAC.02676-15

Cingolani, P., Platts, A., Wang, L. L., Coon, M., Nguyen, T., Wang, L., et al. (2012). A program for annotating and predicting the effects of single

nucleotide polymorphisms, SnpEff. *Fly (Austin)* 6, 80–92. doi: 10.4161/fly.1
9695

Clausen, P. T. L. C., Aarestrup, F. M., and Lund, O. (2018). Rapid and precise alignment of raw reads against redundant databases with KMA. *BMC Bioinformatics* 19:307. doi: 10.1186/s12859-018-2336-6

Colque, C. A., Albarracin Orio, A. G., Feliziani, S., Marvig, R. L., Tobares, A. R., Johansen, H. K., et al. (2020). Hypermutator *Pseudomonas aeruginosa* exploits multiple genetic pathways to develop multidrug resistance during long-term infections in the airways of cystic fibrosis patients. *Antimicrob. Agents Chemother.* 64:e02142-19. doi: 10.1128/AAC.02142-19

Curran, B., Jonas, D., Grundmann, H., Pitt, T., and Dowson, C. G. (2004). Development of a multilocus sequence typing scheme for the opportunistic pathogen *Pseudomonas aeruginosa*. *J. Clin. Microbiol.* 42, 5644–5649. doi: 10. 1128/JCM.42.12.5644-5649.2004

Djordjevic, Z., Folic, M. M., Zivic, Z., Markovic, V., and Jankovic, S. M. (2013). Nosocomial urinary tract infections caused by *Pseudomonas aeruginosa* and *Acinetobacter* species: sensitivity to antibiotics and risk factors. *Am. J. Infect. Control* 41, 1182–1187. doi: 10.1016/j.ajic.2013.02.018

Gomila, A., Carratalà, J., Eliakim-Raz, N., Shaw, E., Wiegand, I., Vallejo-Torres, L., et al. (2018). 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. doi: 10.2147/IDR. S185753
Hocquet, D., Petitjean, M., Rohmer, L., Valot, B., Kulasekara, H. D., Bedel, E., et al. (2016). Pyomelanin-producing *Pseudomonas aeruginosa* selected during chronic infections have a large chromosomal deletion which confers resistance to pyocins. *Environ. Microbiol.* 18, 3482–3493. doi: 10.1111/1462-2920.13336

Hussain, M., Suliman, M., Ahmed, A., Altayb, H., and Elneima, E. (2017). Draft genome sequence of a multidrug-resistant *Pseudomonas aeruginosa* strain isolated from a patient with a urinary tract infection in Khartoum, Sudan. *Genome Announc.* 5:e00203-17. doi: 10.1128/genomeA.00203-17

Klockgether, J., and Tümler, B. (2017). Recent advances in understanding *Pseudomonas aeruginosa* as a pathogen. *F1000Res.* 6:1261. doi: 10.12688/f1000research.10506.1

Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. *Mol. Biol. Evol.* 35, 1547–1549. doi: 10.1093/molbev/msy096

La Rosa, R., Johansen, H. K., and Molin, S. (2019). Adapting to the airways: metabolic requirements of *Pseudomonas aeruginosa* during the infection of cystic fibrosis patients. *Metabolites* 9:234. doi: 10.3390/metabo9100234

Lamas Ferreiro, J. L., Álvarez Otero, J., González González, L., Novoa Lamazares, L., Arca Blanco, A., Bermúdez Sanjurjo, J. R., et al. (2017). *Pseudomonas aeruginosa* urinary tract infections in hospitalized patients: mortality and prognostic factors. *PLoS One* 12:e0178178. doi: 10.1371/journal.pone.0178178

Larsen, M. V., Cosentino, S., Rasmussen, S., Friis, C., Hasman, H., Marvig, R. L., et al. (2012). Multilocus sequence typing of total-genome-sequenced bacteria. *J. Clin. Microbiol.* 50, 1355–1361. doi: 10.1128/JCM.06094-11

Le, S., Yao, X., Lu, S., Tan, Y., Rao, X., Li, M., et al. (2014). Chromosomal DNA deletion confers phage resistance to *Pseudomonas aeruginosa*. *Sci. Rep.* 4:4738. doi: 10.1038/srep04738

Letunic, I., and Bork, P. (2019). Interactive Tree Of Life (iTOL) v4: recent updates and new developments. *Nucleic Acids Res.* 47, W256–W259. doi: 10.1093/nar/gkz239

López-Causapé, C., Cabot, G., Del Barrio-Tofiño, E., and Oliver, A. (2018). The versatile mutational resistome of *Pseudomonas aeruginosa*. *Front. Microbiol.* 9:685. doi: 10.3389/fmicb.2018.00685

Magiorakos, A.-P., Srinivasan, A., Carey, R. B., Carmeli, Y., Falagas, M. E., Giske, C. G., et al. (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. doi: 10.1111/j.1469-0691.2011.03570.x

Marvig, R. L., Sommer, L. M., Jelsbak, L., Molin, S., and Johansen, H. K. (2015a). Evolutionary insight from whole-genome sequencing of *Pseudomonas aeruginosa* from cystic fibrosis patients. *Future Microbiol.* 10, 599–611. doi: 10.2217/fmb.15.3

Marvig, R. L., Sommer, L. M., Molin, S., and Johansen, H. K. (2015b). Convergent evolution and adaptation of *Pseudomonas aeruginosa* within patients with cystic fibrosis. *Nat. Genet.* 47, 57–64. doi: 10.1038/ng.3148

Olivares, E., Badel-Berchoux, S., Provot, C., Prévost, G., Bernardi, T., and Jehl, F. (2019). Clinical impact of antibiotics for the treatment of *Pseudomonas aeruginosa* biofilm infections. *Front. Microbiol.* 10:2894. doi: 10.3389/fmicb.2019.02894

O'Toole, G. A. (2011). Microtiter dish biofilm formation assay. *J. Vis. Exp.* 2437. doi: 10.3791/2437

Pestrak, M. J., Chaney, S. B., Eggleston, H. C., Dellos-Nolan, S., Dixit, S., Mathew-Steiner, S. S., et al. (2018). *Pseudomonas aeruginosa* rugose small-colony variants evade host clearance, are hyper-inflammatory, and persist in multiple host environments. *PLoS Pathog.* 14:e1006842. doi: 10.1371/journal.ppat.1006842

Qin, X., Zhou, C., Zerr, D. M., Adler, A., Addetia, A., Yuan, S., et al. (2018). Heterogeneous antimicrobial susceptibility characteristics in *Pseudomonas aeruginosa* isolates from cystic fibrosis patients. *mSphere* 3:e00615-17. doi: 10.1128/mSphere.00615-17

Rau, M. H., Marvig, R. L., Ehrlich, G. D., Molin, S., and Jelsbak, L. (2012). Deletion and acquisition of genomic content during early stage adaptation of *Pseudomonas aeruginosa* to a human host environment. *Environ. Microbiol.* 14, 2200–2211. doi: 10.1111/j.1462-2920.2012.02795.x

Rodríguez-Rojas, A., Mena, A., Martín, S., Borrell, N., Oliver, A., and Blázquez, J. (2009). Inactivation of the hmgA gene of *Pseudomonas aeruginosa* leads to pyomelanin hyperproduction, stress resistance and increased persistence in chronic lung infection. *Microbiol. Read. Engl.* 155, 1050–1057. doi: 10.1099/mic.0.024745-0

Soares, A., Alexandre, K., and Etienne, M. (2020). Tolerance and persistence of *Pseudomonas aeruginosa* in biofilms exposed to antibiotics: molecular mechanisms, antibiotic strategies and therapeutic perspectives. *Front. Microbiol.* 11:2057. doi: 10.3389/fmicb.2020.02057

Stothard, P., and Wishart, D. S. (2005). Circular genome visualization and exploration using CGView. *Bioinformatics* 21, 537–539. doi: 10.1093/bioinformatics/bti054

Tatusova, T., DiCuccio, M., Badretdin, A., Chetvernin, V., Nawrocki, E. P., Zaslavsky, L., et al. (2016). NCBI prokaryotic genome annotation pipeline. *Nucleic Acids Res.* 44, 6614–6624. doi: 10.1093/nar/gkw569

Tielen, P., Rosin, N., Meyer, A.-K., Dohnt, K., Haddad, I., Jänsch, L., et al. (2013). Regulatory and metabolic networks for the adaptation of *Pseudomonas aeruginosa* biofilms to urinary tract-like conditions. *PLoS One* 8:e71845. doi: 10.1371/journal.pone.0071845

Tielen, P., Wibberg, D., Blom, J., Rosin, N., Meyer, A.-K., Bunk, B., et al. (2014). Genome sequence of the small-colony variant *Pseudomonas aeruginosa* MH27, isolated from a chronic urethral catheter infection. *Genome Announc* 2:e01174-13. doi: 10.1128/genomeA.01174-13

Treangen, T. J., Ondov, B. D., Koren, S., and Phillippy, A. M. (2014). The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. *Genome Biol.* 15:524. doi: 10.1186/s13059-014-0524-x

Vorhölter, F.-J., Tielen, P., Wibberg, D., Narten, M., Schobert, M., Tüpker, R., et al. (2015). Genome sequence of the urethral catheter isolate *Pseudomonas aeruginosa* MH19. *Genome Announc* 3:e115–e115. doi: 10.1128/genomeA.00115-15

Wardell, S. J. T., Rehman, A., Martin, L. W., Winstanley, C., Patrick, W. M., and Lamont, I. L. (2019). A large-scale whole-genome comparison shows that experimental evolution in response to antibiotics predicts changes in naturally evolved clinical *Pseudomonas aeruginosa*. *Antimicrob. Agents Chemother.* 63:e01619-19. doi: 10.1128/AAC.01619-19

Wibberg, D., Tielen, P., Blom, J., Rosin, N., Schobert, M., Tüpker, R., et al. (2014). Genome sequence of the acute urethral catheter isolate *Pseudomonas aeruginosa* MH38. *Genome Announc* 2:e00161-14. doi: 10.1128/genomeA.00161-14

Wibberg, D., Tielen, P., Narten, M., Schobert, M., Blom, J., Schatschneider, S., et al. (2015). Genome sequence of the urethral isolate *Pseudomonas aeruginosa* RN21. *Genome Announc* 3:e00788-15. doi: 10.1128/genomeA.00788-15

Winstanley, C., O'Brien, S., and Brockhurst, M. A. (2016). *Pseudomonas aeruginosa* evolutionary adaptation and diversification in cystic fibrosis chronic lung infections. *Trends Microbiol.* 24, 327–337. doi: 10.1016/j.tim.2016.01.008

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Cottalorda, Leoz, Dahyot, Gravey, Grand, Froidure, Aujoulat, Le Hello, Jumas-Bilak and Pestel-Caron. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.