# WGS Analysis of Clonal and Plasmidic Epidemiology of Colistin-Resistance Mediated by mcr Genes in the Poultry Sector in Lebanon

Hiba Al-Mir, Marwan Osman, Antoine Drapeau, Monzer Hamze, Jean-Yves

Madec, Marisa Haenni

► To cite this version:

Hiba Al-Mir, Marwan Osman, Antoine Drapeau, Monzer Hamze, Jean-Yves Madec, et al.. WGS Analysis of Clonal and Plasmidic Epidemiology of Colistin-Resistance Mediated by mcr Genes in the Poultry Sector in Lebanon. Frontiers in Microbiology, 2021, 12, pp.624194. 10.3389/fmicb.2021.624194. anses-03987553

HAL Id: anses-03987553

https://anses.hal.science/anses-03987553

Submitted on 14 Feb 2023

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

# WGS Analysis of Clonal and Plasmidic Epidemiology of Colistin-Resistance Mediated by *mcr* Genes in the Poultry Sector in Lebanon

Hiba Al-Mir¹,², Marwan Osman¹, Antoine Drapeau², Monzer Hamze¹, Jean-Yves Madec² and Marisa Haenni²*

¹ Laboratoire Microbiologie Santé et Environnement, Doctoral School of Sciences and Technology, Faculty of Public Health, Lebanese University, Tripoli, Lebanon, ² Université de Lyon – ANSES Laboratoire de Lyon, Unité Antibiorésistance et Virulence Bactériennes, Lyon, France

**OPEN ACCESS**

**Edited by:**

Azucena Mora Gutiérrez,
University of Santiago
de Compostela, Spain

**Reviewed by:**

Remy A. Bonnin,
Université Paris-Saclay, France
Patrick Rik Butaye,
Ross University School of Veterinary
Medicine, Saint Kitts and Nevis

*Correspondence:

Marisa Haenni
marisa.haenni@anses.fr

**Specialty section:**

This article was submitted to
Antimicrobials, Resistance
and Chemotherapy,
a section of the journal
Frontiers in Microbiology

**Received:** 30 October 2020

**Accepted:** 11 February 2021

**Published:** 08 March 2021

**Citation:**

Al-Mir H, Osman M, Drapeau A, Hamze M, Madec J-Y and Haenni M (2021) WGS Analysis of Clonal and Plasmidic Epidemiology of Colistin-Resistance Mediated by mcr Genes in the Poultry Sector in Lebanon. *Front. Microbiol.* 12:624194.
doi: 10.3389/fmicb.2021.624194

Poultry and poultry meat are important contributors to the global antimicrobial burden. Unregulated and illegal use of extended-spectrum cephalosporins (ESC) in this sector has long been identified as a major cause of massive spread of ESC-resistant *Escherichia coli*, and colistin usage is considered a main driver of plasmid-mediated *mcr* genes dissemination. In Lebanon, the first *mcr-1*-positive *E. coli* found in poultry dates back to 2015, followed by a few reports of *mcr-1*-positive *E. coli* in poultry, swine, humans, and the environment. On the contrary, a comprehensive picture of the population structure of *mcr-1*-positive *E. coli* and *mcr-1*-bearing plasmids carrying the *mcr-1* gene using whole-genome analysis is largely lacking. This study reports the prevalence of *mcr-1*-positive *E. coli* in poultry originating from 32 farms across three Lebanese governorates and slaughtered in the same place. We report 27/32 (84.4%) *mcr-1* positive farms, leading to a total of 84 non-duplicate *E. coli* collected, of which 62 presented the *mcr-1* gene. Numerous associated resistances were identified, including to ESC through the presence of *bla*<sub>CTX-M</sub> or *bla*<sub>CMY</sub> genes. The *mcr-1* gene was mostly carried by IncX4 (n = 36) and IncI2 (n = 24) plasmids, which are both known for their efficient transfer capacities. A high genetic diversity was detected, arguing for the lack of contamination during the slaughter process. ST744 and ST1011 were the most widely identified clones, which have been both regularly associated to *mcr-1*-carrying *E. coli* and to the poultry sector. The wide dissemination of colistin-resistance, coupled to resistances to ESC and numerous other molecules, should urge authorities to implement efficient guidelines for the use of antibiotics in the poultry sector in Lebanon.

**Keywords:** *mcr-1*, poultry, Lebanon, IncX4, IncI2
# INTRODUCTION

Since the discovery of the plasmid-mediated colistin-resistance gene *mcr-1* in 2015 (Liu et al., 2016), this gene has been extensively described in numerous animal settings, and notably in the poultry sector worldwide (Apostolakos and Piccirillo, 2018). Besides living animals, the *mcr-1* gene has also been detected in retail meat, suggesting a possible transfer to humans through under-cooked meat or cross-contamination (Nishino et al., 2017; Budel et al., 2020). To date, up to 10 *mcr* gene variants (*mcr-1* to *mcr-10*) have been recognized from different sources, of which *mcr-1* and *mcr-3* are the most widespread (Nang et al., 2019). Moreover, most studies have reported *mcr* genes in association with other resistance genes, including to critically important antimicrobials such as extended-spectrum cephalosporins (ESC) (Grami et al., 2016; Maciuca et al., 2019). In the poultry sector, *mcr-1* gene has mostly been found located on IncX4 and IncI2 plasmids, and to a lesser extent on IncH12 plasmids (Perreten et al., 2016). The usually high transfer capacity of both IncX4 and IncI2 most probably explains their wide geographical dissemination and their occurrence in a large variety of hosts, both human and animals. Colistin-resistance has often been studied under the prism of plasmid-mediated resistance, so that only a few studies reported the characterization of *mcr*-negative but colistin-resistant isolates and the role of PmrAB and PhoPQ mutations (Quesada et al., 2015).

