ORIGINAL RESEARCH article
In Vitro and in Silico Evidence of Phosphatase Diversity in the Biomineralizing Bacterium Ramlibacter tataouinensis
- 1Centre National de la Recherche Scientifique, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Universités, UMR 7590, Muséum National d'Histoire Naturelle, Université Pierre et Marie Curie, IRD 206, Paris, France
- 2Department of Geological Sciences, University of Colorado, Boulder, CO, United States
- 3Laboratoire d'Écologie Microbienne de la Rhizosphère et Environnements Extrêmes, UMR 7265, Aix Marseille Univ, Centre National de la Recherche Scientifique, Commissariat à l'Énergie Atomique et aux Énergies Alternatives, Saint-Paul-lez-Durance, France
Microbial phosphatase activity can trigger the precipitation of metal-phosphate minerals, a process called phosphatogenesis with global geochemical and environmental implications. An increasing diversity of phosphatases expressed by diverse microorganisms has been evidenced in various environments. However, it is challenging to link the functional properties of genomic repertoires of phosphatases with the phosphatogenesis capabilities of microorganisms. Here, we studied the betaproteobacterium Ramlibacter tataouinensis (Rta), known to biomineralize Ca-phosphates in the environment and the laboratory. We investigated the functional repertoire of this biomineralization process at the cell, genome and molecular level. Based on a mineralization assay, Rta is shown to hydrolyse the phosphoester bonds of a wide range of organic P molecules. Accordingly, its genome has an unusually high diversity of phosphatases: five genes belonging to two non-homologous families, phoD and phoX, were detected. These genes showed diverse predicted cis-regulatory elements. Moreover, they encoded proteins with diverse structural properties according to molecular models. Heterologously expressed PhoD and PhoX in Escherichia coli had different profiles of substrate hydrolysis. As evidenced for Rta cells, recombinant E. coli cells induced the precipitation of Ca-phosphate mineral phases, identified as poorly crystalline hydroxyapatite. The phosphatase genomic repertoire of Rta (containing phosphatases of both the PhoD and PhoX families) was previously evidenced as prevalent in marine oligotrophic environments. Interestingly, the Tataouine sand from which Rta was isolated showed similar P-depleted, but Ca-rich conditions. Overall, the diversity of phosphatases in Rta allows the hydrolysis of a broad range of organic P substrates and therefore the release of orthophosphates (inorganic phosphate) under diverse trophic conditions. Since the release of orthophosphates is key to the achievement of high saturation levels with respect to hydroxyapatite and the induction of phosphatogenesis, Rta appears as a particularly efficient driver of this process as shown experimentally.
Phosphorus is essential to life but a limiting nutrient in many ecosystems (Paytan and McLaughlin, 2007; Pasek, 2008). In response, cellular adaptations have been selected in microorganisms to catalyze an efficient uptake and recycling of P (Sowell et al., 2008; Martiny et al., 2009; Coleman and Chisholm, 2010; Sebastián et al., 2012; Martin et al., 2014). Under certain conditions, the microbial release of phosphorus generates locally high P concentrations, sometimes triggering the precipitation of metal-phosphate phases, a process called phosphatogenesis (Föllmi, 1996; Arning et al., 2009; Goldhammer et al., 2010; Cosmidis et al., 2014). The implication of bacteria in the biomineralization of metal-phosphates, including calcium-, iron- or uranium-phosphates has been suggested based on (1) the frequent occurrence of bacteria fossil in phosphate-rich sedimentary deposits (Zanin and Zamirailova, 2011; Bailey et al., 2013; Cosmidis et al., 2013a,b) and (2) laboratory biomineralization experiments (Powers et al., 2002; Templeton et al., 2003; Beazley et al., 2009; Miot et al., 2009; Rivadeneyra et al., 2010; Yung and Jiao, 2014; Cosmidis et al., 2015). Deciphering the microbial processes involved will help designing better models for the formation of sedimentary phosphate minerals, a major sink in the global cycle of P (Föllmi, 1996), as well as developing new bioremediation strategies to sequester metal pollutants (e.g., Martinez et al., 2007; Handley-Sidhu et al., 2010; Mondani et al., 2011; Liang et al., 2016).
Several studies suggest that microbial phosphatases are particularly important in phosphatogenesis (Hirschler et al., 1990; Blake et al., 1998; Macaskie et al., 2000; Beazley et al., 2007; Nilgiriwala et al., 2008; Shelobolina et al., 2009). Phosphatases catalyse the hydrolysis of phosphoesters, an important reservoir of P in many environments (Clark et al., 1998; Young and Ingall, 2010). This hydrolysis releases orthophosphates, favoring mineral precipitation in the presence of cations such as Ca2+, Fe2+, Fe3+, or U6+ (Powers et al., 2002; Cosmidis et al., 2015). Several types of phosphatases are known, including the alkaline phosphatase superfamily (Coleman, 1992) composed of at least 3 non-homologous families named PhoA, PhoD, and PhoX (Kim and Wyckoff, 1991; Rodriguez et al., 2014; Yong et al., 2014). Interestingly, this superfamily exhibits a broad diversity in (1) substrate specificity profile: for example, some proteins hydrolyze phosphate monoesters such as glycerophosphate or phosphoproteins, while other are more efficient at hydrolyzing phosphate diesters such as nucleic acids; (2) cofactors such as Zn2+, Ca2+, Mg2+, Fe3+; (3) cellular location; and (4) regulation pathways (Eder et al., 1996; Wojciechowski et al., 2002; Zalatan et al., 2006; Zaheer et al., 2009). Finally, the environmental distribution of these families differs, suggesting an adaptive strategy of microbial communities in response to multiple environmental drivers (Ragot et al., 2017), e.g., the relative availability of inorganic phosphorus and metal co-factors (Kathuria and Martiny, 2011; Moore, 2014). For example, PhoX and PhoD families are more abundant than PhoA under marine Zn- and P-depleted conditions (Luo et al., 2009; Martiny et al., 2009; Sebastian and Ammerman, 2009; Temperton et al., 2011). Yet, we have a poor understanding on how this diversity (at cell or community level) controls the capability of ecosystems to induce Ca-phosphate biomineralization. It can be hypothesized that in a given environment where a certain type of organic P substrate prevails and/or some metal co-factors are limiting, some phosphatases might be more efficient than others at releasing orthophosphates and therefore inducing phosphatogenesis. With the possibility to access the whole genomes of numerous microorganisms, an outstanding question to address is whether a microbe's phosphatase gene content can help predict its mineralization capabilities under particular conditions.
In this context, we studied Ramlibacter tataouinensis strain TTB310 (Rta), a betaproteobacterium isolated from a semi-arid region in Tunisia, where abundant calcification occurs (Gillet et al., 2000; Benzerara et al., 2003; Heulin et al., 2003). Rta is able to precipitate calcium phosphates in vitro (Benzerara et al., 2004). Although the involvement of a phosphatase was speculated, the molecular mechanisms inducing biomineralization by Rta remain unknown. The genome of the strain has been sequenced and annotated (De Luca et al., 2011), offering the possibility to explore the functional repertoire of Rta phosphatases. Here, we combined bioinformatics, molecular biology, biochemistry and mineralogy to characterize Ramlibacter tataouinensis phosphatase properties in relation with its capability to induce phosphatogenesis.
Materials and Methods
All reactants, including p-nitrophenylphosphate (pNPP), bis p-nitrophenylphosphate (bis-pNPP), Thymidine 5′-monophosphate p-nitrophenylester sodium salt (TpNPP), calcium glycerophosphate (CaGP) and sodium glycerophosphate (NaGP) were purchased from Sigma-Aldrich. pNPP and glycerophosphate (GP) are phosphomonoesters while bis-pNPP and TpNPP are phosphodiesters.
Ramlibacter tataouinensis strain TTB310 (Rta) was described by Heulin et al. (2003) and its genome analyzed by De Luca et al. (2011). Here, Rta was cultivated in fivefold diluted LB medium (LB) at 30°C with orbital shaking (100 rpm) in the dark. The stationary phase (cell density between 107 and 108 cells/ml) was achieved after 4 weeks.
Gene Cloning and Recombinant Expression of Phosphatase Genes in E. coli
Phosphatase genes predicted from Rta genome (see below for details on the functional annotation procedure) were recombinantly expressed in E. coli to characterize some properties of these enzymes. Sequences were optimized for heterologous expression by using codons more commonly found in E. coli (Cosmidis et al., 2015). Synthetic genes were synthesized in one step by PCR from long synthetic oligonucleotides (Gencust service) and inserted in a pET 22b vector (Novagen) at MscI/Bamh1 sites. E. coli strain BL21 5(DE3) (Agilent) was transformed by pET vectors with synthetic phosphatase genes. One hundred milliliters of E. coli culture were grown in LB medium with 50 μg/ml ampicillin at 37°C, 180 rpm orbital agitation. Expression of phosphatase genes was induced at 30°C by adding 0.5 mM of isopropyl αD-thiogalactoside (IPTG) when cultures reached an OD of 0.8 at 600 nm. After 4 h of expression, E. coli cells were centrifuged at 4,000 g and cellular pellets were stored at −20°C.
Calcification assays were conducted to assess the capability of Rta cells and extracellular extracts to induce the precipitation of Ca-phosphate minerals by enzymatic hydrolysis of a phosphomonoester (glycerophosphate). For that purpose, 30 ml of a 2-week old culture were harvested by centrifugation for 10 min at 4,000 g. The supernatant (extracellular extract) was concentrated 200 times on a 5,000 daltons polyethersulfone membrane inserted in a Millipore Amicon stirred cell and tested for calcification after a hundred-fold dilution. Cell pellets were washed with a 20 mM HEPES buffer at pH 7.5 and concentrated 6-fold. Cells or extracellular extracts were added to a solution composed of 10 mM calcium glycerophosphate and 20 mM HEPES at pH 7.5 (CaGP). No other source of phosphorus (either organic or inorganic) was added. Control assays, without Ca, used 10 mM sodium glycerophosphate and 20 mM HEPES pH 7.5 (NaGP). In order to test the impact of Ca on phosphatase activity with no interference from Ca-phosphate precipitation, additional assays in a medium composed of 10 mM of sodium glycerophosphate as the sole source of phosphorus, 0.75 mM of CaCl2 and 20 mM of HEPES at pH 7.5 (NaGP+Ca), were performed. The assays, in either CaGP, NaGP, or NaGP+Ca media, were run for 35 days at 30°C in the dark with shaking.
To test calcification by transformed E. coli strains, cells were grown overnight, and harvested by centrifugation for 10 min at 4,000 g. Cell pellets were washed with 20 mM of HEPES buffer at pH 7.5 and resuspended in assay solutions at a concentration of 108 cells/ml. Assays were run for 7 days at 37°C with shaking.
Concentrations of dissolved Ca and inorganic phosphate (orthophosphate) were measured at different times in the assay solutions. For this purpose, half a milliliter of the solutions was sampled at different time intervals. Samples were centrifuged at 6,000 g for 10 min. Supernatants were filtered at 0.2 μm. The obtained dissolved fraction (solute and particles smaller than 0.2 μm) was acidified with 1% HNO3 and stored at 4°C. Dissolved calcium concentrations were determined using the cresolphthalein complexone spectrophotometry procedure (Moorehead and Biggs, 1974). Dissolved orthophosphate concentrations were determined using the ascorbic acid spectrophotometry procedure (Chen et al., 1956).
Phosphatase Activity Assays
Phosphatase activity was measured using a conventional spectrophotometric assay based on the ability of phosphatases to hydrolyse a derivative of pNPP (p-nitrophenylphosphate) to p-nitrophenol, a chromogenic product with a maximal absorbance at 405 nm (Engvall, 1980). Conversion between absorption at 405 nm and the concentration of released p-nitrophenol (hence orthophosphate) was achieved using an extinction coefficient of 18,000 M−1 cm−1 at 405 nm (Zhang and VanEtten, 1991). Three different pNPP derivatives were used as substrates at a concentration of 600 μM: p-nitrophenylphosphate monoester (pNPP), bis p-nitrophenylphosphate diester (bis-pNPP) and Thymidine 5′-monophosphate p-nitrophenylester (TpNPP) as a deoxynucleotide analog.
Assays were performed on total soluble proteins extracted from the bacterial cells in order to ensure an optimal access of the enzymes to the substrates. Total soluble proteins were extracted from stationary cultures of Rta (~2 × 109 cells) and E. coli (~4 × 109 cells), using the “Bug buster protein extraction” reagent (Novagen), following the protocol recommended by the supplier. This reagent allowed gentle disruption of cell walls, resulting in the release of soluble proteins (including periplasmic proteins). The lysate was cleared by centrifugation at 13,000 g to get the soluble fraction. The phosphatase activities of the total soluble protein of Rta and E. coli cultures were measured in a buffer containing 50 mM Tris at pH 8, 5 mM of CaCl2, 5 mM of MgCl2, and 0.6 mM of substrate (pNPP, bis-pNPP, or TpNPP).
The following controls were used. (i) We checked that no spontaneous hydrolysis of the 3 pNPP derivatives occurred in the absence of phosphatases under the same chemical conditions as those used for phosphatase activity assays, i.e., after addition of the derivatives to the “Bug Buster protein extraction” reagent in 50 mM Tris, with 5 mM of CaCl2, 5 mM of MgCl2 and at pH 8.5. (ii) An E. coli strain BL21 5(DE3) transformed by a pET vector containing a gene not encoding a phosphatase but a small heat shock protein (Uniprot accession Q38806) was incubated in LB with ampicillin. IPTG was added at 30°C for 4 h when the culture reached an OD of 0.8 at 600 nm. The lysate had no significant phosphatase activity (data not shown). (iii) Last, we used an E. coli strain BL21 5(DE3) transformed by a pET vector containing the E. coli phoA gene as a positive control. PhoA is a well-known monoester phosphatase (Cosmidis et al., 2015). This strain was processed similarly to the other transformants. We tested the phosphatase activity of the culture lysate after IPTG induction, toward the 3 pNPP derivatives.
Protein extracts were loaded on 10% polyacrylamide gels for a denaturing SDS PAGE electrophoresis. Gels and samples were prepared as described by Laemmli (1970). Molecular weight standards and gel staining solution (EZ blue) were purchased from Biorad and Sigma-Aldrich, respectively.
Calculation of Supersaturation Indices with Respect to Hydroxyapatite
The saturation index is defined as the decadic logarithm of the ratio of the ion activity product (IAP) over the solubility constant (Ks): SI = log (IAP/Ks). Dissolved calcium and orthophosphate concentrations were used to calculate the saturation index of the solutions with respect to hydroxyapatite, assuming a pH of 7.5, using the CHESS code (van der Lee, 1998) and a solubility constant for hydroxyapatite KHA as log(KHA) = −57.74 (Fujita et al., 2007). Cosmidis et al. (2015) measured that pH remained constant at 7.5 in a calcification assay similar as the one used here. Saturation indices were not calculated when orthophosphate or calcium concentrations were below the detection limits.
Fourier-Transform Infrared Spectroscopy
Precipitates formed by Rta cells or recombinant E. coli cells expressing either E. coli PhoA or Rta PhoX1 enzymes in CaGP medium were harvested by centrifugation after 35 days (for Rta) or 7 days of incubation (for recombinant E. coli). Samples were washed in a 20 mM HEPES buffer at pH 7.5 first, then in MilliQ water before drying for 48 h at ~45°C. About 3 mg of sample were gently ground in an agate mortar and mixed with ~300 mg of dried potassium bromide. Samples were pelleted under 9 tons of pressure for ~1 min and dried at ~90°C. After a second pressing, Fourier-Transform InfraRed (FT-IR) spectra were recorded between 400 and 4,000 cm−1 with a resolution of 1 cm−1 using a Nicolet 7600 FT-IR spectrometer. One hundred spectra were averaged for each sample. The baseline subtraction was done with the program Omnic 7.3. FT-IR spectra of E. coli cells grown in LB medium, and a reference hydroxyapatite (Sigma-Aldrich) were measured for comparison. All spectra were normalized to the maximum peak intensity.
Scanning Electron Microscopy
Scanning electron microscopy (SEM) analyses were performed using a Zeiss Ultra 55 field emission gun SEM. Images were acquired using an Everhart-Thornley detector at an accelerating voltage of 10 kV and a working distance of ~7 mm or an in column detector (InLens) at 2 kV and a working distance of 1.5 mm. Samples were prepared following the same protocol as Cosmidis et al. (2015). Briefly, cell suspensions were washed and filtered through a 0.2 μm polycarbonate filter (GTTP, Millipore). Filters were then air-dried and coated with carbon, prior to SEM analysis.
Functional Annotation of Phosphatase Genes
Eight families of non-specific phosphatases were previously described in the literature, with an optimal activity in an alkaline (PhoA, PhoD, PhoK, PhoX) or acid (AcpA, PhoN, AphA, NSAPc) pH range. In order to identify their homologs in the genome of Rta, we first selected a set of reference sequences corresponding to bacterial proteins with available 3D structures, one per family (Table 1). In addition, we described each family by a specific set of pre-calculated sequence profiles (PSSMs, for Position-Specific Scoring Matrices). These PSSMs were retrieved from the NCBI Conserved Domain Database (version 3.14) (Marchler-Bauer et al., 2015) according to their similarity to the family reference sequence, using the RPS-Blast program, a standalone version of CD-search (Marchler-Bauer and Bryant, 2004). For each family of phosphatases, we removed the profile redundancy as follows: among the PSSMs whose significant hits on the given reference sequence overlap with a mutual coverage exceeding 70% of the length of the longest hit, only one profile (the one with the lowest E-value) was added to the dataset characterizing the family. The final sets of reference PSSMs are listed in Table 1.
All translated coding sequences (CDS) of Rta, i.e., the putative proteome as provided by the automatic structural annotation of the Rta genome (De Luca et al., 2011), were then compared by RPS-Blast with all selected sets of reference PSSMs. In this functional annotation procedure, a given CDS was annotated as a member of a given phosphatase family when it shared at least one significant hit of each PSSM of the family reference dataset. E-value threshold was set at 1e-10. In order to remove false positive annotations caused by remote homology, i.e., sequences belonging to the same protein superfamily but laying outside any family of interest, validation was further done by blastp (low stringency search, with 1e-3 as E-value threshold) with phosphatase reference sequences as a query.
The entire procedure of functional annotation was also applied to the NCBI set of bacterial genomes completed before March, 2016.
Prediction of Transcription Units and Detection of Regulatory Genomic Sequences
The Rta transcription units encompassing the genes identified as phosphatases using our functional annotation procedure were defined according to DOOR (Mao et al., 2009), a database of computationally predicted operons covering 2072 bacteria genomes. The operon prediction program used to generate this database has been described by Dam et al. (2007) and was ranked to be the best among 14 methods by Brouwer et al. (2008). The prediction accuracy was estimated to be higher than 80% for genomes where only few operons were already known.
Putative transcription factor (TF) binding sites were searched for in the non-coding genomic regions located upstream of the Rta phosphatases CDS, or upstream of the first gene of the operons to which they were predicted to belong. For that purpose, we detected significant similarities between these non-coding genomic sequences and sequence profiles describing the binding specificity of known TFs of E. coli K-12 (86 PSSMs provided by RegulonDB, Huerta et al., 1998), using the program matrix-scan available on the RSAT web server (Turatsinze et al., 2008). Owing to the AT-richness of the prokaryotic regulatory regions whatever the GC content of the genomes (Cordero and Hogeweg, 2009), the large and experimentally validated dataset provided by RegulonDB constituted a reference for the identification of TF binding sites in bacteria. Pairwise hit p-values were estimated according to an organism-specific Markov model of order 1, precalibrated on all the non-coding upstream-gene sequences of Rta and representing their global oligonucleotide composition (i.e., the background model). The p-value threshold was set at 1e-4 (i.e., 1 false positive prediction expected per 10 kb).
Prediction of Protein Subcellular Location and Export Pathway
The subcellular locations of Rta phosphatases were predicted from their amino acid sequences using the program PSORTb v3.0 (Yu et al., 2010) with Gram-negative bacteria models. In addition, the presence of N-terminal signal peptides, the position of their cleavage site, and the corresponding export pathway were investigated by PRED-TAT (Bagos et al., 2010). Theoretical isoelectric point and molecular weight of the predicted precursor and mature protein forms were computed from their respective sequence by ProtParam tool from ExPASy (http://web.expasy.org/protparam).
Molecular Homology Modeling and Analysis of Protein Surface Properties
The protein structures of Rta phosphatases were modeled using the HHpred server (version 5) (Meier and Söding, 2015). The HMM-HMM similarity search was done against the PDB70 database, as available on the 3rd of January 2015. Models based on optimal multiple templates were automatically built using MODELLER (Šali and Blundell, 1993). The quality of the models was checked by Verify3D (Lüthy et al., 1992). Structures and models were superimposed with Matras (Kawabata, 2003). The 3D coordinate files of the protein models are provided in Supplementary Data Sheet 1.
Electrostatic potential maps of 3D structures and models were computed and drawn using the APBS plugin (Baker et al., 2001) in the PyMOL Molecular Graphics System (version 1.8 Schrödinger, LLC), excluding the charges of the active-site metal ions.
The amino acids involved in noncovalent interactions with P ligand and/or metal ions in the reference 3D structures, according to the CCP4 program contact (Winn et al., 2011), were considered as active site residues. After removing ligand, co-factors and putative signal peptide from the 3D coordinate files, POPS server (Cavallo et al., 2003) was used to compute the solvent-accessible surface area (further described as accessibility) of each amino acid. Their sum over the active site residues (or amino acids at homologous positions in protein models) was considered as the active site accessibility.
Rta Phosphatase Activities and Biomineralization
Rta cells were suspended in three different media (CaGP, NaGP, or NaGP+Ca) at a pH of 7.5 with glycerophosphate (GP) as a sole source of phosphorus. Concentrations of dissolved inorganic orthophosphate ions (Pi) and dissolved Ca2+ (Ca) were monitored during the course of the experiment (Figure 1).
Figure 1. Calcification assay performed in different media on Rta cells and the extracellular fraction. Three media were used: (A,B) CaGP (10 mM calcium glycerophosphate and 20 mM HEPES at pH 7.5); (C,D) NaGP+Ca (10 mM of sodium glycerophosphate, 0.75 mM of CaCl2 and 20 mM of HEPES at pH 7.5), and (E,F) NaGP (10 mM sodium glycerophosphate and 20 mM HEPES pH 7.5). The NaGP medium was used as a control assay without Ca. For each assay, time variations of the concentrations of dissolved orthophosphates (A,C,E) and dissolved Ca2+ (B,D,F) were measured. Abiotic controls (labeled Medium) consisted in the media without cells. Error bars (smaller than the symbol size) represent the instrumental error: ±0.85 μM and ±50 μM for dissolved orthophosphates and dissolved calcium, respectively.
Rta cells as well as to a lower extent, Rta extracellular fraction, were able to hydrolyse glycerophosphate (GP), as shown by the increase of Pi concentration over time in all inoculated experiments (Figures 1A,C,E). The increase of Pi concentration was higher for Rta cells and Rta extracellular fraction incubated in the presence of calcium (CaGP and NaGP+Ca, Figures 1A,C) than without Ca (NaGP, Figure 1E), at least for the 9 first days. Moreover, after 9 days, (i) the concentration of dissolved Pi was significantly higher (267 vs. 149 μM) when cells were incubated with 10 mM Ca (CaGP, Figure 1A) than 0.75 mM Ca (NaGP+Ca, Figure 1C); (ii) for Rta cells in CaGP (Figure 1A), Pi concentration plateaued before a decrease down to 116 μM (day 29) due to the precipitation of calcium phosphate. Precipitation in this assay was supported by the net decrease of Ca concentration from 10 mM (CaGP medium, day 0) to 7.4 mM (day 29) (Figure 1B). In contrast, no significant variation of Ca concentration was observed when cells were incubated with 0.75 mM Ca (NaGP+Ca, Figure 1D) or without Ca (NaGP, Figure 1F). Precipitation of Ca phosphates by Rta cells in CaGP was also detected by Fourier-Transform InfraRed (FT-IR) spectroscopy analyses (Figure 2). Precipitates were identified as poorly crystalline hydroxyapatite based on the detection of several characteristics vibrational bands at 472 cm−1 (P-O-P bending ν2), 563 and 602 cm−1 (P-O-P bending ν4), 960–962 cm−1 (P-O stretching ν1), and 1,035–1,045 cm−1 (P-O stretching ν3).
Figure 2. FT-IR spectra of Rta and recombinant E. coli cells in CaGP medium. The FT-IR spectra of E. coli cells (grown in LB medium) and of a reference hydroxyapatite are shown for comparison. The peaks assigned to hydroxyapatite are numbered as follows. 1: 1,035–1,045 cm−1 corresponding to P-O stretching ν3; 2: 960–962 cm−1 corresponding to P-O stretching ν1; 3: 630–633 cm−1 corresponding to librational mode of OH− groups; 4 and 5: 602 and 563 cm−1, respectively, corresponding to P-O-P bending ν4; 6: 472 cm−1 corresponding to P-OP bending ν2.
In contrast, the concentration of dissolved calcium remained constant for cells and extracellular fractions in NaGP+Ca (Figure 1D) as well as extracellular fractions in CaGP (Figure 1B). Based on calculations of the saturation index (SI) of the solutions (Supplementary Table 1), a minimum SI value of 10 was required for precipitation to occur, a condition that was met with Rta cells in CaGP only.
Furthermore, the phosphatase activity of Rta soluble proteins was tested by colorimetric assays in the presence of other organic P substrates: monoester (pNPP), diester (bis-pNPP), and a diester deoxynucleotide analog (TpNPP). A broad phosphatase activity toward the different substrates was detected (Figure 3), the highest rate of reaction being toward monoester and nucleotide analog, corresponding to 1.7 and 1.6 pmol s−1 of orthophosphate released, respectively (1.1 pmol s−1 toward bis-pNPP).
Figure 3. Phosphatase activity of Rta cells toward different substrates. pNPP, p-nitrophenylphosphate (phosphomonoester); bis-pNPP, bis p-nitrophenylphosphate (phosphodiester); TpNPP, Thymidine 5′-monophosphate p-nitrophenylester (phosphodiester, deoxynucleotide analog). Error bars (±0.22 μM) represent the systematic instrumental error.
In Silico Prediction of Rta Non-specific Phosphatases
We identified 5 putative phosphatases in the genome of Rta, according to their sequence similarity with reference profiles and sequences of known phosphatase families (Supplementary Table 2). They belonged to 2 non-homologous families of alkaline phosphatases: 1 PhoD and 4 PhoX, further described as PhoX1-4. No member of other known alkaline phosphatase families, i.e., PhoA and PhoK, and no acid phosphatase homolog were detected. The phoD gene product was already described as a candidate alkaline phosphatase D precursor in the latest annotation release available for the Rta genome from the NCBI. In contrast, no relevant functional annotation was available for the 4 putative PhoX family members, described as (conserved) hypothetical proteins.
The sequences of the 5 putative Rta phosphatases exhibited the same protein domain composition than the reference of their respective families (Supplementary Table 2). No additional phosphatase gene (or pseudogene) was found in the noncoding part of the Rta genome by similarity search of the reference sequences using tblastn.
Genomic Organization, Regulatory Sequences, and Predicted Cellular Location of Rta Phosphatases
The phoD gene was located on the direct DNA strand, while the 4 phoX genes were located on the reverse strand (Table 2). These 5 genes were distant from each other on the genome, i.e., there was no joint organization as gene clusters or operons. The phoX2 gene was predicted as part (3′-end) of an operon in the DOOR database (Mao et al., 2009), together with two other genes: from 5′ (locus tag Rta_37630) to 3′ (Rta_37620), the protein sequences corresponding to these genes were described as (1) uncharacterized (Uniprot accession F5Y2Z8; similar to non-heme di-iron oxygen transport proteins, CDD profile pfam01814) and (2) cytochromes P450-like (Uniprot accession F5Y2Z7), respectively. The 4 other phosphatase-encoding genes were considered as complete transcription units.
In order to predict the regulatory pathways controlling the expression of phosphatases in Rta, and assess whether they were expressed constitutively or under controlled regulation and in the latest case which environmental parameters may impact their transcriptional regulation, we looked for putative regulatory elements. We detected putative cis-regulatory elements that may control the transcription of the phoX3 and phoX4 genes, based on their similarity with known transcription factor (TF) binding sites in E. coli K-12. The significant hits (Table 2) corresponded to 5 TF, the complete (PhoB, Fur, FadR) or partial (Cra, FhlA) homologs of which were identified in the Rta genome by blastp (Supplementary Table 3). No cis-regulatory element was identified in the upstream region of the other putative phosphatase genes, or in the upstream region of the first gene of the phoX2 operon.
While belonging to the same family, the 4 putative PhoX proteins strongly varied in their sequence length, ranging from 458 a.a. (PhoX4) to 756 a.a. (PhoX3) (Table 3). The 5 Rta phosphatases were predicted by PRED-TAT to share a N-terminal signal peptide (length ranging from 31 to 67 a.a.) with specific motifs recognized by the twin-arginine translocase (tat). This bacterial system exports folded proteins to non-cytoplasmic cell space. In parallel, different cell locations were predicted by PSORT-B: extracellular (PhoD), non-cytoplasmic (PhoX2-4) and cytoplasmic (PhoX1). Therefore, there was a conflict between PSORT-B and PRED-TAT predictions for PhoX1, since cytoplasmic proteins do not contain signal peptides.
The theoretical pI of the mature forms of Pho proteins were acid (PhoX1, PhoX3-4), alkaline (PhoX2) and near-neutral (PhoD) ranging from 4.9 to 9.2. Predicted signal sequences did not significantly affect the theoretical pI of Rta phosphatases (Table 3). Finally, the molecular weight of the predicted mature forms ranged from 45.7 (PhoX4) to 73.7 kDa (PhoX3).
Diversity of the Physicochemical Properties of Rta Phosphatases
A molecular model of Rta mature PhoD (499 amino acids) was built by homology with the 3D structure of Bacillus subtilis (Bsu) mature PhoD (527 amino acids), which measured the exact same size as Rta PhoD plus the 28 amino acids of a C-ter cap, not found in Rta PhoD. Both Rta and Bsu mature PhoD shared 45% sequence identity. The root mean square deviation (RMSD) between the C-alpha atoms of the two superposed structures was 0.48 Å (Figure 4A). Based on this model, 9 of the 10 residues of the active site located in the catalytic domain of Bsu PhoD (Figure 4B) were conserved in Rta PhoD, whose active site was composed as follows (sequence numbering): C154, D181, Y184, D235, D236, H237, N241, D242 (N216 in Bsu PhoD, PDB numbering), D425, H427. Two residues of Bsu PhoD active site, labeled R506 and H510 and belonging to a flexible C-terminal extension (C-ter cap, 28 amino acid length), had no equivalent in Rta PhoD (Supplementary Image 1A).
Figure 4. Structural features of Rta PhoD model. (A) Superimposition of 3D model of Rta PhoD and 3D structure of Bsu PhoD used as template (PDB chain 2yeqA). Bsu PhoD: Ig-like domain and catalytic domain are colored in dark blue and cyan, respectively; C-ter cap in gray. Rta PhoD: Ig-like domain (orange), catalytic domain (yellow). (B) 3D view of the active site of Bsu PhoD. The amino acid residues are labeled according to their PDB numbering. Atom name of co-factors and inorganic phosphate molecule are labeled. Hydrogens are hidden. (C,D) Electrostatic potential maps of Bsu PhoD without C-ter cap and Rta PhoD. The surface potentials are contoured from −5 kT/e (red) to 5kT/e (blue). All structures are shown according to their orientation in the superimposition. The positions of ligand and ions are those from Bsu PhoD; their atoms are represented as spheres of given van der Waals radius (A) or small fixed radius (B–D).
Molecular models of mature Rta PhoX1-4 were built based on the 3D experimental structure of Pseudomonas fluorescens (Pfl) PhoX (Figure 5A). This 3D template was also used to generate the molecular model of Sinorhizobium meliloti (Sme) PhoX (Uniprot accession M4MZI0), whose profile of substrate affinity was experimentally determined by Zaheer et al. (2009). The pairwise RMSD and the pairwise sequence identity of Pfl PhoX with each of the PhoX models were: 0.48 Å, 49.5% (Rta PhoX1); 0.46 Å, 49.9% (Rta PhoX2); 0.69 Å, 40.2% (Rta PhoX3); 1.37 Å, 25.4% (Rta PhoX4); 0.62 Å, 41.5% (Sme PhoX). While Rta PhoX1 and PhoX2 exhibited similar sequences (65.7% identity) and structures (0.36 Å RMSD), large insertions and deletions occurred in loops in Rta PhoX3 and PhoX4, respectively. Based on these models, the 10 residues of the active site of Pfl PhoX (Figure 5B) were conserved in the Rta PhoX proteins and also in Sme PhoX (Supplementary Image 1B).
Figure 5. Structural features of Rta PhoX1-4 models. (A) Superimposition of 3D models of Sme PhoX (pink), Rta PhoX1 (green), Rta PhoX2 (orange), Rta PhoX3 (gray), Rta PhoX4 (purple) and 3D structure of Pfl PhoX used as template (PDB chain 4alfA, cyan). (B) 3D view of the active site of Pfl PhoX. The amino acid residues are labeled according to their PDB numbering. Atom name of co-factors and inorganic phosphate molecule are labeled. Hydrogens are hidden. (C–H) Electrostatic potential maps. The surface potentials are contoured from −5 kT/e (red; −15 kT/e for Sme PhoX) to 5kT/e (blue). All structures are shown according to their orientation in the superimposition. The positions of ligand and ions are those from Pfl PhoX; their atoms are represented as spheres of given van der Waals radius (A) or small fixed radius (B–H).
Based on these 3D structures and models, we computed the solvent-accessible surface area of each amino acid located in the protein active site. For a given protein, the sum of these areas over the active site residues (or amino acids at homologous positions in protein models) is called the active site accessibility. This parameter correlates to the protein surface that can accommodate the substrate during the enzymatic process. The accessibility of the active site of the 4 Rta PhoX proteins strongly varied, ranging from small (PhoX2, 213.5 Å2) to large values (PhoX1 and PhoX4, 302.0 and 333.9 Å2, respectively) with an intermediate value for PhoX3 (270.6 Å2) close to the reference Pfl PhoX (268.7 Å2). By comparison, the accessibility of the active site of Sme PhoX was 267.7 Å2. The active site of Rta PhoD exhibited similar accessibility (273.0 Å2) to that of the references Bsu PhoD (261.6 Å2 in an open state, i.e., without its C-ter cap) and Pfl PhoX.
Finally, all 3D structures studied here showed active sites with a negative surface net charge (Figures 4C,D, 5C–H). This likely controls their ability or requirement to interact with metal cations (i.e., co-factors). However, the regions surrounding the active sites exhibited variable patterns of surface electrostatic potential, corresponding to negatively charged patches for Bsu PhoD, Pfl PhoX, Sme PhoX, and Rta PhoX4 or a mosaic of patches with apolar or distinct charged properties for Rta PhoD and Rta PhoX1-X3.
Recombinant Expression of Rta Phosphatases in E. coli; Phosphatase and Calcification Activities
Recombinant expression was achieved for 3 of the 5 Rta phosphatase genes: phoD, phoX1, and phoX3. E. coli phoA was also cloned in the same vector and tested for expression as a control. On SDS-PAGE, an intense band was systematically observed for each transformed strain in the “IPTG” condition that was absent from the “no-IPTG” condition (Supplementary Image 2). These proteins had apparent molecular masses close to those predicted by in silico analyses for the different phosphatases. For E. coli phoA, the overexpression of PhoA was confirmed using Western blot with a commercial anti E. coli alkaline phosphatase antibody (data not shown). Two proteins were clearly visible for PhoA after an induction time of 4 h, corresponding to the expected molecular weights of the precursor and mature forms of this protein. In contrast, no such precursor forms were observed on SDS-PAGE for the Rta PhoD, PhoX1 and PhoX3.
The phoX3 expression was unstable, i.e., the recombinant protein was rapidly degraded over the duration of expression. Recombinant PhoX2 and PhoX4 proteins were toxic in E. coli host cells, i.e., they inhibited E. coli growth. Therefore, the activities of these 3 proteins could not be measured. Phosphatase assays were carried on total soluble proteins extracted from the phoD, phoX1 and phoA E. coli clones (Figure 6). PhoD and PhoX1 hydrolyzed all the pNPP derivatives (pNPP, bis-pNPP, and TpNPP). However, the activity of PhoD was higher for bis-pNPP diester (28.2 pmol s−1 of orthophosphate released) than pNPP monoester (11.1 pmol s−1) and TpNPP nucleotide analog (3.3 pmol s−1), whereas the activity of PhoX1 was higher for monoester (358.6 pmol s−1) than diester (43.1 and 2.3 pmol s−1 for bis-pNPP and TpNPP, respectively).
Figure 6. Enzymatic activities of Rta phosphatases. The soluble protein fractions extracted from E. coli clones overexpressing (A) Rta phoD or (B) Rta phoX1 were assayed for phosphatase activity toward monoester (pNPP), diester (bis-pNPP) or nucleotide analog (TpNPP) substrates. (C) Control assay of the soluble protein fraction of E. coli clones over-expressing E. coli phoA. Error bars (±0.22 μM, smaller than the symbol size) represent the systematic instrumental error.
E. coli cells expressing Rta phoX1 induced hydroxyapatite precipitation in the presence of 10 mM of CaGP as shown by FT-IR spectroscopy (Figure 2), or SEM coupled with energy dispersive x-ray spectrometry (Figure 7). The precipitation of hydroxyapatite induced by E. coli cells expressing E. coli phoA, previously described by Cosmidis et al. (2015), was confirmed. The main peaks in the FT-IR spectrum of the samples could be unambiguously assigned to hydroxyapatite, with broader bands (e.g., P-O stretching ν3 at 1,035–1,045 cm−1) than in the reference hydroxyapatite spectrum, suggesting more disordered structures (Cosmidis et al., 2015). Moreover, the P-O stretching ν1 (960–962 cm−1) was less defined in the “Rta in CaGP” spectrum than in the “Recombinant” spectra, which may indicate a lower crystallinity as well. In contrast, no hydroxyapatite precipitation was evidenced for E. coli cells expressing phoD.
Figure 7. Mineralization of calcium phosphate induced by over-expression of Rta PhoX1 in E. coli clones. (A,B) SEM images in the secondary electron mode. (A) Image of a cluster of cells encrusted by calcium phosphate. (B) Close-up showing the lumen of a cell in the center. The wall is encrusted by calcium phosphate precipitates. (C) Energy dispersive x-ray spectrometry (EDXS) spectrum measured on the mineral phase precipitated on the cells. Emission lines of O, P, and Ca correspond to hydroxyapatite. The emission line of C is due to the cells and the carbon filter on which the sample was deposited.
Broad Substrate Recognition by Rta Phosphatases
Rta cells showed a significant phosphatase activity toward a broad spectrum of substrates. While the highest reaction rate was observed for a phosphate monoester and for the nucleotide-analog phosphodiester (TpNPP), there was also a significant hydrolysis of bis-pNPP, a phosphate phosphodiester molecule smaller than TpNPP. Five genes of non-specific phosphatases were identified in Rta genome, belonging to the two non-homologous families phoD (1 gene) and phoX (4 genes, phoX1-4). Among the 5 phosphatases, phoD and phoX1 heterologously expressed in E. coli were found to catalyze the hydrolysis of different substrates, including diesters at the exception of the nucleotide analog.
In agreement with our observations, PhoD family members have previously been shown to hydrolyze phosphomonoesters and phosphodiesters with a preference for the latter (Rodriguez et al., 2014).
PhoX1 showed a high phosphomonoesterase activity. Therefore, Rta PhoX1 profile for organic P substrate (mono- and diesters) was intermediate between the monoester-specific activity of P. fluorescens PhoX (Yong et al., 2014), and the affinity of S. meliloti PhoX for most of the metabolites with phosphoester bonds, including nucleotides, carbohydrates and amino acids (Zaheer et al., 2009). The 3D model of the S. meliloti PhoX (Figure 5D) revealed (i) an active site accessibility similar to P. fluorescens PhoX but smaller than Rta PhoX1, and (ii) a surface with a high negative net charge, while P. fluorescens PhoX and Rta PhoX1 both exhibited a surface with a low positive net charge.
While sharing a conserved active site, PhoX proteins showed different patterns of predicted surface charges and accessibility of the active site. The substrates used in the present study corresponded to diverse shapes and occupancies, and were negatively charged: −4 for TpNPP, −3 for bis-pNPP, and −2 for pNPP and GP. It has been shown that positive charges at the surface of the phosphatases may facilitate the loading of negatively charged substrates into their catalytic site (Coleman, 1992). The diversity of the substrate profiles observed within the PhoX family thus seems to rely more on the physicochemical properties of the amino acids surrounding the active site and little on the active site residues themselves.
The capability of Rta to hydrolyze nucleotide-analog phosphodiesters was therefore likely due to one of the 3 other phosphatases identified in the genome. The ability of recombinant PhoX2 and PhoX4 to hydrolyze the nucleotide phosphodiester bonds can be discussed, although it could not be verified experimentally. These proteins were encoded by genes with a complete structure and without genomic features related to pseudogenes. According to 3D models, PhoX4 had an active site with a very large solvent accessibility surface (~330 Å2) and therefore appeared as a good candidate for the accommodation of large substrates such as nucleotides. The predicted high positive net charges at the surface of PhoX2 (Figure 5F), together with the predicted Tat signal, which should induce the expression of a folded protein in the cytoplasm (Teter and Klionsky, 1999), may explain the toxicity of this protein, which may catalyze DNA degradation when overexpressed in E. coli. PhoX2 was organized into a 3-genes transcription unit (operon). This genomic organization suggests that the proteins encoded by these genes may form a macromolecular assembly when expressed altogether in Rta (Wells et al., 2016). In this macromolecular assembly, the catalytic site of PhoX2 may be buried, precluding DNA damage within Rta cells. In the present study, Rta phoX2 only was cloned in E. coli without the two other genes of the predicted operon. Moreover, no homolog of these two genes were found in the genome of E. coli. As a result, the catalytic site of PhoX2 would not be buried in a macromolecular assembly in E. coli and would therefore be toxic to the cells. However, we did not find any report in the literature about such an assembly between homologs of these three proteins and therefore further work will be needed to assess the above assumption.
Overall, the diversity of phosphatases in Rta allows the hydrolysis of a broad range of organic P substrates and therefore the release of Pi under diverse trophic conditions. Since Pi release is key to the achievement of high saturation levels with hydroxyapatite and the induction of phosphatogenesis, Rta appears as a particularly efficient driver of this process as shown experimentally. Hydroxyapatite precipitation was shown to first occur in the periplasm of Rta cells, then in the cytoplasm (Benzerara et al., 2004). The prediction of a non-cytoplasmic location for most of the phosphatases is consistent with these observations, except PhoX1 for which location remains unclear.
Rta Harbors an Unusual Set of Diverse Phosphatases
The diversity of alkaline phosphatases in microorganisms was previously investigated at the community level by metagenomics, with a particular focus on the relationship between gene abundance and the availability of inorganic phosphorus in ocean surface waters (Luo et al., 2009; Martiny et al., 2009; Sebastian and Ammerman, 2009; Temperton et al., 2011) or in soils (Lidbury et al., 2017; Neal et al., 2017). Overall, PhoX and PhoD families were evidenced to be more abundant than PhoA in high-pH soils and under marine oligotrophic conditions. Interestingly, (i) the Tataouine sand from which Rta was isolated showed similar P-depleted and Ca-rich conditions (Gillet et al., 2000), and (ii) the Rta genomic repertoire contained phosphatases of both the PhoX and PhoD families. These data support the existence of a relationship between the availability of Ca and P in the environment and the relative abundance of the phosphatase gene families. However, the taxonomic scope of studies on ocean surface waters was limited by the numerical domination of microbial phototrophs, in particular Cyanobacteria, and heterotrophic Alpha- and Gamma-proteobacteria in these environments.
Here, in order to limit the bias due to a limited taxonomic sampling, we annotated the phosphatase genes in the large amount of available bacterial complete genomes. This dataset comprised 3630 genomes from 33 phyla. Eight families, either alkaline or acid phosphatases, were annotated according to the procedure previously applied to Rta genome. Compared to the statistical trends that emerged from the analysis of bacterial genomes, the functional potential of the Rta genome for phosphatases exhibited atypical features: (i) among the 12 bacterial genomes with at least 4 PhoX genes, Rta was the only member of Proteobacteria, a phylum which was represented by 1794 genomes; (ii) while 62% of the bacterial genomes had at least one gene of alkaline phosphatase (71% of the 261 Betaproteobacteria genomes), only 4% (1% for Betaproteobacteria) had at least 5 genes, similarly as Rta; (iii) although it was absent in Rta, the PhoA family was the most widely represented with 35% of the bacterial genomes exhibiting at least one gene of this family. It was less than 25% for each other family, either alkaline or acid, e.g., 22 and 17% in the case of PhoX and PhoD, respectively. Among the genomic features of bacteria, the linear correlation between the genome size and the global number of genes (abundance) was well known (Hou and Lin, 2009) (Supplementary Image 3; R2 = 0.97). Though the abundance of phosphatase genes was related to the genome size in our dataset (Kruskal-Wallis test, p < 2.2e-16), Rta represented the bacterial genome with the highest density of PhoX genes (Supplementary Image 4). In addition, only 0.9% of the bacterial genomes had a higher density of alkaline phosphatase genes than Rta. The variable physicochemical properties of Rta phosphatases suggested a functional diversification of the assimilation of organic P substrates, driven by P-depleted conditions.
Environmental Parameters Involved in the Regulation of the Phosphatase Activity of Rta
Our results highlighted the influence of several chemical parameters on the phosphatase activity of Rta.
(1) First, phosphomonoester hydrolysis by Rta cells was favored by extracellular Ca2+. Accordingly, the requirement of calcium ions for phosphatase activity was previously described for PhoD or PhoX proteins (Zaheer et al., 2009; Rodriguez et al., 2014; Yong et al., 2014). The molecular modeling of the Rta phosphatases revealed the high-level of conservation of the active site residues within PhoD and PhoX families. Based on this high-level of conservation of the active site residues, the cofactors of the reference proteins, i.e., Ca2+ and Fe3+, were predicted as probable cofactors of the 5 Rta phosphatases. Although there is no experimental or thermodynamic confirmation of this prediction based on sequence conservation, it is interesting to note that Ca was abundant in the sandy soil from which Rta was isolated (Gillet et al., 2000). Calcium can thus have a double role in phosphatogenesis: (i) by directly participating to the precipitation of Ca-phosphates but also (ii) indirectly by affecting the activity of some phosphatases. In the first case, the extracellular Ca matters, while intracellular Ca plays a role in the second case. Calcium is a key regulator in eukaryotes and bacteria that can modify the cell state and the pattern of gene expression (Dominguez, 2004; Domínguez et al., 2011). Therefore, the relationship between Ca availability in the environment and the cell induction of phosphatogenesis may overall be complex and depend on processes regulating Ca homeostasis within cells.
(2) In addition, putative cis-regulatory elements similar to known transcription factor binding sites in E. coli K-12 were found in Rta genome in the upstream region of the PhoX3 and PhoX4 genes. In E. coli K-12, these transcription factors are involved in the regulation of the expression of genes related to (i) P uptake and metabolism (pho regulon), controlled by the extracellular concentration of inorganic P (PhoB) (Hsieh and Wanner, 2010; Santos-Beneit, 2015), (ii) ferric iron homeostasis (Fur) (da Silva Neto et al., 2009; Porcheron et al., 2013), and (iii) fatty acid and phospholipid metabolism (FadR) (Iram and Cronan, 2005; Fujita et al., 2007). We detected complete homologs of these transcription factors in Rta genome. Previously, the alkaline phosphatases of different families have been experimentally identified as members of the pho regulon, in particular in the Gammaproteobacteria E. coli (Torriani, 1960; Wanner, 1990), S. meliloti (Bardin and Finan, 1998; Yuan et al., 2006) and P. fluorescens (Monds et al., 2006), as well as in B. subtilis (Eder et al., 1996; Hulett, 1996; Prágai et al., 2004), Mycobacterium tuberculosis (Torres et al., 2001) and Synechocystis sp. PCC 6803 (Suzuki et al., 2004). However, the predicted relationship between iron or fatty acids/phospholipids and non-specific phosphatases at the transcription level may be a specificity of Rta. Such a connection between iron uptake, fatty acids/phospholipid metabolism and phosphatase expression may be due to the fact that these phosphatases need Fe as a co-factor and use phospholipids as substrates. However, the ability of Rta phosphatases to use phospholipids as organic P substrates in P-depleted conditions will have to be evaluated. Interestingly, in addition to being predicted as a probable cofactor of the Rta phosphatases, Fe was also present in the environment where Rta was isolated from, although mobilization under oxic conditions likely requires the presence of organic chelatants (Benzerara et al., 2005, 2006). Moreover, fatty acid metabolism has been shown to increase slightly extracellular pH around cells and impact CaCO3 biomineralization (Marvasi et al., 2010). Together with Pi release, the same pH effect may play a synergistic positive role on phosphatogenesis, which will have to be further explored experimentally. Overall, the profile of P substrates hydrolyzed by Rta cells may thus be tuned by the expression of a varying combination of phosphatases depending on the environmental parameters.
In conclusion, the diverse set of phosphatases found in Rta suggests that this strain can hydrolyze diverse P-containing organic molecules under diverse environmental conditions. Since the resulting production of orthophosphates can induce the precipitation of metal-phosphates under some conditions, Rta may therefore be particularly prone to phosphatogenesis. Such a relationship between phosphatase diversity within a bacterial strain and its role in phosphatogenesis has not been investigated before. However, the approach developed here combining in silico and in vitro analyses sets the foundations for future studies of more complex environmental communities with the aim of inferring their phosphatogenesis capabilities.
FS-P, KB, and ED designed the experiments, ED and GC performed bioinformatics analyses, KB and JC performed SEM analyses, JC and FS-P acquired and interpreted FT-IR spectra, FS-P and CF performed microbial cultures, heterologous expression, and assays, GD and TH provided Rta strain and contributed to microbial cultures. All authors contributed to manuscript writing and revision.
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 CNRS-INSU grants (programs EC2CO and INTERRVIE) to ED. GC internship was funded by the French Cluster of Excellence MATISSE. The Scanning Electron Microscope (SEM) facility at IMPMC was supported by Region Ile-de-France grant SESAME 2006 I-07-593/R, INSU-CNRS, INP-CNRS, University Pierre et Marie Curie – UPMC. Technical support from Maxime Guillaumet (FT-IR spectroscopy) at IMPMC was greatly appreciated.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2017.02592/full#supplementary-material
Bagos, P. G., Nikolaou, E. P., Liakopoulos, T. D., and Tsirigos, K. D. (2010). Combined prediction of Tat and Sec signal peptides with hidden Markov models. Bioinformatics 26, 2811–2817. doi: 10.1093/bioinformatics/btq530
Bailey, J. V., Corsetti, F. A., Greene, S. E., Crosby, C. H., Liu, P., and Orphan, V. J. (2013). Filamentous sulfur bacteria preserved in modern and ancient phosphatic sediments: implications for the role of oxygen and bacteria in phosphogenesis. Geobiology 11, 397–405. doi: 10.1111/gbi.12046
Baker, N. A., Sept, D., Joseph, S., Holst, M. J., and McCammon, J. A. (2001). Electrostatics of nanosystems: application to microtubules and the ribosome. Proc. Natl. Acad. Sci. U.S.A. 98, 10037–10041. doi: 10.1073/pnas.181342398
Beazley, M. J., Martinez, R. J., Sobecky, P. A., Webb, S. M., and Taillefert, M. (2007). Uranium biomineralization as a result of bacterial phosphatase activity: insights from bacterial isolates from a contaminated subsurface. Environ. Sci. Technol. 41, 5701–5707. doi: 10.1021/es070567g
Beazley, M. J., Martinez, R. J., Sobecky, P. A., Webb, S. M., and Taillefert, M. (2009). Nonreductive biomineralization of uranium(VI) phosphate via microbial phosphatase activity in anaerobic conditions. Geomicrobiol. J. 26, 431–441. doi: 10.1080/01490450903060780
Benzerara, K., Chapon, V., Moreira, D., López-García, P., Guyot, F., and Heulin, T. (2006). Microbial diversity on the Tatahouine meteorite. Meteorit. Planet. Sci. 41, 1249–1265. doi: 10.1111/j.1945-5100.2006.tb00519.x
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., Menguy, N., Guyot, F., Skouri, F., De Luca, G., Barakat, M., et al. (2004). Biologically controlled precipitation of calcium phosphate by Ramlibacter tataouinensis. Earth Planet. Sci. Lett. 228, 439–449. doi: 10.1016/j.epsl.2004.09.030
Benzerara, K., Menguy, N., Guyot, F., Vanni, C., and Gillet, P. (2005). TEM study of a silicate-carbonate-microbe interface prepared by focused ion beam milling. Geochim. Cosmochim. Acta 69, 1413–1422. doi: 10.1016/j.gca.2004.09.008
Bihani, S. C., Das, A., Nilgiriwala, K. S., Prashar, V., Pirocchi, M., Apte, S. K., et al. (2011). X-ray structure reveals a new class and provides insight into evolution of alkaline phosphatases. PLoS ONE 6:e22767. doi: 10.1371/journal.pone.0022767
Blake, R. E., O'Neil, J. R., and Garcia, G. A. (1998). Effects of microbial activity on the δ18O of dissolved inorganic phosphate and textural features of synthetic apatites. Am. Mineral. 83, 1516–1531. doi: 10.2138/am-1998-11-1240
Calderone, V., Forleo, C., Benvenuti, M., Thaller, M. C., Rossolini, G. M., and Mangani, S. (2006). A structure-based proposal for the catalytic mechanism of the bacterial acid phosphatase AphA belonging to the DDDD superfamily of phosphohydrolases. J. Mol. Biol. 355, 708–721. doi: 10.1016/j.jmb.2005.10.068
Cavallo, L., Kleinjung, J., and Fraternali, F. (2003). POPS: a fast algorithm for solvent accessible surface areas at atomic and residue level. Nucleic Acids Res. 31, 3364–3366. doi: 10.1093/nar/gkg601
Coleman, M. L., and Chisholm, S. W. (2010). Ecosystem-specific selection pressures revealed through comparative population genomics. Proc. Natl. Acad. Sci. U.S.A. 107, 18634–18639. doi: 10.1073/pnas.1009480107
Cosmidis, J., Benzerara, K., Gheerbrant, E., Estève, I., Bouya, B., and Amaghzaz, M. (2013a). Nanometer-scale characterization of exceptionally preserved bacterial fossils in Paleocene phosphorites from Ouled Abdoun (Morocco). Geobiology 11, 139–153. doi: 10.1111/gbi.12022
Cosmidis, J., Benzerara, K., Guyot, F., Skouri-Panet, F., Duprat, E., Férard, C., et al. (2015). Calcium-phosphate biomineralization induced by alkaline phosphatase activity in Escherichia coli: localization kinetics and potential signatures in the fossil record. Front. Earth Sci. 3:84. doi: 10.3389/feart.2015.00084
Cosmidis, J., Benzerara, K., Menguy, N., and Arning, E. (2013b). Microscopy evidence of bacterial microfossils in phosphorite crusts of the Peruvian shelf: implications for phosphogenesis mechanisms. Chem. Geol. 359, 10–22. doi: 10.1016/j.chemgeo.2013.09.009
Cosmidis, J., Benzerara, K., Morin, G., Busigny, V., Lebeau, O., Jézéquel, D., et al. (2014). Biomineralization of iron-phosphates in the water column of Lake Pavin (Massif Central, France). Geochim. Cosmochim. Acta 126, 78–96. doi: 10.1016/j.gca.2013.10.037
da Silva Neto, J. F., Braz, V. S., Italiani, V. C. S., and Marques, M. V. (2009). Fur controls iron homeostasis and oxidative stress defense in the oligotrophic alpha-proteobacterium Caulobacter crescentus. Nucleic Acids Res. 37, 4812–4825. doi: 10.1093/nar/gkp509
De Luca, G., Barakat, M., Ortet, P., Fochesato, S., Jourlin-Castelli, C., Ansaldi, M., et al. (2011). The cyst-dividing bacterium Ramlibacter tataouinensis TTB310 genome reveals a well-stocked toolbox for adaptation to a desert environment. PLoS ONE 6:e23784. doi: 10.1371/journal.pone.0023784
Eder, S., Shi, L., Jensen, K., Yamane, K., and Hulett, F. M. (1996). A Bacillus subtilis secreted phosphodiesterase/alkaline phosphatase is the product of a Pho regulon gene phoD. Microbiology 142, 2041–2047. doi: 10.1099/13500872-142-8-2041
Felts, R. L., Reilly, T. J., Calcutt, M. J., and Tanner, J. J. (2006a). Cloning, purification and crystallization of Bacillus anthracis class C acid phosphatase. Acta Crystallogr. Sect. F 62(Pt 7), 705–708. doi: 10.1107/S174430910602389X
Felts, R. L., Reilly, T. J., and Tanner, J. J. (2006b). Structure of Francisella tularensis AcpA: prototype of a unique superfamily of acid phosphatases and phospholipases C. J. Biol. Chem. 281, 30289–30298. doi: 10.1074/jbc.M606391200
Gillet, P., Barrat, J. A., Heulin, T., Achouak, W., Lesourd, M., Guyot, F., et al. (2000). Bacteria in the Tatahouine meteorite: nanometric-scale life in rocks. Earth Planet. Sci. Lett. 175, 161–167. doi: 10.1016/S0012-821X(99)00291-5
Handley-Sidhu, S., Renshaw, J. C., Yong, P., Kerley, R., and Macaskie, L. E. (2010). Nano-crystalline hydroxyapatite bio-mineral for the treatment of strontium from aqueous solutions. Biotechnol. Lett. 33, 79–87. doi: 10.1007/s10529-010-0391-9
Heulin, T., Barakat, M., Christen, R., Lesourd, M., Sutra, L., De Luca, G., et al. (2003). Ramlibacter tataouinensis gen. nov. sp. nov. and Ramlibacter henchirensis sp. nov. cyst-producing bacteria isolated from subdesert soil in Tunisia. Int. J. Syst. Evol. Microbiol. 53, 589–594. doi: 10.1099/ijs.0.02482-0
Hirschler, A., Lucas, J., and Hubert, J.-C. (1990). Apatite genesis: a biologically induced or biologically controlled mineral formation process? Geomicrobiol. J. 8, 47–56. doi: 10.1080/01490459009377877
Hou, Y., and Lin, S. (2009). Distinct gene number-genome size relationships for eukaryotes and non-eukaryotes: gene content estimation for dinoflagellate genomes. PLoS ONE 4:e6978. doi: 10.1371/journal.pone.0006978
Kathuria, S., and Martiny, A. C. (2011). Prevalence of a calcium-based alkaline phosphatase associated with the marine cyanobacterium Prochlorococcus and other ocean bacteria. Environ. Microbiol. 13, 74–83. doi: 10.1111/j.1462-2920.2010.02310.x
Lidbury, I. D. E. A., Fraser, T., Murphy, A. R. J., Scanlan, D. J., Bending, G. D., Jones, A. M. E., et al. (2017). The 'known' genetic potential for microbial communities to degrade organic phosphorus is reduced in low-pH soils. MicrobiologyOpen e474. doi: 10.1002/mbo3.474
Macaskie, L. E., Bonthrone, K. M., Yong, P., and Goddard, D. T. (2000). Enzymically mediated bioprecipitation of uranium by a Citrobacter sp.: a concerted role for exocellular lipopolysaccharide and associated phosphatase in biomineral formation. Microbiology 146, 1855–1867. doi: 10.1099/00221287-146-8-1855
Makde, R. D., Mahajan, S. K., and Kumar, V. (2007). Structure and mutational analysis of the PhoN protein of Salmonella typhimurium provide insight into mechanistic details. Biochemistry 46, 2079–2090. doi: 10.1021/bi062180g
Marchler-Bauer, A., Derbyshire, M. K., Gonzales, N. R., Lu, S., Chitsaz, F., Geer, L. Y., et al. (2015). CDD: NCBI's conserved domain database. Nucleic Acids Res. 43, D222–D226. doi: 10.1093/nar/gku1221
Martin, P., Dyhrman, S. T., Lomas, M. W., Poulton, N. J., and Van Mooy, B. A. S. (2014). Accumulation and enhanced cycling of polyphosphate by sargasso sea plankton in response to low phosphorus. Proc. Natl. Acad. Sci. U.S.A. 111, 8089–8094. doi: 10.1073/pnas.1321719111
Martinez, R. J., Beazley, M. J., Taillefert, M., Arakaki, A. K., Skolnick, J., and Sobecky, P. A. (2007). Aerobic uranium (VI) bioprecipitation by metal-resistant bacteria isolated from radionuclide- and metal-contaminated subsurface soils. Environ. Microbiol. 9, 3122–3133. doi: 10.1111/j.1462-2920.2007.01422.x
Martiny, A. C., Huang, Y., and Li, W. (2009). Occurrence of phosphate acquisition genes in Prochlorococcus cells from different ocean regions. Environ. Microbiol. 11, 1340–1347. doi: 10.1111/j.1462-2920.2009.01860.x
Marvasi, M., Visscher, P. T., Perito, B., Mastromei, G., and Casillas-Martínez, 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
Miot, J., Benzerara, K., Morin, G., Bernard, S., Beyssac, O., Larquet, E., et al. (2009). Transformation of vivianite by anaerobic nitrate-reducing iron-oxidizing bacteria. Geobiology 7, 373–384. doi: 10.1111/j.1472-4669.2009.00203.x
Mondani, L., Benzerara, K., Carrière, M., Christen, R., Mamindy-Pajany, Y., Février, L., et al. (2011). Influence of uranium on bacterial communities: a comparison of natural uranium-rich soils with controls. PLoS ONE 6:e25771. doi: 10.1371/journal.pone.0025771
Monds, R. D., Newell, P. D., Schwartzman, J. A., and O'Toole, G. A. (2006). Conservation of the Pho regulon in Pseudomonas fluorescens Pf0-1. Appl. Environ. Microbiol. 72, 1910–1924. doi: 10.1128/AEM.72.3.1910-1924.2006
Moorehead, W. R., and Biggs, H. G. (1974). 2-amino-2-methyl-1-propanol as the alkalizing agent in an improved continuous-flow cresolphthalein complexone procedure for calcium in serum. Clin. Chem. 20, 1458–1460.
Neal, A. L., Rossmann, M., Brearley, C., Akkari, E., Guyomar, C., Clark, I. M., et al. (2017). Land-use influences phosphatase gene microdiversity in soils. Env Microbiol. 19, 2740–2753. doi: 10.1111/1462-2920.13778
Nilgiriwala, K. S., Alahari, A., Rao, A. S., and Apte, S. K. (2008). Cloning and overexpression of alkaline phosphatase PhoK from Sphingomonas sp. strain BSAR-1 for bioprecipitation of uranium from alkaline solutions. Appl. Environ. Microbiol. 74, 5516–5523. doi: 10.1128/AEM.00107-08
Porcheron, G., Garénaux, A., Proulx, J., Sabri, M., and Dozois, C. M. (2013). Iron, copper, zinc and manganese transport and regulation in pathogenic enterobacteria: correlations between strains site of infection and the relative importance of the different metal transport systems for virulence. Front. Cell. Infect. Microbiol. 3:90. doi: 10.3389/fcimb.2013.00090
Powers, L. G., Mills, H. J., Palumbo, A. V., Zhang, C., Delaney, K., and Sobecky, P. A. (2002). Introduction of a plasmid-encoded phoA gene for constitutive overproduction of alkaline phosphatase in three subsurface Pseudomonas isolates. FEMS Microbiol. Ecol. 41, 115–123. doi: 10.1111/j.1574-6941.2002.tb00972.x
Prágai, Z., Allenby, N. E. E., O'Connor, N., Dubrac, S., Rapoport, G., Msadek, T., et al. (2004). Transcriptional regulation of the phoPR operon in Bacillus subtilis. J. Bacteriol. 186, 1182–1190. doi: 10.1128/JB.186.4.1182-1190.2004
Ragot, S. A., Kertesz, M. A., Mészáros, É., Frossard, E., and Bünemann, E. K. (2017). Soil phoD and phoX alkaline phosphatase gene diversity responds to multiple environmental factors. FEMS Microbiol. Ecol. 93:fiw212. doi: 10.1093/femsec/fiw212
Rivadeneyra, M. A., Martín-Algarra, A., Sánchez-Román, M., Sánchez-Navas, A., and Martín-Ramos, J. D. (2010). Amorphous Ca-phosphate precursors for Ca-carbonate biominerals mediated by Chromohalobacter marismortui. ISME J. 4, 922–932. doi: 10.1038/ismej.2010.17
Rodriguez, F., Lillington, J., Johnson, S., Timmel, C. R., Lea, S. M., and Berks, B. C. (2014). Crystal structure of the Bacillus subtilis phosphodiesterase PhoD reveals an iron and calcium-containing active site. J. Biol. Chem. 289, 30889–30899. doi: 10.1074/jbc.M114.604892
Sebastián, M., Pitta, P., González, J. M., Thingstad, T. F., and Gasol, J. M. (2012). Bacterioplankton groups involved in the uptake of phosphate and dissolved organic phosphorus in a mesocosm experiment with P-starved Mediterranean waters. Environ. Microbiol. 14, 2334–2347. doi: 10.1111/j.1462-2920.2012.02772.x
Shelobolina, E. S., Konishi, H., Xu, H., and Roden, E. E. (2009). U(VI) sequestration in hydroxyapatite produced by microbial glycerol 3-phosphate metabolism. Appl. Environ. Microbiol. 75, 5773–5778. doi: 10.1128/AEM.00628-09
Sowell, S. M., Wilhelm, L. J., Norbeck, A. D., Lipton, M. S., Nicora, C. D., Barofsky, D. F., et al. (2008). Transport functions dominate the SAR11 metaproteome at low-nutrient extremes in the Sargasso Sea. ISME J. 3, 93–105. doi: 10.1038/ismej.2008.83
Suzuki, S., Ferjani, A., Suzuki, I., and Murata, N. (2004). The SphS-SphR two component system is the exclusive sensor for the induction of gene expression in response to phosphate limitation in Synechocystis. J. Biol. Chem. 279, 13234–13240. doi: 10.1074/jbc.M313358200
Temperton, B., Gilbert, J. A., Quinn, J. P., and McGrath, J. W. (2011). Novel analysis of oceanic surface water metagenomes suggests importance of polyphosphate metabolism in oligotrophic environments. PLoS ONE 6:e16499. doi: 10.1371/journal.pone.0016499
Templeton, A. S., Trainor, T. P., Spormann, A. M., Newville, M., Sutton, S. R., Dohnalkova, A., et al. (2003). Sorption versus biomineralization of Pb(II) within Burkholderia cepacia biofilms. Environ. Sci. Technol. 37, 300–307. doi: 10.1021/es025972g
Torres, A., Juárez, M. D., Cervantes, R., and Espitia, C. (2001). Molecular analysis of Mycobacterium tuberculosis phosphate specific transport system in Mycobacterium smegmatis. Characterization of recombinant 38 kDa (PstS-1). Microb. Pathog. 30, 289–297. doi: 10.1006/mpat.2001.0434
Turatsinze, J.-V., Thomas-Chollier, M., Defrance, M., and van Helden, J. (2008). Using RSAT to scan genome sequences for transcription factor binding sites and cis-regulatory modules. Nat. Protoc. 3, 1578–1588. doi: 10.1038/nprot.2008.97
Wanner, B. L. (1990). “Phosphorus assimilation and its control of gene expression in Escherichia coli,” in The Molecular Basis of Bacterial Metabolism, eds G. Hauska, and R. K. Thauer (Berlin; Heidelberg: Springer), 152–163.
Winn, M. D., Ballard, C. C., Cowtan, K. D., Dodson, E. J., Emsley, P., Evans, P. R., et al. (2011). Overview of the CCP4 suite and current developments. Acta Cryst. D67, 235–242. doi: 10.1107/S0907444910045749
Wojciechowski, C. L., Cardia, J. P., and Kantrowitz, E. R. (2002). Alkaline phosphatase from the hyperthermophilic bacterium T. maritima requires cobalt for activity. Protein Sci. 11, 903–911. doi: 10.1110/ps.4260102
Yong, S. C., Roversi, P., Lillington, J., Rodriguez, F., Krehenbrink, M., Zeldin, O. B., et al. (2014). A complex iron-calcium cofactor catalyzing phosphotransfer chemistry. Science 345, 1170–1173. doi: 10.1126/science.1254237
Young, C. L., and Ingall, E. D. (2010). Marine dissolved organic phosphorus composition: insights from samples recovered using combined electrodialysis/reverse osmosis. Aquat. Geochem. 16, 563–574. doi: 10.1007/s10498-009-9087-y
Yu, N. Y., Wagner, J. R., Laird, M. R., Melli, G., Rey, S., Lo, R., et al. (2010). PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics 26, 1608–1615. doi: 10.1093/bioinformatics/btq249
Yuan, Z.-C., Zaheer, R., Morton, R., and Finan, T. M. (2006). Genome prediction of PhoB regulated promoters in Sinorhizobium meliloti and twelve proteobacteria. Nucleic Acids Res. 34, 2686–2697. doi: 10.1093/nar/gkl365
Yung, M. C., and Jiao, Y. (2014). Biomineralization of uranium by PhoY phosphatase activity aids cell survival in Caulobacter crescentus. Appl. Environ. Microbiol. 80, 4795–4804. doi: 10.1128/AEM.01050-14
Zaheer, R., Morton, R., Proudfoot, M., Yakunin, A., and Finan, T. M. (2009). Genetic and biochemical properties of an alkaline phosphatase PhoX family protein found in many bacteria. Environ. Microbiol. 11, 1572–1587. doi: 10.1111/j.1462-2920.2009.01885.x
Zalatan, J. G., Fenn, T. D., Brunger, A. T., and Herschlag, D. (2006). Structural and functional comparisons of nucleotide pyrophosphatase/phosphodiesterase and alkaline phosphatase: implications for mechanism and evolution. Biochemistry 45, 9788–9803. doi: 10.1021/bi060847t
Keywords: microbial phosphatases, biomineralization, phosphatogenesis, metal-phosphate mineral phases, enzymatic activity, genomic repertoire, organic phosphates, hydroxyapatite
Citation: Skouri-Panet F, Benzerara K, Cosmidis J, Férard C, Caumes G, De Luca G, Heulin T and Duprat E (2018) In Vitro and in Silico Evidence of Phosphatase Diversity in the Biomineralizing Bacterium Ramlibacter tataouinensis. Front. Microbiol. 8:2592. doi: 10.3389/fmicb.2017.02592
Received: 29 June 2017; Accepted: 12 December 2017;
Published: 11 January 2018.
Edited by:David Emerson, Bigelow Laboratory For Ocean Sciences, United States
Reviewed by:Clara S. Chan, University of Delaware, United States
Eric D. van Hullebusch, IHE Delft Institute for Water Education, Netherlands
Copyright © 2018 Skouri-Panet, Benzerara, Cosmidis, Férard, Caumes, De Luca, Heulin and Duprat. 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: Elodie Duprat, email@example.com