First Report and Characterization of the mcr-1 Positive Multidrug-Resistant Escherichia coli Strain Isolated from Pigs in Croatia

The emergence and rapid spread of the plasmid-mediated colistin-resistant mcr-1 gene introduced a serious threat to public health. In 2021, a multi-drug resistant, mcr-1 positive Escherichia coli EC1945 strain, was isolated from pig caecal content in Croatia. Antimicrobial susceptibility testing and whole genome sequencing were performed. Bioinformatics tools were used to determine the presence of resistance genes, plasmid Inc groups, serotype, sequence type, virulence factors, and plasmid reconstruction. The isolated strain showed phenotypic and genotypic resistance to nine antimicrobial classes. It was resistant to colistin, gentamicin, ampicillin, cefepime, cefotaxime, ceftazidime, sulfamethoxazole, chloramphenicol, nalidixic acid, and ciprofloxacin. Antimicrobial resistance genes included mcr-1, blaTEM-1B, blaCTX-M-1, aac(3)-IId, aph(3’)-Ia, aadA5, sul2, catA1, gyrA (S83L, D87N), and parC (A56T, S80I). The mcr-1 gene was located within the conjugative IncX4 plasmid. IncI1, IncFIB, and IncFII plasmids were also detected. The isolate also harbored 14 virulence genes and was classified as ST744 and O101:H10. ST744 is a member of the ST10 group which includes commensal, extraintestinal pathogenic E. coli isolates that play a crucial role as a reservoir of genes. Further efforts are needed to identify mcr-1-carrying E. coli isolates in Croatia, especially in food-producing animals to identify such gene reservoirs.


