Longitudinal Monitoring Reveals Persistence of Colistin-Resistant Escherichia coli on a Pig Farm Following Cessation of Colistin Use

Introduction

Colistin (polymyxin E) is one of the World Health Organization's highest priority antimicrobials: it is regarded as a last resort antibiotic, and is the treatment of choice for multidrug-resistant Enterobacteriaceae infections (1). Unfortunately, the emergence of mobile colistin resistance genes of the mcr gene family has jeopardized the efficacy of colistin. The plasmid mediated colistin resistance gene (mcr-1) was firstly identified in E. coli of porcine origin from China (2). Subsequently, other mcr variants including mcr-2 to mcr-10 were discovered mainly from members of the Enterobacteriaceae family (3–6) from different geographical areas (7, 8). The mcr genes encode phosphoethanolamine transferases enzymes which change the lipid A portion of the lipopolysaccharides (LPS), suppressing colistin binding (9). The mcr genes have been reported not only from various livestock origins (pigs, poultry, bovine) (10, 11) and food products (12) but also from the environment as well as from humans (13). Since the extensive usage of colistin in livestock farms played a major role in the occurrence of colistin resistant mcr genes, controlling the dissemination of these resistant genes from farms to the environment has become a critical concern (14). Moreover, mcr genes could be co-located with other important antibiotic resistance genes such as Extended spectrum beta-lactamase (ESBL) and carbapenemase genes (15, 16). These reports raised awareness of colistin usage and the challenge to clinical medicine.

In the swine industry, colistin had been applied therapeutically and/or prophylactically in several countries (17). Because of the importance of colistin usage in clinical infections, many countries have restricted the prophylactic usage of colistin in pig productions (18). Following the first identification of mcr-1 during nation-wide surveillance in Thailand (19), from the start of 2017 the Department of Livestock and Development (DLD) has prohibited prophylactic use of colistin sulfate in pig farms. Data regarding colistin resistance in bacteria from livestock in Thailand is still limited, although mcr positive E. coli have been detected in pig slaughterhouses from the Thailand border areas even after implementation of the colistin withdraw policy (20). Therefore, currently it is still debatable how this withdrawal of colistin usage may have influenced the emergence and spread of mcr genes in pigs and in the farm environment (21). It is thought to be unlikely that resistance plasmids can be entirely eliminated from bacterial populations (22). Previous studies have found withdrawal it was likely to be beneficial in controlling the emergence of mcr-1 in pigs once the selective pressure is removed (21, 23). However, the variations in duration of mcr gene persistence after cessation of colistin usage and the rate of dissemination from pigs to the farm environment are still concerning. Moreover, resistant bacteria from livestock can potentially spread to farmers or the environment resulting in the occurrence of antibiotic resistance genes (ARG). Furthermore, some studies have shown that even after drastic reductions of antibiotic use on farms, antibiotic resistant bacteria could be maintained in the farm by various factors (24–26).

The aim of this study was to determine the occurrence and extent of persistence of MCRPE following the cessation of colistin sulfate use on a representative pig farm which had a history of a high prevalence of MCRPE (19). Representative E. coli isolates that contained mcr genes were obtained from pigs, wastewater and farm workers over the study period and were characterized for antimicrobial susceptibility patterns, virulence factors, plasmid replicons, and clonal relationships.

Materials and Methods

Study Area and Farm Selection

A typical industrial pig farm with more than 1,000 breeder sows located in the central area of Thailand was selected for use in this study. Prior to 2017, colistin sulfate had been administered routinely to all suckling piglets from birth to weaning to prevent and control diarrhea. It was given via the water at a dose of 10 mg/kg body weight. The farm withdrew prophylactic colistin use in piglets from the beginning of 2017, following the guidance of the DLD. The farm management systems were not otherwise altered, and they continued to follow the recommendations of the Thai standard livestock farm criteria. Piglets with diarrhea were separated from healthy piglets by placing them in separate pens until they recovered. In cases of diarrhea in breeding sows and piglets, antibiotic injections including gentamicin, ceftriaxone, and/or penicillin/streptomycin combinations were used for treatment of individual animals.

Sample Collection and Processing

The number and types of samples (from pigs, wastewater and humans) that were obtained are summarized in Table 1. Sample sizes were calculated based on the prevalence of mcr genes detected in the pig farm from our previous study (19) by using Epitools program http://epitools.ausvet.com.au. Samples were collected at five-time points spanning a 3.5 year period from cessation of colistin use: June 2017, September 2018, March 2019, April 2019, and June 2020.

TABLE 1

Table 1. Details of the sample types and numbers collected at five different sampling times between 2017 and 2020, and numbers of samples found positive for MCRPE.

Approximately 25 g of fecal samples were collected from parity 1–6 farrowing sows, aged between 1 to 3 years (n = 50). Rectal swab samples from 21-day old suckling piglets belonging to the sows that were sampled were also collected (n = 20). Each farrowing sow with their respective litters were kept in farrowing pens, and at each visit one or two sows were sampled from each zone of the farrowing house. The same pens were visited at each sampling time, although the same sows were not necessarily sampled because of animal movements. In September 2018, only fecal samples from sows were collected since at the time of sampling the newly weaned piglets had been moved to another farm.

Wastewater samples (n = 50) from the wastewater tanks on the farm were collected, with 10 samples obtained at each visit. The wastewater was composed of pig manure along with the water used to clean the pig housing. Approximately 500 ml volumes were collected from wastewater tanks located before and after-biogas treatment, which were sited close to the sampled pig pens. The biogas process involves anaerobic fermentation by fermentative, acetogenic, and methanogenic bacteria to produce methane, carbon dioxide, hydrogen, and hydrogen sulfide gases. In addition, at the request of the company, at each sample collection time the farm submitted rectal swab samples from the same 10 farm workers for routine diagnostic purposes (n = 50).

