Evaluation of acidogenesis products’ effect on biogas production performed with metagenomics and isotopic approaches
Biotechnology for Biofuels volume 14, Article number: 125 (2021)
During the acetogenic step of anaerobic digestion, the products of acidogenesis are oxidized to substrates for methanogenesis: hydrogen, carbon dioxide and acetate. Acetogenesis and methanogenesis are highly interconnected processes due to the syntrophic associations between acetogenic bacteria and hydrogenotrophic methanogens, allowing the whole process to become thermodynamically favorable. The aim of this study is to determine the influence of the dominant acidic products on the metabolic pathways of methane formation and to find a core microbiome and substrate-specific species in a mixed biogas-producing system.
Four methane-producing microbial communities were fed with artificial media having one dominant component, respectively, lactate, butyrate, propionate and acetate, for 896 days in 3.5-L Up-flow Anaerobic Sludge Blanket (UASB) bioreactors. All the microbial communities showed moderately different methane production and utilization of the substrates. Analyses of stable carbon isotope composition of the fermentation gas and the substrates showed differences in average values of δ13C(CH4) and δ13C(CO2) revealing that acetate and lactate strongly favored the acetotrophic pathway, while butyrate and propionate favored the hydrogenotrophic pathway of methane formation. Genome-centric metagenomic analysis recovered 234 Metagenome Assembled Genomes (MAGs), including 31 archaeal and 203 bacterial species, mostly unknown and uncultivable. MAGs accounted for 54%–67% of the entire microbial community (depending on the bioreactor) and evidenced that the microbiome is extremely complex in terms of the number of species. The core microbiome was composed of Methanothrix soehngenii (the most abundant), Methanoculleus sp., unknown Bacteroidales and Spirochaetaceae. Relative abundance analysis of all the samples revealed microbes having substrate preferences. Substrate-specific species were mostly unknown and not predominant in the microbial communities.
In this experimental system, the dominant fermentation products subjected to methanogenesis moderately modified the final effect of bioreactor performance. At the molecular level, a different contribution of acetotrophic and hydrogenotrophic pathways for methane production, a very high level of new species recovered, and a moderate variability in microbial composition depending on substrate availability were evidenced. Propionate was not a factor ceasing methane production. All these findings are relevant because lactate, acetate, propionate and butyrate are the universal products of acidogenesis, regardless of feedstock.
Anaerobic digestion (AD) of biomass to methane and carbon dioxide is a widespread process in the natural environment or it can be the result of human activity. This complex process involves the interaction of many groups of microorganisms responsible for, respectively, hydrolysis of polymeric compounds to monomers, acidic fermentations, acetogenesis and methanogenesis in anaerobic environments where the concentrations of electron acceptors, such as nitrates, compounds of iron (III) and manganese (IV) or sulfates, are low and the redox potential is below –240 mV [1,2,3,4]. AD is both a promising method for renewable energy production and a way to treat tons of waste generated around the world. AD can be conducted in one-stage or multi-stage systems where hydrolysis and acidogenesis are separated from acetogenesis and methanogenesis. Multistep systems provide optimal conditions for each step, stabilize the processes, and increase energy recovery from biomass [5, 6].
Due to the limited number of substrates for methanogenesis (acetate, carbon dioxide and hydrogen/formate, methylated compounds), the acetogenic and methanogenic steps of AD are tightly connected, the methanogens and acetogens forming syntrophic systems. The essence of the syntrophic interactions is interspecies electron transfer (IET) making the entire syntrophic metabolism efficient and thermodynamically favorable. IET can occur indirectly, mediated by hydrogen and formate (IIET), or directly (DIET) via contact between the microorganisms [7,8,9,10]. The metabolic pathways utilized for syntrophic oxidation of common non-gaseous products of acidogenesis include beta-oxidation for butyrate; the methylmalonyl-CoA pathway recognized in Syntrophobacter [11, 12] or the dismutation pathway recognized in Smithella propionica for propionate ; the Wood–Ljungdahl pathway for acetate; the pathway of ethanol oxidation recognized in the genera Pelobacter and Desulfovibrio in the absence of other electron acceptors ; the lactate oxidation recognized in Desulfovibrio in the absence of sulfate [9, 15]. Studies on Acetobacterium woodii revealed methanogens-independent metabolic pathways of (i) ethanol and carbon dioxide conversion to acetate involving a bifunctional ethanol/acetaldehyde dehydrogenase ; (ii) lactate oxidation to acetate involving a complex composed of the FAD-dependent lactate dehydrogenase (LDH) and the electron transfer flavoprotein (EtfA/B) . Both require the reductive carbon monoxide dehydrogenase/acetyl-CoA synthase pathway. Genes encoding the lactate-oxidizing metabolic machinery homologous to those of A. woodii and D. vulgaris are present in the domain Bacteria. The Syntrophomonadaceae are highly specialized syntrophic microbes that can oxidize butyric, propionic and long-chain fatty acids, with the best recognized species being the butyrate-oxidizing Syntrophomonas wolfei and the propionate-oxidizing Syntrophobacter fumaroxidans, S. wolinii, S. pfennigii, S. sulfatireducens[19,20,21,22]. Another butyrate oxidizer is Syntrophus acidotrophicus (Syntrophobacterales order in Deltaproteobacteria), whereas representatives of the Desulfotomaculum and Pelotomaculum genera from the Peptococcaceae family (Clostridiales)  or Smithella propionica (Syntrophobacterales order in Deltaproteobacteria)  are recognized propionate oxidizers. Known acetate-oxidizing bacteria belong to the following groups: Synergistetes—genera Synergistes ; Clostridia—Thermoacetogenium phaeum, Clostridium ultunense, Clostridium sporomusa and Moorella sp.; and the Deltaproteobacteria—Geobacter spp. ; Spirochaetes . Uncultivable Cloacimonetes, including WWE1 (Waste Water of Evry 1), are capable of propionate and butyrate oxidation [9, 27]. Representatives of Chloroflexi and Plantomycetes are probably also involved in butyrate oxidation . Previously, methane-producing Archaea were thought to belong only to the Euryarchaeota (orders: Methanobacteriales, Methanococcales, Methanomicrobiales, Methanosarcinales, Methanopyrales, Methanocellales and Methanomassiliicoccales) [28,29,30,31], but recent studies indicate that they may also belong to Bathyarchaeota and Verstaraeteachaeota [32,33,34,35,36]. It is noteworthy that many Archaea and most of the Bacteria belong to the group of "viable but non-culturable" (VBNC) microorganisms, and new species are detected on the basis of data obtained from metagenomic studies of environmental or anaerobic digesters’ samples . Methane is a by-product in the reaction of methyl coenzyme M with coenzyme B, catalyzed by archaeal methyl coenzyme M reductase (encoded by the mtrA gene) regardless of the methanogenic pathways, i.e., splitting of acetate (acetoclastic/acetotrophic methanogenesis); reduction of CO2 with H2 or formate and, rarely, ethanol or secondary alcohols as electron donors (hydrogenotrophic methanogenesis); reduction of methyl groups of methylated compounds such as methanol, methylated amines or methylated sulfides (hydrogen-dependent and hydrogen-independent methylotrophic methanogenesis) [38, 39]. Surprisingly, only two known genera Methanosarcina and Methanotrix, formerly Methanosaeta, members of the order Methanosarcinales, are capable of methane production from acetate. However, the recent studies change the commonly accepted viewpoint that Methanothrix species can only utilize acetate as substrate via the acetoclastic pathway. It was shown that the ethanol-fed methanogenic community was dominated by the Methanotrix species which were metabolically active via the carbon dioxide reductive (hydrogenotrophic) pathway, rather than the acetate decarboxylation (acetotrophic) pathway. Methanothrix species are capable of accepting electrons via the direct interspecies electron transfer (DIET) for the reduction of carbon dioxide to methane [40, 41]. The identified methylotrophic methanogens belong to the Methanosarcinales order. All the other known methanogens produce methane by reduction of CO2 [2, 4, 28, 29, 39].
Even though the general scheme of anaerobic digestion is commonly known, it is still not completely characterized at the molecular level. Our understanding of the microbial ecology and physiology associated with AD is crucial for process stabilization and optimization. Since the majority of microorganisms involved in AD are not cultivable, the culture-dependent techniques are insufficient. Thus amplicon-based approaches (e.g., 16S rRNA gene sequencing), metagenomics, and genome-centric techniques were employed to analyze samples from different-scale bioreactors [37, 42,43,44]. The former concentrate on sequencing of the genes encoding 16S rRNA or other selected genes, e.g., the mcrA gene, while metagenomics and genome-centric techniques involve shotgun sequencing of total DNA. All these approaches confirmed that the majority of the species involved were not isolated but can be classified at the genus, family or order level performing a taxonomic assignment of the genes/protein encoded [42, 45,46,47,48,49,50]. Furthermore, the genome-centric approaches allowed obtaining large numbers of metagenome-assembled genomes (MAGs) from anaerobic digesters. In their majority they represent novel microorganisms residing in AD microbiomes whose physiology and potential functions are being recognized. In fact, the use of approaches based on functional prediction and pathways identification allow to putatively assign these species roles in the degradation of organic matter.
To understand metabolic networks in AD microbiomes, the metagenomic studies focus on analysis of the microbiomes’ composition, description of interactions within microbial communities, and assignment of functions in AD steps to specific groups of microbes [37, 51,52,53,54,55]. Metagenomic approaches show interdependencies between microbial communities and the type of feedstock, the C:N ratio, the operating conditions such as temperature, pH, bioreactor construction, organics loading rate, hydraulic retention time, etc. All the factors select and shape the structure of microbiomes, determine the metabolic pathway of methane production in biogas reactors, and influence the efficiency of methane production. It is believed that feedstock is a factor of special significance and it has a great influence on species abundance and on their functional activity [45,46,47,48,49,50, 56, 57].
Although the studies mentioned above explored the influence of feedstock, it is not completely understood how short-chain fatty acids can influence the last step of organic matter conversion to methane. The aim of the present study is to determine the influence of the common products of the acidogenic step (i.e., butyrate, propionate, lactate and acetate) on the metabolic pathways of methane formation and the microbial community composition. In accordance with the current trends, we combined three independent approaches: metagenomics, stable carbon isotope analysis of the fermentation gas, and monitoring of performance of Up-flow Anaerobic Sludge Blanket (UASB) bioreactors, to describe the complexity of acetogenesis and methanogenesis.
Object of the study
The objects of this study were four methane-yielding microbial communities designated M1, M2, M3, M4 in 3.5-L Up-flow Anaerobic Sludge Blanket (UASB) bioreactors (Fig. 1a, b) operating for 128 weeks (896 days) and processing artificial media intended to imitate a mixture of non-gaseous products of different types of acid fermentations with the domination of one of lactate (M1), butyrate (M2), propionate (M3) or acetate (M4). The methanogenic inoculum was anaerobic sludge from a municipal waste treatment plant “Warszawa Południe” in Warsaw, Poland, sampled during the winter (pH = 7.4). The artificial media were based on a modified M9 medium (BD Diffco)  without MgSO4, CaCl2 and glucose (Table 1). The medium was supplemented with the Bacto yeast extract (BD), sodium lactate (VWR Chemicals), butyric acid (Sigma Aldrich)/sodium butyrate (Alfa Aesar), propionic acid (Sigma Aldrich)/sodium propionate (Alfa Aesar), and acetic acid (VWR Chemicals)/sodium acetate trihydrate (Chempur) as shown in Table 1. The concentration of nitrogen-containing compounds in the media were comparable to those in the studies where co-cultures of acetogens and methanogens were examined [59, 60]. At the beginning, all the UASB bioreactors were filled with 1.5 L of the methanogenic inoculum and 2 L of the respective artificial medium (Variant 1) (pH = 7), and incubated for two weeks at room temperature (20–25 °C). Until the 100th day of bioreactors’ operation (week 15), the pH was adjusted with calcium hydroxide (5 g/L). Starting from the 17th day (week 3), the respective media (Fig. 1c, Table 1) were supplied to the UASB bioreactors using a peristaltic pump (ZALIMP, Poland) such that the hydraulic retention time was 7 days. In the 320th day of the process (week 46), 500 mL of the microbial community (anaerobic sludge) was removed and supplemented by methane-yielding sludge from a 50-L UASB bioreactor processing an acidic effluent from molasses fermentation  originally inoculated with the same above-described anaerobic sludge from the municipal waste treatment plant “Warszawa Południe” in Warsaw.
An initial operation period of 667 days (96 weeks) was aimed at adaptation to the substrates, stabilization and optimization of the process, elaboration of optimal concentration of mineral salts and yeast extract in the media, and as a result obtaining substrate-specific microbial communities. The idea of long-term operation of the process was based on (i) studies showing long-term systems producing methane [61, 62]; and (ii) natural and built anaerobic environments with permanent, stable methane production.
From the 610th day (week 88) to the end of bioreactors’ operation, the supplied medium contained 0.5 g/L of yeast extract and Na2HPO4, KH2PO4, NaCl and NH4Cl (Variant 3 in Table 1a) in amounts corresponding to a fourfold dilution of mineral compounds in the M9 medium. From the beginning till the 701st day (week 101) of bioreactors’ operation, the media contained lactate, propionate, acetate and butyrate with a 70% predominance of one of them (Fig. 1c, Table 1).
Starting from the 701st day of bioreactors’ operation, two main experiments (Experiment 1 and Experiment 2) were carried out. Experiment 1 was preceded by a 9-day period (from the 701st to the 709th days, week 102), during which the medium flow was switched off, and the COD of the effluent from bioreactors dropped below 100 mgO2/L, with the aim to minimize the concentrations of organic components with an unknown 13C isotope content. In the Experiment 1 (from the 710th to the 751st days, weeks 102–108), the bioreactors were again supplied with the media containing 70% of sodium salts of lactic (M1), butyric (M2), propionic (M3) or acetic (M4) acids (Table 1b). Carbon isotopic compositions of the media were determined (as described in the section “Isotopic analyses”). Between Experiment 1 and Experiment 2, there was a 47-day period (from the 752nd to the 798th days, weeks 108–114) when again the medium flow was switched off to decrease the COD of the effluent from the bioreactors to a value below 100 mgO2/L. In the Experiment 2 (from the 799th to the 886th days, weeks 115–127), the media contained only one, previously dominant component, i.e., sodium lactate (M1), sodium butyrate (M2), sodium propionate (M3), and sodium acetate (M4) (Table 1b). Carbon isotopic compositions of the media were again determined. Both experiments ended with collecting of methanogenic sludge samples for DNA analyses on the 751st day (week 108) for Experiment 1 and on the 884th day (week 127) for Experiment 2.
During the entire operation, physicochemical parameters (pH, COD—chemical oxygen demand, biogas production rate, CH4 and CO2 content) describing the UASB bioreactor performance were measured. Additionally, in the periods of Experiments 1 and 2, samples of biogas and effluents from bioreactors were collected for stable carbon isotope analyses (see the section “Isotopic analyses”) and to measure the concentrations of short fatty acids, sulfide (S2–), total nitrogen and soluble iron (Fe2+/3+) in the effluent from the bioreactors.
The pH of the media and the effluents from the UASB bioreactors as well as the redox potential inside the bioreactors was measured using a standard pH meter (ELMETRON model CP-502) equipped with a combination ORP (redox, mV) electrode type ERPt-13. The chemical oxygen demand (COD) values of the media and the effluents were determined using a NANOCOLOR COD 1500 kit (Macherey–Nagel) according to ISO 15705:2002.
The total rate of gas production was measured using a MGC-1 MilliGascounter (RITTER). The measurement data were recorded manually. The composition of the fermentation gas was analyzed using a HPR20 mass spectrometer (Hiden, England) with a QGA version 1.37 and a Fisons Gas Chromatograph 8000 series with a thermal conductivity detector.
Short-chain fatty acids were analyzed by HPLC with photometric detection (Waters HPLC system with Waters 2996—Photodiode Array Detector, and 300 × 7.8 mm Aminex HPX-87 H column with guard column at 30 °C). The samples were eluted for 45 min with an isocratic flow (0.6 mL/min) of 4 mM sulfuric acid.
The concentration of sulfide (S2–) in the effluents was determined using a NANOCOLOR SULFID 3 kit (test No. 985 073, Macherey–Nagel) according to DIN 38405-D26/27. The concentration of total nitrogen in the effluents was determined using a NANOCOLOR TNb 60 kit (test No. 985 092, Macherey–Nagel) based on ISO 7890-1. The concentration of soluble iron (Fe2+/3+) in the effluents was determined using a NANOCOLOR IRON 3 kit (test No. 985 037, Macherey–Nagel). Effluents were centrifuged before the analyses to remove microbial cells and debris.
Data from all analyses performed on samples collected from the UASB reactors are presented as mean values ± SD (standard deviation).
XLSTAT software by Addinsoft was used to prepare the box-plot figures and calculate Kolmogorov–Smirnov and Mann–Whitney tests for detailed presentation of COD reduction and methane production over time and testing similarities between the samples for the periods of Experiments 1 and 2.
Total DNA isolation and sequencing
Total DNAs from the microbial communities formed in all UASB reactors were isolated from samples of anaerobic sludge taken on the 751st day (week 108) for Experiment 1 and on the 884th day (week 127) for Experiment 2. There were three independent sample collections from across the whole bioreactor using a glass pipette with a broken tip. The samples of the sludge were immediately frozen in liquid nitrogen and, before DNA isolation, they were combined in one tube and mixed. DNA was extracted and purified using a DNeasy PowerSoil kit (Qiagen) in five replicates according to the manufacturer's protocol with some modifications. Five samples (0.3 g each) were placed into five bead tubes for extraction, incubated at 70 °C for 15 min, and shaken horizontally in a MoBio vortex adapter for 15 min at maximum speed. The remaining steps were performed as directed by the manufacturer. DNA quality was checked by running the sample on 1% agarose gel and template quantity was measured by a fluorimeter using Qubit 3.0 and High Sensitivity Picogreen reagents (Thermo, USA). Final samples of DNA extracted from the five replicates were pooled and stored at − 80 °C.
DNA sequencing: DNA was mechanically sheared using Covaris (Covaris, MA, USA) and sequencing libraries were constructed using KAPA Library Preparation kit (KAPA Biosystems, Wilmington, USA). The libraries were quality-checked using KAPA Library Quantification kit (KAPA-Roche, Basel, Switzerland), pooled in equimolar ratio and sequenced on a NextSeq 550 instrument using the NextSeq HighOutput reagent v2.5 (300 cycle) chemistry kit (Illumina, San Diego, USA).
Illumina DNA sequences were filtered using Trimmomatic software (v0.39) . Analysis of the shotgun reads was performed to determine the microbial composition using MetaPhlAn3 (v3.0.7–0) , while the heat-map visualization was generated with “hclust2.py” software. The relative abundance profiles obtained from MetaPhlAn3 and a Newick tree relating all the species included in MetaPhlAn 3.0 were used to calculate the weighted UniFrac distances with the script “calculate_unifrac.R”. Estimates of microbial community diversity was calculated on the filtered reads using Nonpareil3 software (v3.3.3–1) . Filtered reads were co-assembled with MEGAHIT (v1.2.4-beta) . All the scaffolds shorter than 1 kb were removed from the assembly and assembly statistics were determined using QUAST_1 . N50 and N90 values are statistics of a set of contigs or scaffolds lengths. The N50 value is calculated by first ordering every contig/scaffold by length starting from the longest. Then, starting from the longest contig/scaffold, the length of each is summed until this running sum is equal to one-half (or 90% for N90) of the total length of all contigs/scaffolds in the assembly. Metagenome binning was performed with metabat2 (v2.12) . The quality of the metagenome-assembled genomes (MAGs) was determined using CheckM software (v1.0.3)  and evaluated considering the MAG quality standards proposed by the Genomic Standards Consortium . The average nucleotide identity (ANI) was calculated considering the genomes deposited in the NCBI Reference Sequence Database . The putative taxonomic classification of the unclassified MAGs were further assessed based on ubiquitous proteins with GTDB-Tk (v1.0.2) . The genomes hits having ANI higher than 95% threshold were used for MAGs classification at the species level [73, 74]. The number of recovered MAGs reported in the global AD microbiome database  was estimated calculating the ANI values with dRep software (v2.3.2) . Genes were predicted for each MAG using Prodigal (v2.6.3) and annotated using a combination of strategies based on eggNOG-mapper (v2.0.1–1) , DRAM (v1.1.1–0)  and KEGG pathway analysis . MAGs coverage was determined through sequences recruitment using Bowtie 2 (v2.2.4)  and checkM (v1.0.3) .
Raw sequence data were uploaded to the Sequence Read Archive (SRA, NCBI) with project ID PRJNA680596.
Analyses of stable carbon isotope composition of fermentation gas and substrates
Samples of the fermentation gases were collected from the UASB bioreactors using a syringe and injected into 20-ml glass ampoules filled with a supersaturated NaCl water solution. Stable carbon isotope analyses of carbon dioxide and methane were carried out with an on-line method on a Delta V Advantage Mass Spectrometer coupled with a Trace GC Ultra gas chromatograph with a GC Isolink device (Thermo Scientific). The GC column used for gas analyses was a HP-PLOT/Q (Agilent Technologies, dimensions: 30 m × 0.32 mm × 20 µm). Helium was used as the carrier gas. The GC oven was initially held at 30 °C for 4 min, then heated at a rate of 10 °C/min to 210 °C, and held for 4 min. A CO2 certified gas standard (δ13CVPDB = − 36.2‰, Air Liquide Deutschland, GmbH) was used for calibration. A gas with the known carbon isotopic composition was analyzed regularly to check the accuracy of the measurement.
Stable carbon isotope analyses of substrates (ingredients in the fermentation experiment) were carried out using two analytical techniques: an off-line preparative system and continuous flow. The first analytical technique involved the use of about 2–5 mg or µl (where relevant) of pure substrates, which were combusted using a CuO wire in a sealed ca. 10 cm2 quartz tube, under vacuum at 900 °C. The produced CO2 gas was cryogenically purified off-line (liquid nitrogen and dry ice + ethanol mixture). The purified gas was introduced into an isotope ratio mass spectrometer (IRMS; Delta V Advantage/dual inlet, Thermo Scientific) for the analysis of the stable carbon isotope ratio.
The continuous flow technique was performed with a Thermo Finnigan Elemental Analyzer interfaced via a Conflo IV to a Finnigan Delta V Advantage (EA–IRMS). About 400 µg was weighted and placed into a tin capsule, sealed and packed using a hand-press device. The EA operated at an oxidation furnace temperature of 1020 °C, reduction furnace temperature of 650 °C, and a packed-column temperature of 45 °C .
For the normalization of the δ13C values, international standards (NBS22, USGS24 and USGS40 distributed by the International Atomic Energy Agency, Vienna) were used and then the values were recalculated and reported relative the Vienna Pee Dee Belemnite (V-PDB) scale with ± 0.1‰ precision [81, 82].
Performance of the UASB bioreactors M1–M4
To examine lactate, butyrate, propionate and acetate transformation to methane and carbon dioxide, four methane-producing microbial communities continuously processing, respectively, lactate-, butyrate-, propionate- and acetate-rich artificial media in UASB reactors were studied. Composition of the growth media was designed to imitate the acidic products cocktail of microbial communities in environments where, respectively, lactate, butyrate, propionate and acetate dominate in the acidogenic step of AD. Both the acetogenic and methanogenic steps took place in the bioreactors. Similar to a previous study , methane-yielding microbial communities were used instead of pure cultures of microorganisms. Our system is novel and closer to the natural environments and biogas digesters where microbial communities and not pure cultures exist.
Figure 2 and Additional file 1 show the performance (pH of the effluent, daily biogas production, content of methane in biogas, substrate utilization measured as COD reduction, and methane production per g of substrate’s COD reduction) of bioreactors M1–M4 operated for 128 weeks, whereas Figs. 3, 4 and Table 2 report more detailed data from operation of bioreactors M1–M4 in the periods of Experiment 1 and Experiment 2, including the composition of the effluents from the bioreactors.
Analysis of the physicochemical parameters revealed that pH of the effluents from the UASB bioreactors was close to neutral, which indicates stability of the methanogenesis process. Notice that the pH of the supplied media was neutral till the 100th day (week 15) of bioreactors’ operation. From the 101st day (week 15) of bioreactors’ operation till the beginning of Experiment 2, acidic media (pH ≈ 5) were supplied (Table 1a). The pH of the media during Experiment 2 was 7 (Tables 1 and 2). The optimum pH range for methane production is 6.8—7.2 . Thus, two-step or multi-step anaerobic digestion ensures pH stability, which is a commonly described advantage of such systems [5, 6].
Because the concentration of mineral salts in the substrate subjected to methanogenesis influences the methanogenic process, the periods of stabilization and optimization involved supplying undiluted and diluted (four-, ten- and fivefold) mineral compounds of the M9 medium to the bioreactors. Also, different concentrations of the yeast extract (0.5 and 2 g/L) were tested (Table 1a). In all bioreactors, (i) production of biogas containing 70–90% of methane and (ii) substrate utilization measured as % reduction of the COD of substrates were observed (Fig. 2b–d and Table 2). Interestingly, the highest methane content was found in the biogas formed in the bioreactor fed with the substrate rich in acetate (M4), whereas the lowest one in the biogas produced as the result of processing of the substrate rich in lactate. To take into account differences in CODs of the media and to standardize the results of methane generation, we calculated the amount of methane produced per g of COD reduction (Fig. 2d) in different periods of bioreactors’ operation. We concluded that in the experimental system used, no spectacular differences in methane production were observed for the various medium variants shown in Table 1.
Since the periods of Experiments 1 and 2 were the most important for the study, their precise analysis is presented below with reference to Figs. 3, 4 and Table 3. The usage of the substrate during Experiment 1 was on average ~ 70%. However, over time, the analysis revealed that the acetate-containing medium (M4) was utilized more efficiently in comparison to the butyrate- (M2), and propionate-containing (M3) media (p = 0.005 and p < 0.03, respectively). Differences in the substrate usage were much more evident during Experiment 2. Lactate- (M1) and acetate-containing (M4) media were utilized at a similar level (on average ~ 80%), whereas lower values were observed for propionate (M3) and butyrate (M4) (~ 70% and ~ 60%, respectively). Over time, the analysis revealed differences between M1 vs M2 and M1 vs M3 (p < 0.001) as well as between M4 vs M2 and M4 vs M3 (p ≤ 0.001). Also propionate (M3) was metabolized more efficiently than butyrate (M2), p ≤ 0.025 (Table 2, Fig. 3, Additional file 2).
Methane contents in the produced biogas during Experiment 1 were on average 75%, 82%, 81% and 76% for M1, M2, M3 and M4, respectively. Methane production was calculated per gram of reduction of the substrate COD. In Experiment 1, the results were on average 0.4 ± 0.1, 0.6 ± 0.1, 0.4 ± 0.1 and 0.5 ± 0.0 dm3/g COD reduction for bioreactors M1, M2, M3 and M4, respectively. Interestingly, over time, the analysis revealed that the highest efficiency was achieved for the M2 microbial community processing the butyrate-containing medium (in comparison to M1 and M3, p ≤ 0.03; in comparison to M4, p ≤ 0.005). It is noteworthy that the butyrate-containing medium was utilized in 60% only. In Experiment 2, when the media contained exclusively sodium lactate (M1), sodium butyrate (M2), sodium propionate (M3) and sodium acetate (M4), the methane contents were on average 74%, 79%, 84% and 88%, respectively. Interestingly, the highest number was observed in the bioreactor (M4) processing exclusively acetate. The results of methane production were on average 0.5 ± 0.1, 0.4 ± 0.1, 0.4 ± 0.1, 0.4 ± 0.1 dm3/g COD reduction for bioreactors M1, M2, M3, M4, respectively. In this case, the analysis over time found a statistically significant difference only between the lactate-utilizing (M1) and the acetate-utilizing (M4) microbial communities (p < 0.02) (Table 2, Fig. 3, Additional file 2).
Qualitative and quantitative analyses of the effluents from bioreactors (Fig. 4) showed that lactate was most effectively and almost completely utilized in all bioreactors independent of its initial concentration in the media. The average percent of butyrate and propionate utilization was 88.5 and 72.2, respectively, when they were the dominant components of the media. Butyrate, propionate and lactate were oxidized to acetate. Lactate was also metabolized to propionate (Fig. 4).
In the effluents from bioreactors, the S2− and Fe2+/3+ ions were detected in low concentrations of ≤ 0.2 mg/L and ≤ 1.5 mg/L, respectively (Table 2). It indicated that the processes of sulfate and iron reduction were insignificant, due to the low redox potential in bioreactors, as required for methanogenesis processes (Table 2).
To find differences between the metabolic pathways of transformation of lactate, butyrate, propionate and acetate to methane and carbon dioxide, isotopic analyses were conducted as a method independent of both monitoring of the bioreactor’s performance and metagenomics.
The carbon isotope signatures (δ13CSUB) of the substrates used in Experiment 1 and Experiment 2 (organic salts: sodium lactate, sodium butyrate, sodium propionate, sodium acetate) were determined prior to incubation (Table 3). The range of δ13CSUB values in the substrates used in this study was from − 45.4 to − 23.0‰, which is similar to those in the natural environment [84, 85]. The distribution of carbon isotopes during decomposition of any organic compound is controlled by the isotope signature of the substrate, the isotopic mass balance controlled by the redox processes, and the kinetic isotope effects [86,87,88]—in this study particularly relating to microbial degradation of organic salts. In the natural environment, the δ13C(CH4) values typical for acetate fermentation are in the range from − 60 to − 33‰, while for CO2 reduction in range from − 110 to − 60‰ (e.g., [87, 89, 90]). Combined analyses of δ13C(CH4) and δ13C(CO2) values in biogas from incubation experiments allow to calculate the potential isotopic fractionation factor α13CCO2-CH4 (e.g., [86, 91, 92]) according to the equation: α13CCO2-CH4 = (δ13CCO2 + 1000)/(δ13CCH4 + 1000) [87, 93]. The values of α13CCO2-CH4 in range 1.049–1.095 are typical for CO2 reduction, 1.039–1.058 are typical for acetate fermentation, and 1.005–1.03 are typical for methane oxidation . Therefore, tracing changes of the fractionation factor α13CCO2-CH4 in time is a helpful tool to determine the dominant process responsible for the methanogenesis.
The results of the mean δ13C values in CH4 and CO2 and of the calculated fractionation factor are shown in Table 3. The mean δ13C(CH4) values obtained for the biogas from the bioreactors were in the range from − 57.2 to − 36.2‰. Such a range of δ13C(CH4) values is typical for acetate fermentation and partially for mixing of methane from acetate fermentation and CO2 reduction. In experiments with lactate domination (M1, Experiments 1 and 2), a similar isotopic fractionation of carbon isotopes between CH4 and CO2 occurred with α13CCO2-CH4 equalling 1.043 and 1.042, respectively. The lactate was oxidized immediately to acetate (data from HPLC not shown here), which was the direct precursor for CH4 and CO2 formation. In experiments with butyrate domination (M2, Experiments 1 and 2) and propionate domination (M3, Experiments 1 and 2) we observed an enrichment of light carbon isotopes in CH4. The mean δ13C(CH4) values obtained for the biogas from bioreactors M2 and M3 were in range from − 57.0 to − 46.1‰. The reason for the enrichment of light carbon in CH4 is the influence of CO2 reduction processes during incubation. This observation is confirmed by the increase of the isotopic fractionation factor α13CCO2-CH4 in time in cases of butyrate (M2) and propionate (M3) domination, both in Experiment 1 (Fig. 5) and Experiment 2 (Fig. 6).
In the experiments with acetate domination (M4, Experiments 1 and 2), the mean δ13C(CH4) values were − 53.3 and − 57.2‰, respectively. The δ13C(CO2) values in the experiments with acetate (M4) differed significantly between Experiment 1 and Experiment 2 and were equal to − 20.9 and − 33.4‰, respectively. The isotopic fractionation factor α13CCO2-CH4 for incubations with acetate domination (M4) was equal to 1.034 in Experiment 1 and 1.026 in Experiment 2. Such values are close to those typical for acetate fermentation and/or methane oxidation. Processes of methane oxidation were excluded during incubation, because free oxygen was not present in the headspace gas samples (data from gas chromatography not shown here). Calculation of the isotopic mass balance with the equation δ13Ccalc = 0.5 × δ13CCH4 + 0.5 × δ 13CCO2 is a very useful tool, especially in the case of the methanogenic decomposition of acetate . Assuming that acetate was the main substrate for the methanogenesis in bioreactor M4, and based on the mean values of δ13C(CH4) and δ13C(CO2), the δ13Ccalc value should be similar or close to the δ13C value for the sodium acetate used in the experiments. In the case of the experiments with acetate domination (M4, Experiments 1 and 2), the calculated δ13C value was equal to − 37.1 and − 45.7‰, respectively. The δ13C for the sodium acetate used in this study was − 45.4‰. The calculations of the isotopic mass balance for Experiment 2 confirmed the decomposition of sodium acetate as the only substrate in the bioreactor. The enrichment of light carbon isotopes in the CO2 is an effect of the access of the microorganisms to a large pool of readily available acetate, continuously transferred to the bioreactors. The microorganisms decompose the 13C-depleted molecules of acetate slightly faster, which is typical for most of biodegradation reactions [86, 87, 89, 95]. In Experiment 1 with domination of acetate (M4), lactate was an additional source for decomposition and, therefore, formation of an additional pool of acetate . The δ13C value for the lactate used in this study was equal to -23.0‰ and its decomposition with sodium acetate (δ13C = − 45.4‰) resulted in the shift to a heavier values of the mean δ13C(CH4) in Experiment 1 with acetate domination (M4). The isotopic fractionation factor α13CCO2-CH4 for Experiment 2 with domination of acetate (M4) apparently might indicate methane oxidation, but it is probably not the case, as the incubation took place in an open system with continuous stable supply of new portions of external acetate to the fermentation solution (Additional file 3).
To explain the different contribution of the acetotrophic and hydrogenotrophic pathways of methane production depending on the type of substrate, metagenomic analyses of the microbiomes from M1 to M4 were performed.
Shotgun Illumina sequencing was performed on eight samples (four samples M1–M4 from Experiment 1, four samples M1–M4 from Experiment 2), obtaining from 9.5 to 33.6 Gbases of sequences depending on the experiment. Taxonomic analysis of the shotgun reads (unassembled) was performed to have a global representation of the microbiome composition including the rare components. Weighted UniFrac distances calculation revealed a clear separation between samples supplemented with lactate (M1), and the remaining samples, in particular those supplemented with acetate (M4) (Additional file 4). All the samples were dominated by Euryarchaeota (average 79.2% relative abundance), while bacteria were 20.8% on average (Additional file 4); in all the reactors the relative abundance of Archaea increased even more in the second period. This behavior was particularly evident in samples supplemented with acetate. In total, eight bacterial phyla were identified, and the most abundant are different in the reactors: Actinobacteria are abundant in butyrate M2 (Experiment 1), Bacteroidetes in lactate M1 (Experiment 2) and acetate M4 (Experiment 1), Synergistetes in butyrate M2 (Experiment 2). Analysis of the shotgun reads cannot provide a high level of detail in defining associations between functional pathways and microbial species. For this reason a genome-centric approach was implemented.
Reads belonging to all the experiments were co-assembled obtaining in total 291,272 contigs larger than 1 kb and having a total size of 1,057 Mbp. Contigs larger than 10 kbp accounted for approximately 41% of the total, evidencing the good quality of the assembly process.
Binning process recovered 234 Metagenome Assembled Genomes (MAGs) having completeness higher than 50% and contamination lower than 10% (Fig. 7; Additional file 5); 104 of them were of very high quality with completeness higher than 90% and contamination lower than 5%, other 90 MAGs have medium–high quality with completeness higher than 70% and contamination lower than 10%.
Taxonomy assignment of MAGs was performed taking into account the results from taxonomic informative proteins, and also Average Nucleotide Identity with the genomes present in NCBI database (Fig. 8; Additional file 6); these analyses revealed the presence of 31 archaeal (13.2%) and 203 bacterial species (86.8%). The high number of Archaea is an interesting feature of this metagenomic analysis project.
Reads alignment on the assembled scaffolds highlighted that MAGs accounted for 54–67% of the entire microbial community (depending on the experiment) and evidenced that the microbiome is extremely complex in terms of species number. It was not really expected that a microbiome focused on the last part of the AD food chain was so complex. The metric of sequence diversity calculated on the shotgun not-assembled reads ranged from 18.1 (bioreactor M3, Experiment 2) to 19.8 (bioreactor M2, Experiment 1) and evidenced that the complexity of the AD microbiome is similar to that calculated for human stool samples . The values for all the samples recovered from the first sampling point were higher than the corresponding values for the second sampling point; this is expected since a simplification of the feedstock resulted in a reduction of the number of species present.
By taking into account MAGs relative abundance for each experiment, only 8–15 of them (depending on the sample considered) had an abundance higher than 1% (Fig. 9; Additional file 7). Interestingly, five MAGs were identified at high abundance in all the samples (MAG_138, MAG_62, MAG_42, MAG_33, MAG_218). Considering the average value in all the experiments, MAG_138 (Methanothrix soehngenii AS29adLBPA_138) was the most abundant (12.1%) and reached the highest relative abundance in the reactor fed with acetate (18–19%); this archaeal species is highly similar to Methanothrix soehngenii GP6 (Average Nucleotide Identity-ANI 99%) an acetoclastic methanogen previously named Methanosaeta concilii GP6 (Additional file 6). Two archaeal species were also found with a high average abundance, MAG_62 (7.58%) (Methanoculleus sp. AS29adLBPA_62, 96.9% ANI) and MAG_103 (2.03%) (Methanocorpusculum sp. AS29adLBPA_103). Methanoculleus sp. AS29adLBPA_62 reached a very high abundance when propionate was used as feeding, but it was present at high abundance in all the samples. Methanocorpusculum sp. AS29adLBPA_103 had a more “scattered” distribution since it was more abundant in samples from bioreactor M1 (domination of lactate, Experiment 1) and bioreactor M2 (domination of butyrate, Experiment 2).
Considering bacterial species, the most abundant was MAG_218 (2.15%), an unknown species of the Spirochaetaceae family (Spirochaetaceae sp. AS29adLBPA_218), followed by MAG_33 (1.92%) (Bacteroidales sp. AS29adLBPA_33) which is highly similar to Lentimicrobium sp. 002433025 (98.52% ANI) and by MAG_149 (1.64%) (Synergistales sp. AS29adLBPA_149) belonging to the Thermovirgaceae family. MAG AS29adLBPA_218 reached the highest abundance in reactors fed with propionate and in the one fed with pure lactate, while MAG AS29adLBPA_149 was abundant in the reactors fed with pure acetate and butyrate.
Substrate-specific microbial species
The average relative abundance was calculated for all the MAGs and the values obtained for each sample/condition were compared with it (Additional file 7). This revealed the changes in abundance of each species in different samples and helped to identify microbes having substrate preferences. Due to the difficulty in defining preferences for substrate utilization for the MAGs, those specifically enriched in samples fed with one substrate were selected. It was found that 37 MAGs were enriched when acetate was used as feeding (> 1%, Bacteroidales sp. AS29adLBPA_29 and Methanomassiliicoccaceae sp. AS29adLBPA_85); 47 were enriched on butyrate (> 1%, Bacteroidales sp. AS29adLBPA_70, Firmicutes sp. AS29adLBPA_139, Synergistales sp. AS29adLBPA_149, Syntrophomonas sp. AS29adLBPA_86), 8 on propionate (> 1%, Peptococcaceae sp. AS29adLBPA_161), and 29 on lactate (> 1%, Methanofollis liminatans AS29adLBPA_137; Verrucomicrobia sp. AS29adLBPA_71; Methanocorpusculum sp. AS29adLBPA_103) (Fig. 10, Additional file 7). These microbes are to a certain extent “specialized” since they have preferences for one substrate; however, in most cases, their relative abundance is low, below 1%. Based on these findings, it seems that butyrate and acetate are the feeding substrates favoring more the growth of substrate-specific MAGs. On the contrary, 113 MAGs did not have preferences for one substrate and they had a more “generalistic” behavior and, thus, can be regarded as a core microbiome (Fig. 10, Additional file 7).
To identify the microbial species, for which the increased abundance in specific reactors is due to their ability to utilize or degrade specific compounds, the coverage profile of MAGs was compared with their predicted functional role, which was determined using a combination of predictive methods (including analysis of KEGG pathways).
Metabolic pathways present in the MAGs enriched in reactors fed with butyrate revealed that ten microbes have a complete pathway involved in SCFA and alcohol conversions “Butyrate, pt 1” (Additional file 8). Among these species, two were present at very high abundance in the microbiome: Synergistales sp. AS29adLBPA_149 and Bacteroidales sp. AS29adLBPA_70. The first one is a high-quality MAG (91.5% compl., 0.0% cont.) reaching 3.7% relative abundance when butyrate was the only carbon substrate provided. It has 21 complete KEGG modules including “fatty acid biosynthesis, elongation” (M00083) and “beta-Oxidation, acyl-CoA synthesis” (M00086) (Additional file 8). It is able to use butyrate catalyzing the conversion with butyrate kinase (K00929) and phosphate butyryltransferase (K00634). The second MAG has very similar characteristics (94.5% compl., 1.0% cont.) but the relative abundance with butyrate as feeding was slightly lower (1.2%). This MAG had a higher number of complete KEGG modules (N = 37) than Synergistales sp. AS29adLBPA_149, and this is probably due to its larger genome size (2.56 Mbp in comparison to 2.04 Mbp). These characteristics, suggestive of a more complex metabolism, still included genes for butyrate utilization (buk, K00929; ptb, K00625,K00634) and complete modules like M00083 and M00086. Only two MAGs were predicted to encode genes for a complete Wood–Ljungdahl pathway, Anaerolineales sp. AS29adLBPA_61 (compl. 93.6%, cont. 3.01%) and Deltaproteobacteria sp. AS29adLBPA_53 (compl. 87.8% cont. 3.2%). Both are enriched in samples fed with butyrate, but Deltaproteobacteria sp. AS29adLBPA_53 has a stronger enrichment; however, these species remain at a quite low level in all the samples (maximum 0.16% and 0.37% relative abundance). Both these MAGs have also the ability to convert butyrate to Acetyl-CoA, which can enter the W-L pathway; however, the low relative abundance suggests that their impact on acetate conversion is limited.
Analysis of the functional profile of MAGs enriched in the reactor fed with lactate revealed that five of them have a complete “SCFA and alcohol conversion: lactate D” pathway (Additional file 8). Despite these microbes being enriched in the reactor supplemented with lactate (Additional file 7), all these species have quite low relative abundance (Additional file 6). Among these, Lachnospiraceae sp. AS29adLBPA_171 (compl. 98.7%, cont. 1.5%) is interesting despite the low relative abundance (0.39%). This species encoded the electron transfer flavoprotein subunits A/B (K03521; K03522) (49957_15; 49957_16). The same pathway was described for the anaerobic lactate oxidation pathway used by the acetogen Acetobacterium woodii . This enzyme complex works in a stable complex with FAD-dependent lactate dehydrogenase LDH, possibly encoded in the same transcriptional unit (49957_13; K00104, K03777). Additionally, the RNFA-G subunits were all identified in the genome and present in the same transcriptional unit (277210_85 to 277210_90); this complex can drive ferredoxin reduction with NADH as reductant. The other MAGs having high abundance belong to Archaea (e.g., Methanocorpusculum sp. AS29adLBPA_103; comp. 97.5%, cont. 4.9%, Methanofollis liminatans AS29adLBPA_137; comp. 96.8%, cont. 2.6%) or to species not having functional pathways directly related to lactate utilization (e.g., Verrucomicrobia sp. AS29adLBPA_71) (comp. 93.9%, cont. 2.03%). AS29adLBPA_103, AS29adLBPA_137 and AS29adLBPA_71 reached 6.4%, 3.4% and 1.05% relative abundance, respectively.
Among the eight MAGs enriched in the samples fed with propionate, only Deltaproteobacteria sp. AS29adLBPA_60 (comp. 85.8%, cont. 1.9%) has a complete “SCFA and alcohol conversions: Propionate, pt 2” pathway (Additional file 8). This species has a low relative abundance (maximum 0.64%), and a manual verification of the gene content revealed that it harbors all the genes previously reported and belonging to the complex pathway involved in the conversion from propionate to acetate  (excluding succinyl-CoA synthetase, which absence could be due to limitations related to genome reconstruction). The conversion from propionate to succinyl-CoA is based on the path described in “KEGG propanoate metabolism”, followed by the conversion to pyruvate (KEGG TCA cycle) and finally to acetate (KEGG Carbon metabolism).
Other two dominant MAGs in the reactor fed with propionate are the unclassified Candidatus Cloacimonetes sp. AS29adLBPA_112 and the Peptococcaceae sp. AS29adLBPA_161. Both have a high-quality genome (compl. 98.9%, cont. 0%; compl. 94.3%, cont. 4.7%) and a maximum relative abundance close to 1%. As previously reported, the genomic analysis of the uncultivable “Candidatus Cloacamonas acidaminovorans” revealed the presence of all the genes involved in syntrophic propionate degradation. We can also suggest that Cloacimonetes sp. AS29adLBPA_112 have the same metabolism despite two genes in the pathway (succinate dehydrogenase and malate dehydrogenase) were not identified. Among the Archaea, Methanospirillum sp. AS29adLBPA_21 was enriched when propionate was the dominant substrate (Additional file 7). Previous data obtained from culture collections revealed that different isolated propionate‐degrading bacteria can act syntrophically with Methanospirillum hungatei .
Acetate addition led to the enrichment of many Archaea (Additional file 7). Some of these microbes, for example Methanosarcina mazei AS29adLBPA_30 are well known for their ability to perform acetoclastic methanogenesis  some others, such as Methanomassiliicoccaceae sp. AS29adLBPA_85 (comp. 97.2%, cont. 5.9%, maximum relative abundance 2.4%) have complete pathways for the conversion of methanol to methane and a nearly complete pathway for the conversion of methylamines (Additional file 8). Since there are seven MAGs related to Methanomassiliicoccaceae having an increased abundance in the acetate-enriched medium, it is tempting to speculate that they use methanol generated by other microbes. It is known that acetogens such as Acetobacterium woodii can directly assimilate formate and methanol into the reductive acetyl-CoA pathway producing pyruvate . Under the conditions of high acetate concentration they can possibly reverse the process leading to the generation of methanol. However, this process is still not demonstrated and can be an interesting target for future investigations. Notably, there was an enrichment of Acetoanaerobium noterae AS29adLBPA_56 (compl. 97.2%, cont. 3.8%) in the reactor fed with acetate, but it always remained at very low abundance. It was reported that this species can produce acetate from H2 and CO2 , and it has a nearly complete (83%) Wood–Ljungdahl pathway; however, the high acetate concentration present in the reactor suggests this is not the case. This species can possibly revert the Wood–Ljungdahl pathway using acetate to generate other sub-products.
The most abundant MAG identified among those enriched under these conditions is Bacteroidales AS29adLBPA_29 (relative abundance 2.8%), that, despite the high genome quality (compl. 93.2%, cont. 6.3%) does not have pathways related to acetate utilization (e.g., the Wood–Ljungdahl). The pathways involved in acetate utilization in this highly abundant species remain to be clarified.
Effects of dominant products of acidogenesis: acetate, butyrate, lactate and propionate on methane formation
Using culture-independent techniques and a long-term system, we traced the processing of four dominant non-gaseous products of bacterial acidic fermentations to methane and carbon dioxide for better understanding of the metabolic pathways and syntrophic cooperation between microorganisms in the methane-yielding communities. Our system allows examining the acetogenic and methanogenic stages of anaerobic digestion and helps understanding microbial processes in multi-stage systems processing organic matter to biogas. However, anaerobic digestion of methanogenic sludge inside the bioreactor should also be considered.
The media subjected to methanogenesis were dominated by, or contained exclusively, one of the acidic fermentation products: lactate (M1), butyrate (M2), propionate (M3) or acetate (M4). Utilization of the substrate measured as % reduction of substrates COD revealed that acetate and lactate were used by the microbial communities more efficiently than butyrate and propionate, especially when the media contained exclusively one compound. With regard to the efficiency of methane production, the interpretation of the results is rather ambiguous. The highest methane production in Experiment 1 was achieved for butyrate processing with 60% substrate utilization; whereas in Experiment 2, for lactate processing with 74% substrate utilization. The tested substrates had a lower impact on the final performance of the bioreactors than expected. In contrast, large differences were observed in the results of isotopic analyses clearly showing that domination of acetotrophic or hydrogenotrophic pathways of methane synthesis is substrate dependent.
Acetate is a direct substrate for acetotrophic methanogens. It explains its efficient utilization and domination of the acetotrophic pathway of methanogenesis confirmed by isotopic analyses. Acetate detected in the effluents from bioreactor M4 in both experiments as well as from the other bioreactors could come from non-utilized substrate or from anaerobic digestion of methanogenic sludge inside the bioreactor. Acetate could also be oxidized to carbon dioxide and hydrogen. However, since acetate oxidation is an endoergic reaction (ΔG0' = + 94.9 kJ/reaction), it requires syntrophic cooperation with hydrogenotrophic methanogens, which results in ΔG0' = − 36.3 kJ/reaction. Oxidation of butyrate and propionate are also endoergic reactions, with ΔG0' = + 48.3 kJ/reaction and ΔG0' = + 76.0 kJ/reaction, respectively, that become thermodynamically favorable with hydrogenotrophic methanogens, with ΔG0' = − 17.3 kJ/reaction and ΔG0' = − 22.4 kJ/reaction, respectively . This explains the domination of hydrogenotrophic methanogenesis in M2 and M3 microbial communities revealed by isotopic analyses. Previously, we have concluded that lactate is oxidized mainly to acetate during the acetogenic step of AD and this includes the acetotrophic pathway of methanogenesis . The present study confirms our previous results. Compared to butyrate and propionate, lactate is the most efficiently used substrate. It can be explained by the thermodynamics of lactate oxidation reactions, as it was discussed previously . Briefly, lactate can be oxidized directly to acetate by many bacteria without the contribution of methanogenic Archaea according to the mechanism described for Acetobacterium woodii (ΔG0' = − 61 kJ/mol) . The fermentation of lactate to propionate is also a thermodynamically favorable reaction (ΔG0' = − 169.7 kJ/reaction for Desulfobulbus propionicus) . Propionate formation from lactate was also observed in this study. Formation of propionate probably induces the metabolic pathways of propionate oxidation. To summarize the results of the study, (i) butyrate and propionate are transformed mainly to acetate, while lactate is transformed mainly to acetate and propionate; (ii) oxidation of acetate and lactate determines the acetotrophic pathway of methanogenesis, whereas oxidation of butyrate and propionate determines the hydrogenotrophic pathway of methanogenesis. Propionate is one of the major intermediates in AD. It is estimated that 6–35% of methane can be produced from propionate. Propionate degradation to methane draws the attention of many researchers, because accumulation of propionate is often observed in bioreactors with poor methane production . Interestingly, inhibition of methanogenesis in the bioreactor fed with a propionate-rich medium was not confirmed in this study. It could indicate that accumulation of propionate should be considered an indicator, rather than a cause of disturbances in anaerobic digestion and biogas production. Furthermore, accumulation of propionate is related to changes in the anaerobic digestion process, unstable pH and temperature, overload of feedstock, too high a concentration of short-chain fatty acids and hydrogen partial pressure, improper reactor configuration and hydraulic retention time [22, 103, 104].
Lactate-, butyrate-, propionate- and acetate-selected microbial communities
Metagenomic analysis revealed that the dominant bacteria and archaea identified in this study are probably unknown, uncultivable species. A striking distinctive feature is the high number of Archaea found in the examined microbial communities, 31 archaeal species out of 234 microbial species identified (13.2%). A recent study performed metagenomic binning starting from a range of different biogas reactors but, out of 1635 microbial species identified, only 61 Archaea (3.7%) were recovered . However, other studies revealed a remarkable number of archaeal species, which can also represent a large fraction of the entire microbiome [44, 105]. To determine how many MAGs reported in the present study represented new species not previously reported in the AD biogas database, a comprehensive Average Nucleotide Identity calculation was performed. Interestingly, 75.6% (177) of the MAGs identified here were not present in the global AD biogas database reported by Campanaro and colleagues, evidencing a very high level of new species recovered in this study. It can be speculated that the use of four diverse SCFAs as feedstock promoted the growth of species associated with the terminal part of the anaerobic degradation food chain (methanogenesis) where Archaea play a crucial role in methane production. On the other hand, most of the previous studies performed on the AD microbiome (those investigated by Campanaro and colleagues) did not focus on the acetogenic and methanogenic steps of the biogas food chain, where there are still many unknown microbes that were identified in the present project. As a confirmation of this finding, 20 out of 31 archaeal species identified here were not present in the global AD biogas database.
113 MAGs exhibited a comparable relative abundance in all bioreactors. They can be considered a sort of core microbiome with some dominant microbes such as Methanothrix soehngenii, Methanoculleus sp., representatives of Bacteroidales, Spirochaetaceae and unclassified bacteria. Interestingly, substrate-specific bacteria are not predominant in the microbial communities. Furthermore, the recognized acetate, lactate, butyrate, and propionate oxidizers described in the introduction are not in their majority present in the studied microbial communities form the bioreactors. The exceptions are butyrate-specific Syntrophomonas wolfei AS29adLBPA_159 (MQ) or lactate-specific Desulfovibrio desulfuricans AS29adLBPA_31. This is a confirmation that the vast majority of microorganisms, especially those requiring syntrophic growth and/or those having a very low growth rate, are yet to be isolated and cultivated. However, they represent the orders, families, and genera to which belong some recognized syntrophic bacteria involved in oxidation of products derived from acidic fermentations. These latter species were identified at high relative abundance in all bioreactors (e.g., Spirochaetaceae sp. AS29adLBPA_218, Bacteroidales sp. AS29adLBPA_33 or Synergistales sp. AS29adLBPA_149). Some other microbes were classified as “lactate-specific”, such as for example Deltaproteobacteria sp. AS29adLBPA_110 (MQ), Firmicutes sp. AS29adLBPA_210, Clostridiales sp. AS29adLBPA_49 (MQ); butyrate-specific Syntrophomonas sp. AS29adLBPA_81 (MQ), Synergistales sp. AS29adLBPA_149, Syntrophaceae sp. AS29adLBPA_197; propionate-specific Peptococcaceae sp. AS29adLBPA_161, Desulfobacteraceae sp. AS29adLBPA_204; acetate-specific Clostridiales sp. AS29adLBPA_6, Sphaerochaeta sp. AS29adLBPA_43 (MQ), Geobacter sp. AS29adLBPA_1 (MQ), Anaerolineaceae sp. AS29adLBPA_207.
Comments on substrate-specific and rare methanogenic species are reported in the Results section; additionally, functional analysis of medium-quality (MQ) MAGs has to be considered with caution, but their enrichment in specific samples can provide suggestions on their putative role.
Preference for specific substrates can be determined based on many different indicators, one being for sure the possibility to degrade a specific compound, but also syntrophic behaviors determined by compounds exchange between microbes are involved in shaping the structure of the microbiomes . For this reason, it is expected that only a part of the “specialized microbes” is able to directly utilize the substrate provided; other species probably rely on the chemical compounds released. It can simply indicate a wide capability of one bacterium to use several different substrates, which should also be illustrated by investigation of genes expression (metatranscriptomics). Differences in genes expression seem to be indirectly supported by the results of stable carbon isotope composition of methane and carbon dioxide in the fermentation gas. Metagenomic analysis of the microbial communities clearly shows that dominant hydrogenotrophic or acetotrophic pathways of methane formation elucidated by isotopic analysis are not associated with changes in the contribution of methanogens utilizing hydrogen and carbon dioxide or acetate. The most dominant methanogens in all the bioreactors are Methanothrix soehngenii AS29adLBPA_138 (10–19%) and Methanoculleus AS29adLBPA_62 (5.5–11.8%), the latter being a Medium-Quality (MQ) MAG. Methanoculleus, similarly to other hydrogenotrophic methanogens, produced methane from carbon dioxide and hydrogen generated during syntrophic oxidation of SCFAs, this assumption being based on taxonomic assignment since the recovered genome of AS29adLBPA_62 was MQ and did not undergo functional analysis. It may indicate induction of genes of either the hydrogenotrophic or acetotrophic pathways of methane formation depending on the supplied substrate. Our explanation is as follows. Acetate as well as lactate, which is easily metabolized to acetate via the mechanism described for A. woodii , determines the acetoclastic pathway of methane formation in Methanothrix soehngenii. Butyrate and propionate initiate metabolic pathways found for methane-yielding microbial communities fed with ethanol, where the Methanothrix species reduced carbon dioxide to methane with electrons accepted via DIET . This hypothesis should be confirmed by metatranscriptomic analysis of the microbial communities. Since, in addition to DNA, RNAs have been isolated from the same samples, the work on gene expression in the examined microbial communities is ongoing and the results will be included in a future report.
Recently, DIET is being increasingly mentioned in the context of anaerobic digestion [10, 106]. The microbial community fed with ethanol and dominated by the Methanotrix species exhibited an elevated abundance of the bacterial pilA gene and the methanogenic sludge showed a higher conductivity. In another study Methanospirillum hungatei was shown to form electrically conductive filaments that are analogs of e-pilli in Geobacter species. In the case of M. hungatei, it is the archaellum, whose core consists of tightly packed phenylalanines . In our study the Methanospirillum_AS29adLBPA_21 strain was classified as a propionate-specific methanogen in the microbial community fed with the propionate-rich substrate.
In the study by Barua and co-workers (2018) , addition of conductive carbon fibers to the bioreactors fed with butyrate- and propionate-containing media resulted in (i) an increase of methane production, (ii) a higher efficiency of substrate utilization, (iii) an increased contribution of electroconductive bacteria such as Desulfuromonas, Pseudomonas, Azonexus or Azovibrio that accompanied butyrate and propionate oxidizers, and (iv) a domination of Methanoseta species among the methanogens. The results indicated that DIET is involved in processing of propionate and butyrate by the microbial community.
In the present research, it was shown that the dominant components in the media (lactate, acetate, propionate or butyrate) subjected to methanogenesis moderately modified the final effect of bioreactor performance in terms of methane production and substrate utilization, whereas strongly affected the methanogenic pathways. Isotopic analysis evidenced different contributions of acetotrophic and hydrogenotrophic pathways for methane production, i.e., acetate and lactate favored the acetotrophic pathway, whereas propionate and butyrate favored the hydrogenotrophic pathway. Most of the 234 MAGs (31 archaeal and 203 bacterial species) were identified as new species. The core microbiome is represented by five MAGs present in high relative abundance (two methanogens: Methanothrix soehngenii and Methanoculleus sp., three bacterial MAGs classified only at high taxonomic level) and 108 other MAGs with a low relative abundance. Considering the relative abundance and/or their predicted functional role (determined according to KEGG pathways), three MAGs were found as propionate specific; four MAGs as lactate specific; four MAGs as butyrate-specific; and three MAGs as acetate-specific microbes. Analyzing the core microbiome and the substrate-specific species, we hypothesize the substrate may first of all change the metabolic activity of the bacteria/methanogens, rather than the composition of the microbial community. This requires confirmation using metatranscriptomic analysis.
Interestingly, we did not observe a reduction in methane production in the bioreactor fed with the propionate-rich medium. It may indicate that propionate commonly connected with inhibition of methanogenesis is rather an indicator than the cause of disturbances in anaerobic digestion. All these findings are relevant due to the fact that lactate, acetate, propionate and butyrate are the universal products of acidogenesis, regardless of feedstock.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files]. The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Barua S, Zakaria BS, Dhar BR. Enhanced methanogenic co-degradation of propionate and butyrate by anaerobic microbiome enriched on conductive carbon fibers. Bioresour Technol. 2018;266:259–66.
Liu Y, Whitman W. Metabolic, phylogenetic, and ecological diversity of the methanogenic Archaea. Ann NY Acad Sci. 2008;1125:171–89.
Detman A, Chojnacka A, Błaszczyk M, Kaźmierczak W, Piotrowski J, Sikora A. Biohydrogen and biomethane (Biogas) production in the consecutive stages of anaerobic digestion of molasses. Pol J Environ Stud. 2017;26(3):1023–9.
Thauer R, Kaster A-K, Seedorf H, Buckel W, Hedderich R. Methanogenic archaea: ecologically relevant differences in energy conservation. Nat Rev Micro. 2008;6:579–91.
Rabii A, Aldin S, Dahman Y, Elbeshbishy E. A review on anaerobic co-digestion with a focus on the microbial populations and the effect of multi-stage digester configuration. Energies. 2019;12:1106.
Van DP, Fujiwara T, Leu Tho B, Song Toan PP, Hoang Minh G. A review of anaerobic digestion systems for biodegradable waste: configurations, operating parameters, and current trends. Environ Eng Res . 2020;25(1):1–17.
Schink BSA. Syntrophism among prokaryotes. In: Dworkin M, Falkow S, Rosenberg E, Schleifer K-H, Stackebrandt E, editors. The Prokaryotes. 3rd ed. New York: Springer; 2006. p. 309–35.
Stams AJ, Plugge CM. Electron transfer in syntrophic communities of anaerobic bacteria and archaea. Nat Rev Microbiol. 2009;7(8):568–77.
Sieber JR, McInerney MJ, Gunsalus RP. Genomic insights into syntrophy: the paradigm for anaerobic metabolic cooperation. Annu Rev Microbiol. 2012;66:429–52.
Nozhevnikova AN, Russkova YI, Litti YV, Parshina SN, Zhuravleva EA, Nikitina AA. Syntrophy and interspecies electron transfer in methanogenic microbial communities. Microbiology. 2020;89(2):129–47.
Müller N, Worm P, Schink B, Stams AJ, Plugge CM. Syntrophic butyrate and propionate oxidation processes: from genomes to reaction mechanisms. Environ Microbiol Rep. 2010;2(4):489–99.
Worm P, Koehorst JJ, Visser M, Sedano-Núñez VT, Schaap PJ, Plugge CM, Sousa DZ, Stams AJM. A genomic view on syntrophic versus non-syntrophic lifestyle in anaerobic fatty acid degrading communities. Biochim Biophys Acta. 2014;1837(12):2004–16.
de Bok FA, Stams AJ, Dijkema C, Boone DR. Pathway of propionate oxidation by a syntrophic culture of Smithella propionica and Methanospirillum hungatei. Appl Environ Microbiol. 2001;67(4):1800–4.
Schmidt A, Frensch M, Schleheck D, Schink B, Müller N. Degradation of acetaldehyde and its precursors by Pelobacter carbinolicus and P. acetylenicus. PLoS ONE. 2014;9(12):e115902.
Meyer B, Kuehl J, Deutschbauer AM, Price MN, Arkin AP, Stahl DA. Variation among Desulfovibrio species in electron transfer systems used for syntrophic growth. J Bacteriol. 2013;195(5):990–1004.
Bertsch J, Siemund AL, Kremp F, Müller V. A novel route for ethanol oxidation in the acetogenic bacterium Acetobacterium woodii: the acetaldehyde/ethanol dehydrogenase pathway. Environ Microbiol. 2016;18(9):2913–22.
Weghoff MC, Bertsch J, Müller V. A novel mode of lactate metabolism in strictly anaerobic bacteria. Environ Microbiol. 2015;17(3):670–7.
Detman A, Mielecki D, Pleśniak Ł, Bucha M, Janiga M, Matyasik I, Chojnacka A, Jędrysek M-O, Błaszczyk MK, Sikora A. Methane-yielding microbial communities processing lactate-rich substrates: a piece of the anaerobic digestion puzzle. Biotechnol Biofuels. 2018;11:116.
Hatamoto M, Imachi H, Yashiro Y, Ohashi A, Harada H. Detection of active butyrate-degrading microorganisms in methanogenic sludges by RNA-based stable isotope probing. Appl Environ Microbiol. 2008;74(11):3610–4.
Liu P, Qiu Q, Lu Y. Syntrophomonadaceae-affiliated species as active butyrate-utilizing syntrophs in paddy field soil. Appl Environ Microbiol. 2011;77(11):3884–7.
Struchtemeyer CG, Duncan KE, McInerney MJ. Evidence for syntrophic butyrate metabolism under sulfate-reducing conditions in a hydrocarbon-contaminated aquifer. FEMS Microbiol Ecol. 2011;76(2):289–300.
Li J, Ban Q, Zhang L. Jha ak: syntrophic propionate degradation in anaerobic digestion: a review. Int J Agric Biol. 2012;14:843–50.
Schmidt A, Müller N, Schink B, Schleheck D. A proteomic view at the biochemistry of syntrophic butyrate oxidation in Syntrophomonas wolfei. PLoS ONE. 2013;8(2):e56905.
Liu Y, Balkwill DL, Aldrich HC, Drake GR, Boone DR. Characterization of the anaerobic propionate-degrading syntrophs Smithella propionica gen. nov., sp. Nov. and Syntrophobacter wolinii. Int J Syst Bacteriol. 1999;49:545–56.
Ito T, Yoshiguchi K, Ariesyady HD, Okabe S. Identification of a novel acetate-utilizing bacterium belonging to Synergistes group 4 in anaerobic digester sludge. ISME J. 2011;5(12):1844–56.
Lee SH, Park JH, Kim SH, Yu BJ, Yoon JJ, Park HD. Evidence of syntrophic acetate oxidation by Spirochaetes during anaerobic methane production. Bioresour Technol. 2015;190:543–9.
Limam RD, Chouari R, Mazéas L, Wu T-D, Li T, Grossin-Debattista J, Guerquin-Kern J-L, Saidi M, Landoulsi A, Sghir A, et al. Members of the uncultured bacterial candidate division WWE1 are implicated in anaerobic digestion of cellulose. MicrobiologyOpen. 2014;3(2):157–67.
Hedderich RWW. Physiology and biochemistry of the methane-producing Archaea. In: Dworkin M, Falkow S, Rosenberg E, Schleifer K-H, Stackebrandt E, editors. The prokaryotes. 3rd ed. New York: Springer; 2006. p. 1050–79.
Sakai S, Takaki Y, Shimamura S, Sekine M, Tajima T, Kosugi H, Ichikawa N, Tasumi E, Hiraki AT, Shimizu A, et al. Genome sequence of a mesophilic hydrogenotrophic methanogen Methanocella paludicola, the first cultivated representative of the order methanocellales. PLoS ONE. 2011;6(7):e22898.
Borrel G, Parisot N, Harris HMB, Peyretaillade E, Gaci N, Tottey W, Bardot O, Raymann K, Gribaldo S, Peyret P, et al. Comparative genomics highlights the unique biology of Methanomassiliicoccales, a Thermoplasmatales-related seventh order of methanogenic archaea that encodes pyrrolysine. BMC Genomics. 2014;15(1):679.
Borrel G, O’Toole PW, Harris HMB, Peyret P, Brugère J-F, Gribaldo S. Phylogenomic data support a seventh order of Methylotrophic methanogens and provide insights into the evolution of Methanogenesis. Genome Biol Evol. 2013;5(10):1769–80.
Borrel G, Adam PS, Gribaldo S. Methanogenesis and the wood-ljungdahl pathway: an ancient, versatile, and fragile association. Genome Biol Evol. 2016;8(6):1706–11.
Evans PN, Parks DH, Chadwick GL, Robbins SJ, Orphan VJ, Golding SD, Tyson GW. Methane metabolism in the archaeal phylum Bathyarchaeota revealed by genome-centric metagenomics. Science. 2015;350(6259):434–8.
Vanwonterghem I, Evans PN, Parks DH, Jensen PD, Woodcroft BJ, Hugenholtz P, Tyson GW. Methylotrophic methanogenesis discovered in the archaeal phylum Verstraetearchaeota. Nat Microbiol. 2016;1(12):16170.
Williams TA, Szöllősi GJ, Spang A, Foster PG, Heaps SE, Boussau B, Ettema TJG, Embley TM. Integrative modeling of gene and genome evolution roots the archaeal tree of life. Proc Natl Acad Sci. 2017;114(23):E4602–11.
Berghuis BA, Yu FB, Schulz F, Blainey PC, Woyke T, Quake SR. Hydrogenotrophic methanogenesis in archaeal phylum Verstraetearchaeota reveals the shared ancestry of all methanogens. Proc Natl Acad Sci. 2019;116(11):5037–44.
Campanaro S, Treu L, Rodriguez-R L, Kovalovszki A, Ziels R, Maus I, Zhu X, Kougias P, Basile A, Luo G, et al. The anaerobic digestion microbiome: a collection of 1600 metagenome-assembled genomes shows high species diversity related to methane production. bioRxiv. 2019;531:435.
Dziewit L, Pyzik A, Romaniuk K, Sobczak A, Szczesny P, Lipinski L, Bartosik D, Drewniak L. Novel molecular markers for the detection of methanogens and phylogenetic analyses of methanogenic communities. Front Microbiol. 2015;6:694.
Thauer RK. Biochemistry of methanogenesis: a tribute to Marjory Stephenson. 1998 Marjory Stephenson prize lecture. Microbiology. 1998;144:2377–406.
Zhao Z, Wang J, Li Y, Zhu T, Yu Q, Wang T, Liang S, Zhang Y. Why do DIETers like drinking: metagenomic analysis for methane and energy metabolism during anaerobic digestion with ethanol. Water Res. 2020;171:115425.
Yin Q, Miao J, Li B, Wu G. Enhancing electron transfer by ferroferric oxide during the anaerobic treatment of synthetic wastewater with mixed organic carbon. Int Biodeterior Biodegradation. 2017;119:104–10.
Kleinsteuber S. Metagenomics of methanogenic communities in anaerobic digesters. In: Stams A, editor. Biogenesis of hydrocarbons handbook of hydrocarbon and lipid microbiology. Cham: Springer; 2018. p. 1–23.
Liu Z, Ma A, Mathé E, Merling M, Ma Q, Liu B. Network analyses in microbiome based on high-throughput multi-omics data. Brief Bioinform. 2020. https://doi.org/10.1093/bib/bbaa005.
Wirth R, Böjti T, Lakatos G, Maróti G, Bagi Z, Rákhely G, Kovács KL. Characterization of core microbiomes and functional profiles of mesophilic anaerobic digesters fed with Chlorella vulgaris green microalgae and maize silage. Front Energy Res. 2019;7:111.
Cai M, Wilkins D, Chen J, Ng S-K, Lu H, Jia Y, Lee PKH. Metagenomic reconstruction of key anaerobic digestion pathways in municipal sludge and industrial wastewater biogas-producing systems. Front Microbiol. 2016;7:778.
Granada CE, Hasan C, Marder M, Konrad O, Vargas LK, Passaglia LMP, Giongo A, Oliveira RRD, Pereira LDM, Trindade FdJ et al: Biogas from slaughterhouse wastewater anaerobic digestion is driven by the archaeal family Methanobacteriaceae and bacterial families Porphyromonadaceae and Tissierellaceae. Renew Energy 2018, 118:840–846.
Delforno TP, Lacerda GV Jr, Sierra-Garcia IN, Okada DY, Macedo TZ, Varesche MB, Oliveira VM. Metagenomic analysis of the microbiome in three different bioreactor configurations applied to commercial laundry wastewater treatment. Sci Total Environ. 2017;587–588:389–98.
Campanaro S, Treu L, Kougias PG, De Francisci D, Valle G, Angelidaki I. Metagenomic analysis and functional characterization of the biogas microbiome using high throughput shotgun sequencing and a novel binning strategy. Biotechnol Biofuels. 2016;9(1):26.
Luo G, Fotidis IA, Angelidaki I. Comparative analysis of taxonomic, functional, and metabolic patterns of microbiomes from 14 full-scale biogas reactors by metagenomic sequencing and radioisotopic analysis. Biotechnol Biofuels. 2016;9(1):51.
Guo J, Peng Y, Ni B-J, Han X, Fan L, Yuan Z. Dissecting microbial community structure and methane-producing pathways of a full-scale anaerobic reactor digesting activated sludge from wastewater treatment by metagenomic sequencing. Microb Cell Fact. 2015;14:33–33.
Basile A, Campanaro S, Kovalovszki A, Zampieri G, Rossi A, Angelidaki I, Valle G, Treu L. Revealing metabolic mechanisms of interaction in the anaerobic digestion microbiome by flux balance analysis. Metab Eng. 2020;62:138–49.
Ortseifen V, Stolze Y, Maus I, Sczyrba A, Bremges A, Albaum SP, Jaenicke S, Fracowiak J, Pühler A, Schlüter A. An integrated metagenome and -proteome analysis of the microbial community residing in a biogas production plant. J Biotechnol. 2016;231:268–79.
Ziels RM, Svensson BH, Sundberg C, Larsson M, Karlsson A, Yekta SS. Microbial rRNA gene expression and co-occurrence profiles associate with biokinetics and elemental composition in full-scale anaerobic digesters. Microb Biotechnol. 2018;11(4):694–709.
Bremges A, Maus I, Belmann P, Eikmeyer F, Winkler A, Albersmeier A, Pühler A, Schlüter A, Sczyrba A. Deeply sequenced metagenome and metatranscriptome of a biogas-producing microbial community from an agricultural production-scale biogas plant. GigaScience. 2015. https://doi.org/10.1186/s13742-015-0073-6.
Theuerl S, Klang J, Heiermann M, De Vrieze J. Marker microbiome clusters are determined by operational parameters and specific key taxa combinations in anaerobic digestion. Bioresour Technol. 2018;263:128–35.
Stahel W. Circular economy. Nature. 2016;531:435–8.
Maus I, Koeck DE, Cibis KG, Hahnke S, Kim YS, Langer T, Kreubel J, Erhard M, Bremges A, Off S, et al. Unraveling the microbiome of a thermophilic biogas plant by metagenome and metatranscriptome analysis complemented by characterization of bacterial and archaeal isolates. Biotechnol Biofuels. 2016;9(1):171.
Miller J. Experiments in molecular genetics. New York: Cold Spring Harbor Laboratory; 1972.
Lee MJ, Zinder SH. Isolation and characterization of a thermophilic bacterium which oxidizes acetate in syntrophic association with a methanogen and which grows acetogenically on H(2)-CO(2). Appl Environ Microbiol. 1988;54(1):124.
Houwen FP, Dijkema C, Schoenmakers CHH, Stams AJM, Zehnder AJB. 13C-NMR study of propionate degradation by a methanogenic coculture. FEMS Microbiol Lett. 1987;41(3):269–74.
van der Star WRL, Abma WR, Blommers D, Mulder J-W, Tokutomi T, Strous M, Picioreanu C, van Loosdrecht MCM. Startup of reactors for anoxic ammonium oxidation: experiences from the first full-scale anammox reactor in Rotterdam. Water Res. 2007;41(18):4149–63.
Garcia ML, Angenent LT. Interaction between temperature and ammonia in mesophilic digesters for animal waste treatment. Water Res. 2009;43(9):2373–82.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Beghini F, McIver LJ, Blanco-Míguez A, Dubois L, Asnicar F, Maharjan S, Mailyan A, Thomas AM, Manghi P, Valles-Colomer M, et al. Integrating taxonomic, functional, and strain-level profiling of diverse microbial communities with bioBakery 3. bioRxiv. 2020;11:2500.
Rodriguez-R LM, Gunturu S, Tiedje JM, Cole JR, Konstantinidis KT. Nonpareil 3: fast estimation of metagenomic coverage and sequence diversity. mSystems. 2018;3(3):e00039-e118.
Li D, Luo R, Liu CM, Leung CM, Ting HF, Sadakane K, Yamashita H, Lam TW. MEGAHIT v1.0: a fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods. 2016;102:3–11.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29(8):1072–5.
Kang DD, Froula J, Egan R, Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ. 2015;3:e1165.
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25(7):1043–55.
Bowers RM, Kyrpides NC, Stepanauskas R, Harmon-Smith M, Doud D, Reddy TBK, Schulz F, Jarett J, Rivers AR, Eloe-Fadrosh EA, et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol. 2017;35(8):725–31.
O’Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, Rajput B, Robbertse B, Smith-White B, Ako-Adjei D, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016;44(D1):D733-745.
Chaumeil PA, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the genome taxonomy database. Bioinformatics. 2019. https://doi.org/10.1093/bioinformatics/btz848.
Rinke C, Schwientek P, Sczyrba A, Ivanova NN, Anderson IJ, Cheng JF, Darling A, Malfatti S, Swan BK, Gies EA, et al. Insights into the phylogeny and coding potential of microbial dark matter. Nature. 2013;499(7459):431–7.
Konstantinidis KT, Ramette A, Tiedje JM. The bacterial species definition in the genomic era. Philos Trans R Soc Lond B Biol Sci. 2006;361(1475):1929–40.
Olm MR, Brown CT, Brooks B, Banfield JF. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. Isme j. 2017;11(12):2864–8.
Huerta-Cepas J, Forslund K, Coelho LP, Szklarczyk D, Jensen LJ, von Mering C, Bork P. Fast genome-wide functional annotation through orthology assignment by eggNOG-mapper. Mol Biol Evol. 2017;34(8):2115–22.
Shaffer M, Borton MA, McGivern BB, Zayed AA, La Rosa SL, Solden LM, Liu P, Narrowe AB, Rodríguez-Ramos J, Bolduc B, et al. DRAM for distilling microbial metabolism to automate the curation of microbiome function. Nucleic Acids Res. 2020;48(16):8883–900.
Kanehisa M, Sato Y. KEGG Mapper for inferring cellular functions from protein sequences. Protein Sci. 2020;29(1):28–35.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.
Platzner IT, Habfast K, Walder AJ, Goetz A. Modern isotope ratio mass spectrometry. New York: Wiley; 2015.
Paul D, Skrzypek G, Fórizs I. Normalization of measured stable isotopic compositions to isotope reference scales–a review. Rapid Commun Mass Spectrom. 2007;21(18):3006–14.
Skrzypek G, Sadler R, Paul D. Error propagation in normalization of stable isotope data: a Monte Carlo analysis. Rapid Commun Mass Spectrom. 2010;24(18):2697–705.
Reungsang A, Pattra S, Sittijunda S. Optimization of key factors affecting methane production from acidic effluent coming from the sugarcane juice hydrogen fermentation process. Energies. 2012;5:4746–57.
Conrad R, Claus P, Chidthaisong A, Lu Y, Fernandez Scavino A, Liu Y, Angel R, Galand PE, Casper P, Guerin F, et al. Stable carbon isotope biogeochemistry of propionate and acetate in methanogenic soils and lake sediments. Org Geochem. 2014;73:1–7.
Londry KL, Des Marais DJ. Stable carbon isotope fractionation by sulfate-reducing bacteria. Appl Environ Microbiol. 2003;69(5):2942–9.
Jędrysek MO. Spatial and temporal patterns in diurnal variations of carbon isotope ratios of early-diagenetic methane from freshwater sediments. Chem Geol. 1999;159(1):241–62.
Whiticar MJ. Carbon and hydrogen isotope systematics of bacterial formation and oxidation of methane. Chem Geol. 1999;161(1):291–314.
Szynkiewicz A, Modelska M, Jędrysek M-O, Kurasiewicz M, Mastalerz M. Ageing of organic matter in incubated freshwater sediments; Inferences from C and H isotope ratios of methane. Geol Q. 2008;52:383-396.
Sugimoto A, Wada E. Carbon isotopic composition of bacterial methane in a soil incubation experiment: Contributions of acetate and CO2H2. Geochim Cosmochim Acta. 1993;57(16):4015–27.
Jedrysek MO. Carbon isotope evidence for diurnal variations in methanogenesis in freshwater lake sediments. Geochim Cosmochim Acta. 1995;59(3):557–61.
Krzycki JA, Kenealy WR, Deniro MJ, Zeikus JG. Stable Carbon Isotope Fractionation by Methanosarcina barkeri during Methanogenesis from Acetate, Methanol, or Carbon Dioxide-Hydrogen. Appl Environ Microbiol. 1987;53(10):2597–9.
Games LM, HayesRobert JM, Gunsalus P. Methane-producing bacteria: natural fractionations of the stable carbon isotopes. Geochim Cosmochim Acta. 1978;42(8):1295–7.
Chanton J. The effect of gas transport on the isotope signature of methane in wetlands. Org Geochem. 2005;36:753–68.
Corbett JE, Tfaily MM, Burdige DJ, Cooper WT, Glaser PH, Chanton JP. Partitioning pathways of CO2 production in peatlands with stable carbon isotopes. Biogeochemistry. 2013;114(1):327–40.
Wada E, Nakamura K. Carbon isotopic studies of global methane production with emphasis on paddy fields. In: Wroclaw MOJ, editor. Course-book on Isotope Geology. Wroclaw: Wroclaw University and Committe on Mineralogical Sciences; 1990. p. 141–9.
Schink B. Energetics of syntrophic cooperation in methanogenic degradation. Microbiol Mol Biol Rev. 1997;61(2):262–80.
Muyzer G, Stams AJM. The ecology and biotechnology of sulphate-reducing bacteria. Nat Rev Microbiol. 2008;6(6):441–54.
Cairó JJ, Clarens M, Touzel JP, Bardulet M, París JM. Methanosarcina mazei JC2, a new methanogenic strain isolated from lake sediments, that does not use H2/CO2. Microbiologia. 1992;8(1):21–31.
Cotton CAR, Claassens NJ, Benito-Vaquerizo S, Bar-Even A. Renewable methanol and formate as microbial feedstocks. Curr Opin Biotechnol. 2020;62:168–80.
Sleat R, Mah R, Robinson R. Acetoanaerobium noterae gen. nov., sp. Nov.: an anaerobic bacterium that forms acetate from H2 and CO2. Int J Syst Bacteriol. 1985;35:10–5.
Kamagata Y. Syntrophy in anaerobic digestion. Anaerobic Biotechnol. 2015. https://doi.org/10.1142/9781783267910_0002.
Oyekola OO, van Hille RP, Harrison ST. Study of anaerobic lactate metabolism under biosulfidogenic conditions. Water Res. 2009;43(14):3345–54.
Ahlert S, Zimmermann R, Ebling J, König H. Analysis of propionate-degrading consortia from agricultural biogas plants. MicrobiologyOpen. 2016;5(6):1027–37.
Fukuzaki S, Nishio N, Shobayashi M, Nagai S. Inhibition of the fermentation of propionate to methane by hydrogen, acetate, and propionate. Appl Environ Microbiol. 1990;56(3):719–23.
Kouzuma A, Tsutsumi M, Ishii Si, Ueno Y, Abe T, Watanabe K. Non-autotrophic methanogens dominate in anaerobic digesters. Sci Rep. 2017;7(1):1510.
Morita M, Malvankar NS, Franks AE, Summers ZM, Giloteaux L, Rotaru AE, Rotaru C, Lovley DR. Potential for direct interspecies electron transfer in methanogenic wastewater digester aggregates. MBio. 2011;2(4):e00159-e1111.
Walker DJF, Martz E, Holmes DE, Zhou Z, Nonnenmann SS, Lovley DR. The archaellum of Methanospirillum hungatei is electrically conductive. MBio. 2019;10(2):e00579-e1519.
We would like to thank Dr. Piotr Piela for editorial assistance. We would like to thank Dr. Damian Mielecki and Karolina Żuchniewicz for technical assistance. We would like to thank the director of the Municipal Water and Sewage Enterprise in the capital city of Warsaw, Poland who issued a permission to sample the anaerobic sludge and use it for scientific research.
ASi is a professor, AD, RG, JG are assistants, AG is a PhD student at the Institute of Biochemistry and Biophysics PAS, Warsaw, Poland. SC is a professor and LT is an assistant at the Department of Biology, University of Padova, Padova, Italy. EŁ is a professor and MB is an assistant at the Faculty of Earth Sciences, University of Silesia in Katowice, Poland. AC is an assistant and MJKB is a professor emeritus at the Department of Biochemistry and Microbiology, Institute of Biology, Warsaw University of Life Sciences, Warsaw, Poland. IM is a professor and MJan is an assistant at the Oil and Gas Institute, National Research Institute, Cracow, Poland. MOJ is a professor at the Institute of Geological Sciences, University of Wroclaw, Poland. ŁP, WD, MJak are assistants at the Institute of Geological Sciences, University of Wroclaw and at the Institute of Biochemistry and Biophysics PAS, Warsaw, Poland. ASa is an assistant at the Institute of Agricultural and Food Biotechnology, Warsaw, Poland.
We acknowledge the support of The National Science Centre, Poland, through grant UMO-2015/17/B/NZ9/01718 awarded for the years 2016–2019.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Anna Sikora and Stefano Campanaro contributed equally to this work.
Performance of bioreactors M1–M4: pH of effluents, substrate utilization and methane production—detailed data for Fig. 2.
Experimental data from individual measurements: δ 13C(CH4), δ 13C(CO2) and δ 13Ccalc.
Evaluation of microbial composition based on analysis of unassembled reads performed independently from binning, using MetaPhlAn3.
Quality of recovered MAGs according to results obtained from ckeckM software. Main parameters calculated for MAG genomes are also included.
Taxonomic results obtained for recovered MAGs. Results obtained using CAT/BAT and GTDB-Tk software are reported.
Relative abundance of MAGs calculated for different reactors and time points. MAGs enriched under different conditions are also reported.
Functional annotation of recovered MAGs obtained using DRAM software (Distilled and Refined Annotation of Metabolism). KEGG IDs of annotated genes are also reported.
About this article
Cite this article
Detman, A., Bucha, M., Treu, L. et al. Evaluation of acidogenesis products’ effect on biogas production performed with metagenomics and isotopic approaches. Biotechnol Biofuels 14, 125 (2021). https://doi.org/10.1186/s13068-021-01968-0
- Short-chain fatty acids
- Isotopic analysis
- Pathways of methanogenesis