Sporadic clone Escherichia coli ST615 as a vector and reservoir for dissemination of crucial antimicrobial resistance genes

There is scarce information concerning the role of sporadic clones in the dissemination of antimicrobial resistance genes (ARGs) within the nosocomial niche. We confirmed that the clinical Escherichia coli M19736 ST615 strain, one of the first isolates of Latin America that harbors a plasmid with an mcr-1 gene, could receive crucial ARG by transformation and conjugation using as donors critical plasmids that harbor bla_CTX-M-15, bla_KPC-2, bla_NDM-5, bla_NDM-1, or aadB genes. Escherichia coli M19736 acquired bla_CTX-M-15, bla_KPC-2, bla_NDM-5, bla_NDM-1, and aadB genes, being only blaNDM-1 maintained at 100% on the 10th day of subculture. In addition, when the evolved MDR-E. coli M19736 acquired sequentially bla_CTX-M-15 and bla_NDM-1 genes, the maintenance pattern of the plasmids changed. In addition, when the evolved XDR-E. coli M19736 acquired in an ulterior step the paadB plasmid, a different pattern of the plasmid’s maintenance was found. Interestingly, the evolved E. coli M19736 strains disseminated simultaneously the acquired conjugative plasmids in different combinations though selection was ceftazidime in all cases. Finally, we isolated and characterized the extracellular vesicles (EVs) from the native and evolved XDR-E. coli M19736 strains. Interestingly, EVs from the evolved XDR-E. coli M19736 harbored bla_CTX-M-15 though the pDCAG1-CTX-M-15 was previously lost as shown by WGS and experiments, suggesting that EV could be a relevant reservoir of ARG for susceptible bacteria. These results evidenced the genetic plasticity of a sporadic clone of E. coli such as ST615 that could play a relevant transitional link in the clinical dynamics and evolution to multidrug/extensively/pandrug-resistant phenotypes of superbugs within the nosocomial niche by acting simultaneously as a vector and reservoir of multiple ARGs which later could be disseminated.

1 Introduction

Escherichia coli is common among the aerobic bacteria in the gastrointestinal tract microbiota of both humans and mammals (Denamur et al., 2021). Simultaneously, some lineages have developed into a pathogen well adapted to their host causing different diseases (Geurtsen et al., 2022), including adaptation to the nosocomial niche as high-risk clones or “superbugs” that rapidly evolve to extreme drug resistance (XDR). The lineage represented by sequence type (ST) 131 is the predominant isolate of hospital infections worldwide among E. coli strains that behave like an epidemic clone, also known as pandemic clones (Pitout and DeVinney, 2017; Soncini et al., 2022; Pitout and Chen, 2023). On the other hand, little is known about the role of sporadic clones of this species in the adaptation to multidrug resistance among clinical isolates. Since a few strains unrelated to outbreaks of E. coli ST615 have been reported from Tunisia (Maamar, 2016), Poland (Jamborova et al., 2015), and Spain (Ojer-Usoz et al., 2017) and E. coli M19736 from Argentina (Tijet et al., 2017), this ST may be considered as a sporadic clone. According to the World Health Organization, multidrug resistance (MDR) in Gram-negative bacilli (Gnb) has become a challenge due to its high global incidence and prevalence. In 2019, it was estimated that 4.95 million deaths were due to infections associated with antimicrobial resistance (AMR) (Global antimicrobial resistance and use surveillance system (GLASS) report: 2022, 2022). These pathogens represent a particular threat in nosocomial infections, and among those of greatest clinical interest, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacteriaceae producers of carbapenemases have been identified as a critical priority (Rello et al., 2019). MDR Gnb usually harbor multiple mobile genetic elements such as gene cassettes, transposons, and plasmids that confer their MDR phenotypes. These mobile genetic elements can be transferred by conjugation, transformation, and transduction and by the most recently discovered mechanism known as vesiduction through the extracellular vesicles (EVs).

The rapid increase of carbapenem-resistant Gnb due to the expression of enzymes such as KPC-2 (Klebsiella pneumoniae carbapenemase-2) and NDM-1 (New Delhi metallo-β-lactamase-1) is a global public health concern. Consequently, interest in another family of antibiotics, polymyxins, has recently increased as a last resort used in medical clinics despite its high nephrotoxicity (Rodríguez et al., 2018; Binsker et al., 2022). Furthermore, the intensive use of polymyxins in veterinary medicine not only for the treatment of infections but also as a growth promoter has led to an increase of the isolation of Gnb strains resistant to this antibiotic in clinical settings (Binsker et al., 2022). Resistance to polymyxins includes chromosome-encoded resistance traits, as well as the mobile plasmid-encoded polymyxin resistance determinants such as the mcr-1 gene (Liu et al., 2016; Caniaux et al., 2017). The transferable mcr-1 gene was first detected in E. coli isolates from animals, food, and patients in China (Liu et al., 2016). Aside from this gene, 10 other mcr-like genes (from mcr-2 to mcr-10) as well as several of their variants designated as mcr-1.2, mcr-1.3, mcr-1.12, etc. have been identified (Hussein et al., 2021). The bla_KPC, bla_NDM, and mcr-like genes are usually found in conjugative plasmids of diverse incompatibility groups, which enhance the challenge to combat the bacteria possessing these antibiotic-resistant determinants (Brandt et al., 2019; Cejas et al., 2022; Knecht et al., 2022; Li C. et al., 2022; Molina et al., 2022; Sanz et al., 2022; Boutzoukas et al., 2023; Riccobono et al., 2023). Despite the key role that conjugative plasmids have in nosocomial isolates, understanding how they can persist in bacterial populations in the absence of positive selection is challenging for pandemic and sporadic clones.

The nanosize EV entities secreted by Gnb have been recognized for their cardinal importance in intercellular communication among cells and to be responsible for the modulation of various biological functions (McMillan and Kuehn, 2021; Toyofuku et al., 2023). EVs are now well recognized for their role as long-distance secretion–delivery systems that eliminate the need for cell–cell contact (Toyofuku et al., 2023). EVs transport, harbor, and deliver in a concentrated, protected, and directed way biologically active proteins, lipids, nucleic acids, metabolites, and virulence factors between two bacterial cells (Tran and Boedicker, 2017; McMillan and Kuehn, 2021; Toyofuku et al., 2023). The exchange of DNA mediated by EV has been identified as an additional form of horizontal genetic transfer (HGT) (Chatterjee et al., 2017). It has been reported that they can carry DNA associated with the membrane and protect luminal genetic material against DNases and RNases (Kim et al., 2015). Of the DNA that is transferred, there are genes associated with AMR, which has huge clinical implications since the potential of propagation of a gene depends to a great extent on the competence to transmit (Chatterjee et al., 2017; Wagner et al., 2018). Earlier studies have reported the transfer by EV of a penicillin resistance gene in Neisseria gonorrhoeae (Dorward et al., 1989); the bla_OXA-24 (Rumbo et al., 2011) and the bla_NDM-1 carbapenemase genes in Acinetobacter baumannii (Chatterjee et al., 2017); the bla_CTX-M-15 and bla_TEM-1 genes in E. coli O104:H4 (Bielaszewska et al., 2020); and finally, the bla_NDM, bla_KPC, bla_SHV, bla_CTX-M-9, and aac(6′)-Ib genes in K. pneumoniae strains (Li P. et al., 2022). In the present study, we wonder about the ability of the E. coli M19736 clinical strain belonging to the sporadic clone ST615, which was one of the first isolates harboring the mcr-1 gene in Latin America (Tijet et al., 2017), to adapt to the XDR phenotype (Magiorakos et al., 2012) by conjugation and/or transformation assays. The conservation of ARG in the evolved MDR and XDR-E. coli M19736 strains showed different patterns of maintenance as well as variability in the capacity to transfer the acquired conjugative plasmids including the co-transfer of different plasmids in several combinations. In addition, we identified EVs in the native and in the evolved XDR-E. coli M19736 strains with crucial content for the competence of AMR in the nosocomial niche, including the bla_CTX-M-15 gene though the plasmid containing this gene, pDCAG1-CTX-M-15, was previously lost. Collectively, this scenario reveals the relevant role of sporadic clones as common vectors for the dissemination of acquired AMR within the framework of the HGT processes.

2 Methods

2.1 Bacterial strains and growth conditions

Escherichia coli M19736 was isolated in November 2015, in Argentina, from the blood of a patient with peritonitis secondary to colon cancer (Tijet et al., 2017). The strain was shown by MLST to be a single-locus variant of ST615 (Rapoport et al., 2016; Tijet et al., 2017). Escherichia coli M19736 has the pM19736 plasmid with the size of 63,230 bp that harbors an mcr-1 gene as described previously (Tijet et al., 2017). Escherichia coli SM5 (Sennati et al., 2012), Klebsiella pneumoniae HA7Kp (Knecht et al., 2022), Klebsiella pneumoniae HA31Kp (Álvarez et al., 2024), and Serratia marcescens SM938 (this study) (Table 1) were used as donors for conjugation assays as described below. Also, E. coli J53 was used as a laboratory model for the control of the conjugation experiments.

Table 1

Table 1 Bacterial strains used for the experiments.

The plasmid paadB from the E. coli TOP10::paadB strain was used for the transformation assays (Table 1). The aadB gene cassette from S. marcescens SCH909 (Gambino et al., 2021) was subcloned from a pCR2.1TOPO vector (Invitrogen, Carlsbad, CA) into the commercial vector pACYC184; paadB is resistant to chloramphenicol due to the background of its vector and to gentamicin due to the Pc promoter contained upstream the aadB gene cassette (Gambino et al., 2021).

All the strains were grown in Luria–Bertani (LB) broth at 37°C with shaking (150 rpm). When needed, antibiotics were used at the following concentrations: ceftazidime (8 μ/ml), meropenem (2 μg/ml), or gentamicin (25 μg/ml).

2.2 Antibiotic susceptibility testing, minimum inhibitory concentration, and phenotypic detection of β-lactamases

Receptor, donor, and transconjugant strains were tested by antibiotic susceptibility testing (AST) and minimum inhibitory concentration (MIC). AST was carried out by the agar disk diffusion method; MIC was determined by agar dilution. Both assays were conducted according to the Clinical and Laboratory Standards Institute (CLSI) guidelines (CLSI, 2023), and the results were interpreted under the same guidelines. On the other hand, phenotypic detection of β-lactamases was determined by two tests. First, the detection of carbapenemases was performed by the modified Hodge test as previously described (Pasteran et al., 2016). Secondly, the synergy produced between the extended-spectrum cephalosporins and clavulanic acid was used for the detection of extended-spectrum β-lactamase.

2.3 Polymerase chain reaction assays

Total DNA extraction was done using the boiling technique for all the experiments of conjugation and transformation. All polymerase chain reaction (PCR) reactions were done using 2 U of Taq DNA polymerase (Inbio Highway, Tandil, Argentina) in 0,5× Taq buffer (Inbio Highway) supplemented with 2,5 mM of MgCl₂, 0,2 mM dNTP mix, and 0,4 uM of each primer in a final volume of 25 µl The PCR conditions were 5 min at 95°C, 30 cycles of 30 s at 95°C, 45 s at the appropriate annealing temperature and 1 min at 72°C, followed by a final extension of 5 min at 72°C.

Detection of each gene of interest—bla_CTX-M-15, bla_KPC-2, bla_NDM-5, bla_NDM-1, aadB, and/or mcr-1 genes—was conducted with specific primers listed in Table 2. PCR products were separated on agarose gels by electrophoresis, stained with SYBR green, and visualized by UV transillumination.

Table 2

Table 2 Primers used for the detection of genes of interest.

2.4 Conjugation assays