Sampling from the pigs and the wastewater was conducted by an authorized veterinarian for the farm. The biohazard execution control was approved by the Institutional Biosafety Committee of the Faculty of Veterinary Science, Chulalongkorn University (IBC 2031011). The wastewater sample collection protocol was applied according to HACH water analysis guidelines (27).

Bacterial Isolation and Identification

All samples were enriched in EC broth (Difco) containing 2 μg/ml colistin sulfate at a 1:9 ratio and incubated at 37°C overnight. The sample suspensions were grown on eosin-methylene blue (EMB) (Oxoid) agar containing 2 μg/ml colistin sulfate and were incubated overnight. One to three representative colonies with a characteristic metallic sheen on the EMB plates were randomly chosen and sub-cultured on tryptic soy agar (TSA) (Difco) from the samples from which growth was obtained. The colonies were identified as E. coli using IMViC biochemical tests and Matrix-Assisted Laser Desorption Ionization combined with time of-flight analysis (MALDI Biotyper, Bruker, USA), according to the manufacturer's recommendations (28). For minimal inhibitory concentration (MIC) determinations, antibiotic susceptibility testing, and PCR detection for virulence genes and plasmid replicon types, a single representative isolate from each positive sample was used.

Antimicrobial Susceptibility Testing

The MIC for colistin was determined using the broth microdilution technique following CLSI guidelines (29). An MIC value of >2 μg/ml was considered to indicate colistin resistance (29). The antibiogram for E. coli isolates was determined using the AST-GN 38 test kit in a Vitek2 compact automated susceptibility level detection apparatus (BioMérieux, France). The antimicrobial groups that were included in Vitek2 were synchronized with veterinary guidelines (30). The 18 antimicrobials comprised amikacin (AK), amoxicillin (AMX), amoxicillin/clavulanic acid (AMC), ampicillin (AMP), cefalexin (CEX), cefpodoxime (CPD), cefovecin (INN), ceftiofur (XNL), chloramphenicol (C), enrofloxacin (ENR), gentamicin (GEN), imipenem (IMP), marbofloxacin (MBR), nitrofurantoin (NIT), piperacillin (PIP), tetracycline (TET), tobramycin (TOB), and trimethoprim/sulfamethoxazole (SXT). The MIC interpretations from the Vitek2 machine system (version-9) were made according to the Food and Drug Administration recommendations (31), CLSI guidelines (32) and EUCAST values (33). E. coli ATCC 25922, Pseudomonas aeruginosa ATCC 27853, and Staphylococcus aureus ATCC 25913 were used as the control strains.

Detection of Plasmid-Mediated Colistin Resistance Genes

Genomic DNA was extracted from all available MCRPE isolates using the Thermo Scientific GeneJET Genomic DNA Purification Kit (Thermo Fisher Scientific). Multiplex-PCR was used to detect mcr1-5 genes, following a previously published protocol (34). E. coli strain CUP13 (35) that is positive for mcr-1 and mcr-3 as confirmed by Sanger sequencing was used for the positive control, and E. coli ATCC25922 was the negative control. The PCR conditions for mcr 6, 7, and 8 were adjusted and performed according to a previous description (4).

Plasmid Replicon Typing

The 18 plasmid replicon types of Enterobacteriaceae were investigated by a set of multiplex and simplex PCRs. The primers used and the PCR conditions followed previously described methods (36). Briefly, PCR amplification, except the F-simplex, were conducted at 94°C for 5 min, followed by 30 cycles at 94°C for 1 min, 60°C for 30 s, 72°C for 1 min. The amplification was concluded with an extension program of 1 cycle at 72°C for 5 min. The PCR for F-simplex was performed in the same way except for annealing at 52°C.

Detection of Virulence Genes

The mcr positive E. coli were examined for virulence genes that are commonly present in enterotoxigenic E. coli (ETEC) and enterohaemorrhagic E. coli (EHEC) by using previously described PCRs (37). Previously sequenced ETEC and EHEC strain were used as positive controls (38). The PCR assays were performed with GoTaq^® green master mix (Promega, USA) with the thermocycler conditions being an initial denaturation at 94°C for 10 min, followed by 30 cycles of denaturation at 94°C for 30 s, and annealing at 55°C for 45 s. Extension was at 72°C for 1.5 min increased by 3 s each cycle, followed by a final extension at 72°C for 10 min.

Conjugation Assay

To determine whether mcr genes were located on transmissible plasmids, and their transferability rate, conjugation assays were performed by the broth mating technique (39). All the mcr positive strains detected by PCR were designated as donors, and E. coli J53, resistant to sodium azide, was used as the recipient strain. This recipient E. coli J53 strain is negative for fertility factors, is resistant to sodium azide (MIC >512 μg/ml) and is sensitive to colistin (MIC <2 μg/ml). Briefly, an overnight culture of bacterial colonies was diluted in Lysogeny broth (LB) and adjusted to OD600 value 1. A 1:1 ratio of donor and recipient then was mixed to obtain a final volume of 2 ml which was incubated overnight. Ten-fold serial dilutions of the overnight mixture were plated on LB agar (Oxoid) plates containing colistin (2 μg/ml) and sodium azide (100 μg/ml). The plates were incubated at 37°C for 2 days, and the transconjugant colonies were counted. The rate of conjugal transfer frequencies was calculated by dividing the number of transconjugant colonies by the number of donor colonies (40). Phenotypic colistin susceptibility testing, PCR detection of mcr genes and replicon type detections were repeated on the transconjugants.

Pulsed-Field Gel Electrophoresis (PFGE)

