Emergence and Clonal Spread of CTX-M-65-Producing Escherichia coli From Retail Meat in Portugal

The emergence and dissemination of resistance to third- and fourth-generation cephalosporins among Enterobacteriaceae from different sources impose a global public health threat. Here, we characterized by whole-genome sequencing four Escherichia coli strains harboring the blaCTX–M–65 gene identified among 49 isolates from beef and pork collected at retail. The genomic content was determined using the Center for Genomic Epidemiology web tools. Additionally, the prediction and reconstruction of plasmids were conducted, the genetic platform of the blaCTX–M–65 genes was investigated, and phylogenetic analysis was carried out using 17 other genomes with the same sequence type and harboring the blaCTX–M–65 gene. All strains harbored blaCTX–M–65, blaOXA–1, and blaTEM–1B, and one also carried the blaSHV–12 gene. Other resistance genes, namely, qnrS2, aac(6′)-Ib-c, dfrA14, sul2, tetA, and mphA, were present in all the genomes; the mcr-1.1 gene was identified in the colistin-resistant strains. They belong to sequence type 2179, phylogenetic group B1, and serotype O9:H9 and carried plasmids IncI, IncFIC(FII), and IncFIB. All strains share an identical genetic environment with IS903 and ISEcp1 flanking the blaCTX–M–65 gene. It seems likely that the blaCTX–M–65 gene is located in the chromosome in all isolates based on deep in silico analysis. Our findings showed that the strains are clonally related and belong to two sub-lineages. This study reports the emergence of CTX-M-65-producing E. coli in Portugal in food products of animal origin. The chromosomal location of the blaCTX–M–65 gene may ensure a stable spread of resistance in the absence of selective pressure.


