Characterization and genomic analysis of chromate resistant and reducing Bacillus cereus strain SJ1
© He et al. 2010
Received: 15 April 2010
Accepted: 19 August 2010
Published: 19 August 2010
Skip to main content
© He et al. 2010
Received: 15 April 2010
Accepted: 19 August 2010
Published: 19 August 2010
Chromium is a toxic heavy metal, which primarily exists in two inorganic forms, Cr(VI) and Cr(III). Chromate [Cr(VI)] is carcinogenic, mutational, and teratogenic due to its strong oxidizing nature. Biotransformation of Cr(VI) to less-toxic Cr(III) by chromate-resistant and reducing bacteria has offered an ecological and economical option for chromate detoxification and bioremediation. However, knowledge of the genetic determinants for chromate resistance and reduction has been limited so far. Our main aim was to investigate chromate resistance and reduction by Bacillus cereus SJ1, and to further study the underlying mechanisms at the molecular level using the obtained genome sequence.
Bacillus cereus SJ1 isolated from chromium-contaminated wastewater of a metal electroplating factory displayed high Cr(VI) resistance with a minimal inhibitory concentration (MIC) of 30 mM when induced with Cr(VI). A complete bacterial reduction of 1 mM Cr(VI) was achieved within 57 h. By genome sequence analysis, a putative chromate transport operon, chrIA1, and two additional chrA genes encoding putative chromate transporters that likely confer chromate resistance were identified. Furthermore, we also found an azoreductase gene azoR and four nitroreductase genes nitR possibly involved in chromate reduction. Using reverse transcription PCR (RT-PCR) technology, it was shown that expression of adjacent genes chrA1 and chrI was induced in response to Cr(VI) but expression of the other two chromate transporter genes chrA2 and chrA3 was constitutive. In contrast, chromate reduction was constitutive in both phenotypic and gene expression analyses. The presence of a resolvase gene upstream of chrIA1, an arsenic resistance operon and a gene encoding Tn7-like transposition proteins ABBCCCD downstream of chrIA1 in B. cereus SJ1 implied the possibility of recent horizontal gene transfer.
Our results indicate that expression of the chromate transporter gene chrA1 was inducible by Cr(VI) and most likely regulated by the putative transcriptional regulator ChrI. The bacterial Cr(VI)-resistant level was also inducible. The presence of an adjacent arsenic resistance gene cluster nearby the chrIA1 suggested that strong selective pressure by chromium and arsenic could cause bacterial horizontal gene transfer. Such events may favor the survival and increase the resistance level of B. cereus SJ1.
The wide use of chromium (Cr) in textile, leather tanning and electroplating industries with subsequent sewage disposal causes severe contamination of global soil-water systems [1, 2]. Highly soluble, hexavalent chromium [chromate, CrO4 2-] is very toxic. As an analogue of sulfate, chromate can enter bacterial and mammalian cells readily via sulfate transport systems . The subsequent reduction of Cr(VI) by glutathione, thiols and other metabolites, and coproduction of reactive oxygen species (ROS) that damage DNA and other cellular components are the cause of the carcinogenic, mutational, and teratogenic potential of chromate [4–6]. On the other hand, the trivalent chromium [Cr(III)] is less bioavailable, thermodynamically stable and less toxic . Accordingly, the reduction of toxic Cr(VI) to stable Cr(III) is an efficient way to remove chromate from soil and water systems.
Bioremediation of chromate-contaminated sites, especially when stimulating indigenous microbial communities, is getting more and more attention because of its economical and environmental friendly aspects compared to chemical and physical methods [8–10]. An increasing number of Cr(VI)- reducing bacteria have been detected and studied including a pseudomonad strain CRB5 , Brucella sp. , Bacillus sp. strain QC1-2 , Burkholderia cepacia MCMB-821  and Thermus scotoductus strain SA-01 .
Bacteria have developed different strategies of chromate resistance including chromate efflux and chromate reduction. Efflux of chromate, which is mediated by the chromate transporter protein ChrA, has been confirmed in Pseudomonas aeruginosa [15, 16], Ochrobactrum tritici 5bvl1  and Shewanella sp. ANA-3 . Prior studies have not identified a chromate-responsive regulatory protein. Most chromate reduction studies have focused on soluble enzymes encoded by genes located on chromosomes . However, very few of the proteins responsible for chromate reduction have been purified and characterized because of technical difficulties. When examining induction of chromate resistance and reduction genes, several strains including Shewanella oneidensis MR-1 , Ochrobactrum tritici 5bvl1  and Ralstonia metallidurans strain CH34  have been shown to contain genes induced by chromate.
In this study, a chromate-resistant and reducing strain Bacillus cereus SJ1 was successfully isolated from chromium contaminated wastewater of a metal electroplating factory. Three chromate transporter related genes chrA, a chromate responsive regulator chrI, four nitR genes encoding nitroreductase and one azoreductase gene azoR possibly involved in chromate reduction were identified by the draft genome sequence. Using RT-PCR technology, we found that all of the five genes encoding putative chromate reductases appeared to be expressed constitutively. In contrast, the gene chrA1 encoding a transporter with high homology to other transporters linked to chromate resistance was up-regulated by the addition of Cr(VI) together with the adjacent putative transcriptional regulator chrI. Since chrA1 is probably regulated by chrI, this suggests identification of the first known chromate-responsive regulator.
Strain SJ1 showing both high Cr(VI) resistance and reduction abilities was isolated from industrial wastewater of a metal plating factory. SJ1 was a Gram positive, rod shaped bacterium. The 16 S rDNA sequence was used for bacterial identification. SJ1 showed the highest identity (100%) with B. cereus 03BB102 [GenBank: CP001407] and was hereafter referred to as B. cereus SJ1.
Draft genome sequence analysis of B. cereus SJ1 showed a genome size of about 5.2 Mb distributed in 268 contigs with an average GC content of 35.4%, containing 5,708 putative coding sequences (CDSs). There are 100 tRNA genes representing all 20 amino acids and 6 scattered ribosomal RNA genes identified on the draft genome. The likely origin of replication of the chromosome of B. cereus SJ1 was located in a 9.0 kb region that included co-localization of six genes (rpmH, gyrA, gyrB, recF, dnaN and dnaA). It was localized by comparing its draft genome to complete genomes of several strains of the B. cereus group though MUMmer3.20.
Since B. cereus SJ1 was isolated from industrial wastewater containing various toxic elements in addition to chromium, the MICs of B. cereus SJ1 for these heavy metals were determined. For B. cereus SJ1, the highest resistance was found for As(V), while Hg(II) was the most toxic compared to the other metal ions. When B. cereus SJ1 was incubated with increasing As concentration, no viable cells were recovered at concentrations above 50 mM As(V) and 4 mM As(III). The MICs of B. cereus SJ1 for Cu(II), Co(II), Ni(II), Cd(II), Ag(I) and Hg(II) were 0.9 mM, 0.8 mM, 0.7 mM, 0.2 mM, 0.02 mM and 0.007 mM, respectively. In order to survive in such unfavorable habitat, B. cereus SJ1 must have various determinants to tolerate such harsh conditions. For example, the copper concentration of the wastewater was as high as 0.65 mM and the MIC of B. cereus SJ1 to copper was 0.9 mM in R2A medium. When we analyzed the genome sequence of B. cereus SJ1, several genes related to copper resistance including copper-exporting P-type ATPase CopA, copper export protein CopC, copper resistance protein CopD, copper homeostasis protein CutC and two multicopper oxidases were identified. Furthermore, many other putative heavy metal resistance genes including those for As, Zn, Mn, Co, Cd, Te and Hg were also identified in the B. cereus SJ1 draft genome (Additional file 2).
Using the procedure described in Methods, we found that the uninduced and induced cells grew to similar cell densities in medium containing 5 mM Cr(VI) as determined spectrophotometrically at OD600. However, the induced cells grew to higher cell densities than the uninduced cells at higher Cr(VI) concentrations in the growth medium. The MIC of induced B. cereus SJ1 to K2CrO4 was 30 mM whereas that of the uninduced strain was 20 mM (Figure 4B).
The chrI gene located upstream of chrA1 encodes a protein with 98% amino acid sequence identity to the PadR-family transcriptional regulator from B. thuringiensis serovar konkukian str. 97-27 [GenBank: YP036529]. As chrI was a potential transcriptional regurator, it should be responsive to the inducer (Cr), so we analyzed the transcription of chrI at 5 and 15 min after addition of K2CrO4. A very weak PCR product was detected with cDNA from uninduced cells as shown in Figure 6B. The level of the chrI gene transcript was 16-fold higher (analyzed using BandScan 5.0 program) in cells induced for 15 min compared to the uninduced culture (lane 4 vs 6), confirming substrate-mediated regulation of chrI. To confirm the hypothesis that chrI-chrA1 was transcribed as a single transcription unit, RT-PCR was carried out with mRNA prepared from B. cereus SJ1 grown with and without K2CrO4 (0.3 mM) as described above. PCR products of the expected size (1,130 bp) were obtained with cDNA from both induced and uninduced cultures as the templates (Figure 6B), which indicated chrI and chrA1 were arranged as an operon. No PCR products were amplified using total RNA as the template that was designed to detect DNA contamination.
The arrangement of chrI genes in an operon together with chrA encoding a chromate transporter can be detected in both Gram positive and Gram negative bacteria (Additional file 3). An alignment of ChrI homologs was constructed using ChrI of B. cereus SJ1 and other related proteins encoded in operons having a chrI gene adjacent to a chrA gene (Additional file 4). The more-conserved domains were located in the N- and C-terminal regions. Within the conserved domains, two amino acids, lysine and arginine, were identified that might be involved in chromate binding and recognition.
Chromate-reducing bacteria have been discovered in both contaminated and non-polluted environments [1, 13, 24, 25]. In this study, a chromate-resistant strain B. cereus SJ1 was isolated from chromium contaminated wastewater of a metal plating factory in China. B. cereus SJ1 showed a rapid growth rate in chromate containing medium and efficient chromate-reducing ability under aerobic conditions. Since the isolation site for B. cereus SJ1 was contaminated with as much as 1.89 mg Cr per liter (36.28 μM), we reasoned that genes conferring chromate resistance could be present in this strain. By genome sequencing analysis, we were able to identify a putative chromate transport operon chrIA1 and two additional chrA genes encoding putative chromate transporters that confer chromate resistance. Furthermore, we also found an azoreductase gene azoR and four nitR genes that encode nitroreductases which may catalyze reduction of chromate [19, 23].
The membrane transporter protein ChrA has been shown to be responsible for extrusion of chromate ions across the cytoplasmic membrane in P. aeruginosa [15, 16], Ochrobactrum tritici 5bvl1  and Shewanella sp. ANA3 . It was demonstrated that the chromate transporter ChrA functions as a chemiosmotic pump that extrudes chromate using proton-motive force . ChrA protein belongs to the CHR superfamily which includes dozens of putative homologs from all three domains of life . Cr(VI) induction of B. cereus SJ1 in this study conferred the ability to survive at a higher chromate concentration. Exposure to chromate resulted in the up-regulation of chrA1 and higher chromate resistance. Possibly increased level of ChrA1 is responsible for higher chromate resistance.
The chrI gene product located upstream of chrA1 showed a high homology to PadR-family transcriptional regulators. The padA gene encoding phenolic acid decarboxylase, is a member of the PadR family that has been identified as a transcriptional repressor in Pediococcus pentosaceus  and Lactobacillus plantarum . Although genes encoding PadR homologs located either upstream or downstream of putative chromate transporter gene chrA have been identified in many genera, such as B. thuringiensis serovar konkukian str. 97-27 [GenBank: YP036529], Oceanobacillus iheyensis HTE831 [GenBank: NP694199], B. licheniformis ATCC 14580 [GenBank: YP093604) and Alkaliphilus oremlandii OhILAs [GenBank: YP001512811], the real function of a PadR homolog associated with chromate resistance has never been reported. In this study, this gene encoding a PadR homolog was renamed as chrI since it was induced by chromate. By an alignment of most PadR-like regulators which form an operon with the chromate transporter gene chrA, highly conserved basic amino acids (lysine and arginine) were identified in ChrI and the homologs that might be involved in chromate binding and recognition because they would carry a positive charge under physiological conditions. Possibly the negatively charged oxyanion CrO4 2- would preferentially bind the basic, positively charged amino acids conserved in the putative transcriptional regulator ChrI.
A strong selective pressure for transformation of metal- and metalloid-related resistance genes is present in heavy metal contaminated environments [29, 30]. Horizontal gene transfer (HGT) events driven by mobile genetic elements, such as phages, plasmids, insertion sequences, integrons and transposons, have been shown to provide microbes with a wide variety of adaptive traits for microbial survival under hostile environmental conditions. In this study, B. cereus SJ1 was isolated from wastewater contaminated with multiple heavy metals. The presence of a resolvase gene upstream of chrIA1 and an arsenic resistance gene cluster, and the Tn7-like transposition proteins ABBCCCD gene downstream of the arsenic resistance operon in B. cereus SJ1 but absent in other strains of B. cereus implied the possibility of a recent HGT event. Interestingly, other strains of B. cereus harbor a gene encoding CHRD-domain-containing protein adjacent to the chrA gene. Whether these proteins have a regulatory role is currently unknown . In addition, ChrA1 from B. cereus SJ1 is only distantly related to ChrA proteins from other strains of B. cereus indicating potential horizontal gene transfer from other Gram-positive bacteria as an adaptation to survive in a highly chromate contaminated environment.
Chromate can be reduced nonenzymatically as well as by various bacterial enzymes. Dihydrolipoamide dehydrogenase from Thermus scotoductus SA-01 , azoreductase in Shewanella oneidensis  and flavoproteins from P. putida and E. coli  were previously reported to be associated with Cr(VI) reduction. Compared to the one electron transfer chromate reductase gene chrR from P. putida, yieF from E. coli was proposed to be a more appropriate gene for bioremediation applications because of the three-electron transfer ability of its gene product and consequently, the generation of fewer reactive oxygen species (ROS) . In our study, one azoreductase gene azoR and four nitR genes encoding nitroreductase obtained from B. cereus SJ1 showed high identities with other Cr(VI) reductases and were expressed constitutively. Since Cr(VI) reduction of strain SJ1 was not inducible by chromate, other potential chromate reductases in B. cereus SJ1 must also be constitutively expressed and the enzyme activity is probably adventitious.
This study describes insights into the chromate resistance and reduction capabilities of B. cereus SJ1 using both physiological and molecular techniques. The expression of the chromate transporter gene chrA1 was inducible by Cr(VI) and most likely regulated by chrI. Even though the physiological function of ChrI has not been verified due to the absence of a genetic system for this Gram positive strain, ChrI is most likely the first identified chromate responsive regulator. In addition, genome analysis identified a number of putative genes encoding gene products with possible functions in chromate resistance and reduction which may be the basis for the observed high chromate resistance and reduction ability of this strain. Furthermore, possible horizontal gene transfer events indicated in this study may have enabled B. cereus SJ1 to survive in metal (loid) contaminated environments.
Industrial wastewater samples were obtained from a metal electroplating factory in Guangdong, China. The total concentrations of Cr, Cu, Zn, Mn, Pb, Co, As and Cd in this sample determined by atomic absorption spectrometry were 36.28 μM, 0.65 mM, 24.88 μM, 7.83 μM, 0.49 μM, 0.41 μM, 0.32 μM, and 0.007 μM, respectively. Isolation of chromate-resistant and reducing bacteria was performed as described . The abilities of the chromate-resistant bacteria to reduce Cr(VI) (K2CrO4) were determined using a spectrophotometric method using the reagent 1, 5-diphenylcarbazide (DPC) . Several chromate-resistant bacteria were isolated and strain SJ1 was chosen for this study. The 16 S rDNA of strain SJ1 was obtained from the genome sequence (see below) and analyzed by BlastN searching tools http://www.ncbi.nlm.nih.gov/blast. Cell morphologies were examined under a scanning electron microscope (SEM; JSM-6390, JEOL, Japan) with 20,000 V accelerating voltage and 15,000 times enlargement.
The MIC, defined as the lowest concentration of heavy metals that inhibited growth in R2A broth (Becton Dickinson, MD, USA), was performed with strain SJ1. A 1% inoculum of an overnight culture was introduced into R2A medium amended with different concentrations of CuCl2, NiCl2, Co(NO3)2, Na2HAsO4, NaAsO2, HgCl2, CdCl2 and AgNO3, incubated at 37°C on a rotary shaker at 200 rpm for 3 days. MIC values were determined spectrophotometrically at OD600.
The exponential phase cultures of uninduced, and induced with 1 mM K2Cr2O6 for 8 h, were adjusted to the same OD600. One hundred microliters of each culture was added to 10 ml fresh LB medium with increasing amounts of K2CrO4, and incubated at 37°C with 200 rpm shaking for 3 days. The OD600 values were then determined spectrophotometrically. For chromate reduction, the uninduced and induced cultures were prepared as above and inoculated into 100 ml LB medium amended with 1 mM K2CrO4 and incubated at 37°C on a rotary shaker at 200 rpm for about 60 h. The residual Cr(VI) concentration was monitored as described above. LB medium with 1 mM K2CrO4 without bacterial cells was incubated as a negative control to monitor abiotic chromate reduction.
High-molecular-mass genomic DNA isolated from B. cereus SJ1 using Blood & Cell Culture DNA Mini Kit (Qiagen, MD, USA) was used to construct a 4 kb to 40 kb random genomic library. Whole genome shotgun sequencing was performed by the University of Arizona Genetics Core facility, using a Roche 454 Genome Sequencer FLX instrument. The B. cereus SJ1 DNA sample was loaded onto one region of a standard four-region plate. A local Linux computing cluster was used for signal processing on the images produced by the FLX instrument. The Roche gsassembler software version 2.0.01 was used for de novo assembly of the 271,408 reads. Using the default assembly parameters, 141 contigs of length greater than 500 bp were built, along with 127 shorter contigs. These 268 contigs were submitted to the RAST annotation server  for subsystem classification and functional annotation. Genome comparison was performed through SEED server http://rast.nmpdr.org/seedviewer.cgi and the subsequent results were modified manually. GC content was analyzed using CLC Main Workbench 5 program http://www.clcbio.com. The NCBI Prokaryotic Genomes Automatic Annotation Pipeline was used for gene annotation in preparation for data submission to GenBank. http://www.ncbi.nlm.nih.gov/genomes/static/Pipeline.html.
RT-PCR was used to assess induced expression and co-transcription of the chromate resistance and reduction related genes of strain SJ1. Total RNA was obtained from mid-exponential phase strain SJ1 cells grown from 0 h to 3 h in the presence or absence of 0.3 mM K2CrO4 in LB medium. Total RNA was isolated by the RNeasy Mini Kit (Qiagen) and then digested with DNase I (Fermentas, MD, USA) to remove any DNA. The OD260 values were then determined spectrophotometrically for the total RNA concentration. Equal amounts of total RNA were used to perform cDNA synthesis using iScript™Select cDNA Synthesis Kit (Biorad, CA, USA). Standard PCR programs were used to generate amplicons from 3 μl of the reverse transcription reaction mixture using the specific primer pairs listed in Additional file 5. PCR amplification using RNA as template was served as the control to investigate the potential presence of DNA contamination. The relative levels of the cDNAs of RT-PCR were determined by densitometric analyses using BandScan 5.0 software (GLyko Inc., Novato, CA, USA) using 16 S rRNA genes as references.
B. cereus SJ1 was deposited in The Agricultural Research Service Culture Collection, USA (NRRL http://nrrl.ncaur.usda.gov) under the accession number of NRRL B-59452. The Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank http://www.ncbi.nlm.nih.gov/sites/genome under the accession number of ADFM00000000. The version described in this paper is the first version, ADFM01000000.
MH is supported by the exchanging PhD student scholarship of the Ministry of Education, China. This work is funded by the National Natural Science Foundation of China (30970075).
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.