Front. Cell. Infect. Microbiol. Frontiers in Cellular and Infection Microbiology Front. Cell. Infect. Microbiol. 2235-2988 Frontiers Media S.A. 10.3389/fcimb.2023.1269732 Cellular and Infection Microbiology Original Research Development of a high-throughput platform to measure plasmid transfer frequency Arbé-Carton Kepa 1 Rey-Sogo Ana 1 Santos-Fernández Nagore 1 Altube Oihane 1 Garbisu Carlos 2 Arana Lide 3 Alkorta Itziar 1 * 1 Department of Biochemistry and Molecular Biology, University of the Basque Country (UPV/EHU), Bilbao, Spain 2 Department of Conservation of Natural Resources, NEIKER-Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Spain 3 Department of Applied Chemistry, University of the Basque Country (UPV/EHU), Donostia, Spain

Edited by: Ibrahim Bitar, Charles University, Czechia

Reviewed by: Vittoria Mattioni Marchetti, University of Pavia, Italy; Gloria Del Solar, Spanish National Research Council (CSIC), Spain

*Correspondence: Itziar Alkorta, itzi.alkorta@ehu.eus

11 10 2023 2023 13 1269732 30 07 2023 25 09 2023 Copyright © 2023 Arbé-Carton, Rey-Sogo, Santos-Fernández, Altube, Garbisu, Arana and Alkorta 2023 Arbé-Carton, Rey-Sogo, Santos-Fernández, Altube, Garbisu, Arana and Alkorta

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Antibiotic resistance represents one of the greatest threats to global health. The spread of antibiotic resistance genes among bacteria occurs mostly through horizontal gene transfer via conjugation mediated by plasmids. This process implies a direct contact between a donor and a recipient bacterium which acquires the antibiotic resistance genes encoded by the plasmid and, concomitantly, the capacity to transfer the acquired plasmid to a new recipient. Classical assays for the measurement of plasmid transfer frequency (i.e., conjugation frequency) are often characterized by a high variability and, hence, they require many biological and technical replicates to reduce such variability and the accompanying uncertainty. In addition, classical conjugation assays are commonly tedious and time-consuming because they typically involve counting colonies on a large number of plates for the quantification of donors, recipients, and transconjugants (i.e., the bacteria that have received the genetic material by conjugation). Due to the magnitude of the antibiotic resistance problem, it is critical to develop reliable and rapid methods for the quantification of plasmid transfer frequency that allow the simultaneous analysis of many samples. Here, we present the development of a high-throughput, reliable, quick, easy, and cost-effective method to simultaneously accomplish and measure multiple conjugation events in 96-well plates, in which the quantification of donors, recipients, and transconjugants is estimated from the time required to reach a specific threshold value (OD600 value) in the bacterial growth curves. Our method successfully discriminates different plasmid transfer frequencies, yielding results that are equivalent to those obtained by a classical conjugation assay.

 antibiotic resistance bacterial conjugation conjugation frequency conjugation inhibitors Escherichia coli high-throughput screening platform horizontal gene transfer Ministerio de Ciencia e Innovación10.13039/501100004837 Eusko Jaurlaritza10.13039/501100003086 section-in-acceptance Antibiotic Resistance and New Antimicrobial drugs Introduction

Antibiotics are indispensable tools for the treatment of bacterial infections. Regrettably, the abuse and inappropriate use of antibiotics have caused the emergence and spread of antibiotic resistant bacteria (ARB) that harbor antibiotic resistance genes (ARGs), as a consequence of antibiotic selective pressure, leading to the loss of efficacy of existing antibiotics (Arana et al., 2021). Some ARB show resistance to several or even all known antibiotics, being currently one of the greatest threats to public health. If this worldwide problem continues to grow as projected, it is estimated that by the year 2050 ARB could kill more than 10 million people (O’Neill, 2016).

ARB can transfer ARGs to other bacteria through horizontal gene transfer (HGT) carried out by three main mechanisms: transformation, transduction, and conjugation. Bacterial conjugation, presumably the most relevant mechanism of HGT, is mediated by mobile genetic elements (MGEs), such as conjugative plasmids and integrative conjugative elements (Jechalke et al., 2014). When mediated by plasmids, the conjugation process involves the transfer of a plasmid from a donor to a recipient bacterium, which then acquires new traits (e.g., antibiotic resistance - AR) as well as the ability to transfer the acquired conjugative element to other bacteria, thereby amplifying the risk of AR spread among bacterial populations. The links between environmental bacteria and human pathogenic bacteria are increasingly being recognized as critical to the AR problem (Jauregi et al., 2021a; Jauregi et al., 2021b; Beltrán de Heredia et al., 2023).

Many studies on the bacterial resistome are nowadays focused on the abundance of ARGs and MGE genes in specific settings and under specific conditions (Jauregi et al., 2021c; Beltrán de Heredia et al., 2023), most likely as a consequence of the development of high-throughput (HT) sequencing techniques which have facilitated the realization of metagenomic studies. While this gene-centric approach is undoubtedly generating important scientific knowledge, we must not forget that genetic information reflects only potentiality and that the real risk of AR dissemination lies in the actual transfer of ARGs to human pathogenic bacteria. Then, conjugation assays for the measurement of plasmid transfer frequency (conjugation frequency) between bacterial cells are essential to assess the risk of AR spread among bacterial populations and, specifically, to investigate the transfer of ARGs from environmental bacteria (e.g., bacteria from agro-livestock settings or wastewater treatment plants) to human bacterial pathogens under different conditions and scenarios. Moreover, plasmid transfer assays can be used to determine the capacity of specific compounds to inhibit bacterial conjugation, in an attempt to define new strategies to control AR spread (Álvarez-Rodríguez et al., 2020a).

Conjugation assays performed to study AR spread have traditionally focused on the transfer frequency of ARG-containing conjugative plasmids. In these experiments, a bacterial population harboring an ARG-containing conjugative plasmid (donor) is put into contact with a recipient bacterial population lacking that specific AR under certain conditions of temperature, salinity, pH, nutrients, presence of contaminants or other compounds (e.g., disinfectants, pharmaceuticals), etc. for a given period of time. Subsequently, the recipient bacteria that have acquired the conjugative plasmid (transconjugants) are selected by growing them in the simultaneous presence of the antibiotic for which the conjugative plasmid confers resistance and the antibiotic for which the recipient bacteria are constitutively resistant. In order to quantify plasmid transfer frequency, the ratio of the number of transconjugants per mL to the number of donors or recipients per mL is calculated.

