The Role of CzcRS Two-Component Systems in the Heavy Metal Resistance of Pseudomonas putida X4

The role of different czcRS genes in metal resistance and the cross-link between czcRS and czcCBA in Pseudomonas putida X4 were studied to advance understanding of the mechanisms by which P. putida copes with metal stress. Similar to P. putida KT2440, two complete czcRS1 and czcRS2 two-component systems, as well as a czcR3 without the corresponding sensing component were amplified in P. putida X4. The histidine kinase genes czcS1 and czcS2 were inactivated and fused to lacZ by homologous recombination. The lacZ fusion assay revealed that Cd2+ and Zn2+ caused a decrease in the transcription of czcRS1, whereas Cd2+ treatment enhanced the transcription of czcRS2. The mutation of different czcRSs showed that all czcRSs are necessary to facilitate full metal resistance in P. putida X4. A putative gene just downstream of czcR3 is related to metal ion resistance, and its transcription was activated by Zn2+. Data from quantitative real-time polymerase chain reaction (qRT-PCR) strongly suggested that czcRSs regulate the expression of czcCBA, and a cross-link exists between different czcRSs.


Introduction
Heavy metal pollution leads to serious ecological and health problems due to the toxic effects of metal ions and their accumulation throughout the food chain. Even the essential metal ions, Zn 2+ , Cu 2+ , and Co 2+ could be toxic at high concentrations [1][2][3]. Many microorganisms naturally possess the ability to inhabit metal-polluted areas. An advanced understanding of the mechanisms by which microbes cope with metal stress will facilitate the rational design of strategies for detection and bioremediation of heavy metal-polluted water and soil systems.
To cope with metal stresses, bacteria develop various detoxification mechanisms. The most effective way is the extrusion of excessive metal ions out of the cell by active transport. A number of studies over the past years have been devoted to the protein families that export metal ions. The protein families include P-type ATPases driven by ATP hydrolysis, cation diffusion facilitator family transporters, which act as chemiosmotic ion-proton exchangers, and the resistance nodulation division (RND) family of transporters that mediate proton-driven efflux [4]. Usually, more than one kind of efflux protein is found in metal resistance bacteria [5,6]. CzcCBA is a CBA transporter. This structural gene region encodes outer membrane factors CzcC, membrane fusion protein CzcB, and CzcA protein of the resistance-nodulation-cell division protein family [4]. The transcriptional responses of these exporting genes are controlled by different regulators, including MerR family regulators [7], the ArsR/SmtB family [8], and two-component systems [9][10][11].
Many adaptive bacterial responses to environmental changes are governed by two-component signal transduction systems [12]. In typical two-component systems, environmental variation is first sensed by a histidine sensor kinase capable of autophosphorylation. The phosphoryl group is then transferred to an aspartic acid residue of the response regulator. Usually, the phosphorylated response regulator binds to DNA, resulting in either activation or repression of target genes. Caille et al. [9] observed that in P. aeruginosa PAO1, Zn 2+ and Cu 2+ enhance the expression of czcRS. Phosphorylated CzcR then activates the expression of czcCBA operon encoding an efflux pump specific for zinc, cadmium, and cobalt. The copRS two-component system has been identified as key genes involved in copper resistance [13]. In Cupriavidus metallidurans CH34, the regulatory genes of czcCBA are arranged in an upstream region consisting of czcN and czcI, and a downstream region consisting of czcD, czcR, czcS and czcE [14]. CzcRS and a periplasmic copper-binding protein designated CzcE [15], exert metal-dependent control of czcNICBA expression via regulation of czcNp activity [16]. Grosse et al. [10] demonstrated that uninduced C. metallidurans CH34 with a mutation in czcS contains more czcCBA message and resumes growth faster when challenged. These pieces of evidence demonstrate the participation of two-component systems in the regulation of metal extrusion.
P. putida is a ubiquitous saprophytic bacterium, which has been extensively studied for the biodegradation of organic pollutants [17]. Although the genomes of P. putida contain several czcRS operons [5], no research has been reported regarding the simultaneous function of different czcRSs in the same strain. Our previous study indicated that P. putida X4 has a high tolerance and absorption capacity for metal ions [18]. In this study, PCR with primers designed from the homologous sequence in genome sequenced P. putida were used to detect the presence of czcRS two-component systems and the czcCBA operon. Because two czcRS two-component systems and one czcR3 without its corresponding sensor gene were detected, we examined their different roles and how they interacted with the czcCBA operon.

