Original Research ARTICLE
Genomic and Transcriptomic Insights into Calcium Carbonate Biomineralization by Marine Actinobacterium Brevibacterium linens BS258
- 1Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
- 2Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology, Qingdao, China
- 3College of Earth Science, University of Chinese Academy of Sciences, Beijing, China
- 4College of Biotechnology and Bioengineering, Zhejiang University of Technology, Hangzhou, China
- 5College of Life Sciences, Qingdao University, Qingdao, China
Calcium carbonate (CaCO3) biomineralization has been investigated due to its wide range of scientific and technological implications, however, the molecular mechanisms of this important geomicrobiological process are largely unknown. Here, a urease-positive marine actinobacterium Brevibacterium linens BS258 was demonstrated to effectively form CaCO3 precipitates. Surprisingly, this bacterium could also dissolve the formed CaCO3 with the increase of the Ca2+ concentration. To disclose the mechanisms of biomineralization, the genome of B. linens BS258 was further completely sequenced. Interestingly, the expression of three carbonic anhydrases was significantly up-regulated along with the increase of Ca2+ concentration and the extent of calcite dissolution. Moreover, transcriptome analyses revealed that increasing concentration of Ca2+ induced KEGG pathways including quorum sensing (QS) in B. linens BS258. Notably, most up-regulated genes related to QS were found to encode peptide/nickel ABC transporters, which suggested that nickel uptake and its associated urease stimulation were essential to boost CaCO3 biomineralization. Within the genome of B. linens BS258, there are both cadmium and lead resistance gene clusters. Therefore, the sequestration abilities of Cd2+ and Pb2+ by B. linens BS258 were checked. Consistently, Pb2+ and Cd2+ could be effectively sequestered with the precipitation of calcite by B. linens BS258. To our knowledge, this is the first study investigating the microbial CaCO3 biomineralization from both genomic and transcriptomic insights, which paves the way to disclose the relationships among bacterial metabolisms and the biomineralization.
Carbonates account for 41.9% of the total carbon on Earth and represent an important carbon reservoir (Ehrlich, 1996). At the Earth’s surface, a large fraction of the insoluble carbonate is biogenic and involves bacteria, fungi, algae, and metazoan (Gadd, 2010). A significant portion of carbon exported from the sea surface to deep waters and sediments is in the form of biogenic CaCO3 (Berelson et al., 2007). In a geological perspective, the importance of biogenic calcification is evident from the global existence of calcium carbonate (CaCO3) rocks including limestone and dolomite (Heldal et al., 2012). Therefore, biogenic production and sedimentation of calcite in the ocean have profound implications for the ocean carbon cycle and ocean acidification (Heldal et al., 2012).
Biomineralization of CaCO3 is a phenomenon of biologically controlled or induced mineralization of CaCO3, and it is also an important process resulting in the production of biominerals by a range of taxa including bacteria, archaea and eukarya (Bundeleva et al., 2012; Rodriguez-Navarro et al., 2012; Benzerara et al., 2014). Generally, biomineralization is classified as either biologically controlled mineralization (BCM) or biologically induced mineralization (BIM) (Marvasi et al., 2010). BCM of CaCO3 is controlled through intracellular organic matrices or vesicles, which usually takes place in eukaryotes, such as molluscs, sea urchin spines, and fish otoliths (Lowenstam and Weiner, 1989). In contrast, BIM of CaCO3 does not involve direct control of the biomineralization process by the organisms, and it is regulated by physiological activities and carried out in open environments (Lowenstam and Weiner, 1989; Ries et al., 2008; Rodriguez-Navarro et al., 2012). BIM occurs either passively, due to metabolically driven changes in the chemistry of the bulk solution or around the living organisms, or actively, when the organism and/or its metabolic byproducts provide nucleation sites for mineralization (Mitchell and Ferris, 2006; Rodriguez-Navarro et al., 2012). BIM of CaCO3 typically occurs in the presence of unicellular organisms such as bacteria and is often correlated with metabolic pathways such as photosynthesis, urea hydrolysis and sulfate reduction (Dupraz and Visscher, 2005; Braissant et al., 2007; Marvasi et al., 2010). However, the precipitation of CaCO3 by bacteria through urea hydrolysis is the most widely used method (Anbu et al., 2016).
Urease-based microbially induced carbonate precipitation (MICP) is brought about through the ureolytic pathway, which leads to alkalinization of the surrounding environment (Stocks-Fischer et al., 1999). In addition to 2 M of NH4+, 1 M of carbonate ion is formed (reaction 1). Thus, in the presence of Ca2+ (reaction 2), calcium carbonate (CaCO3) is produced (Kang et al., 2014).
Calcium carbonate precipitation by many bacterial genera has been observed in different natural environments including seawater (Wright and Oren, 2005). The ecological implications of MICP are various (Wright and Oren, 2005), and the ecological relevance of MICP has been exploited in biotechnology, with applications in the bioremediation of metal contaminated soil and groundwater (Mitchell et al., 2010; Li et al., 2015).
It is generally accepted that MICP is correlated to metabolic activities and cell surface structures (Barabesi et al., 2007). A plausible link of CaCO3 production to Bacillus subtilis genes has been reported, and this research suggested that bacterial precipitation of these minerals brings some evolutionary advantage to the microorganisms (Barabesi et al., 2007). For instance, CaCO3 precipitation may help to fixate excess Ca2+ (Obst et al., 2009) or may help bacteria to survive unfavorable conditions in a cryptobiotic state (Yoshida et al., 2010). Altogether, these studies are helpful to explain why bacterial CaCO3 precipitation appears to be a ubiquitous phenomenon. However, the actual roles and corresponding molecular mechanisms especially from genomic and transcriptomic insights of bacteria in carbonate mineralization are still debated.
In this study, a urease positive marine actinobacterium B. linens BS258 was demonstrated to form and dissolve calcite precipitation in conditions at different concentration of Ca2+. Genome sequencing, transcriptome profiling and other assays of B. linens BS258 suggested that calcite precipitation and dissolution were linked to carbonic anhydrase (CA) activities, quorum sensing (QS) and other essential energy metabolisms. Heavy metal resistance and removal abilities related to MICP by B. linens BS258 were further dissected. Thus, the results of this study provided novel understanding of the mechanisms involved in bacterial carbonate biomineralization and proved its potentials for metal bioremediation in the marine environment.
Materials and Methods
Isolation and Identification of Bacterial Strains Producing Calcite
Sediments were collected from the China Yellow Sea (118°25′E, 35°56′N) and maintained in the laboratory at 4°C. The marine bacterial strains used in this study were isolated from the sediments with dilution method described as previously (Schut et al., 1993), and cultured in marine broth 2216E (5 g/L tryptone, 1 g/L yeast extract, one liter filtered seawater, 15 g agar, pH adjusted to 7.4–7.6) at 28°C. Urease-positive strains were preliminary screened with 9219 medium containing peptone (1 g), KH2PO4 (2 g), NaCl (5 g), urea (20 g), d-glucose (1 g), phenol red (0.012 g) and agar (15 g) in one liter filtered seawater, pH 6.8 ± 0.2 (Balan et al., 2012), and the plates would turn red for urease-positive strains when incubated at 28°C. CaCO3-producing strains were then screened from these urease-positive strains using CaCO3-producing medium, after incubation at 28°C for 5 days, the bacterial colonies were checked by microscope to confirm the production of calcite. The CaCO3-producing medium was prepared as following: nutrient broth (3 g), NaHCO3 (2.12 g), NH4Cl (10 g), CaCl2 (16.65 g), urea (20 g) and agar (15 g) in one liter filtered seawater. Genomic DNA was extracted from the isolate, and PCR (polymerase chain reaction) was performed to amplify the 16S rDNA gene sequence as described previously with universal primers 27F (5′-AGAGTTTGATCCTGGCTCAG-3′) and 1492R (5′-AAGGAGGTGATCCACCC-3′) (Wu et al., 2016). To determine the phylogenetic position of the B. linens BS258, the 16S rDNA gene sequence was analyzed by the BLAST programs1. The phylogenetic tree was constructed with MEGA6.0 (Tamura et al., 2013).
Brevibacterium linens BS258 was firstly incubated in 2216E medium to a logarithmic phase, then 1 mL starter was inoculated in 100 mL 9219 base medium for the production of calcite. The 9219 base medium was prepared as following: 10 g peptone, 5 g yeast extract, 3.33 g CaCl2, 2 g CH3COONa, 2 g KH2PO4, 0.005 g NiSO4, 0.005 g MnSO4, 5 g NaCl, 5 g urea and 20 g d-glucose in one liter filtered seawater. After incubating at 28°C for 10 days, the bacterial cultures were centrifuged at 30 g for 4 min (Riding, 2000). The collected mineral precipitate was then washed with Milli-Q water for three times and dried at 60°C for 24 h. Dried minerals were then further analyzed by Scanning Electron Microscopy (SEM), Energy-Dispersive X-ray (EDX), Transmission Electron Microscopy (TEM) and Fourier Transform Infrared Spectroscopy (FTIR).
SEM, TEM, and EDX Analyses of the Minerals Formed by B. linens BS258
Minerals precipitated by B. linens BS258 were mounted on double-sided carbon adhesive tape on aluminum stubs prior to examination by SEM and EDX. For SEM, particulate materials mounted on stubs were sputter coated for 5 min with gold and platinum (10 nm) using a Hitachi MC1000 Ion Sputter (Li et al., 2015). Specimens were examined using a field emission scanning electron microscope (Hitachi S-3400N) operating at an accelerating voltage of 5 kV. For EDX, uncoated samples were used, and operation was at an accelerating voltage of 15 keV for 100 s (Li et al., 2015).
Aliquots of suspensions containing untreated/treated live bacterial cells and precipitated minerals were also examined using TEM with a JEOL JEM 12000 EX (equipped with a field emission gun) at 100 kV. The cell suspension was rinsed using sterile nutrient solution or Milli-Q water, centrifuged 2 min at 7000 × g (Bundeleva et al., 2014). Samples were taken by immersing grids coated with a carbon film for 10 s in prepared bacterial suspension and dried several minutes at room temperature. The possible precipitated carbonate observed with TEM was further confirmed using Selected Area Electron Diffraction (SAED) (Han et al., 2014).
Calcite Precipitation and Dissolution by B. linens BS258 under Different Conditions
To further investigate the calcite precipitation by B. linens BS258 in the plate, 5 μL overnight cultured cell suspension of B. linens BS258 was loaded in the center of each plate with CaCO3 precipitation medium supplemented with different extra concentration of CaCl2 (0 mM, 50 mM, 150 mM, and 500 mM). The plates were incubated at 28°C for 5 days. To study the dynamics of calcite precipitation and dissolution by B. linens BS258, the plates supplemented with additional 150 mM CaCl2 in the CaCO3 precipitation medium were monitored with light microscopy every 12 h. Meanwhile, the urease activities of B. linens BS258 at the same time points were checked as described above.
Changes in pH and Dissolved Calcium Concentrations Over Time during Growth of B. linens BS258 under Different Conditions
Brevibacterium linens BS258 was firstly incubated in 2216E medium to a logarithmic phase, then 1 mL starter was inoculated to 50 mL modified Luria-Bertani (LB) medium (10 g/L peptone, 5 g/L yeast extract, 10 g/L NaCl and 20 g/L urea) prepared with Milli-Q water supplemented with 50, 150, and 500 mM CaCl2 for the production of calcite. After 5 days of incubation, the pH and concentration of Ca2+ in the medium were measured, respectively. The supernatant for measuring pH and the concentration of dissolved calcium was collected by centrifugation at 7000 × g for 2 min. And the supernatant was thoroughly digested with nitric acid and perchloric acid, and diluted with Milli-Q water for calcium concentration detection (Li et al., 2013). The dissolved calcium concentrations were measured with an inductively coupled plasma source mass spectrometer (Optima 7300 DV, PerkinElmer).
Quantitative Real-time PCR (qRT-PCR) Analyses of Carbonic Anhydrases Expression of B. linens BS258 under Different Conditions
For qRT-PCR assay, the cells of B. linens BS258 incubated in modified LB medium amended with three different concentrations of Ca2+ (with additional 0, 50, 150 mM Ca2+) were collected respectively after 6 h incubation. Total RNAs were extracted using the RNApure Bacteria Kit (DNase I) (CWBio, China). Then the total RNAs were reverse transcribed into cDNA and the transcriptional levels of three genes of CAs were determined by qRT-PCR using Sybr Green Premix Low rox (MDbio, China) and the QuantStudioTM 6 Flex (Thermo Fisher Scientific, America). RNA degradation and contamination was monitored on 1% agarose gels. RNA purity was checked using the NanoPhotometer® spectrophotometer (IMPLEN, Westlake Village, CA, USA). RNA concentration was measured using Qubit® RNA Assay Kit in Qubit® 2.0 Flurometer (Life Technologies, Carlsbad, CA, USA). RNA integrity was assessed using the RNA Nano6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA). 16S rDNA was used as an internal reference and the relative gene expression was calculated using the 2-ΔΔCt method (Qin et al., 2015). Specific primers for CAs and 16S rDNA were designed using Primer 6.0 as shown in Supplementary Table S1. All qRT-PCR runs were conducted with three biological and three technical replicates. Specific primers for CAs and 16S rDNA were designed using Primer 6.0 as shown in Supplementary Table S1.
Transcriptional Profiling of B. linens BS258 Challenged with Different Concentrations of Ca2+
Total RNAs of B. linens BS258 incubated in the modified LB medium with three different concentrations of Ca2+ (with additional 0, 50, 150 mM CaCl2) after 24 h of incubation were extracted and checked as described above. The detailed protocols of RNA-seq were described in the Supplementary Information. qRT-PCR verification of the data from transcriptional profiling was performed as described above. Five genes (A2T55_RS10820, A2T55_RS11225, A2T55_RS11230, A2T55_RS01425, A2T55_RS01430) were chosen for verification of RNA-Seq data by qRT-PCR. Total RNAs were extracted and checked as described above. Then the total RNAs were reversely transcribed into cDNA and the transcriptional levels of 5 selected genes were determined by qRT-PCR using Sybr Green Premix Low rox (MDbio, China) and the QuantStudio 6 Flex (Thermo Fisher Scientific, America). 16S rDNA was used as an internal reference and the relative gene expression was calculated using the 2-ΔΔCt method (Qin et al., 2015). Specific primers were designed using Primer 6.0 as shown in Supplementary Table S3.
Co-precipitation of Ca2+, Cd2+ and Pb2+ by B. linens BS258
A modified LB medium (10 g/L peptone, 5 g/L yeast, 10 g/L NaCl, 80 mM urea and 30 mM CaCl2) supplemented with different metal ions was used to investigate the remediation ability of B. linens BS258. Briefly, B. linens BS258 was incubated for 24 h with shaking at 28°C in LB medium supplemented with different concentration of CdCl2 (0.025, 0.05, 0.1, 0.13, and 0.16 mM), or different concentration of PbCl2 (1, 2, 3, 4, and 5 mM). To measure the concentrations of Cd2+ and Pb2+ at the end of the incubation, the supernatant was collected by centrifugation (13400 × g, 2 min), and was digested completely and diluted with Milli-Q water. The concentrations of Cd2+ and Pb2+ in the solution were measured with an inductively coupled plasma source mass spectrometer (Optima 7300 DV, PerkinElmer) (Han et al., 2014).
The significant differences among groups were subjected to one-way analysis of variance (one-way ANOVA) and multiple comparisons by using the SPSS 18.0 program. A statistically significance was defined in our study by P < 0.05 (indicated by ∗ in all figures) or P < 0.01 (indicated by ∗∗ in all figures).
Nucleotide Sequence Accession Numbers
The GenBank accession number for 16S rDNA gene of B. linens BS258 is KU883150. The GenBank accession number for the whole genome of B. linens BS258 is CP014869.1. The GenBank accession numbers for three CAs are A2T55_RS07325, A2T55_RS14335, A2T55_RS17060.
Identification of Calcite-forming Bacteria Isolated from Marine Sediments
To obtain calcite-forming bacteria, ten marine sediment samples were collected from beaches along China Yellow Sea (118°25′E, 35°56′N). After incubated on isolation media, 357 growing bacteria were isolated and further screened by 9219 media to get urease-positive bacteria because calcite-forming is most often correlated with urea hydrolysis (Hammes and Verstraete, 2002). Urea in the media is hydrolyzed by urease produced by bacteria, which, in turn, leads to the increase of pH and color change in media. According to the color change in 9219 urea media, 38 isolates were identified as urease-positive, and strain 258 showed the strongest urease activity. Thereafter, urease-positive bacteria were inoculated in CaCO3-producing media supplemented with 330 mM urea and 150 mM CaCl2, and crystal-precipitating colonies were checked after 5 days cultivation with stereoscopy. Consistent with the urease activity, strain 258 was found to exhibit pronounced calcite biomineralization, which could form 1000s of crystals with different shapes and sizes both on and around the colony (Figures 1A,B and Supplementary Figure S1). In view of the high homology with marine actinobacteria Brevibacterium linens by 16S rDNA phylogenetic analysis (Supplementary Figure S2), this strain was designated as Brevibacterium linens BS258 and was selected for further experiments due to its obvious urease-producing and calcite-precipitating abilities.
FIGURE 1. Microscopic observation of minerals formed by B. linens BS258. Light microscopic observation of minerals formed by B. linens BS258 with inverted microscope (Nikon TS100, Japan) at 12.5 X (A) or 200 X (B). (C) Scanning electron microscopic observation of minerals precipitated by B. linens BS258. (D) Energy dispersive X-ray analysis of the selected area in (C). (E) Transmission electron microscopic observation of minerals precipitated by B. linens BS258. (F) Electron diffraction pattern of the selected area in (E). Typical images (spectra) are shown from many similar examples.
SEM, TEM, EDX and FTIR Analyses of the Calcite Formed by B. linens BS258
To get more insights of the minerals precipitated by B. linens BS258, the precipitates separated from 9219 base liquid media were further chosen for SEM analysis. SEM results revealed that bacteria imbedded in the crystal matrices (Figure 1C and Supplementary Figure S3B), which coincided with the report that the bacteria served as the nucleation sites for the mineralization process (Achal et al., 2009). In addition, imprints of the cells were also observed on the surface of the different minerals (Supplementary Figures S3C–F), which further demonstrated that bacterial cells might provide nucleation sites or structural support for the process of mineral precipitation. The small holes matched the diameters of cells of B. linens BS258 very well, which could also be taken as evidence for the importance of cells in carbonate biomineralization (Supplementary Figures S3A,C,D). Similar results were acquired using TEM on whole single cell and its surroundings. Representative examples of cells in the control (nutrient-rich, without adding CaCl2) and precipitation experiment (nutrient-rich, 30 mM CaCl2 was added) were shown in Figure 1E and Supplementary Figure S4. The cells in the nutrient media without calcium appeared intact without any traces of solid precipitates in the vicinity of the cells (Supplementary Figure S4A). In contrast, significant CaCO3 precipitation occurred in the vicinity of bacterial cells in the media containing calcium (Figure 1E and Supplementary Figures S4B–D).
To ensure that the crystals produced by B. linens BS258 were indeed CaCO3 mineral polymorphs, they were examined using EDX combined with SEM. High amount of calcium in the bacterial sample indicated that calcite was present in the form of CaCO3, also those spectra showed in EDX analysis indicated Ca, C, O peaks associated with calcite crystals (Figure 1D). Moreover, high-resolution images of the crystals from TEM and the associated SAED (Figure 1F) showed that they were well-crystallized calcite single crystals as described previously (Benzerara et al., 2003). Under TEM, the amorphous nature of the mineral precipitates is clearly evidenced (Figure 1E and Supplementary Figure S4). They start to form as thin coatings defining the cell surfaces (Supplementary Figure S4B), and evolve into globular precipitates attached to cell surfaces and within the bacterial organic film between cells (Figure 1E and Supplementary Figure S4C). The globules have a diameter of few 10s to few 100 nm, and show abundant electron-lucent centers resulting from devolatilization caused by electron-beam irradiation (Figure 1F).
Additionally, FTIR was also applied in an attempt to further characterize the crystals formed by B. linens BS258 (Supplementary Figure S5). The FTIR spectra showed broad bands in the high-energy region, which probably belonged to the water-related band (3374 cm-1), indicating the minerals were in a hydrated form. Theoretically, there are four vibrational modes for carbonate ions, ν1 (∼1000 cm-1), ν2 (∼ 873 cm-1), ν3 (1650∼1300 cm-1) and ν4 (∼700 cm-1) (Li et al., 2015). As seen in Supplementary Figure S5, there were also four peaks corresponding to ν1 (1086 cm-1), ν2 (876.5 cm-1), ν3 (1472.6 cm-1) and ν4 (745.5 cm-1) presented in the crystals formed by B. linens BS258, which all indicated the presence of carbonate. Combining the results of SEM, EDX, TEM, SEAD and FTIR, we concluded that the crystals formed by B. linens BS258 were indeed CaCO3.
Calcite Precipitation and Dissolution by B. linens BS258
The concentration of Ca2+ is one of the most important factors determining the formation of calcite (Hammes and Verstraete, 2002). In order to investigate the effects of Ca2+ concentration on the mineralization of B. linens BS258, we checked calcite formation of this strain in the media supplemented with different concentration of Ca2+. As expected, B. linens BS258 could form obvious calcite precipitation on or around the colony but not the area far from the colony in the medium without additional Ca2+ added (Figure 2A). However, this bacterium started to dissolve the formed calcite crystals around the colony when additional 50 mM Ca2+ was added in the medium (Figure 2B). More interestingly, calcite crystals around the colony were completely dissolved when additional 150 mM Ca2+ was added (Figure 2C). It is noteworthy that the size of clear zone gradually became bigger with the concentration of Ca2+ increased up to 500 mM in the medium (Figure 2D). To investigate the details of calcite precipitation and dissolution by B. linens BS258, the plate with 150 mM Ca2+ was checked by light microscopy with high magnification. No any obvious calcite precipitates around the colony were observed (Supplementary Figure S6B) even lots of crystals covered the surface of the colony completely (Supplementary Figure S6A). There was a clear boundary between the clear zone and precipitation layer in the medium, and the density in the clear region margin was even higher than that in the medium (Supplementary Figures S6C,D).
FIGURE 2. Characterization of calcite precipitation of B. linens BS258 in different concentration of Ca2+. Calcite precipitation and dissolution patterns in the CaCO3-producing media supplemented with additional 0 mM (A), 50 mM (B), 150 mM (C) or 500 mM CaCl2 (D) after growth of B. linens BS258 for 5 days at 28°C. Typical images are shown from many similar examples. (E) pH changes of BS258 in the liquid medium amended with additional 0, 50, 150, and 500 mM Ca2+, respectively. (F) Ca2+ removal rate of BS258 in the liquid medium amended with additional 0, 50, 150, and 500 mM Ca2+, respectively.
Changes in pH and Calcium Concentration Over Time during Growth of B. linens BS258
To further understand the changes of pH and Ca2+ concentration in the medium during the process of calcite-producing, B. linens BS258 was incubated in modified LB medium for 5 days, and pH and Ca2+ concentrations were measured at the end of incubation. Correspondingly, B. linens BS258 induced pH increase in media supplemented with different concentrations of Ca2+, which was a key characteristic for CaCO3 bioprecipitation. Meanwhile, ΔpH (pH difference between the final and initial timepoint of incubation) of the medium declined from 2.99 to 1.61 when Ca2+ concentrations increased from 0 to 500 mM (Figure 2E), which indicates that cells produced more H+ with the concentration of Ca2+ increased and is in accordance with the results of calcite dissolution in Figures 2A–D. Notably, as for the media amended with 50 mM and 150 mM Ca2+, more than 99% Ca2+ was removed from the supernatant, and Ca2+ removal rate of the medium with 500 mM Ca2+ was still up to 69.2% (Figure 2F). Therefore, B. linens BS258 will be a good candidate for removal of Ca2+ from industrial wastewater via an effective, eco-friendly and simple method in the future.
Genome Sequencing and Analysis of B. linens BS258
In order to understand the molecular mechanisms of calcite formation by B. linens BS258 and exploit its potentials, whole genome of B. linens BS258 was sequenced. The sequencing results showed that the genome of B. linens BS258 consisted of one circular chromosome of 3,862,244 bp (Table 1). Notably, six genes were found to be involved in urea catabolic process (Supplementary Figure S7), including three urease genes (urease subunit alpha, urease subunit beta, urease subunit gamma) and three urease accessory proteins (UreF, UreH, UreG). To our knowledge, this is the first Brevibacterium sp. with complete genome sequenced, which provides a platform to pursue the molecular mechanisms of CaCO3 biomineralization.
The Involvement of Carbonic Anhydrase in Calcite Dissolution
One of the notable characteristics of calcite precipitation by B. linens BS258 is that it can dissolve formed calcite when Ca2+ concentration changes (Figure 2). Notably, CA derived from Bacillus mucilaginosus was reported to be responsible for the wollastonite dissolution (Xiao et al., 2015). CA is capable of catalyzing the reversible hydration reaction, , of atmospheric and self-generated CO2 (Smith and Ferry, 2000; Xiao et al., 2015). There are three carbonic anhydrases (CA1, CA2, CA3) encoding genes in the genome of B. linens BS258, and the phylogenetic results indicates that CA3 is much closer related to CA2 than that of CA1 (Figure 3A). To investigate whether CA is involved in the calcite dissolution of B. linens BS258, qRT-PCR was performed to check the expression of CA in different conditions where dissolution happened. The results showed that the expressions of three CAs were all dramatically up-regulated when the concentration of calcium increased in the media, which was consistent well with the trend of calcite dissolution (Figure 3B). The consistency between the increase in CA gene expression and calcite dissolution suggests that CAs of B. linens BS258 play a role in the dissolution of calcite at atmospheric CO2 concentrations.
FIGURE 3. Characterization of carbonic anhydrases of B. linens BS258. (A) The consensus phylogenetic tree of three carbonic anhydrases of B. linens BS258 with other related carbonic anhydrases obtained from GenBank (accession numbers are indicated before the species name) constructed by the neighbor-joining method. Numbers above the branches are bootstrap values based on 1000 replicates. (B) The expression change patterns of three carbonic anhydrases of B. linens BS258 with the increase of Ca2+ concentration. All qRT-PCR runs were conducted with three biological and three technical replicates. CA: carbonic anhydrase. P < 0.05 indicated by ∗ and P < 0.01 indicated by ∗∗ versus the control (0 mM CaCl2) in the same group.
Transcriptome Analyses of Calcite Precipitation in Different Concentration of Ca2+
Transcriptome analysis is a powerful tool for revealing the molecular mechanisms underlying the responses of bacteria to environmental stressors. To disclose the overall perspective of biomineralization of B. linens BS258, we investigated the transcriptome of this bacterium under different concentration of Ca2+. The three B. linens BS258 transcriptome examples were Ca0 (without additional Ca2+ amended), Ca50 (with additional 50 mM Ca2+ amended) and Ca150 (with additional 150 mM Ca2+ amended). A total of 1534 unigenes were differentially expressed between any two-way comparison of Ca0, Ca50 and Ca150 (| fold changes| > 2, corrected P < 0.005) (Figure 4A).
FIGURE 4. Transcriptome profiles of B. linens BS258 challenged in different concentration of Ca2+. (A) Differentially expressed genes (DEGs) that are unique or shared among three groups of Ca0 (without additional Ca2+ amended), Ca50 (with additional 50 mM Ca2+ amended) and Ca150 (with 150 additional Ca2+ amended). Numbers of each circle show the number of Ca2+ stress responsive DEGs that are uniquely (inside of non-overlapping part) or commonly (inside of overlapping part) regulated. (B) Heat map shows the Ca2+ stress responsive relatively up-regulated DEGs related to quorum sensing among Ca0, Ca50, and Ca150 enriched by KEGG analysis. (C) Heat map shows the Ca2+ stress responsive relatively up-regulated DEGs related to ribosome biosynthesis among Ca0, Ca50, and Ca150 enriched by KEGG analysis. (D) Heat map shows the Ca2+ stress responsive relatively up-regulated DEGs related to quorum sensing among Ca0, Ca50, and Ca150 enriched by KEGG analysis. All the Genbank accession numbers are provided together with the corresponding heat map. The heat map is made by the software of Heml 18.104.22.168.
KEGG enrichment analysis showed that QS, ribosome biosynthesis, and ABC transporter were enriched among the two-way comparison of Ca0, Ca50, and Ca150 (Supplementary Tables S2A–C). Under the term of “Quorum Sensing,” expression of 19 DEGs was enhanced continuously by the increase of calcium concentration (Figure 4B). Notably, most up-regulated genes related to QS were found to encode peptide/nickel ABC transporters. One of the DEGs in QS was also controlled by “Fatty acid metabolism,” namely long-chain fatty acid-CoA ligase (A2T55_RS10245). Ribosome is an intracellular organelle that serves as the site of biological protein synthesis. There were 39 DEGs involved in ribosome biosynthesis that showed continuously up-regulated trend with the increase of calcium concentrations (Figure 4C). As the calcium concentration increased, the expression of 28 DEGs enriched in “ABC transporter” was progressively up-regulated (Figure 4D). The substrates transported by the ABC transporters of the 28 up-regulated DEGs were sugar, amino acid, spermidine/putrescine, metal, iron and molybdate (Figure 5). The DEG of ABC transporter with the highest fold change was spermidine/putrescine ABC transporter permease (A2T55_RS03660). To disclose the relationship between the QS-related transporters and other 28 DEGs enriched ABC transporters, we analyzed the phylogenetic relationship of all transporters mentioned above. Interestingly, the up-regulated QS-related transporters annotated as ATP binding protein and permease were clustered together (Figure 5). On the other hand, 28 DEGs enriched ABC transporters were clustered into four subgroups, ATP binding protein, two permeases and substrate binding protein (Figure 5), which comprise the full functional counterparts of ABC transporter (Hiron et al., 2010).
FIGURE 5. Phylogenetic analysis of up-regulated DEGs linked to QS-related ABC transporters and other regular ABC transporters among Ca0, Ca50, and Ca150 enriched by KEGG analysis in B. linens BS258. The consensus phylogenetic tree is constructed by the neighbor-joining method. Numbers above the branches are bootstrap values based on 1000 replicates. The star (∗) means the gene is out of the exact scope of ATP binding, permease and substrate binding counterparts of ABC transporter. The arrow indicates the position of QS-related peptide/nickel ABC transporter substrate-binding protein.
Verification of RNA-seq Data
Reliability of the RNA-seq data was verified by qRT-PCR analysis. Totally, five genes were selected for validation. Consistently, similar trends were observed between qRT-PCR and RNA-seq results, which supported the validity of the RNA-seq data (Supplementary Figure S8).
Heavy Metal Sequestration Assays of B. linens BS258
Microbial carbonate precipitation is known as an efficient process for the remediation of heavy metals from contaminated environments (Li et al., 2015). Next, we sought to check the capability of B. linens BS258 to remove heavy metals. Within the genome of B. linens BS258, there are both cadmium and lead resistance gene clusters. Therefore, we checked the sequestration abilities of Cd2+ and Pb2+ by B. linens BS258. The results showed that B. linens BS258 only removed very few heavy metals if Ca2+ and urea were absent in the media. However, heavy metal removal abilities increased dramatically if 30 mM Ca2+ and 80 mM urea were present in the media (Figure 6). Taken together, the above results indicated that B. linens BS258 can only remove these heavy metals by co-precipitation with Ca2+ by B. linens BS258. We proposed that the possible molecular structures of these minerals could be mixtures of CaCO3, XCO3 and (CanX1-n)CO3 (X is Cd or Pb) as described previously (Li et al., 2015). On the other hand, after incubation supplemented with different heavy metals, all media became alkaline (∼pH 9.2), which is similar to that of medium with only Ca2+. These results showed that supplement of heavy metal ions didn’t have obvious impact on the Ca2+ precipitation, while heavy metal ions could be effectively precipitated.
FIGURE 6. Heavy metal ions removal by B. linens BS258. (A) Cd2+ removal test by B. linens BS258 in different concentration of Cd2+ with or without amended with 30 mM Ca2+ and 80 mM urea in the medium. (B) Pb2+ removal test by B. linens BS258 in different concentration of Pb2+ with or without amended with 30 mM Ca2+ and 80 mM urea in the medium.
In this study, the marine derived actinobacterium B. linens BS258 was demonstrated to effectively precipitate CaCO3 minerals in the plates and liquid media combining the results of SEM, EDX, TEM, SAED and FTIR (Figure 1 and Supplementary Figures S1, S3–S5). It is known that MICP is a rather straightforward chemical process governed by four key factors: (i) Ca2+ concentration, (ii) the concentration of dissolved inorganic carbon, (iii) the pH and (iv) the availability of nucleation sites (Hammes and Verstraete, 2002). Correspondingly, in our study, the pH of the liquid cultures (Figure 2E), amount of carbonate ion could be increased significantly through the ureolytic pathway of B. linens BS258 (Figure 7 and Supplementary Figure S7 and reaction 1). Thus, in the environment with free Ca2+, CaCO3 could be effectively formed by using cells as nucleation sites (Figure 1C).
FIGURE 7. The proposed model of calcite precipitation by B. linens BS258 in the ocean ecosystem. Firstly, B. linens BS258 produces QS signals at a basal rate as the cells grow and divide until a threshold concentration of QS signals (like nickel concentration) has accumulated, and then the bacteria start to form calcite precipitation with the synergic action of urease, carbonic anhydrase and other necessary factors in the ocean surface. Secondly, the formed minerals precipitate to the sea floor and form the sediments. Lastly, the sediments are dissolved in the harsh conditions like high concentration heavy metal ions and the free cells go back to the ocean surface for the next round calcification. Therefore, the concentration of heavy metals in the environment is a key factor determining the equilibration between calcite precipitation and dissolution. During the course of calcite precipitation and dissolution, many energy metabolisms close related to tricarboxylic acid cycle, such as sugar biosynthesis and glycolysis, amino acid biosynthesis and degradation, fatty acid metabolism, are up-regulated. Correspondingly, the activities of many transporters of metal ions, amino acid and sugar are also positive regulated. More details are described in the Section “Discussion.” TCA: tricarboxylic acid cycle. ABC transporter: ATP-binding cassette transporters.
One of the important reasons for microbes to form minerals is to enhance their competitiveness in a microbial community via the detoxification of poisonous inorganic species (Ehrlich, 1996). Generally, bacteria keep tight control of their intracellular Ca2+ concentration at 100–300 nM because of the key role of calcium in the regulation of many essential processes and its potential detriment to cell structures (Dominguez, 2004). In this study, for B. linens BS258, in the medium amended with 150 mM CaCl2, the extracellular Ca2+ concentration is about 105 times higher than that of intracellular, thus, passive transport usually accounts for calcium influx as previous reported (Norris et al., 1996; Hammes and Verstraete, 2002). To survive in such harsh conditions, the cells are required to actively export of intracellular excess calcium to reduce intracellular Ca2+ and simultaneously compensate the proton loss, which could be performed via the ATP dependent calcium pumps (Figure 7). With the extrusion of Ca2+ and uptake of protons in the same region, a localized increase in pH would happen, together with the alkalization caused by ureolytic pathway would form an ideal localized precipitation micro-environment (Figure 7).
Another important characteristic of B. linens BS258 toward calcite precipitation is that the formed calcite around the colony could be dissolved to some extent when increasing the concentration of Ca2+ in the media (Figures 2 and Supplementary Figure S6). For B. linens BS258, the expression levels of CAs were dramatically up-regulated when higher calcium concentration existed in the medium (Figure 3). CA is able to catalyze 1 molar CO2 to 1 molar and H+, which suggested that the production of H+ would be up-regulated while calcium concentration increased. Once the concentration of H+ is high enough to dissolve the formed calcite precipitation around the colony, the calcite dissolution happens (Figure 2).
It is reported that microbes are able to dissolve a mineral by using it as a source of energy to enhance competitiveness in a microbial community (Ehrlich, 1996). CA can participate in the formation of malonyl-CoA with bicarbonate and acetyl-CoA as the substrate (Costa et al., 2007), which means CA is an important regulator of fatty acid metabolism. Additionally, the dissolution of calcite generated a certain amount of by CO2 hydration, which is a key substrate for many essential biological pathways including gluconeogenesis, lipogenesis, ureagenesis and synthesis of several amino acids (Figure 7) (Xiao et al., 2015). For B. linens BS258, the expressions of three CAs were all markedly up regulated with the increase of Ca2+ concentration (Figure 3B), which is in accordance with the suggestion that an increase in CA gene expression level is advantageous to the microbe’s survival chances in the harsh environments (Xiao et al., 2015).
To better understand the calcification of B. linens BS258, the whole genome of this bacterium was sequenced and the transcriptome profiling of B. linens BS258 challenged by different condition of Ca2+ was further performed. Surprisingly, KEGG enrichment analysis showed that only three pathways related to QS, ribosome biosynthesis, and ABC transporter were enriched (Figure 4). With the increase of Ca2+ concentration, the expression of ribosome biosynthesis and ABC transporter in charge of sugar and amino acid transport was significantly up-regulated (Figures 4C,D), which guarantees B. linens BS258 to reach higher microbial density. Consistently, the expression of genes related to QS was also markedly up-regulated (Figure 4B). It is reported that QS enhancement of stress tolerance, therefore having a broad impact in bacterial ecology (Garcia-Contreras et al., 2015). Our transcriptome data indicate that B. linens BS258 changes regulation of its biochemical pathways in a manner that is adaptive for a cooperative lifestyle in the presence of calcium to obtain phytoplankton-derived nutrients and higher microbial density.
It is noteworthy that most up-regulated genes related to QS were all peptide/nickel ABC transporters counterparts (Supplementary Table S2B). Interestingly, in bacteria, nickel uptake is mediated by peptide/nickel ABC transporters (Hiron et al., 2010). As soluble Ni2+ in natural environments usually exists in trace amounts, high-affinity uptake of nickel is needed to assure intracellular metalloenzyme activities (Hiron et al., 2010). It is well known that urease is a multimeric nickel-dependent enzyme mainly involved in the neutralization of acidic environments and influence the MICP to a great extent (Achal et al., 2009), and the activity of urease was determined by nickel transporters (Hiron et al., 2010). Hiron et al. (2010) found that a nickel ABC-transporter of Staphylococcus aureus was involved in urinary tract infection by forming infection stones through ureolytic MICP. Therefore, it is understandable that the expressions of QS-related peptide/nickel ABC transporters in B. linens BS258 were all evidently up-regulated with the increase of Ca2+ concentration. More functions of QS and its related peptide/nickel ABC transporters need to investigate further and eventually disclose the full perspective of biomineralization as illustrated in Figure 7.
The ecological relevance of MICP has been exploited in biotechnology with applications in many fields including bioremediation of heavy metal contamination (Li et al., 2015). In our experiments, B. linens BS258, incubated in media supplemented with urea and different metal ions showed very good metal removal abilities through the co-precipitation with calcium (Figure 6), which indicates that this bacterium is powerful and promising in the application of heavy metal pollution management.
YZ, NM, and CS conceived and designed the experiments. YZ finished the calcite production-bacteria screening and genome sequencing. NM performed the transcriptome analyses and heavy metal removal experiments. WJ helped to do the FTIR analysis. YZ, NM, and CS analyzed the data. YZ, NM, and CS prepared the figures and wrote the paper. SW helped to write the manuscript. All authors reviewed the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by Young TaiShan Scholars Foundation of Shandong Province, Natural Science Outstanding Youth Foundation of Shandong Province (JQ201607), AoShan Talents Program supported by Qingdao National Laboratory for Marine Science and Technology (No.2015ASTP), “100-Talent Project” of Chinese Academy of Sciences for Chaomin Sun.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2017.00602/full#supplementary-material
Achal, V., Mukherjee, A., Basu, P. C., and Reddy, M. S. (2009). Strain improvement of Sporosarcina pasteurii for enhanced urease and calcite production. J. Ind. Microbiol. Biotechnol. 36, 981–988. doi: 10.1007/s10295-009-0578-z
Barabesi, C., Galizzi, A., Mastromei, G., Rossi, M., Tamburini, E., and Perito, B. (2007). Bacillus subtilis gene cluster involved in calcium carbonate biomineralization. J. Bacteriol. 189, 228–235. doi: 10.1128/Jb.01450-06
Benzerara, K., Menguy, N., Guyot, F., Dominici, C., and Gillet, P. (2003). Nanobacteria-like calcite single crystals at the surface of the Tataouine meteorite. Proc. Natl. Acad. Sci. U.S.A. 100, 7438–7442. doi: 10.1073/pnas.0832464100
Benzerara, K., Skouri-Panet, F., Li, J. H., Ferard, C., Gugger, M., Laurent, T., et al. (2014). Intracellular Ca-carbonate biomineralization is widespread in cyanobacteria. Proc. Natl. Acad. Sci. U.S.A. 111, 10933–10938. doi: 10.1073/pnas.1403510111
Berelson, W. M., Balch, W. M., Najjar, R., Feely, R. A., Sabine, C., and Lee, K. (2007). Relating estimates of CaCO3 production, export, and dissolution in the water column to measurements of CaCO3 rain into sediment traps and dissolution on the sea floor: a revised global carbonate budget. Global Biogeochem. Cycles 21, 1–15. doi: 10.1029/2006gb002803
Braissant, O., Decho, A. W., Dupraz, C., Glunk, C., Przekop, K. M., and Visscher, P. T. (2007). Exopolymeric substances of sulfate-reducing bacteria: interactions with calcium at alkaline pH and implication for formation of carbonate minerals. Geobiology 5, 401–411. doi: 10.1111/j.1472-4669.2007.00117.x
Bundeleva, I. A., Shirokova, L. S., Benezeth, P., Pokrovsky, O. S., Kompantseva, E. I., and Balor, S. (2012). Calcium carbonate precipitation by anoxygenic phototrophic bacteria. Chem. Geol. 291, 116–131. doi: 10.1016/j.chemgeo.2011.10.003
Bundeleva, I. A., Shirokova, L. S., Pokrovsky, O. S., Benezeth, P., Menez, B., Gerard, E., et al. (2014). Experimental modeling of calcium carbonate precipitation by cyanobacterium Gloeocapsa sp. Chem. Geol. 374, 44–60. doi: 10.1016/j.chemgeo.2014.03.007
Costa, M., Borges, C. L., Bailao, A. M., Meirelles, G. V., Mendonca, Y. A., Dantas, S. F. I. M., et al. (2007). Transcriptome profiling of Paracoccidioides brasiliensis yeast-phase cells recovered from infected mice brings new insights into fungal response upon host interaction. Microbiology 153, 4194–4207. doi: 10.1099/mic.0.2007/009332-0
Garcia-Contreras, R., Nunez-Lopez, L., Jasso-Chavez, R., Kwan, B. W., Belmont, J. A., Rangel-Vega, A., et al. (2015). Quorum sensing enhancement of the stress response promotes resistance to quorum quenching and prevents social cheating. ISME J. 9, 115–125. doi: 10.1038/ismej.2014.98
Han, Z. Z., Yan, H. X., Zhao, H., Zhou, S. X., Han, M., Meng, X. Q., et al. (2014). Bio-precipitation of calcite with preferential orientation induced by Synechocystis sp. PCC6803. Geomicrobiol. J. 31, 884–899. doi: 10.1080/01490451.2014.907379
Hiron, A., Posteraro, B., Carriere, M., Remy, L., Delporte, C., La Sorda, M., et al. (2010). A nickel ABC-transporter of Staphylococcus aureus is involved in urinary tract infection. Mol. Microbiol. 78, 788–788. doi: 10.1111/j.1365-2958.2010.07416.x
Kang, C. H., Choi, J. H., Noh, J., Kwak, D. Y., Han, S. H., and So, J. S. (2014). Microbially induced calcite precipitation-based sequestration of strontium by Sporosarcina pasteurii WJ-2. Appl. Biochem. Biotechnol. 174, 2482–2491. doi: 10.1007/s12010-014-1196-4
Li, C. S., Xu, Y., Jiang, W., Dong, X. Y., Wang, D. F., and Liu, B. J. (2013). Effect of NaCl on the heavy metal tolerance and bioaccumulation of Zygosaccharomyces rouxii and Saccharomyces cerevisiae. Bioresour. Technol. 143, 46–52. doi: 10.1016/j.biortech.2013.05.114
Marvasi, M., Visscher, P. T., Perito, B., Mastromei, G., and Casillas-Martinez, L. (2010). Physiological requirements for carbonate precipitation during biofilm development of Bacillus subtilis etfA mutant. FEMS Microbiol. Ecol. 71, 341–350. doi: 10.1111/j.1574-6941.2009.00805.x
Mitchell, A. C., Dideriksen, K., Spangler, L. H., Cunningham, A. B., and Gerlach, R. (2010). Microbially enhanced carbon capture and storage by mineral-trapping and solubility-trapping. Environ. Sci. Technol. 44, 5270–5276. doi: 10.1021/es903270w
Obst, M., Dynes, J. J., Lawrence, J. R., Swerhone, G. D. W., Benzerara, K., Karunakaran, C., et al. (2009). Precipitation of amorphous CaCO3 (aragonite-like) by cyanobacteria: a STXM study of the influence of EPS on the nucleation process. Geochim. Cosmochim. Acta 73, 4180–4198. doi: 10.1016/j.gca.2009.04.013
Qin, J. Y., Wang, X. W., Wang, L. D., Zhu, B. B., Zhang, X. H., Yao, Q. S., et al. (2015). Comparative transcriptome analysis reveals different molecular mechanisms of Bacillus coagulans 2-6 response to sodium lactate and calcium lactate during lactic acid production. PLoS ONE 10:e0124316. doi: 10.1371/journal.pone.0124316
Ries, J. B., Anderson, M. A., and Hill, R. T. (2008). Seawater Mg/Ca controls polymorph mineralogy of microbial CaCO3: a potential proxy for calcite-aragonite seas in Precambrian time. Geobiology 6, 106–119. doi: 10.1111/j.1472-4669.2007.00134.x
Rodriguez-Navarro, C., Jroundi, F., Schiro, M., Ruiz-Agudo, E., and Gonzalez-Munoz, M. T. (2012). Influence of substrate mineralogy on bacterial mineralization of calcium carbonate: implications for stone conservation. Appl. Environ. Microbiol. 78, 4017–4029. doi: 10.1128/Aem.07044-11
Schut, F., Devries, E. J., Gottschal, J. C., Robertson, B. R., Harder, W., Prins, R. A., et al. (1993). Isolation of typical marine-bacteria by dilution culture - growth, maintenance, and characteristics of isolates under laboratory conditions. Appl. Environ. Microbiol. 59, 2150–2160.
Wu, S. M., Liu, G., Jin, W. H., Xiu, P. Y., and Sun, C. M. (2016). Antibiofilm and anti-Infection of a marine bacterial exopolysaccharide against Pseudomonas aeruginosa. Front. Microbiol. 7:102. doi: 10.3389/fmicb.2016.00102
Xiao, L. L., Lian, B., Hao, J. C., Liu, C. Q., and Wang, S. J. (2015). Effect of carbonic anhydrase on silicate weathering and carbonate formation at present day CO2 concentrations compared to primordial values. Sci. Rep. 5:7733. doi: 10.1038/srep07733
Keywords: genomic, transcriptomic, marine, actinobacterium, biomineralization, heavy metal, bioremediation
Citation: Zhu Y, Ma N, Jin W, Wu S and Sun C (2017) Genomic and Transcriptomic Insights into Calcium Carbonate Biomineralization by Marine Actinobacterium Brevibacterium linens BS258. Front. Microbiol. 8:602. doi: 10.3389/fmicb.2017.00602
Received: 20 February 2017; Accepted: 23 March 2017;
Published: 06 April 2017.
Edited by:Qiang Wang, Institute of Hydrobiology (CAS), China
Reviewed by:Xiaoke Hu, Yantai Institute of Coastal Zone Research (CAS), China
Bruno A. Martinez, Lawrence Berkeley National Laboratory, USA
Copyright © 2017 Zhu, Ma, Jin, Wu and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chaomin Sun, email@example.com
†These authors have contributed equally to this work.