# Spread of the mcr-1 colistin-resistance gene in Escherichia coli through plasmid transmission and chromosomal transposition in French goats

Michaël Treilles, Pierre Châtre, Antoine Drapeau, Jean-Yves Madec, Marisa Haenni

► To cite this version:

Michaël Treilles, Pierre Châtre, Antoine Drapeau, Jean-Yves Madec, Marisa Haenni. Spread of the mcr-1 colistin-resistance gene in Escherichia coli through plasmid transmission and chromosomal transposition in French goats. Frontiers in Microbiology, 2023, 13, pp.1023403. 10.3389/fmicb.2022.1023403. anses-03954125

HAL Id: anses-03954125

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

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

**OPEN ACCESS**

**EDITED BY**

Rustam Aminov,
University of Aberdeen,
United Kingdom

**REVIEWED BY**

Qixia Luo,
Zhejiang University,
China
Wenjuan Yin,
Hebei University,
China
Alberto Quesada,
University of Extremadura,
Spain

*CORRESPONDENCE

Marisa Haenni

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

marisa.haenni@anses.fr

**SPECIALTY SECTION**

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

RECEIVED 19 August 2022

ACCEPTED 09 December 2022

PUBLISHED 04 January 2023

**CITATION**

Treilles M, Châtre P, Drapeau A, Madec J-Y and Haenni M (2023) Spread of the *mcr-1* colistin-resistance gene in *Escherichia coli* through plasmid transmission and chromosomal transposition in French goats.
*Front. Microbiol.* 13:1023403.
doi: 10.3389/fmicb.2022.1023403

**COPYRIGHT**