Conjugation experiments were performed according to a method described previously (Di Noto et al., 2016). Briefly, mating assays were carried out on LB agar plates. Escherichia coli SM5 (Caz^R), HA7Kp (Mem^R), K. pneumoniae HA31Kp (Mem^R), and S. marcescens SM98 (Mem^R) were used as donor strains, while E. coli M19736 (Cmp^R) and J53 (Azi^R) were used as recipient strains. Donor and recipient strains were diluted from saturated overnight cultures into 12 ml and grown until OD_{600 nm} ~0.6 at 37°C. The cells were harvested by centrifugation, mixed together in a ratio of 1:1, and spotted onto LB plates. They were also spotted individually on LB plates as controls. After 18 h of incubation at 37°C, mating spots were washed and resuspended in saline; serial dilutions were plated onto LB agar with the specific antibiotic to select for donor, recipient, or transconjugant cells (meropenem 3 µg/ml or ceftazidime 8 µg/ml and chloramphenicol 50 µg/ml or sodium azide 150 µg/ml). Conjugation frequency was expressed as the number of transconjugant cells per donor cell in the mating mixture at the time of plating. Transconjugants obtained were checked by plating on CROMagar and LB plates with double antibiotics. Then, the different genes of interest (bla_CTX-M-15, bla_KPC-2, bla_NDM-5, and bla_NDM-1) were detected by PCR for each transconjugant. The original recipient strains [E. coli M19736 (Cmp^R) and J53 (Azi^R)] were used as negative controls, and the donor strains [E. coli SM5 (Caz^R), HA7Kp (Mem^R), K. pneumoniae HA31Kp (Mem^R), and S. marcescens SM98 (Mem^R)] were used as positive controls (Table 1).

2.5 Transformation assays

Escherichia coli M19736 was treated with 10% glycerol previously. Afterward, electroporation was performed with a Gene Pulser™ apparatus (Bio-Rad Laboratories, Denver, USA) and conducted using the following parameters: 200 Ω resistance, 25 mF capacitance, and 2 kV voltage, resulting in a time constant between 4.5 and 5.0 ms. Transformed E. coli cells were recovered in a 1-ml LB broth and incubated at 37°C with a 200-rpm shaking for 2 h before being plated on LB agar plates supplied with gentamicin (25 μg/ml) using 100 µl.

2.6 Plasmid maintenance assay

Each strain of interest was cultured in 5 ml of LB broth without antibiotic pressure and incubated at 37°C ON with shaking (200 rpm). Consecutive subcultures were made for 10 days. An aliquot was taken from each experiment on the 1st and 10th days and plated on LB agar without antibiotics. From each of the replicates, 30 colonies were selected and analyzed. DNA was then extracted from each of these colonies using the boiling method. Next, each of the 30 colonies taken per replicate was tested by PCR for the presence or absence of each of the acquired genes of interest (bla_CTX-M-15, bla_KPC-2, bla_NDM-1, bla_NDM-5, mcr-1 and aadb). The maintenance percentage of each replicate was then calculated as follows: $\%Maintenance=\left[\frac{(n_{pc}+100\%)}{n_t}\right]$, where n_pc is the number of positive colonies for the gene of interest and n_t is the total number of colonies tested (n_t = 30). Then, the average $(\overline{X})$ of three replicates performed was calculated with their respective standard deviations (SDs).

2.7 Isolation and purification of EV

EVs were isolated from the late log-phase (16 h) culture of E. coli M19736 and evolved XDR-E. coli M19736 on day 1 of subculture. In brief, cells were cultivated in 600 ml of LB broth with 10 µg/ml of ceftazidime and subinhibitory concentrations of meropenem ~14 h at 37°C. The next day, the cultures were adjusted to OD₆₀₀ ~0.7. The cells were pelleted by centrifugation (9,500 rpm, 4°C for 20 min), and the supernatant was filtered through a 0.22-μm membrane filter (Merck Millipore, Tullagreen, Carrigtwohill, Co. Cork, Ireland) to remove cells and cellular debris. The filtrate was subjected to ultracentrifugation (100,000g) for 2 h at 4°C using a P45AT(RP45T) fixed angle rotor (HIMAC CP80NX). For washing the EV, the pellet suspended in EV buffer (137 mM of NaCl and 20 mM of HEPES [pH 7.5]) was ultracentrifuged (100,000g) for 1 h at 4°C using the same rotor. The pellet was finally resuspended in 1500 µl of buffer, filtered and frozen at −80°C. The EVs were grown in 2 ml of LB broth to test for any bacterial growth.

2.8 Dynamic light scattering

The hydrodynamic diameter (Dh) and the size distribution (polydispersity index, PDI) of different EV sources were assayed by dynamic light scattering (DLS) (DLS, Zetasizer Nano-ZS, Malvern Instruments) at a scattering angle of 173°. The nano-ZS contains a 4-mW He–Ne laser operating at a wavelength of 633 nm, a digital correlator ZEN3600, and non-invasive backscatter (NIBS^®) technology. For the measurement, 200 μl of vesicles suspended in EV buffer (137 mM NaCl and 20 mM HEPES [pH 7.5]) were used. All the samples were analyzed at 25°C. Viscosities were between 0.8880 and 0.8872 cP. Results were expressed as mean ± standard deviation (SD) of three independent samples prepared in identical conditions. Data for each single specimen were the result of at least six runs.

2.9 Transmission electron microscopy

To verify the presence of intact EV, the preparations were analyzed using transmission electron microscopy (TEM). The EV suspension was fixed in nickel grids for TEM with carbon (200 mesh) (Agar Scientific Ltd., Stansted, Essex, UK), with 2% glutaraldehyde, 4% formaldehyde, and 5% sucrose in PBS; washed three times with ultrapure water; stained with 3% uranyl acetate; allowed to dry for at least 30 min; and examined under a transmission electron microscope (Zeiss EM 109T equipped with Gatan ES1000W digital camera).

2.10 Mass spectrometry analysis

EV proteins were quantified by the Micro BCA™ Protein Assay Kit (Thermo Fisher Scientific, Rockford, IL, USA). Protein digestion from lysed EV of native E. coli M19736 was performed. We used 40 μg of protein from EV based on the Micro BCA results. The proteins were reduced and alkylated with 10 mM of DTT and 20 mM of iodoacetamide and then precipitated with 15% trichloroacetic acid/acetone and processed for liquid chromatography–MS/MS (LC-MS/MS) analysis. Mass spectrometry analysis was performed at the Proteomics Core Facility (CEQUIBIEM), University of Buenos Aires/CONICET (National Research Council) by analyzing the digests by nanoLC-MS/MS in a Thermo Scientific Q Exactive Mass Spectrometer coupled to a nanoHPLC EASY-nLC 1000 (Thermo Scientific). For the LC-MS/MS analysis, approximately 2 μg of peptides were loaded onto the column and eluted for 120 min using the reverse phase column (C18, 2 µm, 100 A, 50 µm × 150 mm) EASY-Spray Column PepMap RSLC (P/N ES801) suitable for separating protein complexes with a high degree of resolution. The flow rate used for the nano column was 300 nl min⁻¹, and the solvent ranged from 7% B (5 min) to 35% (120 min). Solvent A was 0.1% formic acid in water, whereas solvent B was 0.1% formic acid in acetonitrile. The injection volume was 2 µl. The MS equipment has a high collision dissociation cell (HCD) for fragmentation and an Orbitrap analyzer (Thermo Scientific, Q-Exactive). A voltage of 3.5 kV was used for electrospray ionization (Thermo Scientific, EASY-Spray).

XCalibur 3.0.63 (Thermo Scientific) software was used for data acquisition and equipment configuration that allows peptide identification and chromatographic separation. Full-scan mass spectra were acquired in the Orbitrap analyzer. The scanned mass range was 400–1,800 m/z, at a resolution of 70,000 at 400 m/z, and the 12 most intense ions in each cycle were sequentially isolated, fragmented by HCD, and measured in the Orbitrap analyzer. Peptides with a charge of +1 or with an unassigned charge state were excluded from fragmentation for MS2.

2.10.1 Analysis of mass spectrometry data

Q Exactive raw data were processed using Proteome Discoverer software (version 2.1.1.21 Thermo Scientific) and searched against the E. coli sequence database with trypsin specificity and a maximum of 1 missed cleavage per peptide. Carbamidomethylation of cysteine residues was set as a fixed modification, and oxidation of methionine was set as a variable modification. Proteome Discoverer searches were performed with a precursor mass tolerance of 10 ppm and product ion tolerance of 0.05 Da. Protein hits were filtered for high-confidence peptide matches with a maximum protein and peptide false discovery rate of 1% calculated by employing a reverse database strategy.

2.11 EV protein analysis

The localization of proteins was mostly acquired using DAVID (https://david.ncifcrf.gov/home.jsp). The biological process of EV proteins was derived from Gene Ontology, UniProt, and KEGG (https://www.genome.jp/kegg/). Proteins associated with antibiotic response/resistance were predicted using DAVID.

2.11.1 Protein–protein interaction network analysis

Protein–protein interaction (PPI) data were downloaded using the STRING v10.5 database (Szklarczyk et al., 2017). A PPI network of EV proteins from STRING was incorporated in Cytoscape 3.9.0 (Shannon et al., 2003), and using the Cytoscape StringApp (Doncheva et al., 2019), we constructed, analyzed, and visualized the PPI network. Gene Ontology and KEEG enrichment analysis was performed for the EV proteins.

2.12 Determination of DNA in EV

Intravesicular DNA was quantified following the method of Rumbo et al. (2011) with a few modifications. Fifty micrograms of EVs were treated with DNAse RQ1-Free RNAse 1U/µl (Promega) at 37°C for 30 min to hydrolyze the free and surface-associated DNA. The reaction was stopped with a stop solution and incubation was conducted at 65°C for 10 min. DNase-treated EVs were then lysed with 0.125% Triton X-100 (Sigma-Aldrich, USA) solution for 30 min at 37°C, and DNA was purified using a QIAamp DNA Mini Kit with the protocol for crude cell lysates and other samples (Qiagen, Maryland, USA), according to the manufacturer’s instructions. The DNA was quantified using the Nano-500 Micro-spectrophotometer. The purified DNA was used for further PCR using the primers listed in Table 2.

2.13 DNA extraction, DNA sequencing, and sequence assembly

DNA of E. coli M19736 and evolved XDR-E. coli M19736 on day 1 of subculture was extracted using the mini kit QIAamp DNA (Qiagen) following the manufacturer’s protocol for Gram-negative bacteria. The concentration and purity were measured using a NanoDrop instrument (Nano-500 Micro-Spectrophotometer, Allsheng, Hangzhou, China).

DNA was sequenced on a MiSeq sequencer (Illumina pair ends). The sequencing was performed in the Genomics and Bioinformatics Unit of ANLIS Malbrán (Argentina), and the library preparation was made according to the manufacturer’s protocol. Read quality metrics were evaluated using FASTQC v0.11.9. To remove the low-quality reads and the remaining adapters from the sequencing, trimmomatic v0.39 was used. The parameters used were as follows: -threads 8 -phred33 ILLUMINACLIP:TruSeq3.fa:2:30:10 TRAILING:20 SLIDINGWINDOW:4:20 MINLEN:50. The trimmed Fastq files were evaluated using fastqc to determine their quality. Finally, we performed a short-read assembly using Unicycler v0.4.8 (Wick et al., 2017) with default options. Unicycler was executed by the command line in the GNU/Linux environment. Subsequently, the quality of the assembled files was evaluated using the QUAST v5.0.2 program, using default parameters.

2.14 Genomic and plasmid analysis

Consensus sequences of the complete assembly were imported into the RAST (Aziz et al., 2008) and PROKKA (Seemann, 2014) databases. The search for all ARGs and efflux pumps associated with AMR was performed using ResFinder (Zankari et al., 2012), RGI 6.0.2, and CARD 3.2.7 online databases with a minimum identity of 95%. In addition, chromosomal mutations associated with AMR were searched using PointFinder available in the ResFinder database. PlasmidFinder v 2.1.6 (Clausen et al., 2018; Carattoli and Hasman, 2020) was used to detect the plasmids in our samples using the default parameters. VRprofile2 v2.0 was used to predict mobilome (Wang et al., 2022). Conjugation systems were searched against the NCBI database, and OriT was detected using the online tool oriTfinder (Li et al., 2018). The search and detection of toxin–antitoxin systems was performed using the online tool TADB v2.0 using the default parameters (Guan et al., 2023).

2.15 Statistical analysis

Statistical analysis of the data obtained was performed using the GraphPad Prism 8.0.2 program (GraphPad, La Jolla, CA, USA). Variables were expressed as median (interquartile range, IQR) and compared by Kruskal–Wallis followed by Dunn’s post-hoc test. We looked for statistically significant differences between conjugation frequencies of the E. coli M19736 strain and the laboratory control strain E. coli J53. p-values<0.05 were regarded as statistically significant.

3 Results

3.1 Ability of Escherichia coli M19736 to acquire plasmids from different species

Escherichia coli M19736 was tested as a receptor to receive crucial ARG by transformation and conjugation using different relevant plasmids as donors (Table 3) including i) pDCAG1-CTX-M-15 (>112.000 bp, IncFII) from the clinical strain E. coli SM5 that harbors bla_CTX-M-15, ii) pDCCK₁-KPC (>77.218 bp, IncM1) from the clinical strain K. pneumoniae HA7Kp that harbors bla_KPC-2, iii) pDCVA3-NDM-5 (>534.520 bp, IncFII) from the clinical strain K. pneumoniae HA31Kp that harbors bla_NDM-5, iv) pDCASG-NDM-1 (137.269 bp, IncC) from the clinical strain S. marcescens SM938 that harbors bla_NDM-1, and v) recombinant plasmid paadB from Escherichia coli TOP 10::paadB (5.877 bp, p15A) that harbors the aadB gene cassette. We identified by bioinformatics analysis that the four clinical plasmids had conjugation systems (Table 3), all of them being conjugative to E. coli J53 and E. coli M19736 in our experimental conditions (Figure 1). Also, E. coli M19736 was able to acquire paadB by chemical transformation (Figure 1A). Conjugation was performed as previously described, showing that E. coli M19736 was able to acquire the four plasmids (Figure 1B). Although all conjugation efficiencies were higher when E. coli M19736 was the receptor strain, a statistical difference between E. coli M19736 and E. coli J53 was only found when S. marcescens SM938 (pDCASG6-NDM-1) was used as a donor (Figure 1B).