Minimum Inhibitory Concentrations (MICs) of the Wild-Type and Mutant Strains
The MIC values of the wild-type and mutant strains are shown in Table 1. Determination of heavy metal resistance showed that the strain X4 displayed a high level of Cd 2+ , Co 2+ , Cu 2+ , and Zn 2+ resistance. The functionalities of czcS1, czcS2, and czcH were studied by determining the resistance of gene-deficient strains to different divalent heavy metals. ΔS1 tended to be sensitive to Cd 2+ , Co 2+ , and Zn 2+ ions. The MIC of Co 2+ and Zn 2+ declined to 0.5 and 2.25 mmol·L −1 . The czcS2 mutation decreased the strain's resistance to Cd 2+ and Co 2+ . Both czcS1 and czcS2 were required for generating full cadmium resistance as the ΔS1S2 strains with both interrupted czcS1 and czcS2 could tolerate only 1.25 mmol·L −1 Cd 2+ , which is lower than that tolerated by single gene mutant strains. Small but reproducible repression was observed when ΔS1 and ΔS2 were exposed to copper ions. A putative open reading frame (ORF), designated as czcH in this study, just downstream of czcR3, was identified as necessary to generate full metal tolerance. The resistance of ΔczcH to Cd 2+ , Co 2+ , and Zn 2+ were lower than that of the wild-type strains. The results presented above demonstrate that czcS1, czcS2, and czcH are necessary for the wild-type strain to generate full resistance to Cd 2+ , Co 2+ , Cu 2+ and Zn 2+ .

Genetic Complementation of Mutant Genes
Complementation studies were carried out to confirm that the decreased resistance of strains ΔS1, ΔS2, and ΔczcH to divalent metal ions was caused by the mutation of target genes. For the complementation analysis, the complete ORF of czcS1, czcS2 and czcH was cloned into the broad-host-range vector pVLT31, resulting in plasmids pVLTS1, pVLTS2 and pVLTH. Strains ΔS1c, ΔS2c and ΔczcHc were obtained by mating the complementation plasmid into P. putida ΔS1, ΔS2 and ΔczcH, respectively. The MICs of ΔS1c, ΔS2c and ΔczcHc were tested in Tris-buffered medium with 0.1 mmol·L −1 IPTG using the method described above. The complementation data are shown in Table 1. More than 85% of the tolerance to the tested metals was recorded for ΔS1c and ΔczcHc. The resistance of ΔS2c to Co 2+ and Cd 2+ were recorded at 100% and 94%, respectively. At the same time, the MICs of control strains did not improve.

