Next Article in Journal
Batch Simultaneous Saccharification and Fermentation of Primary Sludge at Very High Solid Concentrations for Bioethanol Production
Next Article in Special Issue
Specific Organic Loading Rate Control for Improving Fermentative Hydrogen Production
Previous Article in Journal
Engineering the C3N Pathway as a Short Detour for De Novo NAD+ Biosynthesis in Saccharomyces cerevisiae
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dark Fermentation in the Dark Biosphere: The Case of Citrobacter sp. T1.2D-12

by
Violeta Gallego-Rodríguez
1,†,
Adrián Martínez-Bonilla
1,†,
Nuria Rodríguez
1,2 and
Ricardo Amils
1,2,*
1
Centro de Biología Molecular Severo Ochoa (CSIC-UAM), Universidad Autónoma de Madrid, 28049 Madrid, Spain
2
Centro de Astrobiología (CAB), CSIC-INTA, 28850 Torrejón de Ardoz, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fermentation 2023, 9(10), 887; https://doi.org/10.3390/fermentation9100887
Submission received: 22 August 2023 / Revised: 23 September 2023 / Accepted: 27 September 2023 / Published: 30 September 2023
(This article belongs to the Special Issue Fermentative Biohydrogen Production)

Abstract

:
Microbial diversity that thrives in the deep subsurface remains largely unknown. In this work, we present the characterization of Citrobacter sp. T1.2D-1, isolated from a 63.6 m-deep core sample extracted from the deep subsurface of the Iberian Pyrite Belt (IPB). A genomic analysis was performed to identify genes that could be ecologically significant in the IPB. We identified all the genes that encoded the formate–hydrogen lyase and hydrogenase-2 complexes, related to hydrogen production, as well as those involved in glycerol fermentation. This is particularly relevant as some of the substrates and byproducts of this process are of industrial interest. Additionally, we conducted a phylogenomic study, which led us to conclude that our isolate was classified within the Citrobacter telavivensis species. Experimentally, we verified the strain’s ability to produce hydrogen from glucose and glycerol and, thus, of performing dark fermentation. Moreover, we assessed the activity of the nitrate and tetrathionate reductase complexes and the isolate’s ability to tolerate high concentrations of heavy metals, especially Zn. These results suggest that C. telavivensis T1.2D-1 can play a role in the carbon, hydrogen, iron, nitrogen, and sulfur cycles that occur in the deep subsurface of the IPB, making it a candidate worthy of further study for possible biotechnological applications.

1. Introduction

The production of hydrogen via dark fermentation encompasses processes performed by microorganisms under aphotic and anaerobic conditions using organic substrates. The interest in this approach has increased due to the possibility of using waste biomass or wastewater as substrates, which would not only reduce hydrogen’s production costs, but can also offer a sustainable way to generate clean energy [1]. However, the importance of hydrogen is not limited to its potential as a renewable energy source; there are multiple environments where this gas is essential to microbial communities’ well-being, including hot springs, the gut microbiome, soils, groundwaters, or the deep subsurface, among others [2,3,4,5,6]. The deep subsurface, in particular, is characterized by the absence of light and oxygen, as well as higher temperatures and pressure as the depth increases [6]. In this environment, the mineralogical composition and hydrology greatly influence life since the lack of organic matter becomes a constraint, and minerals may be the only substrates available for microbial development [7]. These conditions collectively render the deep subsurface an extreme environment, often referred to as the dark biosphere, where microorganisms develop completely isolated from the surface [8].
Subsurface ecosystems that thrive in the absence of photosynthesis-derived products are known as Subsurface Lithoautotrophic Microbial Ecosystems (SLiMEs). At present, the most accepted model holds that the primary producers, which are part of the dark biosphere, are powered by hydrogen [9,10]. Moreover, the crucial role that this gas plays in underground communities was reported [5,9,11,12], and the existence of hydrogen produced by microorganisms within these habitats was also demonstrated [13]. This is the case, for instance, of the Iberian Pyrite Belt (IPB), one of the world’s largest sulfide deposits. Located in the southeast of the Iberian Peninsula, between Huelva (Spain) and Grándola (Portugal), it is 240 km long by 35 km wide and Rio Tinto rises from its core in Peña de Hierro (Huelva, Spain) [14].
Río Tinto is a river renowned for its high content of heavy metals, mainly ferric iron, acidic pH (with a mean of 2.3), and rich microbial diversity [15]. For many years, these extreme conditions were attributed to mining activity, which has been ongoing in the area for over 5000 years [11]. Nevertheless, there are indications that these characteristics predate human intervention. Recent studies have suggested that the IPB hosts an underground bioreactor. The discharge from this bioreactor contains microbial metabolic byproducts that are believed to be responsible for the extreme conditions observed in the river’s basin [16]. In particular, these characteristics have been associated with a reaction known as Nitrate-Dependent Fe2+ Oxidation (NDFO) [11]. It is still unclear whether enzymes are involved or if it is an abiotic process governed by nitrogen reactive forms [17]. In either case, it is believed that the ferric iron generated in this reaction participates in the dissolution of pyrite under anaerobic conditions, producing protons, sulphate, and ferrous iron, which is oxidized by the NDFO creating a cycle [18,19,20]. This can explain the high concentration of these compounds in Río Tinto, as well as its acidic pH. In order to investigate the geological processes occurring in this ecosystem, the Iberian Pyrite Belt Subsurface Life Detection (IPBSL) drilling project was established. Two boreholes, BH10 and BH11, were drilled and sampled, reaching depths of 620 and 340 m, respectively. As a result of this study, the vast microbial diversity that hosts the subsurface of the IPB was assessed, and several microorganisms were isolated under anaerobic conditions. Some of them have already been characterized through either in vitro or in silico approaches and their adequacy in participating in various biogeochemical cycles within this ecosystem demonstrated [11].
The strain labeledT1.2D-1 was isolated from 63.6 m-deep samples from borehole BH11 using culture media designed for the enrichment of chemolithoautotrophic denitrifying microorganisms and an atmosphere consisting of H2:CO2. The preliminary analysis indicated that the strain was a member of the Citrobacter genus [21], which consisted of rod-shaped, Gram-negative, facultative anaerobic bacteria classified within the family Enterobacteriaceae. Traditionally, these microorganisms are studied for their clinical relevance [22]. However, they have also been isolated from different environmental samples, including sediments [23,24], water [25], and wastewater [26]. More recently, this genus has been further explored for its potential applications in biotechnology and the industry [27,28].
In this study, our focus is to explore the relevant features that can explain the ecological role of Citrobacter sp. T1.2D-1 in the deep subsurface of the IPB, with a particular emphasis on its fermentation potential. To achieve this, we conduct a genomic analysis that allows us to identify genes of ecological or biotechnological significance. The metabolic activities of interest are subsequently verified experimentally, proving the ability of the strain T1.2D-1 to perform dark fermentation using various substrates and to participate in the different biogeochemical cycles occurring in the subsurface of the IPB.

2. Materials and Methods

2.1. Methodology for Citrobacter sp. T1.2D-1 Isolation

The methodology followed for drilling borehole BH11 is explained by Amils et al. (2013) [29]. The isolation of the strain T1.2D-1 was performed using ~6 g of powdered rock from a 63.6 m-deep core sample extracted from borehole BH11 as inoculum. The isolation procedure was performed in an anaerobic chamber using anaerobic culture media for the enrichment of chemolithoautotrophic denitrifying microorganisms composed of 0.3 g/L NH4Cl, 0.3 g/L K2HPO4·3H2O, 0.1 g/L MgSO4·7H2O, 0.4 g/L NaHCO3, 0.01 g/L CaCl2·2H2O, 0.1 g/L yeast extract, 40 mM nitrate, 20 mM thiosulfate, and 1 ml/L of trace element solution [30] at an initial pH of 7 and an atmosphere composed of H2:CO2 (80:20, v/v). The isolation procedure is described by Leandro et al. (2018) [31].

2.2. DNA Extraction and Sequencing

