Molecular Epidemiology of mcr-Encoded Colistin Resistance in Enterobacteriaceae From Food-Producing Animals in Italy Revealed Through the EU Harmonized Antimicrobial Resistance Monitoring

Colistin resistance by mobilisable mcr genes has been described in bacteria of food-animal origin worldwide, which has raised public health concerns about its potential foodborne transmission to human pathogenic bacteria. Here we provide baseline information on the molecular epidemiology of colistin-resistant, mcr-positive Escherichia coli and Salmonella isolates in food-producing animals in Italy in 2014-2015. A total 678, 861 and 236 indicator E. coli, Extended Spectrum Beta-Lactamase (ESBL)/AmpC-producing E. coli, and Salmonella isolates, respectively, were tested for colistin susceptibility. These isolates were collected according to the EU harmonized antimicrobial resistance monitoring program and are representative of at least 90 and 80% of the Italian poultry (broiler chickens and turkeys) and livestock (pigs and bovines < 12 months) production, respectively. Whole genome sequencing by Illumina technology and bioinformatics (Center for Genomic Epidemiology pipeline) were used to type 42 mcr-positive isolates by PCR. Colistin resistance was mainly observed in the ESBL/AmpC E. coli population, and was present in 25.9, 5.3, 6.5, and 3.9% of such isolates in turkeys, broilers, pigs, and bovines, respectively. Most colistin-resistant isolates (141/161, 87.5%) harbored genes of the mcr-1 group. mcr-1 was also detected in a small proportion of Salmonella isolates (3/146, 2.0%) in turkeys. Additional mcr types were mcr-3 in four ESBL-producing E. coli from bovines, and two mcr-4 in ESBL (n = 1) and indicator E. coli (n = 1) from pigs and bovines. We describe notable diversity of mcr variants with predominance of mcr-1.1 and mcr-1.2 on conjugative IncX4 plasmids in E. coli and in Salmonella serovars Typhimurium, Newport, Blockley from turkey. A new variant, mcr-1.13 was detected in the chromosome in E. coli in turkey and pig isolates. Additionally, we describe mcr-3.2 and mcr-4.3 in E. coli from bovines, and mcr-4.2 in E. coli from pigs. These findings elucidate the epidemiology of colistin resistance in food-producing animals in Italy along with its genetic background, and highlight the likelihood of mcr horizontal transfer between commensal bacteria and major food-borne pathogens (Salmonella) within the same type of productions. Thorough action and strategies are needed in order to mitigate the risk of mcr transfer to humans, in a “One Health” perspective.

Colistin resistance by mobilisable mcr genes has been described in bacteria of food-animal origin worldwide, which has raised public health concerns about its potential foodborne transmission to human pathogenic bacteria. Here we provide baseline information on the molecular epidemiology of colistin-resistant, mcr-positive Escherichia coli and Salmonella isolates in food-producing animals in Italy in 2014-2015. A total 678, 861 and 236 indicator E. coli, Extended Spectrum Beta-Lactamase (ESBL)/AmpC-producing E. coli, and Salmonella isolates, respectively, were tested for colistin susceptibility. These isolates were collected according to the EU harmonized antimicrobial resistance monitoring program and are representative of at least 90 and 80% of the Italian poultry (broiler chickens and turkeys) and livestock (pigs and bovines < 12 months) production, respectively. Whole genome sequencing by Illumina technology and bioinformatics (Center for Genomic Epidemiology pipeline) were used to type 42 mcr-positive isolates by PCR. Colistin resistance was mainly observed in the ESBL/AmpC E. coli population, and was present in 25.9, 5.3, 6.5, and 3.9% of such isolates in turkeys, broilers, pigs, and bovines, respectively. Most colistin-resistant isolates (141/161, 87.5%) harbored genes of the mcr-1 group. mcr-1 was also detected in a small proportion of Salmonella isolates (3/146, 2.0%) in turkeys. Additional mcr types were mcr-3 in four ESBL-producing E. coli from bovines, and two mcr-4 in ESBL (n = 1) and indicator E. coli (n = 1) from pigs and bovines. We describe notable diversity of mcr variants with predominance of mcr-1.1 and mcr-1.2 on conjugative IncX4 plasmids in E. coli and in Salmonella serovars Typhimurium, Newport, Blockley from turkey. A new variant, mcr-1.13 was detected in the chromosome in E. coli in turkey and pig isolates. Additionally, we describe mcr-3.2 and mcr-4.3 in E. coli from