Classical conjugation experiments are often highly tedious and time-consuming because they typically involve counting colonies (colony-forming units, CFU) on a large number of plates for the quantification of donors, recipients, and transconjugants. Normally, it is necessary to plate a series of dilutions on selective media to achieve an easily countable number of colonies (i.e., to ensure a countable plate) for the quantification of donors, recipients, and transconjugants. Classical conjugation assays are frequently characterized by a high variability and, hence, they require many biological and technical replicates to reduce such variability and the accompanying uncertainty (i.e., to have statistically reliable results) (Lorenzo-Díaz and Espinosa, 2009). All this hampers the performance of studies dealing with the effect of different parameters (e.g., environmental conditions) on plasmid transfer frequency and/or the evaluation of the conjugation capacity of a large number of bacterial strains, since the process is time-consuming, arduous, and cumbersome. It is then critical to develop methods for the rapid quantification of plasmid transfer frequency that allow the simultaneous analysis of many samples to be able to promptly assess the risk of AR spread in a large number of samples or test the capacity of many compounds to inhibit conjugation. Several research groups have developed methodologies for a more rapid quantification of conjugation. Some of these methods are based on the modification of the conjugative plasmid so that its transfer can be evaluated by a fluorescent reporter gene incorporated into the plasmid (Fernández-López et al., 2005; Lorenzo-Díaz and Espinosa, 2009), or by insertion in the conjugative plasmid of genes that provide specific metabolic features (Johnsen and Kroer, 2007). Other authors (Wan et al., 2011) opted for the PCR quantification of a specific gene present in the plasmid which, after conjugation, will also be present in the transconjugants. Kneis et al. (2019) used a methodology for the determination of the presence of transconjugants over time and then compared it to the probability predicted by a mathematical model. These strategies allow a reliable measurement of plasmid transfer frequency but, in general, they lack the set of attributes required for an efficient analysis of the risk of AR dissemination in a large number of samples from complex bacterial communities, that is to say, simplicity, speed, high-throughput, low cost, and reproducibility.

The aim of this study was to develop a HT-conjugation platform to perform plasmid transfer assays in a quick, easy, and cost-effective way. To do so, we combined the 96-filter-well plate method reported by Lorenzo-Díaz and Espinosa (2009) with a procedure that allows the measurement of plasmid transfer frequency from a simple readout (time required to reach a given OD600 value in the bacterial growth curve). Plasmid transfer experiments were carried out to search for linear correlations between the number of CFUs of donors, recipients, and transconjugants per mL [log10(CFU mL-1)] and that OD600 readout. The suitability of the HT-conjugation assay was validated through the study of the inhibitory effect of linoleic acid, a compound that affects the transfer of some plasmids but not others.

 Materials and methods Strains and culture media

Three Escherichia coli strains harboring different conjugative plasmids were used as donors in the conjugation assays: E. coli DH5α, E. coli MS411, and E. coli K-12, containing plasmid R388, R1, and pKM101, respectively. Each plasmid belongs to a different incompatibility group (Inc) and confers resistance to the following antibiotic and concentration: R388 to trimethoprim (TMP, 10 µg mL-1), R1 to kanamycin (KAN, 40 µg mL-1), and pKM101 to ampicillin (AMP, 50 µg mL-1) ( Table 1 ). These three strains and plasmids have extensively been used in conjugation studies. Plasmid R1, an F-like conjugative plasmid of the IncFII incompatibility group, has long served as a model system for studying conjugative plasmid biology (Cox and Schildbach, 2017). Plasmid R1, as well as plasmids R388 (a broad-host-range IncW plasmid that contains one of the simplest conjugation systems at the genetic level) and pKM101 (an IncN plasmid), were selected by Fernández-López et al. (2005) for their study on the inhibition of conjugation by unsaturated fatty acids, finding out that conjugation mediated by the first two was inhibited by linoleic acid, unlike that mediated by plasmid pKM101 (see below experiment on inhibition by linoleic acid). Escherichia coli plasmids R388 and pKM101, whose conjugation systems are well-characterized (Christie, 2016), are phylogenetically closely related, but more distantly related to E. coli F plasmid (González-Rivera et al., 2019). Escherichia coli HMS174, a strain with intrinsic chromosomal resistance to rifampicin (RIF, 50 µg mL-1), was used as the only recipient in the conjugation experiments ( Table 1 ).

Escherichia coli strains and plasmids used in the conjugation experiments.

E. coli strain Conjugative plasmid Size (kb) Inc Group Antibiotic resistance# (μg mL-1) Reference
Donors DH5α R388 33.9 IncW 10 TMP Álvarez-Rodríguez et al., 2020b
MS411 R1 97.6 IncFII 40 KAN Cox and Schildbach, 2017
K-12 pKM101 35.4 IncN 50 AMP Langer et al., 1981
Recipient HMS174 * 50 RIF Campbell et al., 1978

*Not applicable.

#For donor strains, antibiotic resistance was conferred by their respective conjugative plasmid. The recipient strain contains intrinsic chromosomal resistance to rifampicin.

Lysogeny Broth (LB) culture medium (Bertani, 1951), supplemented with the corresponding antibiotic, was used for E. coli growth. For the donor bacteria, LB medium was supplemented with 10 µg TMP mL-1 (Merck, Darmstadt, Germany), 40 µg KAN mL-1 (PanReac AppliChem, Darmstadt, Germany), and 50 µg AMP mL-1 (PanReac AppliChem) for E. coli DH5α, E. coli MS411, and E. coli K-12, respectively. For the recipient bacteria, LB medium was supplemented with 50 µg RIF mL-1 (PanReac AppliChem) ( Table 1 ). In both the classical and the HT-conjugation method, quality control tests for donors, recipient, and transconjugants were included, based on the expected antibiotic resistances, to exclude the possibility of bacterial contaminants interfering with the results.

 Classical conjugation assay

In order to measure plasmid transfer frequencies, classical bacterial conjugation assays (C-conjugation) ( Figure 1 ) were performed following Álvarez-Rodríguez et al. (2020b) with some modifications. Briefly, donor and recipient strains were grown overnight in 3.5 mL of LB supplemented with the corresponding antibiotic at 37 °C and orbital shaking at 220 rpm ( Table 1 ). Cultures were then diluted with LB to reach an optical density (OD) value of 0.6 at 600 nm (OD600 = 0.6). Next, 100 μL of the donor culture was mixed with 100 μL of the recipient culture. The mixture was centrifuged at 6,000 x g for 2 min, and the pellet was resuspended in 50 µL of LB without antibiotic. This 50 µL suspension was deposited onto a Millipore filter (0.22 μm pore size) previously placed on a pre-warmed LB-agar plate without antibiotic. After 6 h at 37 °C, cells were detached from the filter by adding 3 mL of 0.8% (w/v) NaCl and vortexing for 30 s. Subsequently, 50 µL of this bacterial suspension was plated on LB-agar with the corresponding selective antibiotic(s) for recipients and transconjugants. Plates were incubated at 37 °C for 24 or 48 h depending on the specific growth rate of each strain. These differences in growth rate are due to the different fitness costs induced by the different conjugative plasmids, with a concomitant effect on bacterial growth rates (Prensky et al., 2021) ( Table 2 ). Plasmid transfer frequency was quantified as the ratio of the number of CFU of transconjugants per mL to the number of CFU of recipients per mL (CFUT mL-1/CFUR mL-1). All assays were performed with five biological replicates and three technical replicates.