Genomic DNA was extracted from an aerobic culture of the strain T1.2D-1 grown in Tryptic Soy Broth (TSB) using the Cetyltrimethylammonium Bromide (CTAB)-based method [32]. The final DNA concentration was determined with a Qubit v.2.0 fluorometer (Thermo Fisher, Scientific ES, Madrid, Spain). Library preparation and sample sequencing were performed by MicrobesNG (Birmingham, UK) using the Illumina MiSeq platform.

2.3. Genomic Analysis

The quality of paired-end reads was assessed using the Comprehensive Genome Analysis tool of the BV-BRC database v3.29.20 [33]. Genome assembly and annotation were conducted using Galaxy v23.0.rc1 [34]. In particular, SPAdes v3.15.4 was used to perform de novo assembly in contigs of the sequences previously trimmed by MicrobesNG, and plasmidSPAdes v3.15.4 to search for and assemble putative plasmids [35,36]. For contigs generation, all the parameters were maintained in their default values except k-mers, using 49, 91, 101, 121, and 127. Contigs were merged into scaffolds using SSPACE v3.0 software for Linux [37]. The assembled genome was represented using the Circular Viewer tool of BV-BRC.
Genome annotation was conducted using the RAST platform v2.0 [38], keeping all the parameters in their default values. Gene classification into different subsystems was visualized using the Seed viewer tool. The search for proteins related to metabolisms of interest was conducted using the BLAST protein tool from the NCBI server. The sequences of these proteins were retrieved from the UniProt database. We considered a positive result whenever the query cover was higher than 50%, the E-value less than 10-5, and the percent identity higher than 35%.
A phylogenomic analysis was performed through the calculation of Digital DNA-DNA Hybridization (dDDH), Average Nucleotide Identity (ANI), and Average Amino acid Identity (AAI). The dDDH values were estimated using the Genome-to-Genome Distance Calculator (GGDC) from the Leibniz Institute DMSZ [39]. ANI and AAI values were determined with the NCBI Prok tool from the Microbial Genome Atlas (MiGA) [40]. The phylogenomic tree was built with the BV-BRC Bacterial Genome Tree tool, using 1000 genes, 100 bootstrap replicas, and keeping the rest of the parameters as default values. The genomes of the strains used were retrieved from the NCBI database. Archaeopteryx.js v2.0.0a6 was used to visualize the tree [41].

2.4. Utilization of Different Electron Donors under Aerobic and Anaerobic Conditions

The electron donors used by Citrobacter sp. T1.2D-1 under aerobic and anaerobic conditions were determined in M9 medium, composed of 1 g/L NH4Cl, 11 g/L Na2HPO4·7H2O, 3 g/L KH2PO4, 5 g/L NaCl, 4 g/L carbon source, 0.120 g/L MgSO4, and 10 mg/L CaCl2 [42]. The pH was adjusted to 7, before autoclaving with NaOH 1M. The carbon sources used under aerobic conditions were acetate, arabinose, citrate, fructose, galactose, glucose, glycerol, lactate, mannitol, pyruvate, starch, and sucrose; under anaerobic conditions, they were acetate, citrate, galactose, glucose, glycerol, lactose, malic acid, pyruvate, and sucrose. A total of 1 ml of a culture in exponential growth in medium PaFe2N2 [43] supplemented with 10 mM of acetate and 5 mM of nitrate for the anaerobic assays, or in medium TSB for the aerobic assays, were inoculated into 15 ml of medium M9. The experiments were conducted at room temperature in duplicates plus a non-inoculated negative control. The aerobic cultures were incubated in a shaker at 130 rpm whilst the anaerobic ones were left static.
Growth was determined by measuring the optical density at 600 nm using a Thermo Scientific Evolution 300 UV-Visible spectrophotometer. Durham tubes were used for the detection of gas production.

2.5. Evaluation of the Consumption and Production of Hydrogen

The consumption of hydrogen by Citrobacter sp. T1.2D.1 was tested in M9 medium using fumarate as the carbon source and an atmosphere composed of H2:N2 (80:20, v/v). Its production was evaluated in the same medium; however, glycerol and glucose were added as electron donors in an atmosphere of N2. The production experiments were inoculated with 1 ml of the strain in exponential growth in medium M9 supplemented with glucose or glycerol into a volume of 15 ml. The consumption experiment was inoculated with a culture grown in medium PaFe2N2 using the same volumes as the production assay. Both experiments were performed at room temperature and in triplicates. Additionally, the glucose concentration was increased to 10 g/L and glycerol to 40 g/L. The production of gases was detected with a Bruker Series Bypass 450GC chromatograph. It was equipped with a column CP2056 0.6 m x 1/8‘Ultimetal Cromsorb GHP 100–120 mesh and a column CP81073 0.5 m × 1/8’ Ultimetal Hayesep Q 80–100 mesh; a detector TCD at 200 °C for the detection of H2 and CO2 and a detector FID at 250 °C for the detection of CH4. N2 were used as carrier gases.

2.6. Nitrate-Dependent Fe2+ Oxidation Assay

The NDFO assay was performed using PaFe2N2 medium with 5 mM of nitrate, 10 mM of acetate, and 4 mM of FeCl2. The preparation of the medium, its inoculation, and measurements of iron oxidation were performed as explained by Mateos et al. (2022) [43]. The inoculated bottles were manipulated inside an anaerobic chamber in a nitrogen atmosphere (Coy Laboratory products, EEUU). Nitrate reduction was measured through a N-(1-Naphthyl)ethylenediamine (NEDA)-based method [44].

2.7. Tetrathionate Reduction Assay

To assess the ability of Citrobacter sp. T1.2D-1 to reduce tetrathionate, we used a medium composed of 13.6 g/L KH2PO4, 2 g/L (NH4)2SO4, 0.2 g/L MgSO4·7H2O, 2 g/L galactose, 0.3 g/L yeast extract, 0.3 g/L casamino acids, 1 g/L K2S4O6, and 5 g/L KOH [45]. All the components were sterilized by autoclaving, except for tetrathionate, which was sterilized using syringe filters with a pore size of 0.22 μm (Nuclepore, Whatman, Maidstone, UK). A total of 1 ml of a culture in exponential growth cultured in PaFe2N2 supplemented with 5 mM of nitrate and 10 mM of acetate under anaerobic conditions was used as inoculum placed into 15 ml of the medium. A total of 1 ml of each culture was extracted for the measurement of tetrathionate and centrifuged at 4000 x g for 5 minutes. The supernatant was preserved at –20 °C until measurement. The reduction of tetrathionate was quantified using a modified version of the method proposed by Kelly et al. (1969) [46]. Briefly, 200 μL of sample were mixed with 320 μL of NaH2PO4 0.2 M-NaOH 0.2 N buffer at pH 7.4, 280 μL of water, and 400 μL of 0.05 M KCN at 4 °C. The mixtures were incubated for 20 minutes at 4 °C and tempered for 15 minutes. Finally, 240 μL of 1.5 M Fe(NO3)3·9H2O 0.5 M in 4 N HClO4 and 560 μL of water were added. Spectrophotometric measurements were taken at 460 nm.

2.8. pH Tolerance Assay

The tolerance of Citrobacter sp. T1.2D-1 to different pH values, ranging from 3 to 10, was evaluated in duplicates using Luria Bertani (LB) broth as the basal medium at room temperature and a 130 rpm shaking speed. A total of 1 ml of an aerobic culture in exponential growth in medium TSB was used as inoculum placed into 15 ml of LB at each pH. pH adjustment was conducted after autoclaving using 1 M NaOH and HCl stocks sterilized by filtration and a Thermo Scientific Orion 920a pH meter. Growth was determined by measuring the optical density at 600 nm.

2.9. Heavy Metal-Resistance Assay

The heavy metal-resistance assay was performed in triplicate with a negative control in Tryptic Soy Agar (TSA) medium. After autoclaving and before gelation, the TSA medium was supplemented with different concentrations of each metal using 0.5 M stocks of the different salts sterilized by filtration. The salts used were CrCl3·6H2O, CuSO4·5H2O, CdSO4·8/3H2O, NiCl2·6H2O, CoSO4·7H2O, Pb(NO3)2, ZnCl2, and HgCl2 at concentrations of 0.1 mM, 0.5 mM, 1 mM, 5 mM, 7.5 mM, 10 mM, 12.5 mM, 15 mM, 17.5 mM, 20 mM, and 22.5 mM, respectively. The agar plates were incubated at room temperature for two weeks.