INTRODUCTION
Colistin is a polymyxin classified among the Highest Priority Critically Important Antimicrobials for human medicine by the World Health Organization (WHO), and it is considered a last resort antimicrobial for the treatment of infections by carbapenem-resistant Enterobacteriaceae in humans (Poirel et al., 2017).
In 2015, colistin resistance mediated by mcr-1, a phosphoethanolamine transferase gene located on a transferable plasmid was first reported in Escherichia coli from animals, food and patients from China (Liu et al., 2016). Since then, at least 32 countries from the five continents have found mcr-1 in E. coli isolates from different sources including humans, animals and foods (Xavier et al., 2016).
In Europe, the presence of mcr-1 was first detected in E. coli from poultry meat and humans in Denmark (Hasman et al., 2015). Subsequently, this gene was found in Enterobacteriaceae from different sources in almost all European countries, including Italy (Battisti, 2016;Cannatelli et al., 2016). In 2016, a new mcr gene, mcr-2, was described in E. coli in calves and piglets in Belgium (Xavier et al., 2016), followed by the description of three additional mobile colistin resistance genes in 2017, namely mcr-3 (Yin et al., 2017) in E. coli isolated from pig in China, mcr-4 (Carattoli et al., 2017) in Salmonella and E. coli from pigs in Italy, Spain and Belgium, and mcr-5 (Borowiak et al., 2017) in Salmonella Paratyphi B from poultry and environmental sources from Germany.
It is clear that the epidemiology of transferable mcr-mediated colistin resistance is evolving rapidly and timely information on prevalence and molecular epidemiology of mcr-positive isolates is needed to enhance surveillance, and implement measures to prevention and to control further spread of colistin resistance. In the European Union (EU), the harmonized antimicrobial resistance (AMR) monitoring and reporting in poultry and livestock, which includes colistin susceptibility testing in E. coli and Salmonella, ensures that prevalence of colistin-resistant bacteria in a representative proportion of the food-animal population is reported from each Member State (MS) (Decision 2013/652/EU). However, the lack of molecular data limits the epidemiologic monitoring of colistin resistance (Schrijver et al., 2017).
The aim of this study is to determine the prevalence of colistin resistance, and the molecular epidemiology of mcrmediated colistin resistance genes and their genetic environment in commensal E. coli, Extended Spectrum Beta-Lactamase (ESBL)/AmpC-producing E. coli, and Salmonella spp. in foodanimals in Italy in 2014-2015.

Study Design, Sample Collection, Isolation, and Identification of Bacterial Cultures
Study design and sampling were performed according to Decision 2013/652/EU 1 , which mandates each EU Member State (MS) to collect caecal content samples from different epidemiological units of poultry flocks (broiler chickens, fattening turkeys), fattening pigs and bovines < 12 months at slaughter.
Samples were collected from broiler chickens (n = 300) and fattening turkeys (n = 300) in 2014, and from fattening pigs (n = 304) and bovines < 12 months (n = 223) in 2015 ( Table 1). The regional stratification of samples represented at least 90 and 80% of the Italian poultry (broiler chickens and turkeys) and livestock (pigs and bovines < 12 months) production, respectively.
In addition, 558 and 709 samples from fattening turkey and broiler chicken flocks, respectively, were collected within the voluntary national Salmonella monitoring framework (Decision 2013/652/EU) in 2014 ( Table 1). In compliance with Decision 2013/652/EU, E. coli isolation and identification, were performed according to the EURL-AR protocols 2 , whereas Salmonella spp., isolation, identification and serotyping were performed according to the ISO 6579:2002/Amd 1:2017 protocols.

