INTRODUCTION {#s1}
============

Horizontal gene transfer (HGT) is the lateral movement of genetic material between cells ([@B1]) and has received increased attention due to the rapid emergence of multidrug-resistant bacteria ([@B2], [@B3]). HGT enables the incorporation of exogenous DNA to obtain new virulence and resistance traits. Plasmids are extrachromosomal entities that often confer novel and advantageous traits to the host. The plasmidome refers to the entire plasmid content of a microbial population, representing mobile genetic material that can be subjected to HGT. Plasmid-mediated HGT in bacteria can occur via conjugation (transfer of genetic material between bacterial cells by direct cell-to-cell contact), transduction (transfer of genetic material to a bacterium via a virus), and natural transformation (uptake of exogenously available genetic material by a bacterium) ([@B4]). Of these mechanisms, natural transformation is ostensibly the simplest, needing merely the presence of a competent bacterial cell to take up an exogenously available DNA fragment(s). Natural transformation occurs when a cell is competent, a highly regulated physiological state that reflects a window in which the proteins required for DNA binding, processing, and internalization are produced. The DNA taken up can be used as a nutrient source or a source of nucleotides for DNA synthesis, for genome repair, or to acquire variant alleles ([@B5][@B6][@B8]). Natural transformation has been detected in bacteria from all trophic and taxonomic groups, including *Archaebacteria*, suggesting that transformability evolved early in phylogeny ([@B9]). Thus, it has a significant impact on bacterial population dynamics as well as on bacterial evolution and speciation ([@B9]).

The Oak Ridge Field Research Center (ORFRC) ([@B10][@B11][@B13]) is a well-studied United States Department of Energy site that includes both areas with and without metal contamination, referred to as the contaminated and background sites, respectively. In a recent report ([@B14]), the plasmidome of this background site revealed the presence of a highly abundant 8-kb native plasmid containing putative mercury resistance genes (*mer*), despite a lack of detectable mercury contamination in the groundwater. It is not clear if the *mer* genes found at the background site are functional or how they came to be at that location. To get a deeper understanding, we examined the putative *mer*-encoding plasmid, p5343. We synthesized this plasmid and confirm that it provides functional mercury resistance in the model microbe Escherichia coli. We examined the phylogeny of the mercury reductases at ORFRC to obtain an understanding of its prevalence and distribution along with discussing the implications of our findings.

We used the synthesized native plasmid to better understand its potential for HGT in a number of bacterial isolates from this ORFRC environment. Since groundwater is a low-cell-density environment with fluctuating populations, we used natural transformation to examine the dissemination of the native plasmid-encoded metal resistance trait to bacterial isolates from the ORFRC site. Our results from these HGT assays are presented, and they reveal that natural transformation may play a critical role in the spread of resistance genes to divergent strains.

RESULTS {#s2}
=======

We examined the mercury resistance genes encoded on p5343, an 8-kb native plasmid belonging to incompatibility IncA/C core gene plasmid multilocus sequence type (cgPMLST), MOB group MOB~Q~, also encoding genes involved in plasmid replication and mobilization. The lack of detectable mercury in the source groundwater for p5343 ([@B14]) raised the possibility that the *mer* genes on this plasmid are not functional. Vestigial *mer* operons ([@B15], [@B16]) have been reported, and the closest MerA that has been experimentally confirmed to confer mercury resistance has only 47% amino acid identity ([@B17]) to p5343 MerA (closest functional MerA identified based on PaperBLAST \[[@B18]\]). Therefore, to functionally examine the putative *mer* genes of p5343, E. coli was transformed with the *mer*-containing plasmid p5343_UC57 ([Fig. 1](#fig1){ref-type="fig"}) or with a control plasmid (pUC57) that lacks the *mer* genes. The transformed strains were assayed for mercury resistance at mercury chloride concentrations known to be inhibitory to E. coli ([@B19]). Consistent with a functional *mer* system, we observed higher mercury chloride half-maximal inhibitory concentration (IC~50~) values in E. coli carrying p5343_UC57 relative to the same strain with pUC57 ([Fig. 2](#fig2){ref-type="fig"}). Specifically, at 14.3 μM mercury chloride, the strain with p5343_UC57 had a clear growth advantage. This is similar to the 10 μM mercury chloride concentration at which a growth advantage was reported earlier in E. coli transformed with mercury resistance genes ([@B19]). The proposed mercury reduction schematic based on previous illustrations ([@B15], [@B20], [@B21]), and SPOCTOPUS ([@B22])-predicted topology of the *mer* operon in p5343, is depicted in [Fig. 3a](#fig3){ref-type="fig"}. The p5343 encodes genes providing narrow-spectrum mercury resistance \[to inorganic mercury compounds like HgCl~2~ and Hg(NO~3~)~2~\] while lacking genes encoding broad-spectrum mercury resistance (like methylmercury and phenyl mercuric acetate). When p5343_UC57 is cloned into E. coli, the genes *merP*, *merT*, and *merF* likely aid in the transport of mercury while *merA* reduces Hg^+2^ to elemental mercury which volatilizes, leaving the bacterial environment mercury-free, resulting in improved mercury resistance. Thus, the native plasmid p5343 encodes a functional *mer* system that provides resistance to mercury.

![Plasmid map of p5343_UC57. The p5343 part (depicted in blue) encodes helix-turn-helix domain protein (HTH), mercuric ion reductase (MerA), mercuric ion uptake protein (MerF), hypothetical protein (Hyp2), mercuric transport protein (MerP), mercuric transport protein (MerT), regulator of mercury resistance genes (MerR), mobilization protein A (MobA), mobilization protein C (MobC), hypothetical protein (Hyp3), plasmid replication protein (RepA), hypothetical protein (Hyp1), and RelE toxin (RelE), while the pUC57 part (depicted in gray) codes for promotor sequence for ampicillin resistance gene (AmpR promoter), ampicillin resistance marker (AmpR), origin of replication, Lac promoter, LEU2 selection marker (Leu2), and elements ensuring plasmid maintenance (CEN/ARS).](mSystems.00588-19-f0001){#fig1}

![Improved mercury resistance in Escherichia coli DH10B containing p5343_UC57 (blue) versus only the empty vector pUC57 (green), grown in the presence of various mercuric chloride concentrations.](mSystems.00588-19-f0002){#fig2}

![(a) Schematic model of mercury resistance encoded by the *mer* operon; the genes absent on plasmid p5343 are depicted in gray. (b) The consensus *mer* operon (not drawn to scale) is provided for reference along with the putative *mer* operons from the background and contaminated sites as available from published sources ([@B11], [@B24]).](mSystems.00588-19-f0003){#fig3}

Since elevated levels of mercury have been reported at the floodplain of East Fork Poplar Creek (contaminated site), which is close to the background site studied ([@B23]), the movement/flow of microbes between the two sites may be responsible for the presence of *mer* genes in a location with undetectable mercury. To evaluate this possibility, we used existing metagenome ([@B12]) and whole-genome sequences ([@B24]) to compare the *mer* sequences from the contaminated site with those from the plasmidome of the background site. The *mer* operon structures differed between the background and contaminated sites ([Fig. 3b](#fig3){ref-type="fig"}). Further, we performed phylogenetic analysis on the MerA sequences from the ORFRC along with known sequences available on NCBI for providing context. The overall pattern of MerA distribution ([Fig. 4](#fig4){ref-type="fig"}) was similar to the 16S rRNA-based distribution, as reported earlier ([@B25][@B26][@B27]). Notably, no inference about HGT can be drawn because the tree includes MerA sequences from plasmidome and metagenome analysis where the source microbe encoding the plasmid is not known. Increased diversity in mercury resistance genotypes has been reported in the presence of low mercury concentrations ([@B28]). We confirm the same observation---we find that while the homologs of MerA from the background site are distributed across several bacterial and archaeal phyla, those previously observed from the contaminated site cluster mostly with *Gammaproteobacteria*. Additionally, while we observed the presence of genes encoding both narrow- and broad-spectrum mercury resistance in contaminated sites, only those encoding narrow-spectrum resistance were found in the circular plasmids from background sites ([Fig. 3b](#fig3){ref-type="fig"}). Thus, it is evident that the background sites have a diverse set of *merA* genes, which appear different from the metal-contaminated areas.

![The evolutionary history based on MerA sequences was inferred by the maximum likelihood method using MEGA7 ([@B39]). The bootstrap consensus tree is inferred from 1,000 replicates, and bootstrap values (%) for well-supported nodes are indicated. Coherent phyla are indicated to the right. The MerA sequences found in the background site are indicated by pink triangles, while those reported previously from the contaminated site are indicated by red triangles.](mSystems.00588-19-f0004){#fig4}

Next, we tested if the native plasmid-encoded metal resistance trait can be disseminated into other bacteria from this environment. Since this environment has very low cell density, natural transformation might be a prevalent HGT mechanism. Given that the ecological relevance of natural transformation may be better understood by studying bacteria and plasmids (DNA sequence is known to influence its uptake \[[@B29], [@B30]\]) native to this site, we used p5343 to test a suite of groundwater isolates for HGT of the plasmid-encoded mercury resistance trait. We found that the strains tested fall into three distinct categories ([Fig. 5](#fig5){ref-type="fig"}). The first category includes the negatives, i.e., strains that do not show improved growth in mercury chloride in a plasmid-dependent manner. The second set included false positives, i.e., strains that showed improved plasmid-dependent growth in the presence of mercury chloride but were negative for colony PCR against the *merA* gene. The third category included the positives, i.e., strains that both showed improved growth in mercury chloride in a plasmid-dependent manner and were positive for colony PCR for the presence of *merA.* The positives included 17 strains belonging to six different genera and consisted of both Gram-positive (2 strains) and Gram-negative (15 strains) bacteria. Among these, natural competence had never been reported in the genera *Arthrobacter*, *Dermacoccus*, *Acidovorax*, and *Cupriavidus*. Thus, the screen identified 17 strains that acquired mercury resistance and likely took up the plasmid via natural transformation.

![Results of the HGT assay used to screen environmental isolates for uptake of a plasmid-encoded trait via natural transformation. The 58 bacterial isolates are depicted on a 16S rRNA sequence-based phylogenetic tree (constructed using maximum likelihood). The colored bars indicate different bacterial classes, namely, orange is *Actinobacteria*, red is *Cytophagia* (*Bacteroidetes* phylum), yellow is *Alphaproteobacteria*, peach is *Betaproteobacteria*, and pink is *Gammaproteobacteria*. The heat map plots optical density in the presence of p5343_UC57 compared to pUC57 in the presence of different mercury chloride concentrations normalized by the no-mercury-chloride control. Only the strains showing improved plasmid-based growth in the presence of mercury chloride were subjected to *merA*-based colony PCR. The strains positive for colony PCR are indicated in black circles while those that are negative are depicted by gray circles. The genera of the strains positive for natural transformation as depicted in the figure are *Arthrobacter* (Ar), *Dermacoccus* (De), *Acidovorax* (Ax), *Cupriavidus* (Cu), Acinetobacter (Ac), and *Pseudomonas* (Ps). The figure was constructed using iTOL ([@B47]).](mSystems.00588-19-f0005){#fig5}

To further examine uptake of this plasmid via natural competence, two *Pseudomonas* strains (from the 17 positives), along with model strain E. coli DH10B, were subjected to a more comprehensive natural transformation assay with replicates. This involved incubation of strains with plasmids for natural transformation, followed by growth under mercury stress for two serial transfers. The *Pseudomonas* strains 5 and 12, along with DH10B, demonstrated plasmid-dependent improved growth in the presence of mercury chloride over the two serial growth regimes ([Fig. 6](#fig6){ref-type="fig"}). For strain 12, the improved growth was more evident at a lower mercury chloride concentration than at higher concentrations. Following the second growth regime, all three strains were positive for p5343-encoded *merA* via quantitative PCR (qPCR) after the assay while being negative before the assay. Thus, we confirm HGT and functionality of the p5343-encoded *mer* operon in a model bacterium and field-relevant environmental bacteria.

![Confirmation of HGT in strains 5, 12, and DH10B. Growth curves from second growth regime in the presence of various mercury chloride concentrations are depicted when strains 5, 12, and DH10B were naturally transformed with plasmid p5343_UC57 (blue) and pUC57 (green). Assay design is provided in [Fig. S2](#figS2){ref-type="supplementary-material"}.](mSystems.00588-19-f0006){#fig6}

DISCUSSION {#s3}
==========

The presence of both background and metal-contaminated sites at ORFRC makes it an excellent system to study the prevalence and dissemination of metal resistance determinants in microbial communities. Analysis of mercury reductases found in the contaminated and background sites revealed a significant decrease in the gene diversity in the former compared to the latter. This might be due to better functioning of a particular homolog(s) of *merA* in the presence of high levels of heavy metals (particularly mercury) and/or the inability of certain phyla (encoding diverse *merA* homologs) to survive in such conditions, given the reduced taxonomic diversity in contaminated groundwaters. Alternatively, it is possible that the contaminated site sequences available do not capture the entire *merA* diversity, and a targeted plasmidome analysis from the contaminated groundwater can shed light on this matter. Given that the background site is off-gradient from the contaminant plume and that the mercury resistance determinants based on the phylogenetic tree vary between the two sites, there is a reduced likelihood of bacterial carryover from the contaminated to the background site. The high abundance of the *mer*-encoding plasmid in the absence of detectable mercury in the background site could indicate very effective mercury reduction and volatilization of any trace amounts of mercury that leached into the groundwater, by the resident microbiome community. A study demonstrated that plasmid persistence could be attributed to compensatory adaptation, along with brief periods of positive selection ([@B31]), which might be the most plausible explanation for the persistence of a metal resistance gene(s) in plasmids from the background site.

The spread of mercury resistance via conjugation has been well documented in soil microbial communities ([@B32], [@B33]). In contrast, since aquatic environments are typically characterized by low cell densities and availability of exogenous dissolved DNA (1.7 to 88 μg per liter \[[@B34]\]), we hypothesized that natural transformation could be an HGT mechanism in these systems. Based on our screen studying the dissemination of mercury resistance determinants via natural transformation, several environmental strains tested were either negatives or false positives. This could be attributed to the strains not being naturally competent in the conditions tested or being competent but not able to maintain/express the *mer* genes encoded on the plasmid or the *mer* gene provided resistance at mercury concentrations lower than that tested. The presence of the *relE* toxin gene on the p5343 plasmid could also influence the retention and expression of the plasmid in certain strains, limiting it to strains that carry the corresponding antitoxin gene. Interestingly, we identify 17 strains that acquired an environmental plasmid-encoded trait via natural transformation. These comprised 12 *Pseudomonas* strains, which was notable because prior work showed that linear double-stranded DNA (dsDNA) was effective whereas single-stranded DNA (ssDNA) or plasmid DNA did not work for natural transformation in *Pseudomonas* ([@B35]). In fact, in strains known to be naturally competent, genomic DNA is typically preferred over circular plasmid sequences ([@B35], [@B36]), a trait likely attributed to the requirement of linearizing the circular plasmid prior to uptake. Thus, it is significant that, employing a circular native plasmid, we demonstrate HGT via natural transformation into environmental bacteria, indicating this could be a viable route for plasmid dissemination in groundwater.

Successful HGTs frequently occur between closely related organisms ([@B37]), and the compositional similarity between the donor and the recipient genomes promotes homologous recombination leading to DNA acquisition from close relatives. Interestingly, although all the p5343-encoded genes are closest to *Alphaproteobacteria*, the strains that were transformed belong to *Betaproteobacteria*, *Gammaproteobacteria*, and Actinobacteria. Among the 58 strains, only one belonged to *Alphaproteobacteria*, and testing more strains from this class might result in more positives. Natural competence is a transient physiological cell state which allows DNA uptake under specific conditions ([@B9]) (an exception being Neisseria gonorrhoeae, where competence is constitutive \[[@B38]\]). Thus, the ideal set of conditions for competence is strain specific, and the screen can be further optimized by varying plasmid incubation times, the amount of plasmid added, the cell density at which the plasmid is added, changes in medium constituents, and testing other mercury concentrations. Additionally, in its native state, p5343 may be methylated, and further experiments with different methylated versions of the p5343 plasmid might enable HGT into additional bacterial phyla from this environment.

Overall, this study reveals that the highly abundant *mer-*encoding plasmids are functional in providing mercury resistance and capable of being horizontally transferred into relevant bacterial isolates from the ORFRC site. We demonstrate that natural transformation facilitates the interphylum transfer of genetic elements, suggesting that the transient presence of plasmid DNA in close proximity may be sufficient for HGT to occur in the groundwater communities. This suggests that the microbial community studied is likely robust in tolerating low stresses and possesses a latent ability to swiftly adapt to changes in the environmental stress levels using natural transformation.

MATERIALS AND METHODS {#s4}
=====================

Plasmid p5343. {#s4.1}
--------------

Identification of the 8-kb plasmid p5343 is described in a previous report ([@B14]). Briefly, the plasmidome of background groundwaters at the ORFRC was examined using an optimized alkaline hydrolysis method followed by Illumina sequencing and analysis. Based on sequence coverage, p5343 was determined to be highly abundant in the background wells studied ([@B14]). In the present study, p5343 was synthesized and cloned into the vector pUC57 using the GenBrick synthesis service (GenScript, Piscataway, NJ) for successful propagation in E. coli. The resulting plasmid p5343_UC57 and its sequence are available via Addgene (Addgene ID 126645) and contain an Escherichia coli compatible origin along with an antibiotic marker for propagation and maintenance in the strain E. coli DH10B. Plasmid purifications were performed using the miniprep plasmid extraction kit (Qiagen GmbH, Germany).

Analysis of *mer* genes. {#s4.2}
------------------------

For the phylogenetic analysis of MerA, amino acid sequences were obtained from the plasmidome ([@B14]) of the background site along with the metagenome ([@B11], [@B24]) and whole-genome sequences ([@B11], [@B24]) of the contaminated site at ORFRC (given the lack of availability of plasmidome data from the contaminated site). In addition, we also added publicly available MerA sequences from NCBI (accession numbers provided on the phylogenetic tree). The evolutionary history based on MerA sequences was inferred using MEGA7 ([@B39]). The MerA sequences were aligned using MUSCLE ([@B40]). The alignment was manually curated, and all the MerA sequences were cropped to a common length of 515 amino acids. The JTT matrix-based model ([@B41]) was used to construct maximum likelihood trees with 1,000 bootstrap replicates ([@B42]).

Mercury resistance. {#s4.3}
-------------------

The plasmids pUC57 and p5343_UC57 were transformed into E. coli strain DH10B strain separately. Five freshly transformed colonies were picked for each plasmid, and overnight cultures were prepared at 37°C. To compare the growth of DH10B transformed with p5343_UC57 and pUC57, 5 μl of preculture (final OD~600~ of 0.1) was inoculated into 95 μl of LB medium with carbenicillin (100-μg/ml final concentration) in a 96-well transparent flat-bottom tissue culture plate (Corning, Falcon, catalog number 353003) in the presence of 0, 3.5, 7.1, 14.3, 28.7, and 43.1 μM final mercury chloride concentrations. The plate was sealed with BreathEasy seals (E&K Scientific, Santa Clara, CA, USA) and grown at 37°C in a Tecan F200 microtiter plate reader (Tecan Group Ltd., Männedorf, Switzerland) with shaking, measuring the optical density at a wavelength of 600 nm at 20-min intervals.

HGT via high-throughput screen. {#s4.4}
-------------------------------

A previous study ([@B43]) described isolation of 192 strains from the ORFRC environment. Of these, we rearrayed 58 strains in a 96-well format and subjected them to the high-throughput screen to test for HGT (see [Fig. S1](#figS1){ref-type="supplementary-material"} in the supplemental material). The groundwater source, isolation conditions, GenBank accession numbers, taxonomic assignments, and 16S rRNA sequences of the 58 strains are provided in [Table S1](#tabS1){ref-type="supplementary-material"}. The strains are numbered 1 through 58 for ease of representation in figures and text. The design of the assay to screen for HGT in the bacterial isolates is described in [Fig. S1](#figS1){ref-type="supplementary-material"}. The arrayed isolates were recovered in a 96-well deep-well plate (Costar, Thermo Fisher Scientific, Waltham, MA, USA) by addition of 1 ml of appropriate growth (either LB or R2A) medium. The 96-well deep-well plate was grown in a Multitron plate shaker/incubator (Infors, Bottmingen, Switzerland) at 30°C with shaking at 700 rpm, overnight. The 96-well plate was centrifuged at 4,000 rpm for 10 min in an Eppendorf 5810R centrifuge, and the pellets were resuspended in 1 ml fresh growth medium. Approximately 250 μl of each culture was added to individual wells of two 96-well deep-well plates. Approximately 5 μl of 50-ng/μl p5343_UC57 was added to each well of one plate while pUC57 was added to each well of the second plate. Both the plates were incubated overnight at 30°C without shaking to enable natural transformation. Since improved competence has been reported in cells which are stressed ([@B8]), after 24 h, 40 μl of the undiluted overnight cultures from the two 96-well plates was moved to two 384-well flat-bottom transparent plates (Corning, catalog number 353003) using a Biomek FxP (Beckman Coulter) liquid-handling robot. To each quadrant of the 384-well plate, 40 μl of appropriate medium with 0, 112, 224, and 448 μM mercury chloride concentrations was added such that final concentrations were 0, 56, 112, and 224 μM mercury chloride, respectively. The plates were sealed using BreathEasy seals (E&K Scientific, Santa Clara, CA, USA). The growth of both 384-well plates was monitored over 24 h at 30°C in a Tecan F200 microtiter plate reader with shaking, measuring the optical density at a wavelength of 600 nm at 20-min intervals (growth data provided in [Table S2](#tabS2){ref-type="supplementary-material"}). Since certain strains might already encode a native *merA* gene, we evaluated for improved mercury resistance in the presence of p5343_UC57 versus pUC57. The strains with plasmid-dependent improved growth in the presence of mercury chloride were subjected to colony PCR using primers 5′-CACACCGCCCAAAAGTCTAT-3′ and 5′-CAGAGCTGGCACAGATGATG-3′ designed to amplify the *merA* gene. All the cultures tested were confirmed to be negative for *merA* by colony PCR before the assay. For all strains that were positive for *merA* colony PCR postassay, an additional round of 16S rRNA sequencing was done to confirm the identity of the strains.

10.1128/mSystems.00588-19.1

Design of the HGT assay to screen environmental bacterial isolates for uptake of a plasmid-encoded trait via natural transformation. Download FIG S1, TIFF file, 25.1 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

10.1128/mSystems.00588-19.3

Details of strain names, isolation conditions, taxonomic assignments, and 16S rRNA sequences of 58 strains subjected to HGT assay screen. Download Table S1, XLSX file, 0.03 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

10.1128/mSystems.00588-19.4

Raw and normalized (to no-mercury control) ratios of optical density measured at 600 nm in the presence of p5343_UC57 compared to pUC57 after 24 h of growth in the presence of different mercury concentrations. Download Table S2, XLSX file, 0.02 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

To obtain a phylogenetic tree depicting the 16S rRNA sequences of all 58 strains tested, the 16S rRNA sequences of these strains were aligned using MUSCLE ([@B40]). The alignment was manually curated, and evolutionary history was inferred using the minimum evolution (ME) method ([@B44]). The evolutionary distances were computed using the maximum composite likelihood method ([@B45]) and are in the units of the number of base substitutions per site. The ME tree was searched using the Close-Neighbor-Interchange (CNI) algorithm at a search level of 1. The Neighbor-Joining algorithm ([@B46]) was used to generate the initial tree. Evolutionary analyses were conducted in MEGA7 ([@B39]).

Confirmation of HGT in selected strains. {#s4.5}
----------------------------------------

The design of the natural transformation assay for confirmation of HGT in selected strains is depicted in [Fig. S2](#figS2){ref-type="supplementary-material"}. Frozen stocks of the ORFRC isolate strain 5 (*Pseudomonas* FW305-3-2-15-A-LB2) and strain 12 (*Pseudomonas* FW305-20), along with Escherichia coli DH10B, were recovered overnight in culture tubes with 5 ml LB medium grown at 30°C. The cultures were centrifuged, and fresh 2.5 ml LB medium was added to the pellet, followed by resuspension. About 250 μl of each culture was added to individual wells of a 96-well deep-well plate (Costar, Thermo Fisher Scientific). To each well, 5 μl of the appropriate plasmid (p5343_UC57 or pUC57) at a concentration of 50 ng/μl was added. The plate was incubated at 30°C without shaking to enable natural transformation. After 24 h, 50 μl of the undiluted overnight cultures with plasmids from each well was transferred to four wells of a 96-well transparent flat-bottom tissue culture plate (Corning, Falcon, catalog number 353003) and 50 μl of LB medium with 2× mercury chloride concentrations was added to each well such that the final concentrations were 0, 56, 112, and 224 μM mercury chloride. The growth was monitored over 24 h at 30°C in a Tecan F200 microtiter plate reader with shaking, measuring the optical density at a wavelength of 600 nm at 20-min intervals, after which 5 μl from each plate was added to a fresh 96-well plate. To perform a second growth regime assay in the presence of mercury chloride, to each well 45 μl of LB and 50 μl of LB with 2× mercury chloride were added (to reach the same final mercury chloride concentrations), and the growth was monitored over 24 h at 30°C in a Tecan F200 microtiter plate reader. The 24-h growth data were analyzed to see if the strains with p5343_UC57 had an advantage over those with pUC57 when grown in the presence of mercury chloride. In addition, all strains were subjected to quantitative PCR with primers 5′-CGTCCTTGTCGAAGGTTTGT-3′ and 5′-ATAGACTTTTGGGCGGTGTG-3′ designed to amplify p5343-encoded *merA* both before (using glycerol recovered overnight cultures) and after the assay (qPCR raw data provided in [Table S3](#tabS3){ref-type="supplementary-material"}). The qPCR was performed using EvaGreen real-time PCR (Biotium, Hayward, CA, USA) according to the manufacturer's recommendations using 1 μl EvaGreen real-time PCR mixture (Biotium, Hayward, CA, USA), 1 μl of cell culture, 12.5 μl of Q5 high-fidelity 2× master mix (New England Biolabs, Ipswich, MA, USA), 1.25 μl each of 10 μM forward and reverse primers, and 9 μl of nuclease-free water.

10.1128/mSystems.00588-19.2

Design of the natural transformation assay for the confirmation of HGT in strains 5, 12, and DH10B (indicated by different shades of gray) in replicates. Download FIG S2, TIFF file, 25.1 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

10.1128/mSystems.00588-19.5

Results from qPCR depicting quantitation cycle (*C~q~*) values indicating the presence of *merA* gene in strains 5, 12, and DH10B after the second growth regime in mercury chloride in the comprehensive natural transformation assay. The lower the *C~q~* value, the higher initial copy number of the target gene. Download Table S3, XLSX file, 0.01 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

We thank Megan Garber (LBNL) and Nurgul Kaplan (LBNL) for help with liquid-handling systems. We thank Han Zhang (LBNL) for his help with figures. We thank John-Marc Chandonia (LBNL) for help with the environmental strain sequence data tracking.

The authors do not have any conflict of interest.

This work was part of the ENIGMA---Ecosystems and Networks Integrated with Genes and Molecular Assemblies (<http://enigma.lbl.gov>), a Science Focus Area Program at Lawrence Berkeley National Laboratory, and is supported by the U.S. Department of Energy, Office of Science, Office of Biological & Environmental Research, under contract number DE-AC02-05CH11231 between Lawrence Berkeley National Laboratory and the U.S. Department of Energy. The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

[^1]: **Citation** Kothari A, Soneja D, Tang A, Carlson HK, Deutschbauer AM, Mukhopadhyay A. 2019. Native plasmid-encoded mercury resistance genes are functional and demonstrate natural transformation in environmental bacterial isolates. mSystems 4:e00588-19. <https://doi.org/10.1128/mSystems.00588-19>.