Table 3

Table 3 Features of plasmids used for the experiments.

Figure 1

Figure 1 Conjugation/transformation assays of different multidrug-resistant plasmids with Escherichia coli M19736 and E. coli J53 as laboratory control. Panel (A) shows the acquisition of the different plasmids by conjugation or transformation of E. coli M19736. Each color represents the plasmids that were transferred: pDCAG1-CTX-M-15 (violet), pDCCK1-KPC (pink), pDCVA3-NDM-5 (blue), pDCASG-NDM-1 (orange), and paadB (yellow). Five transconjugants using E. coli M19736 as receptor (E. coli M19736::pDCAG1-CTX-M-15, E. coli M19736::pDCCK1-KPC, E. coli M19736:: pDCVA3-NDM-5, E. coli M19736::pDCASG-NDM-1, and evolved MDR-E. coli M19736) and transformant XDR-E. coli M19736 were generated. Panel (B) shows each donor and transconjugant (TRCJ) obtained in the different experiments using E. coli M19736 and E. coli J53 as receptors. Escherichia coli J53 as control is depicted as light gray with dotted bars, E. coli M19736 as black diagonal line bars, and evolved MDR-E. coli M19736 as dark gray with horizontal line bars. Data show the comparison of conjugation efficiencies of each plasmid between E. coli M19736 and E. coli J53. Conjugation efficiencies were calculated as the quotient between the number of transconjugants (Tc) and the number of donors (D) (Tc/D) in triplicates. Data show the median with range from three independent experiments (n = 3). Significant differences between groups are indicated: *p< 0.05; ns, not significant.

ARG acquisition by E. coli M19736 and E. coli J53 as receptor strains was verified by PCR of conjugation markers (Table 2), phenotypic detection of β-lactamases, AST (Supplementary Table 1), and MIC (Table 4). In the experiments of successive conjugation assays, the evolved MDR-E. coli M19736 was able to harbor simultaneously mcr-1, bla_CTX-M-15, and bla_NDM-1 genes (Figure 1A). The evolved XDR-E. coli M19736 was able to acquire also the aadB gene later by transformation assay (Figure 1A). Susceptibility tests showed that the evolved XDR-E. coli M19736 became resistant to all antibiotics tested except to trimethoprim/sulfamethoxazole (Supplementary Table 1). The MIC for carbapenem antibiotics (ERT and MEM) was slightly higher in E. coli M19736::pDCVA3-NDM-5, evolved MDR-E. coli M19736, and evolved XDR-E.coli M19736 than the respective donor strains and transconjugants of E. coli J53 (Table 4).

Table 4

Table 4 Minimum inhibitory concentration (MIC) of Escherichia coli strains.

3.2 Maintenance of plasmids in native and evolved MDR and XDR-Escherichia coli M19736 strains

Firstly, we evaluated the maintenance of plasmids of each donor cell (Escherichia coli SM5, K. pneumoniae HA7Kp, K. pneumoniae HA31Kp, and S. marcescens SM938) on the 1st and 10th days after being subcultured without antibiotic pressure at 37°C. All plasmids were maintained at 100% of each assay (data not shown). Then, each transconjugant of E. coli M19736, or transformant in the case of XDR-E. coli M19736 (Figure 1), was evaluated for its ability to maintain ARG (bla_CTX-M-15, bla_KPC-2, bla_NDM-5, bla_NDM-1, or aadB) by doing the same assay without antibiotic pressure at 37°C (Figure 2). At the same time, the maintenance of the mcr-1 gene was evaluated in all combinations. Each one of the four clinical plasmids (pDCAG1-CTX-M-15, pDCCK1-KPC, pDCVA3-NDM-5, and pDCASG6-NDM-1) harbored a different toxin/antitoxin system (Table 3). The transconjugants were able to maintain each gene of interest on the first day of subculture at 100%. Each transconjugant or transformant maintained bla_CTX-M-15 (pDCAG1-CTX-M-15), bla_KPC-2 (pDCCK1-KPC), bla_NDM-5 (pDCVA3-NDM-5), bla_NDM-1 (pDCASG6-NDM-1), and aadB (paadB) at 98.9%, 73.3%, 88.9%, 100%, and 0%, respectively, on the 10th day of subculture (Figure 2).

Figure 2

Figure 2 Maintenance of crucial ARG harbored by plasmids. Each color/pattern shows the percentage of maintenance of each ARG that was acquired in its respective plasmid by E. coli M19736 (A) or by E. coli J53 (B). The mcr-1 gene in pM19736-MCR-1, the bla_CTX-M-15 gene in pDCAG1-CTX-M-15, the bla_KPC-2 gene in pDCCK₁-KPC, the bla_NDM-5 gene in pDCVA3-NDM-5, the bla_NDM-1 gene in pDCASG-NDM-1, and the aadB gene in paadB were used as target for PCR detection with primers from Table 2, respectively. Data shows the mean ± SD from three independent experiments performed in triplicate (n = 3).

Interestingly, when the evolved MDR and XDR-E. coli M19736 strains acquired progressively bla_CTX-M-15 (pDCAG1-CTX-M-15) and bla_NDM-1 (pDCASG6-NDM-1) or acquired bla_CTX-M-15 (pDCAG1-CTX-M-15), bla_NDM-1 (pDCASG6-NDM-1), and aadB (paadB) plasmids, respectively, a different pattern of maintenance was found (Figure 3). In the case of evolved MDR-E. coli M19736, pDCAG1-CTX-M-15 and pDCASG-NDM-1 were maintained at 41.1% and 91.1%, respectively, on the 10th day of subculture (Figure 3). When the evolved XDR-E. coli M19736 that harbored bla_CTX-M-15 (pDCAG1-CTX-M-15), bla_NDM-1 (pDCASG6-NDM-1), and aadB (paadB) genes was tested, we found that bla_CTX-M-15 (pDCAG1-CTX-M-15) and aadB (paadB) genes were lost while maintaining the bla_NDM-1 gene (pDCASG6-NDM-1) at 98.9% on the 10th day of subculture. Remarkably, in all cases without antibiotic pressure, E. coli M19736 maintained the mcr-1 gene.

Figure 3

Figure 3 Maintenance of the ARG-plasmid located in the evolved MDR and XDR E. coli M19736 strains. The bars represent an evolved MDR-E. coli M19736 or evolved XDR-E. coli M19736 on the 1st day and 10th day of subculture. Each color/pattern and number within each bar represents the percentage of maintenance for the conjugated plasmids. Data shown are the mean ± SD from three independent experiments performed in triplicate (n = 3).

On the other hand, bioinformatics analysis of the WGS of the evolved XDR-E. coli M19736 from subculture on day 1 confirmed the presence of eight ARGs found in pM19736-MCR-1 and pIncFII-M19736 plasmids and chromosome from the host E. coli M19736 strain, five ARGs from pDCASG-NDM-1 as expected, and two ARGs from paadB (Table 5). The replication origins of these plasmids were also found. In contrast, no antibiotic determinant or replication origins of the pDCAG1-CTX-M-15 plasmid were found in the evolved XDR-E. coli M19736 which could be due to the fact that pDCAG1-CTX-M-15 was rapidly lost on day 1 as shown in our maintenance experiments (see below) (Figure 3).

Table 5

Table 5 ARGs found in the genome of the evolved XDR-E. coli M19736.

3.3 Ability of the evolved MDR and XDR-Escherichia coli M19736 strains to disseminate the acquired conjugative plasmids

The evolved MDR and XDR-E. coli M19736, which were co-infected with several plasmids, were tested as donors of clinical conjugative plasmids using again as receptor E. coli J53 (Figures 4A, B). The selection was performed with ceftazidime. The evolved MDR-E. coli M19736 and XDR-E. coli M19736 strains were able to transfer the bla_NDM-1 gene located in pDCASG-NDM-1 to E. coli J53 in the three independent biological replicates. On the other hand, when the evolved MDR-E. coli M19736 was used as donor, we identified two other genotypes in the transconjugants (Figure 4A). The first one transferred simultaneously the bla_NDM-1 (pDCASG6-NDM-1) and bla_CTX-M-15 (pDCAG1-CTX-M-15) genes in the transconjugants of E. coli J53. The second one harbored only the bla_CTX-M-15 (pDCAG1-CTX-M-15) gene. Furthermore, when the evolved XDR-E. coli M19736 was used as donor, we were able to find another genotype in which bla_NDM-1 (pDCASG6-NDM-1) and mcr-1 (pM19736-MCR-1) genes were detected simultaneously (Figure 4B).

Figure 4