Conjugation Experiments
Four representative E. coli and three Salmonella sp. isolates from turkeys were selected as donors. Conjugation experiments were performed as previously described (Franco et al., 2015), with the only modification regarding the MacConkey agar plates selective for transconjugants which contained 2 mg/L colistin sulfate and 100 mg/L rifampicin in this study.

Whole Genome Sequencing (WGS) and Bioinformatics Analysis
A total 42 isolates which tested mcr-positive by PCR, (28 E. coli and three Salmonella enterica isolates from turkeys, five E. coli from pigs and six E. coli from cattle), and the seven E. coli K-12 isolates result of the conjugation experiments were Whole Genome Sequenced. Genomic DNA was extracted using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) following the manufacturer's protocol. Libraries were prepared for Illumina pair-end sequencing using the Illumina (Illumina, Inc., San Diego, CA, United States) NexteraXT R Guide 150319425031942. Sequencing was performed using an Illumina platform (MiSeq or HiSeq2000). Raw sequence data were submitted to the European Nucleotide Archive 3 under study accession no.: PRJEB23728, PRJEB23778, PRJEB21546, and PRJEB26479.
Raw reads were assembled and analyzed using the pipeline from the Center for Genomic Epidemiology (CGE 4 , Thomsen et al., 2016), with default settings. This pipeline performed de novo assembly (Velvet based), species identification (KmerFinder 2.1), Multilocus Sequence Typing (MLST 1.6), identification of virulence (VirulenceFinder 1.2) and antimicrobial resistance genes (ResFinder 2.1), identification of plasmid incompatibility groups (PlasmidFinder 1.2) and plasmid MLST (pMLST 1.4). When identity values for mcr were less than 100% in the ResFinder output, the sequence was submitted to online BLAST 5 (Zhang et al., 2000) to identify the exact mcr variant.
Manual annotation of the contigs containing selected mcr variants was performed using RAST (Aziz et al., 2008), BLAST (Zhang et al., 2000) and ISfinder (Siguier et al., 2006). The contigs harboring the new mcr variants were compared with reference sequences using BLAST and EasyFig (Sullivan et al., 2011). The references sequences used for comparison were CP016034, KP347127 and KY924928.
Colistin Resistance and mcr in E. coli From Fattening Pigs and Bovines <12