To investigate clonal relatedness, PFGE was performed on all 65 available MCRPE isolates from the 33 positive samples (one to three isolates per sample), following the Centers for Disease Control and Prevention standard protocol (41). Briefly, overnight cultures of E. coli isolates were suspended in cell suspension buffer, and the cells were treated with proteinase K and mixed with the agarose gel solution. The gel plugs then were treated with lysis solution, and DNA in the plugs was digested with restriction enzyme XbaI (Thermo Scientific). Gel electrophoresis was undertaken using a Bio-Rad CHEF-DRIII system, with a 200 V field at an angle of 120° run for 17–20 h, incorporating Salmonella serovar Braenderup H9812 DNA as a standard. Dendrograms were created using the GeneTool program (Syngene, India) and analyzed by the GeneDirectory program (Syngene, India).

Multi-Locus Sequence Typing (MLST)

A representative isolate from each of the 34 PFGE pulsotypes that were identified was randomly selected and included for MLST typing. A simplex PCR was performed for each of the 7 housekeeping genes of E. coli used in the Achtman MLST scheme (42). These genes encoded isocitrate/isopropyl malate dehydrogenase (icd), ATP/GTP binding motif (recA), adenylate kinase (adk), DNA gyrase (gyrB), malate dehydrogenase (mdh), adenyl succinate dehydrogenase (purA) and fumarate hydratase (fumc). The sequences were obtained using the Sanger sequencing platform. The E. coli MLST database at http://mlst.warwick.ac.uk/mlst/dbs/Ecoli was used to determine allele and sequence types (STs).

Data Analysis

The colistin-resistance rates and virulence gene profiles for the representative isolates were described as percentages compared to different sources in each sample collection. The mcr positive rates among the samples and the association between each sample collection time were analyzed using Fischer's Exact Test (p ≤ 0.05).

Results

Detection of Colistin Resistant E. coli

A total of 33 of the 170 samples (20.6%) yielded colistin-resistant E. coli, and their MICs to colistin varied from 4 to 8 μg/ml. These positive samples were isolated from pigs (n = 20/70, 28.6%), wastewater (n = 9/50, 18%) and humans (n = 4/50, 8%). A comparison of the prevalence of MCRPE isolates for each sample type over the 3.5 years since colistin cessation is shown in Figure 1, and detailed information about the isolates is presented in Table 2. In pigs the high prevalence found in 2017 (60%) and 2018 (50%) was followed by only a single isolate recovered in 2019 (3.3%), and then another increase in 2020 (33.3%). In humans, resistant isolates were only found in 2017 (40%), while a comparatively low rate of positivity in wastewater in 2017 (20%) and 2018 (10%) was followed by none in 2019, and a high prevalence in 2020 (60%). The majority (8/50: 16%) of MCRPE isolates recovered from wastewater were obtained from samples taken before biogas treatment, with only one isolate recovered in 2020 being from a sample taken after the biogas treatment plant (Table 1).

FIGURE 1

Figure 1. Comparison of the rate of mcr positive E. coli isolated from pigs, workers and the environment in four sample collection years (*significant difference; p ≤ 0.05).

TABLE 2

Table 2. Characterization of 33 colistin-resistant mcr positive E. coli isolates from different years and sources.

Identification of Plasmid-Mediated Colistin Resistance Genes

Of the colistin resistant isolates obtained in 2017, the mcr-1 gene was detected in eight of the pig isolates, while mcr-1 and mcr-3 were detected together in two of these, and mcr-3 alone in one pig isolate. At the same time, mcr-1 was detected in all four of the isolates from workers and in both the isolates from wastewater samples (Table 2). In 2018, after colistin withdrawal for one and a half years, mcr-1 was detected in three and mcr-3 in two of the five resistant isolates from pigs, and mcr-1 was found in the single resistant isolate from wastewater. In 2019 the single isolate from a breeder pig contained mcr-1. In 2020 mcr-1 positive E. coli isolates were found in all 5 piglets that had recent diarrhea symptoms and in wastewater samples (6/10).

Antimicrobial Susceptibility Determination

The antimicrobial resistance (AMR) profiles detected in the MCRPE isolates are shown in Figure 2. ESBL-producing E. coli were identified, and most MCRPE isolates from the first and second samplings were found to demonstrate extreme pan-drug resistance. Interestingly, besides colistin, the isolate from the positive pig sample in 2019 was phenotypically resistant only to ampicillin. On the other hand, the MCRPE isolates from the last sample collection in 2020 were resistant to aminoglycosides, ampicillin, and ceftiofur, and those antibiotics were used for individual treatments on the farm. The antibiogram results comparing isolates between the 4 sampling years are presented in Supplementary Table 1. High rates of resistance to ampicillin, gentamicin, and tetracycline were detected in almost all MCRPE isolates at each sampling time.

FIGURE 2

Figure 2. Comparison between resistance rates against 18 antimicrobials and extended-spectrum beta-lactamase (ESBL) in MCRPE isolates at four different sample collection times.

Various plasmid replicon types were detected among the MCRPE isolates (Table 3). All the mcr positive isolates from different sources contained more than one replicon type. The incompatibility group IncFIB and IncI type plasmids were found most commonly. Although a variety of plasmid types were detected in pigs in 2017 and 2018, there was a decrease in varieties of plasmid types in later sample collection years. For the conjugation assay, the donor E. coli transferred mcr-1 and mcr-3 genes (as confirmed by PCR) to recipient J53 strains with a frequency of 1.7–2 × 10⁻⁴. Phenotypic colistin resistance of transconjugants identified by broth microdilution showed MIC values of >4 μg/ml (Supplementary Table 2). The IncI, IncX and IncF plasmid types were predominantly detected in mcr-1 transconjugants, while the IncHI2 and IncF plasmid types were detected on mcr-3 transconjugants (data not showed).

TABLE 3

Table 3. Plasmid replicon types detected in 33 colistin-resistant E. coli among the three categories of samples at each sample collection time.