Figure 4 Dissemination of plasmids co-infecting evolved Escherichia coli M19736 strains to E. coli J53. Evolved MDR-E. coli M19736 and XDR-E. coli M19736 were tested as donors of co-infecting clinical conjugative plasmids using as receptor E. coli J53 (A, B). The selection was performed with 8 µg/ml of ceftazidime. The number of colonies detected with each gene (bla_NDM-1, bla_CTX-M-15, and/or mcr-1) in transconjugants is shown. Numbers in parentheses represent the conjugation efficiencies of each experiment for the plasmids. Conjugation efficiencies were calculated as the quotient between the number of transconjugants (Tc) and the number of donors (D) (Tc/D). Both evolved strains were able to transfer the bla_NDM-1 gene located in pDCASG-NDM-1 to E. coli J53 from three independent biological replicates. When the evolved MDR-E. coli M19736 was used as donor, two other genotypes were identified in transconjugants in one replicate (A). The first genotype harbored simultaneously plasmids of the bla_NDM-1 (pDCASG-NDM-1) and bla_CTX-M-15 (pDCAG1-CTX-M-15) genes in E. coli J53 transconjugants. The second one harbored only the bla_CTX-M-15 (pDCAG1-CTX-M-15) gene. When the evolved XDR-E. coli M19736 was used as donor, in one replicate, another genotype was identified in which the bla_NDM-1 (pDCASG-NDM-1) and mcr-1 (pM19736-MCR-1) genes were detected simultaneously (B).

3.4 Isolation and characterization of EV from the native and evolved XDR-Escherichia coli M19736 strains

The native E. coli M19736 and evolved XDR-E. coli M19736 strains actively released EV at the log phase of growth and were isolated and collected from the supernatant broth. The cell-free EVs extracted from both strains were purified by filtration and ultracentrifugation. EVs were characterized in terms of morphology, size, and polydispersity index (PDI). The purified EVs appeared at TEM as electron-dense particles, with a spherical morphology, a bilayer membrane, and heterogeneous nanometer size (Figures 5B, D). The purity of the EV was confirmed as there were no bacteria visualized by TEM, and contamination controls on culture plates did not show any growth. This showed that EVs were purified successfully without contamination with other bacterial components for subsequent cell experiments. The obtained data from DLS showed that the typical diameter was approximately 193.2 ± 1.8 nm with a PDI of 0.199 ± 0.012 for the native E. coli M19736 (Figure 5A) and 174.7 ± 0.52 nm with a PDI of 0.2 ± 0.009 for the evolved XDR-E. coli M19736 (Figure 5C).

Figure 5

Figure 5 Phenotypic characterization of EV from native and evolved XDR-E. coli M19736. Vesicles were purified from broth cultures by ultracentrifugation and filtered through a 0.22-μm filter. DLS results are presented as the mean of three independent measurements ± SD. DLS measurement shows an average size of 193.2 ± 1.8 nm for the native E. coli M19736 (A) and 174.7 ± 0.52 nm for the evolved XDR-E. coli M19736 (C). TEM results show in both strains EV with a double membrane, spherical in shape and heterogeneous in size (B, D).

3.5 Content analysis of EV from both Escherichia coli M19736 and evolved XDR-Escherichia coli M19736 strains

The DNA purified from EVs of the native E. coli M19736 gave a specific amplified product for the mcr-1 gene. Also, EV DNA from the evolved XDR-E. coli M19736 allowed us to detect specific PCR products for the bla_CTX-M-15, mcr-1, and aadB genes. On the other hand, the total vesicular proteins were extracted via lysis buffer and then quantified by the Micro BCA kit. The amount of protein for the native and evolved strains was 1,015 μg/ml and 1,091 μg/ml, respectively. To explore the protein contents of EV from E. coli M19736, LC-MS/MS analysis was applied, and 338 different proteins were identified in EVs from the native strain (Supplementary Table 2). Database protein was included in Vesiclepedia 2024 (http://www.microvesicles.org) (Chitti et al., 2024). Proteins were categorized into different classes including the following: cellular localization site (Figure 6) and biological functions (Figure 7). The localization of EV proteins from E. coli M19736 was found to be distributed as follows: 10% of the proteins were located in the cell membrane, 18% in the inner membrane, 16% in the outer membrane, 50% in the cytoplasm, and 6% in the periplasm. Moreover, among the 338 identified proteins, we were able to characterize 292 of them by biological processes/functions (Figure 7). Some proteins have overlapping functions. The majority were involved in the transport, metabolism, and biosynthesis of molecules such as proteins, lipids, and carbohydrates (12.9%) and biosynthesis of secondary metabolites (10.1%). Moreover, the others were involved in microbial metabolism in diverse environments (8.7%); carbon metabolism and utilization (7.1%); tricarboxylic acid cycle/pyruvate metabolism (5.7%); amino acid biosynthesis, metabolism, and transport (5.4%); cell wall and peptidoglycan biosynthesis, metabolism, and degradation (5.1%); translation (4.4%); ion transport and storage (3.1%); transport and ABC transporters (3%); glycolysis/gluconeogenesis (2.7%); two-component system (2.7%); antibiotic response (2.6%); biosynthesis of cofactors (2.4%); nucleotide biosynthesis and metabolism (2.4%); stress response (4.1%); cell cycle and division (1.7%); oxidative phosphorylation (1.4%); rRNA, tRNA, and mRNA processing and degradation (1.4%); transcription (1.4%); vitamin biosynthesis, metabolism, and transport (1.3%); and DNA replication, recombination, repair, damage, and condensation (1.1%). The main functions and pathways enriched in our EV proteins have been related to the same as other EV protein cargoes of XDR K. pneumoniae (Hussein et al., 2023). Other functions were also found to be represented in less than one percent such as quorum sensing, lipopolysaccharide (LPS) and lipid A biosynthesis, biofilm formation, respiratory electron transport chain, glutathione metabolism, Gram-negative bacterium-type cell outer membrane assembly, nitrogen metabolism, heme and porphyrin biosynthesis and metabolism, glyoxylate and dicarboxylate metabolism, methane metabolism, organic substance metabolism, organic acid catabolism, and chemotaxis.

Figure 6

Figure 6 Predicted localization of EV proteins from the native Escherichia coli M19736. The pie diagram represents the localization of different proteins found inside or on the surface of EVs. The results are represented as the percentage of proteins found in different localizations from the native E. coli M19736.

Figure 7

Figure 7 Biological processes associated with EV proteins from the native E. coli M19736. Bar graphs categorizing the proteins with 28 differential biological functions, with distributions in the number of proteins. The functions of EV proteins are graphed with the most abundant function at the bottom and the least abundant function at the top.

Interestingly, proteins associated with the LPS biosynthesis pathway have been found, among others: the LPS assembly OM complex LptDE β-barrel component LptD and LptE (Supplementary Table 2), proteins that form a hetero-oligomeric complex that translocates LPS to the outer membrane and allows it to anchor to the cell wall surface (Lucena et al., 2023). Proteins that form efflux bombs, such as the three proteins that make up the tripartite efflux system AcrAB-TolC (Supplementary Table 2), have also been found. The AcrB and AcrA proteins, respectively, make up the inner membrane and periplasm-spanning regions, and the TolC protein component is located in the bacterial outer membrane and also pairs with subunits of other membrane pumps. In addition, we found other proteins that are involved in pathways of cationic antimicrobial peptide (CAMP) resistance such as D-transpeptidase linking Lpp to murein, N-acetylmuramoyl-L-alanine amidase AmiC, lipoprotein NlpE, and periplasmic serine endoprotease DegP (Supplementary Table 2). Lastly, we detected the MCR-1 protein (Supplementary Table 2). To the best of our knowledge, the presence of this protein in EV has not been described so far. MCR-1 mediates colistin resistance by transferring phosphoethanolamine to bacterial lipid A, thereby reducing its affinity for colistin (Li H. et al., 2021).

3.5.1 PPI network of EV proteins from Escherichia coli M19736

We constructed a PPI network using the STRING database and analyzed it using the Cytoscape software. Gene Ontology and KEEG enrichment analysis permitted to generate the network diagrams. Each node and line represented a term and the correlation between terms, respectively. The color of the terms indicates the classification of nodes based on their functions. We obtained three different PPI networks based on function/biological process/pathway found in the EV protein from E. coli M19736 (Figure 8). The first one showed cell wall and peptidoglycan biosynthesis and metabolism, represented by 27 nodes and 228 edges, with an average local clustering coefficient of 0.666 and PPI enrichment p-value of 1.15e−03. The second one showed an antibiotic response, represented by 23 nodes and 68 edges, and the average local clustering coefficient was 0.819 and the PPI enrichment p-value was 1.0e−16. The third one showed a pathway related to β-lactam resistance, represented by 7 nodes and 12 edges, with an average local clustering coefficient of 0.762 and PPI enrichment p-value of 1.94e−07.

Figure 8

Figure 8 PPI network of proteins of interest found in EV from the native E. coli M19736. Gene Ontology and KEEG enrichment analysis of genes/proteins from the proteomes of EVs derived from the native E. coli M19736 using the STRING software and visualization using the Cytoscape software. The node colors represent the biological process or cellular functions of the genes/proteins of interest according to significant associations of related Gene Ontology and KEEG terms. The arrow points to the MCR-1 protein.

Furthermore, the finding of the interaction of some specific proteins described above is of particular interest to our work. Some of these proteins could interact to generate antibiotic response/resistance, such as LptD and LptE proteins, which have been suggested to be involved in an antibiotic stress response leading to increased production and accumulation in the outer membrane (Figure 8) (Lucena et al., 2023). The efflux pump AcrAB-TolC is involved in conferring CAMP resistance in different bacteria, although this is controversial in E. coli (Blair et al., 2022). Overexpression of these efflux pumps also confers resistance to a variety of antibiotics (Brindangnanam et al., 2022). We also found several proteins involved in ribosomal and RNA degradation, including the chaperone Hsp70 (DnaK), and 30S ribosomal proteins that have been shown to interact with the MCR-1 protein (Supplementary Table 2) (Li H. et al., 2021). It should be noted that in addition to these proteins, the above mentioned proteins associated with CAMP resistance and the AcrA-TolC system are also important in the MCR-1 protein interactome (Li H. et al., 2021).

4 Discussion

HGT is a powerful force that shapes the evolution, diversification, and adaptation of bacterial communities and provides, for instance, a platform for the spread and persistence of ARGs (Piscon et al., 2023). Today, three canonical mechanisms of HGT are recognized, including transformation, transduction, and conjugation (Dubnau and Blokesch, 2019). Also, there are other non-canonical mechanisms including vesiduction, which involves secretion and uptake of EV (Li P. et al., 2022; Marinacci et al., 2023). In stressing habitats such as the nosocomial habitat, genome evolution is driven by antibiotic selection. However, there are several gaps of knowledge in this field. An area not well understood yet, it is what occurs with multidrug-resistant bacterial communities and plasmids that carry ARG in periods of time in which there is no antibiotic selection pressure. Also, the role of sporadic clones related to the spread of ARG is little studied to date. Here, we exposed the ability of a sporadic clone of E. coli, the E. coli M19736 strain belonging to ST615, harboring a plasmid with the mcr-1 gene, to acquire a wide variety of clinical conjugative plasmids, including sequential co-infection of three conjugative plasmids harboring ARGs of current clinical interest (mcr-1, bla_NDM-1, and bla_CTX-M-15) from different species of bacteria (Figure 1). In turn, its competency to disseminate the conjugative plasmids again to another bacterial host was shown (Figure 3). At the same time, other mechanisms that were tested in this strain such as transformation for the paadB plasmid and vesiduction of native and evolved E. coli M19736 strains were identified that could be relevant reservoirs. bla_CTX-M-15 was found in EV even if the evolved XDR-E. coli M19736 had rapidly lost the plasmid pDCAG1-CTX-M-15 after its acquisition, showing the essential role of EV for the dissemination of ARG in bacterial communities. Interestingly, the MICs for carbapenem antibiotics (ERT and MEM) were slightly higher in E. coli M19736::pDCVA3-NDM-5, evolved MDR-E. coli M19736, and evolved XDR-E. coli M19736 than the respective donor strains and transconjugants of E. coli J53, showing the ability of one sporadic clone to express the crucial ARG after HGT acquisition.

It has been a while since the population structure of E. coli has been identified; long-term stability and wide geographic distribution of individual lineages have been identified (Bobay et al., 2015). Some of these clones are pandemic, such as E. coli ST131, which is the predominant extraintestinal pathogenic E. coli that causes multidrug hospital infections, usually harboring bla_CTX-M-15 and/or carbapenemases (Mathers et al., 2015; Sanz et al., 2022). In a recent retrospective epidemiological study performed with 71 relevant carbapenem-resistant E. coli strains isolated from 2008 to 2017 from Argentina, several pandemic clones including ST10, ST38, ST131, ST155, ST648, and ST1193 were found prevalent (Sanz et al., 2022). From this bacterial population under scrutiny, three carbapenem-resistant E. coli strains harboring the mcr-1 gene were also found belonging to the pandemic clone ST10 and to sporadic clones ST12657 and ST12667. Interestingly, two of these strains (E. coli ECO 37 isolated in 2014 and E. coli ECO 81 isolated in 2015) were isolated before the first description of the mcr-1 gene in isolates from animals, food, and patients in China (Liu et al., 2016). Escherichia coli ST615 was not identified in those epidemiological studies and in other studies performed with E. coli strains from Argentina isolated from the clinic or other environments before or after the COVID-19 pandemic (Sennati et al., 2012; Dominguez et al., 2019; Faccone et al., 2023; Gramundi et al., 2023; Piekar et al., 2023), confirming the sporadic condition of this clone. Escherichia coli M19736 ST615 strain isolated in 2015 was identified at that time as one of the first isolates harboring the mcr-1 gene in human infections caused by E. coli in Latin America (Rapoport et al., 2016). The importance and scope of a wide variety of sporadic clones has not yet been studied in-depth. It is likely that E. coli ST615 lineage could capture crucial ARG such as the mcr-1 gene, with the ability to transfer consequently to other strains as shown in the present study and previously (Rapoport et al., 2016). We also identified that the E. coli M19736 ST615 strain was able to acquire a diversity of plasmids of different incompatibility groups harboring multiple ARGs from three different species (E. coli, K. pneumoniae, and S. marcescens) (Figure 1A). In addition, co-infection of plasmids sharing the same incompatibility group such as pDCVA3-NDM-5 and pIncFII-M19736 in E. coli M19736::pDCVA3-NDM-5 or pDCAG1-CTX-M-15 and pIncFII-M19736 in evolved MDR-E. coli M19736 strain was found (Table 1). Several differences were identified among the three IncFII replicons (Supplementary Table 3) that could be related, in part, to their ability to co-infect and to be maintained together during 10 days of daily subcultures by E. coli M19736.

Our studies showed that transconjugants of E. coli M19736 ST615 maintained pDCAG1-CTX-M-15 (IncFII) and pDCASG-NDM-1 (IncC) plasmids for 10 days at 100% while maintaining its native one, pMCR-M19736 (IncI2) with the mcr-1 gene, without antibiotic pressure (Figure 2A). The co-infection with pDCAG1-CTX-M-15 (IncFII) and pDCASG-NDM-1 (IncC) is interesting since both have shown to possess a pandemic behavior. On one hand, the pDCAG1-CTX-M-15 plasmid has an F2:B10 replicon; IncF-type replicons have shaped the evolution of the main fluoroquinolone-resistant ST131-H30 clades adding an advantage resistance to several families of antibiotics including the presence of the bla_CTX-M-15 gene (Johnson et al., 2016). On the other hand, pDCASG-NDM-1 belongs to the IncC group plasmids that are widely distributed among Gnb, with a large range of hosts in which these plasmids can replicate (Ambrose et al., 2018). IncC plasmids have islands of resistance incorporated in different plasmid locations where different ARGs can be accumulated including bla_NDM and bla_KPC genes (Ambrose et al., 2018; Ambrose, 2020). At first glance, since both pDCAG1-CTX-M-15 and pDCASG-NDM-1 plasmids had toxin–antitoxin systems (CcdA/B and pemK/L and HigA/B, respectively) and different incompatibility groups (IncFII and IncC, respectively), a low percentage of loss was expected during co-infected subcultures. Recent advances in the field profoundly questioned the role of toxin–antitoxin systems in bacterial physiology, stress response, and antimicrobial persistence (Jurėnas et al., 2022). More investigations are needed to evaluate their role in clinical isolates harboring several plasmids.

Recently, it has been shown that plasmids carrying a carbapenemase such as KPC or NDM could be efficiently conjugated to strains carrying the mcr-1 gene and vice versa and that these plasmids could stably co-exist in clinical Enterobacteriaceae strains (Liu et al., 2021). Although the clonality of these strains was not determined, these results are congruent with our experiments of the evolved E. coli M19736 strains (Figure 2), which were able to transfer in turn the conjugative plasmids acquired previously. In addition, our experiments showed that when the transconjugant E. coli M19736::pDCAG1-CTX-M-15 was co-infected with additional plasmids, our evolved MDR and XDR-E. coli M19736 ST615 strains showed different patterns of ARG maintenance with the ability to keep almost at 100% the bla_NDM-1 gene located in the pandemic IncC plasmid pDCASG-NDM-1.

Concerning co-infection of plasmids, it is generally overlooked that bacterial strains frequently harbor multiple plasmids, and understanding them is of utmost importance, especially for those relevant in the clinical context (Dionisio et al., 2019). Other studies have shown that bacteria can carry more than one type of plasmid; for example, it has been shown that 27 strains of E. coli producing extended-spectrum β-lactamase harbored multiple different plasmids (García et al., 2007). Positive epistasis between co-infecting plasmids has been shown which minimizes the cost of plasmid carriage and increases the ability of plasmids to persist in the absence of selection for plasmid-encoded traits, suggesting that epistasis may have an important role in resolving the “plasmid paradox” (San Millan et al., 2013), which is in agreement with our results. We also found that maintenance of plasmids without antibiotic selection varied depending on the plasmids that co-infected the host (Figure 2), such as the case of the evolved MDR-E. coli M19736 and XDR-E. coli M19736 strains related to the maintenance of pDCAG1-CTX-M-15 (Figure 3). The addition of plasmid paadB (Inc15A) triggered the loss of pDCAG1-CTX-M-15 from the evolved XDR-E. coli M19736 strain. Concerning this, a strong co-evolution has been shown between some E. coli ST131 lineages and specific plasmids including F2:B10 (Ny et al., 2019), which could explain in part its ease of getting lost in a sporadic clone such as E. coli ST615.

In order to follow the trajectory of plasmids that were co-infecting our evolved strains, conjugation assays were performed revealing that the evolved MDR and XDR-E. coli-M19736 strains were able to transfer one or two plasmids simultaneously while keeping always the pM19736 plasmid (Figure 4). Previous plasmid co-transfer in E. coli strains showed that when hosts harbor two, three, or four distinct plasmids, the co-transfer of both plasmids tends to be limited by the plasmid exhibiting the lowest conjugation rate (Gama et al., 2017; Darphorn et al., 2022). In disagreement with this, we did not find significant differences between the conjugation efficiencies of these plasmids (Figure 1B), so in the cases where the co-transfer happens, we cannot associate with the conjugation rate of each plasmid (Figure 4). Also, our results are congruent with the hypothesis that suggests that de-repression could occur simultaneously on co-resident plasmids. In fact, this could happen as a response to a common stimulus. The idea that a common stimulus triggers the transfer of one or the other plasmid could explain why one or the other plasmid is co-transferred as proposed by Gama et al. (2017). However, it is not clear in our study what stimulus could have triggered co-transfection in our experiments, but selection with ceftazidime may be studied deeply. The tendency of bacteria to maintain several plasmids simultaneously independently of the antibiotic pressure imposed on the pathogens has been shown (San Millan et al., 2013). There are two predictable ways for a strain to harbor multiple plasmids: i) through sequential plasmid acquisition as we showed in the experiments using the sporadic clone E. coli M19736 as the recipient strain that evolved to MDR-E. coli M19736 and XDR-E. coli M19736 strains (Figure 1A) and ii) by simultaneous transmission of several plasmids as we showed in the experiments that we performed using the evolved MDR and XDR-E. coli M19736 strains as donors to E. coli J53 (Figure 4). Our results showed that one strain belonging to a sporadic clone is able to be co-infected with several plasmids in both ways as previously suggested (Dionisio et al., 2019).

