Bioreactor performance
In this study, two types of electron donors, S0 and Fe0 were used for microbially-mediated chromate reduction. Results revealed that the Cr(VI) concentration in both bioreactors decreased significantly and remained stable for each cycle, herein the steady-state biogeochemical process was achieved in these reactors after 160 days incubation (Figure S1). In comparison with the S-bioreactor, the Fe-bioreactor showed a higher reduction rate of Cr(VI) (Fig. 1a). Two peaks corresponding to Cr 2p3/2 (576.0-578.0 eV) and Cr 2p1/2 (586.0-588.0 eV) in XPS spectrum indicated that the precipitates in both bioreactors were Cr(III) species, typically Cr(OH)3 (Figure S2a). In the S-bioreactor, a decline of Cr(VI) accompanied an increase of SO42+ (Fig. 1a), and no sulfite and thiosulfate were detected in the solution, suggesting that S0 oxidation to sulfate may be directly coupled with Cr(VI) reduction. In the Fe-bioreactor, aqueous Fe(II) and Fe(III) were not detected, likely due to the formation of Fe-bearing minerals at circum-neutral pH condition (Figure S2b) [10, 14], suggesting that Cr(VI) reduction may be coupled with Fe0 oxidation. The XPS analysis showed these precipitates had two typical peaks assigned to Fe 2p1/2 (725.0 − 727 eV) and Fe 2p3/2 (710.0 − 712.5 eV) (Figure S2b), which are similar to the peaks of Fe3O4 oxides. Interestingly, a time-course increase of methane was measured in the Fe-bioreactor (Fig. 1b), hinting methanogenesis or related methane metabolism in this bioreactor. Besides, CO2 was detected in both reactors and likely less abundant in the Fe-bioreactor than that in the S-bioreactor within one cycle, and VFAs (e.g., acetate, propionate, isobutyrate, butyrate, isovalerate, and valerate) were produced in both bioreactors (Fig. 1c), revealing that microbial-mediated C metabolism might play vital roles in both systems.
Abiotic experiments without addition of biomass to the bioreactors demonstrated that Cr(VI) reduction was fully attributed to the biotic process in the S-bioreactor, whereas about 14.8 ± 0.7% Cr(VI) was abiotically removed in the first cycle of Fe-bioreactor. After multiple cycles in the Fe-bioreactor, the steady-state removal efficiency remained around 1.9 ± 0.3% of Cr(VI). Therefore, a high reduction efficiency of Cr(VI) in the Fe-bioreactor (98.1 ± 1.2%) was mainly contributed by the biotic process, indicating Fe0 passivation inhibit the additional abiotic process. Initial reduction of Cr(VI) in the Fe-bioreactor might be driven by the electrons directly released from Fe0 oxidation, producing hydrogen [3, 10]. Altough Fe0 passivation prevented more abiotic processes, Fe0 corrosion could still induce a slow release of hydrogen in Fe-bioreactor [10]. Interestingly, hydrogen concentration was under detection limit in both bioreactors in the final cycle, illustrating that hydrogen may be extensively consumed once released. In brief, distinct geochemical variations in two ‘closed’ ecosystems point to different biogeochemical processes that were occuring.
Microbial community overview
To illuminate the biotic process in both bioreactors, metagenomics was constructed for samples from the original inoculum and bioreactors at the end of Cr(VI) reduction. A total of approximately 104 GB raw data were quality filtered and then co-assembled, yielding 89,774 scaffolds (≥ 500 bp) (Supplementary Table S1). Subquently, these 9,323 scaffolds (≥ 2000 bp) were binned, generating 23 draft genomes with estimated completeness of 52%-100% and < 1.5% contamination (Supplementary Table S2). These genomes ranged in size from 1.3 to 3.5 Mb with a GC content between 36.25% and 70.86%. Besides, 1,361 to 3,536 putative genes were predicted for individual genomes, with 80 to 97%, 58 to 75% and 75 to 93% of the gene annotation against the NCBI-nr, KEGG, and eggNOG databases, respectively. In addition, we detected 16S rRNA gene sequences in 15 of these 23 retrieved genomes.
The taxonomic classifications of these bins were inferred using a phylogenomic tree based on the concatenated alignment of 16 ribosomal proteins (Fig. 2a). Specifically, these 23 recovered genomes were affiliated with the phyla Euryarchaeota (n = 7), Bacteroidetes (n = 4), Spirochaetes (n = 2), Actinobacteria (n = 1) and Candidatus Aminicenantes (n = 1), and the classes β-Proteobacteria (n = 4), γ-Proteobacteria (n = 3) and δ-Proteobacteria (n = 1). The relative abundances of these genomes accounted for 59%, 73% and 74% of the communities in the original inoculum, and S- and Fe-bioreactors, respectively (Fig. 2b). In the original inoculum, the Betaproteobacteria (31.7%), Gammaproteobacteria (10.3%), and Bacteroidetes (6.4%) were dominant. After 160 days of incubation in S- and Fe-bioreactors, microbial communities were predominated by the Euryarchaeota (62.9% and 49.4% in the S- and Fe-bioreactors, respectively), followed by the Bacteroidetes (7.6%) in the S-bioreactor, and the Gammaproteobacteria (15.0%) and Spirochaetota (6.8%) in the Fe-bioreactor (Fig. 2b). Moreover, dominant species changed from Rhodocyclaceae QR_Bin.30 (26.1%) and Xanthomonas sp. QR_Bin.24 (8.2%) in the original inoculum to Methanobacterium sp. QR_Bin.6 (48.8%) and Bacteroidetes QR_Bin.18 (5.8%) in the S-bioreactor, and Methanobacterium sp. QR_Bin.6 (35.7%), Gammaproteobacteria QR_Bin.12 (15.0%), and Spirochaetes QR_Bin.15 (6.8%) in the Fe-bioreactor (Fig. 3). These findings indicated a strong selective pressure on microorganisms caused by different substrates, and also implied different mechanisms of microbial-mediated Cr(VI) reduction.
Metabolic potential of key players in bioreactors
Our experimental design capitalized on the niche differentiation required for carbon, nitrogen, sulfur and iron metabolisms coupling with Cr(VI) reduction in the two distinct ‘closed’ systems. To decipher the potential roles of these microorganisms, their metabolic pathways were reconstructed (Fig. 3).
Sulfur metabolism
The key genes dsrAB encoding dissimilatory sulfite reductase, which is considered as the most critical enzyme involved in dissimilatory sulfite reduction and oxidation [26–28], were identified in five recovered genomes (Table S3). The DsrAB phylogenetic tree demonstrated that three betaproteobacterial species (Thiobacillus sp. QR_Bin.5 and QR_Bin.21, and Rhodocyclaceae QR_Bin.30) possessed the oxidative dsrAB genes (Fig. 4a), and the reductive genes were found in Methanosaeta sp. QR_Bin.3 and Spirochaetes QR_Bin.4. Evidenced by the occurrence of the oxidative dsrAB, dsrEFH, aprAB, and sat genes, the three betaproteobacterial species likely carried out sulfur oxidation to sulfate through the reversed dissimilatory sulfate reduction pathway. Besides, they also encoded for sulfide dehydrogenase (fccAB), which are essential for sulfide oxidation to sulfur under anoxic conditions [29]. These above results illustrated that the three species enabled complete oxidation from sulfide to sulfate. Furthermore, Spirochaetes QR_Bin.4 might oxidize sulfur to sulfite as a part of sulfur disproportionation due to the fact that the reductive dsrAB and dsrD genes were present in the genome without dsrEFH genes [26], and Rhodocyclaceae QR_Bin.19 had a potential to mediate the transformation between sulfate and sulfite as it contained sat and aprAB genes, but they disappeared in both bioreactors. In addition, four betaproteobacterial species (QR_Bin.5, QR_Bin.19, QR_Bin.21 and QR_Bin.30) were capable of converting thiosulfate to sulfate due to the identified SOX system and poly-S (sulfiden−1) transformation (sqr and hyd), which accounted for 31.7% and 2.5% in the original inoculum and S-bioreactor, respectively. The above findings indicated an underlying synergetic and collaborative relationship for these microorganisms to complete the sulfur cycle.
Iron metabolism
Iron is not only a micronutrient for nearly all life on Earth, but also serves as an electron donor or acceptor by iron-oxidizing/reducing microorganisms [30, 31]. This study revealed that Gammaproteobacteria QR_Bin.7, Candidatatus Aminicenantes QR_Bin.23 and Bacteroidetes QR_Bin.26 were able to encode the cluster 3-affiliated Cyc2 homologs (Supplementary Table S2), which have been biochemically characterized to catalyze iron oxidation in Acidithiobacillus ferrooxidans [32] and are also found in neutrophilic, obligate iron-oxidizers [33], suggesting that they were potential iron-oxidizers. Moreover, Thiobacillus sp. QR_Bin.5 was inferred to carry out the potential for iron reduction/oxidation owing to the detection of mtoA and mtrB genes and the absence of the mtrC gene [34]. Previous studies have demonstrated that MtoA is homologous to the iron-reducing enzyme, MtrA, of Shewanella oneidensis MR-1, and in fact, it has been evidenced to rescue ΔmtrA mutants of MR-1, partially recovering the ability to reduce ferric iron [35]. Despite this, the MtoA homologs are frequently encoded by known and suspected iron-oxidizing bacteria [34], and the MtrC homolog, an outer-membrane cytochrome thought to participate in dissimilatory iron reduction [36], was missing in this species. Thus, we can not rule out the possibility that this species may use this enzyme to oxidize iron as previously proposed [35].
Methane metabolism
An interesting phenomenon attracted attention that methanogenesis was observed in the Fe-bioreactor with no methane detected in the S-bioreactor (Fig. 1b), despite the community compositions in the two bioreactors were highly similar (Fig. 2b). The key mcrABG genes encoding the methyl-coenzyme M reductase (MCR) complex were identified in all archaeal bins (Supplementary Table S3), suggesting that these archaea enriched in bioreactors were the important actors in methane metabolism. To unravel evolutionary relationships and characteristics of these MCR complexes detected, we built a phylogenetic tree using the concatenated alignment of McrABG subunits (Fig. 4b). Results revealed that these archaea could code for canonical MCR complexes, as they formed clusters with previously studied methanogens such as Methanomicrobiales, Methanobacterilales, and Methanomassiliicoccales [37]. Even so, we cannot rule out the possibility that they may use this enzyme for anaerobic oxidation of methane (AOM) as previously reported [38, 39]. The above-mentioned difference might be due to that in the Fe-bioreactor, H2 abiotically produced by Fe0 facilitated hydrogen-dependent methanogenesis [3, 10], resulting in a high CH4 concentration in the headspace of this reactor. However, for the S-bioreactor, a trace amount of methane generated by methanogens might be rapidly oxidized by methanotrophs, when coupled to sulfate or metal reduction by themselves or other bacteria [37, 40–42].
According to the substrate use, methanogens are broadly characterized as hydrogenotrophic (H2 and CO2), aceticlastic (acetate), methylotrophic (X-CH3) and H2-dependent methylotrophic (H2 and X-CH3) [37, 39]. Thus, Methanolinea sp. QR_Bin.1 and QR_Bin.22, and Methanobacterium sp. QR_Bin.6, QR_Bin.11 and QR_Bin.13 were inferred to be putative hydrogenotrophic methanogens, resulting from the identification of m-WL, mtrABCDEFGH, mcrABG, and frhA/ehbN genes (encoding the [Ni-Fe] hydrogenase) as previously reported [38, 39]. Remarkably, Methanobacterium sp. QR_Bin.6 and QR_Bin.11, and Methanolinea sp. QR_Bin.22 also coded for the carbon monoxide dehydrogenase-acetyl-CoA synthase (Codh-Acs) complex. This complex could activate acetate to generate reduced ferredoxin and methyl-H4MPT, a marker enzyme for acetoclastic methanogenesis, illustrating that these bins had a potential to perform both hydrogenotrophic and acetoclastic methanogenesis. The result was inconsistent with previous findings of the close lineages Methanobacterium formicicium, Methanolinea sp. SDB and Methanolinea tarda sp. NOBI-1 that are considered to be hydrogenotrophic methanogens [37]. Besides, Methanosaeta sp. QR_Bin.3 was also inferred to carry out acetoclastic methanogenesis, as members of the genus Methanosaeta are specialized on acetate degradation [37]. Moreover, Methanomassiliicoccus sp. QR_Bin.16 had the potential for H2-dependent methylotrophic methanogenesis, evidenced by the fact that this species could code for methyltransferases to support methyl-dependent methanogenesis from methanol (mtaABC), methanethiol (mtsAB) and methylamines (mtmB, mtbB and mtbC, and mttB and mttC), and previous studies showed that Methanomassiliicoccus populations were H2-dependent methylotrophic methanogens [37, 43].
Apart from the anaerobic methanogens, we found that Gammaproteobacteria QR_Bin.12 was a putative methane-oxidizing bacteria (MOB) due to the presence of the key genes pomABC encoding methane monooxygenase. This bin accounted for 15.0% of relative abundance in the Fe-bioreactor but only 0.65% in the S-related one. Previous studies revealed methane monooxygenases in a small number of bacterial species affiliated with the Gamma- and Alpha-proteobacteria and Verrucomicrobia habitating the aerobic-anaerobic interface [44–46], and that these microorganisms could use a variety of terminal electron acceptors (e.g., oxygen, nitrate) for methane oxidation [47–49]. As such, the methane produced by methanogens might be partially oxidized by Gammaproteobacteria QR_Bin.12, particularly in the Fe-bioreactor.
Chromium reduction
The primary purpose of this study was to investigate how Cr(VI) reduction was coupled with other biogeochemical cycles, thus genes involved in Cr(VI) reduction were confirmed in these retrieved genomes (Supplementary Table S3). Methanolinea sp. QR_Bin.1 harbored the potential to encode NfsA, an oxygen-insensitive nitroreductase and also a flavoprotein that can reduce Cr(VI) to less soluble and toxic Cr(III) [8, 50], implying that it was a putative chromate reducer. Notably, this reductive process would produce intracellular reactive oxygen species (ROS) that combine with DNA-protein complexes to cause cell DNA and protein damage [50, 51]. Besides that, the nemA gene was identified in Spirochaetes QR_Bin4, Thiobacillus sp. QR_Bin.5 and QR_Bin.21, Gammaproteobacteria QR_Bin.7 and QR_Bin.12, and Rhodocyclaceae QR_Bin.30. The FMN-dependent NADH-azoreductase gene azoR was found in Rhodocyclaceae Bin.19 and Thiobacillus sp. QR_Bin.21. The identification of these genes, which was confirmed to be involved in chromate reduction [8], suggests that all these microorganisms were potential chromate reducers. Previous research have proven that the nethylmaleimide reductase (NemA) belonging to the old yellow enzyme family of flavoproteins could use one or two electrons from the cofactors to reduce chromate intracellularly [52].
Under many oxygen-depleted conditions, Cr(VI) can be reduced extracellularly by coupling with intracellular oxidation of electron donors such as carbohydrates, proteins, fats, and hydrogen [2, 7, 8]. The cytochrome c3 was identified in most of these recovered genomes in both bioreactors (Supplementary Table S3), suggesting that they were able to reduce Cr(VI) through the electron transfer chain. The cytochrome families are frequently shown to participate in extracellular Cr(VI) reduction under anaerobic conditions, including cytochrome c3 of Desulfovibrio vulgaris [53] and Desulfomicrobium norvegicum [54, 55]. Furthermore, the reduction of Cr(VI) may also take place by chemical reactions associated with organic compounds present in intra/extra cellular, such as amino acids, nucleotides, sugars, vitamins, organic acids and glutathione [2, 7, 8]. In addition to the enzymatic reduction, chromate efflux from cells is an efficient and widespread mechanism to decrease chromium toxicity [56, 57]. The protein ChrA, a hydrophobic membrane protein that can export chromate from cytoplasm or periplasm driven by the proton motive force [8, 58], was inferred to be encoded in Conexibacter sp. QR_Bin.2, Thiobacillus sp. QR_Bin.5 and QR_Bin.21, Gammaproteobacteria QR_Bin.7, and Spirochaetes QR_Bin.15. In brief, mechanisms of microbial-mediated Cr(VI) reduction are extremely complicated with multiple possible pathways and unstable redox intermediates [4, 6–9].
Carbon fixation
After one day incubation, relatively stable CO2 concentrations were observed in both bioreactors (Fig. 1b), thus carbon fixation might exist. Among these potential autotrophs (including four archaea and three bacteria), only Methanobacterium sp. QR_Bin.11 might use the Wood-Ljungdahl (WL) pathway to fix CO2, and the others likely perform the Calvin-Benson-Bassham (CBB) cycle (Fig. 3 and Table S3). Intriguingly, these putative archaeal autotrophs and one potential bacterial autotroph (Methanolinea sp. QR_Bin.1 and QR_Bin.22, Methanosaeta sp. QR_Bin.3, Methanobacterium sp. QR_Bin.11, Methanomassiliicoccus sp. QR_Bin.16 and Thiobacillus sp. QR_Bin.21) were enriched in both bioreactors relative to the original inoculum (Fig. 3), revealing that Bacteria primarily contributed to carbon fixation in the original inoculum, whereas Archaea became the main carbon fixers after a long-term incubation in both bioreactors.
Nitrogen metabolism
Key genes relevant to dissimilarity nitrate reduction were evaluated (Supplementary Table S3), and the results revealed that ten species were likely involved in the process transforming nitrate to ammonia, of which Thiobacillus sp. QR_Bin.21 and Gammaproteobacteria QR_Bin.12 were the main participants in the S- and Fe-bioreactors, respectively. The pentaheme nitrite reductase (nrfAH), essential for respiratory nitrite ammonification, was identified in Spirochaetes QR_Bin.4, suggesting that it might be responsible for nitrite reduction to ammonia. Interestingly, the nifDKH genes encoding nitrogenase for nitrogen fixation were confirmed in Euryarchaeota, including Methanosaeta sp. QR_Bin.3, Methanobacterium sp. QR_Bin.6 and QR_Bin.11, and Methanolinea sp. QR_Bin.22, illustrating that these archaea were potentially responsible for supplying organic nitrogen for cell growth in the bioreactors.
Fermentation
The detection of VFAs in both bioreactors gave evidence for microbial fermentation [59] (Fig. 1c). The present study showed that most of the recovered species were able to produce acetate via the Pta-Ack pathway (the pta-ackA genes) and acylphosphatase (acyP) and/or perform ethanol fermentation due to the identification of adh gene in both bioreactors rather than in the original inoculum. In contrast, only Bacteroidetes QR_Bin.27 harbored the potential to produce lactate due to the LDH gene identified.
VFA oxidation
The generated VFAs, e.g., acetate, butyrate, propionate (Fig. 1c), could be oxidized by heterotrophic microorganisms via the reversed WL pathway [60], beta-oxidation pathway [61, 62], or methylmalonyl-CoA pathway [63], linking to the reduction of sulfate, ferric iron and chromium to generate energy for cells [64–66]. As only Methanobacterium sp. QR_Bin.11 harbored key genes encoding the WL pathway [67], this species might oxidize acetate via the reversed WL pathway when CO2 was depleted. The complete beta-oxidation pathway were found in Spirochaetes QR_Bin.4 and gammaproteobacterium QR_Bin.12 (Supplementary Table S3), unveiling their underlying potentials to catalyze butyrate [62]. The complete conventional propionate oxidation pathway, methylmalonyl-CoA pathway [63], were detected in Conexibacter sp. QR_Bin.2, Candidatatus Aminicenantes QR_Bin.23, and Bacteroidetes QR_Bin.18 and QR_Bin, of which Bacteroidetes QR_Bin.18 and QR_Bin.26 potentially mediated the conversion of propionate to succinate and then to acetate in bioreactors. Accrodingly, these VFA-oxidizing microorganisms may play an imperative role in the tranfromation of Cr(VI) to Cr(III) in bioreactors, analogous to previous findings [10, 14]. In addition, the oxidation of butyrate and propionate produces the methanogenic precursors such as acetate [62, 63], which can further facilitate aceticlastic methanogenesis [37].
Integrating inferred niches and activities in bioreactors
In the two bioreactors, despite different substrates (S0 and Fe0) were used, a highly similar community composition appeared to be involved in efficient Cr(VI) removal (Fig. 2). Distinct biogeochemical processes suggest that not only methanogens but also methanotrophs were involved. Moreover, diverse syntrophic metabolisms were found between Cr(VI) reduction and the oxidation of methane, sulfur, iron and acetate in the two ‘closed’ bioreactors (Fig. 5), indicating a metabolic mutualism among different functional members.
In the S-bioreactor, Cr(VI) reduction may be coupled with three pathways, AOM by anaerobic methane-oxidizing archaea (ANME), VFAs oxidation by VFA-oxidizers, and sulfur oxidation by sulfur-oxidizers. First, methane was likely produced via hydrogenotrophic, acetoclastic and H2-dependent methylotrophic methanogenesis. Nonetheless, as methane flux was below the limit of detection in this bioreactor, the generated methane may be oxidized immediately by ANME and methane-oxidizing bacteria (MOB) that possibly cooperated with sulfate- and metal-reducing bacteria [37, 40–42]. The mcr genes do not merely encode the key enzyme, MCR complexes for methanogenesis, but are also present in ANME to reversely oxidize methane, and the [Ni-Fe] hydrogenases (e.g., VhoA, MvhA, FrhA, and EhbN) are an additional indicator for methane-oxidizing potential [37–39]. In this study, two enriched archaea (Methanosaeta sp. QR_Bin.3 and Methanomassiliicoccus sp. QR_Bin.16) did not code for the [Ni-Fe] hydrogenases, but harbored key genes for dissimilarity sulfate reduction (DsrAB) or Cr(VI) reduction (cytochrome c3), thereby they per se possibly hold the potential to couple AOM with the reduction of SO42− or Cr(VI). Next, under elevated CO2 in the S-bioreactor, the putative autotrophic Methanosaeta sp. QR_Bin.11, Methanolinea sp. QR_Bin.22 and Thiobacillus sp. QR_Bin.21 could fix inorganic carbon to support chemolithoautotrophic growth under anaerobic conditions. The resulting dissolved organic carbon (DOC) would be used to ferment into acetate, lactate, methanol or other VFAs by Bacteroidetes QR_Bin.18 and QR_Bin.27 that might be further utilized to produce methane or oxidized by heterotropic Bacteroidetes QR_Bin.26 when coupled with Cr(VI) reduction. Lastly, S0 could be directly used as an effective reductant for Cr(VI) reduction by sulfur-oxdizer Thiobacillus sp. QR_Bin.21, which coupled sulfur oxidation with Cr(VI) reduction.
In the Fe-bioreactor, the reduction of Cr(VI) could be coupled with methane oxidation by MOB and organic matter oxidation by acetate-oxidizing heterotrophs. Prior studies have demonstrated that Fe0 is readily utilized as a slow-release electron donor for methanogenesis [68]. The electron transfer mechanism from Fe0 to microorganisms is through H2, which is generated by the chemical reaction of Fe0 with H2O under anaerobic conditions. In consideration of the lower CO2 and higher acetate (Figs. 1b and 1c), the generated H2 tended to explain the obviously higher CH4 production in Fe-bioreactor than that in S-bioreactor (Fig. 1b), despite both reactors harbored the similar archaea-dominating communities where the hydrogen-dependent methanogenesis dominated. A putative MOB gammaproteobacterium QR_Bin.12, which harbored both pomABC and nemA, constituted the highest abundance in the Fe-bioreactor, primarily owing to the sufficient methane as electron donor. Inconsistent to ANME-mediated AOM in the S-bioreactor, this MOB would perform AOM coupled with Cr(VI) reduction in this bioreactor, which was possibly one of the major contributors for Cr(VI) reduction and might partly explain the higher reduction rate in the Fe-bioreactor (Fig. 1b). Similar to the S-bioreactor, VFA-oxidizing heterotroph Bacteroidetes QR_Bin.26 and gammaproteobacterium QR_Bin.12 could also catalyze Cr(VI) reduction in the Fe-bioreactor, and the produced acetate via VFA oxidation can support methanogenesis. Although the bioreactor was running under a steady-state condition and low abiotic Fe(II) oxidation (< 2%) was detected, we could not rule out that the generated Fe(II) may abiotically reduce Cr(VI) due to its strong reducing potential [10, 14].