Diversity and role of plasmids in adaptation of bacteria inhabiting the Lubin copper mine in Poland, an environment rich in heavy metals

The Lubin underground mine, is one of three mining divisions in the Lubin-Glogow Copper District in Lower Silesia province (Poland). It is the source of polymetallic ore that is rich in copper, silver and several heavy metals. Black shale is also significantly enriched in fossil organic matter in the form of long-chain hydrocarbons, polycyclic aromatic hydrocarbons, organic acids, esters, thiophenes and metalloporphyrins. Biological analyses have revealed that this environment is inhabited by extremophilic bacteria and fungi. Kupfershiefer black shale and samples of water, bottom and mineral sediments from the underground (below 600 m) Lubin mine were taken and 20 bacterial strains were isolated and characterized. All exhibited multi-resistant and hypertolerant phenotypes to heavy metals. We analyzed the plasmidome of these strains in order to evaluate the diversity and role of mobile DNA in adaptation to the harsh conditions of the mine environment. Experimental and bioinformatic analyses of 11 extrachromosomal replicons were performed. Three plasmids, including a broad-host-range replicon containing a Tn3 family transposon, carried genes conferring resistance to arsenic, cadmium, cobalt, mercury and zinc. Functional analysis revealed that the resistance modules exhibit host specificity, i.e., they may increase or decrease tolerance to toxic ions depending on the host strain. The other identified replicons showed diverse features. Among them we identified a catabolic plasmid encoding enzymes involved in the utilization of histidine and vanillate, a putative plasmid-like prophage carrying genes responsible for NAD biosynthesis, and two repABC-type plasmids containing virulence-associated genes. These findings provide an unique molecular insight into the pool of extrachromosomal replicons and highlight their role in the biology and adaptation of extremophilic bacteria inhabiting terrestrial deep subsurface.


INTRODUCTION
Bacterial plasmids, as extrachromosomal mobile genetic elements, are components of many microbial genomes. They have modular structures, since it is possible to dissect them into several functional genetic modules. The plasmid backbone is composed of a set of conserved modules, coding for replication, stability and conjugal transfer functions, which are crucial for plasmid maintenance and spread. Many plasmids also carry accessory genes determining various phenotypes, e.g., antibiotics resistance and utilization of toxic compounds. Such genetic information is not crucial for host viability, but it may play an important role in the adaptation of bacteria to various environments, including terrestrial deep subsurface habitats (Sobecky and Coombs, 2009;Heuer and Smalla, 2012;Nojiri, 2013).
It was shown that one of the most frequently found phenotypic modules carried by bacterial plasmids are heavy metal resistance genes (Silver, 1996). Moreover, metal resistance genes are often co-localized on plasmids together with antibiotic resistance genes, and they are frequently present within transposable and integrative mobile elements (Rahube et al., 2014). Such plasmidencoded heavy metal resistance systems are usually related to chromosomally-encoded determinants found in other bacteria, which exemplifies the naturally occurring gene flow (Silver, 1996).
Horizontal gene transfer is a major mechanism contributing to bacterial diversification and adaptation, and plasmids are the main players in this process (Heuer and Smalla, 2012). Many plasmids are self-transmissible or mobilizable replicons, which can be transferred from one host to another, together with various "passenger" genes (often embedded within transposons), which can be useful under particular selection pressures (Tamminen et al., 2012). This phenomenon is extremely profitable for bacteria and, in effect, speeds up the process of evolution. To explain the link between environmental conditions and the diversity of bacterial plasmids a complex molecular analyses or, so called, metaanalyses of plasmidomes of bacteria inhabiting various (especially unique) extremophilic environments are needed (Dziewit and Bartosik, 2014).
In our previous work we analyzed several strains of the genus Pseudomonas originating from the Lubin copper mine (Szuplewska et al., 2014). Using a trap plasmid strategy to capture functional transposable elements, we identified several insertion sequences and transposons (autonomous and non-autonomous), some of which carry predicted genetic modules of adaptive value (Szuplewska et al., 2014). In this study we extended our investigation of the mobilome of bacteria inhabiting this deep underground environment. We isolated and analyzed the plasmids of strains representing various taxonomic groups, focusing on replicons conferring heavy metal resistance. Such plasmids (and their bacterial hosts) may be exploited in industrial processes, e.g. bioleaching of metals or biormediation. Moreover, the naturally occurring plasmids are good candidates for the construction of vectors for the genetic manipulations of biotechnologically important bacteria. However, this requires an in-depth understanding of the plasmid biology.