3. Results and Discussion

3.1. Genomic Analysis

3.1.1. Genome Assembly and Annotation

We sequenced the genome of Citrobacter sp. T1.2D-1 to characterize it and identify the metabolic activities of interest. The sequencing run generated 951,331 paired-end reads of 2 × 250 bp resulting in a mean coverage of 85X. The bacterial chromosome was assembled into 81 contigs and 67 scaffolds, with a length of 5,435,771 bp and a GC content of 53.56% (Figure 1). No putative plasmids were identified. These results align with the known characteristics of the Citrobacter genus, which typically exhibits a genome size and GC content ranging from 4 to 5 Mb and 51–56%, respectively [47].
Genome annotation using the RAST platform found 80 tRNA genes, 12 rRNA genes, and 5,312 Coding DNA Sequences (CDSs). Nearly 50% of the genes were classified into the categories of carbohydrate metabolism, protein metabolism, and amino acid metabolism. To further investigate the metabolic diversity of strain T1.2D-1 and its ecological role in the subsurface of the IPB, we focused our study on identifying genes encoding enzymes that enabled the isolate to thrive under the extreme oligotrophic conditions in this ecosystem [11] (Table S1).
Starting with the hydrogen cycle, we found the presence of genes that encoded both the formate–hydrogen lyase (FHL) and hydrogenase-2 (Hyd-2) complexes. This finding positions Citrobacter sp. T1.2D-1 as a significant asset in the context of the subsurface of the IPB, as it has the ability to produce and consume hydrogen (Figure 2). The FHL system consists of a formate dehydrogenase that transfers electrons to a [NiFe] hydrogenase encoded by the hyc genes. This complex has been extensively characterized and is known to be responsible for the hydrogen produced by some facultative anaerobes during mixed-acid fermentation [48]. On the contrary, the Hyd-2 complex can oxidize hydrogen in the periplasm coupled to the reduction of menaquinones, upon which the complex Fumarate Reductase (FRD) is dependent. Furthermore, the reversibility of this system has been reported in another strain of the genus [49], allowing hydrogen production during glycerol fermentation when the quinone pool is excessively reduced [50]. Therefore, Citrobacter sp. T1.2D-1 exhibits the potential to generate H2 through fermentation using two different mechanisms depending on the carbon source.
We also identified several activities that would allow the isolate to utilize different carbon sources under anaerobic conditions. The strain T1.2D-1 has the complete mixed-acid fermentation pathway, suggesting that it has the potential to ferment acetate, formate, oxaloacetate, fumarate, malate, or citrate. We also identified the frc gene, involved in the conversion of oxalate to formate, both substrates found in the IPB [11]. The production of CO2 via this pathway holds ecological significance since it can be used in the production of organic matter, acetate, or methane by microorganisms capable of performing carbon fixation, acetogenesis, or methanogenesis, respectively, contributing to the survival of fermenters or methanotrophs. Likewise, we delved into the ability of the isolate to ferment glycerol, as it is a process that is arousing increasing interest in the industry [51]. Citrobacter sp. T1.2D-1 has all the genes to convert glycerol into 1,3-propanediol (1,3-PDO) via a reductive pathway, which involves glycerol dehydratase and 1,3-PDO oxidoreductase activities, or into hydrogen through an oxidative route that implies the production of glycolytic byproducts and the mixed-acid fermentation route [52]. Therefore, strain T1.2D-1 can produce hydrogen when glycerol is used as a substrate via Hyd-2 or FHL complexes. Furthermore, 1,3-PDO is relevant since it can be used in the synthesis of biopolymers from petrochemical compounds and in different industries, such as cosmetics or food [53]. Moreover, it has been reported that its production by microorganisms is advantageous in comparison to chemical synthesis [54]. In addition, we also identified the genes involved in the conversion of dihydroxyacetone phosphate into methylglyoxal and the posterior synthesis of 1,2-PDO, which is also highly demanded in the industry due to its wide range of applications [55].
Concerning the nitrogen cycle, we found the complete operon of nitrate reductase, indicating that the isolate had the ability to reduce nitrate to nitrite (Figure 3). Nitrite is particularly significant in the subsurface of the IPB as it has been associated with NDFO, which, as mentioned above, is considered one of the main factors contributing to the extreme conditions observed in the Río Tinto basin. Moreover, nitrate has been found in the subsurface of the IPB, as well as nitrite [11], meaning that the reaction would be feasible. Furthermore, the ability of strain T1.2D-1 to perform this reaction involves it in the iron cycle of the IPB.
Lastly, regarding the sulfur cycle, which is of special relevance in this context since the IPB is one of the largest metal sulfide deposits on Earth, we found the complete tetrathionate reductase operon (Figure 3). These preliminary results prove that Citrobacter sp. T1.2D-1 can be metabolically active in the deep subsurface of the IPB, contributing to the sustenance of hydrogen, carbon, nitrogen, iron, and sulfur biogeochemical cycles identified in this ecosystem.

3.1.2. Phylogenomic Analysis

As mentioned above, a previous analysis using the 16S rRNA gene sequence indicated that the strain T1.2D-1 was classified within the Citrobacter genus. To further identify it taxonomically, we conducted a genomic-based taxonomic analysis through the calculation of phylogenomic indices, such as ANI, AAI, and dDDH. This procedure was performed against the NCBI database. The results show that our isolate belongs to the species C. telavivensis, owing to the fact that we obtained the highest values of ANI (99.92%), AAI (95%), and dDDH (99.8%) for the strain C. telavivensis 6105 compared to the rest of species of the genus (Table S2), and that the three parameters exceed the thresholds for species demarcation (ANI ≥ 95%; AAI ≥ 95–96%; dDDH ≥ 70%) [38,56,57]. To further confirm this result, we built a phylogenomic tree using the codon tree method, where C. telavivensis appeared again as the species immediately closest to the isolate (Figure S1).

3.2. Carbon Sources Utilization

Regarding the preliminary results for the carbon cycle, we studied the ability of Citrobacter sp. T1.2D-1 to use various carbon sources under aerobic and anaerobic conditions. The results show that the isolate is able to use the majority of the carbon sources tested for both conditions (Table 1), providing evidence that several reactions of the mixed-acid fermentation route are active, and that the strain T1.2D-1 can potentially use diverse carbon sources detected in the subsurface of the IPB, such as formate or oxalate.
These results confirm that Citrobacter sp. T1.2D-1, in accordance with the characteristics of its genus [22], can ferment a wide range of substrates, including low-value substrates, such as glycerol. At present, high quantities of glycerol are generated as a byproduct of other industrial processes. Consequently, this surplus needs to be disposed of or transformed into products with a higher value, which can be performed through fermentation [58]. The fermentation of glycerol, as previously mentioned, can produce 1,3-PDO and 1,2-PDO, which have a higher industrial value. Therefore, it would be interesting to study the ability of the isolate to produce these two compounds. Furthermore, we can also study if the isolate can also ferment other low-value substrates, such as wastewater, waste biomass, or agricultural biomass, which have already been described as fermentation substrates [59,60,61].
The fermentative ability of Citrobacter sp. T1.2D-1 is also interesting for its products. For instance, since it can perform mixed acid fermentation, Citrobacter sp. T1.2D-1 can produce compounds with industrial interest, such as lactate, used both in the food and non-food industries [62], or ethanol, which is has gained more interest in recent years as it can be used as a biofuel [63].

Hydrogen Production