Virulence Gene Detection

Virulence gene detection was performed on all the 33 mcr positive E. coli isolates. Most of the isolates from pigs contained genes associated with ETEC strains (enterotoxin genes), with StaP and Stb being the most frequent pathotype found in 2017 (Table 2). One strain from a pig in 2017 showed a hybrid ETEC–EHEC genotype. Two of the four colistin-resistant E. coli recovered from farm workers in 2017 contained a combination of StaP and Stb genes. In contrast, the wastewater samples and the piglets' samples obtained after 2017 only contained the Stb enterotoxin gene.

Molecular Genotypic Characterizations

Thirty-four diverse PFGE patterns were obtained for the 65 MCRPE isolates from different sources (Figure 3). No dominant pulsotypes were responsible for mcr gene clonal carriage. Moreover, most of the strains from each sample collection time were dispersed on different branches of the dendrogram and were not closely related genetically. The pulsotypes of the MCRPE from humans were not clonally related to any of those from pigs or wastewater. Strains with high similarity (>80%) occurred rarely and were found mainly in the same set of pig or human samples from the same sampling year. Only the MCRPE strains from piglets and wastewater samples in 2020 showed high clonal relatedness, suggesting that MCRPE strains from the piglets with diarrhea had contaminated the wastewater.

FIGURE 3

Figure 3. Dendrogram generated from pulsed field gel electrophoresis analysis demonstrating the genetic relatedness among 65 mcr positive Escherichia coli strains (one to three isolates per positive sample) that were obtained from different sources at each sampling time.

MLST gave similar results to PFGE, with most isolates belonging to different STs (Supplementary Table 3). Isolates of the common E. coli clonal complex ST10 were detected in 2 pigs and one human sample on the first sampling, and in one wastewater sample on the last sample collection. Isolates belonging to ST 641 were detected in 2 pigs and one wastewater sample in 2018. In 2020, MCRPE isolates belonging to ST3345 (n = 3 in wastewater, n = 2 in pigs) and ST 5218 (n = 2 in pigs) were commonly detected.

Discussion

The geographical distribution and characterization of colistin-resistant E. coli on large-scale pig farms across Thailand has been reported previously (19). The current longitudinal study investigated the persistence and diversity of mcr positive E. coli on a selected commercial Thai pig farm following withdrawal of prophylactic colistin usage. According to the farm history, batches of piglets previously were consistently prescribed colistin up until the time that it was withdrawn from use at the start of 2017. In searching for potential changes in resistance to colistin after its withdrawal, this study focused on examining colistin resistance in E. coli from young sows and their suckling pigs, as well as from wastewater. Prior to 2017 colistin was mainly used for controlling E. coli in suckling pigs in their first 3 weeks of life, so it seemed logical to target this bacterial species and this age group when looking for ongoing resistance. In addition, sows were examined since piglets become colonized by oral exposure from the fecal microbiota of their mothers. The sows were exposed to colistin prior to 2017 and might be persistently colonized and hence transfer resistant bacteria to their piglets. The piglets themselves were destined for slaughter by around 5–6 months of age, and so by definition could not be involved in direct transmission in following years. Accordingly, more sows than piglets were sampled to determine whether they still represented a potential long-term reservoir of MCRPE infection. Wastewater also was sampled, as wastewater tanks on pig farms serve as hotspots for accumulating resistant bacteria since they are composed of pooled fecal discharge from large numbers of pigs housed in the same area. Inclusion of this material in the study increased the likelihood of detecting MCRPE. Although relatively small numbers of samples were examined at each sampling time, they were sufficient to confirm the presence of MCRPE throughout the study.

Even following colistin withdrawal for 21 months, MCRPE that were carrying mcr genes were still quite commonly found in pigs, indicating that this period is insufficient to have a significant effect on reducing the presence of colistin resistant bacteria after the drug's withdrawal. The presence of mcr positive E. coli in the feces of farm workers in mid-2017 is a matter of considerable concern. Bacteria from animals can be transmitted to humans either directly or through food or the environment, and then may transfer resistance genes to pathogenic bacteria that infect humans (43). Farm and food chain workers are likely to be exposed to resistant bacteria throughout the pig production cycle (44). Moreover, Stb and StaP virulence genes were found in MCRPE isolates from pigs and in two of the workers. These enterotoxin genes are linked to neonatal or postweaning diarrhea in pigs, but bacteria carrying the genes also can be shed in feces from healthy animals (45). The Stb enterotoxin is commonly found in E. coli strains from pigs but is rarely found in humans and is not associated with diarrhea in humans (46, 47). Therefore, these findings suggested that subclinical ETEC carriers can be found at various stages in the pig production cycle and may represent a source of transmission to humans. Even though isolates of identical genetic types were not found in humans and pigs, our results highlight possible transmission of mcr genes from bacteria infecting livestock to isolates that are present in humans and in the environment. The failure to recover MCRPE from human samples after 2017 may be associated with reduced exposure to colistin and/or to MRCPE from pigs and the environment. One possibility is that following identification of MCRPE in the workers in 2017 these individuals took greater care of their hygiene to reduce their exposure to MCRPE of pig origin.