BACTERIAL STRAINS, PLASMIDS AND CULTURE CONDITIONS
The bacterial strains and plasmids used in this study are listed in Table S1. The strains were grown in LB (Luria-Bertani) medium (Sambrook and Russell, 2001) at 37 • C (E. coli) or 22 and 30 • C (other strains). Where necessary, the medium was supplemented with sucrose (10%) and the antibiotics, kanamycin (50-1000 μg/ml) and rifampicin (50 μg/ml).

SAMPLE COLLECTION AND BACTERIAL ISOLATION
Samples of black shale, mineral sediment, bottom sediment and water, collected from various sites within the Lubin underground copper mine (below 600 m), were placed in sterile plastic tubes and held at 4 • C until they were processed in the laboratory.
To isolate bacterial strains, samples of solid matter (10 g) were resuspended in 20 ml of 0.85% NaCl (pH 7.0) and shaken at 22 • C for 2 h. Then a series of dilutions were prepared in saline solution and plated onto solid LB medium. Diluted water samples were plated directly onto LB medium. The plates were incubated at 22 • C for 2 weeks. All operations were carried out aseptically.

AMPLIFICATION AND SEQUENCING OF 16S rRNA GENES
A colony PCR method was used for the amplification of 16S rRNA gene fragments (Gathogo et al., 2003). PCR was performed with the primers 27f and 1492r (Lane, 1991). The amplified 16S rDNA fragments were used as templates for DNA sequencing with an ABI Prism 377 automatic sequencer (Applied Biosystems).

PHYSIOLOGICAL ANALYSES OF THE BACTERIAL STRAINS
The temperature, pH and salinity tolerance of bacteria were analyzed by monitoring changes in the optical density of cultures (in comparison with non-inoculated controls) according to the procedures described previously . For the motility assay LB soft agar plates containing 0.3, 0.35, or 0.4% (w/v) agar were inoculated with the bacteria using a sterile toothpick and incubated at 30 • C for 48 h.
The ability to produce siderophores was examined by application of the modified chrome azurol S (CAS) agar plate method (Schwyn and Neilands, 1987). The plates were incubated at 30 • C for 72 h in the dark and the formation of halos around colonies was recorded.

DNA MANIPULATIONS, PLASMID ISOLATION AND INTRODUCTION OF PLASMID DNA INTO BACTERIAL CELLS
Bacterial plasmids were isolated according to the method of Birnboim and Doly (Birnboim and Doly, 1979), and when required, the DNA was further purified by CsCl-ethidium bromide gradient centrifugation (Sambrook and Russell, 2001). Plasmid DNA was also isolated using a Plasmid Mini Kit (A&A Biotechnology), Plasmid Miniprep Kit Gene Matrix (EUR x ) and GeneJET Plasmid Miniprep Kit (Thermo SCIENTIFIC). The visualization of mega-sized replicons was achieved by in-gel lysis and DNA electrophoresis according to a method described previously (Wheatcroft et al., 1990). The common DNA manipulations were performed as described earlier (Sambrook and Russell, 2001). PCR was performed in a Mastercycler (Eppendorf) using Taq DNA polymerase (Qiagen; with supplied buffer), dNTP mixture and appropriate primer pairs (Table S1). Triparental mating was performed as described previously (Bartosik et al., 2001).