Genomic Characterisation
High diversity of Sequences Types (STs) was evident when studying the genome of the 28 E. coli and three Salmonella enterica isolates from turkeys, five E. coli from pigs and six E. coli from bovine <12 months selected for WGS. In E. coli from turkeys, 20 STs were identified. The most represented STs were ST-155 and ST-156 with four isolates each, followed by ST-744 with three isolates and ST-101 with two isolates (Supplementary Table 1 Table 1).
mcr-1.1 (NG_050417.1) was present in 24 isolates: 18 E. coli and two Salmonella isolates from turkeys, three E. coli from pigs and one from bovines. mcr-1.2 (KX236309.1) was detected in eight E. coli and one Salmonella from turkeys. In addition, one E. coli isolate from turkey presented a silent mutation, C801T, and two E. coli isolates, one from turkey and one from fattening pigs showed a new variant (hereafter termed mcr-1.13) of mcr-1.1 with two non-synonymous mutations: M2V and S14G. Four E. coli from bovines <12 months presented a mcr-3 variant, identified as mcr-3.2 (NG_055523.1), and one of them presented both mcr-1.1 and mcr-3.2 (Supplementary Table 1). Also, one ESBL(CTX-M-1)-producing E. coli isolate from pigs and one indicator E. coli from bovines presented mcr-4.2 and mcr-4.3, respectively.
The two E. coli isolates harboring mcr-1.13 allele, ID:14077295 from turkey and ID:15056414 from fattening pig, were genetically different, as shown by MLST (ST-69; ST-5995) and different resistance gene content (Supplementary Table 1). The virulence gene gad was detected in both isolates, but it was the only virulence-associated gene found in the pig isolate (Supplementary Table 1 Table 1). All four were MDR and presented the ESBL gene bla CTX−M−55 , although always found in contigs other than those containing mcr-3.2. The four isolates also shared other resistance accessory genes as sul3, aac(3)-Iid, aadA2, or floR. A number of different plasmids were detected in all four isolates, with only IncF (F46:A-:B20), IncX1 and IncR detected in all of them. The contigs harboring mcr-3.2 did not contain any plasmid replicon type (Supplementary Table 1).
Description of the Contig Harboring mcr-1.13 mcr-1.13 was identified in E. coli strain 14077295 (contig 564; 10,310 bp lenght) from turkeys and in E. coli strain 15056414 (contig 27; 50,840 bp lenght) from pigs. The BLAST alignment of the two contigs showed 99% identity across the entire length of the shortest contig (10,310 bp). No plasmid replicon was detected in any of the two contigs. mcr-1.13 (from nt 33,836 to nt 35,462 in contig 27,) was upstream to a PAP2 superfamily hypothetical protein (from nt 35,509 to nt 36,256). Insertion sequences (IS) were detected flanking the mcr cassette: upstream mcr-1.13 a truncated IS66 (from nt 33,335 to 33,413 bp and from 33,455 to 33,651 bp) and downstream a small fragment of IS66 (from nt 36,261 to nt 36,392) and IS110 (from nt 36,449 to nt 36,816). Upstream the cassette, there was the gene coding for 50S ribosome-binding GTPase family protein (from nt 15,712 to nt 16,627) and downsteam the cassette, there was the coding gene of the subunit YeeA of the methylase (from nt 37,032 to nt 39,849) (Figure 1). In BLAST results, contig sequences showed 93-99% identity to E. coli chromosome sequence CP016034 except for the region flanked by ISs and containing mcr-1.13 and PAP2 superfamily coding genes that had 99 and 100% identity with the same region of the mcr plasmid pHNSHP45(KP347127), for the mcr gene and the PAP2 superfamily coding gene, respectively (Figure 1).

Plasmid Transferability
E. coli harboring mcr-1.1 (IDs:14043377 and 14065450), mcr-1.2 (ID:14047606) and mcr-1.13 (ID:14077295), S. Typhimurium harboring mcr-1.1 (ID:14043372), S. Newport harboring mcr-1.1 FIGURE 1 | Graphical representation of the mcr-1.13 contig. Graphical representation of the BLAST analysis between the contig harboring mcr-1.13 from the pig-origin E. coli isolate (contig27; 15056414) and the mcr plasmid pHNSHP45 (KP347127; region 21000bp-26000bp) and genomic DNA of E. coli Co6114 (CP016034; region 91844bp -99527bp). * Indicated that the gene was manually annotated. The gray area represents the blast identities, the percentage of identity is indicated in the legend. Gene colors: red, mcr-1; blue, transposases or IS elements. FIGURE 2 | Graphical representation of the mcr-3.2 contig. Graphical representation of the BLAST analysis between the contig harboring mcr-3.2 in the bovine-origin E. coli (contig77; 15054212) and the mcr plasmid pWJ1 (KY924928; region 160180 to 163525 bp). * Indicated that the gene was manually annotated. The gray area represents the blast identities, the percentage of identity is indicated in the legend. Gene colors: red, mcr-1; blue, transposases or IS elements.
(ID:14038647) and S. Blockley harboring mcr-1.2 (ID:14085183) were used as donors for conjugation experiments ( Table 2  and Supplementary Table 1). All donors, except the mcr-1.13harboring E. coli, transferred mcr by conjugation to the E. coli K-12 recipient strain. mcr-1.1 and mcr-1.2 were carried by IncX4 plasmids that transferred either alone or in combination with additional plasmids ( Table 2 and Supplementary Table 1). One isolate (ID: 14065450) displayed in the same contig both mcr-1.1 and the replicon of the IncX4 plasmid (contig 68).