Most pigs were still colonized with colistin-resistant E. coli when sampled 21 months after colistin withdrawal; however, by the third year there was a sharp decline in carriage by pigs and neither workers nor wastewater samples were positive for MCRPE. Therefore, the national ban on prophylactic use of colistin was highly likely to have been beneficial for controlling the emergence and dissemination of mcr-1 on this and on other pig farms. Nevertheless, some pigs and wastewater still contained MRCPE when samples 3.5 years after the withdrawal of colistin. Pigs reared for meat production are only kept for around 5–6 months before slaughter, although breeder pigs are retained for up to 3–4 years. Presumably transmission cycles of MCRPE between batches of pigs that are selected for meat production or breeding, and/or exposure to contaminated environments allowed them to remain for at least 3.5 years. A more extended study is required to determine for how long this carriage my persist. The long duration of persistence that was identified contrasted with a previous report from Britain, where mcr-1 was undetectable in isolates from pigs after the cessation of colistin use for approximately 20 months (7). The reason for the re-occurrence and increase in numbers of pigs shedding MCRPE and in isolates recovered from wastewater in the last year of the current study is unclear. These colonized pigs had shown diarrhea symptoms prior to sampling and had been given therapeutic antimicrobial treatments, unlike the situation in previous batches sampled in earlier years. A possible explanation for the re-occurrence without selective pressure applied from colistin exposure may be the existence of cross-resistance between colistin and other therapeutic antibiotics used in the piglets. A similar phenomenon was reported in previous studies where colistin resistance was found when other antimicrobials such as quinolones or cephalosporins were used in livestock farms (48, 49). However, more complete genomic characterization of the MCRPE isolates involved is required to investigate possible reason for this correlation. Nevertheless, these results are of concern because short-term β-lactam (ceftiofur) or gentamicin use in animals may select for mcr-1 in E. coli and maintain persistence on farms.

From the antibiotic susceptibility testing, some of the mcr positive E. coli isolates were found to be ESBL producers and showed extreme pan-drug resistance. A larger number of E. coli isolates with ESBL were observed in the samples from 2018 compared to the first sampling time. In Thailand, the application of antimicrobials in pig farms varies according to the management system and geographical area. In the central area of Thailand, the antimicrobials that are mainly used are colistin, cephalosporins, tiamulin, amoxicillin, tilmicosin, aminoglycosides (gentamycin), and oxytetracycline (50). The use of other antimicrobials during the production cycle of pigs could co-select for colistin resistance (51, 52). Resistances to other potential agents like heavy metals or biocides that may be linked with antibiotics resistance genes also are a matter for concern.

In the conjugation experiment, MCRPE recovered from pigs without selective pressures from colistin use showed a high transfer frequency. Moreover, various replicon types were found in the colistin-resistant E. coli isolates. According to previous reports, mcr-1 and mcr-3 genes have been found on IncI, IncHI2, and IncX4 plasmids (53). Likewise, mcr-1 was predominantly harbored on the IncX4 plasmid in isolates from healthy human beings in China (54). Different AMR genes can be located on the same plasmid or on different plasmids within the same bacterial host, and these represent multidrug resistant clones. Plasmids encoding the mcr genes, which co-exist with other antimicrobial resistance genes, are a problem for public health. To date, the majority of mcr genes have been identified in various plasmid types and are able to locate and/or transfer with other resistance genes by conjugation (55).

A large number of PFGE pulsotypes were observed among the mcr positive isolates. Therefore, no epidemic strains were dominant on the farm over time, and the mcr genes found in E. coli isolates were mainly plasmid-borne. A high diversity of MCRPE isolates from different hosts also was observed in a study from China (56). Similarly, in a Dutch study where ESBL positive E. coli from animals and humans were examined, ESBL transmission did not involve strain transfer but rather plasmid transfer by identical plasmids of the IncI and IncK types (57). Nevertheless, in our study some clonal relatedness was found in MCRPE from piglets and wastewater samples at the last sampling. In this case the resistant bacteria from pigs were likely to be the primary source of mcr genes contaminating wastewater. Thus, despite moderate persistence of mcr genes in pigs and low-level environmental dissemination in tested wastewater, the distribution of diverse strains with virulence potential from different niches across years is worrisome. Genes from these mcr-1 and mcr-3 positive isolates might be transferred to other sources and/or other pathogens.

In this study, E. coli carrying mcr genes were recovered, but with a gradual decline over 3.5 years following cessation of colistin use. Hence, banning colistin for prophylaxis use was effective for reducing the emergence and dissemination of mcr-1 in pigs and the pig farm environment. However, even in the absence of selective pressure exerted by colistin use, the application of other antimicrobials during the production cycle might co-select indirectly for the mcr genes and favor their spread. This study provides an initial insight into the reduction in dissemination of colistin resistant E. coli from pigs and the farm environment. Further long-term genomic investigations are necessary to improve understanding and control of MCRPE and colistin-resistance in the pig industry.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

Ethical review and approval was not required for the animal study because All fecal samples were submitted from veterinarians in pig industrial field to the veterinary diagnostic laboratory as the annual surveillance. However, the biohazard execution control was approved by the Institutional Biosafety Committee of the Faculty of Veterinary Science, Chulalongkorn University (IBC 2031011).

Author Contributions

NK contributed to conception and design, critical revision, analysis and interpretation of data, and drafting of manuscript. KL contributed to conception and design, analysis, and interpretation of data. WN performed analysis and interpretation of data. TP organized acquisition of data. DH contributed to critical revision, analysis and interpretation of data, and drafting of manuscript. NP organized study conception and design and contributed to critical revision, analysis, and interpretation of data. All authors contributed to manuscript revision, read, and approved the submitted version.

Funding

This study was supported financially by Royal Golden Jubilee Ph.D. (RGJPHD) program, from the Agricultural Research Development Agency (ARDA) [CRP6205031110], the CHE-TRF Senior Research Fund (RTA6280013), Thailand Science Research and Innovation, and Pathogen Bank, Faculty of Veterinary Science, Chulalongkorn University, Bangkok, Thailand.

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.

Acknowledgments

The 90th Anniversary of Chulalongkorn University Scholarship.

Supplementary Material

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

References