Since recent research has revealed that EVs may play an important role as a reservoir of genes and proteins associated with AMR contributing to this global threat through several other mechanisms (Ciofu et al., 2000; Rumbo et al., 2011; Lee et al., 2013; Chatterjee et al., 2017; Bielaszewska et al., 2020; Marchant et al., 2021), we also investigated the ability of the EV of the E. coli M19736 strain to harbor AMR determinants. First, we analyzed the EV of E. coli M19736 and evolved XDR-E. coli M19736 strains by PCR searching for the ARG. We found in both strains the presence of the mcr-1 gene which has been also previously described (Li X. et al., 2021). Interestingly, EVs from the evolved XDR-E. coli M19736 strains harbored bla_CTX-M-15 though the pDCAG1-CTX-M-15 was previously lost as shown by WGS analysis and maintenance experiments, suggesting that EVs could be a relevant reservoir of ARG for susceptible bacteria. Unexpectedly, we did not find the bla_NDM-1 gene in the EV content of the evolved strain by DNA PCR. This result was surprising because the bla_NDM gene is the most described gene reported in EV so far (Chatterjee et al., 2017; Li X. et al., 2021). These findings lead us to consider what might play a role in DNA cargo selection in EV, as the mechanism by which DNA cargo is enriched or excluded from EV remains unclear (Zavan et al., 2023). Plasmid type has been shown to influence gene loading (Tran and Boedicker, 2017). However, this is not the only reason explaining why bla_NDM-1 located in pDCASG6-NDM-1 was not found in EV. Other features could be involved in the selection of DNA loading in EV, such as the location of each type of plasmid within the cell (Pogliano, 2002). Currently, there are more studies focusing on the effects of growth stage on EV selection load. For example, analysis of EV produced by Helicobacter pylori reveals that bacterial growth stage affects the size, composition, and selection of protein load (Zavan et al., 2019). Other research showed that plasmids belonging to different incompatibility groups bind to different sites within the cell (Pogliano, 2002). Aktar et al. (2021) also performed some experiments to evaluate the incorporation of plasmid DNA into EV. The results showed that peptidoglycan defects increased plasmid sorting into EV. It is well known that the use of cell wall antibiotics increases EV production (Bauwens et al., 2017; Lucena et al., 2023). In summary, since EV DNA packaging is not a random process, more work is needed to understand the rules of plasmid packaging into vesicles.

Furthermore, in EVs from the E. coli M19736 strain, the proteins from the LPS biosynthesis pathway were found (Supplementary Table 2), and 2-dehydro-3-deoxyphosphooctonate aldolase (KdsA) is responsible for linking lipid A and core oligosaccharides through the synthesis of other proteins (Hussein et al., 2023). The experiments of Strohmaier et al. (1995) demonstrate that KdsA undergoes growth phase-dependent regulation at the transcriptional level, which again supports the idea that the growth phase is related to the selection of cargo in EV. Also, this protein increases following polymyxin B treatment (Hussein et al., 2023) which supports the idea that treatment with each antibiotic could collaborate in the selection of load in EV. Also, in EV from the E. coli M19736 strain, the MCR-1 protein was found (Supplementary Table 2). To our knowledge, the presence of this protein has not been previously described. It has been demonstrated that MCR-1-interacting proteins were mainly involved in ribosome and RNA degradation such as the 30S ribosomal subunit proteins S5 and S10 (Li H. et al., 2021). These proteins were also identified in our EV extractions and among other ribosomal proteins as well (Supplementary Table 2). Also, the two-component AcrA-TolC multidrug efflux pump interacts with the MCR-1 protein and is involved in mcr-1-mediated colistin resistance (Li H. et al., 2021). In addition, EVs have shown to have an effect of protection of colistin’s administration over other strains (Marchant et al., 2021) that could apply to the E. coli M19736 strain under scrutiny in the present work. For relevance in the nosocomial niche, the main functions and pathways enriched in our EV proteins were related to EV protein’s cargo of XDR K. pneumoniae strains as previously shown (Hussein et al., 2023).

