ORIGINAL RESEARCH article
Metabolomic Analysis Reveals Contributions of Citric and Citramalic Acids to Rare Earth Bioleaching by a Paecilomyces Fungus
- 1Department of Civil and Environmental Engineering, University of California, Berkeley, Berkeley, CA, United States
- 2Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, United States
- 3Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, United States
Conventional methods for extracting rare earth elements from monazite ore require high energy inputs and produce environmentally damaging waste streams. Bioleaching offers a potentially more environmentally friendly alternative extraction process. In order to better understand bioleaching mechanisms, we conducted an exo-metabolomic analysis of a previously isolated rare earth bioleaching fungus from the genus Paecilomyces (GenBank accession numbers KM874779 and KM 874781) to identify contributions of compounds exuded by this fungus to bioleaching activity. Exuded compounds were compared under two growth conditions: growth with monazite ore as the only phosphate source, and growth with a soluble phosphate source (K2HPO4) added. Overall metabolite profiling, in combination with glucose consumption and biomass accumulation data, reflected a lag in growth when this organism was grown with only monazite. We analyzed the relationships between metabolite concentrations, rare earth solubilization, and growth conditions, and identified several metabolites potentially associated with bioleaching. Further investigation using laboratory prepared solutions of 17 of these metabolites indicated statistically significant leaching contributions from both citric and citramalic acids. These contributions (16.4 and 15.0 mg/L total rare earths solubilized) accounted for a portion, but not all, of the leaching achieved with direct bioleaching (42 ± 15 mg/L final rare earth concentration). Additionally, citramalic acid released significantly less of the radioactive element thorium than did citric acid (0.25 ± 0.01 mg/L compared to 1.18 ± 0.01 mg/L), suggesting that citramalic acid may have preferable leaching properties for a monazite bioleaching process.
Rare earth elements (REEs) are critical for a variety of technologies (US DOE, 2011; Alonso et al., 2012). Monazite, a rare earth phosphate mineral, is one of the main ores used for commercial REE production (Gupta and Krishnamurthy, 1992; Rosenblum and Fleischer, 1995). However, conventional monazite extraction methods involve chemical leaching at high temperatures, and are energy intensive and environmentally damaging (Merritt, 1990; Gupta and Krishnamurthy, 1992; Alonso et al., 2012; Peelman et al., 2014). Additionally, the radioactive element thorium (Th) is also usually present along with REE-phosphates in monazite ore (Rosenblum and Fleischer, 1995). Conventional monazite extraction methods co-extract Th, which must be separated from REEs in downstream processing (Gupta and Krishnamurthy, 1992).
Bioleaching offers a possible alternative to conventional REE extraction, potentially resulting in a more environmentally sustainable extraction process (Brisson et al., 2016; Fathollahzadeh et al., 2018a). Several previous studies have investigated the potential for bioleaching of REEs from monazite by phosphate solubilizing microorganisms (PSMs), including both fungi (genera: Aspergillus, Paecilomyces, and Penicillium) and bacteria (genera: Acidithiobacillus, Bacillus, Burkholderia, Enterobacter, Klebsiella, Microbacterium, Pantoea, Pseudomonas, and Streptomyces) (Brisson et al., 2016; Corbett et al., 2017, 2018; Fathollahzadeh et al., 2018b, c, 2019). Monazite bioleaching by some fungi has also been shown to preferentially release REEs into solution over Th (Brisson et al., 2016). Several recent studies have investigated bioleaching REEs from substrates other than monazite. This includes leaching REEs from bastnaesite, another important REE bearing ore (Zhang et al., 2018), and bioleaching to recover REEs from waste products including spent catalysts, rare earth magnets, and coal fly ash (Auerbach et al., 2019; Jin et al., 2019; Park and Liang, 2019).
Phosphate solubilizing microorganisms are microorganisms that have the ability to solubilize phosphate ions from otherwise insoluble phosphate compounds and minerals (Rodríguez and Fraga, 1999). Current understanding of the mechanisms of phosphate solubilization by PSMs indicates that three main contributing factors are acidification of the medium, exchange reactions, and the formation of complexes between organic acids or other chelating molecules produced by the PSMs and cations associated with phosphate in the mineral and released during solubilization (Bolan et al., 1994; Rodríguez and Fraga, 1999; Nautiyal et al., 2000; Gyaneshwar et al., 2002; Arcand and Schneider, 2006; Scervino et al., 2011). Previous work on monazite bioleaching indicated that although both acidification and complexation with citric acid were capable of contributing to monazite leaching, these contributions did not account for the levels of leaching seen during bioleaching or when leaching with spent bioleaching medium (Brisson et al., 2016). This indicates that other compounds released into solution, but not identified by the methods used in that study, contribute to the bioleaching process.
Mass spectrometry based metabolomics technologies provide the opportunity to accurately detect a large number of different organic molecules and compare relative concentrations across different conditions, providing insight into biological processes. Metabolomic analyses applied to exuded metabolites are sometimes referred to as exometabolomics or metabolic footprinting (Kell et al., 2005). Metabolic footprinting has been applied to investigate other eukaryotic microbial processes including wine production and microalgae growth in bioreactors (Howell et al., 2006; Sue et al., 2011; Richter et al., 2013). Recent metabolomic studies of fungi have investigated a variety of important processes including the production of mycotoxins, fungal degradation of biomass waste, and metabolic changes and metabolite exchange in fungal-bacterial mutualistic interactions (Karpe et al., 2015, 2018; Shi et al., 2017; Uehling et al., 2019). One recent study did a comparative genomic and exometabolomic analysis of 27 members of Aspergillus section Nigri, a group of fungi with important roles in biotechnology, the environment, and human health (Vesth et al., 2018).
A metabolomics approach could potentially enable better understanding of the bioleaching process by identifying a larger array of small molecules released during bioleaching that might be associated with bioleaching effectiveness. In this study, we analyzed metabolites exuded into the growth medium during monazite bioleaching under two different growth conditions: growth with monazite as the only phosphate source (using soluble phosphate limitation to force monazite solubilization) and growth with the addition of a soluble phosphate source (relieving the phosphate limitation stress). This analysis had two parallel goals. One was to examine the effects of phosphate availability on growth and metabolic processes of a bioleaching microorganism, and the second was to identify metabolites exuded into solution that may contribute to monazite solubilization. Once those metabolites of interest were identified, the ability of laboratory prepared solutions of those metabolites to leach monazite under abiotic conditions was studied to further investigate individual metabolites’ contributions to bioleaching.
Materials and Methods
Organism and Bioleaching Growth Conditions
Bioleaching experiments were conducted with a monazite bioleaching fungal isolate designated WE3-F (GenBank accession numbers KM874779 and KM874781), whose isolation and identification as a Paecilomyces species were described previously (Brisson et al., 2016). This organism was selected for further study based on its consistent bioleaching performance in that study.
Growth conditions were based on those described previously (Brisson et al., 2016) with some modifications. Briefly, bioleaching was conducted in 250 mL glass Erlenmeyer flasks, each containing 0.5 g ground monazite sand (City Chemical LLC, West Haven, CT, United States) (finer than 200 mesh) and 50 mL modified ammonium salts medium (AMS medium) (Parales et al., 1994). AMS medium contained 1.0 g/L MgSO4⋅7H2O, 0.2 g/L KCl, 0.66 g/L (NH4)2SO4, 1.0 mL/L 1000× trace elements stock solution, and 1.0 mL/L stock A. The 1000× trace elements stock solution contained 0.5 g/L FeSO4⋅7H2O, 0.4 g/L ZnSO4⋅7H2O, 0.02 g/L MnSO4⋅H2O, 0.015 g/L H3BO3, 0.01 g/L NiCl2⋅6H2O, 0.25 g/L EDTA, 0.05 g/L CoCl2⋅6H2O, and 0.005 g/L CuCl2⋅2H2O. Stock A contained 5 g/L FeNaEDTA and 2 g/L NaMoO4⋅2H2O. 10 g/L glucose was added as a carbon and energy source and air in the headspace served as oxygen source. Each flask was inoculated with 1 mL of spore suspension containing approximately 107 CFU and sealed with a foam stopper. As in previous bioleaching studies with this organism (Brisson et al., 2016), flasks were stirred continuously at 250 rpm to ensure well mixed conditions, and incubated at 28°C for the duration of the bioleaching experiment.
Two different growth conditions were compared to study the effects of a soluble phosphate source: growth with monazite only and growth with K2HPO4 and monazite. For each K2HPO4 and monazite condition flask, 0.4 g/L K2HPO4 was added.
Quantification of REEs, Th, Phosphate, Glucose, pH, and Biomass
Rare earth element, Th, phosphate, glucose, pH, and biomass were quantified by previously described analytical methods (Brisson et al., 2016). Briefly, REE and Th concentrations were measured using an Agilent Technologies 7700 series ICP-MS. Phosphate concentration was measured using the BioVision Phosphate Colorimetric Assay Kit. Glucose was measured by HPLC on a Waters 2695 HPLC system with a BioRad Aminex HPX-87H carbohydrate/organic acids analysis column and a Waters 2414 refractive index detector. pH was measured using a Hanna Instruments HI 2210 pH meter. Biomass was measured as total volatile solids of filter-collected samples by drying at 105°C and subsequent ashing at 550°C as described previously (Brisson et al., 2016) based on United States Environmental Protection Agency Method 1684 (US EPA, 2001). REE, Th, phosphate, glucose, and pH measurements were taken for six biological replicates for each time point (0, 2, 4, and 6 days after inoculation), while biomass measurements were taken for three biological replicates at time points 2, 4, and 6 days.
Samples of bioleaching supernatant were collected, filtered through 0.2 μm syringe filters to remove cells, and immediately frozen and stored at −80°C. Six replicate samples were collected at each time point (0, 2, 4, and 6 days after inoculation). The t = 0 days samples served as a media controls, since they consisted of medium that had not undergone incubation with fungi. Subsequent samples represent metabolite production (and consumption) during the bioleaching process. Metabolomic analysis was performed by the West Coast Metabolomics Center at the University of California, Davis.
Details of the protocols used for metabolite extraction, derivatization, gas chromatography-mass spectrometry (GC-MS) data acquisition, and data processing have been described in detail previously (Fiehn et al., 2008; Fiehn, 2016). Briefly, the samples were extracted with a 3:3:2 solution of acetonitrile, isopropanol, and water. Extracted samples were derivatized via silylation derivatization with N-Methyl-N-(trimethylsilyl)trifluoroacetamide (MSTFA). Derivatized samples were analyzed on an Agilent 6890 GC coupled with a Pegasus TOF MS. The GC column used was an Rtx-5Sil GC column (30 cm long, 0.25 mm internal diameter). Compound identifications were made based on their retention index and comparisons of their mass spectra to the BinBase metabolomics database using the BinBase algorithm (Skogerson et al., 2011).
Hierarchical clustering (Müllner, 2011) of metabolites based on concentration profiles was performed in Python using the SciPy cluster module. Signal intensity data for each metabolite were first centered by subtracting the mean signal intensity for that metabolite, and normalized by dividing by the standard deviation. Hierarchical clustering was performed using the “complete” method, also called the farthest point algorithm, with Euclidian distances (Müllner, 2011).
Identification of Metabolites of Potential Bioleaching Importance
Metabolites that were potentially relevant to bioleaching performance were identified by three methods. The first method identified metabolites that were released at higher concentrations under the monazite-only condition than under the K2HPO4 plus monazite condition. Signal intensities for each metabolite were compared between the two conditions using a two-tailed T-test for independent samples. p-values were corrected for multiple comparisons using the Benjamini/Hochberg correction for false discovery rate for independent samples (Benjamini and Hochberg, 1995). For this analysis only, all metabolites for which the p-values were marginally significant (p < 0.1) were selected for further study. This less stringent p-value criterion was used at this intermediate stage in order to identify a large number of metabolites for the final set of experiments. This analysis was performed independently for time points 2, 4, and 6 days.
The second approach to selecting metabolites of interest was to identify correlations between metabolite concentration (signal intensity) and REE concentration over all time points. This analysis was performed on data from the monazite-only condition, using measurements of metabolite concentrations and REE concentrations at each time point. A least squares linear regression was performed to identify correlations. p-values were corrected for multiple comparisons using the Šidák correction (Šidák, 1967). Metabolites for which the linear regression had a positive slope and a significant corrected p-value (p > 0.05) were selected for further study.
The final approach was to select metabolites with the highest signal intensities. Metabolites for which the average signal intensity was greater than 105 for any condition and time point were selected for further study.
Abiotic Leaching With Identified Metabolites
Abiotic leaching conditions were conducted as previously described (Brisson et al., 2016), with some modifications. Leaching was conducted in 50 mL flat bottomed polypropylene tubes, each containing 0.1 g ground monazite sand (200 mesh). 10 mL leaching solution was added to autoclaved tubes and stirred for 48 h at 250 rpm at room temperature (25–28°C). All leaching solutions were tested in triplicate.
Leaching solutions contained selected metabolites at a concentration of 10 mM, with the exception of stearic acid. This concentration was selected based on observed concentrations of organic acids identified in our previous study of monazite bioleaching (Brisson et al., 2016). Stearic acid, which has an extremely low solubility in water (0.003 g/L or 0.01 mM at 20°C) (Anneken et al., 2000), was dissolved in water for 20 min with vortexing and filtered to remove undissolved particles. Additionally, a combined leaching solution containing all selected metabolites, each at a concentration of 10 mM (except for stearic acid), was also tested. All leaching solutions were adjusted to pH 2.5 by the addition of HCl in order to mimic the pH observed during bioleaching and to eliminate the effects of variations in pH observed previously (Brisson et al., 2016). By making this adjustment, contributions of individual metabolites could be quantified based on their ability to increase leaching beyond what was expected for that pH level. Leaching solutions were filter sterilized through 0.2 μm syringe filters prior to leaching experiments.
Statistical significance of leaching effectiveness was determined using a two-tailed T-test for independent samples to compare REEs released by each leaching solution to a control solution of HCl at a pH of 2.5. p-values were corrected for multiple comparisons using the Šidák correction (Šidák, 1967).
Rare earth element solubilization was greater for the monazite-only condition, when a soluble phosphate source was not provided, reaching concentrations of 42 ± 15 mg/L total REEs after 6 days of leaching (Figure 1A). However, some solubilization of REEs did occur in the cultures provided with K2HPO4, reaching concentrations of 14 ± 9 mg/L after 6 days of bioleaching. Although release of the radioactive element Th was low for both conditions, it was consistently greater for the monazite-only condition (0.6 ± 0.3 mg/L) than for K2HPO4 plus monazite (0.04 ± 0.02 mg/L) (Figure 1B).
Figure 1. Bioleaching of monazite in the absence or presence of soluble phosphate (K2HPO4). Shown are (A) REE concentrations, (B) Th concentrations, (C) phosphate concentrations, (D) glucose concentrations, (E) pH, and (F) biomass measured as volatile solids. REE, phosphate, glucose, and pH data are for six biological replicates. Biomass data are for three biological replicates. Error bars indicate standard deviations around the mean.
Free phosphate concentrations (Figure 1C) remained very low (maximum observed concentration in a single sample: 0.005 mM) when monazite was the only phosphate source. When K2HPO4 was added to the medium, phosphate levels decreased from their initial concentration but remained high throughout the experiment (minimum observed concentration in a single sample: 0.68 mM). This indicates that the concentration of K2HPO4 provided was sufficient to avoid phosphate limiting conditions during bioleaching for this growth condition.
Glucose consumption (Figure 1D) for the monazite-only growth condition lagged behind glucose consumption when K2HPO4 was provided. Similarly, the pH was reduced at a faster rate when soluble phosphate was provided (Figure 1E), resulting in a slightly lower pH for this condition on day two of bioleaching despite the higher initial pH of the medium with added K2HPO4. However, by the fourth day, both conditions had achieved a similar pH. Biomass production (Figure 1F) for the monazite-only condition also lagged behind growth with K2HPO4 plus monazite. By the sixth day, however, biomass accumulation was comparable under both growth conditions (2.8 ± 0.03 g/L for monazite-only and 2.9 ± 0.2 g/L for K2HPO4 with monazite).
Overall Metabolomic Profile
Metabolomic analyses of the fungal supernatant from the two conditions detected 210 metabolites. Of these 87 could be identified as known compounds. The remaining 123 were identified only with BinBase ID numbers based on their characteristic mass spectra, referring to the BinBase metabolomics database (Skogerson et al., 2011). Details of the compound identifications and mass spectra are given in Supplementary Table 1, and peak heights for all identified compounds in all samples are given in Supplementary Table 2. Concentration profiles of all metabolites identified by chemical name are summarized in Figure 2 for all time points and conditions. Concentration profiles for all 210 metabolites are summarized in Supplementary Figure 1.
Figure 2. Concentration profiles of identified metabolites. Heatmap shows average levels of identified metabolites detected during monazite bioleaching for each growth condition and time point. Columns represent different conditions and time points. Rows represent different metabolites. Metabolites are ordered based on hierarchical clustering, with the clustering dendrogram displayed to the left of the heatmap and metabolite names to the right. Heatmap colors indicate standard deviations below (blue) and above (yellow) the overall mean level for each metabolite. Letters between heatmap and metabolite names indicate KEGG pathways involving each metabolite. E = energy metabolism, C = carbohydrate metabolism, A = amino acid metabolism, L = lipid metabolism, N = nucleic acid metabolism, V = metabolism of cofactors and vitamins, O = other KEGG pathways. Metabolites highlighted in green are part of the TCA cycle. Metabolites highlighted in yellow are disaccharides and alcohols detected early in the growth cycle. Metabolites highlighted in purple are long chain (>8 carbon) fatty acids.
Nine identified compounds were detected at significant levels in the medium control (t = 0 days) samples (Figure 2, bottom nine rows). These included glucose, which was added to the medium as a carbon source. Because these were present in the medium, and did not increase with bioleaching, they were not considered to be fungal metabolites.
The concentration profiles of fungal metabolites were consistent with the lag in growth when monazite was the only phosphate source, observed above in glucose consumption, pH reduction, and biomass growth (see section “Bioleaching Performance,” Figures 1D–F). Early exuded metabolites (detected at day two for growth with K2HPO4 and day for monazite only) included several disaccharides (cellobiose, lactulose, and β-gentiobiose) and alcohols (lactitol, lyxitol, and glycerol) (Figure 2, yellow highlights).
The lag in growth for the monazite only condition was paralleled in the concentration time profiles for most metabolites involved in the tricarboxylic acid cycle (TCA cycle) (Figure 2, green highlights), including citric, isocitric, α-ketoglutaric, succinic, fumaric, and malic acids. The concentrations of these TCA cycle components peaked on day four when K2HPO4 was added and on day six when monazite was the only phosphate source. Aconitic acid, also in the TCA cycle, had a different concentration profile, peaking on day six for both growth conditions, while oxaloacetic acid was not detected.
A large group of metabolites (Figure 2, middle rows), were present at the highest levels on day six when K2HPO4 was provided, and remained at much lower levels when monazite was the only phosphate source. These included several long chain (more than eight carbons) fatty acids including azelaic, capric, lauric, oleic, palmitic, pelargonic, and stearic acids (Figure 2, purple highlights).
Identification of Metabolites of Potential Bioleaching Importance
Metabolites of potential bioleaching importance were identified with three separate approaches: higher concentrations when soluble phosphate was not available (Figure 3), correlation with REE concentrations (Figure 4), and high signal intensity. Of those that could be identified by chemical name, a subset were selected for further study (Table 1).
Figure 3. Metabolites of potential bioleaching importance identified by occurrence at higher concentrations in the monazite-only condition compared to K2HPO4 plus monazite. (A) Metabolites identified at 2 days. (B) Metabolites identified for at 4 days. Signal intensities are normalized as the fraction of the maximum observed signal intensity for each metabolite at that time point. Height and error bars indicate mean and standard deviation for six biological replicates. Metabolites selected based on a marginally significant Benjamini/Hochberg adjusted p-value < 0.1. This less stringent p-value criterion was used at this intermediate stage only, in order to identify a large number of metabolites for the abiotic leaching experiments.
Figure 4. Correlations between metabolite signal intensities and REE concentrations for all samples from the monazite only condition. Only compounds found to have a significant positive correlation and having identification beyond BinBase ID numbers are shown. p-values are Šidák adjusted.
Metabolites Released at Higher Concentrations Without Soluble Phosphate
Direct comparison of metabolite levels for the two growth conditions identified metabolites with higher concentrations for the monazite-only growth condition (Figure 3). This analysis identified 15 and 13 metabolites for the 2 and 4 day time points, respectively. Three metabolites were identified for both time points. However, none of these three had identification beyond BinBase ID numbers (20282, 2044, and 1681). No metabolites were identified from the analysis at the 6 day time point because differences in concentration were not found to be statistically significant, likely due to the high variability at this time point.
Of the eleven metabolites identified by name (eight for 2 days and three for 4 days) seven (ribose, ribitol, nicotinic acid, isothreonic acid, gluconic acid, histidine, and citric acid) were selected for further study of their leaching abilities using laboratory prepared solutions of these metabolites to leach abiotically. Sulfuric acid was not considered because the focus of this analysis was organic metabolites, and because the acidification of the medium, which contained sulfate, produces sulfuric acid. Glucose was not considered because it was the provided substrate rather than a metabolite and its higher concentration in the monazite-only condition at 4 days was already shown (Figure 1D). Fructose was also rejected because our previous study of monazite bioleaching with this organism tested fructose as a carbon source found that it did not have any benefits over glucose in REE solubilization (Brisson et al., 2016). Ribonic acid was not used in the further experiment because of the prohibitive cost and low quantities that could be purchased commercially. This not only made it unsuitable for further testing, but also makes it a less desirable target for a potential leaching technology.
Metabolites With Concentrations That Correlated With REE Concentrations
Fifteen metabolites, eight of which were identified by chemical name, were found to have positive correlations with REE concentrations (Figure 4). Six of these (galactinol, citramalic acid, 4-hydroxybenzoic acid, 3,4-dihydroxybenzoic acid, 1-deoxyerythritol, 2-deoxyerythritol) were selected for further leaching studies, while the other two (erythronic acid and dehydroabietic acid) were rejected for the same reasons at ribonic acid above.
High Signal Intensity Metabolites
Seven named metabolites had high overall signal intensities (average signal intensity >105 for at least one condition and time point). Four of these compounds (sorbitol, glycerol, p-cresol, and stearic acid) were selected for further study, while three (glucose, sulfuric acid, and phosphate), were rejected for previously stated reasons.
Abiotic Leaching Effectiveness of Identified Metabolites
Of the 17 tested metabolites, two (citric acid and citramalic acid) showed statistically significant improvements in REE solubilization greater than the pH 2.5 HCl control, indicating additional leaching contributions beyond what was expected for pH effects (p = 0.008 and 0.04 after Šidák correction for citric and citramalic acid respectively) (Figure 5A). Leaching with a combination of all selected metabolites did not improve solubilization significantly beyond the combined effects of individual metabolites, with increases of only approximately 6.5 and 5.1 mg/L above controls (9.9 mg/L for pH = 2.5 control) for citric and citramalic acids, respectively, and did not approach the REE concentrations achieved by direct bioleaching (42 ± 15 mg/L, see Figure 1A). With regard to Th release during bioleaching, only citric acid, citramalic acid, and the combination of all selected metabolites resulted in detectable levels of Th release, with citric acid releasing significantly more Th than citramalic acid (1.18 ± 0.01 mg/L as opposed to 0.25 ± 0.01 mg/L) (Figure 5B).
Figure 5. Abiotic solubilization of REEs and Th from monazite by selected metabolites. Heights and error bars indicate means and standard deviations for three replicates. The ∗ symbol indicates a statistically significant increase in leaching compared to the control based on an adjusted p-value < 0.05. (A) REE solubilization. (B) Th solubilization.
This study investigated exuded metabolites as potential contributors to bioleaching of REEs from monazite by a Paecilomyces fungus. A previous study found that spent bioleaching medium from this fungus could effectively leach REEs from monazite, indicating that exuded compounds contribute to bioleaching (Brisson et al., 2016). However, the organic acids detected and tested in that study did not account for the leaching capacity of the spent medium. In the study presented here, we employed a metabolomics analysis to identify a larger set of exuded compounds of potential bioleaching interest, and tested their leaching capabilities.
This study focused on the impacts of soluble phosphate availability on the exometabolome of this bioleaching fungus during bioleaching. This comparison was chosen because the fungus may regulate the exudation of bioleaching relevant compounds in response to the need to acquire otherwise unavailable phosphate. Thus, only the soluble phosphate availability was varied, while the presence and amount of monazite was kept constant in order to exclude confounding impacts of the presence or absence of the monazite itself and focus on phosphate availability. Monazite may have other impacts on the fungal exometabolome, potentially in response to the REEs themselves, the radioactive element Th, or the physical presence of the ground sand particles. While these effects could be interesting, they were outside the scope and focus of this study.
The absence of soluble phosphate, and reliance on monazite as the sole phosphate source, resulted in a lag in fungal growth, but was necessary for effective leaching of REEs. By the end of the bioleaching experiment, our data indicate that phosphate was not the limiting factor for growth. This is supported by the comparable glucose, pH, and biomass levels achieved under both growth conditions by the sixth day of bioleaching (Figure 1). REE and Th solubilization were much higher when monazite was the only phosphate source. This is consistent both with forcing the organisms to solubilize phosphate from monazite for growth and with possible re-precipitation of REE-PO4 in the medium that contains K2HPO4 at a relatively high phosphate content. The solubilities of REE phosphates are known to be extremely low (approximately 10–13 M) (Firsching and Brune, 1991).
The depletion of glucose by the end of the experiment under both growth conditions may suggest a glucose growth limitation. However, in a previous study with the same organism (Brisson et al., 2016), increasing the glucose concentration to 100 g/L did not improve bioleaching performance. Nitrogen availability is another possible growth limiting factor. Some estimates of the nitrogen content of fungal mycelia range from 0.2 to 9% of dry weight (Lahoz et al., 1966; Dawson et al., 1989; Watkinson et al., 2006). Assuming a typical value of 5% nitrogen content, the 0.66 g/L of (NH4)2SO4 (i.e., 0.14 g/L N) provided in AMS medium would correspond to the production of approximately 2.8 g/L of dry biomass. This is comparable to the biomass produced under both growth conditions in this study, suggesting that nitrogen availability may be limiting biomass production to this level. Scervino et al. (2011) found nitrogen limitation to enhance phosphate mineral solubilization by Penicillium purpurogenum, another phosphate solubilizing fungus, indicating that nitrogen limitation of growth may be desirable for bioleaching performance. Nitrogen limitation has also been found to enhance citric acid production in some fungi (Cunningham and Kuiack, 1992; Papagianni, 2007).
Hierarchical clustering (Müllner, 2011) of metabolites based on their concentration profiles identified groups of metabolites with similar responses to the two growth conditions and progression over the course of fungal growth. As the fungus progressed from exponential growth (day two for growth with K2HPO4, days two and four for monazite only), through a transition (day four for growth with K2HPO4, day six for monazite only), and into stationary phase (day six for growth with K2HPO4) (see Figure 1F), there were distinct groups of metabolites detected at each stage (Figure 2). Metabolites that clustered together and were detected at higher concentrations in the earlier stages include several disaccharides (cellobiose, lactulose, and β-gentiobiose) and alcohols (lactitol, lyxitol, and glycerol) (Figure 2, yellow highlights).
Most of the components of the TCA cycle (Figure 2, green highlights) clustered together, with concentrations peaking near the transition from exponential growth to stationary phase (day four with K2HPO4 and day six for monazite only). This overall trend is consistent with the earlier depletion of glucose when K2HPO4 was provided (Figure 1G) since glucose, through glycolysis and the TCA cycle, feeds into the production of these metabolites (Madigan et al., 2008). Once the glucose is depleted, these TCA cycle components are consumed and not replenished, resulting in the reduced concentrations by day six when K2HPO4 is provided (Figure 2). The importance of maintaining high sugar concentrations for the production and excretion of citric acid by Aspergillus niger, a commercially important production process, have been well documented (Magnuson and Lasure, 2004; Papagianni, 2007).
Long chain fatty acids (longer than eight carbons) (Figure 2, purple highlights) also clustered together, with concentrations peaking on day six when K2HPO4 was provided, and remaining at much lower levels when monazite was the only phosphate source. Long chain fatty acids observed were azelaic, capric, lauric, oleic, palmitic, pelargonic, and stearic acids. This increase in long chain fatty acid production corresponds with the depletion of glucose and the leveling off in biomass production (Figures 1D,F), and may be related to a transition from exponential growth to stationary phase. Long chain fatty acids and their derivatives have been associated with changes in fungal physiology and morphology, and specifically with the transition from growth to spore formation (Mysyakina and Feofilova, 2011).
Analysis of individual metabolite levels identified several metabolites potentially correlated with bioleaching, and further experiments with 17 of these metabolites showed that citric and citramalic acids contribute significantly to monazite leaching. A combination of all 17 metabolites together did not increase solubilization beyond the expected contributions of the individual metabolites, indicating that there were no synergistic effects of the combination of these metabolites.
Previous studies have shown that citric acid was more effective that oxalic, phthalic, salicylic, gluconic, itaconic, succinic, and acetic acids at leaching REEs from monazite (Goyne et al., 2010; Brisson et al., 2016). Although the effect of citric acid appears to be somewhat larger here than that reported in our previous study of monazite leaching (6.5 mg/L here as opposed to 3 mg/L in our previous study) (Brisson et al., 2016), the experimental protocols were quite different (see Materials and methods). This experiment supports the overall result from that study, indicating that citric acid provides some additional REE solubilization, although not sufficient improvements to account for the majority of the bioleaching effectiveness.
To our knowledge, REE solubilization by citramalic acid has not been previously reported. One study found that phosphate stress induced increased levels if citramailc acid exudation from the roots of sugar beet plants (Khorassani et al., 2011). That study also showed that citramalic acid solubilized phosphate from low phosphate soils amended with monocalcium phosphate dihydrate (Ca(H2PO4)2⋅H2O) (Khorassani et al., 2011).
Notably, the same metabolites that contributed to additional REE solubilization also contributed to Th solubilization. However, citric acid leached significantly more Th than citramalic acid, indicating that citramalic acid may have more desirable leaching characteristics for preferential leaching of REEs. Th, which is radioactive, commonly co-occurs in monazite with REEs-phosphates, and can present concerns for downstream processing and waste disposal in monazite extraction (Gupta and Krishnamurthy, 1992).
The ability of citric and citramalic acids to leach REEs, as well as their differential levels of Th leaching, make them interesting targets for development of an industrial bioleaching process. Fungi are known to produce a wide array of organic acids, and are used industrially for the production of citric, gluconic, itaconic, kojic, and oxalic acids (Liaud et al., 2014; Yang et al., 2017; Hossain et al., 2019). Of the 87 metabolites identified in this study, 37 were organic acids (Supplementary Table 1), and eight of the 17 metabolites selected for abiotic leaching tests were organic acids. Only two metabolites, citric and citramalic acid, significantly increased leaching. This result indicates that the particular organic acids, rather that organic acids in general, are important for effective leaching. As another example, oxalic acid, which is also produced by fungi, forms insoluble complexes with REEs and causes them to precipitate out of solution (Brisson et al., 2016). The differences in Th solubilization between citric and citramalic acid further emphasizes the importance of the particular acids present. A recent study describing a citramalic acid biosynthesis pathway in A. niger showed that overexpression of two genes, cimA and mfsB, reduced citric acid production and increased citramalic acid production (Hossain et al., 2019). This type of control could be useful for developing a process to preferentially recover REEs without the release of Th.
This study investigated 17 metabolites that were potentially associated with bioleaching and could be identified and tested for their leaching abilities. However, a large proportion of detected metabolites (123 out of 210) could not be identified, and 21 of these unidentified metabolites were either correlated with REE concentration or were exuded at higher concentrations when soluble phosphate was not provided. Some of these unknown compounds may also contribute to bioleaching activity, but without identifications their contributions could not be investigated in this study. However, as mass spectrometry databases continue to expand, it may be possible to identify these compounds in the future based on their retention indices and mass spectra (Supplementary Table 1).
In addition to the metabolites detected by the methods used in this study, other products exuded by the fungus, such as siderophores or extracellular enzymes, may also contribute to bioleaching. Although generally studied in the context of iron complexation, siderophores have also been shown to form complexes with REEs, and could also contribute to bioleaching (Bau et al., 2013; Osman et al., 2019). Siderophore production and protein excretion by the strain used in this study have not been characterized. However, the related species Paecilomyces variotii has been shown to produce multiple siderophores (Vala et al., 2000). P. variotii is also known to produce a variety of extracellular enzymes, including phytases, which breaks down phytic acid by removing phosphate groups (de Laguna et al., 2015).
This work has advanced our understanding of the monazite bioleaching process by a Paecilomyces fungus. The metabolomic analysis identified the impacts of soluble vs. insoluble phosphate availability on the dynamics of metabolite exudation. Based on that analysis, we identified and tested a set of metabolites of potential bioleaching importance. We confirmed the contribution of citric acid and identified an additional contribution by citramalic acid to the bioleaching process. Further, we showed that citramalic acid exhibits preferable leaching properties to those of citric acid due to the reduced Th release.
Data Availability Statement
The GC-MS metabolomics data have been deposited in the Global Natural Product Social Molecular Networking (GNPS) data repository with the MassIVE (Mass Spectrometry Interactive Virtual Environment) accession number MSV000084613.
All authors designed and planned the experiments, contributed to the discussion of the results, and wrote the manuscript. VB conducted the experiments and analyzed the data.
This research was supported by the Siemens Corporate Research, a division of Siemens Corporation, through award number UCB_CKI-2012-Industry_IS-001-Doyle. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.
Conflict of Interest
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.
The handling Editor declared a shared affiliation, though no other collaboration, with the author LA-C at the time of review.
We thank Dr. Negassi Hadgu for assistance in the laboratory, and the West Coast Metabolomics Center for their help with the metabolomics analysis. Work presented here originally appeared as part of VB’s doctoral dissertation (Brisson, 2015), which is available online at https://escholarship.org/uc/item/742647jc.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.03008/full#supplementary-material
Alonso, E., Sherman, A. M., Wallington, T. J., Everson, M. P., Field, F. R., Roth, R., et al. (2012). Evaluating rare earth element availability: a case with revolutionary demand from clean technologies. Environ. Sci. Technol. 46, 3406–3414. doi: 10.1021/es203518d
Anneken, D. J., Both, S., Christoph, R., Fieg, G., Steinberner, U., and Westfechtel, A. (2000). Fatty Acids, in Ullmann’s Encyclopedia of Industrial Chemistry. Weinheim: Wiley-VCH Verlag GmbH & Co. KGaA, 73–116.
Arcand, M. M., and Schneider, K. D. (2006). Plant- and microbial-based mechanisms to improve the agronomic effectiveness of phosphate rock: a review. An. Acad. Bras. Ciênc. 78, 791–807. doi: 10.1590/S0001-37652006000400013
Auerbach, R., Bokelmann, K., Stauber, R., Gutfleisch, O., Schnell, S., and Ratering, S. (2019). Critical raw materials – Advanced recycling technologies and processes: recycling of rare earth metals out of end of life magnets by bioleaching with various bacteria as an example of an intelligent recycling strategy. Miner. Eng. 134, 104–117. doi: 10.1016/j.mineng.2018.12.022
Bau, M., Tepe, N., and Mohwinkel, D. (2013). Siderophore-promoted transfer of rare earth elements and iron from volcanic ash into glacial meltwater, river and ocean water. Earth Planet. Sci. Lett. 364, 30–36. doi: 10.1016/j.epsl.2013.01.002
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Bolan, N. S., Naidu, R., Mahimairaja, S., and Baskaran, S. (1994). Influence of low-molecular-weight organic acids on the solubilization of phosphates. Biol. Fertil. Soils 18, 311–319. doi: 10.1007/bf00570634
Corbett, M. K., Eksteen, J. J., Niu, X.-Z., Croue, J.-P., and Watkin, E. L. J. (2017). Interactions of phosphate solubilizing microorganisms with natural rare-earth phosphate minerals: a study utilizing Western Australian monazite. Bioprocess Biosyst. Eng. 40, 929–942. doi: 10.1007/s00449-017-1757-3
Corbett, M. K., Eksteen, J. J., Niu, X.-Z., and Watkin, E. L. J. (2018). Syntrophic effect of indigenous and inoculated microorganisms in the leaching of rare earth elements from Western Australian monazite. Res. Microbiol. 169, 558–568. doi: 10.1016/j.resmic.2018.05.007
Dawson, M. W., Maddox, I. S., and Brooks, J. D. (1989). Evidence for nitrogen catabolite repression during citric acid production by Aspergillus niger under phosphate-limited growth conditions. Biotechnol. Bioeng. 33, 1500–1504. doi: 10.1002/bit.260331119
de Laguna, I. H. B., Marante, F. J. T., and Mioso, R. (2015). Enzymes and biproducts produces by the ascomycete fungus Paecilomyces variotii. J. Appl. Microbiol. 119, 1455–1466. doi: 10.1111/jam.12934
Fathollahzadeh, H., Eksteen, J. J., Kaksonen, A. H., and Watkin, E. L. J. (2018a). Role of bioleaching of rare earth elements from primary and secondary resources. Appl. Microbiol. Biotechnol. 103, 1043–1057. doi: 10.1007/s00253-018-9526-z
Fathollahzadeh, H., Hackett, M. J., Khaleque, H. N., Eksteen, J. J., Kaksonen, A. H., and Watkin, E. L. J. (2018b). Better together: potential of co-culture microorganisms to enhance bioleaching of rare earth elements from monazite. Bioresour. Technol. Rep. 3, 109–118. doi: 10.1016/j.biteb.2018.07.003
Fathollahzadeh, H., Hackett, M. J., Khaleque, H. N., Eksteen, J. J., Kaksonen, A. H., and Watkin, E. L. J. (2018c). Microbial contact enhances bioleaching of rare earth elements. Bioresour. Technol. Rep. 3, 102–108. doi: 10.1016/j.biteb.2018.07.004
Fathollahzadeh, H., Khaleque, H. N., Eksteen, J., Kaksonen, A. H., and Watkin, E. L. J. (2019). Efffect of glycine on bioleaching of rare earth elements from Western Australian monazite by heterotrophic and autotrophic microorganisms. Hydrometallurgy 189, 105137. doi: 10.1016/j.hydromet.2019.105137
Fiehn, O., Wohlgemuth, G., Scholz, M., Kind, T., Lee, D. Y., Moon, S., et al. (2008). Quality control for plant metabolomics: reporting MSI-compliant studies. Plant J. 53, 691–704. doi: 10.1111/j.1365-313X.2007.03387.x
Howell, K. S., Cozzolino, D., Bartowsky, E. J., Fleet, G. H., and Henschke, P. A. (2006). Metabolic profiling as a tool for revealing Saccharomyces interactions during wine fermentation. FEMS Yeast Res. 6, 91–101. doi: 10.1111/j.1567-1364.2005.00010.x
Jin, H., Reed, D. W., Thompson, V. S., Fujita, Y., Jiao, Y., Crain-Zamora, M., et al. (2019). Sustainable bioleaching of rare earth elements from industrial waste materials using agricultural wastes. ACS Sustain. Chem. Eng. 7, 15311–15319. doi: 10.1021/acssuschemeng.9b02584
Karpe, A. V., Beale, D. J., Godhani, N. B., Morrison, P. D., Harding, I. H., and Palombo, E. A. (2015). Untargeted metabolic profiling of winery-derived biomass waste degradation by Penicillium chrysogenum. J. Agric. Food Chem. 63, 10696–10704. doi: 10.1021/acs.jafc.5b04834
Karpe, A. V., Dunn, M. S., Taylor, M. C., Nguyen, T., Ong, C., Karla, T., et al. (2018). Nitrogen deprivation in Fusarium oxysporum promotes mycotoxin production via intermediates in the Krebs cycle and unreported methylmalonyl-CoA mutase activity. Metabolomics 14:160. doi: 10.1007/s11306-01801459-0
Kell, D. B., Brown, M., Davey, H. M., Dunn, W. B., Spasic, I., and Oliver, S. G. (2005). Metabolic footprinting and systems biology: the medium is the message. Nat. Rev. Microbiol. 3, 557–565. doi: 10.1038/nrmicro1177
Khorassani, R., Hettwer, U., Ratzinger, A., Steingrobe, B., Karlovsky, P., and Claassen, N. (2011). Citramalic acid and salicylic acid in sugar beet root exudates solubilize soil phosphorus. BMC Plant Biol. 11:121. doi: 10.1186/1471-2229-11-121
Liaud, N., Giniés, C., Navarro, D., Fabre, N., Crapart, S., Herpoël-Gimbert, I., et al. (2014). Exploring fungal biodiversity:organic acid production by 66 strains of filamentous fungi. Fungal Biol. Biotechnol. 1:1. doi: 10.1186/s40694-014-0001-z
Magnuson, J., and Lasure, L. (2004). “Organic Acid Production by Filamentous Fungi,” in Advances in Fungal Biotechnology for Industry, Agriculture, and Medicine, eds J. Tkacz, and L. Lange (Berlin: Springer), 307–340. doi: 10.1007/978-1-4419-8859-1_12
Nautiyal, C. S., Bhadauria, S., Kumar, P., Lal, H., Mondal, R., and Verma, D. (2000). Stress induced phosphate solubilization in bacteria isolated from alkaline soils. FEMS Microbiol. Lett. 182, 291–296. doi: 10.1016/s0378-1097(99)00605-9
Osman, Y., Gebreil, A., Mowafy, A. M., Anan, T. I., and Hamed, S. M. (2019). Characterization of Aspergillus niger siderophore that mediates bioleaching of rare earth elements from phosphorites. World J. Microbiol. Biotechnol. 35:93. doi: 10.1007/s11274-019-2666-1
Papagianni, M. (2007). Advances in citric acid fermentation by Aspergillus niger: biochemical aspects, membrane transport and modeling. Biotechnol. Adv. 25, 244–263. doi: 10.1016/j.biotechadv.2007.01.002
Peelman, S., Sun, Z. H. I., Sietsma, J., and Yang, Y. (2014). “Leaching of rare earth elements: past and present,” in Proceedings of the 1st European Rare Earth Resources Conference, eds E. Balomenos, D. Panias, and I. Paspaliaris, Milos, 446–456.
Richter, C. L., Dunn, B., Sherlock, G., and Pugh, T. (2013). Comparative metabolic footprinting of a large number of commercial wine yeast strains in Chardonnay fermentations. FEMS Yeast Res. 13, 394–410. doi: 10.1111/1567-1364.12046
Scervino, J. M., Papinutti, V. L., Godoy, M. S., Rodriguez, M. A., Della Monica, I., Recchi, M., et al. (2011). Medium pH, carbon and nitrogen concentrations modulate the phosphate solubilization efficiency of Penicillium purpurogenum through organic acid production. J. Appl. Microbiol. 110, 1215–1223. doi: 10.1111/j.1365-2672.2011.04972.x
Shi, Y., Pan, C., Wang, K., Chen, X., Wu, X., Chen, C.-T. A., et al. (2017). Synthetic multispecies microbial communities reveal shifts in secondary metabolism and facilitates cryptic natural product discovery. Environ. Microbiol. 19, 3606–3618. doi: 10.1111/1462-2920.13858
Sue, T., Obolonkin, V., Griffiths, H., and Villas-Bôas, S. G. (2011). An Exometabolomics approach to monitoring microbial contamination in microalgal fermentation processes by using metabolic footprint analysis. Appl. Environ. Microbiol. 77, 7605–7610. doi: 10.1128/aem.00469-11
Uehling, J. K., Entler, M. R., Neredith, H. R., Millet, L. J., Timm, C. M., Aufrecht, J. A., et al. (2019). Microfluidics and metabolomics reveal symbiotic bacterial-fungal interactions between Mortierella elongata and Burkholderia include metabolite exchange. Front. Microbiol. 10:2163. doi: 10.3389/fmicb.2019.02163
Vesth, T. C., Nybo, J. L., Theobald, S., Frisvad, J. C., Larsen, T. O., Nielsen, K. F., et al. (2018). Investigation of inter- and intraspecies variation through genome sequencing of Aspergillus section Nigri. Nat. Genet. 50, 1688–1695. doi: 10.1038/s41588-018-0246-1
Watkinson, S., Bebber, D., Darrah, P., Fricker, M., Tlalka, M., and Boddy, L. (2006). “The role of wood decay fungi in the carbon and nitrogen dynamicas of the forest floor,” in Fungi in Biogeochemical Cycles, ed. G. M. Gadd (Cambridge: Cambridge University Press), 151–181. doi: 10.1017/cbo9780511550522.008
Keywords: bioleaching, rare earth elements, thorium, metabolomics, Paecilomyces, phosphate solubilizing
Citation: Brisson VL, Zhuang W-Q and Alvarez-Cohen L (2020) Metabolomic Analysis Reveals Contributions of Citric and Citramalic Acids to Rare Earth Bioleaching by a Paecilomyces Fungus. Front. Microbiol. 10:3008. doi: 10.3389/fmicb.2019.03008
Received: 05 August 2019; Accepted: 13 December 2019;
Published: 14 January 2020.
Edited by:Aindrila Mukhopadhyay, Lawrence Berkeley National Laboratory, United States
Reviewed by:M. Oves, King Abdulaziz University, Saudi Arabia
Young-Mo Kim, Pacific Northwest National Laboratory (DOE), United States
Erica L.-W. Majumder, SUNY College of Environmental Science and Forestry, United States
Copyright © 2020 Brisson, Zhuang and Alvarez-Cohen. 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) and the copyright owner(s) 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.