In the previous section, we detected, through Durham tubes, the production of gases in certain fermentations. To identify the gas and quantify its production, we performed gas chromatography using anaerobic cultures, where glucose and glycerol were used as sole carbon sources. The Citrobacter genus is known to be one of the most studied and promising hydrogen producers due to its ability to produce this gas from a wide range of substrates and various environmental conditions [64,65]. We observed that strain T1.2D-1 produced hydrogen at room temperature at an initial pH of 7 using glucose and glycerol as electron donors (Figure 4).
Hydrogen is of great significance in the deep subsurface because the absence of light impedes the development of photosynthetic microorganisms capable of generating organic matter for heterotrophs to consume; thus, alternative energy sources are needed. Thus, it might be valuable to examine this ability further and evaluate if, in similar conditions as those found in the subsurface, Citrobacter sp. T1.2D-1 can also produce this gas, since in IPB, sugars, such as glucose, mannose, and sorbitol, have been detected (V. Parro, personal communication) and, as mentioned above, biohydrogen production has been detected in deep subsurface samples [11,13]. Furthermore, it would be interesting to perform an optimization of the conditions of hydrogen production and calculate the maximum yield of the isolate, in order to compare it with the ones described for other strains of the genus and to explore the possible applications.

3.3. Other Metabolisms of Interest

The genomic analysis showed that the strain T1.2D-1 could be metabolically active in the IPB. This is supported by the identification of multiple genes known to be involved in the biogeochemical cycles that occur in this system [11]. With the aim of testing the metabolic activities identified in the genome and how Citrobacter sp. T1.2D-1 can thrive in an extreme environment, such as the deep subsurface, several assays were designed and performed.

3.3.1. Nitrate-Dependent Fe2+ Oxidation Assay

We aimed to confirm the activity of the nitrate reductase operon and the ability of the strain to perform the NDFO. A medium supplemented with nitrate as the electron acceptor, acetate as the electron donor, and FeCl2 was used. We observed that, in a two-week period, Citrobacter sp. T1.2D-1 reduced nearly half of the nitrate, producing nitrite. Concurrently, we confirmed the oxidation of ferrous iron (Figure 5).
This is the second time that it has been reported that a strain of the Citrobacter genus is capable of performing the NDFO [66]. These results emphasize that the deep subsurface of the IPB harbors microorganisms capable of performing reactions associated with the extreme conditions detected in the Río Tinto basin and iron cycling, as explained in Section 3.1. To explore this matter more thoroughly, it would be valuable to investigate how Citrobacter sp. T1.2D-1 interacts with pyrite. Studying whether the isolate can participate in its dissolution [18,19,20], as well as its ability to interact with the mineral surface [67], would help to gain a further insight into the ecology of the deep subsurface of the IPB.

3.3.2. Tetrathionate Reduction Assay

As previously mentioned, in the genomic analysis, we identified the complete tetrathionate reductase operon. In order to prove the activity of this enzyme, we cultured Citrobacter sp. T1.2D-1 in a medium supplemented with galactose as the electron donor and tetrathionate as the electron acceptor under anaerobic conditions. We measured the concentration of tetrathionate in the medium and observed that the strain was capable of effectively reducing all of the tetrathionate present within a two-day period (Figure S3). This makes Citrobacter sp. T1.2D-1 a likely candidate for the sulfur cycle of the IPB. The production of thiosulfate through these mechanisms holds ecological relevance since this element can be utilized in the production of sulfates or sulfides. Notably, both elements have already been detected at various depths within the IPB and can serve as substrates for some of the microorganisms identified in this ecosystem [11].

3.3.3. Hydrogen Consumption Assay

With the aim of further investigating the strain’s metabolic diversity, we studied whether Citrobacter T1.2D-1 had the ability to use hydrogen in the presence of fumarate. The isolate was cultured in medium supplemented with fumarate and an H2:N2 atmosphere. We observed that, within two weeks, the strain used nearly half of the hydrogen present in the headspace of the serum bottles (Figure 6), proving the activity of the Hyd-2 complex and the ability of the isolate to use this gas as an electron donor.
Notably, recent articles have reported Citrobacter strains capable of performing carbon fixation under heterotrophic conditions [68,69]. This process is associated with the presence of a homolog of the rnfC gene of Acetobacterium woodii, which is also identified in the genome of Citrobacter sp. T1.2D-1, and with the ability to use hydrogen as an electron donor via Hyd-2 [69]. Therefore, it would be interesting to study whether strain T1.2D-1 is also capable of performing this reaction, since in an environment, such as the deep subsurface, this can be of great importance when carbon sources are scarce.

3.3.4. pH Tolerance Assay

Within the first 150 m of borehole BH11, the pH was determined as acid, ranging from 3.5 to 6.5 (R. Amils, personal communication). In particular, at a 63.6 m-depth from which Citrobacter sp. T1.2D-1 was isolated, the pH was nearly 4. We determined that the isolate was capable of growing at pH 4.0–9.0 in LB (Figure S2), and thus the strain T1.2D-1 was capable of tolerating the acidic conditions described for some of the depths of the BH11.

3.3.5. Heavy Metal Tolerance Assay

As stated previously, Río Tinto is known for the high concentration of heavy metals in its waters, elements that have also been reported in the subsurface of the IPB [11]. Additionally, members of the Citrobacter genus have been isolated from environments where these compounds are present and have been studied with bioremediation purposes [70]. For these reasons, we considered it an interesting proposition to test whether strain T1.2D-1 was capable of tolerating high concentrations of heavy metals. We determined that the isolate could grow in medium supplemented with 1 mM of Hg2+, 7.5 mM of Co2+ and Cu2+, 10 mM of Ni2+, 12.5 mM of Cd2+, Cr3+, and Pb2+, and 20 mM of Zn2+. The ability of Citrobacter sp. T1.2D-1 to resist high concentrations of heavy metals is clearly superior to those reported for other enterobacteria, such as Escherichia coli-type strains DH5α and K-12 [71,72], as well as other strains of the Citrobacter genus [70]. In particular, we found a remarkable tolerance to Ni2+ and Zn2+. Only one Citrobacter strain, belonging to the Citrobacter intermedius species, was described as being able to grow in higher concentrations of Ni2+ (with a minimal inhibitory concentration of 15 mM) [73], and none were reported with a tolerance to Zn2+ as high as the one determined for T1.2D-1. Zn has been referred to as a soil contaminant of concern. Therefore, taking into account the tolerance of Citrobacter sp. T1.2D-1 to Zn, it might be interesting to look at this ability in more detail, as this isolate can also be potentially used in the bioremediation of Zn [74].
Moreover, this noteworthy tolerance also has a genetic basis. We found multiple genes associated with resistance to these types of compounds, such as cusRS and cusCFBA, the former being a response system that induces the expression of the latter, which functions as a pump that retires the excess of Cu (I) from the periplasmic space of E. coli. We also identified the transcriptional regulator cueR, which enhances the expression of the also present copA, whose function is to transport Cu (I) from the cytosol to the periplasm, where it is reduced to Cu (II) (a less toxic form) by CueO [75]. We also identified nikABCDE operon, which encodes an Ni import system [76], and rcnA, which is translated into a membrane-associated protein whose presence has been related to an increase in resistance to Co and Ni. Lastly, focusing on the mechanisms of tolerance to Zn, we found zntA and zntR, which coded for a type-P ATPase that was predicted as the first Zn-specific transporter, responsible for the release of Zn ions out of the cell. In addition, ZntA could also be activated by Cd, Ag, and Pb, translocating this last element [77]. Conversely, we identified two Zn-uptake systems, another P-type ATPase and a high-affinity periplasmic-binding protein encoded by the genes zosA and znuABC [74,77], respectively, as well as zupT, which were not only involved in importing this metal into the cell but also of Cd, Fe, Co, and possibly Mn, given its broad specificity [78]. These results show that the strain T1.2D-1 is adapted to the presence of heavy metals, as expected from the environment in which it lives.

4. Conclusions