© 2023 Treilles, Châtre, Drapeau, Madec and Haenni. This is an open-access article distributed under the terms of the [Creative Commons Attribution License (CC BY)](https://creativecommons.org/licenses/by/4.0/). 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.

# Spread of the *mcr-1* colistin-resistance gene in *Escherichia coli* through plasmid transmission and chromosomal transposition in French goats

Michaël Treilles¹,², Pierre Châtre³, Antoine Drapeau³, Jean-Yves Madec³ and Marisa Haenni³*

¹Laboratoire d'Analyse Qualyse, Champdeniers Saint-Denis, France, ²Association Régionale de Prévention contre la résistance aux Antimicrobiens, Champdeniers Saint Denis, France, ³Unité Antibiorésistance et Virulence Bactériennes, Agence Nationale de Sécurité Sanitaire (ANSES) – Université de Lyon, Lyon, France

**Introduction:** Colistin-resistance widely disseminated in food-producing animals due to decades of colistin use to treat diarrhea. The plasmid-borne *mcr-1* gene has been extensively reported from bovine, swine and chicken worldwide, but smaller productions such as the goat farming sector were much less surveyed.

**Methods:** We looked for colistin-resistant isolates presenting plasmid-borne genes of the mcr family in both breeding (n=80) and fattening farms (n=5). Localization of the *mcr-1* gene was performed using Southern blot analysis coupled to short-read and long-read sequencing.

**Results:** Only the *mcr-1* gene was identified in 10% (8/80) of the breeding farms and four over the five fattening farms. In total, 4.2% (65/1561) of the animals tested in breeding farms and 60.0% (84/140) of those tested in fattening farms presented a *mcr-1*-positive *E. coli*. The *mcr-1* gene was located either on the chromosome (32.2%) or on IncX4 (38.9%) and IncHI2 (26.8%) plasmids. As expected, both clonal expansion and plasmidic transfers were observed in farms where the *mcr-1* gene was carried by plasmids. Tn6330 transposition was observed in the chromosome of diverse *E. coli* sequence types within the same farm.

**Discussion:** Our results show that the *mcr-1* gene is circulating in goat production and is located either on plasmids or on the chromosome. Evidence of Tn6330 transposition highlighted the fact that chromosomal insertion does not impair the transmission capability of the *mcr-1* gene. Only strict hygiene and biosecurity procedures in breeding farms, as well as a prudent use of antibiotics in fattening farms, can avoid such complex contamination pathways.

**KEYWORDS**

colistin, *mcr-1*, transmission, goat, transposon, IncX4, IncHI2
# 1. Introduction

Colistin is an antibiotic compound that has been widely used in animals during decades whereas it has long been disregarded to treat humans – except as last-resort molecule – due to its nephro- and neurotoxicity. Colistin is active against several bacterial species/ genus among Enterobacterales (mainly *Escherichia coli*, *Salmonella enterica*, *Klebsiella pneumoniae*, or *Enterobacter* spp.) but also non fermenters, such as *Pseudomonas aeruginosa* or *Acinetobacter baumannii*. Colistin-resistance went globally unnoticed until late 2015, when the plasmid-borne *mcr-1* gene was described in *E. coli* and *K. pneumoniae* isolates from food-producing animals, and also from humans, in China. This discovery raised concerns that colistin-resistance in animals could be transmitted to humans and impair last-line treatments in healthcare settings (Liu et al., 2016). Since then, the number of reports on colistin-resistance due to *mcr* genes grew exponentially, mostly in livestock animals, and new *mcr* genes (to date, up to *mcr-10*) and *mcr* gene variants were identified, albeit with very different epidemiological successes (Valiakos and Kapna, 2021). These variants differ by their degree of homology, their frequency (*mcr-1* is reported globally while *mcr-2* has only been identified a couple of times), their geographical distribution and the bacterial host in which they are circulating (*mcr-1* is widespread in Enterobacterales, while *mcr-4* and *mcr-5* are more restricted to *Salmonella* spp.; Hussein et al., 2021).

The *mcr-1* gene is still by far the most widely disseminated gene among the *mcr* family and has principally been associated with IncI2, IncX4, and IncHI2 plasmids worldwide (Matamoros et al., 2017). The selective advantage of these plasmids may rely on increased fitness, together with high transfer efficiency for IncI2 and IncX4 (Liu et al., 2016; Wu et al., 2018), or on capacity to carry multiple antimicrobial resistance (AMR) genes for IncHI2. Altogether, these three *mcr-1*-carrying plasmids have efficiently disseminated in animals since the 2010s, in livestock as well as in pets and wildlife (Migura-Garcia et al., 2020; Moon et al., 2020; Valiakos and Kapna, 2021), but also in humans and the environment (Anyanwu et al., 2020; Gagliotti et al., 2021). In addition to the widely reported plasmid-borne diffusion, the *mcr-1* gene has also been identified on the chromosome, where it is mostly located within the composite transposon Tn6330 (Li et al., 2018).

In livestock in France, the *mcr-1* gene has been detected through various AMR surveillance programs in the pig, poultry and bovine sectors (Haenni et al., 2016; Perrin-Guyomard et al., 2016; Delannoy et al., 2017; Haenni et al., 2018; Hamame et al., 2022) where colistin has long been prescribed to treat diseased animals (ANSES, 2020). In 2016, in the frame of the Resapath network (Mader et al., 2021), we also identified the first *mcr-1*-carrying *E. coli* from a goat. Since small ruminants are usually overlooked in AMR surveys despite a strong usage of antibiotics, including colistin in some breeding practices, this finding prompted us to set up a specific study and search for colistin-resistance in breeding and fattening goat farms. The aim of this study was to determine the proportion of *mcr-1*-positive goat farms and to characterize the collected isolates molecularly, first

to get better insights into the transmission routes of the colistin-resistant clones and plasmids in this food-producing sector, and second to make the breeders aware of their own breeding/fattening practices in order to improve them.

# 2. Materials and methods

## 2.1. Sampling

Between 2018 and 2019, fecal samples were collected from goats in the Poitou-Charentes area, one of the main goat rearing area in France covering 30% and 70% of the total national goat milk herd and milk production, respectively. Two types of farms were sampled. Each farmer provided written consent for sampling and completed a questionnaire on their husbandry and antibiotic use practices. First, 80 goat breeding farms were visited once. In each farm, rectal swabs from 10 adults and 10 kids (<1 month of age) were collected using e-swabs (Biomerieux, Marcy-l'Etoile, France). In a few cases, the expected number of samples could not be obtained, so that a total of 1,561 animals including 775 kids were sampled. Second, in the same period of time, four fattening farms were visited, corresponding to five fattening batches (two successive batches in the same farm). Over the one-month fattening period, four batches were sampled three times (10 goat kids (1-7 days old) upon arrival, 15 and 30 days after arrival) whereas the fifth batch was sampled only twice (upon arrival and 15 days later) due to unexpected earlier slaughtering.

## 2.2. Bacterial isolation and identification

Samples were sent to the laboratory within 48 h after sampling. A total of 50 µl of the fecal suspension were inoculated in an enrichment broth extemporaneously prepared with 9 ml BHI (Brain Heart Infusion) broth supplemented with one colistin disk (10 µg; Mast Diagnostics, Amiens, France). After an overnight incubation at 35 ± 2°C, 50 µl of the broth were streaked onto a ChromID® Colistin agar (Biomerieux). After 18-24 h incubation at 35 ± 2°C, the identification of one pink colony (suspected to be *E. coli*) was confirmed by MALDI-TOF MS (Bruker, Wissembourg, France). All confirmed *E. coli* were kept at -80°C in peptone water plus glycerol (25%) for further investigations. On all isolates, the presence of the *mcr-1* to *mcr-5* genes was tested using a multiplex PCR (Rebelo et al., 2018), and the presence of the *mcr-9* gene was tested using a simplex PCR (Haenni et al., 2020). *mcr-6* to *mcr-8* and *mcr-10*, which are rarer and more restricted to other Enterobacterales species, were not included in the screening.

## 2.3. Antimicrobial susceptibility testing

Susceptibility testing was performed using disc diffusion on Mueller-Hinton agar (BioRad, Marne-la-Coquette, France),
according to the guidelines and clinical breakpoints (human and veterinary) of the Antibiogram Committee of the French Society for Microbiology.<fnref n="1" /> The following discs (Mast Diagnostics) of human and/or veterinary interest were tested: amoxicillin 25 µg, amoxicillin + clavulanic acid 20/10 µg, piperacillin 30 µg, ticarcillin 75 µg, piperacillin + tazobactam 30/6 µg, ticarcillin + clavulanic acid 75/10 µg, cefalotin 30 µg, cefuroxime 30 µg, cefotaxime 30 µg, ceftiofur 30 µg, ceftazidime 30 µg, cefoxitin 30 µg, cefepime 30 µg, cefquinome 30 µg, ertapenem 10 µg, aztreonam 30 µg, streptomycin 10 µg, gentamicin 15 µg, tobramycin 10 µg, netilmicin 30 µg, chloramphenicol 30 µg, florfenicol 30 µg, tetracycline 30 µg, sulfonamides 200 µg, trimethoprim 5 µg, nalidixic acid 30 µg, enrofloxacin 5 µg and ofloxacin 5 µg. Minimum inhibitory concentrations (MICs) were determined by broth microdilution for colistin, according to the European Committee for Antimicrobial Susceptibility Testing (EUCAST). *Escherichia coli* ATCC 25922 was included for quality control both for disc diffusion and broth micro-dilution. All isolates presenting resistance to three or more antibiotic family were considered as Multi-Drug Resistant (MDR; Magiorakos et al., 2012).

## 2.4. Short-read sequencing

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-6,000 instrument (Illumina, San Diego, 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. Quality control statistics of all sequenced isolates are provided as Supplementary Table S1. STs were determined using MLSTFinder v2.0.4 and point mutation in the *pmrAB* genes were identified using PointFinder (CGE online tools, http://www.genomicepidemiology.org/), while resistance genes and replicon content were inferred using ABRicate v1.0.1.<fnref n="2" /> Virulence factors (VFs) were determined using VirulenceFinder and serotypes were determined using SeroType Finder.<fnref n="3" />

## 2.5. Phylogenetic analysis

MLVA characterization was performed prior to sequencing (Caméléna et al., 2019). The cgMLST (core genome multi-locus sequence type) 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 using snp-dists.<fnref n="4" /> The resulting tree was visualized using iTol v.6⁵ (Letunic and Bork, 2019).

## 2.6. Plasmid characterization

The plasmidic content was determined from the WGS data using PlasmidFinder 2.0.1.<fnref n="6" /> Plasmids carrying the *mcr-1* gene were assigned *in silico* when *mcr-1* 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 S1-PFGE 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 for the *mcr-1* gene as well as for the IncX4 and IncH12 replicons, as previously described (Saidani et al., 2019). When plasmidic location could not be evidenced, the chromosomal location was looked for by PFGE on *I-Ceu1*-digested DNA, followed by Southern blot hybridization using a 16S rDNA probe and probes corresponding to the *mcr-1* gene.

## 2.7. Long-read sequencing

To get better insight into a subset of plasmids, MinION long-read sequencing libraries were prepared according to the Oxford Nanopore Technologies, United Kingdom, using the native barcoding expansion kit (catalog number EXP-NBD104; Oxford Nanopore Technologies) and the ligation sequencing kit (SQK-LSK109). Sequencing was performed on a MinION sequencer using a SpotON Mk 1 R9 version flow cell (FLO-MIN106D; Supplementary Table S1). Assembly of both Illumina short reads and Nanopore long reads was performed using Unicycler (Wick et al., 2017). Annotation was done using RASTtk on the Patric 3.25.0 platform<fnref n="7" /> (Brettin et al., 2015) and the genetic environment of the *mcr-1* gene was explored using BLASTn.

## 2.8. Treatment ratios

The treatment ratios was calculated as the number of treatment administered over 12 months divided by the number of breeding animals in the farm.

<footnote n="1">1 www.sfm-microbiologie.org</footnote>

<footnote n="2">2 https://github.com/tseemann/abricate</footnote>

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

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

<footnote n="5">5 http://itol.embl.de/itol.cgi</footnote>

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

<footnote n="7">7 https://www.bv-brc.org/</footnote>
## 2.9. Statistics

The Mann Whitney-Wilcoxon test was used to compare the presence of *mcr-1* and non-*mcr-1* colistin-resistant isolates versus the treatment ratio.

## 2.10. Accession number(s)

All genomic sequences were deposited in DDBJ/EMBL/GenBank under the BioProject accession number PRJNA857909.

# 3. Results

## 3.1. Proportion of colistin-resistance and additional resistance genes

All 80 breeding farms were sampled once (20 animals each) and at least one colistin-resistant (col-R) *E. coli* isolate was identified in 36 of them (36/80, 45.0% of col-R positive farms). Among those 36 col-R positive farms, the *mcr-1* gene was detected in eight of them, leading to an overall 10% (8/80)

proportion of *mcr-1*-positive farms (Supplementary Table S2, Figure 1A). The number of *mcr-1*-positive animals per farm varied from 3 to 15 over 20 goats sampled (Figure 1A). In all, over the 1,561 animals tested, 65 (4.2%) were *mcr-1*-positive, of which 57 were goat kids (57/775, 7.4%). *mcr-1*-positive *E. coli* were thus mostly found in kids and, in a given farm, we never found a *mcr-1*-positive adult alone, in the absence of *mcr-1*-positive kids.

In fattening farms, four over the five batches presented *mcr-1* colistin-resistance at least at one sampling time (Figure 1B), whereas no colistin-resistant isolate was found in the last farm (farm #57) at any of the three sampling times. Similarly, no colistin-resistance was observed in farm #45 at the first sampling time of the second batch; however, colistin-resistant *E. coli* were then identified at the second and third time points. Finally, farm #32 was only sampled twice since goats were slaughtered earlier than expected. The number of *mcr-1*-positive animals per time point varied from 3 to 10 over 10 samples, and *mcr-1*-positive *E. coli* were identified in 60.0% (84/140) of the tested animals.

In both breeding and fattening farms, although we looked for the presence of the *mcr-1* to *mcr-5* genes as well as *mcr-9*, no gene or gene variant other than *mcr-1.1* was identified.

A

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

B

<figure><img src="image_5.png" /><figcaption>FIGURE 1</figcaption></figure>

Schematic representation of *mcr-1*-positive *Escherichia coli* in breeding (A) or fattening (B) farms. The *mcr-1* gene was found on IncX4 (blue colors) or IncHI2 (green colors) plasmids, or on the chromosome (grey colors). (A) Each colored box represents a *mcr-1*-positive goat. All boxes correspond to kids, except those marked by a black dot. Empty boxes indicate that no colistin-resistant isolate found. (B) In ST17 (orange color), *mcr-1* could not be located. Black boxes indicate that the time point was not sampled. Empty boxes indicate that no colistin-resistant isolate found.
## 3.2. Genetic background

All 149 *mcr-1*-positive isolates were characterized by PCR using the published MLVA scheme (Caméléna et al., 2019), and one representative per profile, per farm and per sampling time was further sequenced. In total, 46 *E. coli* were short-read sequenced, among which 25 different STs were identified (Figure 2; Supplementary Table S2).

In the breeding farms, a total of 13 different STs were found (Figures 1A, 2, Supplementary Table S2). The epidemiological success of identified STs varied within the same farm, with some clones identified only once (for example ST215 in farm #12) and others that had spread locally (for example ST591 – identified in 12 of the 15 *mcr-1* positive isolates from farm #47, or ST167 – identified in all 11 isolates from farm #6). Four STs were found in two different breeding farms (ST10, ST57, ST69, ST362). The virulence genes largely varied depending on the identified clones (Supplementary Table S2), but ST301 has particularly caught our attention. This clone, which was found in four isolates of breeding farm #52, was an O80:H2 enterohemorrhagic EHEC possessing, among others, the *stx2*, *eae* and *ehxA* virulence genes (Cointe et al., 2021).

In the fattening farms a total of 18 different STs were found (Figures 1B, 2; Supplementary Table S2). Of note, NT1 was a

single locus variant (SLV) of ST1011, so that the five ST1011 and NT1 isolates only differed by a maximum of 16 SNPs. As two to three sampling campaigns were conducted in each fattening farm, the same ST could be identified over time. When this happened, WGS data proved that this was the same clone in all but one case (2–13 SNPs differences for ST57, ST10, ST744, ST1011, ST3014, non-typable STs NT1 and NT2). The only exception was ST167, which was found in the same farm #45 in 2018 and 2019, but also in two different farms (#32 and #45) 2 months apart. The identified isolates differed, respectively, by 90 and 125 SNPs.

Identical STs were also identified in two to three unrelated breeding or fattening farms on seven other occasions. In five cases (ST10, ST69, ST167, ST362, and ST744), isolates were not genetically related, differing by 29-855 SNPs. In the last two cases, two ST57 (from breeding farm #20 and #37) and two ST1720 (from farm #19 and #32) isolates were identical (0–1 SNP difference). Of note, breeding farm #19 and fattening farm #32 were located on the same geographical site, with specific premises (called breeding farm #19) rearing only adult and kid goats originating from the farm, and other distinct premises (called fattening farm #32) receiving kid goats from farm #19 and other farms of the Poitou-Charentes area.

<figure><img src="image_6.png" /></figure>
### 3.3. Genetic localization of the *mcr-1* gene

Based on the sequenced isolates, the *mcr-1* gene was located either on the chromosome (*n* = 48; 48/149, 32.2%) or on the IncX4 (*n* = 58; 38.9%) and IncHI2 (*n* = 40; 26.8%) plasmids. The *mcr-1* gene could not be reliably located on the chromosome or on a plasmid in three ST17 isolates (fattening farm #32; Supplementary Table S2; Figures 1, 2). In most farms, the *mcr-1* gene was identified exclusively on one genetic determinant.

In breeding farms, the *mcr-1* gene was more frequently found on an IncHI2 plasmid (*n* = 5 farms) than on the chromosome (*n* = 3 farms) or on an IncX4 plasmid (*n* = 1 farm). The two plasmids IncX4 and IncHI2 never cohabitated in a same farm, and the *mcr-1* gene was found located either on the chromosome or on an IncHI2 plasmid only in farm #52.

The situation was similar on the fattening farms, with *mcr-1* located on the chromosome in farm #32, exclusively on an IncX4 plasmid in farm #32 (both in 2018 and 2019), and either on the chromosome or on an IncHI2 plasmid in farm #57.

In four cases (ST167 and ST162 with *mcr-1* located on the chromosome in farms #6 and #32; ST57 with *mcr-1* located on IncHI2 in farms #20 and #37), the *mcr-1* gene had obviously spread in one farm through clonal transmission of one unique clone. On the other hand, plasmid transfer between different *E. coli* isolates was hypothesized for *mcr-1*-bearing IncHI2 plasmids within three farms (#12, #52, #55) and for *mcr-1*-bearing IncX4 within two farms (farm #45, spanning over two batches and five time points, and #47). Long-read sequencing of three isolates (Supplementary Table S2) was thus used to clarify this issue in farm #12, where a *mcr-1*-bearing IncHI2 plasmid was found in three different genetic *E. coli* backgrounds. Sequence analysis showed 99% identity over the 235,257–235,266 bp of the three IncHI2/ST4 plasmids of isolates #49775, #49777, and #49780, strongly arguing for the within-farm spread of a unique *mcr-1*-bearing IncHI2 plasmid between three different *E. coli* isolates. In these IncHI2 plasmids, the *mcr-1* gene alone (without the *pap2* gene) was only preceded by one *ISApl1* element and inserted in a putative kinase. Illumina sequences showed that, in IncX4 plasmids, the *mcr-1-pap2* genes were inserted between hypothetical proteins without any *ISApl1* element (Figure 3).

Chromosomal insertion of *mcr-1* was also observed in several genetic *E. coli* backgrounds in farms #19, #32, and #55. Long-read sequencing of nine isolates (Supplementary Table S2) showed that the *mcr-1* gene was carried on the prototypic Tn6330 transposon additionally presenting the *pap2* gene in seven isolates (Figure 3). Tn6330 was inserted in different locations in the *E. coli* genome, except for isolates #49815 (breeding farm #19) and #49863 (fattening farm #32) which were clonal and presented the Tn6330 insertion at the identical site. Only isolates #49742 (breeding farm #6) and #49870 (fattening farm #32) displayed a particular and identical genetic arrangement, with the absence of the *pap2* gene and the *ISApl1* element. Both isolates belonged to ST167, but differed by 54 SNPs.

### 3.4. Additional resistance genes

While all but one sequenced *E. coli* isolate presented genes conferring resistance to aminopenicillins (TEM-like genes), only two of them (from fattening farm #45) presented the Extended-Spectrum Beta-Lactamase (ESBL) *bla*<sub>CTX-M-15</sub> gene, in coherence with the observed phenotype (Supplementary Table S2). Resistances to tetracyclines (*n* = 45), aminoglycosides (*n* = 43), sulfonamides (*n* = 41), trimethoprim (*n* = 36) and phenicol (*n* = 32) were widespread; the most frequently identified genes were *tet(A)* (*n* = 37), *ant(3’)-Ia* (*n* = 31, conferring resistance to streptomycin/spectinomycin), *sul2* (*n* = 37), *dfrA1* (*n* = 24) and *catA1* (*n* = 19), but several variants of these genes were detected in many isolates (Supplementary Table S2; Figure 2).

In addition to the *mcr-1* gene, seven isolates (7/46, 15.2%) presented mutations in the *pmrA* gene and 33 isolates (71.7%) presented mutations in the *pmrB* gene (Supplementary Table S2). In the *pmrB* gene, the most common mutations were D283G (*n* = 27) and H2R (*n* = 19). Among the mutations that have previously been shown to confer resistance or decreased susceptibility to colistin, the following ones were identified here: G53R and R81S in *pmrA*, each found in one isolate, as well as R93P and P94A in *pmrB*, both found in two isolates.

### 3.5. Antibiotic usage on farms

In France, it is of common knowledge that colistin is widely prescribed to 1-7 days old goat kids at arrival in fattening farms, to avoid digestive disorders related to transportation and rearing poorly immunocompetent young animals of mixed origins. Then, animals may also receive antibiotics to treat inter-current diseases over the fattening period. Here, the different treatments and molecules administered per animal in each farm were unfortunately not available. Nonetheless, available data were the treatment ratios per breeding farm, i.e., the number of treatments administered over 12 months divided by the number of adults in the farm (Figure 4). The vast majority of the farms displayed a ratio between 0 (no treatment) to 1 (mean of one treatment per animal over the breeding period; Figure 4). However, when confronting antibiotic exposure ratios in each farm to the proportion of colistin-resistance or *mcr-1*-positive *E. coli* isolates in the same farm, no statistically relevant correlation (*p* > 0.05) was evidenced.

## 4. Discussion

In this study, we investigated to what extend *mcr* genes conferring colistin resistance may have spread within the caprine sector in France. We thus looked for the presence of *mcr-1* to *mcr-5* genes as well as *mcr-9*. The sole *mcr* gene found here was *mcr-1* (*mcr-1.1* variant), which is among the most frequently reported ones in livestock in Europe. In 34/46 isolates, this *mcr-1*
<figure><img src="image_7.png" /><figcaption>FIGURE 3</figcaption></figure>

Genetic environment of *mcr-1* genes based on Illumina-ONT hybrid assemblies. RM: restriction-modification; hp: hypothetical protein. *Analyses based on Illumina sequences only.

gene was accompanied by point mutations in the *pmrA* and *pmrB* genes. Several of these mutations have been reported in other publications including in colistin-susceptible isolates, so that they most probably play no or little role in colistin-resistance (Kuang et al., 2020; Al-Mir et al., 2021). Among the mutation identified here, only a few had previously been associated with colistin-resistance, namely the G53R and R81S mutation in *pmrA* (Sun et al., 2009; Quesada et al., 2015; Rhouma et al., 2016) as well as the R93P and P94A mutations in *pmrB* (Bourrel et al., 2019; Kuang et al., 2020).

The *mcr-1* gene was detected in 10% of the breeding farms, representing 4.2% of all animals and 7.4% of goat kids tested in those farms, respectively. The early treatment of goat (colistin is mostly prescribed to 1–7 days old kids), together with the maturation of the intestinal microbiota, probably explain the fact

that resistant bacteria are much more frequently identified in young animals. Nevertheless, this 10% proportion of *mcr-1* positive animals is quite low compared to fattening farms where four over the five tested batches were proved to be *mcr-1*-positive, representing 60.0% of the tested animals. Even if comparison between different contexts and sampling protocols can be hazardous, it is undeniable that the proportion of *mcr-1*-positive *E. coli* is more important in fattening farms than in breeding farms. This difference may reflect different antibiotic selective pressures in the two types of farms, a hypothesis that we were unfortunately not capable to confirm since we had no access to individual prescription records. This hypothesis was also not confirmed by more global antibiotic usage estimates based on number of treatments on breeding farms. However, another important driver of colistin resistance in fattening farms may rely
<figure><img src="image_8.png" /></figure>

on *mcr-1* gene horizontal spread, which was supported in this study by identical *mcr-1*-positive clones and/or *mcr-1*-bearing plasmids found in different animals of the same farm.

A major feature of the caprine sector also refers to successive steps in the production chain, where breeding farms are the source of male goat kids in fattening farms. Therefore, *mcr-1*-positive animals in a breeding farm can contaminate several fattening farms, where *mcr-1* prevalence may further amplify through close contacts between individuals. The same snowball effect was observed in France for *E. coli* producing ESBLs in the veal calves sector, where ESBL-positive *E. coli* imported from breeding farms were selected by the antibiotic treatment administered to all animals upon arrival on fattening farms (Gay et al., 2019). These ESBL-positive *E. coli* then tended to disappear over the 5-6 months of fattening; a decrease in *mcr-1*-positive *E. coli* was not observed here, either because the goat fattening period may be too short or because the decrease of *mcr-1* positive isolates would not follow the same dynamic as the decrease of ESBL-positive isolates. In any case, as for ESBL-positive isolates, the vast majority of positive animals were kids, and the proportion of resistance dropped drastically with the passage to adulthood. An issue could come from the persistence of *mcr-1*-positive clones in the environment. Indeed, genomic analyses have shown the presence of identical clones (in the case of ST1011 and NT2) 1 year apart, even though the animal batches had been renewed. Also, some STs - such as ST57 or ST1011 - might be particularly prone to carry the *mcr-1* gene. Both STs have recurrently been reported in the poultry

production (Maciuca et al., 2019; Dandachi et al., 2020; Al-Mir et al., 2021; Jamin et al., 2021; Sadek et al., 2021) but also in humans (Elnahriry et al., 2016; Al-Mir et al., 2019; Hassan et al., 2021).

Most isolates of this study, which were collected from healthy animals, displayed multiple resistance genes (including third and fourth generation cephalosporins) and virulence factors. This indicates on the one hand that these isolates may be co-selected by the use of different classes of antibiotics, and on the other hand that they may become pathogenic either for the animals or for humans if the context allows. In particular, one adult goat and three kids in breeding farm #52 presented an O80:H2 enterohemorrhagic EHEC possessing the *stx2*, *eae* and *ehxA* virulence genes (Cointe et al., 2021). This serotype was among the three major ones causing pediatric haemolytic uraemic syndrome (HUS) in 2015 in France (Bruyand et al., 2019). Knowing that consumption of raw milk and raw milk cheese is one of the risk factors of EHEC transmission, this finding warrants further surveillance.

In this study, isolates presented the *mcr-1* gene either on the chromosome (32.2%) or on the IncX4 (38.9%) and IncH12 (26.8%) plasmids. IncI2, which is more widespread in Asian countries was not identified here (Matamoros et al., 2017).

The IncX4 plasmid, even though it was carried by the largest number of *mcr-1*-positive isolates, was found only in one breeding and one fattening farm. This indicates that IncX4 might not be the most widespread vector of *mcr-1*, but suggests that it is highly successful and persistent once established on a farm. Interestingly,
the only farm in which *mcr-1* was carried on the IncX4 plasmid (breeding farm #47) was selling goat kids to the only fattening farm (#45) in which *mcr-1* was also circulating on the IncX4 plasmids. Even though we did not identify one identical ST between these two farms, the hypothesis of an IncX4/mcr-1 plasmid spread is very likely since this plasmid was found in numerous different genetic backgrounds, showing its high transfer ability. In this farm #45, the first sampling point presented either few (in 2018) or no (in 2019) *mcr-1*-positive *E. coli*, while 90–100% of the sampled animals were positive at the second and third time points, probably indicating an on-farm dynamic of *mcr-1* acquisition over time. We also observed the persistence of the ST1011 *E. coli* over 1 year.

The *mcr-1* gene was carried by the IncHI2 plasmid in five out the eight farms sampled, and in one of the fattening farms (#55). We observed both clonal diffusion and plasmidic transfers. The ST57 was found in two breeding farms (#20 and #37) that apparently share no epidemiological link, and at two time points of the fattening farm #55 which did not receive goat kids from farms #20 and #37. The common source for these three farms, if any, remains unknown. Long-read sequencing proved plasmidic transfer in farm #12 between the ST69, ST744 and ST215 genetic backgrounds. The plasmid was an IncHI2 plasmid of ~235kbp additionally carrying the *bla*~TEM-1~A, *aac*(3)-IIa, *sul*1, *sul*2, *dfr*A1, and *tet*(A) genes conferring resistances to penicillins, sulfonamides-trimethoprim and tetracyclines.

Finally, *mcr-1* was found on the chromosome in 32.2% (n = 48) of the isolates originating from three breeding and two fattening farms. Even though *mcr-1* is usually described as plasmid-borne, this proportion is very close to what has been described in healthy residents in Vietnam (36.8%; Yamaguchi et al., 2020). The presence of a unique *mcr-1*-carrying ST was observed in breeding farms #6 and #52, and in fattening farms #32 (first time point) and #55 (first and third time points). Nevertheless, chromosomal insertions of the *mcr-1* gene were also identified in different genetic backgrounds in breeding farm #19 and fattening farms #32 and #55 (both on the second time point). Among the nine isolates that were long-read sequenced, seven carried the IS*Apl1*-mcr-1-pap2-IS*Apl1* prototypic Tn6330 transposon, which is considered as the principal vehicle of mcr-1 spread, inserted in diverse locations on the *E. coli* chromosome. Several studies suggested that IS*Apl1* enables efficient mcr-1 transposition through the generation of a circular intermediate (Li et al., 2018; He et al., 2019), so that transposon-mediated excision/insertion might be as successful as plasmid transfer. This composite transposon might also lose one or both IS*Apl1* sequences: this was the case in our two last sequenced isolates, which displayed a truncated chromosomally-encoded IS*Apl1*-mcr-1-IS1-like element. This was also the case in the IncHI2 and IncX4plasmid sequenced in this study, which, respectively, presented either the truncated IS*Apl1*-mcr-1 element or only the mcr-1-pap2 genes with no trace of the original IS*Apl1* sequences. Our results reinforce the hypothesis that the loss of the IS*Apl1* sequences might favor the stabilization of the *mcr-1* in a plasmidic backbone (Strepis et al., 2021).

## 5. Conclusion

Our results show that the *mcr-1* gene is actively circulating in goat production, a sector that is often overlooked. The low proportion of *mcr-1*-positive *E. coli* in goat breeding farms observed in this study (4.2% for both adults and kids; 7.4% for goat kids only) appeared sufficient to maintain this resistance gene at a high level (60.0%) in fattening farms. Only strict hygiene and biosecurity procedures in breeding farms, as well as a prudent use of antibiotics in fattening farms, can avoid such a contamination. Between- and within-farm spread of the *mcr-1* gene was also strongly suggested. This spread was due for one third each to the transmission of the IncX4 and IncHI2 plasmids as well as to the Tn6330 transposition in diverse locations of the *E. coli* chromosome, highlighting the fact that chromosomal insertion does not impair the transmission capability of the *mcr-1* gene.

## 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

Written informed consent was obtained from the owners of the animal(s) for the participation in the study.

## Author contributions

MT and MH obtained the EcoAntibio2 grant. MT, J-YM, and MH designed the study. MT conducted the sampling and bacterial isolation. PC and MT performed the bench work. AD performed the WGS analyses. MH drafted the manuscript. All authors contributed to the article and approved the submitted version.

## Funding

This work was supported by internal funding of the French Agency for Food, Environmental and Occupational Health & Safety (ANSES) and by the National EcoAntibio2 Action Plan funded by the Ministry in charge of Agriculture (grant number 2017-389).

## Acknowledgments

We thank all farmers and veterinarians that participated to this study.
## 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.

## Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated

organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Supplementary material

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

## References

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

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

ANSES (2020). *Suivi des ventes de médicaments vétérinaires contenant des antibiotiques en France en 2019*. (Fougères, France: Agence nationale de sécurité sanitaire de l'alimentation, de l'environnement et du travail - Agence nationale du médicament vétérinaire).

Anyanwu, M. U., Jaja, I. F., and Nwobi, O. C. (2020). Occurrence and characteristics of mobile colistin resistance (*mcr*) gene-containing isolates from the environment: a review. *Int. J. Environ. Res. Public Health* 17:1028. doi: 10.3390/ijerph17031028

Bourrel, A. S., Poirel, L., Royer, G., Darty, M., Vuillemin, X., Kieffer, N., et al. (2019). Colistin resistance in Parisian inpatient faecal *Escherichia coli* as the result of two distinct evolutionary pathways. *J. Antimicrob. Chemother.* 74, 1521–1530. doi: 10.1093/jac/dkz090

Brettin, T., Davis, J. J., Disz, T., Edwards, R. A., Gerdes, S., Olsen, G. J., et al. (2015). RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. *Sci. Rep.* 5:8365. doi: 10.1038/srep08365

Bruyand, M., Mariani-Kurkdjian, P., Le Hello, S., King, L. A., Van Cauteren, D., Lefevre, S., et al. (2019). Paediatric haemolytic uraemic syndrome related to Shiga toxin-producing *Escherichia coli*, an overview of 10 years of surveillance in France, 2007 to 2016. *Euro Surveill.* 24:24. doi: 10.2807/1560-7917. ES.2019.24.8.1800068

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–e02818. doi: 10.1128/AEM.02812-18

Cointe, A., Bizot, E., Delannoy, S., Fach, P., Bidet, P., Birgy, A., et al. (2021). Emergence of new ST301 Shiga toxin-producing *Escherichia coli* clones harboring extra-intestinal virulence traits in Europe. *Toxins (Basel)* 13:686. doi: 10.3390/toxins13100686

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

Delannoy, S., Le Devendec, L., Jouy, E., Fach, P., Drider, D., and Kempf, I. (2017). Characterization of colistin-resistant *Escherichia coli* isolated from diseased pigs in France. *Front. Microbiol.* 8:2278. doi: 10.3389/fmicb.2017.02278

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

Gagliotti, C., Bolzoni, L., Carretto, E., Sarti, M., Ricchizzi, E., Ambretti, S., et al. (2021). Reduction trend of *mcr-1* circulation in Emilia-Romagna region, Italy. *Eur. J. Clin. Microbiol. Infect. Dis.* 40, 2585–2592. doi: 10.1007/s10096-021-04318-y

Gay, E., Bour, M., Cazeau, G., Jarrige, N., Martineau, C., Madec, J.-Y., et al. (2019). Antimicrobial usages and antimicrobial resistance in commensal *Escherichia coli* from veal calves in France: evolution during the fattening process. *Front. Microbiol.* 10:792. doi: 10.3389/fmicb.2019.00792

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 *bla*<sub>CTX-M-55</sub> in cattle in France. *J. Antimicrob. Chemother.* 73, 533–536. doi: 10.1093/jac/dkx418

Haenni, M., Metayer, V., Gay, E., and Madec, J. Y. (2016). Increasing trends in *mcr-1* prevalence among extended-spectrum-beta-lactamase-producing *Escherichia coli* isolates from French calves despite decreasing exposure to colistin. *Antimicrob. Agents Chemother.* 60, 6433–6434. doi: 10.1128/AAC.01147-16

Haenni, M., Métayer, V., Jarry, R., Drapeau, A., Puech, M. P., Madec, J. Y., et al. (2020). Wide spread of *bla*<sub>CTX-M-9</sub>/<i>mcr-9</i> IncH12/ST1 plasmids and CTX-M-9-producing *Escherichia coli* and *Enterobacter cloacae* in rescued wild animals. *Front. Microbiol.* 11:601317. doi: 10.3389/fmicb.2020.601317

Hamame, A., Davoust, B., Rolain, J.-M., and Diene, S. M. (2022). Genomic characterisation of an *mcr-1* and *mcr-3*-producing *Escherichia coli* strain isolated from pigs in France. *J. Global Antimicrob. Resist.* 28, 174–179. doi: 10.1016/j.jgar.2022.01.014

Hassan, J., Mann, D., Li, S., Deng, X., and Kassem, I. I. (2021). Emergence of the mobile colistin resistance gene *mcr-1* in multidrug-resistant *Escherichia coli* isolated from the fecal matter of toddlers in a community. *Antimicrob. Agents Chemother.* 65, e02434–e02211. doi: 10.1128/AAC.00243-21

He, Y.-Z., Li, X.-P., Miao, Y.-Y., Lin, J., Sun, R.-Y., Wang, X.-P., et al. (2019). The *ISAPl2* dimer circular intermediate participates in *mcr-1* transposition. *Front. Microbiol.* 10:15. doi: 10.3389/fmicb.2019.00015

Hussein, N. H., Al-Kadmy, I. M. S., Taha, B. M., and Hussein, J. D. (2021). Mobilized colistin resistance (*mcr*) genes from 1 to 10: a comprehensive review. *Mol. Biol. Rep.* 48, 2897–2907. doi: 10.1007/s11033-021-06307-y

Jamin, C., Sanders, B. K., Zhou, M., Costessi, A., Duijsings, D., Kluytmans, J. A. J. W., et al. (2021). Genetic analysis of plasmid-encoded *mcr-1* resistance in Enterobacteriaceae derived from poultry meat in the Netherlands. *JAC- Antimicrob. Resist.* 3:dlab156. doi: 10.1093/jacamr/dlab156

Kuang, Q., He, D., Sun, H., Hu, H., Li, F., Li, W., et al. (2020). R93P substitution in the PmrB HAMP domain contributes to colistin heteroresistance in *Escherichia coli* isolates from swine. *Antimicrob. Agents Chemother.* 64, e01509–e01520. doi: 10.1128/AAC.01509-20

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

Li, R., Yu, H., Xie, M., Chen, K., Dong, N., Lin, D., et al. (2018). Genetic basis of chromosomally-encoded *mcr-1* gene. *Int. J. Antimicrob. Agents* 51, 578–585. doi: 10.1016/j.ijantimicag.2017.11.015

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

Maciuca, 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

Mader, R., Jarrige, N., Haenni, M., Bourély, C., Madec, J. Y., and Amat, J. P. (2021). OASIS evaluation of the French surveillance network for antimicrobial resistance in diseased animals (RESAPATH): success factors underpinning a well-performing voluntary system. *Epidemiol. Infect.* 149:e104. doi: 10.1017/S0950268821000856

Magiorakos, A. P., Srinivasan, A., Carey, R. B., Carmeli, Y., Falagas, M. E., Giske, C. G., et al. (2012). Multidrug-resistant, extensively drug-resistant and
pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. *Clin. Microbiol. Infect.* 18, 268–281. doi: 10.1111/j.1469-0691.2011.03570.x

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. doi: 10.1038/s41598-017-15539-7

Migura-Garcia, L., González-López, J. J., Martinez-Urtaza, J., Aguirre Sánchez, J. R., Moreno-Mingorance, A., Perez De Rozas, A., et al. (2020). *Mcr*-colistin resistance genes mobilized by IncX4, IncH12, and IncI2 plasmids in *Escherichia coli* of pigs and white stork in Spain. *Front. Microbiol.* 10:3072. doi: 10.3389/fmicb.2019.03072

Moon, D. C., Mechesso, A. F., Kang, H. Y., Kim, S. J., Choi, J. H., Kim, M. H., et al. (2020). First report of an *Escherichia coli* strain carrying the colistin resistance determinant *mcr-1* from a dog in South Korea. *Antibiotics (Basel)* 9:768. doi: 10.3390/antibiotics9110768

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

Perrin-Guyomard, A., Bruneau, M., Houee, P., Deleurme, K., Legrandois, P., Poirier, C., et al. (2016). Prevalence of *mcr-1* in commensal *Escherichia coli* from French livestock, 2007 to 2014. *Euro Surveill.* 21:pii=30135. doi: 10.2807/1560-7917. ES.2016.21.6.30135

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. doi: 10.1093/jac/dku320

Rebelo, A. R., Bortolaia, V., Kjeldgaard, J. S., Pedersen, S. K., Leekitcharoenphon, P., Hansen, I. M., et al. (2018). Multiplex PCR for detection of plasmid-mediated colistin resistance determinants, *mcr-1*, *mcr-2*, *mcr-3*, *mcr-4* and *mcr-5* for surveillance purposes. *Euro Surveill.* 23:pii=17-00672. doi: 10.2807/1560-7917.ES.2018.23.6.17-00672

Rhouma, M., Beaudry, F., Thériault, W., and Letellier, A. (2016). Colistin in pig production: chemistry, mechanism of antibacterial action, microbial resistance

emergence, and one health perspectives. *Front. Microbiol.* 7:1789. doi: 10.3389/fmicb.2016.01789

Sadek, M., Ortiz De La Rosa, J. M., Abdelfattah Maky, M., Koraše Dandrawy, M., Nordmann, P., and Poirel, L. (2021). Genomic features of MCR-1 and extended-spectrum β-lactamase-producing Enterobacterales from retail raw chicken in Egypt. *Microorganisms* 9:195. doi: 10.3390/microorganisms9010195

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. *Microb. 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

Strepis, N., Voor Int Holt, A. F., Vos, M. C., Zandijk, W. H. A., Heikema, A. P., Hays, J. P., et al. (2021). Genetic analysis of mcr-1-carrying plasmids from gram-negative bacteria in a Dutch tertiary care hospital: evidence for intrapatient and interspecies transmission events. *Front. Microbiol.* 12:727435. doi: 10.3389/fmicb.2021.727435

Sun, S., Negrea, A., Rhen, M., and Andersson, D. I. (2009). Genetic analysis of colistin resistance in *salmonella enterica* Serovar typhimurium. *Antimicrob. Agents Chemother.* 53, 2298–2305. doi: 10.1128/AAC.01016-08

Valiakos, G., and Kapna, I. (2021). Colistin resistant *mcr* genes prevalence in livestock animals (swine, bovine, poultry) from a multinational perspective. A systematic review. *Vet. Sci.* 8:8. doi: 10.3390/vetsci8110265

Wick, R. R., Judd, L. M., Gorrie, C. L., and Holt, K. E. (2017). Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. *PLoS Comput. Biol.* 13:e1005595. doi: 10.1371/journal.pcbi.1005595

Wu, R., Yi, L. X., Yu, L. F., Wang, J., Liu, Y., Chen, X., et al. (2018). Fitness advantage of *mcr-1*-bearing IncI2 and IncX4 plasmids *in vitro*. *Front. Microbiol.* 9:331. doi: 10.3389/fmicb.2018.00331

Yamaguchi, T., Kawahara, R., Hamamoto, K., Hirai, I., Khong, D. T., Nguyen, T. N., et al. (2020). High prevalence of colistin-resistant *Escherichia coli* with chromosomally carried *mcr-1* in healthy residents in Vietnam. *mSphere* 5, e00117–e00120. doi: 10.1128/mSphere.00117-20