Scheme of the plasmid transfer frequency assays. Classical conjugation assay: one single conjugation event can be carried out in each Petri dish. For quantification of plasmid transfer frequencies by colony counting, serial dilutions must be plated in separate LB-agar plates supplemented with the antibiotic used for selection. High-throughput conjugation assay: 96 conjugation events can be carried out simultaneously in a 96-filter-well plate. Quantification of plasmid transfer frequencies can be inferred by the linear correlation between log10(CFU mL-1) values of donors, recipients, and transconjugants and a simple readout, i.e., the time needed to reach a specific OD600 value (time to threshold).

Experimental conditions for donors, recipient, and transconjugants.

Strain Antibiotic for selection(μg mL-1) Growthtime# (h)
Donors E. coli DH5α (R388) 10 TMP 48
E. coli MS411 (R1) 40 KAN 24
E. coli K-12 (pKM101) 50 AMP 24
Recipient E. coli HMS174 50 RIF 24
Transconjugants* TR388 10 TMP+50 RIF 48
TR1 40 KAN+50 RIF 24
TpKM101 50 AMP+50 RIF 24

*TR388, transconjugant that acquired plasmid R388; TR1, transconjugant that acquired plasmid R1; TpKM101, transconjugant that acquired plasmid pKM101.

#Bacteria were grown for 24 or 48 hours depending on their growth rate.

 High-throughput conjugation assay

For the high-throughput bacterial conjugation assays (HT-conjugation) ( Figure 1 ), donor and recipient bacteria were grown overnight in 3.5 mL of LB supplemented with the corresponding antibiotic at 37 °C and orbital shaking at 220 rpm ( Table 1 ). Cultures were then diluted with LB to an OD600 = 0.6. Next, 100 μL of the donor culture was mixed with 100 μL of the recipient culture. The mixture was centrifuged at 6,000 x g for 2 min, and the pellet was resuspended in 50 µL of LB without antibiotic. This 50 µL suspension was deposited in a well of a sterile 96-filter-well plate (0.22 µm) (Merck, Darmstadt). The plate was adapted to a vacuum manifold (MultiScreenHTS) and vacuum (between 5-10 mm Hg) was applied to attach the cells to the filter, following Lorenzo-Díaz and Espinosa (2009). Then, 150 μL of LB without antibiotic was added to the well and the plate was incubated at 37°C for 6 h onto a LB-agar rectangular plate. For those experiments aimed at determining the effect of linoleic acid on plasmid transfer frequencies, the corresponding amounts of linoleic acid were added to the wells to reach the concentrations to be tested. After the 6-hour conjugation time, 150 μL of each well was removed by vacuum filtration as described above. For cell detachment from the filter, 200 μL of 0.8% (w/v) NaCl was added to each well and bacteria were resuspended by pipette mixing 10 times. Next, 10 μL of this suspension was inoculated into a sterile 96-well flat-bottomed microplate (TC-Plate 96 Well, Standard, R, Sarstedt, Nümbrecht, Germany), and 140 μL of LB supplemented with the corresponding antibiotic for the selection of recipients and transconjugants was added ( Table 2 ). Plates were incubated at 37°C for 24 h in an Infinite MNano plate reader (Tecan, Männedorf, Switzerland) and OD600 readings were taken every 10 min to obtain the growth curves for recipients and transconjugants. All assays were performed with five biological replicates and three technical replicates. Plasmid transfer frequency was also quantified as the ratio of the number of CFU of transconjugants per mL to the number of CFU of recipients per mL (CFUT mL-1/CFUR mL-1). In this case, instead of counting colonies on LB-agar plates, CFUT mL-1 and CFUR mL-1 values were estimated from the abovementioned growth curves.

Linear correlations between the time needed to reach a given OD600 value in the growth curves of recipients and transconjugants and the number of CFU mL-1 of recipients and transconjugants (on a logarithmic scale) ( Figures 2B ; 3D–F ) were searched for. To obtain the growth curves, serial dilutions of the recipient and transconjugant suspensions described above were prepared with LB medium ( Figures 2A ; 3A–C ). Depending on the growth rate of each strain, dilutions ranged from 10-1 to 10-5 for recipient and from 10-1 to 10-6 for transconjugant suspensions. Dilutions for plasmid R388-containing transconjugants (TR388) ranged from 10-1 to 10-3; dilutions for plasmid R1-containing transconjugants (TR1) ranged from 10-1 to 10-5; and dilutions for plasmid pKM101-containing transconjugants (TpKM101) ranged from 10-1 to 10-4. Growth curves were obtained in an Infinite MNano plate reader as described above. Dilutions of recipients and transconjugants were plated on LB-agar medium supplemented with the corresponding selective antibiotics ( Table 2 ) to determine the number of CFU mL-1. When required to ensure a countable plate, suspensions were further diluted to the appropriate concentration at which CFU mL-1 could be easily counted in the LB-agar plates.

Correlation between log10(CFU mL-1) values and the time needed to reach an OD600 = 0.3 (time to threshold) for recipient E. coli HMS174 strain. (A) Growth curves of the different bacterial dilutions. The dotted line marks the threshold value. (B) Linear regression between log10(CFU mL-1) values obtained by colony counting in LB-agar plates supplemented with 50 μg RIF mL-1 and the time at which each growth curve reaches an OD600 = 0.3 (time to threshold). The solid line represents the best data fit. Dashed lines correspond to the margin of error of the data with a 95% confidence interval.

Correlation between log10 (CFU mL-1) values and the time needed to reach an OD600 = 0.3 (time to threshold) for transconjugants. (A) TR388; (B) TR1; (C) TpKM101) Growth curves of the different bacterial dilutions. The dotted line marks the threshold value. (D) TR388; (E) TR1; (F) TpKM101) Linear regression between log10 (CFU mL-1) values obtained by colony counting in LB-agar plates supplemented with the appropriate antibiotic for transconjugant selection and the time at which each growth curve reaches an OD600 = 0.3 (time to threshold). The solid line represents the best data fit. Dashed lines correspond to the margin of error of the data with a 95% confidence interval.

Calibration curves that correlate the time required to reach a predetermined threshold value (“time to threshold”) with inoculum size can be constructed, so that unknown initial cell densities can be predicted from time to threshold values (Bethke et al., 2020). In our study, it was observed that, at an OD600 = 0.3, all cultures had reached the exponential phase ( Figures 2A ; 3A–C ). Then, the CFUT mL-1 and CFUR mL-1 values obtained in the LB-agar plates inoculated with the abovementioned dilutions were represented against the time (in hours) needed to reach that threshold value (OD600 = 0.3) in the growth curves, i.e., “time to threshold” ( Figures 2 ; 3D–F ). The resulting plots can be used as calibration plots to calculate the number of recipients and transconjugants of unknown samples according to the time required to reach the threshold value.