1. Falagas ME, Kasiakou SK. Colistin: the revival of polymyxins for the management of multidrug-resistant gram-negative bacterial infections. Clin Infect Dis. (2005) 40:1333–41. doi: 10.1086/429323

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Liu YY, Wang Y, Walsh TR, Yi LX, Zhang R, Spencer J, et al. 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. (2016) 16:161–8. doi: 10.1016/S1473-3099(15)00424-7

PubMed Abstract | CrossRef Full Text | Google Scholar

3. AbuOun M, Stubberfield EJ, Duggett NA, Kirchner M, Dormer L, Nunez-Garcia J, et al. mcr-1 and mcr-2 (mcr-6.1) variant genes identified in Moraxella species isolated from pigs in Great Britain from 2014 to 2015. J Antimicrob Chemother. (2018) 73:2904. doi: 10.1093/jac/dky272

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Wang X, Wang Y, Zhou Y, Li J, Yin W, Wang S, et al. Emergence of a novel mobile colistin resistance gene, mcr-8, in NDM-producing Klebsiella pneumoniae. Emerg Microbes Infect. (2018) 7:122. doi: 10.1038/s41426-018-0124-z

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Yang YQ, Li YX, Lei CW, Zhang AY, Wang HN. Novel plasmid-mediated colistin resistance gene mcr-7.1 in Klebsiella pneumoniae. J Antimicrob Chemother. (2018) 73:1791–5. doi: 10.1093/jac/dky111

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Carroll LM, Gaballa A, Guldimann C, Sullivan G, Henderson LO, Wiedmann M. Identification of novel mobilized colistin resistance gene mcr-9 in a multidrug-resistant, colistin-susceptible Salmonella enterica Serotype Typhimurium isolate. mBio. (2019) 10:e00853-19. doi: 10.1128/mBio.00853-19

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Duggett NA, Randall LP, Horton RA, Lemma F, Kirchner M, Nunez-Garcia J, et al. Molecular epidemiology of isolates with multiple mcr plasmids from a pig farm in Great Britain: the effects of colistin withdrawal in the short and long term. J Antimicrob Chemother. (2018) 73:3025–33. doi: 10.1093/jac/dky292

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Wang C, Feng Y, Liu L, Wei L, Kang M, Zong Z. Identification of novel mobile colistin resistance gene mcr-10. Emerg Microbes Infect. (2020) 9:508–16. doi: 10.1080/22221751.2020.1732231

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Andrade FF, Silva D, Rodrigues A, Pina-Vaz C. Colistin update on its mechanism of action and resistance, present and future challenges. Microorganisms. (2020) 8:1716. doi: 10.3390/microorganisms8111716

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Soliman AM, Ramadan H, Zarad H, Sugawara Y, Yu L, Sugai M, et al. Coproduction of Tet(X7) conferring high-level tigecycline resistance, fosfomycin fosa4, and colistin mcr-1.1 in Escherichia coli strains from chickens in Egypt. Antimicrob Agents Chemother. (2021) 65:e02084-20. doi: 10.1128/AAC.02084-20

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Irrgang A, Roschanski N, Tenhagen BA, Grobbel M, Skladnikiewicz-Ziemer T, Thomas K, et al. Prevalence of mcr-1 in E. coli from livestock and food in Germany, 2010-2015. PLoS ONE. (2016) 11:e0159863. doi: 10.1371/journal.pone.0159863

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Tartor YH, Gharieb RMA, Abd El-Aziz NK, El Damaty HM, Enany S, Khalifa E, et al. Virulence determinants and plasmid-mediated colistin resistance mcr genes in gram-negative bacteria isolated from bovine milk. Front Cell Infect Microbiol. (2021) 11:761417. doi: 10.3389/fcimb.2021.761417

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Prim N, Rivera A, Rodriguez-Navarro J, Espanol M, Turbau M, Coll P, et al. Detection of mcr-1 colistin resistance gene in polyclonal Escherichia coli isolates in Barcelona, Spain, 2012 to 2015. Euro Surveill. (2016) 21:30183. doi: 10.2807/1560-7917.ES.2016.21.13.30183

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Huang H, Dong N, Shu L, Lu J, Sun Q, Chan EW, et al. Colistin-resistance gene mcr in clinical carbapenem-resistant Enterobacteriaceae strains in China, 2014-2019. Emerg Microbes Infect. (2020) 9:237–45. doi: 10.1080/22221751.2020.1717380

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Aguirre L, Vidal A, Seminati C, Tello M, Redondo N, Darwich L, et al. Antimicrobial resistance profile and prevalence of extended-spectrum beta-lactamases (ESBL), AmpC beta-lactamases and colistin resistance (mcr) genes in Escherichia coli from swine between 1999 and 2018. Porcine Health Manag. (2020) 6:8. doi: 10.1186/s40813-020-00146-2

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Rhouma M, Beaudry F, Theriault W, Letellier A. Colistin in pig production: chemistry, mechanism of antibacterial action, microbial resistance emergence, and one health perspectives. Front Microbiol. (2016) 7:1789. doi: 10.3389/fmicb.2016.01789

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Walsh TR, Wu Y. China bans colistin as a feed additive for animals. Lancet Infect Dis. (2016) 16:1102. doi: 10.1016/S1473-3099(16)30329-2

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Khine NO, Lugsomya K, Kaewgun B, Honhanrob L, Pairojrit P, Jermprasert S, et al. Multidrug resistance and virulence factors of Escherichia coli harboring plasmid-mediated colistin resistance: mcr-1 and mcr-3 genes in contracted pig farms in Thailand. Front Vet Sci. (2020) 7:582899. doi: 10.3389/fvets.2020.582899

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Lay KK, Jeamsripong S, Sunn KP, Angkititrakul S, Prathan R, Srisanga S, et al. Colistin resistance and ESBL production in Salmonella and Escherichia coli from pigs and pork in the Thailand, Cambodia, Lao PDR, and Myanmar border area. Antibiotics. (2021) 10:657. doi: 10.3390/antibiotics10060657

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Xia X, Wang Z, Fu Y, Du XD, Gao B, Zhou Y, et al. Association of colistin residues and manure treatment with the abundance of mcr-1 gene in swine feedlots. Environ Int. (2019) 127:361–70. doi: 10.1016/j.envint.2019.03.061

