- 1Unit of Microbiology and Clinical Microbiology, Department of Clinical-Surgical, Diagnostic and Pediatric Sciences, University of Pavia, Pavia, Italy
- 2Department of Biomedical and Clinical Sciences, Romeo and Enrica Invernizzi Pediatric Research Center, University of Milan, Milan, Italy
- 3Faculty of Medicine, Charles University, Plzen, Czechia
- 4Molecular Virology Unit, Microbiology and Virology Department, Fondazione IRCCS Policlinico San Matteo, Pavia, Italy
- 5Department of Earth and Environment Sciences, University of Pavia, Pavia, Italy
To evaluate the water compartment antibiotic-resistance contamination rates, 11 wells, five streams, and four treatment plants located in the Oltrepò Pavese area were screened for the presence of third generation cephalosporins resistant Gram-negative bacteria. Enterobacteriaceae were also characterized for the Extended-Spectrum-β-Lactamases (ESBLs), carbapenemases, and mcr-1 genes presence. From December 2014 to November 2015, 246 water samples were filtered, plated on Plate Count Agar, MacConkey Agar, and MacConkey Agar with cefotaxime. Isolates were species identified using AutoSCAN-4-System and ESBLs, carbapenemases, and colistin resistance determinants were characterized by PCR, sequencing, and microarray. Plasmid conjugative transfer experiments, PCR-based Replicon typing, Pulsed-Field Gel Electrophoresis, Multi-Locus-Sequence-Typing, and in-silico plasmid characterization were performed. A total of 132 enterobacteria isolates grew on MacConkey agar with cefotaxime: 82 (62.1%) were obtained from streams, 41 (31.1%) from treatment plants, and 9 (6.8%) from wells. Thirty out of 132 (22.7%) isolates, mainly belonging to Escherichia coli (n = 15) species, showed a synergic effect with piperacillin-tazobactam. A single ESBL gene of blaCTX−M-type was identified in 19/30 isolates. In further two E. coli strains, a blaCTX−M−1 gene co-existed with a blaSHV-type ESBL determinant. A blaSHV−12 gene was detected in two isolates of E. coli (n = 1) and Klebsiella oxytoca (n = 1), while any ESBL determinant was ascertained in seven Yersinia enterocolitica strains. A blaDHA-type gene was detected in a cefoxitin resistant Y. enterocolitica from a stream. Interestingly, two Klebsiella pneumoniae strains of ST307 and ST258, collected from a well and a wastewater treatment plant, resulted KPC-2, and KPC-3 producers, respectively. Moreover, we report the first detection of mcr-1.2 ST10 E. coli on a conjugative IncX4 plasmid (33.303 bp in size) from a stream of Oltrepò Pavese (Northern Italy). Both ESBLs E. coli and ESBLs/carbapenemase-producing K. pneumoniae strains showed clonal heterogeneity by Pulsed-Field Gel Electrophoresis and Multi-Locus-Sequence-Typing. During one-year study and taking in account the whole Gram-negative bacterial population, an average percentage of cefotaxime resistance of 69, 32, and 10.3% has been obtained for the wastewater treatment plants, streams, and wells, respectively. These results, of concern for public health, highlight the need to improve hygienic measures to reduce the load of discharged bacteria with emerging resistance mechanisms.
Introduction
Antibiotic resistance, in particular to third generation cephalosporins (3GCs) and carbapenems, threatens healthcare globally. Drug resistance has traditionally been viewed as a clinical problem, but recently natural ecosystems have been recognized as an important reservoir of antibiotic resistance genes (ARGs) (Berglund, 2015).
In aquatic environments the prevalence of antibiotic-resistant bacteria, which may originate from anthropogenic sources such as hospital and municipal effluents, is constantly growing (Baquero et al., 2008; Bouki et al., 2013). The large amounts of antibiotics or their active metabolites released into wastewater with treated urine and feces increase the selective pressure in bacterial populations, allowing the development of antibiotic-resistant microbes' generations (Davies and Davies, 2010).
Furthermore, the surface water could act as resistance hotspots where ARGs disseminate favored by bacteriophages or integrons and new resistant strains are created by horizontal gene transfer (Berglund, 2015; Colombo et al., 2017).
β-lactamase genes have been identified in bacteria isolated from both surface waters and Wastewater Treatment Plants (WWTPs) (Schwartz et al., 2003; Poppe et al., 2006; Bouki et al., 2013). The β-lactamases now include >2,000 naturally occurring amino acid sequences. Some of the clinically most important of these are the Class A penicillinases, the Extended-Spectrum β-Lactamases (ESBLs) (TEM, SHV, VEB, or CTX-M the most spread), the AmpC cephalosporinases (usually CMY, FOX, DHA, ACT, or MOX), and the carbapenem-hydrolyzing enzymes in both the serine (i.e., KPC, IMI, SME, or GES) and metallo-enzyme groups (mainly NDM, VIM, or IMP) (Bonomo, 2017). Because of the versatility of these enzymes to evolve as new β-lactams are used therapeutically, combined approaches to antimicrobial therapy may be required (i.e., carbapenems plus colistin). These genes, frequently located on mobile genetic elements, often coexist with other resistance determinants working on other classes of antibiotics such as aminoglycosides and fluoroquinolones (Marti et al., 2014).
A priority to be currently aware and vigilant of is the worldwide recently emerged plasmid-mediated resistance to colistin, another “last-resort” antibiotic used for the treatment of critical infections caused by MDR Gram-negative pathogens (Olaitan et al., 2016).
The Oltrepò Pavese Plan is located in the province of Pavia with an area of ~1,097 km2 and a population of 146,579 inhabitants; its name comes to the peculiarity of being south of the river Po, in the full northern Apennines. The Oltrepò Pavese area supports many agricultural and industrial activities causing pollution in the shallower aquifers; for this reason, it represents a good model to evaluate the role of aquatic environment as a potential reservoir for spread and evolution of ARGs and their vectors. The routes by which humans may come in contact with these bacteria include the consumption of crops grown by contaminated sludge used as fertilizer, and/or drinking of water drawn from contaminated ground or surface water. When these resistant bacteria enter humans, they have the opportunity to spread their ARGs to the human microbiome (Wellington et al., 2013).
Only few European studies on surface waters are available; furthermore, due to the lack of local and/or large-scale epidemiological studies, our knowledge on the environmental antibiotic-resistance reservoirs in Italy is still very poor (Perilli et al., 2013; Zanotto et al., 2016).
The purpose of this study was (i) to investigate the occurrence of cefotaxime resistant Enterobacteriaceae in Oltrepò Pavese water ecosystem (ii) to determine Gram-negative bacterial counts and percentage of 3GCs resistance, during all the study period, (iii) to characterize the underlying determinants and to assess the clonal features of the most alarming isolates.
Materials and Methods
Design Study and Sample Collection
The aim of the present 1-year study was to monitor the Gram-negative Oltrepò Pavese area bacterial counts during different months per sampling site, and the corresponding 3GCs resistance percentages. Water samples from wells, streams, and WWTPs of the Oltrepò Pavese area were evaluated for the presence of ESBLs/AmpC, carbapenemases, and mcr-1 genes in Enterobacteriaceae.
The sampling sites are located at south of Po river in the Italian Northern Apennines, in the zone between Voghera, Stradella, and Staffora Valley near Varzi.
This area is characterized by the presence of a sandy-gravelly alluvial deposits covering the marine sediments, which are essentially impermeable (Pilla et al., 2007). Local well have a recharging time of 6 months on average and their structure does not provide a complete isolation from the surface water.
Sampling was carried out once a month—during the period December 2014–November 2015. In the study period, a total of 246 water samples were collected. Each of the single stream sampling site (1T−11T) was visited from 9 up to 11 times, while the 11 wells (1P−13P) were sampled from 10 up to 12 times. The four WWTPs included in the study were less frequently inspected (from three to six times), due to the need for a special authorization to get samples.
All the sampled wells are part of the Po Plan Public-Supply System, serving freshwater for human consumption in the area. The outflow of four selected WWTPs, located upstream rivers of the same area, was screened starting from April 2015 only (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/272364/fmicb-08-02232-HTML-r1/image_m/fmicb-08-02232-g001.jpg)
Figure 1. Map of the Oltrepò Pavese Plain showing the sampling sites. P, Well; T, Stream; TP, Wastewater Treatment Plant.
All water samples were collected in aseptic plastic bottles, stored at 4°C, and analyzed within 24 h after harvesting.
The stream sampling sites have been chosen with different location, upstream, and downstream urban centers, including factories, hospitals, and Long Term Care Facilities (LTCFs).
Filtration Method
A water volume of 100 ml from wells and 1 ml from up/down streams or at the outflow of treatment plants were passed through 0.45 μm-pore size membranes by the filtration method. Filter membranes were placed on Plate Count Agar (PCA), MacConkey Agar (MCA), and selective MCA containing 8 mg/l of CTX (MCA+CTX), to select the potential ESBL/AmpC producers. The bacterial count was estimated after 24 h incubation of the plates at 37°C.
The 3GCs resistance percentage (% R) was calculated by directly comparing the colonies count on the antibiotic plate (MCA+CTX) with the corresponding count on the control plate (MCA) (Watkinson et al., 2007).
Bacterial Identification and Susceptibility Testing
Each colony grown on selective MCA+CTX and distinguishable for the unique morphotype, was identified at species level. Species identification and susceptibility testing were performed using the semi-automated system MicroScan autoSCAN-4 (Beckman Coulter). Colistin (CO) MIC was confirmed by broth microdilution method. Susceptibility results were interpreted according to the EUCAST 2015 (http://www.eucast.org) clinical guidelines. Escherichia coli ATCC 25922 was used as control strain.
Phenotypic ESBL detection was performed by the use of chromogenic media ChromArt ESBL (Biolife Italiana S.r.l., Milano, Italy) and the Double-Disc Synergy (DD) test (Jarlier et al., 1988), using piperacillin-tazobactam, CTX, cefepime, ceftazidime, and aztreonam. A phenotype of both cefoxitin and 3GCs resistance without synergistic activity after DD, was regarded as AmpC-production marker. The strains were phenotypically confirmed for the carbapenemases production by the KPC/MBL Confirm kit (Rosco Diagnostic).
Molecular Characterization of Resistance Genes
The genomic DNA was extracted using NucleoSpin Tissue (Macherey-Nagel) kit. PCR was used to detect the presence of ESBLs, carbapenemases-, and MCR-1-encoding genes using primers and the annealing temperatures described in Table 1 (Rasheed et al., 1997; Yigit et al., 2001; Perilli et al., 2002; Pagani et al., 2003; Eckert et al., 2004; Liu et al., 2016). In the case of ESBLs and carbapenemases, the above annealing temperatures were always preceded by an initial 95°C × 5 min denaturation step. The 35 cycling conditions comprised denaturation and extension steps at 94°C × 3 min and 72°C × 1 min respectively; the final extension was at 72°C 5 min. mcr-1 PCR approach was performed using an initial 94°C × 5 min denaturation step; the 35 cycling conditions comprised the denaturation and extension steps at 94°C and 72°C × 30 s, respectively; the final extension was at 72°C 5 min. PCR products were purified using the quantum PrepPCR Kleen Spin Columns kit (ThermoFisher Scientific) and subjected to double-strand sequencing using the automatic sequencer ABI PRISM 3100 genetic analyzer DNA Sequencer (Applied Biosystems, Foster City, CA, USA) and the BigDye Terminator v1.1 Cycle Sequencing kit (Applied Biosystems, Foster City, CA, USA).
The sequences were analyzed according to the BLAST software program (http://blast.ncbi.nlm.nih.gov/Blast.cgi).
In the case of Y. enterocolitica Check-MDR CT103 XL array (Check points Health B.V., Wageningen, The Netherlands) has been used to investigate the bla genes content.
Conjugation Transfer Experiments and Plasmid Characterization
Conjugation transfer of resistance determinants was performed in liquid medium using the E. coli K12 strain J62 (F−, pro, his, trp, lac, Smr) and J53 (F−, met, pro, Rifr) as recipients. To evaluate the transferability of CTX- or CO- resistance, the E. coli transconjugants were screened using MCA plates supplemented with streptomycin (1,000 mg/l) or rifampicin (100 mg/l) plus CTX (8 mg/l) or CO (2 mg/l), respectively. Azide-resistant E. coli J53AzR strain was used as a recipient for carbapenem-resistant Klebsiella pneumoniae isolates. MCA supplemented with sodium azide (200 mg/l) and ertapenem (0.5 mg/l) was used to select for transconjugants. Conjugation frequency per recipient was expressed by dividing the number of transconjugants by the initial number of recipients.
Plasmids were typed according to their incompatibility group using the PCR based replicon typing scheme PBRT 2.0 Kit (Diatheva), as described previously (Carattoli, 2009).
Molecular Typing
PFGE was performed on all ESBL-, carbapenemase-, and MCR- positive E. coli and K. pneumoniae isolates. All the obtained pulsotypes were than compared with the hyper-epidemic clones previously collected in the same area (e.g., 20LOM E. coli), Genomic DNA of isolates was analyzed after digestion with XbaI restriction enzyme. Fragments were separated on a CHEF DRII system (Bio-Rad, Hercules, CA) at 14°C at 6 V/cm for 20 h with an initial pulse time of 0.5 s and a final pulse time of 30 s. Lambda 48.5 kb concatamers (New England BioLabs, Beverly, MA, USA) were used as molecular size markers.
Dendrograms of strain relatedness were created with Fingerprinting II version 3.0 software (Bio-Rad) using UPGMA. The Dice correlation coefficient was used with a 1.0% position tolerance to analyze the similarities of the banding patterns. The restriction patterns of the genomic DNA from the isolates were analyzed and interpreted according to the criteria described previously (Tenover et al., 1995).
Bacterial strains were cultivated in MCA; one single colony per strain was used for DNA extraction using NucleoSpin Tissue (Macherey-Nagel) kit. MLSTs genes were sequenced using the automatic sequencer ABI PRISM 3100 genetic analyzer DNA Sequencer (Applied Biosystems, Foster City, CA, USA) and the BigDye Terminator v1.1 Cycle Sequencing kit (Applied Biosystems, Foster City, CA, USA).
MLST for E. coli was done according to the MLST Databases at University of Warwick (http://mlst.warwick.ac.uk/mlst/dbs/Ecoli) and for K. pneumoniae according to the Institute Pasteur MLST (http://bigsdb.pasteur.fr/klebsiella/klebsiella.html). MLST of Enterobacter cloacae was done according to the E. cloacae MLST website (https://pubmlst.org/ecloacae/) developed by Keith Jolley and sited at the University of Oxford.
Plasmid Characterization
7TE.coli plasmid sequencing was done by Illumina MiSeq technology, using Nextera XT kits for library preparation. Reads were assembled using Spades 3.8 (http://cab.spbu.ru/software/spades/).
The detection of resistance genes and plasmid replication sites have been performed on J53R7TE. coli using the ResFinder and PlasmidFinder tools in the DTU database (http://www.genomicepidemiology.org/).
Contigs containing plasmid sequences were detected, analyzed, and closed using the software Bandage (https://rrwick.github.io/Bandage/). Open Reading Frames (ORFs) and the relative amino acids were predicted using Artemis (http://www.sanger.ac.uk/science/tools/artemis). The annotation was performed manually using the online blast tool on the ncbi database (https://blast.ncbi.nlm.nih.gov/Blast.cgi?CMD=Web&PAGE_TYPE=BlastHomeNew).
Genbank files were formatted and uploaded using the Sequin software (https://www.ncbi.nlm.nih.gov/Sequin/).
Statistical Analysis
Data were expressed as means (± Standard Error). The difference between means of the bacterial count was tested by the Analysis of Variance (ANOVA). A P value below 0.05 was deemed significant.
Results
Identification of Gram-Negative Bacteria Isolates
A total of 246 water samples was collected; 48% (118/246) from five streams, 43.9% (108/246) were from 11 wells, and 8.1% (20/246) from four WWTPs.
Two hundred and sixty-four non-duplicate microorganisms, one-half (50%; n = 132) belonging to enterobacterial species, were overall obtained on MCA+CTX/ChromArt ESBL. The remaining isolates included Pseudomonas spp. (19%; n = 49), Acinetobacter spp. (11%; n = 29), Vibrio fluvialis (7%; n = 18), Aeromonas hydrophila (9%; n = 25), and other Gram-negative bacteria (4%; n = 11).
Bacterial Counts and Percentage of 3GCs Resistance
The number in average of CFU ml−1 of bacteria identified after membrane filtration on PCA ranged from zero to 10,000 (data not shown). High bacterial densities were always observed in the case of WWTPs (5,100–10,000 CFU ml−1), while a month to month fluctuating trend have been observed in streams samples, with number in average bacterial counts ranging between 1,023 and 8,272 CFU ml−1(Table 2).
Concerning the wells results, the bacterial detection frequencies varied between 12.9 and 9,005 CFU ml−1; highest densities being found for the three wells named 6P, 9P, and 10P (Table 2).
The percentages of 3GCs-resistant Gram-negatives, obtained comparing the bacterial growth on MCA and MCA+CTX, ranged between 0 and 100%, depending on the single visit and sampling site (Table 2).
Statistical Analysis
Seasonal changes of total bacterial count were statistically evaluated. No significant seasonal effects on total vital count were found analyzing samples from wells and WWTPs.
The lowest bacterial count (as average value) was registered from the streams in October (126.4 ± 34.1). A significant increase (p < 0.05) in the streams mean bacterial count was observed, starting from the beginning of winter toward the end of spring. If the highest value of the mean bacterial count (10,000 ± 202.5) has been recorded in April, a significant decrease (p < 0.05) was observed from the beginning of summer to the end of autumn.
The highest percentage of CTX-resistant Gram-negative microorganisms (65.65 ± 13.1) was recorded in May. A significant decrease (p < 0.05) in the same percentage was observed since the beginning of summer toward the end of autumn. The lowest value was found in September (10.2 ± 7.8). The mean values of both total bacterial count and percentage of CTX-resistant Gram-negatives are shown in Table 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/272364/fmicb-08-02232-HTML-r1/image_m/fmicb-08-02232-t003.jpg)
Table 3. Mean values of total bacterial count and percentage of 3GCs-resistance (±SE) from stream samples during different seasons.
Phenotypic Testing and Antibiotic Resistance Patterns
A total of 132 enterobacteria grew on MCA+CTX: 82 (62.1%) were from streams, 41 (31.1%) from treatment plants, and 9 (6.8%) from wells samples. Only 30/132 (22.7%) Enterobacteriaceae showed a synergic effect by DD phenotypic test using 3GCs and piperacillin-tazobactam; 15/30 were E. coli, 7/30 Yersinia enterocolitica, 5/30 K. pneumoniae, 2/30 Klebsiella oxytoca, and 1/30 E. cloacae. Two (1.5%) out of 132 Enterobacteriaceae, were phenotypically confirmed as carbapenemases producers by KPC/MBL kit, while 1/132 (0.8%) was identified as AmpC-producer.
All the collected Y. enterocolitica strains showed clinical resistance to 3GCs; the 12.5% were also chloramphenicol and fosfomycin resistant. All the E. coli and Klebsiella spp. isolates showed resistance to 3GCs; the 13.3% (n = 2) of E. coli and 22.2% (n = 2) of Klebsiella spp. resulted, in addition, resistant to the three main carbapenems (ertapenem, meropenem, imipenem). Notably, 60% (n = 9) ESBLs-producing E. coli isolates showed fluoroquinolones resistance (Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/272364/fmicb-08-02232-HTML-r1/image_m/fmicb-08-02232-t004.jpg)
Table 4. Molecular characteristics of ESBLs/AmpC/carbapenemases-producing isolates recovered from streams, wells, and WWTPs.
Only 3/15 (20%) of the E. coli (7TE. coli, 1TPE. coli Stradella, 3TPE. coli Broni) and 1/7 (14.2%) of the K. pneumoniae (9P) out of the 33 ESBLs (n = 30)/AmpC (n = 1)/carbapenemases (n = 2) phenotypically suspected enterobacteria, resulted CO-resistant (MIC > 4 mg/l) by MicroScan-A4 System.
Clinical CO resistance was confirmed only for the strain named 7TE. coli (MIC = 8 mg/l) by MIC broth microdilution method.
KPC/MBL confirm kit test showed a positive result for the 9P and 5TPK. pneumoniae strains, obtained on June 2015 and July 2015 from a well and a WWTP in the Varzi area, respectively. The boronic acid synergistic effect suggested the production of a KPC-type enzyme.
Detection and Characterization of Resistance Genes
A blaCTX−M-type gene, was identified in a total of 21/33 (63.6%) isolates (n = 14 E. coli; n = 5 K. pneumoniae; n = 1 K. oxytoca; n = 1 E. cloacae) from all water compartments (Table 4).
In four cases (3TE. coli; 3TPK. pneumoniae; 7T E. coli; 6TP E. coli Varzi) the blaCTX−M-type gene co-existed with other β-lactamase determinants, as shown in Table 4. The most identified blaCTX−M variant was blaCTX−M−1 (8/21; 38%), followed by the blaCTX−M−15 (5/21; 23.8%), the blaCTX−M−14 and blaCTX−M−28 (3/21, each; 14.3%, each), and blaCTX−M−138 (2/21; 9.5%).
A blaSHV−12 gene was present in three isolates, named 11TE. coli, 7TE. coli, and 6TK. oxytoca.
Only 1/8 Chromogenic ESBL-agar selected Y. enterocolitica strains (the cefoxitin resistant 11T) harbored a blaDHA-type gene, while the remaining seven resulted negative for any bla resistance gene by Check-MDR CT103 XL array.
PCR and sequencing confirmed the presence of the blaKPC−2 and blaKPC−3 determinants in the 5TP and 9PK. pneumoniae strains, respectively. A co-occurrence of mcr-1.2, blaCTX−M−1, and blaSHV−12 genes was detected by in-silico plasmid characterization in the 7TE. coli collected on November 2015 from Coppa stream (Table 4).
Conjugation Experiments Results
Conjugation experiments were performed on all the 25 ESBLs/KPC positive strains: 15 E. coli, seven K. pneumoniae, two K. oxytoca, one E. cloacae.
Lateral transfer of resistance genes was observed in 15/25 (60%) isolates: 11 E. coli, three K. pneumoniae, and one K. oxytoca, as shown in Table 5. The resistance profiles of donors and transconjugants confirmed lateral transfer of 3GCs, carbapenems (ertapenem), and colistin resistance. PCR analysis confirmed the presence of the resistance genes in all the transconjugants (Table 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/272364/fmicb-08-02232-HTML-r1/image_m/fmicb-08-02232-t005.jpg)
Table 5. Phenotipical and molecular characteristics of 3GCs-resistant transconjugants; comparison with E. coli and Klebsiella spp. donors resistance and β-lactamases.
The transfer of CTX resistance was observed at a frequency of ~10−3 transconjugants per recipient. Compared to the E. coli J53 and J62 strains used as recipients, the transconjugants exhibited a decreased susceptibility to 3GCs (data not shown).
While in the case of KPC-3-producing 9PK. pneumoniae the conjugation experiment failed, the transfer of the CO resistance trait from MCR-1.2-producing 7TE. coli to J53 (met, pro-, rifr) E. coli K12 was possible, at a frequency of ~10−2 transconjugants per recipient. Susceptibility testing by autoSCAN-4 System revealed a MIC ≥ 4 mg/l for the mcr-1-harboring J53R7T E. coli.
PCR and sequencing analysis on J53R7T E. coli yielded positive results for the presence of the mcr-1.2 gene. Inc group plasmid analysis performed on J53R7TE. coli using PBRT 2.0 kit, showed that mcr-1.2 resistance determinant was located in plasmid belonging to the IncX4 Group. In addition to the IncX4, also IncX3 incompatibility group was observed in the donor 7T E. coli.
Molecular Typing
PFGE analysis carried out on all 15 ESBLs-producing E. coli showed clonal diversity (Figure 2). The 1TE. coli strain resulted identical to the 20LOM E. coli ST 131 (100% similarity between pulsotypes), previously collected in the same area (Figure 2B). Five out of seven ESBLs- and/or carbapenemases- producing K. pneumoniae were unique by PFGE (Figure 3). Multiple MLST lineages were identified, coherently with the previously obtained PFGE clonal relationships. Among these, were detected the hyperepidemic clones ST131 E. coli and ST258 K. pneumoniae (Table 4). The ST1601 was shared by the two PFGE clonally related strains (but showing different resistance profiles and blaCTX−M variants) collected in May 2015 from Varzi WWTP.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/272364/fmicb-08-02232-HTML-r1/image_m/fmicb-08-02232-g002.jpg)
Figure 2. (A) Cluster analysis of the ESBL-producing E. coli isolates. (B) Comparison between the profiles of 1T E.coli and 20LOM E. coli. Date of isolation (month, year), resistance determinants and ST are also included. The scale bar at the top (left) indicates similarity coefficient (%).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/272364/fmicb-08-02232-HTML-r1/image_m/fmicb-08-02232-g003.jpg)
Figure 3. Cluster analysis of the seven ESBLs- or carbapenemases- producing K. pneumoniae. Date of isolation (month, year) and resistance determinants are also included. The scale bar at the top (left) indicates similarity coefficient (%).
Plasmid Characterization
In-silico analysis of the NGS data of the transconjugant J53R7TE. coli, revealed the presence of the CO resistance gene mcr-1.2 on an IncX4 plasmid. Using Bandage, the plasmid was extracted and closed resulting in a 33,303 bp circular plasmid named pIBMC-MCR1.2 (Figure 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/272364/fmicb-08-02232-HTML-r1/image_m/fmicb-08-02232-g004.jpg)
Figure 4. Linear map of the pIBMC_MCR1.2 plasmid, showing the genetic structures surrounding mcr-1.2 gene.
pIBMC-MCR1.2 showed 99% similarity with the recently reported pMCR1.2-IT plasmid detected in a KPC-3-producing ST512 K. pneumoniae clinical isolate collected in Italy (Di Pilato et al., 2016). Other similar IncX4 plasmids carrying mcr-1 gene were reported in China from a clinical K. pneumoniae isolate, in Estonia from a pig sludge E. coli isolate, and in South Africa from an E. coli clinical isolate (Li et al., 2016; Poirel et al., 2016). In all these plasmids, the mcr neighboring regions were highly similar, with the downstream presence of pap2 gene (encoding for a putative PAP family transmembrane protein), and the absence of the mobilization region ISApl1.
The plasmid pIBMC-MCR1.2 has an average G/C-content of 42%; it comprises 26 ORFs, of which 21 encoding proteins with known functions and five hypothetical proteins. Any additional resistance genes are carried on this plasmid. The plasmid contains a region responsible for mobilization (14.1 kb) and another for replication (1 kb) and maintenance (2.7 kb) (Figure 4).
Discussion
The occurrence of antibiotic-resistant bacteria from European aquatic sites is increasingly described (Zarfel et al., 2017).
The dissemination in such settings of ESBLs-producing Enterobacteriaceae could be particularly worrisome, due to the ability of nosocomial pathogens to transfer ARGs among different hosts and environments (Machado et al., 2009).
This study describes the occurrence of ESBLs, acquired cephalosporinases and carbapenemases among Enterobacteriaceae from surface, ground water and WWTP aquatic compartments of the Oltrepò Pavese Plain, Italy. Moreover, this is the first report on the presence of a CTX-M-1, SHV-12 ST10 E. coli strain harboring a mcr-1.2 gene from an Italian stream. CTX-M-type producing ST10 E. coli was already reported from farm animals, healthy, and hospitalized humans (Hansen et al., 2014).
The Oltrepò Pavese area is densely populated (355 inhabitants/km2 in average), supports many agricultural (as vineyards, orchards), intensive livestock, and industrial activities. Although the study was undertaken in a restricted area, the microbiological results here presented can be regarded as representative, since this hydrogeological site is typical of many other Po plain zones.
During one-year period, a total of 33 CTX-M-, SHV-, DHA-, KPC-type producing Enterobacteriacae were identified from 11 wells, five streams, and four WWTPs. A β-lactamase production mechanism was not detected in a large amount of the selected 3GCs resistant enterobacterial species (n = 99/132; 75%).
High levels of bacterial contamination and CTX-resistance rates were constantly observed in WWTPs, while seasonal changes—with highest values in spring—were recorded from stream samples. Moreover, fluctuations occurred depending on the sampling sites location, up- or downstream the rivers. The above trend, that well highlights how the local anthropogenic activities directly affect the water quality, was observed also in a recent research from Guadeloupe (Guyomard-Rabenirina et al., 2017).
Surprisingly, high bacterial counts were detected sampling 2/11 wells located nearby the heavily contaminated Varzi WWTP; these data confirm the huge pollution risks existing in this short term (about 60 days) ground water recharge area.
The CTX-M-type ESBLs resulted the most spread 3GCs hydrolyzing enzymes among the collected isolates. The overall β-lactamases detected were mainly of CTX-M-type, followed by SHV-type, TEM-1, and KPC-type.
The high (60%) fluoroquinolones non-susceptibility level detected in 3GCs resistant E. coli represents a particular concern, as this antibiotic class is frequently employed in treatment of urinary tract infections.
Fifteen out of 33 (45.4%) β-lactamases-producing Enterobacteriacae were collected from WWTPs.
Although in wells 3GCs-resistance rates resulted very low, a ST258 KPC-3-producing 9PK. pneumoniae strain has been found in Varzi area. To our knowledge, only one description of an environmental strain ST512 KPC-3-producing K. pneumoniae was previously reported in Central Italy from a WWTP (Perilli et al., 2013). In 2011, the European Antimicrobial Resistance Surveillance Network (EARS-Net) described an ongoing worrisome increasing trend of carbapenemase-producing K. pneumoniae from invasive samples in Italian Intensive Care Units. Moreover, recent National data show that KPC-producing K. pneumoniae Clonal Complex 258 is also frequently found in patients of geriatric or medicine wards (Cristina et al., 2016).
Seven out of eight Y. enterocolitica, although initially suspected as ESBLs-positive using phenotypic tests, resulted bla genes negative by microarray. The detection of Y. enterocolitica showing CTX and cefepime clinical resistance is not common in our country. Thereby, as the intrinsic AmpC production is not enough to explain this phenotype, other mechanisms as porin loss and/or over-expression of efflux systems could be responsible in the resistance detected in our isolates (Frazão et al., 2017).
The plasmid-mediated mcr-1 gene has recently been reported from food animals, vegetables, asymptomatic carriers, and hospitalized patients in China (Liu et al., 2016). Since then, it has been across five continents (Wang et al., 2017). Regarding the environmental compartment, mcr-1 was reported in Malaysa (Zurfuh et al., 2016), Switzerland river (Li et al., 2016), and Spain sewage water (Ovejero et al., 2017). In the present study, a mcr-1.2 variant was detected in the 7TE. coli on a conjugative IncX4 plasmid named pIBMC-MCR1.2. Other similar IncX4 plasmids carrying mcr-1 determinat were reported in China, Estonia, South Africa, and very recently Italy, in different Enterobacteriaceae from both clinical and animal origin (Di Pilato et al., 2016; Li et al., 2016; Poirel et al., 2016). In the Italian strain, in particular, the mcr1.2 gene is carried on a plasmid showing 99% similarity with the here reported pIBMC-MCR1.2.
The spreading potential of mcr-1-harboring plasmids in MDR microorganisms poses significant challenges for clinical treatment and infection control strategies, as when β-lactams, aminoglycosides or quinolones are ineffective, colistin serve as the final alternative. The advent of transmissible colistin resistance indicates that the evolution of Enterobacteriaceae from extensively to pan-drug resistant is inevitable (Biswas et al., 2012).
The occurrence of ESBLs-, DHA-, KPCs-, MCR-1.2- environmental Enterobacteriaceae highlights the importance to improve surveillance and remediation actions on surface and ground waters of Oltrepò Pavese area. In this regard, and out of the aims of this report, we observed clonal relatedness between the 1TE. coli and the 20LOME.coli strain isolated in October 2013 from a Long Term and Rehabilitation Facility located near the Versa stream (Figure 2B). Appropriate measures urgently need to be enforced in order to reduce the anthropogenic burden of antibiotic resistance. Improvement of water status is of major concern: new strategies for the treatment of wastewaters, i.e., the use of sand filters or more-stringent chlorine disinfection, need to be taken into consideration to prevent resistant bacteria from being released into the aquatic environment.
The current results support the hypothesis that environmental water might represent an important hidden resistance reservoir.
Nucleotide Sequence Accession Number
The complete sequence of the pIBMC_MCR1.2 plasmid was deposited into GenBank under the accession number MF449287.
Author Contributions
MC and EN played an important role in interpreting the results and in writing the manuscript. MS and GP helped to acquired data. FZ, SP, FN, VM, AP, IB, and MD carried out experimental work. LP and RM supervised the experiments and revised the final manuscript, which was approved by all authors.
Funding
The study was supported by Banca Del Monte Foundation (Grant F12I14000180007). The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Conflict of Interest Statement
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.
Acknowledgments
We would like to thank Dr. R. Brerra for technical support and Mr. G. Chiesa (ACAOP Spa-Stradella), and Mr. C. Passarella (ASM-Voghera) for the possibility to reach the wells sampling sites.
References
Baquero, F., Martínez, J. L., and Cantón, R. (2008). Antibiotics and antibiotic resistance in water environments. Curr. Opin. Biotechnol. 19, 260–265. doi: 10.1016/j.copbio.2008.05.006
Berglund, B. (2015). Environmental dissemination of antibiotic resistance genes and correlation to anthropogenic contamination with antibiotics. Infect. Ecol. Epidemiol. 5:28564. doi: 10.3402/iee.v5.28564
Biswas, S., Brunel, J. M., Dubus, J. C., Reynaud-Gaubert, M., and Rolain, J. M. (2012). Colistin: an update on the antibiotic of the 21st century. Expert Rev. Anti Infect. Ther. 8, 917–934. doi: 10.1586/eri.12.78
Bonomo, R. A. (2017). β-Lactamases: a focus on current challenges. Cold Spring Harb. Perspect. Med. 7:a025239. doi: 10.1101/cshperspect.a025239
Bouki, C., Venieri, D., and Diamadopoulos, E. (2013). Detection and fate of antibiotic resistant bacteria in wastewater treatment plants: a review. Ecotoxicol. Environ. Saf. 91, 1–9. doi: 10.1016/j.ecoenv.2013.01.016
Carattoli, A. (2009). Resistance plasmid families in Enterobacteriaceae. Antimicrob. Agents Chemother. 53, 2227–2238. doi: 10.1128/AAC.01707-08
Colombo, S., Arioli, S., Neri, E., Della Scala, G., Gargari, G., and Mora, D. (2017). Viromes as genetic reservoir for the microbial communities in aquatic environments: a focus on antimicrobial-resistance genes. Front. Microbiol. 8:1095. doi: 10.3389/fmicb.2017.01095
Cristina, M. L., Sartini, M., Ottria, G., Schinca, E., Cenderello, N., Crisalli, M. P., et al. (2016). Epidemiology and biomolecular characterization of carbapenem-resistant Klebsiella pneumoniae in an Italian hospital. J. Prev. Med. Hyg. 3, 149–156. doi: 10.15167/2421-4248/jpmh2016.57.3.639
Davies, J., and Davies, D. (2010). Origins and evolution of antibiotic resistance. Microbiol. Mol. Biol. Rev. 74, 417–433. doi: 10.1128/MMBR.00016-10
Di Pilato, V., Arena, F., Tascini, C., Cannatelli, A., Henrici De Angelis, L., Fortunato, S., et al. (2016). mcr-1.2, a new mcr variant carried on a transferable plasmid from a colistin-resistant KPC carbapenemase-producing Klebsiella pneumoniae strain of sequence type 512. Antimicrob. Agents Chemother. 60, 5612–5615. doi: 10.1128/AAC.01075-16
Eckert, C., Gautier, V., Saladin-Allard, M., Hidri, N., Verdet, C., Ould-Hocine, Z., et al. (2004). Dissemination of CTX-M-type beta-lactamases among clinical isolates of Enterobacteriaceae in Paris, France. Antimicrob. Agents Chemother. 48, 1249–1255. doi: 10.1128/AAC.48.4.1249-1255.2004
Frazão, M. R., Andrade, L. N., Darini, A. L. C., and Falcão, J. P. (2017). Antimicrobial resistance and plasmid replicons in Yersinia enterocolitica strains isolated in Brazil in 30 years. Braz. J. Infect. Dis. 21, 477–480. doi: 10.1016/j.bjid.2017.04.006
Guyomard-Rabenirina, S., Dartron, C., Falord, M., Sadikalay, S., Ducat, C., Richard, V., et al. (2017). Resistance to antimicrobial drugs in different surface waters and wastewaters of Guadeloupe. PLoS ONE 12:e0173155. doi: 10.1371/journal.pone.0173155
Hansen, F., Olsen, S. S., Heltberg, O., Justesen, U. S., Fuglsang-Damgaard, D., Knudsen, J. D., et al. (2014). Characterization of third-generation cephalosporin-resistant Escherichia coli from bloodstream infections in Denmark. Microb. Drug Resist. 20, 316–324. doi: 10.1089/mdr.2013.0157
Jarlier, V., Nicolas, M. H., Fournier, G., and Philippon, A. (1988). Extended broad-spectrum β-lactamases conferring transferable resistance to newer β-lactam agents in Enterobacteriaceae: hospital prevalence and susceptibility patterns. Rev. Infect. Dis. 10, 867–878 doi: 10.1093/clinids/10.4.867
Li, A., Yang, Y., Miao, M., Chavda, K. D., Mediavilla, J. R., Xie, X., et al. (2016). Complete sequences of mcr-1-harboring plasmids from extended-spectrum-β-lactamase- and carbapenemase-producing Enterobacteriaceae. Antimicrob. Agents. Chemother. 7, 4351–4354. doi: 10.1128/AAC.00550-16
Liu, Y. Y., Wang, Y., Walsh, T. R., Yi, L. X., Zhang, R., Spencer, J., et al. (2016). Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study. Lancet Infect. Dis. 16, 161–168. doi: 10.1016/S1473-3099(15)00424-7
Machado, E., Coque, T. M., Cantón, R., Sousa, J. C., Silva, D., Ramos, M., et al. (2009). Leakage into Portuguese aquatic environments of extended-spectrum-β-lactamase-producing Enterobacteriaceae. J. Antimicrob. Chemother. 63, 616–618. doi: 10.1093/jac/dkn510
Marti, E., Huerta, B., Rodríguez-Mozaz, S., Barceló, D., Jofre, J., and Balcázar, J. L. (2014). Characterization of ciprofloxacin-resistant isolates from a wastewater treatment plant and its receiving river. Water. Res. 61, 67–76. doi: 10.1016/j.watres.2014.05.006
Olaitan, A. O., Morand, S., and Rolain, J. M. (2016). Emergence of colistin-resistant bacteria in humans without colistin usage: a new worry and cause for vigilance. Int. J. Antimicrob. Agents 47, 1–3. doi: 10.1016/j.ijantimicag.2015.11.009
Ovejero, C. M., Delgado-Blas, J. F., and Calero-Caceres, W. (2017). Spread of mcr-1-carrying Enterobacteriaceae in sewage water from Spain. J. Antimicrob. Chemother. 72, 1050–1053. doi: 10.1093/jac/dkw533
Pagani, L., Dell'Amico, E., Migliavacca, R., D'Andrea, M. M., Giacobone, E., Amicosante, G., et al. (2003). Multiple CTX-M-type extended-spectrum β-lactamases in nosocomial isolates of enterobacteriaceae from a hospital in Northern Italy. J. Clin. Microbiol. 41, 4264–4269. doi: 10.1128/JCM.41.9.4264-4269.2003
Perilli, M., Bottoni, C., Pontieri, E., Segatore, B., Celenza, G., Setacci, D., et al. (2013). Emergence of bla(KPC-3)-Tn4401a in Klebsiella pneumoniae ST512 in the municipal wastewater treatment plant and in the university hospital of a town in central Italy. J. Glob. Antimicrob. Resist. 1, 217–220. doi: 10.1016/j.jgar.2013.07.002
Perilli, M., Dell'Amico, E., Segatore, B., de Massis, M. R., Bianchi, C., Luzzaro, F., et al. (2002). Molecular characterization of extended-spectrum beta-lactamases produced by nosocomial isolates of Enterobacteriaceae from an Italian nationwide survey. J. Clin. Microbiol. 40, 611–614. doi: 10.1128/JCM.40.2.611-614.2002
Pilla, G., Sacchi, E., and Ciancetti, G. (2007). “Hydrochemical and isotopic groundwater investigation in the Oltrepo region (Po valley, Northern Italy),” in Proceedings of Symposium: “Advances in Isotope Hydrology and its Role in Sustainable Water Resources Management (HIS-2007),” (Vienna: IAEA, Proceedings Series; IAEA–CN−151/41), 49–58.
Poirel, L., Kieffer, N., Brink, A., Coetze, J., Jayol, A., and Nordmann, P. (2016). Genetic features of MCR-1-producing colistin-resistant Escherichia coli isolates, South Africa. Antimicrob. Agents Chemother. 60, 4394–4397. doi: 10.1128/AAC.00444-16
Poppe, C., Martin, L., Muckle, A., Archambault, M., McEwen, S., and Weir, E. (2006). Characterization of antimicrobial resistance of Salmonella Newport isolated from animals, the environment, and animal food products in Canada. Can. J. Vet. Res. 70, 105–114.
Rasheed, J. K., Jay, C., Metchock, B., Berkowitz, F., Weigel, L., Crellin, J., et al. (1997). Evolution of extended-spectrum beta-lactam resistance (SHV-8) in a strain of Escherichia coli during multiple episodes of bacteremia. Antimicrob. Agents Chemother. 41, 647–653.
Schwartz, T., Kohnen, W., Jansen, B., and Obst, U. (2003). Detection of antibiotic-resistant bacteria and their resistance genes in wastewater, surface water, and drinking water biofilms. FEMS Microbiol. Ecol. 43, 325–435. doi: 10.1111/j.1574-6941.2003.tb01073.x
Tenover, F. C., Arbeit, R. D., Goering, R. V., Mickelsen, P. A., Murray, B. E., Persing, D. H., et al. (1995). Interpreting chromosomal DNA restriction patterns produced by pulsed-field gel electrophoresis: criteria for bacterial strain typing. J. Clin. Microbiol. 33, 2233–2239.
Wang, Y., Tian, G. B., Zhang, R., Shen, Y., Tyrrell, J. M., Huang, X., et al. (2017). Prevalence, risk factors, outcomes, and molecular epidemiology of mcr-1-positive Enterobacteriaceae in patients and healthy adults from China: an epidemiological and clinical study. Lancet Infect. Dis. 17, 390–399. doi: 10.1016/S1473-3099(16)30527-8
Watkinson, A. J., Micalizzi, G. R., Bates, J. R., and Costanzo, S. D. (2007). Novel method for rapid assessment of antibiotic resistance in Escherichia coli isolates from environmental waters by use of a modified chromogenic agar. Appl. Environ. Microbiol. 73, 2224–2229. doi: 10.1128/AEM.02099-06
Wellington, E. M., Boxall, A. B., Cross, P., Feil, E. J., Gaze, W. H., Hawkey, P. M., et al. (2013). The role of the natural environment in the emergence of antibiotic resistance in gram-negative bacteria. Lancet Infect. Dis. 13, 155–165. doi: 10.1016/S1473-3099(12)70317-1
Yigit, H., Queenan, A. M., Anderson, G. J., Domenech-Sanchez, A., Biddle, J. W., Steward, C. D., et al. (2001). Novel carbapenem-hydrolyzing b-lactamase, KPC-1, from a carbapenem-resistant strain of Klebsiella pneumoniae. Antimicrob. Agents Chemother. 45, 1151–1161. doi: 10.1128/AAC.45.4.1151-1161.2001
Zanotto, C., Bissa, M., Illiano, E., Mezzanotte, V., Marazzi, F., Turolla, A., et al. (2016). Identification of antibiotic-resistant Escherichia coli isolated from a municipal wastewater treatment plant. Chemosphere 164, 627–633. doi: 10.1016/j.chemosphere.2016.08.040
Zarfel, G., Lipp, M., Gürtl, E., Folli, B., Baumert, R., and Kittinger, C. (2017). Troubled water under the bridge: screening of river mur water reveals dominance of CTX-M harboring Escherichia coli and for the first time an environmental VIM-1 producer in Austria. Sci. Total Environ. 593–594, 399–405. doi: 10.1016/j.scitotenv.2017.03.138
Zurfuh, K., Poirel, L., Nordmann, P., Nüesch-Inderbinen, M., Hächler, H., and Stephan, R. (2016). Occurrence of the plasmid-borne mcr-1 colistin resistance gene in extended-spectrum-β-lactamase-producing Enterobacteriaceae in river water and imported vegetable samples in Switzerland. Antimicrob. Agents. Chemother. 60, 2594–2595. doi: 10.1128/AAC.00066-16
Keywords: water ecosystems, Gram-negative bacteria, carbapenemases, colistin resistance, molecular characterization
Citation: Caltagirone M, Nucleo E, Spalla M, Zara F, Novazzi F, Marchetti VM, Piazza A, Bitar I, De Cicco M, Paolucci S, Pilla G, Migliavacca R and Pagani L (2017) Occurrence of Extended Spectrum β-Lactamases, KPC-Type, and MCR-1.2-Producing Enterobacteriaceae from Wells, River Water, and Wastewater Treatment Plants in Oltrepò Pavese Area, Northern Italy. Front. Microbiol. 8:2232. doi: 10.3389/fmicb.2017.02232
Received: 06 April 2017; Accepted: 31 October 2017;
Published: 10 November 2017.
Edited by:
Manuela Caniça, Istituto Nacional de Saúde, PortugalReviewed by:
Alberto Quesada, Universidad de Extremadura, SpainMonika Dolejska, University of Veterinary and Pharmaceutical Sciences Brno, Czechia
Copyright © 2017 Caltagirone, Nucleo, Spalla, Zara, Novazzi, Marchetti, Piazza, Bitar, De Cicco, Paolucci, Pilla, Migliavacca and Pagani. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Roberta Migliavacca, r.miglia@unipv.it