The suitability of the HT-conjugation assay was validated through the study of the inhibitory effect of linoleic acid, a compound that has been reported to affect the transfer of some conjugative plasmids but not others (Fernández-López et al., 2005). Using the HT-conjugation assay, plasmid transfer frequencies were measured in the presence of different concentrations of linoleic acid (0, 1, 10, 50, 100, 200, 300, and 500 µM). Not all of these concentrations were tested for all cases (concentrations of linoleic acid were adjusted depending on the observed degree of inhibition for each case). According to preliminary studies (data not shown), none of these linoleic acid concentrations affected the growth of either the donor or recipient cells.

 Graphical representations and statistical analysis

Bacterial growth curves ( Figures 2A ; 3A–C ) and calibration plots ( Figures 2B ; 3D–F ) were graphed using GraphPad Prism 9 (GraphPad Software, San Diego, CA). Statistical analyses, carried out to search for statistically significant differences in plasmid transfer frequencies between C-conjugation and HT-conjugation assays, were performed using the two-tailed paired Student’s t-test using GraphPad Prism 9 (GraphPad Software), where p < 0.05 was considered to be significant.

 Results Development of the high-throughput platform

In conjugation assays, plasmid transfer frequency is often quantified as the ratio of the number of CFU of transconjugants per mL to the number of CFU of recipients (or donors) per mL, which usually requires counting colonies on a large number of plates. Here, instead of counting colonies on selective medium plates, CFUT mL-1 and CFUR mL-1 values were estimated from the growth curves of recipient and transconjugant bacteria. Such estimations were based on linear correlations between the number of CFU mL-1 (on a logarithmic scale) and the time needed to reach a specific OD600 value in the growth curves ( Figures 2 , 3 ). Initially, serial dilutions of recipient and transconjugant suspensions were prepared and plated on LB-agar medium supplemented with the corresponding selective antibiotics to determine the number of CFU mL-1 of recipient and transconjugants ( Table 2 ). These recipient and transconjugant suspensions were used to obtain the growth curves (OD600). As shown in Figures 2 (2A for recipient) and 3 (3A–C for transconjugants), at an OD600 = 0.3, all cultures had reached the exponential phase. Then, the time required to reach an OD600 = 0.3 was set as threshold value. As shown in Figures 2 (2B for recipient) and 3 (3D–F for transconjugants), very good inverse linear correlations between log10(CFU mL-1) values and time to threshold values were detected: r2 values ranged from 0.95 to 0.99. Maximum OD600 values reached by the different bacterial cultures were around 0.75-0.80 and the duration of the lag phase ranged from close to zero for larger inoculum sizes to 10-12 h for smaller inoculum sizes. The good regression parameters obtained for all cases indicate that the calibration plots can be used to calculate the number of recipients ( Figure 2B ) and transconjugants ( Figures 3D–F ) from “time to threshold” values.

The HT-conjugation assay was validated by comparison with the results obtained with the classical conjugation assay (C-conjugation). To this purpose, once the HT-conjugation assay had been designed and the calibration plots [based on linear regressions between log10(CFU mL-1) values and time to threshold values] for recipients and transconjugants had been established, it was evaluated whether the plasmid transfer frequencies obtained with the HT-conjugation assay were comparable to those obtained with the C-conjugation assay. To this goal, transfer frequencies of plasmids R388, R1, and pKM101 were measured by both assays, finding out similar values for the three plasmids. Conjugation frequency values for plasmid R388 obtained with the C-conjugation and HT-conjugation assays were 2.27 x 10-4 ± 4.12 x 10-5 and 4.69 x 10-4 ± 5.22 x 10-5, respectively. For plasmid R1 conjugation frequency, the values obtained by both assays were 4.49 x 10-3 ± 5.95 x 10-4 (C-conjugation) and 1.59 x 10-3 ± 2.15 x 10-4 (HT-conjugation), and for plasmid pKM101 were 2.03 x 10-1 ± 3.05 x 10-2 (C-conjugation) and 1.04 x 10-1 ± 3.35 x 10-2 (HT-conjugation) ( Figure 4 ). According to the two-tailed paired Student’s t-test, the observed differences between the values obtained from both assays for the three plasmids were not statistically significant (differences were in the range of the variability observed within each assay). In other words, the observed differences in transfer frequencies of the three plasmids between C-conjugation and HT-conjugation assays were not biologically relevant, pointing out to the suitability of the HT-conjugation assay presented here.

Plasmid transfer frequencies of plasmids R388 (A), R1 (B), and pKM101 (C) measured by the classical conjugation (C-conjugation) versus the high-throughput conjugation (HT-conjugation) assay. In the C-conjugation assay, plasmid transfer frequencies were obtained by colony counting (CFU mL-1) in LB-agar plates. In the HT-conjugation assay, plasmid transfer frequencies were obtained according to the calibration plots used to calculate the number of recipients and transconjugants (CFU mL-1) inferred from the time to threshold value. Plasmid transfer frequency was defined as the ratio between the number of CFU mL-1 of transconjugants and the number of CFU mL-1 of recipients (CFUT mL-1/CFUR mL-1). Results are the mean ± SE of five independent experiments performed in triplicate. Differences between the two methods were not statistically significant (p > 0.05).

 Validation of the high-throughput platform with linoleic acid as inhibitor

Once it had been verified that the HT-conjugation assay could quantify the transfer frequency of the three plasmids, its capacity to differentiate expected differences in conjugation frequencies, according to previously reported knowledge, was tested. To do so, a compound that has been reported to inhibit the transfer of some conjugative plasmids but not others, i.e., linoleic acid (Fernández-López et al., 2005), was selected. Linoleic acid does inhibit conjugation mediated by plasmid R388 and F plasmids (such as plasmid R1), but not by plasmids pKM101, R6K, and RP4 (Fernández-López et al., 2005). R388, R1, and pKM101 plasmid transfer frequencies were measured in the presence of different concentrations of linoleic acid using the HT-conjugation assay, finding out that it did detect the inhibitory effect of linoleic acid on the conjugation frequency of plasmids R388 and R1 ( Figure 5 ). Furthermore, the transfer frequency of plasmid pKM101 was not inhibited by linoleic acid ( Figure 5 ). The values of plasmid transfer frequencies for each plasmid in the presence of different concentrations of linoleic acid are shown in Table 3 . The transfer frequency of plasmid R388 decreased from 5.84 x 10-6 ± 2.94 x 10-7 in the absence of linoleic acid to 2.72 x 10-6 ± 4.44 x 10-8 in the presence of 100 µM linoleic acid; thus, linoleic acid inhibited the transfer of plasmid R388 by approximately 50% (statistical significance: ***p < 0.005). For plasmid R1, the frequency of conjugation decreased from 1.81 x 10-3 ± 9.07 x 10-5 in the absence of linoleic acid to 4.66 x 10-4 ± 4.94 x 10-5, 2.41 x 10-4 ± 2.28 x 10-5, 1.03 x 10-4 ± 7.86 x 10-6, and 1.67 x 10-4 ± 1.07 x 10-5 in the presence of 50, 100, 200, and 300 µM linoleic acid, respectively. All the differences were statistically significant. Then, an inhibition of almost 75% was observed at 50 µM linoleic acid and above. Finally, plasmid pKM101 conjugation frequencies measured with the HT-conjugation assay showed no statistically significant differences in the presence versus absence of linoleic acid, reflecting the fact that this compound does not inhibit the transfer of plasmid pKM101, in agreement with Fernández-López et al. (2005). It was concluded that the HT-conjugation assay can differentiate expected differences (according to previously reported knowledge) in conjugation frequencies.