Cd 2+ and Zn 2+ Repressed the Transcription of czcRS1, Induced the Transcription of czcRS2 and czcH, Respectively
Expression patterns and inducer specificities were determined by measuring the changes of LacZ activity and verified by qRT-PCR. All promoters were tested for their inducibility in the presence of 0-180 μmol·L −1 Zn 2+ , Co 2+ , Cu 2+ , and Cd 2+ . The expression patterns and inducer specificities of czcS1 are shown in Figure 1a. The czcS1 is transcribed without any metal ions. The transcription of czcS1 dropped to 40% in the presence of 180 μmol·L −1 Zn 2+ , then to 34% after a 180 μmol·L −1 Cd 2+ treatment. Compared with Zn 2+ and Cd 2+ , Cu 2+ and Co 2+ moderately decreased the transcription of czcS1. The transcription of czcS2 was only induced by Cd 2+ ; LacZ activity in ΔS2 increased about 24-fold after exposure to 180 μmol·L −1 Cd 2+ (Figure 1b). The czcH gene was significantly induced by Zn 2+ , but was not affected by Co 2+ , Cu 2+ and Cd 2+ . The highest level of czcH transcription was found at 100 μmol·L −1 Zn 2+ and it was decreased when the Zn 2+ concentration was higher than 100 μmol·L −1 (Figure 1c). These observations reveal that transcription from the promoters of czcS1, czcS2 and czcH are different, and that the main inducers or repressors are Zn 2+ and Cd 2+ . qRT-PCR data confirmed the inducibility presented above ( Figure 2); in the wild-type strain, addition of 300 μmol·L −1 Cd 2+ or Zn 2+ to the culture resulted in down-regulation of czcRS1 and up-regulation of czcRS2 and czcH, respectively (30-, and 35-fold). Besides, the transcription of czcCBA was strongly induced by Cd 2+ (90-fold) and Zn 2+ (47-fold) to keep the equilibrium of these metal ions in the cells.

The Transcription Variation of Mutant Strains
After treatment with Cd 2+ or Zn 2+ , transcription of czcA was determined by qRT-PCR in the wild-type X4 strain and mutant strains. As shown in Figure 3, the mutation of czcS1 and czcS2 reduced the mRNA amount of czcA to 14%, and 23% that of the wild-type, respectively, in the presence of Cd 2+ . A similar trend was detected upon induction with Zn 2+ . The mutation in czcH also caused repressed transcription of the czcA. These findings show that all these three genes are required to maintain the expression of czcA at a high level. The transcript variation of czcR1 and czcR2 is shown in Figure 4. Compared with the wild-type, interruption of czcS1 enhanced the transcription of czcR1 by 3.22-and 2.5-fold, in the presence of Cd 2+ and Zn 2+ , respectively. Mutation of czcS2 reduced the czcR1 and czcR2 mRNA to 40% and 2% of the wild type, while cultured with Cd 2+ . These transcription variations indicated that czcS1 in wild-type act negatively on the transcription of czcR1, and that functional czcS2 is needed to keep transcription of czcR1 and czcR2 at a relatively high level.   Table S1.

Discussion
In the present study, we demonstrate the roles played by the czcRS two-component systems in the control of czcCBA transcription as well as their different inducibility. To our knowledge, this is the first report on the functional difference of czcRS two-component systems in P. putida. Our results clearly show that all czcRS two-component systems are needed for P. putida X4 to generate full metal resistance. The reductions in MICs of Zn 2+ , Co 2+ , and Cd 2+ in the czcS mutant strains were less than eight-fold. It cannot be excluded that czcCBA is also induced by other regulators, or that other transporters present in X4 may contribute to the intrinsic resistance of this organism to heavy metals. The amplified czcA in this study which is induced by Zn 2+ , and Cd 2+ in wild type revealed a >99% nucleotide identity with the czcCBA1 operon in P. putida KT2400. The czcCBA1 operon in P. putida KT2440 plays a major role in generating metal resistance to Zn 2+ , Cd 2+ , and possibly Pb 2+ [6]. The czcCBA found in C. metallidurans CH34 confers resistance to Zn 2+ , Co 2+ , and Cd 2+ [10]. Similar systems have been found in P. aeruginosa PAO1 and P. aeruginosa CMG103 [9,19].
Plasmid pVIK112 was constructed by Kalogeraki et al. [20] to simultaneously disrupt and create fusions to target genes of diverse bacteria. The sequence of the fragment amplified in this study and the complete genome analysis of P. putida KT2440 [5] revealed that the czcRS in P. putida form separate operons and czcSs is located at the end of the operons. Therefore, a minimal polar effect was generated after single cross-over with pVIK112. The remarkable increase in metal resistance in the complementation strains supports this idea. Combining the data from the β-galactosidase activity assay and the transcription variation of czcRs and czcCBA, the cross-interaction pathway among Zn 2+ , Cd 2+ , czcRS and czcCBA was depicted ( Figure 5). The control of czcCBA expression through czcRS has been studied in P. aeruginosa and C. metallidurans [9][10][11]. Compared with previous studies, the regulation of czcCBA in P. putida was different. P. putida evolved two czcRS operons. Simultaneously, other regulation genes, such as czcN and czcI were not found in the P. putida sequenced genome. The expression of the CzcCBA efflux pump was induced by two czcRS operons with different inducibilities. The czcH gene also influenced the transcription of czcCBA. The function of czcRS1, czcRS2 and czcH in P. putida X4 may be complementary. The transcription of czcRS1 was inhibited in the presence of Zn 2+ and Cd 2+ , which may caused the reduced transcription of czcCBA. However, a high level of efflux pump CzcCBA expression is essential for coping with high concentrations of metal ions (Figure 2). The low expression trend of czcCBA resulting from the CzcR1 may be compensated for by the induction effect caused by czcRS2 and czcH, which were induced with Cd 2+ and Zn 2+ , respectively. Complete czcRS1 negatively feedbacks on its expression, which creates a common regulating pathway in prokaryotes for maintaining the expression of transcription factors at appropriate levels [21]. Complete czcRS2 has a positive effect on the transcription of itself. This positive autoregulation loop could allow the cells to respond rapidly to the presence of small amounts of heavy metals in the environment [11]. Based on the data from the LacZ fusion assay and qRT-PCR, the cross-talk among Zn 2+ , Cd 2+ , czcRS, and czcA was depicted. Zn 2+ and Cd 2+ repressed the transcription of czcRS1, and induced the transcription of czcRS2 and czcH, respectively. The czcRs and czcH gene could induce czcCBA transcription, leading to resistance to cobalt, zinc, and cadmium. Data from qRT-PCR also revealed a cross-interaction between czcRS2 and czcRS1.
In conclusion, there was a clear distinction between the different czcRSs. czcRS1 is constitutively expressed and possibly acts as a housekeeping resistance mechanism against Zn 2+ , Co 2+ , and Cd 2+ in P. putida X4. The transcription of czcRS1 was repressed by Cd 2+ and Zn 2+ , whereas the transcription of czcRS2 was induced only by Cd 2+ . Both czcS1 and czcS2 act positively on the transcription of czcA.
Although czcH plays an important role in the metal resistance of P. putida X4, the underlying mechanism merits further studies. Different czcRS genes, along with other regulators, are needed to make precise adjustments to the gene expression of P. putida to adapt to changes in metal ions in the environment.

Bacterial Strains, Plasmids, and Growth Conditions
Detailed information regarding the bacterial strains and plasmids are listed in Table 2. P. putida X4 was isolated from the waste soil around the Academy of Hubei Agricultural Sciences, Wuhan, China [22]. For construction and maintenance of plasmids, Escherichia coli strains DH5α and C118λpir were used. Bacteria were routinely grown in Luria-Bertani medium at 37 °C (E. coli) or 30 °C (P. putida). When required, antibiotics were added at the following final concentrations: for E. coli, 100 µg/mL ampicillin, 50 μg/mL kanamycin, and 12.5 μg/mL tetracycline; and for P. putida, 50 μg/mL kanamycin and 12.5 μg/mL tetracycline. Tris-buffered medium [26] containing 0.1% glucose was used for testing metal resistance and promoter inducibility. Analytical grade reagents CuSO4·5H2O, CdCl2·H2O, ZnCl2, and CoCl2·6H2O were prepared into 1.0 M stock solutions and sterilized by filtration.

Determination of MIC Values
The minimum inhibitory concentrations (MICs) of trace metals were determined as previously described [26]. Briefly, bacteria cultures were spread onto the Tris-buffered medium containing metal salts from 0.5 to 10 mmol·L −1 (0.25 mmol·L −1 interval). The MICs were defined as the lowest concentration of metal salts at which no CFU were observed after 60 h of incubation at 30 °C.