PubMed Abstract | CrossRef Full Text | Google Scholar

22. de Toro M, Garcillaon-Barcia MP, De La Cruz F. Plasmid diversity and adaptation analyzed by massive sequencing of Escherichia coli plasmids. Microbiol Spectr. (2014) 2:2–6. doi: 10.1128/microbiolspec.PLAS-0031-2014

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Randall LP, Horton RA, Lemma F, Martelli F, Duggett NAD, Smith RP, et al. Longitudinal study on the occurrence in pigs of colistin-resistant Escherichia coli carrying mcr-1 following the cessation of use of colistin. J Appl Microbiol. (2018) 125:596–608. doi: 10.1111/jam.13907

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Yates CM, Shaw DJ, Roe AJ, Woolhouse ME, Amyes SG. Enhancement of bacterial competitive fitness by apramycin resistance plasmids from non-pathogenic Escherichia coli. Biol Lett. (2006) 2:463–5. doi: 10.1098/rsbl.2006.0478

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Brolund A, Sundqvist M, Kahlmeter G, Grape M. Molecular characterisation of trimethoprim resistance in Escherichia coli and Klebsiella pneumoniae during a two year intervention on trimethoprim use. PLoS ONE. (2010) 5:e9233. doi: 10.1371/journal.pone.0009233

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Sundqvist M, Geli P, Andersson DI, Sjolund-Karlsson M, Runehagen A, Cars H, et al. Little evidence for reversibility of trimethoprim resistance after a drastic reduction in trimethoprim use. J Antimicrob Chemother. (2010) 65:350–60. doi: 10.1093/jac/dkp387

PubMed Abstract | CrossRef Full Text | Google Scholar

27. HACH. Water Analysis Guide (2013).

Google Scholar

28. Singhal N, Kumar M, Kanaujia PK, Virdi JS. MALDI-TOF mass spectrometry: an emerging technology for microbial identification and diagnosis. Front Microbiol. (2015) 6:791. doi: 10.3389/fmicb.2015.00791

PubMed Abstract | CrossRef Full Text | Google Scholar

29. CLSI. Clinical and Laboratory Standards Institute. Performance Standards for Antimicrobial Susceptibility Testing. Wayne, PA (2020).

Google Scholar

30. Plumb. Plumb's Veterinary Drug Handbook. Indianapolis, IN: Wiley-Blackwell Press (2015).

Google Scholar

31. FDA. Microbiology Data for Systemic Antibacterial Drugs — Development, Analysis, and Presentation Guidance for Industry. Food and Drug Administration Center for Drug Evaluation and Research (CDER) (2018).

Google Scholar

32. CLSI. Clinical and Laboratory Standards Institute. Performance Standards for Antimicrobial Susceptibility Testing. Wayne, PA (2017).

Google Scholar

33. EUCAST. The European Committee on Antimicrobial Susceptibility Testing: Breakpoint Tables for Interpretation of MICs (2018).

Google Scholar

34. Rebelo AR, Bortolaia V, Kjeldgaard JS, Pedersen SK, Leekitcharoenphon P, Hansen IM, et al. 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. (2018) 23:17-00672. doi: 10.2807/1560-7917.ES.2018.23.6.17-00672

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Lugsomya K, Chanchaithong P, Tribudharat C, Thanawanh N, Niyomtham W, Prapasarakul N. Prevalence of plasmid-mediated colistin resistance genes (mcr-1) in commensal Escherichia coli from fattening pigs in Thailand. Thai J Vet Med. (2016) 46:2016.

Google Scholar

36. Carattoli A, Bertini A, Villa L, Falbo V, Hopkins KL, Threlfall EJ. Identification of plasmids by PCR-based replicon typing. J Microbiol Methods. (2005) 63:219–28. doi: 10.1016/j.mimet.2005.03.018

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Casey TA, Bosworth BT. Design and evaluation of a multiplex polymerase chain reaction assay for the simultaneous identification of genes for nine different virulence factors associated with Escherichia coli that cause diarrhea and edema disease in swine. J Vet Diagn Invest. (2009) 21:25–30. doi: 10.1177/104063870902100104

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Prapasarakul N, Tummaruk P, Niyomtum W, Tripipat T, Serichantalergs O. Virulence genes and antimicrobial susceptibilities of hemolytic and nonhemolytic Escherichia coli isolated from post-weaning piglets in central Thailand. J Vet Med Sci. (2010) 72:1603–8. doi: 10.1292/jvms.10-0124

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Gray KJ, Wilson LK, Phiri A, Corkill JE, French N, Hart CA. Identification and characterization of ceftriaxone resistance and extended-spectrum beta-lactamases in Malawian bacteraemic Enterobacteriaceae. J Antimicrob Chemother. (2006) 57:661–5. doi: 10.1093/jac/dkl037

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Apiwatsiri P, Pupa P, Yindee J, Niyomtham W, Sirichokchatchawan W, Lugsomya K, et al. Anticonjugation and antibiofilm evaluation of probiotic strains Lactobacillus plantarum 22F, 25F, and Pediococcus acidilactici 72N against Escherichia coli harboring mcr-1 gene. Front Vet Sci. (2021) 8:614439. doi: 10.3389/fvets.2021.614439

PubMed Abstract | CrossRef Full Text | Google Scholar