Effect of linoleic acid on plasmid transfer frequency measured by the high-throughput conjugation assay. Data are expressed as percentage of inhibition relative to control value (0% = no linoleic acid). Results are the mean ± SE of five independent experiments performed in triplicate. The statistical significance of the observed differences with respect to the control (no linoleic acid) was analyzed by Student’s t test (***p <0.005, ****p < 0.001).

Effect of different concentrations of linoleic acid on conjugation frequencies of plasmids R388, R1, and pKM101.

Linoleic acid(µM) Conjugation frequency(CFUT mL-1/CFUR mL-1)
R388 0 5.84 x 10-6 ± 2.94 x 10-7
1 6.42 x 10-6 ± 6.99 x 10-7
10 5.94 x 10-6 ± 1.79 x 10-7
100 2.72 x 10-6 ± 4.44 x 10-8
R1 0 1.81 x 10-3 ± 9.07 x 10-5
1 1.27 x 10-3 ± 9.49 x 10-5
50 4.66 x 10-4 ± 4.94 x 10-5
100 2.41 x 10-4 ± 2.28 x 10-5
200 1.03 x 10-4 ± 7.86 x 10-6
300 1.67 x 10-4 ± 1.07 x 10-5
pKM101 0 1.42 x 10-2 ± 8.01 x 10-4
200 1.22 x 10-2 ± 6.36 x 10-4
500 1.68 x 10-2 ± 6.81 x 10-4

Mean values ± SE of five independent experiments performed in triplicate.

 Discussion

Antibiotic resistance phenotypes rely on the presence of ARGs in bacterial chromosomes and/or extrachromosomal genetic elements. ARG-harboring bacteria can transfer ARGs to other bacteria by HGT mediated by MGEs. Then, we must closely monitor the possibility of ARG transfer from non-pathogenic environmental bacteria to human pathogenic bacteria (Perry and Wright, 2013; von Wintersdorff et al., 2016). Conjugation, the main mechanism responsible for AR dissemination (Frost et al., 2005), requires a physical contact between the donor cells carrying the plasmid and the recipient cells to then generate the transconjugants (Ochman et al., 2000), which can later contribute to plasmid transfer to new recipients. When conjugative plasmids enter a bacterium, they often trigger its SOS response, driving the swapping of integron cassettes and stimulating the generation of multiresistant strains (Baharoglu et al., 2010). Moreover, antibiotic treatment has been suggested to promote HGT, but this fact is still a matter of much debate (Møller et al., 2017).

Due to its crucial role in AR spread, conjugation has been extensively studied from different points of view and using a variety of experimental approaches (Huisman et al., 2022). The quantification of plasmid transfer frequency has been widely applied to assess the risk of AR spread (González-Plaza et al., 2019; Macedo et al., 2022). Nonetheless, there is no generally accepted method to quantify conjugation frequency (Huisman et al., 2022). Being aware of this limitation, Huisman et al. (2022) compared a variety of methods, finding out that conjugation frequency estimates strongly depend on time of measurement, initial population densities, and the ratio of donor to recipient populations (differences in growth rates also affected conjugation frequency estimates). The quantification of plasmid transfer has mostly focused on a few laboratory strains, which might be poor models for natural bacterial populations (Hobman et al., 2007; Dimitriu et al., 2019). For instance, E. coli K-12 has been questioned as a model organism since the current versions of this strain appear to be quite distinct from any original, wild-type, free-living ancestor and, besides, it is usually studied under conditions far removed from its natural habitats (Hobman et al., 2007).

In order to carry out experimental work on a large number of environmental and/or clinical samples, it is crucial to have a method that can accurately, quickly, and cost-effectively quantify plasmid transfer frequency. Classical laboratory assays for the determination of conjugation frequency are commonly highly tedious and time-consuming. Besides, they are poorly reproducible and, hence, it is often necessary to increase the number of replicates for statistical robustness. In consequence, the high number of biological and technical replicates that are required to provide statistical robustness substantially complicates the development of experiments involving a large number of samples and/or experimental conditions (Kneis et al., 2019). In the case of environmental samples, methodological strategies based on the modification of the conjugative plasmid to assess its transfer (Fernández-López et al., 2005; Johnsen and Kroer, 2007; Lorenzo-Díaz and Espinosa, 2009) are often not feasible, as the very nature of the samples and their diversity make these approaches a daunting task. Within this context, Alalam et al. (2020) designed a high-throughput platform for the measurement of AR plasmid conjugation and the identification of all the E. coli genes that control F-plasmid donation, as well as novel conjugation genes not previously linked to AR transmission.

Inspired by the work of Lorenzo-Díaz and Espinosa (2009), conjugation assays were carried out here in 96-filter-well plates. Filter mating provides high cell density and close proximity for donor and recipient cells, facilitating the formation of mating pairs for plasmid transfer (Lim et al., 2008; Aviv et al., 2016). After conjugation, instead of quantifying the number of recipients and transconjugants by plating a series of dilutions on LB-agar plates supplemented with the corresponding selective antibiotics, sterile 96-well flat-bottomed microplates were used to record the growth curves of the diluted cell suspensions in the appropriate selective media. The quantification of recipients and transconjugants was carried out by defining a “time required to reach a predetermined threshold value of an easily measurable parameter” (OD600), following Bethke et al. (2020). From the growth curves, it was observed that at an OD600 = 0.3 all cultures had reached the exponential phase, and the time taken for each culture to reach this OD600 value was established as “time to threshold”. A linear negative relationship between “time to threshold” and inoculum size of recipients and transconjugants, measured as number of CFU per mL of inoculum in LB-agar plates [expressed as log10(CFU mL-1)], was found. This negative linear correlation between the log of the initial inoculum size [expressed as log10(CFU mL-1)] and the “time to threshold” value, observed for the range of serial dilutions used here, is mostly, if not entirely, due to the decrease in the initial number of cells growing exponentially (i.e., the increase in the “time to threshold” value observed in the diluted cultures is a direct consequence of the lower initial concentration of cells). The negative linear correlation observed in exponentially growing cells between the log of the initial inoculum size [log10(CFU mL-1)] and the “time to threshold” value was then used to determine conjugation frequencies. The HT conjugation procedure presented here is valid as long as the experimental values of the “time to threshold” are within the linear range.