Since a few clinical sporadic clones have shown to be prone to accept plasmids by transformation (Ramirez et al., 2011), deeply analyzing the genomic plasticity of E. coli M19736, artificial competency was also achieved. This assay evidenced the ability of this strain to accept the paadB plasmid by chemical competency of the native and evolved MDR-E. coli M19736 strains, which could be maintained by E. coli M19736 during at least 58 generations (Figure 2). Furthermore, sporadic clones were found to capture relevant features of environmental plasmids that later could spread within the nosocomial niche such as the case of pDCPR1 which has a mosaic structure with parts of successful phytopathogen plasmids such as pIII from pathogenic Xanthomonas albilineans and pPSV48C from another pathogenic species Pseudomonas savastanoi (Vilacoba et al., 2014). This ability to capture genetic platforms from the bacterial communities of the open environment adds a relevant role to sporadic clones in the framework necessary for the success of AMR in the nosocomial niche.

In summary, our studies showed that E. coli M19736 strains belonging to the sporadic clone ST615 have a range of valuable tools for survival under antibiotic pressure as we identified its genomic plasticity to acquire plasmids from various species with crucial ARGs, with proficiency to keep them and disseminate in turn by conjugation to other strains. Furthermore, its competence to produce EVs containing several valuable AMR determinants with the potential to be disseminated and its ability to receive plasmids by transformation as demonstrated in the present study indicated a great adaptation to the in-hospital environment. Based on this, sporadic clones could play a crucial role in a clinical context by adapting to different selection pressures and at the same time by spreading plasmids to other bacteria which guarantees an important link in the chain of events that leads bacteria to resist in the nosocomial niche.

Data availability statement

