# Dynamics of extended-spectrum beta-lactamase-producing Enterobacterales colonization in long-term carriers following travel abroad

Laurence Armand-Lefèvre, Emilie Rondinaud, Dimitri Desvillechabrol, Jimmy

Mullaert, Olivier Clermont, Marie Petitjean, Etienne Ruppe, Thomas

Cokelaer, Christiane Bouchier, Olivier Tenaillon, et al.

► To cite this version:

Laurence Armand-Lefèvre, Emilie Rondinaud, Dimitri Desvillechabrol, Jimmy Mullaert, Olivier Clermont, et al. Dynamics of extended-spectrum beta-lactamase-producing Enterobacterales colonization in long-term carriers following travel abroad. Microbial Genomics, 2021, 7 (7), pp.000576. 10.1099/mgen.0.000576. pasteur-03509051

HAL Id: pasteur-03509051

https://pasteur.hal.science/pasteur-03509051

Submitted on 3 Jan 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>

OPEN
DATA

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

OPEN
ACCESS

# Dynamics of extended-spectrum beta-lactamase-producing *Enterobacterales* colonization in long-term carriers following travel abroad

Laurence Armand-Lefèvre¹,²*, Emilie Rondinaud¹,², Dimitri Desvillechabrol³, Jimmy Mullaert², Olivier Clermont², Marie Petitjean², Etienne Ruppe¹,², Thomas Cokelaer³,⁴, Christiane Bouchier³, Olivier Tenaillon², Laurence Ma³, Yasmine Nooroya², Sophie Matheron²,⁵, The VOYAG-R Study Group, Antoine Andremont¹,², Erick Denamur²,⁶ and Sean P. Kennedy⁷

## Abstract

Travel to tropical regions is associated with high risk of acquiring extended-spectrum beta-lactamase-producing *Enterobacterales* (ESBL-E) that are typically cleared in less than 3 months following return. The conditions leading to persistent carriage that exceeds 3 months in some travellers require investigation. Whole-genome sequencing (Illu-mina MiSeq) was performed on the 82 ESBL-E isolates detected upon return and 1, 2, 3, 6 and 12 months later from the stools of 11 long-term (>3 months) ESBL-E carriers following travel abroad. One to five different ESBL *Escherichia coli* strains were detected per traveller upon return, and this diminished to one after 3 months. Long-term carriage was due to the presence of the same ESBL *E. coli* strain, for more than 3 months, in 9 out of 11 travellers, belonging to epidemic sequence type complexes (STc 10, 14, 38, 69, 131 and 648). The mean carriage duration of strains belonging to phylogroups B2/D/F, associated with extra-intestinal virulence, was higher than that for commensal-associated A/B1/E phylogroups (3.5 vs 0.5 months, P=0.021). Genes encoding iron capture systems (*fyuA*, *irp*), toxins (*senB*, *sat*), adhesins (*flu*, *daaF*, *afa/nfaE*, *pap*, *ecpA*) and colicin (*cjrA*) were more often present in persistent strains than in transient ones. Single-nucleotide polymorphism (SNP) analysis in persistent strains showed a maximum divergence of eight SNPs over 12 months without signs of adaptation. Genomic plasticity was observed during the follow-up with the loss or gain of mobile genetic elements such as plasmids, integrons and/or transposons that may contain resistance genes at different points in the follow-up. Long-term colonization of ESBL-E following travel is primarily due to the acquisition of *E. coli* strains belonging to epidemic clones and harbouring 'virulence genes', allowing good adaptation to the intestinal microbiota.

## DATA SUMMARY

Raw sequence data and assemblies are available
in the European Nucleotide Archive (EMBL-EMI)