For donors, a linear relationship was also observed (data not shown) but, since the quantification of plasmid transfer frequencies was referred here to a common recipient strain (E. coli HMS174), only the calibration plots for recipients and transconjugants are shown. The linear regression parameters of all the calibration plots were good, with r2 values ranging from 0.95 to 0.99. All the experimental points of the calibration plots were within a margin of error below the 95% confidence interval. The calibration plots were used to infer the number of recipients and transconjugants from the “time to threshold” values obtained from the growth curves, thus allowing the simultaneous assessment of 96 conjugation events.

Donor and recipient cells were put into contact for 6 h (conjugation time = 6 h). When comparing data among conjugation studies, the duration of the conjugation time must be taken into consideration, as conjugation frequency curves often resemble bacterial growth curves, with a lag phase occurring after the initial mixing of donor and recipient cells, followed by a period of increasing conjugation that normally ends in a plateau (Andrup and Andersen, 1999; Fox et al., 2008; Wan et al., 2011; Headd and Bradford, 2020). Although most studies on conjugation frequency are focused on cultures of bacteria growing for short periods of time, under natural conditions bacteria often grow episodically and at much lower rates than those observed in the laboratory. Headd and Bradford (2020) studied the dynamics of conjugation of E. coli cells harboring plasmid pUUH239.2 (an IncFII plasmid containing multiantibiotic- and metal-resistance genes) over a 24-h period in the presence and absence of antibiotics, finding out that (i) conjugation occurred only in a narrow time frame when E. coli cells were transitioning from a growing to a non-growing phase; and (ii) the conjugation plateau developed because of a lack of capable donor cells. Their results suggest that episodic growth in nutrient-depleted environments (as it is frequently the case in many environmental settings) could result in more conjugation than sustained growth in a nutrient rich environment (Headd and Bradford, 2020). Conjugation dynamics are basically governed by two kinetic processes: the rate of gene transfer (or conjugation efficiency) and the relative growth rate of transconjugants (or the growth dynamics) (Prensky et al., 2021). The hours immediately following conjugation can play a critical role in both short- and long-term plasmid prevalence; actually, newly generated transconjugants can exhibit transient reduced growth rates and/or prolonged lag times, pointing out to the presence of an increased burden immediately after conjugation (plasmid acquisition cost) which occurs independently of long-term fitness effects (Prensky et al., 2021). It has been reported (Newbury et al., 2022) that the degree to which a plasmid spreads through a host community (both the number of host species in which it is found and the evenness of its prevalence between host species) can be explained as a direct result of the effect the plasmid has on its hosts’ population growth rate. In any event, many more studies on the dynamics of bacterial conjugation under conditions that mimic real settings are needed to properly assess the risk of AR dissemination.

Huisman et al. (2022) emphasized that conjugation methods, based on the ratio between population densities, e.g., transconjugants/recipients or transconjugants/donors, vary as a function of initial population densities, initial donor to recipient ratios, and the length of the conjugation assay (Simonsen et al., 1990; Zhong et al., 2012). In a study on plasmid transfer among E. coli isolates harboring one to five plasmids, Darphorn et al. (2022) found that, under standardized conditions, plasmid transfer rates varied by a factor of almost 106, depending on the specific recipient/donor strain combination. In consequence, to allow for comparisons among studies, comprehensive information on the experimental conditions and strain growth rates must always be included (Huisman et al., 2022). On the other hand, conjugation efficiency is influenced by the taxonomic relatedness between donor and recipient bacteria (Hunter et al., 2008). Dimitriu et al. (2019) studied the importance of genetic similarity between naturally co-occurring E. coli isolates for plasmid transfer, finding out that plasmid transfer was strongly biased toward clone-mates, but not correlated to genetic distance when donors and recipients were not clone-mates. Alderliesten et al. (2020) carried out a meta-analysis on conjugation frequencies from E. coli donors to various recipient species and observed that taxonomic relatedness was limiting conjugation in liquid matings, but not in filter matings, suggesting that taxonomic relatedness might not be a limiting factor for conjugation in those environments where bacteria are fixed in space.

Although a wide variety of biotic and abiotic factors can affect plasmid transfer (e.g., temperature, nutrient concentration, pH, moisture, cell physiology, population densities, type of plasmid, type of donor and recipient, growth on surfaces versus growth in well-mixed liquids), the transfer of conjugative plasmids in biofilms is of paramount importance (Águila-Arcos et al., 2017), as most bacteria live in biofilms, microcolonies or similar forms of clumped cells with an explicit spatial structure (Stalder and Top, 2016). Biofilm growth can promote the initial transfer of some plasmids, but also limit further plasmid invasion into the population or community (Stalder and Top, 2016). A much less explored factor that can also affect plasmid transfer is the ecological growth strategy (K-strategists, r-strategists, oligotrophs, copiotrophs) of donors and recipients (Varner et al., 2022). In any case, the monitoring of environmental HGT remains challenging mainly due to cultivation bias (Pallares-Vega et al., 2021).

The high-throughput platform presented here has been validated according to two different approaches. First, the transfer frequencies of three plasmids belonging to different incompatibility groups (R388, R1, pKM101) were quantified by both the classical (C-conjugation) and the high-throughput (HT-conjugation) assay. In all cases, similar values were obtained by the C-conjugation versus the HT-conjugation assay. Due to the variability of the conjugation process itself, the observed differences between both methods were smaller than those found in different experiments carried out with the same method. Secondly, it was tested whether the HT-conjugation assay could differentiate expected differences in conjugation frequencies. For this purpose, the effect of linoleic acid, a compound that inhibits the transfer of some conjugative plasmids but not others (Fernández-López et al., 2005), was determined. When HT-conjugation experiments were performed with plasmids R388 and R1, for which the inhibition of conjugation by linoleic acid has been reported (Fernández-López et al., 2005), the expected inhibition was observed. When the frequency of conjugation was measured with plasmid pMK101, linoleic acid showed no effect on such frequency, as previously reported (Fernández-López et al., 2005).

From all of the above, it was concluded that the high-throughput conjugation platform developed here is a suitable, reliable, time and money saving, methodological tool for the simultaneous assessment of plasmid transfer frequency in a wide number of environmental and/or clinical samples. This high-throughput platform can, for instance, be used to simultaneously screen a large number of molecules or compounds as potential inhibitors of bacterial conjugation, in an attempt to minimize the risk of AR dissemination. Despite these promising results, it must be taken into consideration that this study has been carried out using only three E. coli strains (each harboring a different plasmid) as donor cells and one E. coli strain as recipient. For a more comprehensive validation, the HT-conjugation platform needs to be further put to the test with other strains of the same or different bacterial species, other types of conjugative plasmids, and other experimental conditions, a requirement which will be fulfilled as it is being used in other conjugation studies.

 Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

 Author contributions

KA-C: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. AR-S: Data curation, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing – review & editing. NS-F: Methodology, Validation, Writing – review & editing. OA: Methodology, Validation, Writing – review & editing. CG: Conceptualization, Investigation, Supervision, Validation, Writing – original draft, Writing – review & editing. LA: Conceptualization, Formal Analysis, Investigation, Supervision, Writing – original draft, Writing – review & editing. IA: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.

 Funding