PLASMID DNA SEQUENCING AND ASSEMBLY
The complete nucleotide sequences of plasmids pLM16A1, pLM20P1-5, pLM8P1, and pLM12P1 were determined in the Laboratory of DNA Sequencing and Oligonucleotide Synthesis (oligo.pl) at the Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Warsaw, Poland. Plasmids pLM19O1, pLM19O2, and pLM21S1 were sequenced at the Center for Biotechnology (CeBiTec), Bielefeld University, Germany. Highthroughput sequencing of the MID-tagged shotgun plasmidlibraries was performed using a Genome Sequencer FLX system (Roche/454 Life Sciences). The GS de novo assembler software (Roche) was applied for the sequence assemblies. Primer walking and polymerase chain reaction (PCR) were used to close physical gaps between assembled contigs. The amplified DNA fragments were sequenced using an ABI3730xl DNA Analyzer (Applied Biosystems).

IDENTIFICATION OF ACTIVE TRANSPOSABLE ELEMENTS (TEs)
The identification of functional TEs using trap plasmid pMAT1 was performed as described previously (Szuplewska and Bartosik, 2009;Dziewit et al., 2012).

PHAGE INDUCTION
Phage of Sinorhizobium sp. LM21 was induced using mitomycin C as described previously (Dziewit et al., 2014b).

BIOINFORMATIC ANALYSES
Plasmid nucleotide sequences were analyzed using Clone Manager Professional, version 9.0 (Sci-Ed software), Artemis (Carver et al., 2008) and GenDB 2.0 (Meyer et al., 2003). Similarity searches were performed using the BLAST programs (Altschul et al., 1997) and the PRIAM tool (Claudel-Renard et al., 2003). Putative tRNA genes were identified with the tRNAscan-SE program (Lowe and Eddy, 1997). Comparative genomic analyses were performed with the application of the ACT: the Artemis comparison tool (Carver et al., 2008). The reference data set for the computational prediction of metabolic pathways was obtained from the MetaCyc database (Caspi et al., 2014). Phylogenetic analysis was performed using MEGA6 (Tamura et al., 2013), applying the neighbor-joining algorithm with 1000 bootstrap replicates. The initial alignment obtained with ClustalW (Chenna et al., 2003) was manually refined. For the analysis 100 homologs (best BLAST hits) of proteins ArsB, ArsC, CzcD, and MerA retrived from the UniProt database (Apweiler et al., 2004) were used.
Preliminary physiological characterization of all identified bacterial strains revealed that 13 could grow at temperatures ranging from 15 to 37 • C, or even 42 • C (this growth pattern is typical for mesophilic bacteria), while 7 strains were psychrotrophs, able to grow at temperatures between 4 and 37 • C (the optimum temperature for all strains was either 21 or 30 • C) (Table S2). All strains grew in LB medium at pH values close to 7, which is characteristic for neutrophilic bacteria. Three strains (LM7, LM16, and LM19) displayed properties of facultative alkaliphiles, since they could tolerate pH values of up to 11 (Table S2) (Slonczewski et al., 2009). Salinity tolerance testing revealed that the majority (18) of the strains were non-halophilic bacteria, while LM14 and LM26 could tolerate a higher NaCl concentration (6%); thus they were classified as halotolerant (Table S2) (Larson, 1986). The motility assay revealed that 15 of the tested strains were motile (Table S2). Moreover, application of the universal chrome azurol S (CAS) agar plate assay indicated that 16 of the analyzed strains produced iron-chelating siderophores (Table S2).

GENERAL CHARACTERIZATION OF THE PLASMIDOME OF THE IDENTIFIED BACTERIA
Plasmid screening of all bacterial isolates revealed the presence of 12 circular replicons ranging in size from approx. 1.7 to 120 kb. Six strains contained plasmids: Achromobacter sp. LM16, Ochrobactrum sp. LM19, P. yeei LM20, Pseudomonas sp. LM8, Pseudomonas sp. LM12 and Sinorhizobium sp. LM21. The nucleotide sequences of all but one of these replicons [pLM8P2 of Pseudomonas sp. LM8 was described in our previous study (Szuplewska et al., 2014)], were obtained and analyzed. This revealed that the plasmids contained from 2 to 150 predicted genes and their average GC content varied between 53.5 and 67.3%. The plasmids carried as many as 22 types of predicted genetic modules, involved in various functions. The results of the overall characterization of the plasmids are shown in Table 1. A summary of genes identified in each plasmid, including their position, the size of the putative encoded proteins and their closest homologs, is presented in Tables S3-S7. For each plasmid, a conserved backbone, composed of the maintenance and conjugal transfer modules (including modules responsible for mobilization for conjugal transfer) was distinguished. In silico analysis of the plasmid-encoded accessory modules showed plasmids containing genes of direct adaptive value: (i) heavy metal resistance plasmids (pLM16A1, pLM20P1 and pLM20P2), (ii) a catabolic plasmid (pLM20P5), and other replicons, including a (iii) plasmid-like prophage (pLM21S1), (iv) putative virulence plasmids (pLM19O1 and pLM19O2) and (v) cryptic plasmids (pLM20P3, pLM20P4, pLM8P1 and pLM12P1). In our analyses we mainly focused on heavy metal resistance plasmids.