In this study, we reported the first Citrobacter strain isolated from hard-rock samples extracted from the deep subsurface. The genomic analysis indicates that strain T1.2D-1 belongs to the Citrobacter telavivensis species and helps identify the enzymatic activities of ecological significance. The experimental verification of features, such as hydrogen production and consumption, the fermentation of various sugars, nitrate reduction coupled to iron oxidation, and tetrathionate reduction, demonstrates that Citrobacter telavivensis T1.2D-1 can be metabolically active in the deep subsurface of the IPB, and that it may play a role in the hydrogen, carbon, nitrogen, iron, and sulfur biogeochemical cycles. Furthermore, the possibility of the isolate being metabolically active in this ecosystem is supported by its ability to tolerate high concentrations of heavy metals and a low pH, conditions characteristic of Río Tinto and the IPB. These results suggest that Citrobacter telavivensis T1.2D-1 is worthy of further study, as potential biotechnological applications for this strain can emerge. In particular, delving into some of its most remarkable features, such as its capacity to produce hydrogen and substrates of biotechnological interest, as well as its ability to withstand high concentrations of heavy metals, may very well reveal the potential uses of the strain in the industry and bioremediation, respectively.
The deep subsurface hosts, such enormous microbial diversity, remain, at present, uncharacterized. Since hydrogen production via dark fermentation becomes increasingly plausible for clean energy production, it is imperative to study different aspects with the aim of finding candidates that can take part in this process. In this respect, it is not unreasonable to think that some of them are members of the dark biosphere, giving us yet another reason to continue exploring the life beneath our feet.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/fermentation9100887/s1, Figure S1: phylogenomic tree built with the codon tree method using representative strains of the Citrobacter genus; Figure S2: growth of Citrobacter sp. T1.2D-1 at different pH values; Figure S3: tetrathionate reduction by Citrobacter sp. T1.2D-1; Table S1: genes related to metabolisms of interest identified in the genome of Citrobacter sp. T1.2D-1; Table S2: values of the phylogenomic analysis conducted comparing Citrobacter sp. T1.2D-1 with the closest isolates of the genus.

Author Contributions

Conceptualization, V.G.-R. and A.M.-B.; methodology, V.G.-R., A.M.-B. and N.R.; formal analysis, V.G.-R. and A.M.-B.; investigation, V.G.-R., A.M.-B. and N.R.; data curation, V.G.-R.; writing—original draft preparation, V.G.-R.; writing—review and editing, A.M.-B. and R.A.; supervision, R.A.; funding acquisition, R.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Ministry of Science and Innovation (projects TED2021-129563B-I00 and PID2022-136607NB-I00). A.M.B. acknowledges the Grant FPU19/01743 from the Spanish Ministry of Universities.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Reads have been deposited at DDBJ/ENA/GenBank under the accession number ERR11267709, and the complete genome sequences and annotations have been deposited under the accession numbers CATNHH010000000 for this chromosome and GCA_950098555.1 for the annotation.

Acknowledgments