The authors declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by MCIN/AEI/10.13039/501100011033 (PID2020-116495RB-I00) and the Basque Government (IT1578-22). AR-S received an UPV/EHU Training of pre-doctoral Research Staff Grant. NS-F received an Ikasiker grant from the Basque Government.

 Acknowledgments

The authors are grateful to Drs. Manuel Espinosa, Alicia Bravo, Sofia Ruiz-Cruz, and Ana Moreno-Blanco for their pertinent scientific comments. The authors wish to thank support from Euskampus – JRL Environmental Antibiotic Resistance.

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 References Águila-Arcos S. Álvarez-Rodríguez I. Garaiyurrebaso O. Garbisu C. Grohmann E. Alkorta I. (2017). Biofilm-forming clinical Staphylococcus isolates harbor horizontal transfer and antibiotic resistance genes. Front. Microbiol. 8. doi: 10.3389/fmicb.2017.02018 Alalam H. Graf F. E. Palm M. Abadikhah M. Zackrisson M. Boström J. . (2020). A highthroughput method for screening for genes controlling bacterial conjugation of antibiotic resistance. mSystems 5, e01226-20. doi: 10.1128/mSystems.01226-20 Alderliesten J. B. Duxbury S. J. Zwart M. P. De Visser J. Stegeman A. Fischer E. A. (2020). Effect of donor-recipient relatedness on the plasmid conjugation frequency: a meta-analysis. BMC Microbiol. 20, 1–10. doi: 10.1186/s12866-020-01825-4 Álvarez-Rodríguez I. Arana L. Ugarte-Uribe B. Gómez-Rubio E. Martín-Santamaría S. Garbisu C. . (2020a). Type IV coupling proteins as potential targets to control the dissemination of antibiotic resistance. Front. Mol. Biosci. 7. doi: 10.3389/fmolb.2020.00201 Álvarez-Rodríguez I. Ugarte-Uribe B. de la Arada I. Arrondo J. L. R. Garbisu C. Alkorta I. (2020b). Conjugative coupling proteins and the role of their domains in conjugation, secondary structure and in vivo subcellular location. Front. Mol. Biosci. 185. doi: 10.3389/fmolb.2020.00185 Andrup L. Andersen K. (1999). A comparison of the kinetics of plasmid transfer in the conjugation systems encoded by the F plasmid from Escherichia coli and plasmid pCF10 from Enterococcus faecalis . Microbiology 145, 2001–2009. doi: 10.1099/13500872-145-8-2001 Arana L. Gallego L. Alkorta I. (2021). Incorporation of antibiotics into solid lipid nanoparticles: a promising approach to reduce antibiotic resistance emergence. Nanomaterials 11, 1251. doi: 10.3390/nano11051251 Aviv G. Rahav G. Gal-Mor O. (2016). Horizontal transfer of the Salmonella enterica serovar infantis resistance and virulence plasmid pESI to the gut microbiota of warm-blooded hosts. mBio 7, e01395–e01316. doi: 10.1128/mBio.01395-16 Baharoglu Z. Bikard D. Mazel D. (2010). Conjugative DNA transfer induces the bacterial SOS response and promotes antibiotic resistance development through integron activation. PloS Genet. 6, e1001165. doi: 10.1371/journal.pgen.1001165 Beltrán de Heredia I. Garbisu C. Alkorta I. Urra J. González-Gaya B. Ruiz-Romera E. (2023). Spatio-seasonal patterns of the impact of wastewater treatment plant effluents on antibiotic resistance in river sediments. Environ. Pollut. 319, 120883. doi: 10.1016/j.envpol.2022.120883 Bertani G. (1951). Studies on lysogenesis. I. The mode of phage liberation by lysogenic Escherichia coli . J. Bacteriol. 62, 293–300. doi: 10.1128/jb.62.3.293-300.1951 Bethke J. H. Davidovich A. Cheng L. Lopatkin A. J. Song W. Thaden J. T. . (2020). Environmental and genetic determinants of plasmid mobility in pathogenic Escherichia coli . Sci. Adv. 6, eaax3173. doi: 10.1126/sciadv.aax3173 Campbell J. L. Richardson C. C. Studier F. W. (1978). Genetic recombination and complementation between bacteriophage T7 and cloned fragments of T7 DNA. Proc. Natl. Acad. Sci. 75, 2276–2280. doi: 10.1073/pnas.75.5.2276 Christie P. J. (2016). The mosaic type IV secretion systems. EcoSal. Plus 7 (1). doi: 10.1128/ecosalplus.ESP-0020-2015 Cox K. E. L. Schildbach J. F. (2017). Sequence of the R1 plasmid and comparison to F and R100. Plasmid 91, 53–60. doi: 10.1016/j.plasmid.2017.03.007 Darphorn T. S. Koenders-Van Sintannel 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, e0270205. doi: 10.1371/journal.pone.0270205 Dimitriu T. Marchant L. Buckling A. Raymond B. (2019). Bacteria from natural populations transfer plasmids mostly towards their kin. Proc. R. Soc. B. 286, 20191110. doi: 10.1098/rspb.2019.1110 Fernández-López R. Machón C. Longshaw C. M. Martin S. Molin S. Zechner E. L. . (2005). Unsaturated fatty acids are inhibitors of bacterial conjugation. Microbiology 151, 3517–3526. doi: 10.1099/mic.0.28216-0 Fox R. E. Zhong X. Krone S. M. Top E. M. (2008). Spatial structure and nutrients promote invasion of IncP-1 plasmids in bacterial populations. ISME J. 2, 1024–1039. doi: 10.1038/ismej.2008.53 Frost L. S. Leplae R. Summers A. O. Toussaint A. (2005). Mobile genetic elements: the agents of open source evolution. Nat. Rev. Microbiol. 3, 722–732. doi: 10.1038/nrmicro1235 González-Plaza J. J. Blau K. Milaković M. Jurina T. Smalla K. Udiković-Kolić N. (2019). Antibiotic-manufacturing sites are hot-spots for the release and spread of antibiotic resistance genes and mobile genetic elements in receiving aquatic environments. Environ. Int. 130, 104735. doi: 10.1016/j.envint.2019.04.007 González-Rivera C. Khara P. Awad D. Patel R. Li Y. G. Bogisch M. . (2019). Two pKM101-encoded proteins, the pilus-tip protein TraC and Pep, assemble on the Escherichia coli cell surface as adhesins required for efficient conjugative DNA transfer. Mol. Microbiol. 111, 96–117. doi: 10.1111/mmi.14141 Headd B. Bradford S. A. (2020). The conjugation window in an Escherichia coli K-12 strain with an IncFII plasmid. Appl. Environ. Microbiol. 86, e00948–e00920. doi: 10.1128/AEM.00948-20 Hobman J. L. Penn C. W. Pallen M. J. (2007). Laboratory strains of Escherichia coli: model citizens or deceitful delinquents growing old disgracefully? Mol. Microbiol. 64, 881–885. doi: 10.1111/j.1365-2958.2007.05710.x Huisman J. S. Benz F. Duxbury S. J. De Visser J. A. G. Hall A. R. Fischer E. A. . (2022). Estimating plasmid conjugation rates: a new computational tool and a critical comparison of methods. Plasmid. 121, 102627. doi: 10.1016/j.plasmid.2022.102627 Hunter P. R. Wilkinson D. C. Catling L. A. Barker G. C. (2008). Meta-analysis of experimental data concerning antimicrobial resistance gene transfer rates during conjugation. Appl. Environ. Microbiol. 74, 6085–6090. doi: 10.1128/AEM.01036-08 Jauregi L. Epelde L. Alkorta I. Garbisu C. (2021a). Agricultural soils amended with thermally-dried anaerobically-digested sewage sludge showed increased risk of antibiotic resistance dissemination. Front. Microbiol. 12, 666854. doi: 10.3389/fmicb.2021.666854 Jauregi L. Epelde L. Alkorta I. Garbisu C. (2021b). Antibiotic resistance in agricultural soil and crops associated to the application of cow manure-derived amendments from conventional and organic livestock farms. Front. Vet. Sci. 8, 633858. doi: 10.3389/fvets.2021.633858 Jauregi L. Epelde L. González A. Lavín J. L. Garbisu C. (2021c). Reduction of the resistome risk from cow slurry and manure microbiomes to soil and vegetable microbiomes. Environ. Microbiol. 23, 7643–7660. doi: 10.1111/1462-2920.15842 Jechalke S. Schreiter S. Wolters B. Dealtry S. Heuer H. Smalla K. (2014). Widespread dissemination of class 1 integron components in soils and related ecosystems as revealed by cultivation-independent analysis. Front. Microbiol. 4. doi: 10.3389/fmicb.2013.00420 Johnsen A. R. Kroer N. (2007). Effects of stress and other environmental factors on horizontal plasmid transfer assessed by direct quantification of discrete transfer events. FEMS Microbiol. Ecol. 59, 718–728. doi: 10.1111/j.1574-6941.2006.00230.x Kneis D. Hiltunen T. Heß S. (2019). A high-throughput approach to the culture-based estimation of plasmid transfer rates. Plasmid 101, 28–34. doi: 10.1016/j.plasmid.2018.12.003 Langer P. J. Shanabruch W. G. Walker G. C. (1981). Functional organization of plasmid pKM101. J. Bacteriol. 145, 1310–1316. doi: 10.1128/jb.145.3.1310-1316.1981 Lim Y. M. Groof A. J. C. D. Bhattacharjee M. K. Figurski D. H. Schon E. A. (2008). Bacterial conjugation in the cytoplasm of mouse cells. Infect. Immun. 76, 5110–5119. doi: 10.1128/IAI.00445-08 Lorenzo-Díaz F. Espinosa M. (2009). Large-scale filter mating assay for intra- and inter-specific conjugal transfer of the promiscuous plasmid pMV158 in Gram-positive bacteria. Plasmid 61, 65–70. doi: 10.1016/j.plasmid.2008.09.005 Macedo G. Olesen A. K. Maccario L. Hernandez Leal L. van der Maas P. Heederik D. . (2022). Horizontal gene transfer of an IncP1 plasmid to soil bacterial community introduced by Escherichia coli through manure amendment in soil microcosms. Environ. Sci. Technol. 56, 11398–11408. doi: 10.1021/acs.est.2c02686 Møller T. S. B. Liu G. Boysen A. Thomsen L. E. Lüthje F. L. Mortensen S. . (2017). Treatment with cefotaxime affects expression of conjugation associated proteins and conjugation transfer frequency of an IncI1 plasmid in Escherichia coli . Front. Microbiol. 8. doi: 10.3389/fmicb.2017.02365 Newbury A. Dawson B. Klümper U. Hesse E. Castledine M. Fontaine C. . (2022). Fitness effects of plasmids shape the structure of bacteria–plasmid interaction networks. Proc. Natl. Acad. Sci. 119, e2118361119. doi: 10.1073/pnas.2118361119 Ochman H. Lawrence J. G. Groisman E. A. (2000). Lateral gene transfer and the nature of bacterial innovation. Nature 405, 299–304. doi: 10.1038/35012500 O’Neill J. (2016). Review on antimicrobial resistance: tackling drug-resistant infections globally: final report and recommendations. United Kingdom: Wellcome Trust and HM Government. Pallares-Vega R. Macedo G. Brouwer M. S. Hernandez Leal L. van der Maas P. Van Loosdrecht M. C. . (2021). Temperature and nutrient limitations decrease transfer of conjugative IncP-1 plasmid pKJK5 to wild Escherichia coli strains. Front. Microbiol. 12, 656250. doi: 10.3389/fmicb.2021.656250 Perry J. A. Wright G. D. (2013). The antibiotic resistance “mobilome”: searching for the link between environment and clinic. Front. Microbiol. 4, 138. doi: 10.3389/fmicb.2013.00138 Prensky H. Gomez-Simmonds A. Uhlemann A. C. Lopatkin A. J. (2021). Conjugation dynamics depend on both the plasmid acquisition cost and the fitness cost. Mol. Syst. Biol. 17, e9913. doi: 10.15252/msb.20209913 Simonsen L. Gordon D. M. Stewart F. M. Levin B. R. (1990). Estimating the rate of plasmid transfer: an end-point method. Microbiology 136, 2319–2325. doi: 10.1099/00221287-136-11-2319 Stalder T. Top E. (2016). Plasmid transfer in biofilms: a perspective on limitations and opportunities. NPJ Biofilms Microbiomes 2, 1–5. doi: 10.1038/npjbiofilms.2016.22 Varner P. M. Allemann M. N. Michener J. K. Gunsch C. K. (2022). The effect of bacterial growth strategies on plasmid transfer and naphthalene degradation for bioremediation. Environ. Technol. Innovation 28, 102910. doi: 10.1016/j.eti.2022.102910 von Wintersdorff C. J. Penders J. Van Niekerk J. M. Mills N. D. Majumder S. Van Alphen L. B. . (2016). Dissemination of antimicrobial resistance in microbial ecosystems through horizontal gene transfer. Front. Microbiol. 7, 173. doi: 10.3389/fmicb.2016.00173 Wan Z. Varshavsky J. Teegala S. Mclawrence J. Goddard N. L. (2011). Measuring the rate of conjugal plasmid transfer in a bacterial population using quantitative PCR. Biophys. J. 101, 237–244. doi: 10.1016/j.bpj.2011.04.054 Zhong X. Droesch J. Fox R. Top E. M. Krone S. M. (2012). On the meaning and estimation of plasmid transfer rates for surface-associated and well-mixed bacterial populations. J. Theor. Biol. 294, 144–152. doi: 10.1016/j.jtbi.2011.10.034