Regarding animals in Lebanon, *mcr-1*-mediated colistin-resistance was first reported in 2015 in poultry (Dandachi et al., 2018). A single *mcr-1*-positive *Escherichia coli* was recovered from one rectal swab over 982 samples (0.1%) taken in 49 farms for surveillance purposes. In 2017, a follow-up study was performed in this first *mcr-1*-positive farm, which showed that colistin-resistance had widely disseminated since 181/200 chicken and 6/6 workers carried *mcr-1*-producing *E. coli* or *Klebsiella pneumoniae*, while litter and feed were less heavily contaminated (6 and 20%, respectively) (Dandachi et al., 2020). In 2017 as well, 23 over 114 fecal samples from swine were resistant to colistin due to the presence of the *mcr-1* gene, and four of these isolates co-harbored resistances to ESC (Dandachi et al., 2019). Between 2017 and 2018, 88/93 (94.6%) fecal samples collected from the three major poultry farms in Lebanon presented the *mcr-1* gene, among which 35.5% co-harbored a *bla*<sub>CTX-M</sub> gene (Hmede and Kassem, 2018). The presence of the *mcr-1* gene was also detected in water samples, from either irrigation water or the Mediterranean Sea (Hmede et al., 2019; Sourenian et al., 2020), suggesting environmental contamination. Finally, *mcr-1* has also been reported in human clinical isolates in Lebanon that had been collected as early as 2011 (Al-Mir et al., 2019). The plasmidic location of the *mcr-1* gene has, however, rarely been investigated in those studies: it was found located on an IncX4 plasmid in human isolates and in a rainbow trout isolate, and on IncX4 and IncI2 in water isolates (Al-Mir et al., 2019; Hassan et al., 2020; Sourenian et al., 2020). Likewise, data on the bacterial population structure hosting *mcr-1*-mediated colistin-resistance remain largely unknown.

The goal of this study was thus to look for the presence of colistin-resistance in chicken fecal samples collected

from 32 chicken farms located in three governorates of Lebanon. Based on both phenotypic and molecular analyses (including next-generation sequencing), we characterized the population structure of colistin-resistant *E. coli*, the genetic support of *mcr*-dependent or *mcr*-independent colistin-resistance, and the plasmid types carrying the *mcr* genes in order to lay the foundations for a better understanding of colistin-resistance spread in poultry, but also in humans and the environment in Lebanon.

# MATERIALS AND METHODS

## Ethics

This investigation was approved by the Azm Center/Lebanese University ethical committee (document CE-EDST-3-2018), authorized by the Lebanese Ministry of Public Health.

## Bacterial Isolation and Identification