DISCUSSION
The high diversity of transferable colistin resistance mediated by mcr genes and alleles, quickly spreading across pathogenic enterobacteria globally (Kluytmans, 2017), is an emerging challenge for treatment of Gram-negative infections due to increased occurrence of Healthcare-Associated Extremely Drug Resistant bacterial pathogens (EARS-Net, 2015).
In the present study, we found high prevalence (∼25%) of colistin resistance in both indicator commensal E. coli and ESBL/AmpC-producing E. coli in turkeys in Italy. In other Italian primary productions such as broilers, fattening pigs and bovines <12 months, colistin-resistant E. coli occurred at relatively lower, though still relevant levels.
From a (Veterinary) Public Health perspective, the spread of colistin resistance across production types is of major concern also considering that it was often associated with resistance to multiple drugs and notably also to fluoroquinolones. For instance, MDR E. coli exhibiting co-resistance to three Critically Important Antimicrobials for humans were detected in approximately 23% of turkey flocks investigated.
By using WGS, the vast diversity of mcr types and variants was evident within and across the major production types. At least three mcr-1 alleles were identified in turkeys, pigs and bovines <12 months in a 2-year-period, one of which has not been previously described. The majority of isolates harbored mcr-1.1, but mcr-1.2, mcr-1.13, mcr-3.2, mcr-4.2 and mcr-4.3 were also detected.
The core genomes of the E. coli and Salmonella isolates harboring mcr-1 were quite variable as shown by diversity of MLST types, in agreement with what has been observed in different contexts previously (Hasman et al., 2015;Bernasconi et al., 2016;Cannatelli et al., 2016;Zhao et al., 2017).
Conjugation experiments and WGS data analysis confirmed that mcr-1.1 and mcr-1.2 were located on IncX4 plasmids in E. coli and Salmonella spp. from turkey. In Europe, IncX4 plasmids have been described in association with mcr-1 in isolates of animal origin in Belgium (Xavier et al., 2016;Garcia-Graells et al., 2017) and from humans in Italy (Di Pilato et al., 2016).
The mcr-1.2 gene was first described on an IncX4 plasmid in a Klebsiella pneumoniae isolated from a human sample in the Tuscany Region (Italy) (Di Pilato et al., 2016). In our study, this variant was found on IncX4 plasmids in intestinal Salmonella and E. coli spread throughout the Italian fattening turkey production. Our observations suggest the transmission of mcr-positive IncX4 plasmids between different bacterial species, with the possibility of transmission from animals to humans, or vice versa.
The mcr-3 gene was first described on an IncHI plasmid (Yin et al., 2017), whereas the mcr-3.2 variant has been described from the whole genome shotgun sequence of a Shigella sonnei strain (NG_055523; NPZH01000177) isolated in the United States, but to our understanding, no formal description has been published to date. In our study, four ESBL/AmpC-producing E. coli isolates from bovines <12 months presented mcr-3.2. The WGS data analysis performed on these isolates did not allow us to determine the genomic location of mcr-3.2. However, only one of these isolates harbored an IncHI plasmid. A ST heterogeneity of E. coli isolates harboring the mcr-3.2 gene was evident. Further conjugation or transformation experiments are needed in order to elucidate if the mcr-3.2 found in E. coli in Italian cattle is plasmid-borne and transferable.
One E. coli from cattle was simultaneously positive for mcr-3.2 and mcr-1.1. In this isolate, mcr-1.1, was located on an IncX4 plasmid. An E. coli carrying both mcr-1.1 and mcr-3.2 and also isolated from cattle has been recently described in Spain (Hernández et al., 2017). Interestingly, all four mcr-3.2-positive E. coli isolated detected in bovines <12 months in Italy carried the ESBL bla CTX−M−55 gene, apparently an emerging variant in veal calves in Italy. The same observation has been made by other authors (Hernández et al., 2017;Roer et al., 2017). Taken together, these features could be suggestive of a genetic linkage between the two genes. Our study, however, provides preliminary data that the mcr-3.2 and the bla CTX−M−55 genes may not be linked on the same genetic element, and that this aspect needs further investigation.
The mcr-4.2 gene was described recently from two S. enterica Typhimurium (monophasic variant) isolated from human samples in Italy (Carretto et al., 2018). So far, this variant has not been associated to any plasmid. In our study, the mcr-4.2 gene identified in a ESBL(CTX-M-1)-producing E. coli isolated from a fattening pig, was located on a ColE plasmid, exactly the same type of plasmid in which mcr-4.1 was described by Carattoli et al., 2017. The mcr-4.3 variant was first described from a Salmonella Kedougou isolated from pigs in Spain (Rebelo et al., 2018). In our study, mcr-4.3 was harbored by an indicator E. coli from a bovine caecal sample. All the mcr-4 variants described so far have been found in Salmonella or E. coli from animals and from humans, despite their localisation on non-conjugative plasmids, similarly to what has been found with mcr-4.1 by Carattoli et al. (2017).
The mcr-1.13 gene was located on a mobile genetic element inserted in the chromosome. The mcr cassette found is different from the chromosomal mcr cassettes described so far, because of the presence of ISs other than ISplA (Snesrud et al., 2016;Li et al., 2017). By considering the high rate of self-transferability of the IncX plasmids (Sun et al., 2017) and that IncX4 is dominant in the Italian primary productions surveyed, this finding suggests that the mcr cassette may have been acquired from an IncX4 plasmid which lacks ISplA insertion sequence (Bernasconi et al., 2016;Matamoros et al., 2017;Sun et al., 2017). It is likely that the insertion occurred in different moments, since the two E. coli isolates derived from two different animal productions and presented different genome content. The ability of the mcr cassette to jump into several types of plasmids (IncI2, IncX4, IncHI2, and IncP) or into the chromosome is highly concerning also in view of the possibility of insertion into plasmids already present in MDR isolates, or even of the creation of a new megaplasmid as occurred in Salmonella Infantis (Franco et al., 2015). This would imply that colistin resistance could be co-selected by the use of others antimicrobials. Indeed, multidrug resistance in colistin-resistant, mcr-positive isolates is a constant feature we have observed in these population-based studies. These findings highlight how rapidly mobile genetic elements can be acquired by Enterobacteriaceae and how genes can mutate in environments with high selective pressure as that occurring in the intensive farming systems. In this respect, the selection pressure exerted by the use of colistin and the complex co-selection mechanisms triggered by the overall high use of several antibiotic classes in the meat-producing industry is likely to have played a major role. For instance, the overall exposure to antibiotics in food-producing animals in Italy in 2014 was estimated around 360 mg/Population Correction Unit (PCU, i.e., mg per kilogram of biomass of farmed animals), and the exposure to colistin only was estimated around 29 mg/PCU (ESVAC, 2016). The mean, median and range of total sales in Europe (29 countries) were 108, 66, and 3 -419 mg/PCU, respectively, while 22 out of 29 countries reported sales of colistin ≤5 mg/PCU (range of the 29 EU countries: 0.06-36.10 mg/PCU). It is well known that the intensively farmed, meat-producing animals (turkeys, broilers, pigs, bovines <12 months) are the population categories at risk of exposure to higher amounts of colistin in modern farming systems, and the results of this study strongly support the general concept that the spread of mcr-mediated colistin resistance has been favored by the semicontinuous and high exposure to colistin in these production chains.
Whether animals are an important source for human extraintestinal pathogenic E. coli (ExPEC) infections in humans is still a matter of debate (Bélanger et al., 2011). Few studies provide information and comparison on STs and plasmids harboring ESBLs or transferable AmpCs (Leverstein-van Hall et al., 2011) for both human and food-producing animal isolates. Although conclusive epidemiological evidence is still lacking, it has been proposed that some human ExPEC infections could arise from poultry and pig ExPEC reservoirs (Jakobsen et al., 2010). That being said, it is interesting to notice that at least one-fourth of E. coli isolates described in our study belong to the same STs as isolates associated with human ExPEC infections in Europe. In some cases (Leverstein-van Hall et al., 2011;Mavroidi et al., 2012;Brolund et al., 2014) they also share the same ESBL (e.g., ST69 and CTX-M-14, ST410 and CTX-M-1). As for colistinresistant E. coli, one of the mcr-1.1 positive isolates from turkeys here described belongs to ST131, although it lacks the virulence gene markers and ESBL generally associated with the globally spread human clinical clone. Indeed, a mcr-1-positive, ESBLnegative ST-131 E. coli was also described in human bloodstream infections in Italy in 2017 (Corbella et al., 2017). Overall, these observations should be taken very cautiously, since methods for genome analysis and parameters for assessing relatedness among both core genomes and accessory genomes are quickly evolving, which implies that it may be necessary to re-evaluate any earlier conclusions on relatedness or source attribution based on partial molecular characteristics, as previously shown (De Been et al., 2014).