INTRODUCTION
Bacteria harboring antimicrobial resistance genes can be spread to humans, constituting a global public health concern (O' Neill, 2016). Of particular importance is the resistance to β-lactam antibiotics, enzymatic inactivation being the most common mechanism of resistance through which β-lactamases cause the cleavage of the β-lactam ring (Blair et al., 2015). Extended-spectrum β-lactamase (ESBL) is a group of enzymes that hydrolyze oxymino-beta-lactam antibiotics, conferring resistance to a wide variety of β-lactams, including penicillins, first-, second-, third-, and fourth-generation cephalosporins, and monobactams (Bush, 2018). These enzymes can be inhibited by β-lactam inhibitors such as clavulanic acid, sulbactam, and tazobactam through the covalent link to the serine residue active site (Tooke et al., 2019). Mobile genetic elements, particularly plasmids, are involved in the spread of ESBL genes, resulting in the rapid increase of ESBL-producing bacteria among different sources (Rozwandowicz et al., 2018;Silva et al., 2019). ESBLs are classified into several enzymatic groups, of which CTX-M, OXA, SHV, and TEM are frequently observed, CTX-M being the most prevalent (Silva et al., 2019;Palmeira and Ferreira, 2020). Currently, 230 CTX-M variants have been described according to GenBank records (last accessed on January 4, 2021). In Europe, the most frequently reported variants of Enterobacteriaceae species isolated from foodproducing animals and food products are CTX-M-1, CTX-M-14, CTX-M-15, and CTX-M-2, CTX-M-15 being associated with outbreaks of severe extraintestinal infections in humans caused by the multidrug-resistant (MDR) Escherichia coli ST131 (Silva et al., 2019;Palmeira and Ferreira, 2020). In Portugal, the CTX-M-1, CTX-M-14, and CTX-M-32 variants were reported as the most prevalent in commensal and pathogenic E. coli isolated from food-producing animals and food products (Clemente et al., 2019;Silva et al., 2019).
In 2008, a new variant of the CTX-M family identified as CTX-M-65, belonging to the CTX-M-9 group and cluster 14, was described for the first time in E. coli isolated from a human urine sample in the United States (Doi et al., 2008). CTX-M-65 differs from CTX-M-14 by two amino acid substitutions, namely, alanine by valine at position 77 (A77V) and serine by arginine at position 274 (S274R). In 2017, Tate and colleagues reported the emergence of a CTX-M-65 Salmonella Infantis isolated from food animals, retail chickens, and humans in the United States that is highly resistant to most of the studied β-lactam antimicrobials (Tate et al., 2017). In Europe, the CTX-M-65 variant has also been reported from Italy in Salmonella Infantis isolated from broilers and humans (Franco et al., 2015) and from the Netherlands in E. coli from cattle (Palmeira and Ferreira, 2020). Globally, the enzyme is widely distributed, with reports from China, Korea, and South America in E. coli isolated from humans, foodproducing animals, retail chickens, and in a giant anteater from a zoo (Zheng et al., 2012;Bartoloni et al., 2013;Rao et al., 2014;Tate et al., 2017;Furlan et al., 2019;Park et al., 2019;Vinueza-Burgos et al., 2019;Wang et al., 2020). It is noteworthy to mention that a multidrug carbapenemase strain of Klebsiella pneumoniae co-producing CTX-M-65 was the causative agent of a severe nosocomial outbreak in China (Zhan et al., 2017).
In the present study, four CTX-M-65-producing E. coli isolated from beef and pork samples collected at retail in 2017 were characterized by whole-genome sequencing (WGS). Here, we described the molecular epidemiology of the resistance genes, the identification of the plasmids, and the bla CTX−M−65 genetic environment. Moreover, we determined and described the genetic relatedness with other E. coli genomes for better insights into the public health impact of an ESBL producer rarely found in Europe.

Bacterial Isolates and Antimicrobial Susceptibility Testing
Two hundred and twenty beef and 220 pork samples were collected at retail stores across mainland Portugal in compliance with the European Commission Implementing Decision of November 12, 2013 to monitor and report antimicrobial resistance in zoonotic and commensal bacteria (Commission Decision 652/2013) in 2017. The isolation and identification of extended-spectrum β-lactamase/plasmid-mediated AmpC (ESBL/PMAβ) E. coli producers from meat samples were performed according to the laboratory protocols defined by the European Union Reference Laboratory for antimicrobial resistance (EURL-AR) 1 . Briefly, 25 g of each meat sample was mixed with 225 ml of buffered peptone water, followed by incubation at 37 • C for 18-22 h. Enriched samples were plated onto MacConkey agar supplemented with 1 mg/L of cefotaxime (Glentham, Corsham, United Kingdom) and incubated at 44 • C for 18-22 h. Presumptive E. coli colonies were selected for biochemical identification on ChromID R coli agar (bioMérieux, Marcy-l'Étoile, France), and after confirmation, the isolates were sub-cultured and stored at −80 • C before further analyses.
The isolates were tested for antimicrobial susceptibility through the determination of the minimum inhibitory concentrations (MICs) using commercially available 96well microplates, EUVSEC and EUVSEC2 (Sensititre R , Trek Diagnostic Systems, East Grinstead, United Kingdom) panels, and the results were interpreted according to EUCAST epidemiological breakpoints 2 .

Molecular Characterization of β-Lactam Resistance
Resistance mechanisms associated with ESBL/PMAβ enzymes were screened by PCR using primers targeting bla TEM , bla SHV , bla OXA , bla CTX−M , bla ACC , bla FOX , bla MOX , bla DHA , bla CIT , and bla EBC (Dallenne et al., 2010). Amplified products were purified with ExoSAP-IT TM (Applied Biosystems TM , Warrington, United Kingdom), followed by Sanger sequencing using the BigDye R Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems). The sequencing of fragments was performed in an automatic sequencer ABI3100 (Applied Biosystems), and the identification of resistance genes was determined using the Basic Local Alignment Search Tool (BLAST) from the NCBI website (Altschul et al., 1990) and The Comprehensive Antibiotic Resistance Database (CARD) (Alcock et al., 2019).

Whole-Genome Sequencing and Bioinformatics Analysis of the CTX-M-65 E. coli Producers
Four isolates identified as positive for bla CTX−M−65 by PCR and Sanger sequencing were further characterized by whole-genome sequencing. The isolates were recovered from pork and beef samples, collected on the same day from the same retail store in the north of Portugal (INIAV_ECX027 and INIAV_ECX036) and on different days from different retail stores in Lisbon and Tejo Valley (INIAV_ECX016 and INIAV_ECX035).
The genome of two isolates (INIAV_016ECX and INIAV_027ECX) was sequenced at EURL-AR, DTU, Lyngby, Denmark, under the scope of the European Food Safety Authority (EFSA) confirmatory testing. Genomic DNA was extracted using an Invitrogen Easy-DNA Kit TM (Invitrogen, Carlsbad, CA, United States) and the DNA concentrations determined using the Qubit dsDNA BR assay kit (Invitrogen). Genomic DNA was prepared for Illumina pair-end sequencing using the Illumina (Illumina, Inc., San Diego, CA, United States) Nextera XT R Guide following the protocol revision C1. A sample of the pooled Nextera XT Libraries was loaded onto an Illumina MiSeq reagent cartridge using MiSeq Reagent Kit v3. The libraries were sequenced using an Illumina MiSeq platform (Illumina). The raw reads were de novo assembled using the assembler pipeline (version 1.4) available from the Center for Genomic Epidemiology (CGE) 3 . Raw sequence data from these two isolates were submitted to the European Nucleotide Archive (ENA) 4 under study accession numbers: ERS3535656 and ERS3535669.
DNA extraction of the remaining two isolates (INIAV_035ECX and INIAV_036ECX) was carried out at INIAV, Oeiras, Portugal, using the PureLink R Genomic DNA kit (Invitrogen) according to the manufacturer's instructions, with minor modifications. Briefly, the incubation period at 55 • C was performed for 90 min and the DNA eluted with 50 µl of Tris-HCl buffer, pH 8.5. The DNA quality and quantity were assessed using a spectrophotometer [NanoDrop R 2000, Thermo Scientific, emergency use authorization (EUA), Waltham, MA, United States] and sequenced using the Illumina HiSeq sequencing technology (NovaSeq 6000 S2 PE150 XP sequencing mode, Eurofins Genomics Europe Sequencing GmbH, Ebersberg, Germany). The raw sequence data from these two isolates were submitted to the ENA under study accession numbers: ERS5493675 and ERS5493676. Raw data quality was assessed by FastQC 5 . BBDuk from the BBTools package 6 was used to remove possible contamination by adapter sequences and for trimming/removing low-quality reads, all performed with a minimum quality of Q20 using the Phred algorithm, with a minimum read length of 50 and with a k-mer length parameter of 19. All pre-processed reads were assembled with SPAdes 3.12.0 (Nurk et al., 2013).
Additionally, PLACNETw (Plasmid Constellation Network) was used to predict and reconstruct the plasmids (Vielva et al., 2017) in specifically assembled contigs of the bacterial genomes. For identifying the genetic platform of the CTX-M-65 enzyme, the contigs containing the bla CTX−M−65 gene were annotated using Prokka version 1.14.6 (Seemann, 2014), followed by analysis with Artemis (Carver et al., 2012) and EasyFig version 2.2.5 (Sullivan et al., 2011). A plasmid database, PLSDB, was also used to search for the plasmid nucleotide sequences contained in each of the selected contigs using the "mash screen" option (Galata et al., 2018).
A phylogenetic analysis based on the single nucleotide polymorphisms (SNPs) present in the genomes using CSI Phylogeny version 1.4 (10 reads of minimal depth at SNP positions, 10% minimal relative depth at SNP positions, 10 bp of minimal distance between SNPs, minimal SNP quality of 30, minimal read mapping quality of 25, and a minimal Z-score of 1.96)  from the CGE website was conducted with the E. coli isolates from this study and 17 E. coli genomes with the same sequence type (ST2179) and harboring the bla CTX−M−65 gene, retrieved from EnteroBase (Zhou et al., 2020). The tree was visualized using FigTree version v1.4.3 7 .

Characterization of Antimicrobial Resistance
From 49 isolates of E. coli phenotypically resistant to thirdgeneration cephalosporins, a total of 10 different bla genes were identified in 42 ESBL and seven PMAβ producers: The bla CTX−M−65 gene was detected in four isolates (8.2%) recovered from beef (3/26) and pork (1/23) retail meat samples, and all exhibited a MDR phenotype, being resistant to ciprofloxacin, nalidixic acid, azithromycin, chloramphenicol, tetracycline, trimethoprim, and sulfamethoxazole. Resistance to colistin was observed in all, except isolate INIAV_ECX035 (Table 1).

Genome Analysis of the CTX-M-65 E. coli Producers
The assemblies of the reads originated between 70 and 120 contigs and a genome size of about 5 Gb. The genotypic traits of the four isolates, their resistome and mobilome, are summarized in Table 2. According to the ResFinder tool, three β-lactam-encoding genes were found in all isolates, namely, bla CTX−M−65 , bla OXA−1 , and bla TEM−1B . Moreover, one isolate (INIAV_ECX035) also carried bla SHV−12 . The mcr-1.1 gene was identified in the three isolates resistant to colistin, and plasmidmediated quinolone resistance (PMQR) genes, namely, qnrS2 and aac(6 )-Ib-cr, were detected in all isolates. Additional genes conferring resistance to trimethoprim/sulfamethoxazole (dfrA14 and sul2), tetracycline (tetA), and azithromycin (mphA) were also present. Point mutations were identified at the gyrA and parC subunits of DNA, namely, serine by lysine at position 83 (S83L) and serine by isoleucine at position 80 (S80I), respectively, conferring resistance to quinolones. Based on these results, the phenotype and genotype were in accordance.
According to the PlasmidFinder tool, IncFIC (FII) and IncFIB replicons were identified in all isolates, and p0111 was also predicted in one isolate (INIAV_ECX016). Isolate INIAV_ECX035 also harbored the IncI1-I plasmid, while the remaining isolates carried the IncI2 plasmid. All strains have identical profiles regarding MLST (ST2179), fumC65/fimH32 alleles, serotype (O9:H9), and phylogroup (B1). The somatic antigen of the INIAV_ECX036 isolate was not typable. VirulenceFinder predicted the presence of five virulence factors (lpfA, iroN, iss, cma, and gad) in all isolates except one (INIAV_ECX016), which carries four as gad is absent.
Using PLACNETw, a new assembly was generated from the raw reads, and a network was produced representing the contigs defined as belonging to the chromosome and plasmids, according to the reference genomes from the software database. The contigs were identified as belonging to plasmids based on the recognition of relaxase and/or replicon protein sequences. A manual pruning of the original network ( Figure 1A) was performed according to Lanza et al. (2014) to reconstruct the graphical representation of the genome (Figure 1B). With the FASTA files obtained from this assembly, a new analysis with PlasmidFinder and ResFinder was performed to identify in which contigs were the antimicrobial resistance genes and plasmids located. Hereafter, it was possible to pinpoint the exact position of the resistance genes in the genome's graphical representation. Thus, based on the analysis of all isolates, the bla CTX−M−65 gene, along with floR and mdf A, was identified in contigs belonging to the chromosome, and the remaining resistance genes, including the other bla genes, were located in contigs belonging to plasmids ( Figure 1B). In isolates INIAV_ECX016, INIAV_ECX027, and INIAV_ECX036, genes bla CTX−M−65 and floR were located in the same contig, and mdf A was in a different contig, while in isolate INIAV035, the three genes were in different contigs. Most resistance genes were located in IncFIC (FII) and IncFIB plasmids, except bla SHV−12 and mcr1.1, which were carried in IncI-1 and IncI2, respectively.
The analysis of the genetic platform of the bla CTX−M−65 genes using the EasyFig tool revealed that all isolates have the IS903 (IS5 family) and ISEcp1 (IS1380 family) flanking the bla CTX−M−65 gene, and all but one isolate also harbored the transposon TnAs3 (Tn3 family), IS1006 (IS6 family) and ISVsa3 (IS91 family). In this isolate (INIAV_ECX035), the contig starts  in a different position, and its genetic platform seems to be not fully represented (Figure 2A). Consequently, the genetic platform is split and the downstream elements (TnAs3, IS1006, and ISVsa3) are not represented in this contig, being in a different one. Looking at a larger region of the contigs (Figure 2B), it is possible to realize that all four contigs have high homology between each other. The alignment of each contig using EasyFig with the plasmid sequences identified by PlasmidFinder revealed no homology between those sequences. Comparison of the sizes of the contigs containing the bla CTX−M−65 gene (ranging from 227,681 to 340,574 bp, except for the INIAV_ECX035 isolate with 82,631 bp in length) with the sizes of the plasmids identified in this study (ranging from 64,015 to 99,159 bp) revealed that the contigs are mostly wider than the plasmids. Furthermore, no plasmids were also identified, using the "mash screen" option of the PLSDB tool, in the contig's input sequence containing the bla CTX−M−65 gene from each isolate.
An unrooted phylogenetic tree (radial cladogram option) using E. coli ATCC 25922 as the reference genome was constructed using 21 genomes from different geographic regions worldwide and isolated from multiple sources (Figure 3). Multiple phylogenetic groups can be distinguished in the tree: one group includes two strains, one from Nigeria and one from Colombia isolated from poultry and human samples; another phylogenetic group including European strains from the United Kingdom, Hungary, and Italy, isolated mainly from poultry; other groups formed with one strain from the United States and one strain from China isolated from human samples; another group includes strains from China, Singapore, and Canada isolated from human and livestock samples; and two groups with only one strain each from Portugal and Japan isolated from bovine and environmental samples. All shared about 43,430 SNPs against the reference genome. The SNP analysis showed FB7969_China_Human_2018 as the strain having more SNPs, between 1,636 and 1,746. Most shared less than 112 SNPs with the four strains from this study, except for HB6973_China_Livestock_2015 and MA2351_Japan_Environment_2015 separated by between 219 and 394 SNPs. Three Portuguese strains were grouped with one strain from the UK (RA3201_UK_2018), separated by 37 and 45 SNPs. The fourth Portuguese strain (INIAV_ECX035) is separated from the remaining three from this study by 85 and 101 SNPs, being in a different phylogenetic group. The most closely related strains from this study are INIAV_ECX016 and INIAV_ECX036, sharing only 8 SNPs between each other; these two strains were obtained from beef samples collected from distinct geographic regions in Portugal.

DISCUSSION
In this study, four CTX-M-65-producing E. coli were identified among 49 ESBL/PMAβ producers isolated from meat. Here, bla CTX−M−65 is identified in Portugal for the first time in food of animal origin, with no reports linking to human infections; a complete characterization by WGS of the four MDR E. coli harboring this gene is described. E. coli harboring the bla CTX−M−65 gene is commonly found in food-producing animals and meat from Southwest Asian and South American countries (Rao et al., 2014;Na et al., 2019;Vinueza-Burgos et al., 2019;Zurita et al., 2019). Although rarely occurring in Europe, previous studies have reported CTX-M-65-producing E. coli from Dutch beef calves (Ceccarelli et al., 2019) and wild birds from Switzerland (Zurfluh et al., 2019). In our study, three isolates were from beef and one from pork; therefore, we cannot confirm the source of the CTX-M-65 E. coli detected in meat. There are several potential sources of bacteria in meat, including the animals from which the meat was derived, cross-contamination from other products, equipment and the environment, and the workers who are producing and handling the meat (EFSA and ECDC, 2019).
Based on PLACNETw, the bla CTX−M−65 genes were located in contigs with homology to the chromosome. Moreover, the contigs containing this gene are wider than the plasmids identified in this study, and none of the four contigs showed homology to any sequence of the identified plasmids, reinforcing that the bla CTX−M−65 gene is in the chromosome. The chromosomal location of bla CTX−M genes has already been reported, namely, bla CTX−M−2 (Zhao and Hu, 2013;Ferreira et al., 2014), bla CTX−M−14 and bla CTX−M−15 (Hamamoto and Hirai, 2019), bla CTX−M−55 (Zhang et al., 2019), and bla CTX−M−65 (He et al., 2017).
The successful spread of ESBL enzymes is based on their ability to disseminate their resistance genes on mobile genetic elements to other bacteria, the skill to expand their spectrum of activity, and the acquisition of point mutations (Bush, 2018;Galal et al., 2018;Tooke et al., 2019). Nevertheless, the chromosomal location of the resistance genes can benefit the stable propagation of resistance, regardless of the bacterial host's habitat (Yoon et al., 2020). The genetic platform analysis revealed that all share an identical genetic environment with IS903 and ISEcp1 flanking the bla CTX−M−65 gene downstream and upstream, respectively. ISEcp1 is one of the most important insertion sequences associated with bla CTX−M genes (Canton et al., 2012;Zhao and Hu, 2013) and Tn3, a conjugative transposon. These mobile genetic elements found in the chromosome may promote excision and reintegration in a new chromosome or transference to other bacteria through a conjugative plasmid (Canton et al., 2012).
Particularly worrying is the co-occurrence of genes encoding resistance to other critically important antimicrobials, namely, fluoroquinolones, macrolides, and polymyxins, except for isolate INIAV_ECX035, which was susceptible to colistin. This isolate also carried the bla SHV−12 gene, exhibiting a higher MIC to ceftazidime (MIC = 32 µg/ml) compared to the remaining isolates (MIC = 1-2 µg/ml), which is in accordance with the previously described (Maina et al., 2011).
Of note is that all isolates were resistant to azithromycin and harbored the mph(A) gene, confirming this gene's relevant role in macrolide susceptibility. As previously reported (Gomes et al., 2019), most of the mph(A)-carrying isolates show a MIC > 32 µg/L, as observed in our study, where the MIC values were between 32 and 64 µg/L. Moreover, other antimicrobial resistance genetic determinants were found in the isolates, including those also associated with resistance to antimicrobials frequently used in the rearing of food-producing animals, such as sulfamethoxazole (sul1 and sul2), trimethoprim (dfrA14 and dfrA17), phenicols (catB3 and floR), tetracycline (tetA), and aminoglycosides [aph(3 )-Ib and aph(6)-Id] (ESVAC, 2020).
In all isolates, different replicon-typing plasmids were identified [IncFIC(FII) and IncFIB] carrying most of the resistance genes. IncF plasmids are frequently described from human and animal sources and are considered epidemic resistance plasmids, bearing the greatest variety of resistance genes in Enterobacteriaceae (Rozwandowicz et al., 2018). Although the precise gene location on plasmids was not determined in this study, based on PLACNETw, mcr-1.1 was located on IncI2 replicon-typing plasmid. IncI2 plasmids have been associated with the mobilization of mcr genes widely spread in Europe in E. coli isolates from animals and humans (Rozwandowicz et al., 2018;Migura-Garcia et al., 2020). The IncI-1 plasmid predominantly described in Europe was identified in the INIAV_ECX035 strain carrying bla SHV−12 , in accordance with previous reports (Rozwandowicz et al., 2018).
WGS analysis also revealed that all isolates belong to ST2179, phylogroup B1, and fumC65/fimH32 type, suggesting that all have a common clonal origin, with minor differences. Recently, E. coli ST2179 bearing the bla CTX−M−65 gene but belonging to phylogroup A was reported from ducks in South Korea (Na et al., 2019). Isolates belonging to phylogroup B1 are commonly associated with non-pathogenic commensal E. coli reported from humans, animals, and food products (Bailey et al., 2010;Coura et al., 2015;Scheinberg et al., 2017;Belaynehe et al., 2018;Zurita et al., 2019).
The phylogenetic analysis revealed three isolates (INIAV_ECX016, INIAV_ECX027, and INIAV_ECX036) grouped in the same cluster, showing high genetic homology between each other. The UK strain (RA3201_UK_2018) was closely related to the Portuguese strains with 37-45 SNPs. The in silico analysis of the UK strain revealed the bla CTX−M−65 gene to also be located in the chromosome. These findings suggest the clonal spread of CTX-M-65-producing E. coli isolates in Europe. INIAV_ECX035 is separated from INIAV_ECX027 by 101 SNPs and from the remaining two strains by 85 SNPs. Although from different animal species, the isolates from the same retail store (INIAV_ECX027 and INIAV_ECX036) are closely related, pointing out the hypothesis of cross-contamination. Moreover, the phylogenetic tree and the WGS analysis suggest that the four isolates belong to two sub-lineages, one composed of the three strains that grouped and the second with the fourth strain. INIAV_ECX035 showed some differences regarding the resistance genes and plasmids compared with the other three strains: the presence of the bla SHV−12 gene conferring a higher resistance to ceftazidime and the absence of resistance to colistin, also the presence of the IncI1-I plasmid instead of the IncI2 plasmid carrying the mcr-1.1 gene. These genetic differences may justify the higher number of SNPs found and the existence of two evolutionary sub-lineages suggesting that the CTX-M-65 variant is emerging in our country.
The emergence and clonal spread of E. coli harboring bla CTX−M -65 can be a problem for the livestock industry and human health, given their multidrug resistance profile to critically important antimicrobials and their presence in the food supply chain. To our knowledge, this is the first time that CTX-M-65 is identified in Portugal in food products of animal origin. The chromosomal addition of the bla CTX−M−65 gene may ensure the spread of resistance in the absence of selective pressure. A better understanding of the factors that contribute to the emergence and dissemination of ESBL genes rarely seen in Europe, but highly prevalent in Southwest Asian and South American countries, is strongly advisable and is worthy of close monitoring. Tourism and migration flow to Portugal and the trade treaties established between countries by importing meat and meat products may be sources of cross-contamination with uncommon MDR strains, facilitating their dissemination to the community.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below. https://www.ebi.ac.uk/ ERS3535656, ERS3535669, ERS5493675, and ERS5493676.

AUTHOR CONTRIBUTIONS
CL contributed to the whole-genome sequencing (WGS) experiments, bioinformatics analysis, and interpretation of the data and wrote the manuscript. LM contributed to the laboratory experiments. AS, IH, and RH contributed to the WGS and bioinformatics analysis. LC and AA designed the study, interpreted the data, and reviewed and edited the manuscript. All authors read and approved the manuscript.

FUNDING
This work was supported by the project PTDC/CVT-CVT/28469/2017 "CIAinVET-Food-producing animals as reservoirs of resistance to Critically Important Antibiotics" financed by the "Fundação para a Ciência e Tecnologia" (FCT), Portugal. CL was also funded by National Funds through FCT-Foundation for Science and Technology under the Project UIDB/05183/2020.