Between May and August 2018, poultry samples were collected in one big slaughterhouse in Chekka, North Lebanon. Animals were originated from 32 individual farms hosting from 4500 to 195,000 chicken individuals, and located in seven districts from Akkar, North Lebanon, and Mount Lebanon governorates (Figure 1). At slaughterhouse, five different samples from each farm were collected: three caeca and two necks corresponding to five different chicken carcasses were randomly sampled with all precautions needed to avoid inter-sample contamination. Each farm was sampled once, except farm 2 that was sampled twice. All samples were put in plastic bags, conserved at 4°C, and rapidly transported to the Laboratoire Microbiologie Santé et Environnement (LMSE) in Tripoli, Lebanon. Resistant Enterobactérales were isolated on MacConkey agar (Bio-Rad, Hercules, CA, United States) supplemented with colistin (3.5 mg/L) (Sigma-Aldrich, St. Louis, MO, United States). Selective plates were incubated at 37°C for 24 h. One presumptive *E. coli* colony per morphology was arbitrarily selected from each selective plate. Identification of isolates was performed using matrix-assisted laser desorption/ionization time-of-flight MALDI TOF VITEK MS Version 3.0 (bioMérieux, Marcy L'Etoile, France).

## Antibiotic Susceptibility Testing and Phenotypic Characterization

Susceptibility testing was performed by the disc diffusion method on Mueller-Hinton agar according to the guidelines of the Antibiogram Committee of the French Society for Microbiology (CA-SFM)<fnref n="1" />. The *E. coli* ATCC 7624 strain was used as quality control. A total of 16 beta-lactam (amoxicillin, piperacillin, ticarcillin, amoxicillin-clavulanic acid, piperacillin-tazobactam, ticarcillin-clavulanic acid, cefalotine, cefuroxime, cefotaxime, ceftiofur, ceftazidime, cefoxitin, cefepime, cefquinome, aztreonam, and ertapenem) and 14 non-beta-lactam (tetracycline, kanamycin, tobramycin, gentamicin, amikacin, apramycin, netilmicin, streptomycin, florfenicol,

<footnote n="1">¹www.sfm-microbiologie.org</footnote>
<figure><img src="image_3.png" /></figure>

chloramphenicol, sulfonamides, trimethoprim, nalidixic acid, and enrofloxacin) antibiotics of both veterinary and human interest were tested. Minimum inhibitory concentrations (MICs) were determined by broth microdilution for colistin, according to the European Committee for Antimicrobial Susceptibility Testing (EUCAST).

## Molecular Typing of the Isolates

The detection of the major *E. coli* phylogenetic groups (A, B1, B2, or D) was performed as described by Doumith et al. (2012). Duplicate isolates collected from the same farm were detected by multiple-locus variable-number tandem-repeat analysis (MLVA) using the multiplex-based PCR described by Caméléna et al. (2019).

## Whole-Genome Sequencing (WGS)

DNA was extracted using the NucleoSpin Microbial DNA extraction kit (Macherey-Nagel, Hoerdt, France) according

to the manufacturer's instructions. Library preparation was performed using the Nextera XT technology, and sequencing was performed on a NovaSeq instrument (Illumina, San Diego, CA, United States). After sequencing, reads were quality trimmed and *de novo* assembled using Shovill v1.0.4, and the quality of assemblies was assessed using QUAST v5.0.2. STs were determined using the CGE online tool MLSTFinder v2.0.4², while resistance genes and replicon content were inferred using ABRicate v1.0.1³. Mutations in *gyrA* and *parC* were searched using PointFinder (see text footnote 2). The PmrA, PmrB, PhoP, and PhoQ amino acid sequences were extracted from the assemblies of all isolates and compared using Clustal Omega with the *E. coli* K12 reference strain (NP_418537.1 for PmrA, NP_418536.1 for PmrB, NP_415648.1 for PhoP, and NP_415647.1 for PhoQ). Virulence factors (VFs) were determined using VirulenceFinder, and serotypes

<footnote n="2">²http://www.genomicepidemiology.org/</footnote>

<footnote n="3">³https://github.com/tseemann/abricate</footnote>
were determined using SeroType Finder (see text footnote 2). Avian pathogenic *E. coli* (APEC) was defined according to Johnson et al. (2008).

## Characterization of the *mcr-1*-Carrying Plasmids

The replicon content was determined from the whole-genome sequencing (WGS) data using PlasmidFinder 2.0.1 (see text footnote 2). Plasmids carrying the *mcr-1* gene were assigned *in silico* when the *mcr-1* gene was located on the same contig as the plasmidic marker. When *in silico* data showed no co-occurrence on the same contig, plasmids carrying the *mcr-1* gene were detected using PFGE-S1 gels (6 V/cm for 20 h with an angle of 120° at 14°C with pulse times ranging from 1 to 30 s) followed by Southern blot using adequate probes as previously described (Saidani et al., 2019). Plasmid co-localization was assessed by comparison between the bands corresponding to the resistance gene and those corresponding to the Inc type of the plasmid. When Southern blots did not lead to interpretable results and for all isolates that could not be typed by PFGE (smearing profile), the plasmid of interest was transferred by conjugation in an *E. coli* J53 rifampicin-resistant recipient strain. Conjugation was performed in liquid medium using rifampicin and cefotaxime (5 mg/L) or colistin (2 mg/L) to select transconjugants (TC). Only TC presenting a unique plasmid, as assessed by S1-PFGE, were further characterized by PBRT and Southern blots as described above.

## Phylogenetic Analysis

The core genome multi-locus sequence type (cgMLST) was extracted from the WGS data. The pan-genome was determined, and core gene alignments were generated, for each collection, with Roary v. 3.13.0 (Page et al., 2015) using a Protein BLAST identity of 80% and a core definition of 90%. In the first step, all assemblies were annotated *de novo* with Prokka v1.14.6 using default settings (Seemann, 2014). The Prokka annotations were provided to Roary as input. Subsequently, recombination was removed with gubbins v2.4.1 and a maximum likelihood tree was constructed from the core gene alignment produced by Roary using RAxML v.8.2.12 using default parameters. Pairwise single nucleotide polymorphism (SNP) distances were calculated from core genome alignments generated by Roary using snp-dists<fnref n="4" />. The SNP distance matrix is provided as **Supplementary Table S1**. The resulting tree for both analysis/collections was visualized using iTol v.5.5.1⁵ (Letunic and Bork, 2019).

## Statistical Analyses

Statistical analyses were performed with GraphPad Prism 6.0 (GraphPad Software, Inc., San Diego, CA, United States) using chi-square test in order to search for an association between *mcr-1* gene and ESBL or AmpC genes among colistin-resistant isolates. The tests were two-sided, with a type I error set at $\alpha = 0.05$.

## Accession Number(s)

The whole genome shotgun project was deposited in DDBJ/EMBL/GenBank under the BioProject accession number PRJNA671785.

## RESULTS

### Detection of Colistin-Resistance

All presumptive colistin-resistant *E. coli* isolated on selective plates was characterized according to their resistance phenotype, phylogroup, and MLVA profiles, and only non-duplicate isolates were kept for further studies. Among the 32 farms tested, 27 presented at least one colistin-resistant isolate (27/32, 84.4%), with MICs ranging from 4 to 64 mg/L. Since several different isolates were retrieved from one farm (up to 12), and also from one animal (up to three), a total of 84 colistin-resistant isolates were collected (**Supplementary Table S1**).

### E. coli Characterization and Virulence Patterns

The 84 *E. coli* isolates belonged to phylogroups A (*n* = 44, 52.4%), B1 (*n* = 22, 26.2%), B2 (*n* = 1, 1.2%), and D (*n* = 17, 20.2%). Thirty-one different serogroups were identified, but only O21 (*n* = 14), O102 (*n* = 8), O101 (*n* = 6), O38 (*n* = 5), O1 (*n* = 3), and O109 (*n* = 3) were found in more than two isolates. Twenty-four isolates (24/84, 28.6%) could be considered as APEC according to the definition by Johnson et al. (2008), i.e., concomitantly presenting the *iss*, *iutA*, *hlyF*, *iroN*, and *ompT* virulence genes (**Supplementary Table S1**).

### Plasmid-Mediated Colistin-Resistance

Whole-genome sequencing data revealed that no other *mcr* gene than *mcr-1* was present in the *E. coli* genomes. The *mcr-1* gene was detected in 62 *E. coli* isolates (62/84, 73.8%) originating from 21 different farms. Two copies of *mcr-1* were found in one isolate, located on the same contig. The plasmidic location of the *mcr-1* gene was proved in 61/62 of the isolates, mostly by deduction from WGS data when *mcr-1* and Inc genes were located on the same contig (37 isolates) or next to conjugation and Southern blot experiments (14 isolates). For one isolate, the plasmidic or chromosomal location of *mcr-1* was not resolved despite numerous attempts (absence of conjugation and smearing PFGE profile). In all, the *mcr-1* gene was mostly carried by an IncX4 plasmid (*n* = 36), followed by an IncI2 (*n* = 24) and IncH12 (*n* = 1) plasmid.

### Chromosome-Mediated Colistin-Resistance

Amino-acid variations in the PhoP, PhoQ, PmrA, and PmrB were extracted from the WGS data of all *mcr-1*-positive and *mcr-1*-negative isolates (**Supplementary Table S1** and **Table 1**). Modifications were found at one site in PhoP, 11 sites in PhoQ, two sites in PmrA, and six sites in PmrB. Four modifications (H2R and D283G in PmrB, I44L in PhoP, and A482T in PhoQ) were found in ≥25 isolates. Eleven modifications (one in PmrA

<footnote n="4">⁴https://github.com/tseemann/snp-dists</footnote>

<footnote n="5">⁵http://itol.embl.de/itol.cgi</footnote>
<table><caption>TABLE 1 | Amino-acid modifications (PmrA, PmrB, PhoP, and PhoQ) in mcr-1-positive and mcr-1-negative isolates.</caption><thead><tr><th></th><th>PmrA</th><th>PmrB</th><th>PhoP</th><th>PhoQ</th></tr></thead><tbody><tr><td>G29S</td><td>V129L</td><td>H2R</td><td>A242V</td><td>A360V</td></tr><tr><td>mcr-1-positive isolates</td><td>1</td><td>2</td><td>61</td><td>15</td></tr><tr><td>mcr-1-negative isolates</td><td>1</td><td>0</td><td>21</td><td>1</td></tr><tr><td>Total</td><td>2</td><td>2</td><td>82</td><td>16</td></tr></tbody></table><table><thead><tr><th></th><th>I44L</th><th>I175F</th><th>V228I</th><th>A341T</th><th>N346K</th><th>T348I</th><th>T348N</th><th>V373I</th><th>V386L</th><th>A390T</th><th>E464D</th><th>L467M</th><th>A482T</th></tr></thead><tbody><tr><td>mcr-1-positive isolates</td><td>1</td><td>1</td><td>0</td><td>1</td><td>1</td><td>0</td><td>2</td><td>1</td><td>1</td><td>1</td><td>19</td><td>2</td><td>30</td></tr><tr><td>mcr-1-negative isolates</td><td>1</td><td>0</td><td>0</td><td>3</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Total</td><td>2</td><td>2</td><td>44</td><td>17</td><td>18</td><td>25</td><td>4</td><td>2</td><td>3</td><td>1</td><td>1</td><td>19</td><td>2</td><td>30</td></tr></tbody></table>

and 10 in PhoQ) were found only in *mcr-1*-positive isolates (Table 1), while two modifications in PhoQ (three isolates with the V228I modification and one isolate with the T348I modification) were exclusively found in *mcr-1*-negative isolates. However, no pattern and not even one specific mutation were found to be associated with all 22 *mcr-1*-negative colistin-resistant isolates.

## Associated Resistance Genes

Among the 62 *mcr-1*-positive *E. coli* isolates, 30 presented an associated ESBL gene, including *bla*<sub>CTX-M-65</sub> (*n* = 15), *bla*<sub>CTX-M-3</sub> (*n* = 9), *bla*<sub>CTX-M-14</sub> (*n* = 2), *bla*<sub>CTX-M-15</sub> (*n* = 2), *bla*<sub>CTX-M-64</sub> (*n* = 1), and *bla*<sub>SHV-12</sub> (*n* = 1) (Figure 2). No ESBL/*mcr-1* co-localization could be evidenced, neither by conjugation nor by analyzing WGS data. Plasmidic AmpC genes were also detected in 17 isolates, with *bla*<sub>CMY-2</sub> in 16 of them. All but one CTX-M-3-positive *E. coli* isolate also displayed the CMY-2 enzyme. On the contrary, only one CTX-M-65 and nine CMY-2-producing *E. coli* were identified among *mcr-1*-negative colistin-resistant isolates. The presence of the ESBL genotype was 20-fold more common in *mcr-1*-positive than in *mcr-1*-negative colistin-resistant isolates (OR: 19.7, CI: 2.5–155.6, *P* = 0.0003). On the contrary, no significant difference related to AmpC phenotype was observed between the two aforementioned colistin-resistant groups.

In addition to genes responsible for colistin-resistance, *E. coli* isolates mostly presented resistance genes to aminoglycosides (75/84, 89.3%), tetracyclines (71/84, 84.5%), sulfonamides (66/84, 78.6%), trimethoprim (62/84, 73.8%) and phenicol (53/84, 63.1%). The corresponding genes were mainly *tet(A)*, *sul2*, *dfrA12*, and *floR*, while aminoglycoside resistance was largely multifactorial (Supplementary Table S1). Quinolone-resistance was observed in 70 isolates, of which 56 were also resistant to fluoroquinolones. Variants of the plasmid-mediated resistance gene *qnrB* (*n* = 4), as well as the *qnrS1* (*n* = 14) gene were detected in 14 different isolates. The mutations S83L and D87N in *GyrA* and S80I in *ParC* were concomitantly detected in 47/56 fluoroquinolone-resistance isolates (Supplementary Table S1). Only one apramycin-resistant isolate was identified, which presented the *aac(3′)-IV* gene, while none of the isolates was resistant to ertapenem and amikacin (Table 2).

## Population Structure of the mcr-1-Positive and mcr-1-Negative E. coli Isolates

The 84 isolates belonged to 38 different STs (Figure 2 and Supplementary Table S1). Two STs were predominant [ST93 (*n* = 14) and ST1011 (*n* = 11, including one SLV)], while two were each found in five different isolates (ST744 and ST1140). The remaining isolates were identified in three or less isolates. Within one ST, colistin-resistance was usually mediated by the same mechanism (either *mcr-1* or chromosomal mutations). However, six STs (ST10, ST48, ST155, ST744, ST1140, and ST1635) encompassed isolates that were resistant to colistin through either plasmid- or chromosome-mediated mechanisms.
<figure><img src="image_4.png" /></figure>

Among the STs presenting five isolates or more, only ST93 formed a homogeneous group, with isolates differing by 14–72 SNPs. On the other side, ST744, ST1011, and ST1140 were split

into two groups, each differing by > 900, > 2700 SNPs, and > 600 SNPs, respectively (Supplementary Table S2). The number of SNPs was also detailed for all isolates sharing the same ST. ST359
<table><caption>TABLE 2 | Associated resistance phenotypes.</caption><thead><tr><th rowspan="2"></th><th colspan="2">mcr-1-positive isolates (n = 62)</th><th colspan="2">mcr-1-negative isolates (n = 22)</th><th colspan="2">Total (n = 84)</th></tr><tr><th>No.</th><th>Percentage</th><th>No.</th><th>Percentage</th><th>No.</th><th>Percentage</th></tr></thead><tbody><tr><td>Streptomycin</td><td>31</td><td>50.0</td><td>17</td><td>77.3</td><td>48.0</td><td>57.1</td></tr><tr><td>Kanamycin</td><td>44</td><td>71.0</td><td>15</td><td>68.2</td><td>59.0</td><td>70.2</td></tr><tr><td>Amikacin</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>Apramycin</td><td>1</td><td>1.6</td><td>0</td><td>0.0</td><td>1.0</td><td>1.2</td></tr><tr><td>Gentamicin</td><td>26</td><td>41.9</td><td>9</td><td>40.9</td><td>35.0</td><td>41.7</td></tr><tr><td>Tobramycin</td><td>24</td><td>38.7</td><td>9</td><td>40.9</td><td>33.0</td><td>39.3</td></tr><tr><td>Netilmicin</td><td>11</td><td>17.7</td><td>6</td><td>27.3</td><td>17.0</td><td>20.2</td></tr><tr><td>Chloramphenicol</td><td>42</td><td>67.7</td><td>11</td><td>50.0</td><td>53.0</td><td>63.1</td></tr><tr><td>Florfenicol</td><td>35</td><td>56.5</td><td>9</td><td>40.9</td><td>44.0</td><td>52.4</td></tr><tr><td>Tetracyclines</td><td>51</td><td>82.3</td><td>17</td><td>77.3</td><td>68.0</td><td>81.0</td></tr><tr><td>Sulfonamides</td><td>46</td><td>74.2</td><td>18</td><td>81.8</td><td>64.0</td><td>76.2</td></tr><tr><td>Trimethoprim</td><td>47</td><td>75.8</td><td>17</td><td>77.3</td><td>64.0</td><td>76.2</td></tr><tr><td>Nalidixic acid</td><td>53</td><td>85.5</td><td>17</td><td>77.3</td><td>70.0</td><td>83.3</td></tr><tr><td>Enrofloxacin</td><td>47</td><td>75.8</td><td>9</td><td>40.9</td><td>56.0</td><td>66.7</td></tr></tbody></table>

presented two highly similar isolates (25 SNPs) originating from two different animals from the same farm, while ST616 and ST1196 each comprised two isolates from different farms that only differed by 3 SNPs. On the contrary, ST48 and ST6856 each included two isolates from different farms, respectively, differing by > 5000 and 391 SNPs. ST4162 comprised three isolates from different farms, of which two (farms F1 and F32) were fully identical, while the third one (F8) differed by 27 SNPs. ST155 encompassed three isolates from three different farms: isolates from farms F16 and F29 displayed the *mcr-1* gene and differed by 34 SNPs, while the isolate from farm F18 was *mcr-1*-negative and differed from the two others by > 5000 SNPs. ST1638 comprised three isolates: two were from the same farm F9 and were highly similar (5 SNPs), while the third one from farm F22 clearly diverged (>5000 SNPs). Finally, ST10 included three isolates from three different farms that were not similar, isolates from farm F2 and F20 differing by > 200 SNPs, and further differing by > 6000 SNPs from the isolate from farm F26.

## DISCUSSION

This study first reveals a massive spread of colistin-resistance in poultry farms in Lebanon (27/32, 84.4%), with the wide dissemination of the *mcr-1* gene (21/27 positive farms). The poultry sector is the very first one where *mcr-1* was detected in Lebanon (Dandachi et al., 2018), and a recent study in three poultry farms reported high *mcr-1* prevalence in broilers (Hmede and Kassem, 2018). The present work on 32 farms therefore expands our knowledge on the magnitude of *mcr-1* spread in the poultry sector in Lebanon. This high prevalence can most probably be explained by the unregulated use of colistin in this country in the poultry sector, without implementation of antibiotic stewardship programs (Kassem et al., 2019). Even if we have no specific information on colistin use in the farms from which our samples originated, it is known

that colistin is easily available without the requirement of a veterinarian's prescription. Around 12 different drug brands that contain colistin are legally available over the counter in agriculture stores in Lebanon and particularly advised for the treatment and prevention of infections in poultry. Consequently, there is a high probability that colistin had been used in several of the tested farms. Additional ESBL/pAmpC genes were detected in 49/84 *E. coli* isolates, and ESBL genes were statistically more associated with *mcr-1*-positive isolates. These 49 ESBL/pAmpC-positive isolates, as well as 28 of the 35 remaining isolates, were multidrug-resistant (resistance to three or more antibiotic families). This suggests that colistin-resistance would be easily co-selected in the poultry gut by the use of most of the veterinary-licensed antibiotics. The *bla*<sub>ESBL</sub> genes found in these isolates (*bla*<sub>CTX-M-3</sub> and *bla*<sub>CTX-M-65</sub>) do not betray a human origin, where the *bla*<sub>CTX-M-15</sub> gene (which was only detected in two poultry isolates) is clearly dominating.

The mechanism of colistin-resistance in *mcr*-negative isolates was not elucidated here. Numerous mutations were observed compared to the K12 reference strain, most of which were detected in both *mcr*-positive and *mcr*-negative isolates, if not only in *mcr*-positive ones. The rare mutations that have been associated with colistin-resistance (R81S in PmrA; T156K, A159V, G161V in PhoP; E375K in PhoQ) were not detected in our dataset (Jeannot et al., 2017). Likewise, the N346K modification highlighted by Luo et al. (2017) as a possible colistin-resistance-related modification was found here associated to *mcr-1*-positive isolates. We might hypothesize that colistin-resistance-related mutations in PmrAB and PhoPQ may also arise in *mcr-1*-positive isolates and contribute to an increased MIC to colistin, but the highest MICs to colistin (≥16mg/L) observed in this study did not correlate with a specific modification pattern. Moreover, several modifications detected here (S29G in PmrA, H2R in PmrB, and D283G in PhoQ) have also been described in susceptible isolates (Luo et al., 2017). Consequently, further studies are needed to identify other genes that may be associated to colistin-resistance.

The *E. coli* population structure described here appeared very diverse, with 38 different STs detected, as it is mostly the case when samples are originating from different farms. Since all animals were sampled in the same slaughterhouse (even though poultry originated from 32 different farms), we looked for a potential contamination related to the slaughterhouse. Such a massive contamination is very unlikely because collected *E. coli* were genetically diverse, which rather indicates multiple origins. Nevertheless, a one-source contamination cannot be excluded in the cases of ST616 and ST1196 (two isolates each, originating from different farms) and ST93 (recovered in 14 isolates collected from 10 different farms), which only differed by a few SNPs. Clonal spread of *mcr-1*-carrying ST93 clinical isolates has been described in companion animals attending a veterinary hospital in China (Wang et al., 2018), suggesting that this clone may survive in the environment (a clinic or a slaughterhouse) before further dissemination. Interestingly, the *mcr-1* gene has also been detected in ST93
*E. coli* isolated from a human patient in Uruguay and from one healthy person in Finland (Grondahl-Yli-Hannuksela et al., 2018; Papa-Ezdra et al., 2020).

Besides ST93, and despite the high clonal diversity, two other main STs (ST744 and ST1011) were detected. These STs have already been both concomitantly reported in *mcr-1*-positive isolates of poultry origin in Czechia and in colistin-susceptible poultry isolates from Algeria (Belmahdi et al., 2016; Gelbicova et al., 2019). ST1011 has also been reported in *mcr-1*-positive environmental samples of swine farms in Germany, in pigs in China and Belgium, in poultry in Egypt, and in companion animals in China, as well as in colistin-susceptible poultry meat isolates in Egypt (Elnahriry et al., 2016; El Garch et al., 2017; Guenther et al., 2017; Wang et al., 2018; Ramadan et al., 2020; Shen et al., 2020). Interestingly, *mcr-1*-positive ST1011 *E. coli* isolates have also been recently identified from poultry farm workers in Lebanon (Dandachi et al., 2020). Consequently, the *mcr-1*-positive ST1011 *E. coli* isolate that was identified in 2013 in a Lebanese patient may well have a poultry (or at least an animal) origin (Al-Mir et al., 2019). ST744 has also been reported in *mcr-1*-positive isolates from poultry in Romania and from swine in China, and from *mcr-3*-producing *E. coli* in veal calves in France (Haenni et al., 2018; Maciua et al., 2019; Shen et al., 2020), but also from *mcr-1*-positive clinical *E. coli* in Portugal (Tacao et al., 2017), suggesting that this clone may be particularly prone to harbor colistin-resistance. Nonetheless, it should be kept in mind that WGS data nowadays strongly challenge any lineage distribution based on MLST only, as also highlighted in some occasions in the present work. Therefore, WGS approaches are required for any further comprehensive pictures of the cross-sectorial distribution of *mcr-1*-positive *E. coli*.

A strength of our work refers to the identification of *mcr-1*-bearing plasmids, which was almost absent from other studies in Lebanon (Al-Mir et al., 2019; Hassan et al., 2020; Sourenian et al., 2020). With the single exception of one isolate where plasmid or chromosomal location was not clarified, all *mcr-1* genes were detected on plasmids, mostly on IncX4 (*n* = 36) but also on IncI2 (*n* = 24). IncX4 is the main plasmid spreading *mcr* genes worldwide, and notably *mcr-1* (Matamoros et al., 2017). IncX is a family of small and narrow-range plasmids, and experiments proved that it has a very weak fitness cost and high transfer frequencies at 30°C, allowing its wide spread in environmental settings (Lo et al., 2014). IncI2 plasmids are also spreading efficiently globally, and a recent study proved *in situ* in a mouse model that this plasmid family had a particularly high capacity to transfer DNA in the gut (Neil et al., 2020). This high transfer capacity of both IncX4 and IncI2 plasmids may explain the occurrence of *mcr-1* in such a high number of genetic backgrounds of *E. coli*, and the relative absence of clonal transmission on farm, since animals from the same farm mostly carried different *E. coli* clones.

## CONCLUSION

We report a high prevalence and a massive spread of *mcr-1*-positive *E. coli* in poultry farms in Lebanon. Based on

WGS analysis, we deciphered that the colistin-resistance gene *mcr-1* has widely disseminated in the poultry sector in diverse genetic backgrounds of *E. coli*, and principally on IncX4 and IncI2 plasmids. Interestingly, no other *mcr* variant than *mcr-1* was found in this sector, while other still unknown non-*mcr* genes most likely also contribute to colistin-resistance. Also, a comprehensive WGS-based picture of the global clonal and plasmidic epidemiology of *mcr-1*-positive *E. coli* in a One Health perspective still lacks in all sectors to further conclude or hypothesize on major sources and routes of transmission in Lebanon. Nonetheless, some situations would warrant further investigations, such as the occurrence of ST1011 in poultry (where it is a major *E. coli* clone), in a human patient and a poultry farm worker in the same country. These results should be used to inform and increase breeders' awareness of the consequences of uncontrolled use of antibiotics in their daily practices. Overall, the wide dissemination of colistin-resistance, coupled to a low-level awareness of antibiotic stewardship in the Lebanese community (Al Omari et al., 2019) and high resistance rates to ESC and numerous other molecules, should urge authorities to implement efficient guidelines for the use of antibiotics in the poultry—and more globally the Agri-food sector in Lebanon.

## 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 in the article/Supplementary Material.

## ETHICS STATEMENT

The animal study was reviewed and approved by the Azm Center/Lebanese University Ethical Committee (document CE-EDST-3-2018), authorized by the Lebanese Ministry of Public Health.

## AUTHOR CONTRIBUTIONS

MH and MO designed the experiments and supervised the sampling campaign. HA-M performed the experiments. AD performed all WGS analyses. MH, HA-M, and J-YM analyzed the data. MH drafted the manuscript. All authors approved the final version of this manuscript.

## FUNDING

This work was supported by internal funding of the French Agency for Food, Environmental and Occupational Health and Safety (ANSES), and by a grant from the Hamidi Medical Center in Tripoli, Lebanon. HA-M was supported by a fellowship from Association AZM and Saadeh,
Lebanese Association for Scientific Research (LASER), and
Université de Lyon.

ACKNOWLEDGMENTS

The authors would like to thank Taha Abdou, Mariam Yehya,
Anas Al-Mir, Adel Al-Mir, Basel Halabieh, Veronique Metayer,
Raquel Garcia Fierro, and Wilco PM poultry slaughterhouse

and processing plant workers for their assistance in sample
collection and processing.

SUPPLEMENTARY MATERIAL

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

REFERENCES

Al Omari, S., Al Mir, H., Wrayde, S., Merhabi, S., Dhaybi, I., Jamal, S., et al. (2019). First Lebanese antibiotic awareness week campaign: knowledge, attitudes and practices towards antibiotics. *J. Hosp. Infect.* 101, 475-479. doi: 10.1016/j.jhin. 2018.07.009

Al-Mir, H., Osman, M., Azar, N., Madec, J. Y., Hamze, M., and Haenni, M. (2019). Emergence of clinical *mcr-1*-positive *Escherichia coli* in Lebanon. *J. Glob. Antimicrob. Resist.* 19, 83–84. doi: 10.1016/j.jgar.2019.08.019

Apostolakos, I., and Piccirillo, A. (2018). A review on the current situation and challenges of colistin resistance in poultry production. *Avian Pathol.* 47, 546–558. doi: 10.1080/03079457.2018.1524573

Belmahdi, M., Bakour, S., Al Bayssari, C., Touati, A., and Rolain, J. M. (2016). Molecular characterisation of extended-spectrum beta-lactamase- and plasmid AmpC-producing *Escherichia coli* strains isolated from broilers in Bejaia, Algeria. *J. Glob. Antimicrob. Resist.* 6, 108–112. doi: 10.1016/j.jgar.2016. 04.006

Budel, T., Kuenzli, E., Campos-Madueno, E. I., Mohammed, A. H., Hassan, N. K., Zinsstag, J., et al. (2020). On the island of Zanzibar people in the community are frequently colonized with the same MDR Enterobacterales found in poultry and retailed chicken meat. *J. Antimicrob. Chemother.* 75, 2432–2441. doi: 10.1093/jac/dka198

Caméléna, F., Birgy, A., Smail, Y., Courroux, C., Mariani-Kurkdjian, P., Le Hello, S., et al. (2019). Rapid and simple universal *Escherichia coli* genotyping method based on multiple-locus variable-number tandem-repeat analysis using single-tube multiplex PCR and standard gel electrophoresis. *Appl. Environ. Microbiol.* 85:e02812-18.

Dandachi, I., Fayad, E., El-Bazzal, B., Daoud, Z., and Rolain, J. M. (2019). Prevalence of extended-spectrum beta-lactamase-producing Gram-negative bacilli and emergence of *mcr-1* colistin resistance gene in Lebanese swine farms. *Microb. Drug Resist.* 25, 233–240. doi: 10.1089/mdr.2018.0110

Dandachi, I., Fayad, E., Sleiman, A., Daoud, Z., and Rolain, J. M. (2020). Dissemination of multidrug-resistant and *mcr-1* Gram-negative bacilli in broilers, farm workers, and the surrounding environment in Lebanon. *Microb. Drug Resist.* 26, 368–377. doi: 10.1089/mdr.2019.0137

Dandachi, I., Leangapichart, T., Daoud, Z., and Rolain, J. M. (2018). First detection of *mcr-1* plasmid-mediated colistin-resistant *Escherichia coli* in Lebanese poultry. *J. Glob. Antimicrob. Resist.* 12, 137–138. doi: 10.1016/j.jgar.2018. 01.004

Doumith, M., Day, M. J., Hope, R., Wain, J., and Woodford, N. (2012). Improved multiplex PCR strategy for rapid assignment of the four major *Escherichia coli* phylogenetic groups. *J. Clin. Microbiol.* 50, 3108–3110. doi: 10.1128/jcm. 01468-12

El Garch, F., Sauget, M., Hocquet, D., Lechaudee, D., Woehrle, F., and Bertrand, X. (2017). *mcr-1* is borne by highly diverse *Escherichia coli* isolates since 2004 in food-producing animals in Europe. *Clin Microbiol Infect* 23, 51.e1–51.e4.

Elnahriry, S. S., Khalifa, H. O., Soliman, A. M., Ahmed, A. M., Hussein, A. M., Shimamoto, T., et al. (2016). Emergence of plasmid-mediated colistin resistance gene *mcr-1* in a clinical *Escherichia coli* isolate from Egypt. *Antimicrob. Agents Chemother.* 60, 3249–3250. doi: 10.1128/aac.00269-16

Gelbicova, T., Barakova, A., Florianova, M., Jamborova, I., Zelendova, M., Pospisilova, L., et al. (2019). Dissemination and comparison of genetic determinants of mcr-mediated colistin resistance in *Enterobacteriaceae* via retailed raw meat products. *Front. Microbiol.* 10:2824. doi: 10.3389/fmicb.2019. 02824

Grami, R., Mansour, W., Mehri, W., Bouallegue, O., Boujaafar, N., Madec, J. Y., et al. (2016). Impact of food animal trade on the spread of *mcr-1*-mediated colistin resistance, Tunisia, July 2015. *Euro Surveill.* 21:30144.

Grondahl-Yli-Hannuksela, K., Lonnqvist, E., Kallonen, T., Lindholm, L., Jalava, J., Rantakokko-Jalava, K., et al. (2018). The first human report of mobile colistin resistance gene, *mcr-1*, in Finland. *APMIS* 126, 413–417. doi: 10.1111/apm. 12834

Guenther, S., Falgenhauer, L., Semmler, T., Imirzalioglu, C., Chakraborty, T., Roesler, U., et al. (2017). Environmental emission of multiresistant *Escherichia coli* carrying the colistin resistance gene *mcr-1* from German swine farms. *J. Antimicrob. Chemother.* 72, 1289–1292.

Haenni, M., Beyrouthy, R., Lupo, A., Chatre, P., Madec, J. Y., and Bonnet, R. (2018). Epidemic spread of *Escherichia coli* ST744 isolates carrying *mcr-3* and *blaCTX-M-55* in cattle in France. *J. Antimicrob. Chemother.* 73, 533–536. doi: 10.1093/jac/dkx418

Hassan, J., Eddine, R. Z., Mann, D., Li, S., Deng, X., Saoud, I. P., et al. (2020). The mobile colistin resistance gene, *mcr-1.1*, is carried on IncX4 plasmids in multidrug resistant *E. coli* isolated from rainbow trout aquaculture. *Microorganisms* 8:1636. doi: 10.3390/microorganisms8111636

Hmede, Z., and Kassem, I. I. (2018). The colistin resistance gene *mcr-1* is prevalent in commensal *Escherichia coli* isolated from preharvest poultry in Lebanon. *Antimicrob. Agents Chemother.* 62:e01304-18.

Hmede, Z., Sulaiman, A. A., Jaafar, H., and Kassem, I. I. (2019). Emergence of plasmid-borne colistin resistance gene *mcr-1* in multidrug-resistant *Escherichia coli* isolated from irrigation water in Lebanon. *Int. J. Antimicrob. Agents* 54, 102–104. doi: 10.1016/j.ijantimicag.2019.05.005

Jeannot, K., Bolard, A., and Plesiat, P. (2017). Resistance to polymyxins in Gram-negative organisms. *Int. J. Antimicrob. Agents* 49, 526–535. doi: 10.1016/j. ijantimicag.2016.11.029

Johnson, T. J., Wannemuehler, Y., Doetkott, C., Johnson, S. J., Rosenberger, S. C., and Nolan, L. K. (2008). Identification of minimal predictors of avian pathogenic *Escherichia coli* virulence for use as a rapid diagnostic tool. *J. Clin. Microbiol.* 46, 3987–3996. doi: 10.1128/jcm.00816-08

Kassem, I. I., Hijazi, M. A., and Saab, R. (2019). On a collision course: the availability and use of colistin-containing drugs in human therapeutics and food-animal farming in Lebanon. *J. Glob. Antimicrob. Resist.* 16, 162–164. doi: 10.1016/j.jgar.2019.01.019

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

Liu, Y. Y., Wang, Y., Walsh, T. R., Yi, L. X., Zhang, R., Spencer, J., et al. (2016). Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study. *Lancet Infect. Dis.* 16, 161–168. doi: 10.1016/s1473-3099(15)00424-7

Lo, W. U., Chow, K. H., Law, P. Y., Ng, K. Y., Cheung, Y. Y., Lai, E. L., et al. (2014). Highly conjugative IncX4 plasmids carrying *blaCTX-M* in *Escherichia coli* from humans and food animals. *J. Med. Microbiol.* 63, 835–840. doi: 10.1099/jmm.0.074021-0

Luo, Q., Yu, W., Zhou, K., Guo, L., Shen, P., Lu, H., et al. (2017). Molecular epidemiology and colistin resistant mechanism of *mcr*-positive and *mcr*-negative clinical isolated *Escherichia coli*. *Front. Microbiol.* 8:2262. doi: 10.3389/fmicb.2017.02262

Maciuka, I. E., Cummins, M. L., Cozma, A. P., Rimbu, C. M., Guguianu, E., Panzaru, C., et al. (2019). Genetic features of *mcr-1* mediated colistin resistance in CMY-2-producing *Escherichia coli* from Romanian poultry. *Front. Microbiol.* 10:2267. doi: 10.3389/fmicb.2019.02267
Matamoros, S., Van Hattem, J. M., Arcilla, M. S., Willemse, N., Melles, D. C., Penders, J., et al. (2017). Global phylogenetic analysis of *Escherichia coli* and plasmids carrying the *mcr-1* gene indicates bacterial diversity but plasmid restriction. *Sci. Rep.* 7:15364.

Nang, S. C., Li, J., and Velkov, T. (2019). The rise and spread of *mcr* plasmid-mediated polymyxin resistance. *Crit. Rev. Microbiol.* 45, 131–161. doi: 10.1080/1040841x.2018.1492902

Neil, K., Allard, N., Grenier, F., Burrus, V., and Rodrigue, S. (2020). Highly efficient gene transfer in the mouse gut microbiota is enabled by the Incl2 conjugative plasmid TP114. *Commun. Biol.* 3:523.

Nishino, Y., Shimojima, Y., Suzuki, Y., Ida, M., Fukui, R., Kuroda, S., et al. (2017). Detection of the *mcr-1* gene in colistin-resistant *Escherichia coli* from retail meat in Japan. *Microbiol. Immunol.* 61, 554–557. doi: 10.1111/1348-0421.12549

Page, A. J., Cummins, C. A., Hunt, M., Wong, V. K., Reuter, S., Holden, M. T., et al. (2015). Roary: rapid large-scale prokaryote pan genome analysis. *Bioinformatics* 31, 3691–3693. doi: 10.1093/bioinformatics/btv421

Papa-Ezdra, R., Grill Diaz, F., Vieytes, M., Garcia-Fulgueiras, V., Caiata, L., Avila, P., et al. (2020). First three *Escherichia coli* isolates harbouring *mcr-1* in Uruguay. *J. Glob. Antimicrob. Resist.* 20, 187–190. doi: 10.1016/j.jgar.2019.07.016

Perreten, V., Strauss, C., Collaud, A., and Gerber, D. (2016). Colistin resistance gene *mcr-1* in avian-pathogenic *Escherichia coli* in South Africa. *Antimicrob. Agents Chemother.* 60, 4414–4415. doi: 10.1128/aac.00548-16

Quesada, A., Porrero, M. C., Tellez, S., Palomo, G., Garcia, M., and Dominguez, L. (2015). Polymorphism of genes encoding PmrAB in colistin-resistant strains of *Escherichia coli* and *Salmonella enterica* isolated from poultry and swine. *J. Antimicrob. Chemother.* 70, 71–74.

Ramadan, H., Jackson, C. R., Frye, J. G., Hiott, L. M., Samir, M., Awad, A., et al. (2020). Antimicrobial resistance, genetic diversity and multilocus sequence typing of *Escherichia coli* from humans, retail chicken and ground beef in Egypt. *Pathogens* 9:357. doi: 10.3390/pathogens9050357

Saidani, M., Messadi, L., Chaouechi, A., Tabib, I., Saras, E., Soudani, A., et al. (2019). High genetic diversity of *Enterobacteriaceae* clones and plasmids disseminating

resistance to extended-spectrum cephalosporins and colistin in healthy chicken in Tunisia. *Microbial. Drug Resist.* 25, 1507–1513. doi: 10.1089/mdr.2019.0138

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. *Bioinformatics* 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

Shen, C., Zhong, L.-L., Yang, Y., Doi, Y., Paterson, D. L., Stoesser, N., et al. (2020). Dynamics of *mcr-1* prevalence and *mcr-1*-positive *Escherichia coli* after the cessation of colistin use as a feed additive for animals in China: a prospective cross-sectional and whole genome sequencing-based molecular epidemiological study. *Lancet Microbe* 1, e34–e43.

Sourenian, T., Mann, D., Li, S., Deng, X., Jaafar, H., Kassem, et al. (2020). Dissemination of multidrug-resistant *Escherichia coli* harboring the mobile colistin resistance gene *mcr-1.1* on transmissible plasmids in the Mediterranean Sea. *J. Glob. Antimicrob. Resist.* 22, 84–86. doi: 10.1016/j.jgar.2020.05.007

Tacao, M., Tavares, R. D. S., Teixeira, P., Roxo, I., Ramalheira, E., Ferreira, S., et al. (2017). *mcr-1* and *bla*<sub>KPC-3</sub> in *Escherichia coli* sequence type 744 after meropenem and colistin therapy. *Portugal. Emerg. Infect. Dis.* 23, 1419–1421. doi: 10.3201/eid2308.170162

Wang, J., Huang, X. Y., Xia, Y. B., Guo, Z. W., Ma, Z. B., Yi, M. Y., et al. (2018). Clonal spread of *Escherichia coli* ST93 carrying *mcr-1*-harboring IncN1-IncH12/ST3 plasmid among companion animals, China. *Front. Microbiol.* 9:2989. doi: 10.3389/fmicb.2018.02989

**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 Al-Mir, Osman, Drapeau, Hamze, Madec and Haenni. 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.