41. CDC. Centers for Disease Control and Prevention, Standard Operating Procedure for Pulsenet pfge of Escherichia coli O157:H7, Escherichia coli non-O157 (STEC), Salmonella serotypes, Shigella sonnei and Shigella flexneri. Centers for Disease Control and Prevention (2013).

Google Scholar

42. Adiri RS, Gophna U, Ron EZ. Multilocus sequence typing (MLST) of Escherichia coli O78 strains. FEMS Microbiol Lett. (2003) 222:199–203. doi: 10.1016/S0378-1097(03)00295-7

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Archawakulathep A, Kim CTT, Meunsene D, Handijatno D, Hassim HA, Rovira HR, et al. Perspectives on antimicrobial resistance in livestock and livestock products in ASEAN countries. Thai J Vet Med. (2014) 44:5–13. Available online at: https://he01.tci-thaijo.org/index.php/tjvm/article/view/17308

Google Scholar

44. You Y, Leahy K, Resnick C, Howard T, Carroll KC, Silbergeld EK. Exposure to pathogens among workers in a poultry slaughter and processing plant. Am J Industrial Med. (2016) 59:453–64. doi: 10.1002/ajim.22594

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Moredo FA, Pineyro PE, Marquez GC, Sanz M, Colello R, Etcheverria A, et al. Enterotoxigenic Escherichia coli subclinical infection in pigs: bacteriological and genotypic characterization and antimicrobial resistance profiles. Foodborne Pathog Dis. (2015) 12:704–11. doi: 10.1089/fpd.2015.1959

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Echeverria P, Seriwatana J, Patamaroj U, Moseley SL, McFarland A, Chityothin O, et al. Prevalence of heat-stable II enterotoxigenic Escherichia coli in pigs, water, and people at farms in Thailand as determined by DNA hybridization. J Clin Microbiol. (1984) 19:489–91. doi: 10.1128/jcm.19.4.489-491.1984

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Casey TA, Herring CJ, Schneider RA, Bosworth BT, Whipp SC. Expression of heat-stable enterotoxin STb by adherent Escherichia coli is not sufficient to cause severe diarrhea in neonatal pigs. Infect Immun. (1998) 66:1270–2. doi: 10.1128/IAI.66.3.1270-1272.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Nguyen NT, Nguyen HM, Nguyen CV, Nguyen TV, Nguyen MT, Thai HQ, et al. Use of colistin and other critical antimicrobials on pig and chicken farms in Southern Vietnam and its association with resistance in commensal Escherichia coli bacteria. Appl Environ Microbiol. (2016) 82:3727–35. doi: 10.1128/AEM.00337-16

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Cao YP, Lin QQ, He WY, Wang J, Yi MY, Lv LC, et al. Co-selection may explain the unexpectedly high prevalence of plasmid-mediated colistin resistance gene mcr-1 in a Chinese broiler farm. Zool Res. (2020) 41:569–75. doi: 10.24272/j.issn.2095-8137.2020.131

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Pokhrel N, Kasemsuwan S, Goutard F, Boonsoongnern A, Chumsing S, Tulayakul P. Practices and factors influencing the use of antibiotics in swine farms in Central region of Thailand. J Kasetsart Vet. (2019) 29:33–55. Available online at: https://kukr2.lib.ku.ac.th/kukr_es/kukr/search_detail/result/408382

PubMed Abstract | Google Scholar

51. Migura-Garcia L, González-López JJ, Martinez-Urtaza J, Aguirre Sánchez J, Moreno-Mingorance A, Perez de Rozas A, et al. mcr-colistin resistance genes mobilized by IncX4, IncHI2, and IncI2 plasmids in Escherichia coli of pigs and White Stork in Spain. Front Microbiol. (2020) 10:3072. doi: 10.3389/fmicb.2019.03072

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Vines J, Cusco A, Napp S, Alvarez J, Saez-Llorente JL, Rosas-Rodoreda M, et al. Transmission of similar mcr-1 Carrying plasmids among different Escherichia coli lineages isolated from livestock and the farmer. Antibiotics. (2021) 10:313. doi: 10.3390/antibiotics10030313

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Kieffer N, Nordmann P, Moreno AM, Zanolli Moreno L, Chaby R, Breton A, et al. Genetic and functional characterization of an MCR-3-Like enzyme-producing Escherichia coli isolate recovered from swine in Brazil. Antimicrob Agents Chemother. (2018) 62:e00278–18. doi: 10.1128/AAC.00278-18

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Shen Y, Zhou H, Xu J, Wang Y, Zhang Q, Walsh TR, et al. Anthropogenic and environmental factors associated with high incidence of mcr-1 carriage in humans across China. Nat Microbiol. (2018) 3:1054–62. doi: 10.1038/s41564-018-0205-8

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Rhouma M, Theriault W, Rabhi N, Duchaine C, Quessy S, Fravalo P. First identification of mcr-1/mcr-2 genes in the fecal microbiota of Canadian commercial pigs during the growing and finishing period. Vet Med. (2019) 10:65–7. doi: 10.2147/VMRR.S202331

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Shen C, Zhong L-L, Yang Y, Doi Y, Paterson DL, Stoesser N, et al. Dynamics of mcr-1 prevalence and mcr-1-positive Escherichia coli after the cessation of colistin use as a feed additive for animals in China: a prospective cross-sectional and whole genome sequencing-based molecular epidemiological study. Lancet Microbe. (2020) 1:e34–43. doi: 10.1016/S2666-5247(20)30005-7

CrossRef Full Text | Google Scholar

57. de Been M, Lanza VF, de Toro M, Scharringa J, Dohmen W, Du Y, et al. Dissemination of cephalosporin resistance genes between Escherichia coli strains from farm animals and humans by specific plasmid lineages. PLoS Genet. (2014) 10:e1004776. doi: 10.1371/journal.pgen.1004776

PubMed Abstract | CrossRef Full Text | Google Scholar