The authors would like to thank D. Ruano for providing valuable insights during the development of this work, T. Leandro for the isolation of the strain, V. Parro for the unpublished results, the IPBSL team for obtaining the data from the drilled borehole, and R. Samalot and N. Aldaba for improving the English of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wong, Y.M.; Wu, T.Y.; Juan, J.C. A review of sustainable hydrogen production using seed sludge via dark fermentation. Renew. Sustain. Energy Rev. 2014, 34, 471–482. [Google Scholar] [CrossRef]
  2. Otaki, H.; Everroad, R.C.; Matsuura, K.; Haruta, S. Production and consumption of hydrogen in hot spring microbial mats dominated by a filamentous anoxygenic photosynthetic bacterium. Microb. Environ. 2012, 27, 293–299. [Google Scholar] [CrossRef]
  3. Hylemon, P.B.; Harris, S.C.; Ridlon, J.M. Metabolism of hydrogen gases and bile acids in the gut microbiome. FEBS Lett. 2018, 592, 2070–2082. [Google Scholar] [CrossRef] [PubMed]
  4. Fan, X.; Zhang, X.; Zhao, G.; Zhang, X.; Dong, L.; Chen, Y. Aerobic hydrogen-oxidizing bacteria in soil: From cells to ecosystems. Rev. Environ. Sci. Bio/Technol. 2022, 21, 877–904. [Google Scholar] [CrossRef]
  5. Ruff, S.E.; Humez, P.; de Angelis, I.H.; Diao, M.; Nightingale, M.; Cho, S.; Connors, L.; Kuloyo, O.O.; Seltzer, A.; Bowman, S. Hydrogen and dark oxygen drive microbial productivity in diverse groundwater ecosystems. Nat. Commun. 2023, 14, 3194. [Google Scholar] [CrossRef]
  6. Kieft, T.L. Microbiology of the Deep Continental Biosphere. In Advances in Environmental Microbiology; Springer: Berlin/Heidelberg, Germany, 2016; pp. 225–249. [Google Scholar] [CrossRef]
  7. Escudero, C.; Oggerin, M.; Amils, R. The deep continental subsurface: The dark biosphere. Int. Microbiol. 2018, 21, 3–14. [Google Scholar] [CrossRef] [PubMed]
  8. Meyer-Dombard, D.A.R.; Malas, J. Advances in defining ecosystem functions of the terrestrial subsurface biosphere. Front. Microbiol. 2022, 13, 891528. [Google Scholar] [CrossRef] [PubMed]
  9. Stevens, T.O.; McKinley, J.P. Lithoautotrophic microbial ecosystems in deep basalt aquifers. Science 1995, 270, 450–455. [Google Scholar] [CrossRef]
  10. Pedersen, K. Exploration of deep intraterrestrial microbial life: Current perspectives. FEMS Microbiol. Lett. 2000, 185, 9–16. [Google Scholar] [CrossRef]
  11. Amils, R.; Escudero, C.; Oggerin, M.; Puente Sánchez, F.; Arce Rodríguez, A.; Fernández Remolar, D.; Rodríguez, N.; García Villadangos, M.; Sanz, J.L.; Briones, C.; et al. Coupled C, H, N, S and Fe biogeochemical cycles operating in the continental deep subsurface of the Iberian Pyrite Belt. Environ. Microbiol. 2022, 25, 428–453. [Google Scholar] [CrossRef]
  12. Wu, X.; Pedersen, K.; Edlund, J.; Eriksson, L.; Åström, M.; Andersson, A.F.; Bertilsson, S.; Dopson, M. Potential for hydrogen-oxidizing chemolithoautotrophic and diazotrophic populations to initiate biofilm formation in oligotrophic, deep terrestrial subsurface waters. Microbiome 2017, 5, 37. [Google Scholar] [CrossRef] [PubMed]
  13. Sanz, J.L.; Rodriguez, N.; Escudero, C.; Carrizo, D.; Amils, R. Biological production of H2, CH4 and CO2 in the deep subsurface of the Iberian Pyrite Belt. Environ. Microbiol. 2021, 23, 3913–3922. [Google Scholar] [CrossRef] [PubMed]
  14. Tornos, F. Environment of formation and styles of volcanogenic massive sulfides: The Iberian Pyrite Belt. Ore Geol. Rev. 2006, 28, 259–307. [Google Scholar] [CrossRef]
  15. Amils, R.; González-Toril, E.; Fernández-Remolar, D.; Gómez, F.; Rodríguez, N.; Durán, C. Interaction of the sulfur and iron cycles in the Tinto River ecosystem. Rev. Environ. Sci. Biotechnol. 2002, 1, 299–309. [Google Scholar] [CrossRef]
  16. Amils, R.; Fernández-Remolar, D. IPBSL Team Río Tinto: A geochemical and mineralogical terrestrial analogue of Mars. Life 2014, 4, 511–534. [Google Scholar] [CrossRef] [PubMed]
  17. Bryce, C.; Blackwell, N.; Schmidt, C.; Otte, J.; Huang, Y.; Kleindienst, S.; Tomaszewski, E.; Schad, M.; Warter, V.; Peng, C. Microbial anaerobic Fe(II) oxidation–ecology, mechanisms and environmental implications. Environ. Microbiol. 2018, 20, 3462–3483. [Google Scholar] [CrossRef] [PubMed]
  18. Edwards, K.J.; Goebel, B.M.; Rodgers, T.M.; Schrenk, M.O.; Gihring, T.M.; Cardona, M.M.; Mcguire, M.M.; Hamers, R.J.; Pace, N.R.; Banfield, J.F. Geomicrobiology of pyrite (FeS2) dissolution: Case study at Iron Mountain, California. Geomicrobiol. J. 1999, 16, 155–179. [Google Scholar]
  19. Schippers, A.; Jørgensen, B.B. Biogeochemistry of pyrite and iron sulfide oxidation in marine sediments. Geochim. Cosmochim. Acta 2002, 66, 85–92. [Google Scholar] [CrossRef]
  20. Juncher Jørgensen, C.; Jacobsen, O.S.; Elberling, B.; Aamand, J. Microbial oxidation of pyrite coupled to nitrate reduction in anoxic groundwater sediment. Environ. Sci. Technol. 2009, 43, 4851–4857. [Google Scholar] [CrossRef]
  21. Leandro, T. Microbial Diversity and Anaerobic Metabolisms in the Subsurface of The Iberian Pyritic Belt. Ph.D. Dissertation; Universidad de Coimbra: Coimbra, Portugal, 2018. Available online: https://cab.inta-csic.es/en/tesis/microbial-diversity-and-anaerobic-metabolisms-in-the-subsurface-of-the-iberian-pyritic-belt/ (accessed on 1 January 2018).
  22. Borenshtein, D.; Schauer, D. The genus Citrobacter. Prokaryotes 2006, 6, 90–98. [Google Scholar]
  23. Furlan, J.P.R.; Ribeiro, T.G.; Novais, Â.; Cabral, H.; Peixe, L.; Stehling, E.G. Citrobacter meridianamericanus sp. nov., isolated from a soil sample. Int. J. Syst. Evol. Microbiol. 2023, 73, 005720. [Google Scholar] [CrossRef]
  24. Wang, H.; Hou, H.; Huang, J. Citrobacter arsenatis sp. nov., an arsenate-reducing bacterium isolated from freshwater sediment. Antonie Van Leeuwenhoek 2021, 114, 1285–1292. [Google Scholar] [CrossRef] [PubMed]
  25. Ribeiro, T.G.; Goncalves, B.R.; da Silva, M.S.; Novais, A.; Machado, E.; Carrico, J.A.; Peixe, L. Citrobacter portucalensis sp. nov., isolated from an aquatic sample. Int. J. Syst. Evol. Microbiol. 2017, 67, 3513–3517. [Google Scholar] [CrossRef] [PubMed]
  26. Jiang, X.; Cui, X.; Liu, W.; Xu, H.; Zheng, B. Genetic characterization of a novel sequence type of multidrug-resistant Citrobacter freundii strain recovered from wastewater treatment plant. Infect. Drug Res. 2019, 12, 2775–2779. [Google Scholar] [CrossRef] [PubMed]
  27. Metsoviti, M.; Zeng, A.; Koutinas, A.A.; Papanikolaou, S. Enhanced 1, 3-propanediol production by a newly isolated Citrobacter freundii strain cultivated on biodiesel-derived waste glycerol through sterile and non-sterile bioprocesses. J. Biotechnol. 2013, 163, 408–418. [Google Scholar] [CrossRef]
  28. Unuofin, J.O.; Moubasher, H.A.; Okoh, A.I.; Nwodo, U.U. Production of polyextremotolerant laccase by Achromobacter xylosoxidans HWN16 and Citrobacter freundii LLJ16. Biotechnol. Rep. 2019, 22, e00337. [Google Scholar] [CrossRef]
  29. Amils, R.; Fernández-Remolar, D.; Parro, V.; Rodríguez-Manfredi, J.A.; Timmis, K.; Oggerin, M.; Sánchez-Román, M.; López, F.J.; Fernández, J.P.; Puente, F.; et al. Iberian Pyrite Belt Subsurface Life (IPBSL), a drilling project of biohydrometallurgical interest. Adv. Mater. Res. 2013, 825, 15–18. [Google Scholar]
  30. Sanz, J.L.; Rodríguez, N.; Amils, R. Effect of chlorinated aliphatic hydrocarbons on the acetoclastic methanogenic activity of granular sludge. Appl. Microbiol. Biotechnol. 1997, 47, 324–328. [Google Scholar] [CrossRef]
  31. Leandro, T.; Rodriguez, N.; Rojas, P.; Sanz, J.L.; Da Costa, M.S.; Amils, R. Study of methanogenic enrichment cultures of rock cores from the deep subsurface of the Iberian Pyritic Belt. Heliyon 2018, 4, e00605. [Google Scholar] [CrossRef]
  32. Kate, W. Preparation of genomic DNA from bacteria. Curr. Protoc. Mol. Biol. 2001, 56, 2–4. [Google Scholar] [CrossRef]
  33. Olson, R.D.; Assaf, R.; Brettin, T.; Conrad, N.; Cucinell, C.; Davis, J.J.; Dempsey, D.M.; Dickerman, A.; Dietrich, E.M.; Kenyon, R.W. Introducing the Bacterial and Viral Bioinformatics Resource Center (BV-BRC): A resource combining PATRIC, IRD and ViPR. Nucleic Acids Res. 2023, 51, 678–689. [Google Scholar] [CrossRef]
  34. The Galaxy Community. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2022 update. Nucleic Acids Res. 2022, 50, 345–351. [Google Scholar] [CrossRef] [PubMed]
  35. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [PubMed]
  36. Antipov, D.; Hartwick, N.; Shen, M.; Raiko, M.; Lapidus, A.; Pevzner, P.A. plasmidSPAdes: Assembling plasmids from whole genome sequencing data. Bioinformatics 2016, 32, 3380–3387. [Google Scholar] [CrossRef] [PubMed]
  37. Boetzer, M.; Henkel, C.V.; Jansen, H.J.; Butler, D.; Pirovano, W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics 2011, 27, 578–579. [Google Scholar] [CrossRef]
  38. Aziz, R.K.; Bartels, D.; Best, A.A.; DeJongh, M.; Disz, T.; Edwards, R.A.; Formsma, K.; Gerdes, S.; Glass, E.M.; Kubal, M. The RAST Server: Rapid annotations using subsystems technology. BMC Genom. 2008, 9, 75. [Google Scholar] [CrossRef]
  39. Meier-Kolthoff, J.P.; Carbasse, J.S.; Peinado-Olarte, R.L.; Göker, M. TYGS and LPSN: A database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes. Nucleic Acids Res. 2022, 50, 801–807. [Google Scholar] [CrossRef]
  40. Rodriguez-R, L.M.; Gunturu, S.; Harvey, W.T.; Rosselló-Mora, R.; Tiedje, J.M.; Cole, J.R.; Konstantinidis, K.T. The Microbial Genomes Atlas (MiGA) webserver: Taxonomic and gene diversity analysis of Archaea and Bacteria at the whole genome level. Nucleic Acids Res. 2018, 46, 282–288. [Google Scholar] [CrossRef]
  41. Han, M.V.; Zmasek, C.M. phyloXML: XML for evolutionary biology and comparative genomics. BMC Bioinform. 2009, 10, 356. [Google Scholar] [CrossRef]
  42. Zhou, G.; Li, L.; Shi, Q.; Ouyang, Y.; Chen, Y.; Hu, W. Effects of nutritional and environmental conditions on planktonic growth and biofilm formation of Citrobacter werkmanii BF-6. J. Microbiol. Biotechnol. 2013, 23, 1673–1682. [Google Scholar] [CrossRef]
  43. Mateos, G.; Bonilla, A.M.; de Francisco de Polanco, S.; Martínez, J.M.; Escudero, C.; Rodríguez, N.; Sánchez-Andrea, I.; Amils, R. Shewanella sp. T2. 3D-1.1 a Novel Microorganism Sustaining the Iron Cycle in the Deep Subsurface of the Iberian Pyrite Belt. Microorganisms 2022, 10, 1585. [Google Scholar] [CrossRef] [PubMed]
  44. Bricio, C. Producción de óxidos de Nitrógeno Gaseosos en Thermus Thermophilus. Ph.D. Dissertation; Universidad Autónoma de Madrid: Madrid, Spain, 2012. Available online: http://hdl.handle.net/10486/11954 (accessed on 26 October 2012).
  45. Novotný, Č Instability of the tetrathionate respiratory chain of Citrobacter freundii. Folia Microbiol. 1978, 23, 428–432. [CrossRef] [PubMed]
  46. Kelly, D.P.; Chambers, L.A.; Trudinger, P.A. Cyanolysis and spectrophotometric estimation of trithionate in mixture with thiosulfate and tetrathionate. Anal. Chem. 1969, 41, 898–901. [Google Scholar] [CrossRef]
  47. Yuan, C.; Yin, Z.; Wang, J.; Qian, C.; Wei, Y.; Zhang, S.; Jiang, L.; Liu, B. Comparative genomic analysis of Citrobacter and key genes essential for the pathogenicity of Citrobacter koseri. Front. Microbiol. 2019, 10, 2774. [Google Scholar] [CrossRef] [PubMed]
  48. Sinha, P.; Roy, S.; Das, D. Role of formate hydrogen lyase complex in hydrogen production in facultative anaerobes. Int. J. Hydrogen Energy 2015, 40, 8806–8815. [Google Scholar] [CrossRef]
  49. Kim, S.; Seol, E.; Raj, S.M.; Park, S.; Oh, Y.; Ryu, D.D. Various hydrogenases and formate-dependent hydrogen production in Citrobacter amalonaticus Y19. Int. J. Hydrogen Energy 2008, 33, 1509–1515. [Google Scholar] [CrossRef]
  50. Evans, R.M.; Beaton, S.E. A novel overproduction system for the structural determination of a proton-pumping hydrogen-producing [NiFe]-hydrogenase. In Methods in Enzymology; Academic Press: Cambridge, MA, USA, 2018; Volume 613, pp. 91–116. [Google Scholar]
  51. Poblete-Castro, I.; Wittmann, C.; Nikel, P.I. Biochemistry, genetics and biotechnology of glycerol utilization in Pseudomonas species. Microb. Biotechnol. 2020, 13, 32–53. [Google Scholar] [CrossRef]
  52. Clomburg, J.M.; Gonzalez, R. Anaerobic fermentation of glycerol: A platform for renewable fuels and chemicals. Trends Biotechnol. 2013, 31, 20–28. [Google Scholar] [CrossRef]
  53. Ju, J.; Wang, D.; Heo, S.; Kim, M.; Seo, J.; Kim, Y.; Kim, D.; Kang, S.; Kim, C.; Oh, B. Enhancement of 1, 3-propanediol production from industrial by-product by Lactobacillus reuteri CH53. Microb. Cell Fact. 2020, 19, 6. [Google Scholar] [CrossRef]
  54. Lee, J.H.; Jung, M.; Oh, M. High-yield production of 1, 3-propanediol from glycerol by metabolically engineered Klebsiella pneumoniae. Biotechnol. Biofuels 2018, 11, 104. [Google Scholar] [CrossRef]
  55. Tao, Y.; Bu, C.; Zou, L.; Hu, Y.; Zheng, Z.; Ouyang, J. A comprehensive review on microbial production of 1, 2-propanediol: Micro-organisms, metabolic pathways, and metabolic engineering. Biotechnol. Biofuels 2021, 14, 216. [Google Scholar] [CrossRef] [PubMed]
  56. Jain, C.; Rodriguez-R, L.M.; Phillippy, A.M.; Konstantinidis, K.T.; Aluru, S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 2018, 9, 5114. [Google Scholar] [CrossRef] [PubMed]
  57. Konstantinidis, K.T.; Tiedje, J.M. Towards a genome-based taxonomy for prokaryotes. J. Bacteriol. 2005, 187, 6258–6264. [Google Scholar] [CrossRef] [PubMed]
  58. Yazdani, S.S.; Gonzalez, R. Anaerobic fermentation of glycerol: A path to economic viability for the biofuels industry. Curr. Opin. Biotechnol. 2007, 18, 213–219. [Google Scholar] [CrossRef] [PubMed]
  59. Hong, C.; Wu, J.; Huang, R.; Zhang, W.; He, W.; Deng, Z.; Han, Y.; Xiao, B.; Luo, H.; Qu, W. Effects of Temperature and Total Solid Content on Biohydrogen Production from Dark Fermentation of Rice Straw: Performance and Microbial Community Characteristics. Chemosphere 2022, 286, 131655. [Google Scholar] [CrossRef]
  60. Corsino, S.F.; Di Trapani, D.; Capodici, M.; Torregrossa, M.; Viviani, G. Optimization of acetate production from citrus wastewater fermentation. Water Resour. Ind. 2021, 25, 100140. [Google Scholar] [CrossRef]
  61. Liu, Z.; Zhang, C.; Lu, Y.; Wu, X.; Wang, L.; Wang, L.; Wang, L.; Han, B.; Xing, X.H. States and challenges for high-value biohythane production from waste biomass by dark fermentation technology. Bioresour. Technol. 2013, 135, 292–303. [Google Scholar] [CrossRef]
  62. Abedi, E.; Hashemi, S.M.B. Lactic Acid Production—Producing Microorganisms and Substrates Sources-State of Art. Heliyon 2020, 6, 04974. [Google Scholar] [CrossRef]
  63. Tse, T.J.; Wiens, D.J.; Reaney, M.J. Production of bioethanol—A review of factors affecting ethanol yield. Fermentation 2021, 7, 268. [Google Scholar] [CrossRef]
  64. Patel, S.K.; Kumar, P.; Mehariya, S.; Purohit, H.J.; Lee, J.; Kalia, V.C. Enhancement in hydrogen production by co-cultures of Bacillus and Enterobacter. Int. J. Hydrogen Energy 2014, 39, 14663–14668. [Google Scholar] [CrossRef]
  65. Cabrol, L.; Marone, A.; Tapia-Venegas, E.; Steyer, J.; Ruiz-Filippi, G.; Trably, E. Microbial ecology of fermentative hydrogen producing bioprocesses: Useful insights for driving the ecosystem function. FEMS Microbiol. Rev. 2017, 41, 158–181. [Google Scholar] [CrossRef] [PubMed]
  66. Li, B.; Pan, X.; Zhang, D.; Lee, D.; Al-Misned, F.A.; Mortuza, M.G. Anaerobic nitrate reduction with oxidation of Fe(II) by Citrobacter freundii strain PXL1—A potential candidate for simultaneous removal of As and nitrate from groundwater. Ecol. Eng. 2015, 77, 196–201. [Google Scholar] [CrossRef]
  67. Escudero, C.; Del Campo, A.; Ares, J.R.; Sánchez, C.; Martínez, J.M.; Gómez, F.; Amils, R. Visualizing microorganism-mineral interaction in the Iberian Pyrite Belt subsurface: The Acidovorax case. Front. Microbiol. 2020, 11, 572104. [Google Scholar] [CrossRef] [PubMed]
  68. Amulya, K.; Mohan, S.V. Fixation of CO2, electron donor and redox microenvironment regulate succinic acid production in Citrobacter amalonaticus. Sci. Total Environ. 2019, 695, 133838. [Google Scholar] [CrossRef] [PubMed]
  69. Jin, T.; Wang, Y.; Yao, S.; Hu, C.; Ma, T.; Xia, W. Bioconversion of carbon dioxide to succinate by Citrobacter. Chem. Eng. J. 2023, 452, 139668. [Google Scholar] [CrossRef]
  70. Ajmal, A.W.; Saroosh, S.; Mulk, S.; Hassan, M.N.; Yasmin, H.; Jabeen, Z.; Nosheen, A.; Shah, S.M.U.; Naz, R.; Hasnain, Z. Bacteria isolated from wastewater irrigated agricultural soils adapt to heavy metal toxicity while maintaining their plant growth promoting traits. Sustainability 2021, 13, 7792. [Google Scholar] [CrossRef]
  71. Hassan, S.; Abskharon, R.; Gad El-Rab, S.; Shoreit, A. Isolation, characterization of heavy metal resistant strain of Pseudomonas aeruginosa isolated from polluted sites in Assiut city, Egypt. J. Basic Microbiol. 2008, 48, 168–176. [Google Scholar] [CrossRef]
  72. Hassen, A.; Saidi, N.; Cherif, M.; Boudabous, A. Resistance of environmental bacteria to heavy metals. Bioresour. Technol. 1998, 64, 7–15. [Google Scholar] [CrossRef]
  73. Keramati, P.; Hoodaji, M.; Tahmourespour, A. Multi-metal resistance study of bacteria highly resistant to mercury isolated from dental clinic effluent. Afr. J. Microbiol. Res. 2011, 5, 831–837. [Google Scholar]
  74. Hussain, S.; Khan, M.; Sheikh, T.M.M.; Mumtaz, M.Z.; Chohan, T.A.; Shamim, S.; Liu, Y. Zinc essentiality, toxicity, and its bacterial bioremediation: A comprehensive insight. Front. Microbiol. 2022, 13, 900740. [Google Scholar] [CrossRef]
  75. Boyd, S.M.; Rhinehardt, K.L.; Ewunkem, A.J.; Harrison, S.H.; Thomas, M.D.; Graves, J.L., Jr. Experimental Evolution of Copper Resistance in Escherichia coli Produces Evolutionary Trade-Offs in the Antibiotics Chloramphenicol, Bacitracin, and Sulfonamide. Antibiotics 2022, 11, 711. [Google Scholar] [CrossRef] [PubMed]
  76. Diep, P.; Leo Shen, H.; Wiesner, J.A.; Mykytczuk, N.; Papangelakis, V.; Yakunin, A.F.; Mahadevan, R. Engineered nickel bioaccumulation in Escherichia coli by NikABCDE transporter and metallothionein overexpression. Eng. Life Sci. 2023, 23, 2200133. [Google Scholar] [CrossRef] [PubMed]
  77. Choudhury, R.; Srivastava, S. Zinc resistance mechanisms in bacteria. Curr. Sci. 2001, 81, 768. Available online: https://www.jstor.org/stable/24106396 (accessed on 21 July 2023).
  78. Roberts, C.S.; Ni, F.; Mitra, B. The zinc and iron binuclear transport center of ZupT, a ZIP transporter from Escherichia coli. Biochemistry 2021, 60, 3738–3752. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Circular representation of the genome of Citrobacter sp. T1.2D-1. From the outside to the inside, the following are shown: (1) the size of each contig in Mbp; (2) contigs in which the genome is assembled; (3) coding regions in the forward strand; (4) coding regions in the reverse strand; (5) non-coding regions; (6) GC content; (7) GC skew.