Introduction
The global rise of antimicrobial resistance (AMR) and the emergence and spread of multidrug-resistant (MDR) bacteria represent a huge public health threat.Against such Gram-negative pathogens, polymyxins have been regarded as last-resort antibiotics [1].Polymyxins are a group of cationic polypeptide antibiotics that consist of five different compounds (polymyxin A-E).Currently, only polymyxin E (colistin) and polymyxin B have been put into clinical use.This group was discovered in 1947 as a metabolite produced by the bacterium Paenibacillus polymyxa (formerly known as Bacillus polymyxa var.colistinus) and was originally named 'aerosporin' [2,3].Due to its neurotoxicity and nephrotoxicity, colistin was rarely used in human medicine and was initially restricted to ophthalmic and topical use [4][5][6].However, it was widely used in veterinary medicine to treat infections caused by Enterobacteriaceae.The main indication for colistin use in veterinary medicine was the treatment of gastrointestinal tract infections caused by non-invasive Escherichia (E.) coli in pigs, poultry, cattle, sheep, goats, and rabbits.Typically, colistin products were administered orally, as a drench, in feed, in drinking water, or through milk replacer diets [7].Colistin was thought to be a 'safe' antimicrobial because the prevalence of resistance remained very low.Colistin resistance was traditionally linked to the chromosomal-encoded mechanisms that involved two-component systems, phoPQ and pmrAB, and mutation in the mgrB regulator, resulting in increasing the positive charge on lipopolysaccharides (LPSs) and decreasing colistin binding isolates, leading to colistin resistance [8][9][10][11].
The first plasmid-mediated mobilize-colistin-resistant (mcr) gene, the mcr-1 gene, was detected in 2015 in China [12].The mcr-1 gene confers colistin resistance due to the ability to encode a phosphoethanolamine transferase that modifies lipid A, reducing its attraction to colistin [13].This discovery introduced a serious threat to public health associated with the fact of horizontal mcr gene transfer and the emergence of bacteria resistant to all available antimicrobials [10,14,15].As a result, polymyxins are classified as Highest Priority Critically Important Antimicrobials (HPCIA) for human medicine by the World Health Organization (WHO) [16], and as Veterinary Highly Important Antimicrobial agents (VHIA) by the World Organization for Animal Health (WOAH) [17].According to the Antimicrobial Advice Ad Hoc Expert Group (AMEG) categorization, polymyxins are classified as Category B with the indication 'Restrict' use in veterinary medicine within the EU [18].
Multiple studies have found that the mcr-1 gene spreads rapidly in animals, humans, and the environment [15,19].Until now, nine additional mcr genes encoding colistin resistance have been identified (mcr-2 to mcr-10) [20][21][22][23][24][25][26][27][28].By mid-2022, mcr genes in E. coli have been found on five continents, with a total of 54 countries reporting publications from 2018.The mcr-1 variant is the most common and widely distributed across all continents and hosts.The total crude prevalence of worldwide distribution of mcr in E. coli in nonclinical isolates was 6.52% and 1.76% in clinical isolates [29].The mcr-1 gene has been detected on plasmids of various incompatibility (Inc) groups, with IncX4, IncI2, and IncHI2 being the most common types [30,31].There is evidence that this plasmid-borne gene is spreading rapidly among the pig population.[32].
In this study, we report the first detection of mcr-1 plasmid gene in the MDR E. coli strain isolated from pig caecal content in Croatia.The aim of this research includes describing the resistance profile against several classes of antimicrobials, as well as the virulence potential of E. coli EC1945 isolate.
Bacterial identification was performed using the VITEK 2 system (bioMérieux, Marcyl'Etoile, France).The isolate was stored at −80 • C in the Tryptic Soy broth medium containing 30% glycerol for further study.
Reference strain E. coli ATCC 25922 was used for quality control.

Polymerase Chain Reaction (PCR)
Multiplex PCR for the detection of mcr-1-mcr-5 variants was performed using the protocol recommended by EURL-AR [35].

Whole Genome Sequencing (WGS)
The E. coli EC1945 isolate genome was sequenced at EURL-AR, DTU, Lyngby, Denmark, under the scope of the EFSA confirmatory testing [36].Genomic DNA was extracted using an Invitrogen Easy-DNA KitTM (Invitrogen, Carlsbad, CA, United States), and the DNA concentrations were 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, USA) Nextera XT ® 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) (https://www.genomicepidemiology.org/) (accessed on 25 April 2022).
The Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession JAUTEE000000000.
The core genome MLST (cgMLST) method was used to compare alleles and complete linkage cluster analysis using BioNumerics software (version 8.1, Applied Maths, Belgium).The cgMLST results of our strain were compared with the results of strains present in the PubMLST E. coli database (https://pubmlst.org/organisms/escherichia-spp)(accessed on 4 May 2023) [40] based on categorical differences in alleles and complete linkage clustering analysis using BioNumerics (version 8.1, Applied Maths, Belgium).

Antimicrobial Susceptibility Testing and Polymerase Chain Reaction (PCR)
The minimum inhibitory concentrations for EC1945 E. coli isolate are listed in Table 1.This strain exhibited MDR phenotype and was found to be resistant to colistin, gentamicin, ampicillin, cefepime, cefotaxime, ceftazidime, sulfamethoxazole, chloramphenicol, nalidixic acid, and ciprofloxacin.However, it remained susceptible to azithromycin, amikacin, tigecycline, cefoxitin, temocillin, ertapenem, meropenem, imipenem, tetracycline, and trimethoprim.Considering the EUCAST guidelines [41], this isolate was identified as an Extended Spectrum Beta-Lactamase (ESBL) producer due to its resistance to cefotaxime and ceftazidime, susceptibility to cefepime, and synergism between cefotaxime/clavulanic acid and ceftazidime/clavulanic acid.Multiplex PCR analysis confirmed the presence of the mcr-1 gene.

Whole Genome Sequencing and Bioinformatic Analysis
Based on the WGS analysis, it was determined that the E. coli EC1945 isolate had a length of 5,095,464 base pairs.ResFinder 4.1 results revealed the presence of the mcr-1.1 gene, which showed a 100% identity to E. coli strain SHP45 plasmid pHNSHP45 [12].
p.A56T), and parC (p.S80I), were located on plasmids.However, quinolone and fluoroquinolone resistance were chromosomally mediated.Based on our findings, it can be concluded that the mcr-1 gene was present on the IncX4 plasmid, while the bla CTX-M-1 and aadA5 genes were located on the IncI1 plasmid, as they were identified on the same contigs.Among the virulence genes, only traT, sitA, and cea were identified on contigs that belonged to plasmids.
Using the Proksee v1.0.0a6 tool, we were able to reconstruct and annotate the IncX4 plasmid from all the plasmid contigs, as presented in Figure 1.Unfortunately, we were unable to confidently assemble the other plasmids due to significant overlap on multiple contigs.
quinolone resistance were chromosomally mediated.Based on our findings, it can be concluded that the mcr-1 gene was present on the IncX4 plasmid, while the blaCTX-M-1 and aadA5 genes were located on the IncI1 plasmid, as they were identified on the same contigs.Among the virulence genes, only traT, sitA, and cea were identified on contigs that belonged to plasmids.
Using the Proksee v1.0.0a6 tool, we were able to reconstruct and annotate the IncX4 plasmid from all the plasmid contigs, as presented in Figure 1.Unfortunately, we were unable to confidently assemble the other plasmids due to significant overlap on multiple contigs.
Through analysis of the complete plasmid sequences, it was found that the backbone of EC1945_IncX4 plasmid is highly similar to previously sequenced mcr-1-carrying IncX4 plasmids, with a similarity of over 98% (Figure 1).Within the EC1945_IncX4 plasmid, we identified the mcr-1/PAP2 casse e (2513 bp) as well as putative conjugal transfer components such as the auxiliary factor TaxA, relaxase TaxC, type IV secretion system genes (T4SS), pilX1-pilX11, and type IV coupling protein gene, taxB (T4CP).Additionally, the ori-T-like region contained a pair of 14-bp insert repeats (GCAGGTGAGCAAAG…CTTTGTTCACCTGA) (coordinates 27,161-27,194 bp).Through analysis of the complete plasmid sequences, it was found that the backbone of EC1945_IncX4 plasmid is highly similar to previously sequenced mcr-1-carrying IncX4 plasmids, with a similarity of over 98% (Figure 1).

Multilocus Sequence Typing and Core Genome Multilocus Sequence Typing
Further analysis using MLST and serotyping revealed that the E. coli strain EC1945 belonged to sequence type 744 (ST744) and O101:H10 serotype (Table 4).Figure 2 displays a dendrogram that illustrates the connection between the identified strain and the most closely related strains present in the database utilized in this study.This relationship includes information such as the source, country of origin, species, and year of isolation.The phylogenetic tree, constructed using allele numbers obtained from cgMLST analysis, is shown in Figure 3.The dendrogram classifies 60 patterns of strains, with numbers on the branches indicating the number of allelic differences multiplied by 100.The cgMLST analysis indicates that the E. coli strain ST744, which was isolated in Scotland, United Kingdom, is the closest strain to the EC1945 strain from this study, with only a 20 allele difference.Unfortunately, the isolate's origin is unknown.The entire cluster of ST744 isolates, originating from food, environment, and animals, is within 120 allele differences.

Discussion
The overuse and misuse of antibiotics have resulted in the resistance of bacteria to all known antimicrobial classes, contributing to the emergence of multi-drug resistant (MDR) bacteria in both human medicine and veterinary practices [42][43][44].Food-producing animals also play a significant role in the spread of antimicrobial resistance (AMR) as they act as a reservoir and intermediary for AMR between humans, animals, and the environment [45,46].Additionally, food-producing animals serve as an important source of AMR genes, which have a direct impact on humans and the environment [47].Colistin has been used in food-producing animals for decades, not only for the treatment of infections, but also as a growth promoter, so it is not surprising that widespread use of colistin has aided the spread of colistin resistance [10].
In recent years, there has been an increase in reports of plasmid-mediated colistinresistant E. coli isolates in both livestock and humans.These isolates, in addition to being resistant to colistin, were found to be MDR [19,29].More than 40 studies from 15 countries around the world have reported the presence of E. coli containing the mcr-1 gene in pigs [19].
In this study, we analyzed the first mcr-1 positive, MDR E. coli strain isolated in pig caecal content in Croatia.
In 2021, MDR producing, O101:H10 E. coli isolate, belonging to ST744, that displayed a colistin MIC of 8 mg/L, was revealed.It was also resistant to gentamicin, ampicillin, cefepime, cefotaxime, ceftazidime, sulfamethoxazole, chloramphenicol, nalidixic acid, and ciprofloxacin (Table 1).Antimicrobial susceptibility testing confirmed ESBL phenotype.Resistance gene analysis provided clear evidence of the relationship between phenotype and genotype resistance.
The mcr-1 gene was located on an IncX4-type plasmid, which is very common and part of an emerging plasmid expansion [32].A total of 14 different plasmid incompatibility groups capable of carrying the mcr-1 gene were identified.Over 90% of the worldwide identified plasmids belonged to IncX4, IncI2, and IncHI2 groups.The IncX4 group was the predominant Inc plasmid group carrying the mcr-1 gene in Europe (48.9%),while the IncI2 group dominated in Asia (52.1%) [30].The plasmid-mediated mcr-1 gene confers colistin resistance by encoding a phosphoethanolamine transferase, which catalyzes the addition of a phosphoethanolamine moiety to the lipid-A lipopolysaccharide of the bacterium.A structural change in lipidA reduces its affinity for polymyxin [48].
Conjugative plasmids play a critical role in the dissemination of the mcr-1 gene, and the IncI2 and IncX4 are the two leading plasmid types responsible for the global spread of colistin resistance [15,30].Previous reports have shown that the IncX4 plasmid is highly transmissible, showing 10 2 -10 5 -fold higher transfer frequencies than IncFII plasmids [20].It was also noted that IncX4 plasmids carrying the mcr-1 gene are strikingly similar and show very high architectural conservation [49].
We identified typical conjugal modules on the EC1945_IncX plasmid, including auxiliary factor TaxA, relaxase TaxC, type IV secretion system genes (T4SS), pilX1-pilX11, and type IV coupling protein gene, taxB (T4CP) (Figure 1).Conjugative T4SS was found in other mcr-1-harbouring E. coli strains and could be responsible for horizontal gene transfer in conjugative plasmids [50].Furthermore, a pair of 14-bp insert repeats were also found in the ori-T-like region, confirming that the EC1945_IncX4 plasmid is conjugative.
Although the rapid transmission of IncX4 plasmids containing mcr-1 is a major concern, the mechanisms that have enabled IncX4 plasmids to become successful vectors for the global spread of the mcr-1 gene are mostly unclear.A recent study identified a novel transfer activator PixR, which is specific to IncX4 and IncX7 plasmids.pixR directly activates the expression of transfer genes, increases transfer capability, and results in the successful dissemination of IncX4 plasmids harboring the mcr-1 gene [51].
Previous studies have shown that IncF plasmids have a conjugation transfer rate approximately 400 times lower than IncX4 plasmids, but a transfer rate 2.5 times higher than that of IncI plasmids [52].Our study on E. coli EC1945 strain revealed the presence of four plasmids (Table 3).It has been noted that transfer rates can be influenced by the number of plasmids in the bacterial cell, with strains possessing two or more plasmids displaying faster transfer rates [53].Therefore, we can confidently state that the plasmids detected in E. coli EC1945 are transferable.
The mcr-1 gene was located on the IncX4 plasmid, while the bla CTX-M-1 and aadA5 genes were found on the IncI1 plasmid.Only traT, sitA, and cea virulence genes were found on plasmid contigs.However, the specific location of the other AMR and virulence genes is unclear as there is significant overlap on multiple plasmid contigs.
The E. coli EC1945 isolate was found to be resistant to cefotaxime, ceftazidime, and cefepime with MICs of 4 mg/L, 2 mg/L, and 8 mg/L, respectively.Resistance genes analysis confirmed that this bla CTX-M-1 gene was found on the IncI1 plasmid.This gene belongs to CTX-M-type extended-spectrum β-lactamase, which confers resistance to extendedspectrum cephalosporins and is inhibited by clavulanic acid [54].The bla CTX-M-1 gene is the most frequently detected ESBL-type coded gene in animals in Europe, particularly in pigs [14,[55][56][57].
We also detected the bla TEM-1b plasmid gene in E. coli EC1945 isolate, which is the most prevalent gene in the world associated with resistance to narrow-spectrum β-lactamases and inactivation of penicillin and aminopenicillins [58].Previous studies from Australia, Denmark, Switzerland, the United States, and South Korea have described the bla TEM-1 gene as the most common gene in clinical E. coli isolates responsible for β-lactam resistance in pigs [59][60][61][62][63].
Regarding fluoroquinolones resistance, we detected known double point mutations in gyrA (S83L, D87N) and parC (A56T, S80I).These have been closely associated with resistance to fluoroquinolones.This explains the MICs of 8 mg/L and >64 mg/L to ciprofloxacin and nalidixic acid, respectively.Previously, it was pointed out that fluoroquinolone resistance in the Enterobacteriaceae family has been mainly caused by point mutations in genes encoding DNA gyrase and topoisomerase IV (gyrA, gyrB, and parC) in quinolone resistance-determining regions (QRDRs) [64][65][66].Mutations within gyrA and parC genes have been described in E. coli isolates originating from pigs and other food-producing animals [67][68][69].
In the E. coli EC1945 isolate, sulfamethoxazole resistance was mediated by the plasmid sul2 gene.Resistance to sulphonamides, which are longstanding antimicrobials and the most commonly used class of antimicrobials in animal production, is prevalent globally [70].Sulphonamide resistance in bacteria is mediated by sul1, sul2, and sul3 genes, among which the sul2 gene is most widely distributed in porcine, avian, or human E. coli sulphonamideresistant isolates [61,63,71,72].
The E. coli EC1945 strain harbored the plasmid catA1 gene, which was responsible for chloramphenicol resistance.Chloramphenicol is a broad-spectrum antibiotic that was extensively used in veterinary medicine in all major food-producing animals until concerns over its toxicity emerged [73].Since 1994, it has been banned for use in food-producing animals in the European Union (EU) [74] and in many other countries, including the USA, Canada, Australia, Japan, and China.Phenicol resistance in E. coli of animal origin is mostly mediated by enzymatic inactivation of phenicols by chloramphenicol acetyltransferases encoded by cat genes [75].The catA1 gene was isolated worldwide in E. coli isolates from pigs [61,[76][77][78].
Aminoglycosides are often used for treating complicated infections such as sepsis, pneumonia, meningitis, and urinary tract/abdominal infections, and are extremely important in both human and veterinary medicine [79].Resistance to aminoglycosides in Gram-negative bacteria is mainly due to the production of aminoglycoside-modifying enzymes or modification of the ribosome by acquired 16S rRNA methyltransferases.Aminoglycoside-modifying enzymes catalyze the modification at −OH or −NH 2 groups of the 2-deoxystreptamine nucleus or the sugar moieties and can be acetyltransferases (AACs), nucleotidyltranferases (ANTs), or phosphotransferases (APHs) [80].Mutations in the genes responsible for the synthesis of these enzymes lead to the emergence of new enzyme variants and, consequently, resistance to a greater number of antibiotics of this class.Their ability to transfer at the molecular level as part of integrons, gene cassettes, transposons, or integrative conjugative elements results in the ability of this resistance mechanism to reach virtually all bacterial types [81].It was reported that E. coli from pigs may be an important reservoir for the transfer of gentamicin resistance genes or bacteria to humans [82].In our study, we identified two plasmid genes responsible for aminoglycoside resistance: aac(3)-IId and aph (3 )-Ia.AAC-3 enzymes catalyze the acetylation of the -NH2 group in an aminoglycoside antibiotic at the third position [80] and are responsible for resistance to apramycin, gentamicin, netilmicin, tobramycin, sisomicin, and dibekacin [83].APHs catalyze the transfer of a phosphate group to the aminoglycoside molecule [84].The APH(3 )-I subclass shows a resistance profile including kanamycin, neomycin, paromomycin, ribostamycin, and lividomycin, and is composed of three enzymes that are widely distributed mainly among gram-negatives [85].
The Plasmid sitABCD gene was also found in E. coli EC1945 isolate.The SitABCD system mediates the transport of iron and manganese.Its ability to obtain manganese contributes to the resistance to oxidative stress and protection against agents such as hydrogen peroxide [86].
The E. coli strain EC1945 screened in this study belonged to ST744.It is a single-locus variant of ST10, belonging to clonal complex (CC) 10, and A phylogenetic E. coli group.This clonal complex is widely disseminated [87] and ranks as the third most common extraintestinal pathogenic E. coli (ExPEC) in a systematic review of human studies [88].ExPEC strains are associated with a variety of infections, including urinary tract infections (UTI), neonatal meningitis, septicemia, diverse intraabdominal infection, pneumonia, osteomyelitis, and soft-tissue infection [89].It has been previously emphasized that the existence and high prevalence of MDR E. coli isolates harboring the mcr-1 gene in the ST10-related population is due to their intrinsic ability to acquire AMR genes.This population of commensal E. coli isolates then plays a crucial role as a reservoir for these genes [30].MDR E. coli ST744 isolates have been found in multiple animal species, environments, and diseased humans worldwide.ST744, mcr-1 positive, isolates have been obtained from human bloodstream infection in Denmark [90], human fecal samples in China and Brazil [91,92], UTI in Portugal and Brazil [93,94], and sputum/body fluid in China [95].E. coli ST744 strains have also been reported in mcr-1-positive isolates from poultry in Romania and Lebanon [96,97], weaning pigs in Japan [98], while mcr-3-producing E. coli ST744 isolate was reported in veal calves in France [99].Additionally, previous studies have emphasized that E. coli ST744 isolates consistently had mutations in QRDRs of gyrA and parC genes [100], as demonstrated in the isolate in our study.
According to the cgMLST analysis, the E. coli strain ST744 found in Scotland is very similar to the EC1945 strain from our study, with only a difference of 20 alleles.These 2 strains form a separate cluster and are 40 alleles different from the main linkage cluster that consists of 2 strains from Ecuador and Spain, and one strain from Kenya (Figure 3).The entire cluster of ST744 isolates is within 120 allele differences, indicating quite a similarity of ST744 strains.
We analyzed the virulence factor-encoding genes (VFGs) present in the sequenced E. coli EC1945 isolate to determine its potential to cause disease.The virulence factors (VFs) present in bacteria, such as adhesins, toxins, siderophores, capsules, hemolysins, and invasins, largely determine the virulent potential of bacteria.These VFs help the microorganism to avoid host defenses, invade, colonize, and cause disease [101].According to the virulence profile (Table 2), E. coli EC1945 had extra-intestinal pathogenic potential, primarily due to the presence of fimH, gad, hlyE, hlyF, iss, and TraT genes.Previously, it was described that FimH adhesion located at the tip of the bacterial type 1 fimbrium mediates the binding to urothelial cells and prevents bacterial washout by micturition [102,103].The gad gene encodes the enzyme glutamate decarboxylase and is responsible for E. coli survival of low pH exposure and passage through the stomach after ingestion [104].Hemolysin E is one of the numerous virulence factors of E. coli strains responsible for ExPEC infections.It is a pore-forming toxin that lyses mammalian erythrocytes.It is also toxic toward cultured mammalian cells and induces apoptosis in macrophages [105,106].Hemolysin F was supposed to be an avian hemolysin associated with avian pathogenic E. coli [107], and it has been shown that plays an important role in the virulence of ExPEC [108].Another very important ExPEC virulence gene is iss, which is responsible for surviving the bactericidal effects of the complement system [109,110].It is associated with the avian pathogenic E. coli (APEC) subpathotype of ExPEC and is found to occur in around 60% of uropathogenic and neonatal meningitis-associated E. coli strains [111,112].The TraT is another gene that was described as significant VF in neonatal meningitis E. coli (NMEC), sepsis-associated E. coli (SEPEC), uropathogenic E. coli (UPEC), and APEC in the ExPEC group.It expresses a transfer protein that inhibits the classical pathway of complement activation [113].

Conclusions
We characterized the first mcr-1 positive E. coli strain isolated in Croatia.The present study provides a detailed analysis of the MDR E. coli EC1945 strain harboring the mcr-1 gene isolated from pig caecal content.The mcr-1 gene was located on the IncX4 plasmid, which is the predominant Inc plasmid group carrying the mcr-1 gene in Europe.Apart from IncX4, the E. coli EC1945 strain also harbored three other transferable Inc plasmid groups, namely IncI1, IncFIB, and IncFII.However, due to overlapping plasmid contigs, we could not reconstruct them in this study.Our data indicate that the E. coli EC1945 isolate exhibits phenotypic and genotypic resistance to nine antimicrobial classes, and possesses 14 virulence genes, which classify it as MDR and ExPEC isolate.Regarding the Achtman scheme, our strain was classified as ST744, which is a single locus variant of ST10, commensal E. coli isolates then play a crucial role as a gene reservoir.Further efforts are needed to identify mcr-1-carrying E. coli isolates in Croatia, especially in food-producing animals to identify such gene reservoirs, as we think that this strain is not an isolated case.

Microorganisms 2023 , 20 Figure 2 .
Figure 2. Dendrogram showing relationship of identified strain to the most related strains present in database used in this study according to MLST 7 Achtman scheme.Strains are also identified according to ST, source, country of origin species, and year of isolation.STs are designated by different colors.Isolate from this study is in rectangle.

Figure 2 .
Figure 2. Dendrogram showing relationship of identified strain to the most related strains present in database used in this study according to MLST 7 Achtman scheme.Strains are also identified according to ST, source, country of origin species, and year of isolation.STs are designated by different colors.Isolate from this study is in rectangle.

Figure 3 .
Figure 3.Dendrogram showing relationship of identified strain to the most related strains present in database used in this study calculated from cgMLST results.cgMLST profiles were generated using the 2513 core genes.Strains are also identified according to ST, source and country of origin species, and year of isolation.Complete linkage cluster analysis was used in the calculation with scaling factor of 100.Numbers indicated on the branches correspond to the number of allelic differences between strains multiplied by 100.Different colors indicate similarities between strains within same linkage cluster.Isolate from this study is in rectangle.

Figure 3 .
Figure 3.Dendrogram showing relationship of identified strain to the most related strains present in database used in this study calculated from cgMLST results.cgMLST profiles were generated using the 2513 core genes.Strains are also identified according to ST, source and country of origin species, and year of isolation.Complete linkage cluster analysis was used in the calculation with scaling factor of 100.Numbers indicated on the branches correspond to the number of allelic differences between strains multiplied by 100.Different colors indicate similarities between strains within same linkage cluster.Isolate from this study is in rectangle.