The E. coli M19736 genome data presented in the study is deposited in GenBank under the accession numbers PRJNA1080650 and PRJNA1102182. The assembled plasmids pDCASGNDM-1 and pDCAG1-CTX-M-15 were deposited in GeneBank under the accession numbers PP622803 and PP632099. The data base protein of EV was included in Vesiclepedia 2024 (http://www.microvesicles.org) study ID 3592.

Author contributions

LCP: Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. MO: Conceptualization, Investigation, Project administration, Writing – review & editing. AG: Data curation, Writing – review & editing. TP: Investigation, Methodology, Writing – review & editing. AA: Data curation, Investigation, Writing – review & editing. MQ: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Resources, Supervision, Writing – review & editing. DC: Conceptualization, Formal analysis, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing, Validation.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by grant ISID/Pfizer 2019 given to DC, grant PUE 2522 from CONICET given to IMPaM (UBA/CONICET), and PICT 2021 (1120) given to DC.

Acknowledgments

We are very grateful to Dr. Alejandro Petroni from INEI-ANLIS “Dr. Carlos G. Malbrán” for providing us with the Escherichia coli M19736. Also, we would like to thank Claudia Vanessa Vargas and Nicolas Donis for their help in the experimental assays. DC and MQ are career investigators of CONICET. LCP is a recipient of a doctoral fellowship from CONICET.

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.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

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.

Supplementary material

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

Abbreviations

ARG, antimicrobial resistance gene; EV, extracellular vesicle; ST, sequence type; XDR, extreme drug resistance; MDR, multidrug resistance; Gnb, Gram-negative bacilli; AMR, antimicrobial resistance; LB, Luria–Bertani; PDI, polydispersity index; LPS, lipopolysaccharide; CAMP, cationic antimicrobial peptide.

References

Aktar, S., Okamoto, Y., Ueno, S., Tahara, Y. O., Imaizumi, M., Shintani, M., et al. (2021). Incorporation of plasmid DNA into bacterial membrane vesicles by peptidoglycan defects in Escherichia coli. Front. Microbiol. 12. doi: 10.3389/fmicb.2021.747606

PubMed Abstract | CrossRef Full Text | Google Scholar

Ambrose, S. J. (2020). Properties and evolution of IncA and IncC plasmids and their interactions with Salmonella genomic island 1. Available at: https://ses.library.usyd.edu.au/handle/2123/24101 (Accessed 6 January 2024).

Google Scholar

Ambrose, S. J., Harmer, C. J., Hall, R. M. (2018). Evolution and typing of IncC plasmids contributing to antibiotic resistance in Gram-negative bacteria. Plasmid 99, 40–55. doi: 10.1016/j.plasmid.2018.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Aziz, R. K., Bartels, D., Best, A. A., DeJongh, M., Disz, T., Edwards, R. A., et al. (2008). The RAST Server: Rapid annotations using subsystems technology. BMC Genomics 9, 1–15. doi: 10.1186/1471-2164-9-75/TABLES/3

PubMed Abstract | CrossRef Full Text | Google Scholar

Álvarez, V. E., Carrera Páez, L., Piekar, M., García Allende, N., Campos, J., Mendiondo, N., et al. (2024). Genomic characterization of two NDM-5-producing isolates of Klebsiella pneumoniae ST11 from a single patient. Curr. Microbiol. 81 (3), 76. doi: 10.1007/s00284-023-03596-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Bauwens, A., Kunsmann, L., Karch, H., Mellmann, A., Bielaszewska, M. (2017). Antibiotic-mediated modulations of outer membrane vesicles in enterohemorrhagic Escherichia coli O104:H4 and O157:H7. Antimicrob. Agents Chemother. 61 (9), e00937–17. doi: 10.1128/AAC.00937-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Bielaszewska, M., Daniel, O., Karch, H., Mellmann, A. (2020). Dissemination of the blaCTX-M-15 gene among Enterobacteriaceae via outer membrane vesicles. J. Antimicrob. Chemother. 75, 2442–2451. doi: 10.1093/jac/dkaa214

PubMed Abstract | CrossRef Full Text | Google Scholar

Binsker, U., Käsbohrer, A., Hammerl, J. A. (2022). Global colistin use: a review of the emergence of resistant Enterobacterales and the impact on their genetic basis. FEMS Microbiol. Rev. 46, 1–37. doi: 10.1093/femsre/fuab049

CrossRef Full Text | Google Scholar

Blair, J. M. A., Zeth, K., Bavro, V. N., Sancho-Vaello, E. (2022). The role of bacterial transport systems in the removal of host antimicrobial peptides in Gram-negative bacteria. FEMS Microbiol. Rev. 46 (6), fuac032. doi: 10.1093/femsre/fuac032

PubMed Abstract | CrossRef Full Text | Google Scholar

Bobay, L. M., Traverse, C. C., Ochman, H. (2015). Impermanence of bacterial clones. Proc. Natl. Acad. Sci. United States America 112, 8893. doi: 10.1073/pnas.1501724112

CrossRef Full Text | Google Scholar

Boutzoukas, A. E., Komarow, L., Chen, L., Hanson, B., Kanj, S. S., Liu, Z., et al. (2023). International epidemiology of carbapenemase-producing Escherichia coli. Clin. Infect. Dis. 77, 499–509. doi: 10.1093/cid/ciad288

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandt, C., Viehweger, A., Singh, A., Pletz, M. W., Wibberg, D., Kalinowski, J., et al. (2019). Assessing genetic diversity and similarity of 435 KPC-carrying plasmids. Sci. Rep. 9, 11223. doi: 10.1038/S41598-019-47758-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Brindangnanam, P., Sawant, A. R., Prashanth, K., Coumar, M. S. (2022). Bacterial effluxome as a barrier against antimicrobial agents: structural biology aspects and drug targeting. Tissue Barriers 10 (4), 2013695. doi: 10.1080/21688370.2021.2013695

PubMed Abstract | CrossRef Full Text | Google Scholar

Caniaux, I., van Belkum, A., Zambardi, G., Poirel, L., Gros, M. F. (2017). MCR: modern colistin resistance. Eur. J. Clin. Microbiol. Infect. Dis. 36, 415–420. doi: 10.1007/s10096-016-2846-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Carattoli, A., Hasman, H. (2020). PlasmidFinder and in silico pMLST: identification and typing of plasmid replicons in whole-genome sequencing (WGS). Methods Mol. Biol. (Clifton NJ) 2075, 285–294. doi: 10.1007/978-1-4939-9877-7_20

CrossRef Full Text | Google Scholar

Cejas, D., Elena, A., González-Espinosa, F. E., Pallecchi, L., Vay, C., Rossolini, G. M., et al. (2022). Characterisation of blaKPC-2-harbouring plasmids recovered from Pseudomonas aeruginosa ST654 and ST235 high-risk clones. J. Global Antimicrob. Resist. 29, 310–312. doi: 10.1016/j.jgar.2022.04.017

CrossRef Full Text | Google Scholar

Chatterjee, S., Mondal, A., Mitra, S., Basu, S. (2017). Acinetobacter baumannii transfers the blaNDM-1 gene via outer membrane vesicles. J. Antimicrob. Chemother. 72, 2201–2207. doi: 10.1093/jac/dkx131

PubMed Abstract | CrossRef Full Text | Google Scholar

Chitti, S. V., Gummadi, S., Kang, T., Shahi, S., Marzan, A. L., Nedeva, C., et al. (2024). Vesiclepedia 2024: an extracellular vesicles and extracellular particles repository. Nucleic Acids Res. 52, D1694. doi: 10.1093/nar/gkad1007

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciofu, O., Beveridge, T. J., Kadurugamuwa, J., Walther-Rasmussen, J., Høiby, N. (2000). Chromosomal beta-lactamase is packaged into membrane vesicles and secreted from Pseudomonas aeruginosa. J. Antimicrob. Chemother. 45, 9–13. doi: 10.1093/jac/45.1.9

PubMed Abstract | CrossRef Full Text | Google Scholar

Clausen, P. T. L. C., Aarestrup, F. M., Lund, O. (2018). Rapid and precise alignment of raw reads against redundant databases with KMA. BMC Bioinf. 19 (1), 307. doi: 10.1186/s12859-018-2336-6

CrossRef Full Text | Google Scholar

CLSI. (2023). Performance Standards for Antimicrobial Susceptibility Testing, supplement M100. 33nd ed. (USA: Clinical and Laboratory Standards Institute). Available at: https://clsi.org/media/tc4b1paf/m10033_samplepages-1.pdf.

Google Scholar

Darphorn, T. S., Koenders-van Sintanneland, B. B., Grootemaat, A. E., van der Wel, N. N., Brul, S., Ter Kuile, B. H. (2022). Transfer dynamics of multi-resistance plasmids in Escherichia coli isolated from meat. PloS One 17 (7), e0270205. doi: 10.1371/journal.pone.0270205

PubMed Abstract | CrossRef Full Text | Google Scholar

Denamur, E., Clermont, O., Bonacorsi, S., Gordon, D. (2021). The population genetics of pathogenic Escherichia coli. Nat. Rev. Microbiol. 19, 37–54. doi: 10.1038/s41579-020-0416-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dionisio, F., Zilhão, R., Gama, J. A. (2019). Interactions between plasmids and other mobile genetic elements affect their transmission and persistence. Plasmid 102, 29–36. doi: 10.1016/j.plasmid.2019.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Dominguez, J. E., Faccone, D., Tijet, N., Gomez, S., Corso, A., Fernández-Miyakawa, M. E., et al. (2019). Characterization of Escherichia coli Carrying mcr- 1-Plasmids Recovered From Food Animals From Argentina. Front. Cell. Infect. Microbiol. 9. doi: 10.3389/fcimb.2019.00041

CrossRef Full Text | Google Scholar

Doncheva, N. T., Morris, J. H., Gorodkin, J., Jensen, L.J. (2019). Cytoscape StringApp: network analysis and visualization of proteomics data. J. Proteome Res. 18, 623–632. doi: 10.1021/acs.jproteome.8b00702

PubMed Abstract | CrossRef Full Text | Google Scholar

Dorward, D. W., Garon, C. F., Judd, R. C. (1989). Export and intercellular transfer of DNA via membrane blebs of Neisseria gonorrhoeae. J. Bacteriol. 171, 2499–2505. doi: 10.1128/jb.171.5.2499-2505.1989

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubnau, D., Blokesch, M. (2019). Mechanisms of DNA uptake by naturally competent bacteria. Annu. Rev. Genet. 53, 217–237. doi: 10.1146/annurev-genet-112618-043641

PubMed Abstract | CrossRef Full Text | Google Scholar

Faccone, D., Gomez, S. A., de Mendieta, J. M., Sanz, M. B., Echegorry, M., Albornoz, E., et al. (2023). Emergence of hyper-epidemic clones of enterobacterales clinical isolates co-producing KPC and metallo-beta-lactamases during the COVID-19 pandemic. Pathog. (Basel Switzerland) 12 (3), 479. doi: 10.3390/pathogens12030479

CrossRef Full Text | Google Scholar

Gama, J. A., Zilhão, R., Dionisio, F. (2017). Co-resident plasmids travel together. Plasmid 93, 24–29. doi: 10.1016/j.plasmid.2017.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Gambino, A. S., Déraspe, M., Álvarez, V. E., Quiroga, M. P., Corbeil, J., Roy, P. H., et al. (2021). Serratia marcescens SCH909 as reservoir and source of genetic elements related to wide dissemination of antimicrobial resistance mechanisms. FEMS Microbiol. Lett. 368 (14), fnab086. doi: 10.1093/femsle/fnab086

PubMed Abstract | CrossRef Full Text | Google Scholar

García, A., Navarro, F., Miró, E., Villa, L., Mirelis, B., Coll, P., et al. (2007). Acquisition and diffusion of blaCTX-M-9 gene by R478-IncHI2 derivative plasmids. FEMS Microbiol. Lett. 271, 71–77. doi: 10.1111/fml.2007.271.issue-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Geurtsen, J., de Been, M., Weerdenburg, E., Zomer, A., McNally, A., Poolman, J. (2022). Genomics and pathotypes of the many faces of Escherichia coli. FEMS Microbiol. Rev. 46 (6), fuac031. doi: 10.1093/femsre/fuac031

PubMed Abstract | CrossRef Full Text | Google Scholar

Global antimicrobial resistance and use surveillance system (GLASS) report: 2022. (2022). Global antimicrobial resistance and use surveillance system (GLASS) report: 2022 (World Health Organization). Available at: https://www.who.int/publications/i/item/9789240062702 (Accessed 5 March 2023).

Google Scholar

Gramundi, I., Albornoz, E., Boutureira, M., Rapoport, M., Gomez, S., Corso, A., et al. (2023). Characterization of third generation cephalosporin-resistant Escherichia coli clinical isolates from Ushuaia, Argentina. Rev. Argent. Microbiol. 55, 43–48. doi: 10.1016/j.ram.2022.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Gravel, A., Messier, N., Roy, P. H. (1998). Point mutations in the integron integrase IntI1 that affect recombination and/or substrate recognition. J. Bacteriol. 180, 5437. doi: 10.1128/JB.180.20.5437-5442.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

Guan, J., Chen, Y., Goh, Y. X., Wang, M., Tai, C., Deng, Z., et al. (2023). TADB 3.0: an updated database of bacterial toxin-antitoxin loci and associated mobile genetic elements. Nucleic Acids Res. 52 (D1), D784–D790. doi: 10.1093/nar/gkad962

CrossRef Full Text | Google Scholar

Hussein, N. H., Al-Kadmy, I. M. S., Taha, B. M., Hussein, J. D. (2021). Mobilized colistin resistance (mcr) genes from 1 to 10: a comprehensive review. Mol. Biol. Rep. 48, 2897–2907. doi: 10.1007/s11033-021-06307-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Hussein, M., Jasim, R., Gocol, H., Baker, M., Thombare, V. J., Ziogas, J., et al. (2023). Comparative proteomics of outer membrane vesicles from polymyxin-susceptible and extremely drug-resistant Klebsiella pneumoniae. mSphere 8 (1), e0053722. doi: 10.1128/msphere.00537-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Jamborova, I., Dolejska, M., Vojtech, J., Guenther, S., Uricariu, R., Drozdowska, J., et al. (2015). Plasmid-mediated resistance to cephalosporins and fluoroquinolones in various Escherichia coli sequence types isolated from rooks wintering in Europe. Appl. Environ. Microbiol. 81, 648–657. doi: 10.1128/AEM.02459-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, T. J., Danzeisen, J. L., Youmans, B., Case, K., Llop, K., Munoz-Aguayo, J., et al. (2016). Separate F-type plasmids have shaped the evolution of the H30 subclone of Escherichia coli sequence type 131. mSphere 1 (4), e00121–16. doi: 10.1128/mSphere.00121-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Jurėnas, D., Fraikin, N., Goormaghtigh, F., Van Melderen, L. (2022). Biology and evolution of bacterial toxin-antitoxin systems. Nat. Rev. Microbiol. 20, 335–350. doi: 10.1038/S41579-021-00661-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J. H., Lee, J., Park, J., Gho, Y. S. (2015). Gram-negative and Gram-positive bacterial extracellular vesicles. Semin. Cell Dev. Biol. 40, 97–104. doi: 10.1016/j.semcdb.2015.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Knecht, C. A., García Allende, N., Álvarez, V. E., Prack McCormick, B., Massó, M. G., Piekar, M., et al. (2022). Novel insights related to the rise of KPC-producing Enterobacter cloacae complex strains within the nosocomial niche. Front. Cell. Infect. Microbiol. 12. doi: 10.3389/fcimb.2022.951049

CrossRef Full Text | Google Scholar

Lee, J., Lee, E. Y., Kim, S. H., Kim, D. K., Park, K. S., Kim, K. P., et al. (2013). Staphylococcus aureus extracellular vesicles carry biologically active β-lactamase. Antimicrob. Agents Chemother. 57, 2589–2595. doi: 10.1128/AAC.00522-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Xie, Y., Liu, M., Tai, C., Sun, J., Deng, Z., et al. (2018). oriTfinder: a web-based tool for the identification of origin of transfers in DNA sequences of bacterial mobile genetic elements. Nucleic Acids Res. 46, W229. doi: 10.1093/nar/gky352

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Wang, Y., Chen, Q., Xia, X., Shen, J., Wang, Y., et al. (2021). Identification of functional interactome of colistin resistance protein MCR-1 in Escherichia coli. Front. Microbiol. 11. doi: 10.3389/fmicb.2020.583185

CrossRef Full Text | Google Scholar

Li, X., Sun, L., Li, C., Yang, X., Wang, X., Hu, X., et al. (2021). The Attenuated Protective Effect of Outer Membrane Vesicles Produced by a mcr-1 Positive Strain on Colistin Sensitive Escherichia coli. Front. Cell. Infect. Microbiol. 11. doi: 10.3389/fcimb.2021.701625

CrossRef Full Text | Google Scholar

Li, C., Jiang, X., Yang, T., Ju, Y., Yin, Z., Yue, L., et al. (2022). Genomic epidemiology of carbapenemase-producing Klebsiella pneumoniae in China. Genom. Proteomics Bioinf. 20, 1154–1167. doi: 10.1016/j.gpb.2022.02.005

CrossRef Full Text | Google Scholar

Li, P., Luo, W., Xiang, T. X., Jiang, Y., Liu, P., Wei, D. D., et al. (2022). Horizontal gene transfer via OMVs co-carrying virulence and antimicrobial-resistant genes is a novel way for the dissemination of carbapenem-resistant hypervirulent Klebsiella pneumoniae. Front. Microbiol. 13, 945972. doi: 10.3389/fmicb.2022.945972

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Chan, E. W. C., Chen, S. (2021). Transmission and stable inheritance of carbapenemase gene (blaKPC-2 or blaNDM-1)-encoding and mcr-1-encoding plasmids in clinical Enterobacteriaceae strains. J. Global Antimicrob. Resist. 26, 255–261. doi: 10.1016/j.jgar.2021.05.022

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucena, A. C. R., Ferrarini, M. G., de Oliveira, W. K., Marcon, B. H., Morello, L. G., Alves, L. R., et al. (2023). Modulation of Klebsiella pneumoniae Outer Membrane Vesicle Protein Cargo under Antibiotic Treatment. Biomedicines 11 (6), 1515. doi: 10.3390/biomedicines11061515

PubMed Abstract | CrossRef Full Text | Google Scholar

Maamar, E., Ferjani, S., Jendoubi, A., Hammami, S., Hamzaoui, , Mayonnove-Coulange, Z. L., et al. (2016). High prevalence of gut microbiota colonization with broad-spectrum cephalosporin resistant enterobacteriaceae in a Tunisian intensive care unit. Front. Microbiol. 7. doi: 10.3389/fmicb.2016.01859

PubMed Abstract | CrossRef Full Text | Google Scholar

Magiorakos, A. P., Srinivasan, A., Carey, R. B., Carmeli, Y., Falagas, M. E., Giske, C. G., et al. (2012). Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infect.: Off. Publ. Eur. Soc. Clin. Microbiol. Infect. Dis. 18, 268–281. doi: 10.1111/j.1469-0691.2011.03570.x

CrossRef Full Text | Google Scholar

Marchant, P., Carreño, A., Vivanco, E., Silva, A., Nevermann, J., Otero, C., et al. (2021). “One for All”: Functional Transfer of OMV-Mediated Polymyxin B Resistance from Salmonella enterica sv. Typhi ΔtolR and ΔdegS to Susceptible Bacteria. Front. Microbiol. 12. doi: 10.3389/fmicb.2021.672467

PubMed Abstract | CrossRef Full Text | Google Scholar

Marinacci, B., Krzyżek, P., Pellegrini, B., Turacchio, G., Grande, R. (2023). Latest update on outer membrane vesicles and their role in horizontal gene transfer: A mini-review. Membranes 13, 860. doi: 10.3390/membranes13110860

PubMed Abstract | CrossRef Full Text | Google Scholar

Mathers, A. J., Peirano, G., Pitout, J. D. D. (2015). Escherichia coli ST131: The quintessential example of an international multiresistant high-risk clone. Adv. Appl. Microbiol. 90, 109–154. doi: 10.1016/bs.aambs.2014.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsumura, Y., Peirano, G., Pitout, J. D. D. (2018). Complete genome sequence of Escherichia coli J53, an azide-resistant laboratory strain used for conjugation experiments. Genome Announc. 6. doi: 10.1128/genomeA.00433-18

CrossRef Full Text | Google Scholar

McMillan, H. M., Kuehn, M. J. (2021). The extracellular vesicle generation paradox: a bacterial point of view. EMBO J. 40. doi: 10.15252/embj.2021108174

CrossRef Full Text | Google Scholar

Molina, K. C., Bajrovic, V., Bonomo, R. A. (2022). Carbapenemase-producing Enterobacterales in solid organ transplantation: Tip of the iceberg? Transplant. Infect. Dis. 24 (2), e13780. doi: 10.1111/tid.13780

CrossRef Full Text | Google Scholar

Ny, S., Sandegren, L., Salemi, M., Giske, C. G. (2019). Genome and plasmid diversity of Extended-Spectrum β-Lactamase-producing Escherichia coli ST131 – tracking phylogenetic trajectories with Bayesian inference. Sci. Rep. 9, 1–9. doi: 10.1038/s41598-019-46580-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ojer-Usoz, E., González, D., Vitas, A. I. (2017). Clonal diversity of ESBL-producing Escherichia coli isolated from environmental, human and food samples. Int. J. Environ. Res. Public Health 14 (7), 676. doi: 10.3390/ijerph14070676

PubMed Abstract | CrossRef Full Text | Google Scholar

Parmeciano Di Noto, G., Jara, E., Iriarte, A., Centrón, D., Quiroga, C. (2016). Genome analysis of a clinical isolate of Shewanella sp. Uncovered an active hybrid integrative and conjugative element carrying an integron platform inserted in a novel genomic locus. Microbiol. (United Kingdom) 162, 1335–1345. doi: 10.1099/MIC.0.000310/CITE/REFWORKS

CrossRef Full Text | Google Scholar

Pasteran, F., Gonzalez, L. J., Albornoz, E., Bahr, G., Vila, A. J., Corso, A. (2016). Triton hodge test: improved protocol for modified Hodge test for enhanced detection of NDM and other carbapenemase producers. J. Clin. Microbiol. 54, 640. doi: 10.1128/JCM.01298-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Piekar, M., Álvarez, V. E., Knecht, C., Leguina, C., García Allende, N., Carrera Páez, L., et al. (2023). Genomic data reveals the emergence of the co-occurrence of blaKPC-2 and blaCTX-M-15 in an Escherichia coli ST648 strain isolated from rectal swab within the framework of hospital surveillance. J. Global Antimicrob. Resist. 32, 108–112. doi: 10.1016/j.jgar.2022.12.012

CrossRef Full Text | Google Scholar

Piscon, B., Pia Esposito, E., Fichtman, B., Samburski, G., Efremushkin, L., Amselem, S., et al. (2023). The effect of outer space and other environmental cues on bacterial conjugation. Microbiol. Spectr. 11 (3), e0368822. doi: 10.1128/spectrum.03688-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Pitout, J. D. D., Chen, L. (2023). The significance of epidemic plasmids in the success of multidrug-resistant drug pandemic extraintestinal pathogenic Escherichia coli. Infect. Dis. Ther. 12, 1029. doi: 10.1007/s40121-023-00791-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Pitout, J. D. D., DeVinney, R. (2017). Escherichia coli ST131: A multidrug-resistant clone primed for global domination. F1000Research 6, 1–7. doi: 10.12688/f1000research.10609.1

CrossRef Full Text | Google Scholar

Pogliano, J. (2002). Dynamic cellular location of bacterial plasmids. Curr. Opin. Microbiol. 5, 586–590. doi: 10.1016/S1369-5274(02)00370-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Quiroga, M. P. (2012). Estudio molecular y funcional de los sitios de recombinación atípicos de los integrones. Tesis doctoral (Universidad de Buenos aires). Available at: http://repositoriouba.sisbi.uba.ar/gsdl/collect/afamaster/index/assoc/HWA_5917.dir/5917.PDF (Accessed 9 January 2024).

Google Scholar

Ramirez, M. S., Adams, M. D., Bonomo, R. A., Centrón, D., Tolmasky, M. E. (2011). Genomic analysis of Acinetobacter baumannii A118 by comparison of optical maps: identification of structures related to its susceptibility phenotype. Antimicrob. Agents Chemother. 55, 1520–1526. doi: 10.1128/AAC.01595-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramírez, D. G., Nicola, F., Zarate, S., Relloso, S., Smayevsky, J., Arduino, S. (2013). Emergence of Pseudomonas aeruginosa with KPC-type carbapenemase in a teaching hospital: an 8-year study. J. Med. Microbiol. 62, 1565–1570. doi: 10.1099/jmm.0.059923-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Rapoport, M., Faccone, D., Pasteran, F., Ceriana, P., Albornoz, E., Petroni, A., et al. (2016). First Description of mcr-1-Mediated Colistin Resistance in Human Infections Caused by Escherichia coli in Latin America. Antimicrob. Agents Chemother. 60, 4412 LP–4413. doi: 10.1128/AAC.00573-16.Address

CrossRef Full Text | Google Scholar

Rebelo, A. R., Bortolaia, V., Kjeldgaard, J. S., Pedersen, S. K., Leekitcharoenphon, P., Hansen, I. M., et al. (2018). Multiplex PCR for detection of plasmid-mediated colistin resistance determinants, mcr-1, mcr-2, mcr-3, mcr-4 and mcr-5 for surveillance purposes. Eurosurveillance 23, 17–00672. doi: 10.2807/1560-7917.ES.2018.23.6.17-00672/CITE/PLAINTEXT

CrossRef Full Text | Google Scholar

Rello, J., Kalwaje Eshwara, V., Lagunes, L., Alves, J., Wunderink, R. G., Conway-Morris, A., et al. (2019). A global priority list of the TOp TEn resistant Microorganisms (TOTEM) study at intensive care: a prioritization exercise based on multi-criteria decision analysis. Eur. J. Clin. Microbiol. Infect. Dis. 38, 319–323. doi: 10.1007/s10096-018-3428-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Riccobono, E., Salvetti, S., Coppi, M., Montenora, I., Di Pilato, V., Rossolini, G. M. (2023). Citrobacter freundii resistant to novel β-lactamase inhibitor combinations and cefiderocol, co-producing class A, B and D carbapenemases encoded by transferable plasmids. J. Antimicrob. Chemother. 78, 1677–1682. doi: 10.1093/jac/dkad150

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodríguez, C., Brengi, S., Cáceres, M. A., Mochi, S., Viñas, M. R., Rizza, C. A., et al. (2018). Successful management with fosfomycin + ceftazidime of an infection caused by multiple highly-related subtypes of multidrug-resistant and extensively drug-resistant KPC-producing Serratia marcescens. Int. J. Antimicrob. Agents 52, 737–739. doi: 10.1016/J.IJANTIMICAG.2018.07.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Rumbo, C., Fernández-Moreira, E., Merino, M., Poza, M., Mendez, J. A., Soares, N. C., et al. (2011). Horizontal transfer of the OXA-24 carbapenemase gene via outer membrane vesicles: A new mechanism of dissemination of carbapenem resistance genes in Acinetobacter baumannii. Antimicrob. Agents Chemother. 55, 3084–3090. doi: 10.1128/AAC.00929-10

PubMed Abstract | CrossRef Full Text | Google Scholar

San Millan, A., Heilbron, K., MacLean, R. C. (2013). Positive epistasis between co-infecting plasmids promotes plasmid survival in bacterial populations. ISME J. 8, 601–612. doi: 10.1038/ismej.2013.182

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanz, M. B., De Belder, D., de Mendieta, J. M., Faccone, D., Poklepovich, T., Lucero, C., et al. (2022). Carbapenemase-producing extraintestinal pathogenic Escherichia coli from Argentina: clonal diversity and predominance of hyperepidemic clones CC10 and CC131. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.830209

CrossRef Full Text | Google Scholar

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinf. (Oxford England) 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

CrossRef Full Text | Google Scholar

Sennati, S., Santella, G., Di Conza, J., Pallecchi, L., Pino, M., Ghiglione, B., et al. (2012). Changing epidemiology of extended-spectrum β-lactamases in Argentina: emergence of CTX-M-15. Antimicrob. Agents Chemother. 56, 6003–6005. doi: 10.1128/AAC.00745-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Soncini, J. G. M., Cerdeira, L., Sano, E., Koga, V. L., Tizura, A. T., Tano, Z. N., et al. (2022). Genomic insights of high-risk clones of ESBL-producing Escherichia coli isolated from community infections and commercial meat in southern Brazil. Sci. Rep. 12 (1), 9354. doi: 10.1038/S41598-022-13197-Y

PubMed Abstract | CrossRef Full Text | Google Scholar

Strohmaier, H., Remler, P., Renner, W., Högenauer, G. (1995). Expression of Genes kdsA and kdsB Involved in 3-Deoxy-D-manno-Octulosonic Acid Metabolism and Biosynthesis of Enterobacterial Lipopolysaccharide Is Growth Phase Regulated Primarily at the Transcriptional Level in Escherichia coli K-12. J. Bacteriol. 177, 4488–4500. doi: 10.1128/jb.177.15.4488-4500.1995

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Morris, J. H., Cook, H., Kuhn, M., Wyder, S., Simonovic, M., et al. (2017). The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 45, D362. doi: 10.1093/nar/gkw937

PubMed Abstract | CrossRef Full Text | Google Scholar

Tijet, N., Faccone, D., Rapoport, M., Seah, C., Pasterán, F., Ceriana, P., et al. (2017). Molecular characteristics of mcr-1-carrying plasmids and new mcr-1 variant recovered from polyclonal clinical Escherichia coli from Argentina and Canada. PloS One 12, 1–13. doi: 10.1371/journal.pone.0180347

CrossRef Full Text | Google Scholar

Toyofuku, M., Schild, S., Kaparakis-Liaskos, M., Eberl, L. (2023). Composition and functions of bacterial membrane vesicles. Nat. Rev. Microbiol. 21, 415–430. doi: 10.1038/s41579-023-00875-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Tran, F., Boedicker, J. Q. (2017). Genetic cargo and bacterial species set the rate of vesicle-mediated horizontal gene transfer. Sci. Rep. 7, 1–10. doi: 10.1038/s41598-017-07447-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Vilacoba, E., Quiroga, C., Pistorio, M., Famiglietti, A., Rodríguez, H., Kovensky, J., et al. (2014). A blaVIM-2 plasmid disseminating in extensively drug-resistant clinical Pseudomonas aeruginosa and Serratia marcescens isolates. Antimicrob. Agents Chemother. 58, 7017–7018. doi: 10.1128/AAC.02934-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, T., Joshi, B., Janice, J., Askarian, F., Škalko-Basnet, N., Hagestad, O. C., et al. (2018). Enterococcus faecium produces membrane vesicles containing virulence factors and antimicrobial resistance related proteins. J. Proteomics 187, 28–38. doi: 10.1016/j.jprot.2018.05.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Goh, Y. X., Tai, C., Wang, H., Deng, Z., Ou, H. Y. (2022). VRprofile2: detection of antibiotic resistance-associated mobilome in bacterial pathogens. Nucleic Acids Res. 50, W768–W773. doi: 10.1093/nar/gkac321

PubMed Abstract | CrossRef Full Text | Google Scholar

Wick, R. R., Judd, L. M., Gorrie, C. L., Holt, K. E. (2017). Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PloS Comput. Biol. 13, 1–22. doi: 10.1371/journal.pcbi.1005595

CrossRef Full Text | Google Scholar

Zankari, E., Hasman, H., Cosentino, S., Vestergaard, M., Rasmussen, S., Lund, O., et al. (2012). Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 67, 2640–2644. doi: 10.1093/jac/dks261

PubMed Abstract | CrossRef Full Text | Google Scholar

Zavan, L., Bitto, N. J., Johnston, E. L., Greening, D. W., Kaparakis-Liaskos, M. (2019). Helicobacter pylori growth stage determines the size, protein composition, and preferential cargo packaging of outer membrane vesicles. Proteomics 19 (1-2), e1800209. doi: 10.1002/PMIC.201800209

PubMed Abstract | CrossRef Full Text | Google Scholar

Zavan, L., Fang, H., Johnston, E. L., Whitchurch, C., Greening, D. W., Hill, A. F., et al. (2023). The mechanism of Pseudomonas aeruginosa outer membrane vesicle biogenesis determines their protein composition. Proteomics 23 (10), e2200464. doi: 10.1002/pmic.202200464

PubMed Abstract | CrossRef Full Text | Google Scholar