CONCLUSION
In conclusion, harmonized cross sectional studies at slaughter like the ones implemented by the EU represent a very important tool for a deep insight into trends and emergence of antimicrobial resistance traits and patterns in major food-borne pathogens and commensal opportunistic bacteria. Especially when occurring at high prevalence, the spread of transferable colistin resistance in E. coli (both indicator commensal and ESBL/AmpC-producing isolates) is to be considered a concern per se. Additionally, as a general principle, the high spread of resistance increases the probability of transfer of specific resistance traits also to major zoonotic pathogens, such as Salmonella spp. The hypothesis that horizontal transfer, so far, has played a major role in spread of colistin resistance among bacteria in Italian meat-producing animals is supported by the observed heterogeneity of mcrpositive E. coli. Indeed, at least in the Italian turkey productions, we demonstrated that the same transferable determinant of colistin-resistance is being carried on the same conjugative plasmid in both E. coli and major Salmonella serotypes detected in the same intensive-farming industry.
For the above reasons, quick and thorough action should be taken by the farming industry and by the Competent Authorities to drastically reduce the use of colistin in food-producing animals, especially in turkeys, following the recommendations of the European Medicines Agency (≤5 mg/PCU). EU Member States were encouraged to set stricter national targets, ideally below 1 mg/PCU colistin. We also strongly recommend reducing the overall use of all other classes of antibiotics at primary production level, in order to mitigate the effects of the complex mechanisms behind co-selection and multidrug resistance toward Critically Important Antimicrobials, in a "Consumer Protection" and a "One Health" perspective.

AUTHOR CONTRIBUTIONS
PA, PL, AF, RH, VB, and AB conceived and designed the experiments. PA, PL, FF, AI, and FS performed the experiments. PA, AF, FF, AI, RH, VB, and AB analyzed the data. PA, PL, AF, AC, RH, VB, and AB contributed reagents, materials, and analysis tools. PA, AF, RH, VB, and AB wrote the paper.

FUNDING
This work was supported by the project 'Establishing Next Generation sequencing Ability for Genomic analysis in Europe' (ENGAGE) co-funded by the European Food Safety Authority (EFSA, GP/EFSA/AFSCO/2015/01/CT1).