Figure 1. Circular representation of the genome of Citrobacter sp. T1.2D-1. From the outside to the inside, the following are shown: (1) the size of each contig in Mbp; (2) contigs in which the genome is assembled; (3) coding regions in the forward strand; (4) coding regions in the reverse strand; (5) non-coding regions; (6) GC content; (7) GC skew.
Fermentation 09 00887 g001
Figure 2. Complexes identified in the genomic analysis, which are involved in the production and consumption of hydrogen in Citrobacter sp. T1.2D-1. The consumption of hydrogen is represented as an interplay between hydrogenase-2 and fumarate reductase complexes. Hydrogen production is represented as a result of the activity of the formate–hydrogen lyase complex and of glycerol-3-phosphate dehydrogenase. Subunits of the enzymes are indicated. Electron flux is indicated with dashed arrows. Abbreviations: DHAP, dihydroxyacetone phosphate; FHL, formate–hydrogen lyase; FRD, fumarate reductase; G-3-P, glycerol-3-phosphate; G3P-DH, glycerol-3-phosphate dehydrogenase; Hyd-2, hydrogenase-2.
Figure 2. Complexes identified in the genomic analysis, which are involved in the production and consumption of hydrogen in Citrobacter sp. T1.2D-1. The consumption of hydrogen is represented as an interplay between hydrogenase-2 and fumarate reductase complexes. Hydrogen production is represented as a result of the activity of the formate–hydrogen lyase complex and of glycerol-3-phosphate dehydrogenase. Subunits of the enzymes are indicated. Electron flux is indicated with dashed arrows. Abbreviations: DHAP, dihydroxyacetone phosphate; FHL, formate–hydrogen lyase; FRD, fumarate reductase; G-3-P, glycerol-3-phosphate; G3P-DH, glycerol-3-phosphate dehydrogenase; Hyd-2, hydrogenase-2.
Fermentation 09 00887 g002
Figure 3. Overview of the biogeochemical cycles in which Citrobacter sp.T1.2D-1 can participate in the subsurface of the IPB based on the results of the genomic analysis. The name of the enzymatic activities involved is indicated in every reaction. Created with Biorender.
Figure 3. Overview of the biogeochemical cycles in which Citrobacter sp.T1.2D-1 can participate in the subsurface of the IPB based on the results of the genomic analysis. The name of the enzymatic activities involved is indicated in every reaction. Created with Biorender.
Fermentation 09 00887 g003
Figure 4. Hydrogen production from glucose (A) and glycerol (B) by Citrobacter sp. T1.2D-1. The variation in hydrogen concentration in the inoculated replicas (blue) is represented as a function of the negative control (orange). The standard deviation is indicated with error bars.
Figure 4. Hydrogen production from glucose (A) and glycerol (B) by Citrobacter sp. T1.2D-1. The variation in hydrogen concentration in the inoculated replicas (blue) is represented as a function of the negative control (orange). The standard deviation is indicated with error bars.
Fermentation 09 00887 g004
Figure 5. Results of the NDFO assay. (A) Graphic representation of the production of nitrite in the inoculated replicas (blue) versus the negative control (orange); (B) graphic representation of the Fe2+ oxidized throughout the experiment in the inoculated replicas. In both graphs, the error bars represent the standard deviation.
Figure 5. Results of the NDFO assay. (A) Graphic representation of the production of nitrite in the inoculated replicas (blue) versus the negative control (orange); (B) graphic representation of the Fe2+ oxidized throughout the experiment in the inoculated replicas. In both graphs, the error bars represent the standard deviation.
Fermentation 09 00887 g005
Figure 6. Hydrogen consumption by Citrobacter sp. T1.2D-1. The variation in hydrogen concentration in the inoculated replicas is represented as a function of the negative control. The typical deviation for each measure is indicated with error bars.
Figure 6. Hydrogen consumption by Citrobacter sp. T1.2D-1. The variation in hydrogen concentration in the inoculated replicas is represented as a function of the negative control. The typical deviation for each measure is indicated with error bars.
Fermentation 09 00887 g006
Table 1. Use of different electron donors by Citrobacter sp. T1.2D-1 under aerobic or anaerobic conditions. Carbon sources in which we detected the proliferation of the strain are marked with “+” and those in which there is no growth with “−“. The absence of data is indicated as “ND” and the electron donors in which we detected gas as fermentation product are marked with (G).
Table 1. Use of different electron donors by Citrobacter sp. T1.2D-1 under aerobic or anaerobic conditions. Carbon sources in which we detected the proliferation of the strain are marked with “+” and those in which there is no growth with “−“. The absence of data is indicated as “ND” and the electron donors in which we detected gas as fermentation product are marked with (G).
Carbon SourceGrowth
Aerobic ConditionsAnaerobic Conditions
Acetate++
Arabinose+ND
Citrate+
Fructose+ND
FumarateND+
Galactose++ (G)
Glucose++ (G)
Glycerol++ (G)
Lactate+ND
LactoseND+
Malic acidND+
Mannitol+ND
Pyruvate++
Starch+ND
Sucrose++ (G)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gallego-Rodríguez, V.; Martínez-Bonilla, A.; Rodríguez, N.; Amils, R. Dark Fermentation in the Dark Biosphere: The Case of Citrobacter sp. T1.2D-12. Fermentation 2023, 9, 887. https://doi.org/10.3390/fermentation9100887

AMA Style

Gallego-Rodríguez V, Martínez-Bonilla A, Rodríguez N, Amils R. Dark Fermentation in the Dark Biosphere: The Case of Citrobacter sp. T1.2D-12. Fermentation. 2023; 9(10):887. https://doi.org/10.3390/fermentation9100887

Chicago/Turabian Style

Gallego-Rodríguez, Violeta, Adrián Martínez-Bonilla, Nuria Rodríguez, and Ricardo Amils. 2023. "Dark Fermentation in the Dark Biosphere: The Case of Citrobacter sp. T1.2D-12" Fermentation 9, no. 10: 887. https://doi.org/10.3390/fermentation9100887

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop