0SCIENTIFIC REPpRTS
Received: 08 July 2016 Accepted: 21 November 2016 Published: 14 December 2016
Arsenic resistance strategy in Pantoea sp. IMH: Organization, function and evolution of ars genes
Liying Wang1,2, Xuliang Zhuang1,2, Guoqiang Zhuang1,2 & Chuanyong Jing1,2
Pantoea sp. IMH is the only bacterium found in genus Pantoea with a high As resistance capacity, but its molecular mechanism is unknown. Herein, the organization, function, and evolution of ars genes in IMH are studied starting with analysis of the whole genome. Two ars systems - arsl (arsR1B1C1H1) and ars2 (arsR2B2C2H2) - with low sequence homology and two arsC-like genes, were found in the IMH genome. Both arsl and ars2 are involved in the As resistance, where arsl is the major contributor at 15 °C and ars2 at 30 °C. The difference in the behavior of these two ars systems is attributed to the disparate activities of their arsR promoters at different temperatures. Sequence analysis based on concatenated ArsRBC indicates that arsl and ars2 clusters may be acquired from Franconibacter helveticus LMG23732 and Serratia marcescens (plasmid R478), respectively, by horizontal gene transfer (HGT). Nevertheless, two arsC-like genes, probably arising from the duplication of arsC, do not contribute to the As resistance. Our results indicate that Pantoea sp. IMH acquired two different As resistance genetic systems by HGT, allowing the colonization of changing ecosystems, and highlighting the flexible adaptation of microorganisms to resist As.
Arsenic (As) is a toxic element present in many environmental biotopes. Inorganic As exists primarily in two valence states: As(III) and As(V). To resist the disruptive effects of As, microbes have evolved a variety of mechanisms, including As(III) oxidation through the activity of As(III) oxidase and As methylation by methyltrans-ferase. Microorganisms can also utilize As in metabolism either as a terminal electron acceptor in dissimilatory As(V) respiration or as an electron donor in chemoautotrophic As(III) oxidation. Nevertheless, the most universal and well-characterized As resistance mechanism is induced by the ars system1.
The content and organization of the ars system vary greatly between strains. Most of the core genes in ars operons contain arsR, arsB and arsC, and other genes are also reported, such as arsA, arsD, arsT, arsX, arsH, and arsN2. A common organization of the ars cluster is arsRBC3, whereas duplicate ars operons can also be found in a single strain, such as Pseudomonas putida KT24404. Interestingly, in this case more than one ars cluster with different structures is observed in the same strain. In all, the content and organization of ars operons exhibit great diversity and complexity, and subsequently contribute to the As resistance capability of strains as summarized in Table S1.
The complexity of the ars system in diverse bacteria raises the question of its origin and evolution. Different evolution theories have been advanced for the evolutionarily old proteins, efflux pump protein (ArsB) and As(V) reductase (ArsC). For example, ArsB and ArsC may have evolved convergently, as evidenced by sequence analyses5. In contrast, arsC genes are reported to have a common origin and may have been transferred to other domains by HGT in early times, followed by subsequent divergence to the current phylogeny6. Follow-up studies suggest that the HGT events of ars genes may be common in nature7.
Pantoea is a genus of Gram-negative, facultative anaerobic bacteria. This genus belongs to gamma Proteobacteria, family Enterobacteriaceae, and was recently separated from the genus Enterobacter. Currently, the genus contains twenty-six species (http://www.bacterio.net/pantoea.html). Members of this genus are found in various environmental matrices8,9. In 2013, strain Pantoea sp. IMH was isolated for the first time from As-polluted groundwater and reported to resist high concentrations of As, up to 150 mM As(V) and 20 mM As(III)10. However, the hyper As-resistance strategy employed by Pantoea sp. IMH remains unclear.
■"■State Key Laboratory of Environmental Chemistry and Ecotoxicology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, P.O. Box 2871, Beijing 100085, China. 2University of Chinese Academy of Sciences, Beijing 100049, China. Correspondence and requests for materials should be addressed to C.J. (email: cyjing@rcees.ac.cn)
Figure 1. Analysis of co-transcript unit in the ars clusters of Pantoea sp. IMH by RT-PCR. (a) Map position of ars genes and the primers for RT-PCR analysis. Primers used and amplified products (numbered) are given below the schematic representation of the genes. (b) Result of RT-PCR reactions with RNA from IMH grown in 1 mM As(V) condition. The numbering on the top of the gels corresponds to the product numbers drawn schematically in the outline given above. M, DNA mark; (+), positive control in which genomic DNA was used as template in the RT-PCR; RT, standard RT-PCR reaction; (—), negative control in which no reverse transcriptase was added to the RT reaction.
Herein, we present the first study of the molecular mechanism of As resistance in strain Pantoea sp. IMH. Two different ars systems - arsl (arsR1B1C1H1) and ars2 (arsR2B2C2H2) - were identified as being responsible for As resistance, contributing in different ways under changing temperature. In addition, we determined that the ars genes in IMH were probably acquired by HGT. The insights gained in this study improve our understanding of the flexible adaptation of microorganisms to resist As.
Results
As resistance systems in Pantoea sp. IMH. Strain Pantoea sp. IMH was able to resist up to 150 mM As(V) and 20 mM As(III), whereas E. coli W3110 with an arsRBC operon did not survive at concentrations above 50 mM As(V) and 5 mM As(III) (Fig. S1). To explore the molecular basis for its hyper-resistance to As, we determined the genome sequence of IMH and identified eight ars genes, including two arsR encoding a self-repressed transcriptional regulator, two arsB encoding a membrane-bound transporter that extrudes As(III) out of the cell, two arsC encoding a cytoplasmic As(V) reductase, and two arsH encoding an NADPH-dependent FMN reductase with an unknown biological function. These ars genes were organized as an ars1 cluster (arsR1B1C1H1) and ars2 cluster (arsR2B2C2H2) scattered on the chromosome (Fig. 1a). The genes in each ars cluster were separated by a short sequence of only a few nucleotides, suggesting they were organized in the same operon. To justify this hypothesis, we performed RT-PCR experiments using primers across intergenic regions (Table S3). The results indicate that the genes arsRBC within the arsRBCH cluster were organized as a co-transcribed operon, whereas arsH was in another operon (Fig. 1b).
The degree of DNA sequence identity between homologous genes underscores the appreciable differences between these two ars clusters. Specifically, as shown in Fig. S2, arsR1 and arsR2 shared 50% sequence identity, arsB1 and arsB2 shared 75%, arsC1 and arsC2 shared 60%, and arsH1 and arsH2 shared 70%. Moreover, two arsC-like genes with just 25% homology (arsC1 -like and arsC2-like) were found in the genome. It is an exceptional circumstance that IMH contains two ars systems and two As resistance molecular bases, considering that most bacteria have just one such cluster11. Therefore, we were motivated to investigate the functional contributions of each ars resistance system and molecular base to the As resistance in Pantoea sp. IMH.
Contribution of two ars systems and two arsC-like genes to As resistance. We first examined the transcription levels of ars genes in each ars operon and arsC-like gene by performing reverse transcription quantitative PCR (RT-Q-PCR), using 16 S rRNA as an internal control. As shown in Fig. 2, all genes of ars1 and ars2 clusters were completely transcribed. Notably, ars2 genes exhibited about 2-4 fold higher transcription levels than ars1 genes.
On the other hand, arsC1- and arsC2-like genes resulted in almost no expression (Fig. 2), indicating that these genes do not contribute to the As resistance. To justify this conclusion, we further analyzed the residues of ArsC in IMH. In contrast to previous observations that four residues of ArsC (Cys-12, Arg-60, Arg-94, and Arg-107) are required for As resistance12, Arg-60 and Arg-107 were not conserved in ArsC1-like and ArsC2-like proteins, respectively, in IMH (Fig. S3).
Figure 2. Transcription level of ars and arsC-like genes of strain IMH in the presence of As(V) with a concentration of 1 mM by RT-Q-PCR. The 16 S rRNA gene was used as an endogenous non-changing control. Data are shown as the means of three replicates, with the error bar illustrating one standard deviation.
Figure 3. As resistance capability of recombinant bacteria (E. coli-ars1 and E. coli-ars2) and ars cluster deleted strains (Aars1 and Aars2) at 30 °C and 15 °C. Growth of strains for 12 h in stationary phase in LB medium, tested with As(V) and As(III), E. coli W3110 and Pantoea sp. IMH were used as controls, respectively. IMH: wild type Pantoea sp. IMH; Aars1: ars1 cluster mutant strain; Aars2: ars2 cluster mutant strain; E. coli-ars1: E. coli W3110 with recombinant plasmid pUC18-ars1; E. coli-ars2: E. coli W3110 with recombinant plasmid pUC18-ars2. (a) Growth of E. coli-ars1 and E. coli-ars2 at 30 °C. (b) Growth of E. coli-ars1 and E. coli-ars2 at 15 °C. (c) Growth of Aars1 and Aars2 at 30 °C. (d) Growth of Aars1 and Aars2 at 15 °C. Data are shown as the means of three replicates, with the error bar illustrating one standard deviation.
To further identify the role of each ars cluster, we generated strains Aars1 and Aars2 lacking ars1 and ars2 clusters, respectively, based on the genome of IMH as described in the Methods section. The growth of Aars1, Aars2, and the wild type strain IMH was then monitored in the LB medium with 50 mM As(V) and 5 mM As(III). As shown in Fig. 3c, the growth of both Aars1 and Aars2 was substantially suppressed compared to the wild type strain IMH. Specifically, the deletion of ars1 resulted in more suppression of As resistance than that of ars2, implying that ars2 contributes to a greater extent to the overall As resistance. Moreover, we constructed the functional complementary plasmids pLGM1-ars1 and pLGM1-ars2 (see Methods for details) and then introduced them into the Aars1 and Aars2 strains, respectively. As shown in Fig. S4, the As resistance capabilities of complementary strains Aars1/pLMG1-ars1 and Aars2/pLMG1-ars2 were appreciably improved compared with the deleted strains (Aars1 and Aars2), which confirms the origin of As resistance in the corresponding operons.
Heterologous expression experiments were carried out to study the functional contributions of the two ars clusters to the As resistance. Recombinant plasmids for expression of ars1 and ars2 clusters were constructed according to the procedure described in the Methods. These plasmids were separately introduced in E. coli AW3110 (lacking any As resistance system), yielding the recombinant E. coli AW3110-ars1 and E. coli AW3110-ars2 strains. The growth of these strains together with the wild type strain E. coli W3110 (containing
Figure 4. Activity of Pars1 and Pars2 promoters with transcriptional lacZ fusions in E. coli AW3110. (a,b) The
organizations of Parsl and Pars2. -35 and —10 sequences are marked in green, transcription start sites (TSS) in red, RBSs in yellow, and start codons (ATG) in blue. (c) (3-galactosidase activity driven by Par1l and Pars2 promoters. Cultures were grown in LB medium at either 15 °C or 30 °C in the presence or absence of 1 mM As(III). Data are shown as the means of three replicates, with the error bar illustrating one standard deviation.
one ars operon) was then examined in LB medium containing 5 mM As(V) and 1 mM As(III). The results show that the heterologous host with the ar12 system acquired a higher resistance to As than that with the ar1l system (Fig. 3a).
In sum, the two ars clusters of strain IMH together contribute to its As resistance, in which the ars2 cluster is the major contributor. This observation raises a follow-up question: why did Pantoea sp. IMH evolve two ars systems to resist As? When bacteria species survive under changing environmental circumstances, some proteins may not function under all physicochemical conditions. Bacteria may meet this challenge by having two or more copies of genes to realize the same function under different conditions13. An example includes strain Pseudomonas putida KT2440, possessing two copies of equivalent ars operons to expand its functional scope14. Similarly, we propose that possessing two ars systems is an evolved strategy for IMH to survive in different ecological niches. To validate our speculation, the function of these two ars systems in different ecological conditions was further investigated.
Functioning of two ars systems under different environmental conditions. Major environmental factors, including concentrations of As(V) and As(III), pH, and temperature, may affect As resistance. We therefore tested the performance of each ars system in IMH at different pH values (pH 5, 7, and 9), As(V) concentrations (1 and 10 mM), As(III) concentrations (1 and 10 mM), and temperatures (15 and 30 °C). We first examined the transcriptional levels of arsCl and arsC2 using RT-Q-PCR as proxies for the expression of the ars systems. Figure S5 shows that the transcription level of arsC2 was higher than that of arsCl for different pH values and concentrations of As(V) and As(III). In contrast, an opposite result (i.e., arsCl > arsC2) was obtained at a lower temperature (15 °C). This observation suggests that pH, As concentration, and speciation have no influence on the major contribution of ars2 to As resistance, but a low environmental temperature enables arsl to be the predominant contributor.
Furthermore, the growth of E. coli AW3110-ars1, E. coli AW3110-ars2, Aars1, and Aars2 in LB medium with As at 15 °C was tested, where E. coli W3110 and IMH were used as controls. Figure 3b,d shows that E. coli-ars1 and Aars2 grew better than E. coli-ars2 and Aars1, respectively, at 15 °C, contrary to the results at 30 °C (Fig. 3a,c). In agreement with the RT-Q-PCR results, the above observations confirm that arsl contributed more than ars2 at the lower temperature (15 °C).
Why is this phenotype endowed with two ars systems dominant at different temperatures? We hypothesized that the activity of the arsR promoter should regulate its expression at different temperatures. To test this hypothesis, we first determined the transcription start site (TSS) using 5'-RACE, and predicted the —35 and —10 regions, as well as ribosomal binding site (RBS) sequence, of each of the promoter regions of arsR using SoftBerry software. As shown in Fig. 4a,b, the TSS was located 21 bp upstream of the translational start site of arsRl in arsl, and 136 bp in ars2. The distance between the —10 region and start codon was 28 bp in the arsRl promoter and
99 Raoultella planticola (plasmid pPA35) Klebsiella variicola DX120E Klebsiella pneumoniae JKC1 Hafnia alvei FBI
QO Escherichia coli KTE234
Enterobacter hormaechei ATCC 49162 —jRaoultella ornithinolytica 10-5246
100 Klebsiella oxytoca MGH41 —U" Panoea sp. IMH (Arsl)
100 Franconibacter helveticus LMG23732 Erwinia billingiae EB661 Erwinia gerundensis EM595 Enterobacter cloacae MGH 53 Cronobacter sakazakii ATCC BAA-894 I—P Serratia marcescens Dbl 1
Serratia marcescens (plasmid R478) Panoea sp. IMH (Ars2) Yersinia enterocolitica Y11
Citrobacter amalonaticus Y19 Salmonella enterica CFSAN000968 Citrobacter freundii CFNIH1
Pseudomonas putida KT2440 (Arsl) ~ Pseudomonas putida KT2440 (Ars2)
Figure 5. Neighbor-joining tree based on ArsRBCH protein sequences of IMH and representative microorganisms. The numbers at the nodes indicate levels of bootstrap support (%) based on analysis of 100 assembled datasets. Only values at or above 50% are given, bar 0.1 substitutions per amino acid position.
164 bp in the arsR2 promoter. The pronounced difference in the organization of these two arsR promoters may contribute to their different functions.
To examine the activity of the arsR promoter at different temperatures, we assembled equivalent reporter gene fusions between the predicted promoter regions of each ars cluster and a lacZ reporter gene without its promoter, producing plasmids pPR9TT-Pars1 and pPR9TT-Pars2 according to the procedure described in the Methods. The two constructed plasmids were transferred into E. coli AW3110, and their (-galactosidase levels were measured at 15 and 30 °C in LB media with and without As(III) (1 mM), which is the effective trigger for the ars operon. As shown in Fig. 4c, at the higher temperature (30 °C), the (3-galactosidase activity of pPR9TT-Pars2 was higher than that of pPR9TT-Pars1 with and without As(III) induction. In contrast, the activity of Pars2 was noticeably inhibited, while Pars1 showed comparatively high activity with and without As(III), even higher than that of Pars2 at the low temperature (15 °C). These results explain the different behaviors of ars1 and ars2 expression regarding their temperature dependence. The different performance of the two ars promoters was one of the important factors influencing the function of the two ars clusters in response to different temperatures. On this basis, we propose that the evolutionary reason for maintaining two ars systems in IMH is that their combination facilitates the survival of this strain over an extended range of temperatures in arsenic-polluted niches. This result raises a further question: what is the evolutionary origin of these two ars systems and two arsC-like genes?
Evolution of ars clusters and arsC-like genes in Pantoea sp. IMH genome. Substantial differences were observed between the homologous genes, ars1 and ars2, in their sequence identities and As resistance capabilities. We proposed that the differential origins of the two ars clusters may have derived from HGT. To validate our speculation, we first compared the genome of IMH to those of other bacterial strains including P. agglomerans Tx10, P ananatis LMG20103, P dispersa EGD AAK13, P rwandensis ND04, P stewartii DC283 and P. vagans C9-1. A summary of the features for each genome is shown in Table S4. Notably, ars genes exist only in Pantoea sp. IMH, P. agglomerans Tx10, and P. ananatis LMG, which belong to a pan-genome (Fig. S6). This result suggests that the ars systems in Pantoea strains may have been acquired by HGT.
Furthermore, we identified 25 IS elements in the genome of IMH (Table S5) using the IS finder database (http://www-is.biotoul.fr/), and found that some IS elements exist in the flanking region of ars clusters (Fig. S7). The IS elements are responsible for transferring genetic information between different cells15, and their abundance is positively correlated with HGT16.
Moreover, the variation of G + C content between clusters and the genome can be used as an indicator for HGT17. Compared to the G + C content of the IMH genome (54.74%), ars1 (56.28%) and ars2 (51.86%) clusters exhibited a great difference, supporting the HGT hypothesis.
In addition, the phylogenies constructed based on the sequences of ArsRBCH showed that ars1 and ars2 clusters of IMH have a sister-group relationship with the ars clusters of Franconibacter helveticus LMG23732 and Serratia marcescens (plasmid R478), respectively (Fig. 5). This result suggests that the ars1 cluster may originate from a common ancestor with F. helveticus LMG23732, and ars2 from S. marcescens (plasmid R478).
Phylogenies derived from each of the individual ArsR, B, C, and H were congruent with the phylogeny of the concatenated ArsRBCH (Fig. S8-11). Our result is in agreement with previous reports, where microorganisms can obtain the same functional genes from different sources18.
To evaluate the evolution of the two arsC-like genes, we constructed phylogenetic trees with ArsC and ArsC-like sequences. The result clearly shows that ArsC-like sequences were clustered together themselves and
divergent from the sequence of ArsC (Fig. S12). Homology between arsC and arsC-like genes showed that DNA sequence identity between arsCl and arsCl-like was 46%, and arsC2 and arsC2-like was 42%, suggesting that the arsC-like genes may have resulted from the duplication of arsC. Such a phenomenon is not rare in bacterial evolution19.
Discussion
Pantoea sp. IMH has been reported to have high resistance to As10. However, the molecular mechanism remains unknown. Our study determined that IMH has evolved two different As resistance systems, arsl and ars2 operons, by HGT. At higher temperature (30 °C), the ars2 operon plays the major role in As resistance, but its function is inhibited at lower temperature (15 °C), where the arsl system takes on most of the As resistance function. This kind of genetic constitution for As resistance in IMH is unique compared with other Pantoea strains whose genomes have been recently sequenced and annotated. More than one As resistance system in a strain can elevate its As resistance, which explains why IMH can survive in such high As concentrations.
Our genome sequencing showed that there are two arsC-like genes in the genome of IMH. Transcription levels of the two arsC-like genes were not detected when the strain encountered As (Fig. 2). This result suggests that arsC-like genes did not contribute to the As resistance. It is rare for arsC-like genes to show no As resistance capability. We further investigated the reasons for this phenomenon. It was reported that Cys-12, Arg-60, Arg-94, and Arg-107 were four conserved residues of the ArsC protein in the process of As resistance20. Cys-12 was identified as a catalytic residue and was activated by nearby residues Arg-60, Arg-94, and Arg-10712. Alignment analysis showed that Cys-12 and Arg-94 residues were conserved, but residues Arg-107 and Arg-60 in two ArsC-like proteins were not conserved respectively (Fig. S3). These changes in the amino acid sequence further verified that the two ArsC-like proteins did not contribute to As resistance. Interestingly, using phylogenetic analysis, we found that arsC-like sequences fell into distinct groups when compared to arsC genes. This suggests that multiple arsC-like genes may have resulted from arsC duplication and had already evolved with deviance. To clarify their relations with As resistance, the structures of ArsC-like proteins are worthy of further study.
It is thought that variants of a core arsRBC operon are common in the genomes of various bacteria, and it is rare that more than one ars operon appears in the same genome. Bacterial species usually adapt to changing environments by evolving two or more copies of genes, each one performing the same function under differential conditions13. Thus, the composition and gene distribution of a genome usually reflect the capacity for adaptation to different ecological niches14. In the IMH strain, we indeed found two ars systems with different patterns of expression and efficiency at different temperatures. This is a strategy for strain IMH to expand the scope of the encoded function to a wider range of physicochemical settings. Interestingly, our result was consistent with the report on Pseudomonas putida KT244014. We speculate that temperature is the most important environmental factor in the evolutionary history of ars clusters. Of course, if more strains with two or more copies of ars clusters are found, further mechanistic research should be carried out to support this hypothesis.
HGT is an important adaptation strategy to efficiently obtain 'alien' DNA21. To readily adapt to diverse and stringent growing conditions, IMH obtained two ars clusters by HGT. To identify the transfer of genetic information between genomes, we applied three commonly used methods including identification of IS elements and deviant G + C content, and phylogenetic analysis16'19. The existence of IS elements in the flanking region of ars clusters (Fig. S7) together with the greatly different G + C content in the IMH genome (54.74%), arsl (56.28%) and ars2 (51.86%) clusters, suggests that the two ars clusters (arsRlBlClHl and arsR2B2C2H2) may have been acquired by HGT. Phylogenetic analysis further revealed that the arsl cluster may have been acquired via HGT from a source related to Franconibacter helveticus LMG23732 in its early evolution, and the ars2 cluster from Serratia marcescens (plasmid R478). Our result is in agreement with previous reports where microorganisms can obtain the same functional genes from different sources18.
Methods
Genome sequencing, genome annotation and analysis. The genome of strain IMH was sequenced using the IlluminaHiSeq 2000 sequencing platform at the Beijing Genomics Institute (BGI) (Shenzhen, China). Genes were predicted from the assembled result using Glimmer 3.0222. The rRNA and tRNA genes were identified with RNAmmer and tRNAscan-SE23, respectively. Genome annotation was accomplished by analyzing protein sequences. The resulting translations were aligned with databases, including KEGG 5 9 24, GO 1.41925 and Swiss-Prot 20120626. The draft genome has been deposited in GenBank and the accession number used is JFGT01000000.
Strains, plasmids and culture conditions. The stains and plasmids used in this work are summarized in Table S2. E. coli and Pantoea strains were grown in LB medium (per liter contains: 10 g tryptone, 5 g yeast and 10 g NaCl) or LB plates (LB medium with w/v 1.5% agar) at either 15 °C or 30 °C as indicated in each case. When appropriate, antibiotics were added in the following concentrations: 100 |ig/mL ampicillin, 100 |ig/ mL kanamycin, and 100 |ig/mL streptomycin. For testing of minimal inhibitory concentrations (MICs), strains were incubated in LB medium with a series of concentrations of As(V) and As(III) as shown in Fig. S1. For monitoring the growth ofE. coli AW3110-arsl, E. coli AW3110-ars2 and E. coli AW3110, strains were cultured in LB medium with 5 mM As(V) and 1 mM As(III) in 96-microwell plates at either 15 °C or 30 °C, and OD600 was evaluated at 12 h. For monitoring the growth of Aars1, Aars2 and IMH, strains were cultured in LB medium with 50 mM As(V) and 5 mM As(III) in 96-microwell plates at either 15 °C or 30 °C, and OD600 was evaluated at 12 h. When detecting the As resistance under differential pH conditions, the pH of the LB medium was adjusted to pH 5.0, 7.0 and 9.0.
Construction of recombinant plasmids for expression in E. coli. In order to construct the plasmids used in the heterologous expression experiments, genomic DNA of Pantoea sp. IMH was used as a template for
cloning the two ars clusters. A 3.9 kb BamHI-XbaI DNA fragment containing the complete arsl cluster (promoter region, 360 bp upstream of the start codon ATG of arsR, the contiguous four genes arsRlBlClHl and 310 bp upstream of the start codon ATG of arsH) was PCR-amplified with primers Ars1-F and Ars1-R (Table S3). A 3.6 kb BamHI-XbaI DNA fragment containing the complete ars2 cluster (a 301 bp region downstream of the stop codon TAA of arsR2 and the contiguous ten genes arsR2B2C2H2 and 361 bp downstream of the stop codon TAA of arsH2) was PCR-amplified with primers Ars2-F and Ars2-R (Table S3). The PCR products were ligated to the BamHI-XbaI site of plasmid pUC18, yielding plasmids pUC18-ars1 and pUC18-ars2. Then the plasmids were transferred to E. coli AW3110, yielding the recombinant E. coli AW3110-ars1 and E. coli AW3110-ars2 strains, respectively.
Construction of Aars1 and Aars2. To obtain the deleted mutants of arsl and ars2 clusters in Pantoea sp. IMH, the suicide vector pARS10 was constructed by inserting the Invitrogen Gateway attR-CmR cassette into the backbone of SmaI-SmaI digested plasmid pKNG101, where E. coli DH5a (Xpir) was used as the host of pARS10. The Aars1 and Aars2 mutated stains were created using a modified Gateway method described by Choi27. To delete the arsl cluster, the flanking regions of the arsl cluster were amplified by PCR using primers as summarized in Table S3. A kanamycin resistance cassette derived from plasmid pKD4 was inserted between the flanking regions of the arsl cluster using a PCR overlap technique with the primers in Table S3. The resulting PCR products containing the Km-resistance cassette flanked by arsl cluster were cloned into the Gateway Entry vectors pDNOR221. The construct was transferred into the suicide vector pARS10, obtaining plasmid pARS1-1. The plasmid pARS1-1 was transferred into E. coli S17-1 and conjugally introduced into Pantoea sp. IMH. An allelic replacement event was selected based on double resistance. PCR with primers listed in Table S3 was used for the verification of the allelic replacement. Generation of the Aars2 strain followed the same method.
Construction of plasmids for complementation studies. In order to verify the As resistance function of the two ars clusters, plasmids pLGM1-ars1 and pLGM1-ars2 were constructed for complementation studies. As described in the previous section titled "Construction of recombinant plasmids for expression in E. coli", arsl and ars2 were amplified and then ligated to the BamHI- EcoRI site of plasmid pLGM1, yielding pLGM1-ars1 and pLGM1-ars2.
RT-PCR analysis. In order to determine the operons in ars clusters, an RT-PCR experiment with primers designed to span across intergenic regions (Table S3) was carried out. A culture of Pantoea sp. IMH was grown in LB medium with 1 mM As(V). After 8 h, the IMH strains were harvested by centrifugation at 4 °C, and the total RNA was isolated using the PrimeScript® RT reagent Kit with gDNA Eraser (Takara Bio) according to the manufacturer's instructions. The possibility of contamination of genomic DNA was eliminated by digestion with RNase-free DNase I (Takara Bio). The integrity and size distribution of the RNA were verified by agarose gel electrophoresis, and the concentration was determined spectrophotometrically. Synthesis of cDNA was carried out using RT Prime Mix according to the manufacturer's specification (Takara Bio). 1.0 |ig of cDNA was used for the template of RT-PCR.
RT-Q-PCR analysis. In order to understand the differences in each gene's transcription level in arsl and ars2 clusters under differential environmental factors, RT-Q-PCR analysis was used. Pantoea sp. IMH was grown in LB medium with different As(V) or As(III) concentrations (1 mM and 10 mM), in LB medium with 1 mM As(V) at different temperatures (15 °C and 30 °C), and in LB medium with different pH (5, 7 and 9). Then the cDNA was obtained as described in the RT-PCR analysis. Specific cDNA was employed to quantify the transcriptional signals of the ars genes and arsC-like genes, where 16 S rRNA gene was used as an internal reference. Primers used are listed in Table S3. RT-Q-PCR reactions were performed with three replicates using the ABI applied Biosystems vii A7.
Transcription start site identification. To determine the transcription start site (TSS) of the two ars operons, the 5'-RACE method was employed using the SMARTer™ RACE cDNA Amplification Kit (Clontech). Gene-specific primers are listed in Table S3. The PCR product was cloned into the pMD18-T Vector and then sequenced.
Pars-lacZ transcriptional fusions and |3-galactosidase assays. To explain the reason for the distinctly different performance of arsl and ars2 clusters at 15 °C and 30 °C, the promoter activities of the two ars clusters were determined. The promoter of the arsl cluster - a 110 bp DNA fragment (Parsl) (from —107 to +3 relative to the arsRl transcription start codon) - and the promoter of the ars2 cluster - a 241 bp DNA fragment (Pars2) (from —238 to +3 relative to the arsR2 transcription start codon) - were amplified from the total DNA of Pantoea sp. IMH using primers listed in Table S3. The amplified fragments were then cloned into the promoter vector pPR9TT, generating transcriptional fusions between the inserted promoter regions and a pro-moterless, complete lacZ gene, pPR9TT-Parsl and pPR9TT-Pars2. The plasmids were then transformed into E. coli DH5a, yielding E. coli DH5a/Parsl::lacZ and E. coli DH5a/Pars2::lacZ, respectively. For (-galactosidase activity assays, strains E. coli DH5a/Parsl::lacZ and E. coli DH5a/Pars2::lacZ were grown in LB medium with 1 mM As(III) for 12 h at 15 °C and 30 °C with shaking. (-galactosidase activity was measured according to the method described by Miller28. Briefly, a 100 |iL sample was mixed with 900 |iL Z buffer and shaken for 20 sec. Then, 200 |iL o-nitrophenyl-(-D-galactopyranoside (ONPG) (4 mg/mL) was added and incubated for 20 min at 30 °C. To stop the above reaction, 500 |iL of 1 M Na2CO3 solution was used. Finally, the OD420 and OD550 values were measured after the mixture was centrifuged. (-galactosidase (1 unit) = [1000 x (OD420 — 1.7 x OD550)]/[Time (min) x Vol (mL) x OD600].
Comparative genomics. All orthologous pairs between tests of Pantoea genomes were identified by Pan Genome Analysis Pipeline29. The common dataset of shared genes among test strains was defined as their core genome. The total set of genes with test genomes was defined as the pan-genome. The unshared genes in each strain were defined as unique genes. The genomes used in this study are listed in Table S4.
Phylogenetic analysis. In order to analyze the evolution of the two ars clusters of strain IMH, phylogenetic trees were constructed using the neighbor-joining method. Evolutionary distances were calculated according to Kimura's two-parameter model. Bootstrap analysis was performed on the basis of 1000 replications. The software package MEGA version 5.0 was used.
References
1. Oremland, R. S. & Stolz, J. F. The ecology of arsenic. Science 300, 939-944 (2003).
2. Páezespino, D., Tamames, J., De, L. V. & Cánovas, D. Microbial responses to environmental arsenic. Biometals 22, 117-130 (2009).
3. Silver, S. & Phung, L. T. Genes and enzymes involved in bacterial oxidation and reduction of inorganic arsenic. Appl. Environ. Microbiol. 71, 599-608 (2005).
4. Cuebas, M., Villafane, A., McBride, M., Yee, N. & Bini, E. Arsenate reduction and expression of multiple chromosomal ars operons in Geobacillus kaustophilus A1. Microbiology 157, 2004-2011 (2011).
5. Mukhopadhyay, R., Rosen, B. P., Phung, L. T. & Silver, S. Microbial arsenic: from geocycles to genes and enzymes. FEMS Microbiol. Rev. 26, 311-325 (2002).
6. Jackson, C. R. & Dugas, S. L. Phylogenetic analysis of bacterial and archaeal arsC gene sequences suggests an ancient, common origin for arsenate reductase. BMC Evol. Biol. 3, 18 (2003).
7. Fernandez, M., Udaondo, Z., Niqui, J. L., Duque, E. & Ramos, J. L. Synergic role of the two ars operons in arsenic tolerance in Pseudomonasputida KT2440. Environ. Microbiol. Rep. 6, 483-489 (2014).
8. Zhang, L. & Birch, R. G. The gene for albicidin detoxification from Pantoea dispersa encodes an esterase and attenuates pathogenicity of Xanthomonas albilineans to sugarcane. Proc. Natl. Acad. Sci. USA 94, 9984-9989 (1997).
9. Rezzonico, F., Smits, T. H., Montesinos, E., Frey, J. E. & Duffy, B. Genotypic comparison of Pantoea agglomerans plant and clinical strains. BMC Microbiol. 9, 204 (2009).
10. Wu, Q., Du, J., Zhuang, G. & Jing, C. Bacillus sp. SXB and Pantoea sp. IMH, aerobic As (V) -reducing bacteria isolated from arsenic-contaminated soil. J. Appl. Microbiol. 114, 713-721 (2013).
11. Shen, W. et al. Microbial responses to environmental arsenic. Biometals 22, 117-130 (2009).
12. Martin, P. et al. Insights into the structure, solvation, and mechanism of ArsC arsenate reductase, a novel arsenic detoxification enzyme. Structure 9, 1071-1081 (2001).
13. Sanchez-Perez, G., Mira, A., Nyiro, G., Pasic, L. & Rodriguez-Valera, F. Adapting to environmental changes using specialized paralogs. Trends Genet. 24, 154-158 (2008).
14. Páez-Espino, A. D., Durante-Rodríguez, G. & Lorenzo, V. Functional coexistence of twin arsenic resistance systems in Pseudomonas putida KT2440. Environ. Microbiol. 17, 229-238 (2015).
15. Touchon, M. & Rocha, E. P. Causes of insertion sequences abundance in prokaryotic genomes. Mol. Biol. Evol. 24, 969-981 (2007).
16. Francino, M. P. An adaptive radiation model for the origin of new gene functions. Nat. Genet. 37, 573-578 (2005).
17. Yan, Y. et al. Nitrogen fixation island and rhizosphere competence traits in the genome of root-associated Pseudomonas stutzeri A1501. Proc. Natl. Acad. Sci. USA 105, 7564-7569 (2008).
18. Arsene-Ploetze, F. et al. Structure, function, and evolution of the Thiomonas spp. genome. PLoS Genet. 6, e1000859 (2010).
19. Xie, J. B. et al. Comparative genomic analysis of N2-fixing and non-N2-fixing Paenibacillus spp.: organization, evolution and expression of the nitrogen fixation genes. PLoS Genet. 10, e1004231 (2014).
20. Gladysheva, T., Liu, J. & Rosen, B. P. His-8 lowers the pKa of the essential Cys-12 residue of the ArsC arsenate reductase of plasmid R773. J. Biol. Chem. 271, 33256-33260 (1996).
21. Zhaxybayeva, O. & Doolittle, W. F. Lateral gene transfer. Curr. Biol. 21, R242-R246 (2011).
22. Delcher, A. L., Bratke, K. A., Powers, E. C. & Salzberg, S. L. Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics 23, 673-679 (2007).
23. Lagesen, K. et al. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35, 3100-3108 (2007).
24. Kanehisa, M. et al. From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res. 34, D354-D357 (2006).
25. Ashburner, M. et al. Gene Ontology: tool for the unification of biology. Nat. Genet. 25, 25-29 (2000).
26. Magrane, M. & Consortium, U. UniProt Knowledgebase: a hub of integrated protein data. Database 2011, bar009 (2011).
27. Choi, K.-H. & Schweizer, H. P. An improved method for rapid generation of unmarked Pseudomonas aeruginosa deletion mutants. BMC Microbiol. 5, 30 (2005).
28. Miller, J. H. Experiments in molecular genetics. (Cold Spring Harbor Laboratory, 1972).
29. Zhao, Y. et al. PGAP: Pan-genomes analysis pipeline. Bioinformatics 28, 416-418 (2011).
Acknowledgements
We acknowledge the financial support of the National Basic Research Program of China (2015CB932003), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB14020302), and the National Natural Science Foundation of China (41373123, 41425016, 41503094 and 21321004). We thank Yongguan Zhu for the strain E. coli AW3110, and Sanfeng Chen for ppR9TT vector and strain E. coli S17-1.
Author Contributions
L.W. and C.J. conceived and designed the study. L.W. performed the laboratory work and data analysis, with assistance from X.Z. and G.Z., L.W. and C.J. drafted the tables and figures, and prepared the main manuscript.
Additional Information
Supplementary information accompanies this paper at http://www.nature.com/srep Competing financial interests: The authors declare no competing financial interests.
How to cite this article: Wang, L. et al. Arsenic resistance strategy in Pantoea sp. IMH: Organization, function and evolution of ars genes. Sci. Rep. 6, 39195; doi: 10.1038/srep39195 (2016).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
l/JJv ® I This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
© The Author(s) 2016