DNA Manipulation
All primer pairs used in this study were shown in Table 3. Axygene mini-prep, PCR purification, and gel purification kits (Axygene, Union City, CA, USA) were chosen for DNA separation and purification. PCRs were carried out using Taq DNA polymerase. The plasmids were constructed using standard recombinant DNA techniques and introduced into E. coli by transformation. Plasmids derived from pVIK112 [20] were introduced into P. putida by triparental conjugation. The helper strain was E. coli MM294 with the helper plasmid pRK2073 [24]. The essential DNA fragment was verified by DNA sequencing.

Homologous Recombination for Construction of czcS Mutants and lacZ Fusion Reporter Strains
The suicide fusion vector pVIK112 was selected for the mutation and construction of a lacZ fusion reporter system. The plasmid pVIKS1 containing a czcS1 fragment was constructed as follows: the czcS1 fragment was cloned into pTA2 and its orientation was confirmed by sequencing. After digestion with XbaI and KpnI, the restriction fragment was cloned into pVIK112 to create pVIKS1. The same procedure was then employed to create plasmids pVIKS2 and pVIKH. These plasmids were introduced into P. putida X4 by triparental conjugation. To generate double czcS mutation strain ΔS1S2, another suicide plasmid pRRT was constructed by replacing lacZYA and the kanamycin resistance gene cassette in pVIK112 into a tetracycline resistance gene by PstI digestion. After ligating with czcS2 restriction fragment from pTA2, the plasmid was applied to generate double mutation based on the ΔS1.

β-Galactosidase Activity Assay
The activity of β-galactosidase was measured according to the method by Nies [27] using o-nitrophenyl-D-galactoside (4 mg/mL) as the substrate. The unit 1 U was defined as the activity that forms 1 nmol of o-nitrophenol per minute at 30 °C.

Complementation Experiment
Complementation experiments were carried out to confirm that the phenotype variation is caused by the czcS mutant. Two czcS and the czcH open reading frames (ORFs) were amplified using the total DNA from P. putida X4 as the template. After digestion with XbaI and KpnI, the amplified fragments were ligated into the broad-host-range expression vector pVLT31, which contains the lacI and tac promoter [25], yielding pVLTS1, pVLTS2, and pVLTH. The yielded plasmids were mated into mutant strains ΔS1, ΔS2, and ΔczcH respectively via conjugation. The MICs of the different transconjugants were checked as described above with IPTG induction. Control strains were constructed by transformed pVLT31 into ΔS1, ΔS2, and ΔczcH.

RNA Isolation and Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) Analysis
After Cd 2+ or Zn 2+ treatment, the total RNA for qRT-PCR was extracted using the RNA bacteria protect solution and the RNeasy kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. The residual DNA was eliminated by treatment with 20 U RQ1 RNase-free DNase (Promega, Madison, WI, USA) followed by phenol-chloroform extraction. The RNA was precipitated with ethanol and resuspended in RNase-free water. For cDNA synthesis, 3 μg of RNA was reverse transcribed using random hexamer primers and MLV reverse transcriptase (Promega) according to the supplier's instructions. Reverse transcriptase was inactivated by incubation at 70 °C for 15 min and the obtained cDNAs were stored at −20 °C until utilization. The cDNAs were quantitatively measured with a Bio-Rad iCycler machine (Bio-Rad, Berkeley, CA, USA) using a Sybr Green Quantitect kit (Qiagen). The primer pairs used for qRT-PCR were designed using the Primer3 [28] program. To check for contaminated DNA, control reactions without reverse transcriptase were analyzed using rpsL primer set. The pure RNA samples were obtained when no amplification of the non-template control was detected. The cDNA samples were diluted 10-fold and 1 μL of the dilution served as the PCR template, which