(http://www.ebi.ac.uk/ena) under project accession
PRJEB41147 (sample accession numbers ERS5275910
to ERS5275994). The full list and characteristics of
these strains are presented in Table S1 (available in the

Received 11 December 2020; Accepted 01 April 2021; Published 19 July 2021

**Author affiliations:** ¹Laboratoire de Bactériologie, Hôpital Bichat-Claude Bernard, AP-HP Nord-Université de Paris, F-75018 Paris, France; ²Université de Paris, IAME, INSERM UMR 1137, F-75018 Paris, France; ³Plate-forme Technologique Biomics - Centre de Ressources et Recherches Technologiques (C2RT), Institut Pasteur, F-75015 Paris, France; ⁴Hub de Bioinformatique et Biostatistique - Département Biologie Computationnelle, Institut Pasteur, USR 3756 CNRS, F-75015 Paris, France; ⁵Service de Maladies Infectieuses et Tropicales, Hôpital Bichat-Claude Bernard, AP-HP Nord-Université de Paris, F-75018 Paris, France; ⁶Laboratoire de Génétique Moléculaire, Hôpital Bichat-Claude Bernard, AP-HP Nord-Université de Paris, F-75018 Paris, France. ⁷Département Biologie Computationnelle, Institut Pasteur, USR 3756 CNRS, F-75015 Paris, France.

*Correspondence: Laurence Armand-Lefèvre, laurence.armand@aphp.fr

Keywords: *E. coli*; ESBL; carriage; long-term carriage; persistence; whole-genome sequencing; travel.

**Abbreviations:** BHI, brain heart infusion; CGA, clonal group A; ESBL, extended-spectrum beta-lactamase; ESBL-E, extended-spectrum beta-lactamase-producing *Enterobacteriales*; HPI, high pathogenicity island; IS, insertion sequence; ND, not determined; pMLST, plasmid multilocus sequence typing; SNP, single-nucleotide polymorphism; ST, sequence type; STc, sequence types complexe; VBNC, viable but non culturable; WGS, whole-genome sequencing.

**Data statement:** All supporting data, code and protocols have been provided within the article or through supplementary data files. Three supplementary tables are available with the online version of this article.

000576 © 2021 The Authors

<figure><img src="image_4.png" /></figure>
online version of this article) with the genome accession numbers.

## INTRODUCTION

The worldwide dissemination of multi-drug-resistant *Enterobacterales*, especially those producing extended-spectrum beta-lactamase (ESBL-E), is alarming. For nearly two decades, community carriage of ESBL-E has been steadily increasing [1]. However, this increase has been uneven over different regions of the world. The proportion of people carrying ESBL-E in their gut is estimated at 70% in Asia and 45% in Africa, but only 10% in Europe [1]. It is recognized that travel from low-prevalence to endemic areas increases the risk of acquiring ESBL-E, with acquisition rates of ESBL-E varying from 20–50%, depending on the geographical areas visited [2]. However, ESBL-E appear to be cleared quickly, as the median duration of carriage upon return is approximately 1 to 3 months [2–4]. Nevertheless, 5–10% of travellers who acquire an ESBL-E remain carriers 6 or 12 months after return [2, 3]. The analysis of such long-term carriage is of major interest, as it certainly contributes to the enhanced spread of ESBL-E in communities of low-prevalence countries [3]. Previous studies have shown that longer carriage of travel-acquired ESBL-E is associated with several factors, including microbiota composition, high intestinal concentration of ESBL-E, Asian destination, previous antibiotic therapy, presence of *bla*<sub>CTX-M9</sub> and colonization by *Escherichia coli* [2, 3, 5]. Indeed, among *Enterobacterales*, *E. coli* is the species best adapted species to the human intestinal microbiota [6]. Several studies, focusing rather on susceptible strains, have identified genetic factors, such as haemolysin, iron sensors and adhesins/fimbriae associated with long-term *E. coli* carriage [7, 8]. Persistent strains are also more likely to belong to the extra-intestinal virulent B2 phylogenetic group [7, 9]. However, these factors were inconsistently found in the various studies and have not been explored in multi-drug-resistant strain carriage. In a previous study aiming to assess the acquisition rates of multi-drug-resistant *Enterobacterales* in 574 subjects travelling in tropical areas, we showed that 51% of travellers had acquired at least 1 ESBL-E and that the median duration of carriage after return was 1 month. However, we found that 11 of these travellers were still carriers 6 months after return [2]. Here, we describe the dynamics of ESBL-E colonization in these long-term carriers and the bacterial genomic factors associated with persistence.

## METHODS

### General presentation of the study

This study is an ancillary study of the VOYAG-R clinical trial (clinicaltrials.gov number NCT01526187), which aimed to determine the acquisition rate of multi-drug-resistant bacteria after a trip to the tropics and the duration of carriage after return [2]. In 2013, 574 carriage-free travellers provided

### Impact Statement

Travel to the tropics has long been associated with a high risk of intestinal acquisition of extended-spectrum beta-lactamase-producing *Enterobacterales* (ESBL-E). Because ESBL-E are typically eliminated rapidly upon return, long-term carriage remains an underexamined area. The present study explores the conditions giving rise to long-term ESBL-E carriage following travel as well as the underlying mechanisms. Our findings suggest that long-term carriage is primarily due to the carriage of a single persistent strain. This strain, we observe, is the result of neither the transfer of a plasmid-carrying ESBL to a commensal strain nor of a major adaptation to the host. Rather, we find that long-term carriage is a consequence of the acquisition of epidemic clones of *E. coli* mainly belonging to the phylogenetic groups B2/D/F and harbouring specific genetic traits (adhesins, iron capture systems, colicin and toxins) that make them immediately well adapted and allow them to successfully colonize their host.

faecal samples upon return (M0) and in the case of ESBL-E acquisition, 1 (M1), 2 (M2), 3 (M3), 6 (M6) and 12 (M12) months later. Stool samples were plated directly on ChromID ESBL agar (bioMérieux, Marcy-l'Etoile, France) and bi-plate ESBL agar (AES Chemunex, Ivry-sur-Seine, France). Concurrently, an enrichment step with a brain heart infusion (BHI) broth supplemented with 1.5 mg/l<sup>-1</sup> cefotaxime and incubated overnight was performed and plated on ChromID ESBL agar. All colonies with distinct aspects (size, colour, shape) on each medium and at each time point were identified by mass spectrometry (MALDI biotyper, Bruker-Daltonics, Bremen, Germany) and tested for antibiotic susceptibility (amoxicillin, co-amoxiclav, ticarcillin, cefotaxime, ceftazidime, cefepime, cefoxitin, ertapenem, imipenem, gentamicin, amikacin, nalidixic acid, ofloxacin, co-trimoxazole, tetracycline, fosfomycin) using the disc diffusion method (http://www.eucast.org). Isolates with distinct susceptibility patterns were stored in glycerol at -80 °C [2].

### Definitions

A long term-carrier was defined as a traveller carrying an ESBL-E for more than 3 months after return.

A strain was considered persistent if found in the same traveller for more than 3 months. Transient strains were those observed only at return and never afterwards (<1 month).

Except for strains isolated at M0, the duration of carriage of the strains detected at a single follow-up point could not be calculated and was assigned the status 'not determined' (ND).

### Strain selection and characterization

All the ESBL-E isolates identified in the long-term carriers at each time point were studied. DNA was extracted using
the EZ1 DNA tissue kit (Qiagen, Courtaboeuf, France). The phylogenetic group of *E. coli* was determined as described previously [10], resulting in their classification into eight groups (A, B1, B2, C, D, E, F and *Escherichia cryptic clade I*).

## Whole-genome sequencing (WGS) and subsequent analyses

DNA libraries were prepared using the Illumina Nextera Kit (Illumina, San Diego, USA). Pooled libraries were sequenced on an Illumina MiSeq instrument. *De novo* assemblies were created with a pipeline provided within the Sequana project [11], which uses SPAdes (v3.9.0) software [12]. The quality of the *de novo* assemblies was estimated using standard metrics computed using QUAST (v4.3) [13]. Reads used to construct the assemblies were remapped against the assembly contigs to visualize coverage and detect assembly errors.

Parsnp (v1.2) [14] and PhyML (v3.3.20180621) [15] software were used to construct a maximumlikelihood phylogenetic tree. The Center for Genomic Epidemiology website was consulted to identify sequence types (STs) and ST complexes (STc), serotypes, antimicrobial resistance genes, plasmid replicon and ST [plasmid multilocus sequence typing (pMLST)] (http://www.genomicepidemiology.org). Virulence factor genes were detected using Abricate software (v0.9.8) [16]. Localization of the ESBL-encoding genes was predicted by PlaScope (v1.3.1) [17]. Strains were thus characterized using the clone definition, based on the unique combination of phylogroup, Warwick University scheme ST, Pasteur Institute scheme ST, serotype, fimH allele and ESBL enzyme.

Within the same clone, the entire genomes of the isolates were compared. First, single-nucleotide polymorphisms (SNPs) were detected on the common genome of each clone using the variant calling pipeline from the Sequana project designed according to the best practices recommended by the Broad Institute [18]. Insertions/deletions were realigned using GATK3 (v3.8) Indel Realigner software, duplicated reads were ignored, and aligned reads with a mapping quality score <30 were removed (BWA score) [19]. Variant calling was performed with Freebayes (v1.0.2) using the default options except for the ploidy option, which was set at 1. Variants were filtered, and those with a Freebayes score >100, a frequency of 80%, a minimal depth of 10, and a minimal forward and reverse strand ratio of 0.2 were retained. Mutations were confirmed using the Integrative Genomics Viewer (IGV) tool (v2.3.8). Second, a global genome comparison was performed using a k-mer approach. This approach takes into account both SNPs and the presence/absence of genes. Genetic distances between all the genomes were calculated using 'mash dist' function in Mash software (v2.27.1) with a k-mer size of 32 and a sampling size of 5000 [20]. A tree based on the distance matrix was constructed.

Annotations of the genome were performed using RAST software [21] via the PATRIC Platform (v3.0) (https://www.patricbrc.org/) [22] and using the Microbial Genome Annotation and Analysis Platform (MaGe) (v3.5) [23]. Genes functions were checked on Uniprot (https://www.uniprot.org/).

## Mutation analysis

To compute the ratio of synonymous to non-synonymous mutation rates, we counted the number of synonymous sites (S) and non-synonymous sites (N) in two diverged genomes of *E. coli*: strain 536 (B2 phylogroup) and REL606 (A phylogroup) [24]. We computed the rates as the total number of mutations observed (s for the synonymous, n for the nonsynonymous) divided by the number of possible sites, and the ratio of rates was computed as $R = \frac{n/N}{s/S}$.

## Statistical analysis

Based on the carriage status at each time point, the carriage duration was estimated as the number of months during which the strain was detected. Strains isolated on return (M0) were assigned a carriage duration of 0.5 months. As an additional sensitivity analysis, we also compared (alternative method) the carriage durations, considering that all strains isolated up to M3 were present at M0 (i.e. acquired during the journey).

Carriage durations were compared between strains belonging to the commensalism-associated A/B1/E phylogroups and the extra-intestinal virulence-associated B2/D/F phylogroups using the non-parametric Wilcoxon rank sum test. The rationale of this grouping is based on epidemiological and animal model experimental data [6, 25, 26].

The proportion of the presence/absence of virulence genes was compared between persistent and transient strains in a targeted approach. The *P*-value of the Fisher exact test was reported for all 149 genes of the panel. Due to the relatively low number of strains in each category, no multiple testing corrections were performed, and it was determined that the reported *P*-value should be used for gene ranking and not for the evaluation of statistical significance.

The same analysis was performed for the non-targeted approach, using the genes identified by the RAST and MaGe annotation tools. All genes that were not part of the core genome (genes of the variable genome) were included in the comparison analysis. In case of discrepancies, the presence of the gene was confirmed by BLAST on MaGe software.

## RESULTS

Among the 292 travellers who acquired a multi-drug-resistant *Enterobacteriales* during their trip, 11 (3.8%) were long-term carriers (6 travellers were still carriers at M6 and 5 at M12). Eight individuals reported having travelled to Asian countries, while three travelled to Latin America. The median duration of travel was 18 days [interquartile range (IQR), 15.5–22 days]. Six of them reported taking antibiotics once or twice during the follow-up. (Table S1).

## Dynamic of colonization

During the follow-up period, 82 multi-drug-resistant *Enterobacteriales* isolates, all ESBL *E. coli*, were isolated in the 11 long-term carriers (Table S1). The mean (range) number of
<figure><img src="image_5.png" /><figcaption>Fig. 1. Follow-up of ESBL E. coli carriage in travellers after return.</figcaption></figure>

isolates detected per traveller was 2 (1–5) at M0, 1.5 (1–3) at M1, 1.3 (1–2) at M2, 1.1 (1–2) at M3 and 1 at M6 and M12.

The analysis of genomes of the 82 isolates identified 35 genetically different strains (clones) (Fig. 1, Table S1). The duration of carriage of the strains ranged from 0.5 to 12 months (median 1 month). Eleven strains were considered transient (duration of carriage: 0.5 months) and nine persistent (>3 months). The remaining had, for 10, an estimated duration of carriage between 1 and 3 months and for 5, a duration of carriage undetermined because isolated only at a single point of the follow-up (Fig. 1).

In 9 out of 11 travellers, long-term carriage was due to the presence of the same strain for more than 3 months (persistent strain). Indeed, in two travellers (T1 and T2), a single strain was isolated upon return and persisted for 12 months. In four (T3 to T6), several strains were isolated upon return but only one persisted. In three (T7 to T9), several strains were isolated upon return, but the persistent strain was isolated at M1. Last, in two travellers (T10 and T11) the long-term carriage was due to the acquisition at M6 of a new ESBL-E that had not previously been detected, and no persistent strains were observed.

## Global characteristics of strains

The 35 strains were distributed in 16 STc; the most frequent was STc10 (n=10), followed by STc131 and STc648 (n=4), STc38 (n=3), and STc 165 and 467 (n=2), with the others being singletons.

The genes encoding ESBL enzymes were predominantly bla<sub>CTX-M</sub> genes (n=33, 94%) and were mostly predicted to be plasmidic (n=25, 71%). In all, 66 plasmids were detected

in the 35 strains, and these were primarily IncF type (49%, n=32), followed by IncI type (12%, n=9) (Fig. 1).

Fifteen (43%) of these 35 different strains belonged to B2/D/F extra-intestinal virulent phylogroups, the median carriage duration of which was higher than that of those belonging to commensal A/B1/E phylogroups (3.5 vs 0.5 months, P=0.021). The significance of the difference remained when the alternative method for the estimation of carriage (i.e. from the return) duration was used (4.5 vs 0.5, P=0.018) (Fig. 2).

## Transient and persistent strain comparison

In order to investigate genetic traits potentially associated with persistence, the 9 persistent strains were compared to the 11 transient strains. Due to a risk of misclassification of the other strains resulting from intermittent carriage or lack of sensitivity of the culture, they were excluded from the comparison analysis.

## Phylogenetic background

The persistent strains belonged primarily (n=7/9, 78%) to extra-intestinal virulence-associated phylogenetic groups B2 (n=3), D (n=2) and F (n=2), whereas the transient strains fell mostly (8/11, 72.7%) within commensal phylogenetic groups A (n=6), B1 (n=1) and E (n=1). Furthermore, the persistent strains mainly belonged to extra-intestinal pathogenic E. coli (ExPEC) epidemic clonal complexes, often reported as ESBL-E, including STc131 (n=2), STc648 (n=2), STc14/ST1193, STc69 [clonal group A (CGA)] and STc38. Transient strains were found distributed across 10 different STc (Fig. 1).
<figure><img src="image_6.png" /><figcaption>**Fig. 2.** Estimated duration of carriage depending on the phylogenetic group of the ESBL *E. coli* strains. The phylogenetic tree was generated using the maximum-likelihood method. The tree was rooted on *Escherichia fergusonii*. The bootstrap values (1000 replicates) are shown near the nodes. Scale bar represents the number of nucleotide substitutions per site.</figcaption></figure>

**Targeted approach: plasmid, antibiotic resistance and virulence genes**

No differences were found between the 12 types of plasmids isolated in the 9 persistent and the 11 transient strains. Of 42 different detected resistance genes, only the *qnrS* gene, coding for the quinolone resistance, was found more frequently in transient compared to persistent strains (Table 1).

Of 149 genes known to be involved in virulence, 7 genes or operons/gene clusters were found more frequently in persistent strains than in transient ones. Four of them are involved in adhesion: the *papG* and *papC* genes belonging to the *pap* operon encoding P fimbriae, *nfaE* encoding a non-fimbrial adhesin, *daaF* encoding a diffuse adherence adhesin and the *ecpA*, *B*, *C*, *D*, *E* genes encoding *E. coli* common pilus. Two genes are involved in iron capture, the genes *fuyA* and *irp2* encoding the yersiniabactin siderophore and its transporter and belonging to the high pathogenicity island (HPI). The two remaining genes *sat* and *senB* encode toxins (Table 1).

**Non-targeted approach**

We also employed a non-targeted approach by searching for associations between any gene and the characteristics of the strains, persistent or transient. The genomes of 20 strains (9 persistent and 11 transient) were annotated. The RAST tool found a total of 6724 annotated genes, of which 3020 were part of the variable genome and were included in the comparison analysis; the MaGe tool found a total of 4535 annotated genes, of which 1835 (variable genome)

were included. Hypothetical or poorly characterized genes, and genes from phage genomes, were then excluded from further analysis.

The genes that were found to have the greatest difference in the two types of strains are listed in Table 1. The untar- geted method retrieved all of the genes previously identi- fied with the targeted method (resistance and virulence genes). The *nfaE* gene was identified as *afaB* gene (97% homology) encoding for an afimbrial adhesin. This method also identified additional genes that were more frequent in persistent strains, such as the *flu* gene encoding the antigen 43 involved in adhesion, type VI secretion system genes, an unnamed operon involved in iron transport, the *cjrA*/*cjrB* encoding a colicin and four genes encoding metabolic functions (Table 1).

**Evolution of persistent strains during follow-up**

**SNP-based comparison**

Compared to the first isolate (M0 or M1), from zero to a maximum of eight SNPs were observed in persistent strains at each time point of the follow-up. Of the 52 genetic events observed, 48 (92%) were in coding regions (11 frameshifts and 37 mutations). These genes were affected with very few convergences, as only the *htrE* gene, encoding an outer-membrane usher protein, was mutated in two different persistent strains. Thirteen (27%) genetic events were detected at more than one follow-up point and thus considered to be
<table><caption>**Table 1.** List of genes for which a different occurrence was found between persistent and transient *E. coli* strains by the targeted approach (among 47 resistance genes and 149 virulence genes) and the non-targeted approach (whole-genome comparison). Genes detected by both approaches are underlined and in bold. The reported *P*-value should be used for gene ranking and not for the evaluation of statistical significance</caption><thead><tr><th>Genes</th><th>Function</th><th>Persistent strains, n=9 (%)</th><th>Transient strains, n=11 (%)</th><th>P*</th></tr></thead><tbody><tr><td><strong>Resistance gene</strong></td><td></td><td></td><td></td><td></td></tr><tr><td><strong>qnr</strong></td><td>Fluoroquinolone resistance</td><td>0</td><td>5 (45)</td><td>0.038</td></tr><tr><td><strong>Virulence genes</strong></td><td></td><td></td><td></td><td></td></tr><tr><td><em>fyuA</em>, <em>irp1</em>, <em>irp2</em>, <em>irp3</em>, <em>irp4</em>, <em>irp5</em></td><td>Yersiniabactin iron transporter gene cluster (HPI)</td><td>9 (100)</td><td>4 (36)</td><td>0.005</td></tr><tr><td><em>ydeA</em>, <em>ydfA</em></td><td>Active transporter</td><td>5 (56)</td><td>0</td><td>0.008</td></tr><tr><td><em>shiA</em></td><td>Shikimate transporter</td><td>8 (89)</td><td>3 (27)</td><td>0.010</td></tr><tr><td><em>gfcA</em>, <em>gfcB</em>, <em>gfcC</em>, <em>gfcD</em>, <em>gfcE</em>, <em>etp</em>, <em>etk</em></td><td>Capsule biosynthesis gene cluster</td><td>1 (11)</td><td>8 (73)</td><td>0.010</td></tr><tr><td><em>flu</em></td><td>Antigen 43</td><td>9 (100)</td><td>5 (45)</td><td>0.014</td></tr><tr><td><em>sgcA</em>, <em>sgcB</em>, <em>sgcC</em>, <em>sgcE</em>, <em>scgQ</em>, <em>sgcR</em>, <em>sgcX</em>, <em>yjhF</em>, <em>yjhG</em>, <em>yjhH</em>, <em>yjhL</em>, <em>yjhP</em>, <em>yjhU</em></td><td>D-xylonate dehydratase gene cluster</td><td>6 (67)</td><td>1 (9)</td><td>0.017</td></tr><tr><td><em>sat</em></td><td>Serine protease autotransporter toxin</td><td>6 (67)</td><td>1 (9)</td><td>0.017</td></tr><tr><td>Unnamed</td><td>Fe2 +ABC transporter</td><td>7 (78)</td><td>2 (18)</td><td>0.022</td></tr><tr><td><em>scsC</em>, <em>scsD</em></td><td>Suppressor for copper tolerance or thioredoxin</td><td>7 (78)</td><td>2 (18)</td><td>0.022</td></tr><tr><td><em>papG</em>, <em>papF</em>, <em>papE</em>, <em>papK</em>, <em>papI</em>, <em>papD</em>, <strong><em>papC</em></strong>, <em>papH</em>, <em>papA</em>, <em>papI</em>, <em>papB</em></td><td><em>pap</em> gene cluster</td><td>4 (44)</td><td>0</td><td>0.026</td></tr><tr><td><em>afaA</em>, <em>afaB</em>, <em>afaC</em>, <em>afaD</em>, <em>afaF</em> / <strong><em>nfaE</em></strong></td><td>Afimbrial adhesin</td><td>4 (44)</td><td>0</td><td>0.026</td></tr><tr><td><em>daaF</em></td><td>F1845 diffuse adherence adhesin DaaF</td><td>4 (44)</td><td>0</td><td>0.026</td></tr><tr><td><em>senB</em></td><td>Enterotoxin TieB protein</td><td>4 (44)</td><td>0</td><td>0.026</td></tr><tr><td><em>cjrA</em>, <em>cjrB</em>, <em>crjC</em></td><td>Colicin Ia</td><td>4 (44)</td><td>0</td><td>0.026</td></tr><tr><td><em>impG</em>/vasA, <em>impH</em>/vasB</td><td>Type VI secretion system</td><td>8 (89)</td><td>4 (36)</td><td>0.028</td></tr><tr><td><strong><em>ecpR</em></strong>, <em>ecpA</em>, <em>ecpB</em>, <em>ecpC</em>, <em>ecpD</em>, <em>ecpE</em></td><td><em>E. coli</em> common pilus operon</td><td>9 (100)</td><td>6 (55)</td><td>0.038</td></tr><tr><td><em>mtnA</em>, <em>mtnK</em></td><td>Methylthioribose phosphatase and isomerase</td><td>5 (56)</td><td>1 (9)</td><td>0.050</td></tr></tbody></table>

fixed. Among the 37 genetic mutations, 26 (70%) were non-synonymous and 11 (30%) synonymous (Table 2, Fig. 3). We computed a ratio of non-synonymous to synonymous rate of ~0.81 using all detected mutations. This value, close to 1, suggests that non-synonymous mutations accumulate at the same rate as synonymous mutations and are therefore not subject to active selection.

## K-mer-based comparison

The global genome comparison based on the k-mer approach showed a very high similarity between the different isolates of the persistent strain of four travellers (T1, T3, T5 and T9) and some slight differences between the isolates of the persistent strain of five travellers (T2, T4, T6, T7 and T8) (Fig. 4). As very few SNPs have been evidenced, these differences may be explained by the loss or gain of total or partial genomic mobile elements, prophage in T2, integron

or transposon (containing multiple resistance genes) in T4 and T7, and plasmid in T4, T6 and T8 (Table 3).

## Potential ESBL plasmid transmission in individual gut microbiota

Plasmids were characterized (ESBL gene type, replicons, pMLST type) with regard to potential plasmid transmission of ESBL genes between strains isolated from the gut of each traveller. Plasmids were characterized with regard to potential plasmid transmission of ESBL genes between strains isolated from the gut of each traveller. In the 10 travellers (T2 to T11) from whom several strains were detected, neither identical ESBL enzymes (T2 to T4 and T10) nor identical plasmids (T6 to T7 and T11) nor similar location (chromosomic/plasmidic) of the ESBL gene (T9) were observed between the different strains isolated in a single traveller (Fig. 1). This suggests that no emergence of

<footnote>*The reported *P*-value should be used for gene ranking and not for the evaluation of statistical significance.</footnote>
<table><caption>**Table 2.** Mutation and frameshift events observed in the genome of the persistent ESBL E. coli strains during the follow-up and their impact</caption><thead><tr><th>Traveller</th><th>Gène</th><th>Synonym</th><th>Function of the encoded protein</th><th>Event</th><th>Mutation type</th><th>Fixed</th><th>M0</th><th>M1</th><th>M2</th><th>M3</th><th>M6</th><th>M12</th></tr></thead><tbody><tr><td>T1</td><td>mcsS</td><td></td><td>Mechanosensitive channel</td><td>C&gt;T</td><td>Synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>T2</td><td>intergenic</td><td></td><td></td><td>Deletion AC</td><td>Intergenic</td><td>No</td><td>Ref</td><td>X</td><td></td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>mxID</td><td>yscC,hrcC, imG</td><td>Type III secretion outer-membrane pore-forming protein</td><td>C&gt;G</td><td>Non-synonymous</td><td>Yes</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>pdeC</td><td>ylaB, yjcC</td><td>Cyclic-di-GMP phosphodiesterase</td><td>A&gt;T</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td>X</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>fisH</td><td></td><td>ATP-dependent zinc metalloprotease</td><td>T&gt;A</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td>X</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>intergenic</td><td></td><td></td><td>Insertion A</td><td>Intergenic</td><td>Yes</td><td>Ref</td><td></td><td>X</td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>ccmH</td><td>yejP</td><td>Cytochrome c-type biogenesis protein</td><td>T&gt;A</td><td>Non-sense</td><td>No</td><td>Ref</td><td>X</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>uup</td><td>ycbH, ycbI</td><td>ABC transporter ATP-binding protein</td><td>T&gt;A</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td>X</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>purL</td><td></td><td>Phosphoribosylformylglycinamidine synthase</td><td>C&gt;T</td><td>Synonymous</td><td>No</td><td>Ref</td><td>X</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>kdpC</td><td>kac</td><td>Potassium-transporting ATPase C chain</td><td>A&gt;T</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td>X</td><td></td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>aroP</td><td>ansP</td><td>Aromatic amino acid transport protein</td><td>C&gt;T</td><td>Non-synonymous</td><td>Yes</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>chaA</td><td></td><td>Sodium-potassium/proton antiporter</td><td>T&gt;A</td><td>Synonymous</td><td>Yes</td><td>Ref</td><td></td><td>X</td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>acfD</td><td>yghJ</td><td>Putative lipoprotein accessory colonization factor</td><td>T&gt;A</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td>X</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>rfaQ</td><td>waaQ, rfaQ, hepIII</td><td>Lipopolysaccharide core heptosyltransferase</td><td>T&gt;A</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>flu</td><td>yzX, yefQ, agn</td><td></td><td>Deletion TTC</td><td>Frameshift</td><td>No</td><td>Ref</td><td></td><td>X</td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>gmd</td><td>yefN, yefA</td><td>GDP-mannose 4,6-dehydratase</td><td>A&gt;C</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>T3</td><td>focC</td><td></td><td>Chaperone FimC</td><td>Insertion T</td><td>Frameshift</td><td>Yes</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>aidA</td><td>yfaL</td><td>Adhesive DNA-3-methyladenine glycosylase</td><td>C&gt;T</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>queG</td><td></td><td>Epoxyqueuosine reductase</td><td>C&gt;T</td><td>Non-sense</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td><td></td></tr><tr><td></td><td>fisH</td><td></td><td>Cell division protein</td><td>C&gt;A</td><td>Non-synonymous</td><td>Yes</td><td>Ref</td><td>X</td><td>X</td><td></td><td></td><td>X</td><td>X</td></tr><tr><td></td><td>unknown</td><td></td><td>Hypothetical protein</td><td>Deletion G</td><td>Frameshift</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td><td>X</td></tr></tbody></table>

Continued
<table><caption>**Table 2. Continued**</caption><thead><tr><th>Traveller</th><th>Gene</th><th>Synonym</th><th>Function of the encoded protein</th><th>Event</th><th>Mutation type</th><th>Fixed</th><th>M0</th><th>M1</th><th>M2</th><th>M3</th><th>M6</th><th>M12</th></tr></thead><tbody><tr><td rowspan="7">T4</td><td>hsdS</td><td>hss</td><td>Type I restriction-modification system, subunit S</td><td>Insertion AG</td><td>Frameshift</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>htrE</td><td>yehB</td><td>Outer-membrane usher protein</td><td>Insertion G</td><td>Frameshift</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>oxc</td><td>yfdU</td><td>Oxalyl-CoA decarboxylase</td><td>C>T</td><td>Synonymous</td><td>No</td><td>Ref</td><td></td><td>X</td><td></td><td></td><td></td></tr><tr><td>intergenic</td><td></td><td></td><td></td><td>C>T</td><td>Intergenic</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td>X</td></tr><tr><td>pedC</td><td>ycjG</td><td>L-alanine-DL-glutamate epimerase</td><td>C>T</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>apxIB</td><td>lapB</td><td>Type I secretion system ATPase</td><td>Deletion G</td><td>Frameshift</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>pedC</td><td>lapB</td><td>Diguanylate cyclase</td><td>T>A</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td rowspan="7">T5</td><td>dmsC</td><td>yfiH</td><td>Anaerobic dimethyl sulfoxide reductase, chain C</td><td>C>G</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td>X</td><td></td></tr><tr><td>adrA</td><td>dgcC, yaiC</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>adrA</td><td>dgcC, yaiC</td><td>1,2-phenylacetyl-CoA epoxidase, subunit C</td><td>A>T</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td>X</td><td></td><td></td><td></td></tr><tr><td>gutB</td><td>ycjQ srlB</td><td>Zinc-type alcohol dehydrogenase</td><td>G>A</td><td>Synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>unknown</td><td></td><td>Hypothetical protein</td><td>A>T</td><td>Non-sense</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>nanK</td><td>yhcI</td><td>N-acetylmannosamine kinase</td><td>Deletion TT</td><td>Frameshift</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>mdfA</td><td>cmlA, cmr</td><td>Multidrug transporter</td><td>T>C</td><td>Synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td rowspan="7">T6</td><td>intE</td><td></td><td>Prophage integrase</td><td>T>A</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>slt</td><td></td><td>Soluble lytic murein transglycosylase</td><td>C>T</td><td>Synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>intergenic</td><td></td><td></td><td></td><td>T>A</td><td>Intergenic</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td>X</td></tr><tr><td>yghW</td><td></td><td>Uncharacterized protein</td><td>C>T</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>thiC</td><td></td><td>Hydroxymethylpyrimidine phosphate synthase</td><td>A>G</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td>X</td><td></td><td></td><td>X</td></tr><tr><td>pqIC</td><td>ymbA</td><td>Polypeptide intermembrane transport lipoprotein</td><td>A>G</td><td>Non-synonymous</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr><tr><td>lacA</td><td>wbb1</td><td>Galactosyl O acetyl transferase</td><td>Insertion</td><td>Frameshift</td><td>No</td><td>Ref</td><td></td><td></td><td></td><td></td><td>X</td></tr></tbody></table>

Continued
<table><caption>**Table 2.** Continued</caption><thead><tr><th>Traveller</th><th>Gene</th><th>Synonym</th><th>Function of the encoded protein</th><th>Event</th><th>Mutation type</th><th>Fixed</th><th>M0</th><th>M1</th><th>M2</th><th>M3</th><th>M6</th><th>M12</th></tr></thead><tbody><tr><td rowspan="7">T8</td><td>bcr</td><td>bicA, bicR, suxA</td><td>Multidrug pump efflux</td><td>C>T</td><td>Non-synonymous</td><td>Yes</td><td></td><td>Ref</td><td></td><td></td><td>X</td><td>X</td></tr><tr><td>mutS</td><td>ant, pln, fdv</td><td>DNA mismatch repair protein</td><td>C>T</td><td>Non-synonymous</td><td>Yes</td><td></td><td>Ref</td><td></td><td></td><td>X</td><td>X</td></tr><tr><td>htrE</td><td>yehB</td><td>Outer-membrane usher protein</td><td>G>C</td><td>Non-synonymous</td><td>Yes</td><td></td><td>Ref</td><td></td><td></td><td>X</td><td>X</td></tr><tr><td>citT</td><td></td><td>Citrate/succinate antiporter</td><td>G>A</td><td>Synonymous</td><td>Yes</td><td></td><td>Ref</td><td></td><td></td><td>X</td><td>X</td></tr><tr><td>angR</td><td>irp2, mbtP</td><td>Iron acquisition yersiniabactin synthesis enzyme</td><td>T>C</td><td>Non-synonymous</td><td>Yes</td><td></td><td>Ref</td><td>X</td><td>X</td><td>X</td><td>X</td></tr><tr><td>unknown</td><td></td><td>Hypothetical protein</td><td>T>C</td><td>Synonymous</td><td>Yes</td><td></td><td>Ref</td><td>X</td><td></td><td>X</td><td>X</td></tr><tr><td>cusC</td><td>ibeB, ylcB</td><td>Copper/silver export system outer-membrane channel protein</td><td>G>C</td><td>Non-synonymous</td><td>No</td><td></td><td>Ref</td><td>X</td><td></td><td></td><td></td></tr><tr><td rowspan="4">T9</td><td>evgS</td><td></td><td>Hybrid sensory histidine kinase</td><td>G>A</td><td>Synonymous</td><td>No</td><td></td><td>Ref</td><td></td><td></td><td>X</td><td></td></tr><tr><td>ypjD</td><td></td><td>Uncharacterized protein</td><td>C>T</td><td>Synonymous</td><td>Yes</td><td></td><td>Ref</td><td></td><td></td><td>X</td><td>X</td></tr><tr><td>murG</td><td></td><td>N-acetylglucosaminyl transferase</td><td>C>G</td><td>Non-synonymous</td><td>Yes</td><td></td><td>Ref</td><td>X</td><td>X</td><td>X</td><td></td></tr><tr><td>gltR</td><td>hdfR, gltC</td><td>HTH-type transcriptional regulator</td><td>G>A</td><td>Non-synonymous</td><td>No</td><td></td><td>Ref</td><td>X</td><td></td><td></td><td></td></tr></tbody></table>

Hatched cases: the persistent strain has not been detected at this follow-up point.
<figure><img src="image_7.png" /><figcaption>**Fig. 3.** Mutations or frameshift events observed in the genome of persistent ESBL E. coli strains during follow-up. We created, for each traveller, a maximum-likelihood unrooted tree relating all isolates sampled during the follow-up period. For each traveller, the length of the branches is proportional to the number of mutations separating two samples. The diameters of circles are proportional to the numbers of isolates that are identical. Genes in red had non-synonymous mutation, genes in black had synonymous mutation and genes in blue had nonsense mutation or frameshift. This representation does not follow any time scale.</figcaption></figure>

ESBL-E appears to result from gene transfer from a travel-
acquired strain to a resident strain within the gut of each
carrier, but does not exclude ESBL gene transfer via other
genetic mobile elements, such as insertion sequences (ISs),
integrons, or transposons.

DISCUSSION

It is now well established that travel to the tropics is associ-
ated with a high rate of ESBL-E acquisition and that carriage
is usually short-lived. However, some travellers display
<figure><img src="image_8.png" /><figcaption>**Fig. 4.** Genetic similarity between the different isolates of the nine persistent ESBL E. coli strains. The tree was constructed from the similarity matrix generated by the k-mer comparison approach. The tree was unrooted. Scale bar represents the distance between isolates. Phylogroup and sequence type are mentioned in parentheses for each persistent strain.</figcaption></figure>

prolonged carriage [2]. Here, we showed that long-term ESBL-E carriage following travel is associated with the acquisition of epidemic clones of E. coli belonging mainly to the phylogenetic groups B2/D/F and harbouring specific genetic traits (adhesins, iron capture systems and toxins) that allow them to successfully colonize their host.

There is no consensual definition of long-term carriage, which varies from several weeks to several years depending on the studies [7, 27, 28]. Because in our original study 96% of the travellers were carriage-free 3 months following return, we considered long-term carriers to be those still carrying ESBL-E after more than 3 months. Only 11 of the 292 returning carriers met this criterion. Although multiple strains could be acquired by a single traveller, we observed that in 9 out of 11 the long-term carriage was due to colonization by a single persistent strain. Follow-up was interrupted

if a new trip to a tropical region was reported, but several subjects declared travel outside this destination. Thus, some subjects may have acquired new strains in France or Europe. Even if this risk is low, given the low prevalence of ESBL-E carriage in France (5–10%) [29, 30], this event is possible, especially as travellers represent a population at higher risk (non-French people, in contact with other people frequently travelling abroad and who may consume food products coming from high-prevalence areas).

All persistent strains were *E. coli*, the natural host of the human microbiota colonizing the microbiota at high concentrations (10⁷–10⁸) [5]. Several hypotheses might explain this prolonged carriage. One possibility is that a travel-acquired ESBL-E strain transferred its ESBL gene to a commensal *E. coli*, which was already adapted to its host. Indeed, in three travellers, the persistent strain was
<table><caption>**Table 3.** Regions of genomic plasticity observed in the genome of the persistent ESBL *E. coli* strains during the follow-up</caption><thead><tr><th>Traveller strains</th><th>Mobile elements</th><th>Resistance genes</th><th>Virulence genes</th><th>Size (bp)</th><th>M0</th><th>M1</th><th>M2</th><th>M3</th><th>M6</th><th>M12</th></tr></thead><tbody><tr><td>T1-1</td><td>None</td><td>None</td><td>None</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>T2-1</td><td>Phage</td><td>None</td><td>None</td><td>94369</td><td>P</td><td>P</td><td></td><td></td><td></td><td></td></tr><tr><td>T3-2</td><td>None</td><td>None</td><td>None</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>T4-2</td><td>Incl1 plasmid</td><td>None</td><td>None</td><td>143151</td><td></td><td>V</td><td>V</td><td>V</td><td>V</td><td>V</td></tr><tr><td>T5-2</td><td>Integron/ transposon</td><td>dfrA17, sul1, aadA5, mph(A), ermB</td><td>None</td><td>9845</td><td>P</td><td>P</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>None</td><td>None</td><td>None</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>T6-1</td><td>IncF plasmid (part)</td><td>None</td><td>None</td><td>61078</td><td></td><td>P</td><td></td><td></td><td>P</td><td></td></tr><tr><td>T7-5</td><td>Integron</td><td>dfrA17, aadA5, sul1, mph(A), tet(A), strB, strA, sul2, acc(3)-IId</td><td>None</td><td>15452</td><td></td><td></td><td>P</td><td>P</td><td>P</td><td></td></tr><tr><td>T8-5</td><td>InclL/M plasmid</td><td>None</td><td>None</td><td>53063</td><td></td><td></td><td></td><td></td><td>P</td><td>P</td></tr><tr><td>T9-3</td><td>None</td><td>None</td><td>None</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

None, absence of variation; P, presence; V, presence with variable size between isolates; the absence of indication signifies absence of the element.

Hatched cases: the persistent strain has not been detected at this follow-up point.

not isolated upon return but 1 month later. However, this hypothesis is hardly conceivable, as we found no evidence of plasmid transfer between strains isolated in a single traveller. Transfer of ESBL plasmids is known to occur between *Enterobacteriales* species within the intestinal microbiota, but the frequency of this phenomenon remains unknown [31, 32]. Further, low initial concentrations could explain why the persistent strain was not initially detected at M0 in three subjects.

A second hypothesis explaining long-term carriage is that the travel-acquired ESBL *E. coli* strain subsequently evolved and adapted to the host microbiota. Investigation along this path revealed that persistent strains in the guts of travellers accumulated a maximum of eight SNPs per genome over a 12-month period. This mutation rate is consistent with commensal *E. coli* and in line with a previous study that reported up to six SNPs in clones isolated 4 or 7 months apart [33]. Analysis of the mutated genes showed a ratio of non-synonymous to synonymous mutations of 0.81 and very few convergences or fixed mutations, which is more consistent with a genetic drift rather than with an adaptation process and thus poorly supports a significant adaptation of the acquired strain to the gut of the hosts [34]. In addition, we observed genomic plasticity in some strains in which mobile genetic elements, such as integron, transposon or plasmid, may be absent or present in isolates at different follow-up points. It was not possible to determine the exact chronology of gain/loss of these elements as we studied only one isolate per sampling. Furthermore, it cannot be excluded that the observed diversity is the result of laboratory subculture and/or storage. Metapopulation analyses

[35] coupled with *in vitro* fitness assays [36], should be performed in order to arrive at a better understanding of the evolutionary forces at play.

Thus, the most likely hypothesis is that the ESBL *E. coli* acquired during the trip already harboured the genetic factors required for successful adaptation to the host. We took advantage of the simultaneous carriage within the same microbiota of strains with different carriage durations, mimicking a natural competitive environment. We observed that ESBL *E. coli* belonging to extra-intestinal-associated B2/D/F phylogenetic groups have a longer duration of carriage than those belonging to the commensal groups A/B1/E. *E. coli* belonging to group B2 and, to a lesser extent, group D are known to be more pathogenic but also to colonize more efficiently [8, 27, 37]. In addition, persistent *E. coli* belonged to highly epidemic clonal complexes such as B2-STc131, B2-STc14/ST1193, D-STc69 (CGA), D-STc38 and F-STc648, all of them being frequently reported as ESBL producers [26, 38–40]. Genetic traits associated with persistent carriage were investigated in a targeted manner by comparing 42 resistance and 149 virulence genes, and with a non-targeted approach encompassing the entire complement of the genes identified after sequencing. Both methods confirmed that persistent strains harboured more virulence factors in their genome than transient strains. While 50–60% of the commensal *E. coli* [8, 41] carried *fyuA/irp* genes, all persistent ESBL *E. coli* had these genes irrespective of their phylogenetic group. *fyuA/irp* encode the yersiniabactin siderophore, an iron capture system belonging to the HPI, which has already been associated with longer duration of carriage, although not at such a high
frequency [8]. Genes coding for molecules involved in the adhesion process were found much more often in persistent strains. These included adhesins belonging to the Afa/Dr family (daaF and afaA/nfaE) [42], the self-recognizing adhesin, antigen 43 involved in cell aggregation (flu) [43], the P fimbriae involved in adherence in pyelonephritis (pap operon) [44] and the common E. coli pilus (ecpA) [44]. This pattern of prevalence also held for genes coding for the type VI secretion system (vasA/vasB), which could be involved in cell interaction and biofilm formation [45], and toxins such as sat encoding the secreted autotransporter toxin, a vacuolating cytotoxin [46] or senB encoding the enterotoxin TieB. As already described, the senB gene was associated with the crj operon encoding a colicin [47]. Extra-intestinal virulence genes coding for adhesins, iron capture systems and toxins have previously been correlated with successful gut colonization in humans [6]. Our results confirm that virulence factors are not only involved in strain pathogenicity, but also in commensalism [48, 49]. An important result is that the 'virulence factors' associated with persistence were found in both the widely described B2 strains as well as in the D, F and A persistent strains. These so-called virulence factors improve adaptability and competitiveness and may promote longer intestinal colonization [6, 9, 50]. It is a matter of concern when these factors are present in ESBL-E as, by facilitating longer carriage, they also contribute to a higher risk of transmission and therefore dissemination of ESBL-E. These results concerning genetic factors specific to the bacteria must be understood in a global context, integrating the role and the complex interactions between the environment, the host and its microbiota, as they are inextricably linked [5].

Our study has several limitations. First, the modest size of our study, including the number of travellers with long-term carriage, as well as the number of persistent and transient strains, limits the statistical power of the analyses. Nevertheless, we have been able to correlate our findings with similar observations in the literature. Second, we took only one morphologically different colony per plate, which is a limitation in investigating the diversity of ESBL-E carriage. This would require between 5 and 15 colonies per plate to detect subdominant clones, or a metagenomic approach [35]. Third, even though we used an enrichment step, which increases the sensitivity of the culture, strains colonizing the gut microbiota at low concentrations may have been missed. Moreover, since no molecular methods were performed directly in the stool samples, we could not exclude the possibility that the ESBL-E genes could be carried and transmitted from viable but nonculturable (VBNC) bacteria. Fourth, as no samples were taken during travel, we could not determine the exact timing of the acquisition and could not exclude the possibility that rapid strain adaptations (plasmid or mobile element transfer, mutations) may have occurred during the trip. Finally, no household members or pets have been investigated for ESBL-E colonization. Even if the risk seems low, we cannot exclude the possibility that the persistence of carriage may

be due to iterative acquisition of the same strain circulating within the household. Thus, the association between the acquisition of preadapted multi-resistant *E. coli* clones and long-term carriage after travel abroad should be confirmed in further and larger studies.

In conclusion, long-term colonization by ESBL-E following travel appears to be associated with the acquisition of *E. coli* strains belonging to epidemic clones and harbouring specific traits, including extra-intestinal virulence genes, allowing good adaptation to the intestinal microbiota.

**Funding information**

This work was partially supported by the 'Fondation pour la Recherche Médicale' (Equipe FRM 2016, grant number DEQ20161136698). The Biomics Unit is a member of the France Génomique consortium (ANR10-INBS-09-08).

**Acknowledgements**

We thank all the travellers who participated in this study. The results of this study have been reported in part at the European Congress of Clinical Microbiology and Infectious Diseases (ECCMID 2018, Madrid, Spain). The VOYAG-R Study Group: Antoine Andremont, Laurence Armand-Lefèvre, Olivier Bouchaud, Yacine Boussadia, Pauline Campa, Bruno Coignard, Paul-Henri Consigny, Assiya El Mniai, Marina Esposito-Farèse, Candice Estellat, Pierre-Marie Girard, Catherine Goujon, Isabelle Hoffmann, Guillaume Le Loup, JeanChristophe Lucet, Sophie Matheron, Nabila Moussa, Marion Perrier, Gilles Pialoux, Pascal Ralaimazava, Etienne Ruppé, Daniel Vittecoq, Ingrid Wieder and Benjamin Wyplosz.

**Author contributions**

L.A.L. was involved in conceptualization of the project, project management, methodology, resources, data curation, data analysis, visualization and writing of the original draft of the manuscript. E.Ro. was involved in the conceptualization of the project, methodology and formal analysis. D.D. and T.C. were involved in the development of bioinformatics pipelines and in bioinformatic analysis. J.M. and J.N. were involved in statistical analysis. O.C. was involved in data analysis. M.P. was involved in bioinformatics analysis. C.B. and L.M. were involved in methodology and formal analysis of whole-genome sequencing. O.T. was involved in mutation analysis and visualization. S.M. was PI of the VOYAG-R clinical trial and was involved in the generation of metadata. E.R. was involved in the VOYAG-R clinical trial, data analysis, visualization and reviewing the drafts of the manuscript. The VOYAG-R study group performed the princeps's study. A.A. was involved in the VOYAG-R clinical trial, conceptualization of the project, data analysis, visualization and reviewing the drafts of the manuscript. E.D. and S.P.K. were involved in conceptualization of the project, methodology, data analysis, visualization and reviewing the drafts of the manuscript.

**Conflicts of interest**

The authors declare that there are no conflicts of interest.

**Ethical statement**

The VOYAG-R study (clinicaltrials.gov number NCT01526187) was approved by the Ile de France IV ethics committee on 14 November 2011. All individuals have given written consent to participate in the study.

**References**

1. Woerther PL, Burdet C, Chachaty E, Andremont A. Trends in human fecal carriage of extended-spectrum β-Lactamases in the community: toward the globalization of CTX-M. *Clin Microbiol Rev* 2013;26:744–758.

2. Ruppé E, Armand-Lefèvre L, Estellat C, Consigny PH, El Mniai A. High rate of acquisition but short duration of carriage of multidrug-resistant *Enterobacteriaceae* after travel to the tropics. *Clin Infect Dis* 2015;61:593–600.

3. Arcilla MS, van Hattem JM, Haverkate MR, Bootsma MCJ, van Genderen PJJ. Import and spread of extended-spectrum
β-lactamase-producing *Enterobacteriaceae* by international travellers (COMBAT study): a prospective, multicentre cohort study. *Lancet Infect Dis* 2017;17:78–85.

4. Schaumburg F, Sertic SM, Correa-Martinez C, Mellmann A, Köck R, et al. Acquisition and colonization dynamics of antimicrobial-resistant bacteria during international travel: A prospective cohort study. *Clin Microbiol Infect* 2019;25:e1-1287.e7:1287.

5. Leo S, Lazarevic V, Gaïa N, Estellat C, Girard M. The intestinal microbiota predisposes to traveler's diarrhea and to the carriage of multidrug-resistant *Enterobacteriaceae* after traveling to tropical regions. *Gut Microbes* 2019;10:631–641.

6. Tenaillon O, Skurnik D, Picard B, Denamur E. The population genetics of commensal *Escherichia coli*. *Nat Rev Microbiol* 2010;8:207–217.

7. Nowrouzian F, Hesselmar B, Saalman R, Strannegård IL, Åberg N. *Escherichia coli* in infants' intestinal microflora: colonization rate, strain turnover, and virulence gene carriage. *Pediatr Res* 2003;54:8–14.

8. Ostblom A, Adlerberth I, Wold AE, Nowrouzian FL. Pathogenicity island markers, virulence determinants malX and usp, and the capacity of *Escherichia coli* to persist in infants' commensal microbiotas. *Appl Environ Microbiol* 2011;77:2303–2308.

9. Nowrouzian FL, Adlerberth I, Wold AE. Enhanced persistence in the colonic microbiota of *Escherichia coli* strains belonging to phylogenetic group B2: role of virulence factors and adherence to colonic cells. *Microbes Infect* 2006;8:834–840.

10. Clermont O, Christenson JK, Denamur E, Gordon DM. The Clermont *Escherichia coli* phylo-typing method revisited: improvement of specificity and detection of new phylo-groups. *Environ Microbiol Rep* 2013;5:58–65.

11. Desvillechabrol D, Bouchier C, Kennedy S, Cokelaer T. (n.d.) Sequana coverage: detection and characterization of genomic variations using running median and mixture models. *Gigascience*;7.

12. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. *J Comput Biol* 2012;19:455–477.

13. Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. *Bioinformatics* 2013;29:1072–1075.

14. Treangen TJ, Ondov BD, Koren S, Phillippy AM. The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. *Genome Biol* 2014;15:524.

15. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. *Syst Biol* 2010;59:307–321.

16. Bourrel AS, Poirel L, Royer G, Darty M, Vuillemin X. Colistin resistance in Parisian inpatient faecal *Escherichia coli* as the result of two distinct evolutionary pathways. *J Antimicrob Chemother* 2019;74:1521–1530.

17. Royer G, Decousser JW, Branger C, Dubois M, Médigue C. (n.d.) PlaScope: a targeted approach to assess the plasmidome from genome assemblies at the species level. *Microb Genom*;4.

18. Cokelaer T, Desvillechabrol D, Legendre R, Cardon M. Sequana': a Set of Snakemake NGS pipelines. *JOSS* 2017;2:352.

19. Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. *Bioinformatics* 2010;26:589–595.

20. Ondov BD, Treangen TJ, Melsted P, Mallonee AB, Bergman NH. Mash: fast genome and metagenome distance estimation using MinHash. *Genome Biol* 2016;17:132.

21. Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). *Nucleic Acids Res* 2014;42:D206–214.

22. Brettin T, Davis JJ, Disz T, Edwards RA, Gerdes S. RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. *Sci Rep* 2015;5:8365.

23. Vallenet D, Belda E, Calteau A, Cruveiller S, Engelen S. Micro-Scope--an integrated microbial resource for the curation and comparative analysis of genomic and metabolic data. *Nucleic Acids Res* 2013;41:D636–647.

24. Ghalayini M, Magnan M, Dion S, Zatout O, Bourguignon L. Long-term evolution of the natural isolate of *Escherichia coli* 536 in the mouse gut colonized after maternal transmission reveals convergence in the constitutive expression of the lactose operon. *Mol Ecol* 2019;28:4470–4485.

25. Galardini M, Clermont O, Baron A, Busby B, Dion S, et al. Major role of iron uptake systems in the intrinsic extra-intestinal virulence of the genus *Escherichia* revealed by a genome-wide association study. *PLoS Genet* 2020;16:e1009065.

26. Denamur E, Clermont O, Bonacorsi S, Gordon D. The population genetics of pathogenic *Escherichia coli*. *Nat Rev Microbiol* 2021;19:37–54.

27. van Duijkeren E, Wielders CCH, Dierikx CM, van Hoek AHAM, Hengeveld P, et al. Long-term carriage of extended-spectrum β-lactamase-producing *Escherichia coli* and *Klebsiella pneumoniae* in the general population in the Netherlands. *Clin Infect Dis* 2018;66:1368–1376.

28. Löhr IH, Rettedal S, Natås OB, Naseer U, Øymar K. Long-term faecal carriage in infants and intra-household transmission of CTX-M-15-producing *Klebsiella pneumoniae* following a nosocomial outbreak. *J Antimicrob Chemother* 2013;68:1043–1048.

29. Birgy A, Cohen R, Levy C, Bidet P, Courroux C. Community faecal carriage of extended-spectrum beta-lactamase-producing Enterobacteriaceae in french children. *BMC Infect Dis* 2012;12:315.

30. Nicolas-Chanoine MH, Gruson C, Bialek-Davenet S, Bertrand X, Thomas-Jean F. 10-Fold increase (2006–11) in the rate of healthy subjects with extended-spectrum β-lactamase-producing *Escherichia coli* faecal carriage in a Parisian check-up centre. *J Antimicrob Chemother* 2013;68:562–568.

31. Nilsen E, Haldorsen BC, Sundsfjord A, Simonsen GS, Ingebretsen A. Large IncHI2-plasmids encode extended-spectrum β-lactamases (ESBLs) in *Enterobacter* spp. bloodstream isolates, and support ESBL-transfer to *Escherichia coli*. *Clin Microbio Infect* 2013;19:E516–E518.

32. Rashid H, Rahman M. Possible transfer of plasmid mediated third generation cephalosporin resistance between *Escherichia coli* and *Shigella sonnei* in the human gut. *Infect Genet Evol* 2015;30:15–18.

33. Ghalayini M, Launay A, Bridier-Nahmias A, Clermont O, Denamur E, et al. Evolution of a dominant natural isolate of *Escherichia coli* in the human gut over the course of a year suggests a neutral evolution with reduced effective population size. *Appl Environ Microbiol* 2018;84.

34. Tenaillon O, Barrick JE, Ribeck N, Deatherage DE, Blanchard JL. Tempo and mode of genome evolution in a 50,000-generation experiment. *Nature* 2016;536:165–170.

35. Good BH, McDonald MJ, Barrick JE, Lenski RE, Desai MM. The dynamics of molecular evolution over 60,000 generations. *Nature* 2017;551:45–50.

36. Mahérault A-C, Kemble H, Magnan M, Gachet B, Roche D, et al. Advantage of the F2:A1:B-IncF pandemic plasmid over IncC plasmids in *in vitro* acquisition and evolution of blaCTX-M gene-bearing plasmids in *Escherichia coli*. *Antimicrob Agents Chemother* 2019;63.

37. Rogers BA, Kennedy KJ, Sidjabat HE, Jones M, Collignon P. Prolonged carriage of resistant *E. coli* by returned travellers: clonality, risk factors and bacterial characteristics. *Eur J Clin Microbiol Infect Dis* 2012;31:2413–2420.

38. Schaufler K, Semmler T, Wieler LH, Trott DJ, Pitout J, et al. Genomic and functional analysis of emerging virulent and multidrug-resistant *Escherichia coli* lineage sequence type 648. *Antimicrob Agents Chemother* 2019;63.

39. Dautzenberg MJD, Haverkate MR, Bonten MJM, Bootsma MCJ. Epidemic potential of *Escherichia coli* ST131 and *Klebsiella pneumoniae* ST258: a systematic review and meta-analysis. *BMJ Open* 2016;6:e009971.
40. Hu YY, Cai JC, Zhou HW, Chi D, Zhang XF, et al. Molecular typing of CTX-M-producing *Escherichia coli* isolates from environmental water, swine feces, specimens from healthy humans, and human patients. *Appl Environ Microbiol* 2013;79:5988-5996.

41. Massot M, Daubié A-S, Clermont O, Jauréguy F, Couffignal C, et al. Phylogenetic, virulence and antibiotic resistance characteristics of commensal strain populations of *Escherichia coli* from community subjects in the Paris area in 2010 and evolution over 30 years. *Microbiology (Reading)* 2016;162:642-650.

42. Servin AL. Pathogenesis of human diffusely adhering *Escherichia coli* expressing Afa/Dr adhesins (Afa/Dr DAEC): current insights and future challenges. *Clin Microbiol Rev* 2014;27:823-869.

43. Ulett GC, Valle J, Beloin C, Sherlock O, Ghigo JM. Functional analysis of antigen 43 in uropathogenic *Escherichia coli* reveals a role in long-term persistence in the urinary tract. *Infect Immun* 2007;75:3233-3244.

44. Blackburn D, Husband A, Saldaña Z, Nada RA, Klena J. Distribution of the *Escherichia coli* common pilus among diverse strains of human enterotoxigenic E. *J Clin Microbiol* 2009;47:1781-1784.

45. Navarro-Garcia F, Ruiz-Perez F, Cataldi Á, Larzábal M. Type VI secretion system in pathogenic *Escherichia coli*: Structure, role in virulence, and acquisition. *Front Microbiol* 2019;10:1965.

46. Dautin N. Serine protease autotransporters of enterobacteriaceae (SPATEs): biogenesis and function. *Toxins (Basel)* 2010;2:1179-1206.

47. Smajs D, Weinstock GM. The iron- and temperature-regulated cjrBC genes of *Shigella* and enteroinvasive *Escherichia coli* strains code for colicin Js uptake. *J Bacteriol* 2001;183:3958-3966.

48. Le Gall T, Clermont O, Gouriou S, Picard B, Nassif X. Extraintestinal virulence is a coincidental by-product of commensalism in B2 phylogenetic group *Escherichia coli* strains. *Mol Biol Evol* 2007;24:2373-2384.

49. Diard M, Garry L, Selva M, Mosser T, Denamur E, et al. Pathogenicity-associated islands in extraintestinal pathogenic *Escherichia coli* are fitness elements involved in intestinal colonization. *J Bacteriol* 2010;192:4885-4893.

50. Leimbach A, Hacker J, Dobrindt UE. coli as an all-rounder: the thin line between commensalism and pathogenicity. *Curr Top Microbiol Immunol* 2013;358:3-32.

**Five reasons to publish your next article with a Microbiology Society journal**

1. The Microbiology Society is a not-for-profit organization.

2. We offer fast and rigorous peer review – average time to first decision is 4–6 weeks.

3. Our journals have a global readership with subscriptions held in research institutions around the world.

4. 80% of our authors rate our submission process as 'excellent' or 'very good'.

5. Your article will be published on an interactive journal platform with advanced metrics.

**Find out more and submit your article at microbiologyresearch.org.**