GENOMICS OF HEAVY METAL RESISTANCE PLASMIDS
Three plasmids were found to carry predicted genetic modules involved in heavy metal resistance: pLM20P1 and pLM20P2 of P. yeei LM20, and pLM16A1 of Achromobacter sp. LM16 (Figure 1, Tables S3, S4). The average GC content of the nucleotide sequences of the plasmids is 62.7, 63.4, and 64.2% (Table 1), respectively, which is lower than the mean values determined for the total DNA of the bacteria belonging to Achromobacter and Paracoccus genus-66.4% (mean result for 11 genomes) and 67.3% (mean result for 27 genomes), respectively.
Two plasmids of P. yeei LM20 are small replicons (6 and 6.2 kb, respectively) encoding replication initiation proteins with homology to related proteins of plasmids occurring in several carotenoid-producing strains of Paracoccus spp. (Maj et al., 2013). Plasmid pLM20P1 carries three arsenic-tolerance genes (ARS module; pLM20P1_p3-pLM20P1_p5) encoding arsenate reductase (ArsC), an efflux pump (ArsB) and an ArsR family transcriptional repressor. The small cytoplasmatic arsenate reductase ArsC is responsible for the conversion of arsenate (As 5+ ) into arsenite (As 3+ ), while the membrane protein ArsB is an effective chemiosmotic efflux system mediating the removal of As 3+ from the cell (Silver and Phung Le, 2005a,b). Interestingly, the ARS module of pLM20P1 shares 80% nucleotide sequence identity with the appropriate gene cluster of chromid pAMI5 of Paracoccus aminophilus JCM 7686 (Dziewit et al., 2014a). The second plasmid, pLM20P2 carries a CZC module, encoding a MerR family transcriptional regulator (pLM20P2_p3) and a predicted CzcD cation transport membrane protein (pLM20P2_p4) [member of the cation diffusion facilitator (CDF) protein family], which mediates cobalt (Co 2+ ), zinc (Zn 2+ ) and cadmium (Cd 2+ ) resistance (Anton et al., 2004). This DNA region shares high level (69%) of nucleotide sequence identity with appropriate CZC module of a large, extrachromosomal, alphaproteobacterial replicon-NT26_p1 of Rhizobium sp. NT-26 (Andres et al., 2013). The third plasmid, pLM16A1 of Achromobacter sp. LM16 (25 kb), has a mosaic structure. More than half of the plasmid genome shares at least 84% nucleotide sequence identity with transposon TNCP23 of plasmid pKLC102, coexisting as an autonomous replicon and a genomic island in Pseudomonas aeruginosa (Klockgether et al., 2004;Bonnin et al., 2013).
The MER module consists of 6 genes encoding proteins responsible for the enzymatic conversion of Hg 2+ ions to the less toxic form Hg 0 (metallic mercury). The key enzyme is the mercuric reductase (MerA), which reduces Hg 2+ to Hg 0 . The other five MER-encoded proteins include two mercury ion transporters (MerT and MerP), two transcriptional regulators (MerR and MerD) and an accessory membrane protein (MerE) of unknown function. The pLM16A1 MER module is embedded within a transposon (TnAO22a; 8240 bp), which is an isoform (99% nucleotide sequence identity) of TnAO22, previously identified in Achromobacter sp. AO22 (Ng et al., 2009). Using a positiveselection trap vector (pMAT1) to identify functional transposable elements, we demonstrated that TnAO22a is capable of transposition and therefore may contribute to the dissemination of mercury resistance.
We also performed a phylogenetic analysis of the heavy metal resistance modules of plasmids pLM20P1, pLM20P2, and pLM16A1. For the analysis ArsB, ArsC, CzcD, and MerA proteins were used. We found that the close homologs of ArsB and ArcC proteins are encoded mostly within the genomes of Alphaproteobacteria, including several strains of Paracoccus spp. (Figures S1, S2). Interestingly, proteins related to CzcD are encoded not only by Alphaproteobacteria, but also by members of the phylum Thaumarchaeota (Archaea), which may suggest the inter-domain gene transfer ( Figure S3). The fourth analyzed protein, MerA, has close relatives encoded within the genomes of various gram-negative (Beta-and Gammaproteobacteria), as well as gram-positive bacteria (Mycobacterium spp. and Bacilli) ( Figure S4).
Then, the functionality of the resistance modules (from plasmids pLM20P1, pLM20P2, and pLM16A1, and transposable elements ISPpu12a and Tn5563a) was tested in various hosts and the obtained values were compared with those of the wild-type strains (Figure 2). This analysis revealed that the introduction of the plasmids to the recipient strains resulted in a significant increase (at least 2-fold) in the MICs in 30 (20%) of the transconjugants, while decreases in the MICs of at least 2-fold were observed in 22 (14.7%) (Figure 2, Table S8).
The putative CZC module of plasmid pLM20P2 (CZC-LM20) was active, but produced a resistance phenotype in only a limited number of strains: Pseudomonas spp. LM5, LM6 and LM8 (increased resistance to Cd 2+ , Co 2+ , Zn 2+ ; MIC increases of at least 3-fold) and A. tumefaciens LBA288 (elevated resistance to Cd 2+ and Zn 2+ ; MIC increases of 2-and 2.5-fold, respectively). In contrast, the CZC-IS module of ISPpu12a did not confer resistance to Cd 2+ , Co 2+ , or Zn 2+ in any of the tested transconjugant strains and hence probably was inactive. In a few cases, presence of the CZC-LM20 and CZC-IS modules resulted in a decrease in the MIC values (Figure 2, Table S8). The arsenic resistance module of pLM20P1 (ARS-LM20) was tested for its ability to increase the tolerance of bacteria to arsenite and arsenate ions. Only 5 strains carrying pBBR-ARSLM20 (Pseudomonas spp. LM5, LM6, LM12, Sinorhizobium sp. LM21 and A. tumefaciens LBA288) exhibited at least 2-fold higher resistance to arsenite or arsenate compared to their parental wild-type strains (Figure 2, Table S8).
The MER modules analyzed in this study originated from plasmid pLM16A1 (MER-LM16 module) and transposon Tn5563a (MER-TN module) and were predicted to confer resistance to mercury ions. As mentioned above, the MER-LM16 module carries a complete set of genes commonly found in other well defined mercury resistance operons (Silver and Phung Le, 2005b), while MER-TN is a partial module, carrying only three genes (merRTP). Introduction of the MER-TN module resulted in a significant increase in the MIC for Hg 2+ in only two strains, while the MICs of six other strains were at least 2-fold decreased. In contrast, presence of the MER-LM16 module increased tolerance to Hg 2+ in 11 (73.3%) strains. Interestingly, introduction of plasmid pBBR-MERLM16 into two strains (LM6 and LM16) originally exhibiting high levels of resistance to mercury, resulted in an unexpected decrease in tolerance (Figure 2, Table S8).
We also tested the host range of plasmids pLM20P1, pLM20P2, and pLM16A1. The replication system of pLM16A1 was functional in all tested strains, which indicates its broad host range. In contrast, plasmids replication modules of pLM20P1 and pLM20P2 were found to have a relatively narrow host range, limited to a few strains of Alphaproteobacteria. Both plasmids were able to replicate in A. tumefaciens LBA288, and the latter one could also replicate in Sinorhizobium sp. LM21.
Another plasmid, pLM21S1 of Sinorhizobium sp. LM21, was found to carry numerous phage-related genes, encoding proteins required for DNA packaging, capsid and tail assembly, and cell lysis. This predicted plasmid-like prophage lacks an integrase gene, which suggests that it is not able to integrate into chromosomal DNA. Interestingly, pLM21S1 carries a RepC-like replication system, commonly found within large alphaproteobacterial plasmids (Cevallos et al., 2008). To check whether pLM21S1 is an active phage, we treated cells of Sinorhizobium sp. LM21 with mitomycin C, a classical inducer of lambdoid prophages. This approach did not cause phage induction, although it resulted in the induction of another indigenous prophage residing in the bacterial genome (Dziewit et al., 2014b). We speculate that pLM21S1 may be an inactive prophage, or alternatively it may require specific, as yet unidentified, environmental factors for induction.
Plasmid/phage pLM21S1 carries 147 putative genes and 3 tRNA-encoding sequences (Table S5). In total 83 of its genes are conserved in the genome of a related phage, RHEph10 of Rhizobium etli CFN42 (Santamaria et al., 2014) (Figure 3B).  NAD biosynthesis (B). The predicted genetic modules are indicated by black rectangles: ABC, ABC-type transporter system; ALN, allantoate amidohydrolase; COB, part of a cobalamine biosynthesis module; HIS, histidinol-phosphate aminotransferase; HUT, histidine utilization system; LCT, D-lactate dehydrogenase; MLR, microcystin LR degradation protein; MRS, multimer resolution system; NAD, NAD + biosynthesis module; PAR, partitioning system; REP, replication system; TA, toxin-antitoxin system; VAN, vanillate utilization system. Arrows indicate genes and their transcriptional orientation. The tRNA-encoding sequences are marked by black dots. The gray-shaded area connects genes of plasmid pLM21S1 and phage RHEph10 that encode homologous proteins.
These two replicons carry homologous replication and partitioning systems, as well as putative genetic modules for cobalamine biosynthesis (cobTS) (Figure 3B). The CobT and CobS proteins form a complex, which together with an additional subunit, CobN (not encoded by the prophages) catalyzes cobalt incorporation into the corrin ring during the biosynthesis of coenzyme B12 (Debussche et al., 1992). Interestingly, pLM21S1 carries also genes involved in nicotinamide adenine dinucleotide (NAD) biosynthesis, which are not present in RHEph10. Of these, pLM21S1_p57 encodes a bifunctional nicotinamide mononucleotide adenylyltransferase/ADP-ribose pyrophosphatase (NadM, EC 2.7.7.1), while pLM21S1_p58 encodes nicotinamide phosphoribosyltransferase (NadV, EC 2.4.2.12). NadV catalyzes the transformation of nicotinamide (NAm) to β-nicotinamide ribonucleotide (NMN), which is then converted to NAD by the NadM enzyme ( Figure 3B) (Martin et al., 2001). It has been suggested that the nadV-nadM module is not only responsible for NAD biosynthesis, but also for overall recycling of endogenous nicotinamide, which may be generated by the hydrolysis of NAD (Gerdes et al., 2006).
Next two replicons identified in the course of this study originated from Ochrobactrum sp. LM19. The plasmids pLM19O1 and pLM19O2 contain repABC modules (replication and partitioning functions) and carry 75 and 100 genes, respectively ( Figure S5, Table S6). The bioinformatic analyses revealed that both replicons encode several virulence-associated proteins, including: (i) outer membrane autotransporter barrel domain-containing proteins of the type V secretion system, (ii) three invasion-associated locus B (IalB-like) proteins, whose homologs determine the erythrocyteinvasive phenotype of Bartonella bacilliformis (Coleman and Minnick, 2001), (iii) a YadA-like autotransporter adhesin of the non-fimbrial adhesins family, whose representatives are responsible for adhesion in infected tissues and protection against lysis (Hoiczyk et al., 2000), and (iv) an ATP-binding ABC-type transporter (HlyB) and hemolysin D (HlyD) of the hemolysin A type I secretion system (Fath and Kolter, 1993) (Figure S5, Table S6).
The remaining four plasmids (pLM20P3 and pLM20P4 of P. yeei LM20, pLM8P1 of Pseudomonas sp. LM8 and pLM12P1 of Pseudomonas sp. LM12) identified in bacteria from the Lubin mine were designated as cryptic replicons, since they lack accessory modules (Tables S3, S7).

DISCUSSION
In this study we analyzed the plasmid content of 20 bacterial strains (representing eight genera of three classes of Proteobacteria and one genus of the Bacteroidetes phylum) isolated from the Lubin underground copper mine. The plasmidome of these strains comprised 12 replicons, whose genetic structure was revealed. We demonstrated the diversity of the plasmids and defined their adaptive value, focusing on heavy metal resistance, since these toxic elements are major contaminants of the Lubin copper mine environment influencing indigenous microoganisms.
The black shale horizon is exposed to the activity of oxygen, water, and microorganisms. These factors cause the chemical and biological weathering of the rock and strongly influences the geochemical cycles of heavy metals, as well as organic carbon, which results in further redistribution of the elements within the environment (Matlakowska et al., 2010(Matlakowska et al., , 2012Matlakowska and Sklodowska, 2011). The geochemical examination of environmental samples of water and bottom sediment used in presented study confirmed the presence of heavy metals. For example, the concentrations of Cu, Pb, Zn, As and Ni in sample of bottom sediment originated from underground pond were 6839 mg kg −1 , 1128 mg kg −1 , 115 mg kg −1 , 110 mg kg −1 and 14 mg kg −1 , respectively. Similarly, the mineral sediment taken from the surface of black shale contained copper at concentration 721 mg kg −1 .
High concentrations of heavy metals act as a specific selection pressure permitting the survival of only well adapted indigenous strains expressing multi-resistant and hypertolerant phenotypes. Our analysis confirmed that all the strains analyzed in this work exhibited such characteristics. The majority of them were resistant to high levels of copper, arsenic and nickel, which correlates with the composition of Kupferschiefer black shale [35000 mg kg −1 for Cu, 398 mg kg −1 for As, 479 mg kg −1 for Ni (Matlakowska et al., 2012)]. We showed that some of the resistance phenotypes were confirmed by the mobile genetic elements.
Among the mobile genetic elements identified in the analyzed bacteria we found three plasmids (pLM16A1, pLM20P1, pLM20P2) and three transposable elements [TnAO22a, Tn5563a, ISPpu12a] that mediate resistance to arsenate, arsenite, cadmium, zinc, cobalt, and mercury. In-depth analyses demonstrated that the resistance phenotypes conferred by the particular elements are highly dependent on the host strain. Surprisingly, these results revealed that acquisition of a predicted resistance module is not always beneficial for the host, and may paradoxically lead to increased sensitivity. We speculate that this phenomenon may be the result of an unfavorable influence of the introduced genetic module on the overall homeostasis of the cell, e.g., by altering intracellular ion concentrations. All of the analyzed resistance modules encode influx/efflux pumps that are likely to have a relaxed substrate specificity [e.g., ArsB recognizes and functions with antimony, as well as arsenite (Silver and Phung Le, 2005a,b)], which may significantly influence ion fluxes.
An interesting example of a resistance module leading to increased sensitivity is the MER module of Tn3 family transposon Tn5563a [deletion derivative of Tn5044 of Xanthomonas campestris; (Kholodii et al., 2000)]. Transposon Tn5563a lacks merA (which is present in Tn5044), but it carries merP and merT genes encoding transporters responsible for the uptake of toxic mercuric ions. In six of the tested strains, the acquisition of this genetic module resulted in significant decreases in their MICs, which is probably the result of an increased intracellular concentration of toxic ions plus the inability of the strains to inactivate Hg 2+ caused by the lack of MerA mercuric reductase. In contrast, increased tolerance was observed in two Pseudomonas spp. strains, possibly due to the presence of a chromosomal copy of the merA gene. Therefore it seems that the phenotypes determined by this module depend on the genetic background of the host strain.
The analyzed resistance modules are present on mobile genetic elements, so they may be readily transferred to other bacteria. We found that the replication system of the mobilizable plasmid pLM16A1 (containing MER) is functional in a broad range of hosts, which makes this replicon ideal for the horizontal dissemination of the mercury resistance phenotype. Moreover, the MER module of pLM16A1 (determining resistance in the great majority of tested strains) is located within an active Tn3 family transposon (TnAO22a), which enables its spread among different replicons co-residing in a single cell.
The host range of the tested plasmids correlates with the results of phylogenetic analyses, which showed that the close homologs of ArsB, ArcC, and CzcD proteins (whose genes are located within narrow host range plasmids, pLM20P1 and pLM20P2) are encoded mostly within the genomes of Alphaproteobacteria, while MerA protein (encoded by the broad host range plasmid pLM16A1) is conserved among various gram-negative and gram-positive bacteria. Interestingly, it was also found that the homologs of the MER module are frequently colocalized with antibiotic resistance genes in various plasmids of Gammaproteobacteria [e.g., pACM1 (Preston et al., 2014), pHCM1 (Parkhill et al., 2001), pKOX_R1 (Huang et al., 2013)]. Those findings reflect the pivotal role of broad host range replicons in the dissemination of various resistance genes.
In the course of this study other replicons which may contribute to adaptation of their bacterial hosts to harsh environmental conditions of the Lubin mine were also identified. Amongst them, there was a catabolic replicon pLM20P5 of P. yeei LM20. It carries genetic modules, enabling the utilization of histidine and vanillate. These substrates can be used as alternative sources of carbon, nitrogen and energy. Since there is a limited amount of easily degradable carbon sources in the Lubin mine, it is likely that acquisition of catabolic plasmids such as pLM20P5 may be beneficial to the host.
An intriguing finding was also the identification of an unusual plasmid-like prophage element in Sinorhizobium sp. LM21. The replicon pLM21S1 has features typical of both, plasmids (replication and stable maintenance modules) and phages (complete set of proteins necessary for the phage "life cycle"). Plasmid pLM21S1 is related to the phage RHEph10, which is able to infect various Rhizbium etli strains (Santamaria et al., 2014). Accordingly, both replicons may be considered as the archetypes of a novel group of plasmid-like prophages. Interestingly, pLM21S1 carries the nadV-nadM genes responsible for NAD biosynthesis. The introduction of such genes into bacteria that neither possess the NAD de novo biosynthesis pathway, nor the NAD salvage pathway (e.g., Haemophilus influenzae), may transform them into "NAD-independent" strains, able to synthesize NAD from nicotinamide (Martin et al., 2001). We postulate that the presence of pLM21S1 may be beneficial to the host because NAD production is likely to improve the overall fitness of the bacterium, since NAD (and its derivative, NADP) is the most important coenzyme in cellular redox reactions (Martin et al., 2001).
Current knowledge concerning the direction, frequency and range of DNA transfer among microorganisms living in extreme environments, and especially deep terrestrial habitats, is sparse. We believe that the approach followed in the present study, linking geochemical data with physiological characterization of bacteria and detailed analyses of their plasmidome, may considerably increase our understanding of the influence of mobile DNA and horizontal gene transfer on the biology of extremophilic bacteria.

ACKNOWLEDGMENT
This work was supported by the National Science Centre, Poland [grant number N N303 579238].

SUPPLEMENTARY MATERIALS
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fmicb. 2015.00152/abstract Figure S1 | Phylogenetic tree of the ArsB proteins homologous to ArsB efflux pomp encoded within pLM20P1. The analysis was based on 100 sequences retrieved from the UniProt database. The unrooted tree was constructed using the neighbor-joining algorithm and statistical support for the internal nodes was determined by 1000 bootstrap replicates. Values of ≥50% are shown. UniProt accession numbers of the protein sequences used for the phylogenetic analysis are given in parentheses.
The underlined accession numbers indicate Paracoccus proteins.