Statement of Novelty

To introduce novel biomass types for biogas production, an integrated evaluation approach has to be applied. The present study delivers a combination of technological and microbial data about the process of anaerobic digestion of the commonly used substrates Zea mays L. silage and Miscanthus sacchariflorus silage, the latter of which has a high content of lignin and cellulose. The correlations between the abundance of methanogenic microorganisms and the concentrations of individual VFAs produced during digestion of these substrates were estimated. The results reveal novel associations between the type of substrate used, the rate of methane production and the structure of the methanogenic community. The use of materials with a high lignin content limits biodiversity in the anaerobic digesters, which can result in decreases in anaerobic digestion performance and biogas production.

Introduction

Anaerobic digestion (AD) is an important technology for simultaneously utilizing waste and producing biogas. Food waste, municipal solid waste and animal manure are commonly transformed in this manner. Recently, energy crops and crop residues have gained much attention as substrates, mainly due to their abundance. By using perennial crops that are cultivated on marginal and agro-technical lands that are not suitable for producing food, it may be possible to reduce the competition for land between food and energy crops, and mitigate the currently observed increase in demand from the bioenergy sector for vegetable food feedstocks.

AD of organics to methane is a complex process involving at least four different microbial groups: primary fermenting bacteria, secondary fermenting bacteria, and two types of methanogens. Hydrolysis and fermentation are the key stages in reducing the amount of the solid phase and supplying volatile fatty acids (VFAs), which are favorable precursors for the next steps. The pathways leading to formation of individual VFAs depend on several factors, such as the chemical composition of substrate and the operational conditions of digestion [1]. It is commonly believed that an increase in the concentration of VFAs can lead to process failure. One of the most common reasons for process deterioration is reactor acidification due to reactor overload. For this reason, many studies have focused on the relationship between VFAs accumulation and organic loading rate [2].

The effect of the chemical composition of biomass on the build-up of VFAs, particularly with regard to the presence of lignocellulosic materials in the biomass, has received little attention. It is commonly known that acetate and propionate are the main intermediates in the degradation of polysaccharides such as starch or cellulose. Valerate may form as intermediate of l-arabinose and d-xylose fermentation, which are present in the largest amount in hemicelluloses. Branched-chain fatty acids, such as isobutyrate, isovalerate and 2-methylbutyrate, are products of protein degradation [2].

In a complex methanogenic community, some volatile fatty acids can be synthesized from two acids that are present together, or reversible isomerization between isoforms can take place. For example, Wu et al. [3] found that, during conversion to methane, reversible isomerisation between butyrate and isobutyrate took place, and that when butyrate was degraded in the presence of propionate, isobutyrate was synthesized from butyrate and 2-methylbutyrate from propionate. A literature review shows that, in many cases, the appearance of acid isoforms in the digestate leads to the uncoupling of acid formation and conversion of the acids to methane.

In anaerobic digesters, methane production is favored by syntrophic interactions between fermenting bacteria and methanogens. Franke-Whittle et al. [4] found that, when a mixture of cow manure (46%), corn silage (36%), vegetable waste (9%) and potato (9%) was digested under mesophilic conditions, the quantity of Methanosarcina correlated positively with the concentrations of all VFAs, whereas the abundance of Methanobacterium was not affected by VFA concentration. However, when a mixture of cow manure (52%) and food waste (48%) was digested in thermophilic conditions, the abundance of Methanosarcina only had a significant positive correlation with the concentrations of acetate, butyrate and isovalerate, while Methanobacterium abundance had a significant negative correlation with the same VFAs. According to the authors, the concentrations of VFAs are not the only criteria that determine the dominance of a particular genera of methanogens. However, the results of those authors do suggest the need for further research that uses statistical analyses to estimate the association between the type and concentration of VFAs and the structure of the methanogenic community.

Therefore, the aim of this study was to investigate how two kinds of lignocellulosic biomass that differ in their content of cellulose and lignin influence the concentrations of individual VFAs and the structure of the methanogenic community. Silages of Z. mays L. and M. sacchariflorus were anaerobically digested under constant operational conditions. Z. mays L. was chosen as a model crop because it is commonly used in agricultural biogas plants, and M. sacchariflorus because of its high content of fibrous materials, including lignin. The following variables were measured: (i) the removal of individual non-fibrous and fibrous materials, (ii) the concentrations of individual VFAs that are fermentation intermediates, (iii) the Archaeal community structure in digesters, using the polymerase chain reaction coupled with denaturing gradient gel electrophoresis (PCR-DGGE), and (iv) the rate of biogas and methane production. The observed and theoretical methane yield coefficients were calculated, and Pearson product-moment correlation coefficients were calculated to estimate the associations between the concentrations of individual VFAs and the abundance of methanogens.

Methods

Feed Material

Zea mays L. and M. sacchariflorus were used as substrates for biogas production. They were obtained during field experiments at the Production and Experimental Station in Bałcyny (53°35′49″N, 19°51′20.3″E), belonging to the University of Warmia and Mazury in Olsztyn (Poland).

Zea mays L. (mid-early variety LG 3232) was harvested at the BBCH 89 stage (fully ripe: kernels hard and shiny), and M. sacchariflorus (clone variety) at the beginning of the ejection phase of the panicle in the fourth year of vegetative growth. The crops were collected by self-propelled harvesters equipped with cutting drums, which chopped the crops into 2–3 cm long pieces. Both kinds of harvested biomass were then ensiled by concentrating them in 200 L silos lined with foil for 90 days. Formic acid (85%) was added at a ratio of 5 g of acid per 1 kg of biomass.

Feedstock Preparation

For feedstock standardization, the silages were dried at 60 °C and ground in a cutting mill (Retsch SM100, Germany), then passed through a 1-mm screen and stored in plastic containers at room temperature. In the standardized silages, the concentrations of total solids (TS) and volatile solids (VS) were 968.0 (± 5.0) g kg−1 and 902.0 (± 4.0) g kg−1, respectively for Z. mays L., and 971.0 (± 3.0) g kg−1 and 864.0 (± 3.0) g kg−1, respectively, for M. sacchariflorus. The elemental composition of Z. mays L. silage included the following: C, 43.9%; H, 6.01%; N, 1.57%; and O, 41.7%; that of M. Sacchariflorus was as follows: C, 45.9%; H, 5.61%; N, 1.46%; and O, 36.0%. To prepare the feedstocks, 12.8 g of standardized Z. mays L. silage and 13.1 g of standardized M. sacchariflorus were weighed and mixed with 0.133 L distilled water to obtain a 9% VS concentration in both feedstocks.

Digester Performance

The experiments were conducted simultaneously over 75 days in continuously-stirred anaerobic tank reactors (CSTRs) with a working volume of 6 L, equipped with a stirrer and a water jacket. Valves located on the top and bottom of the reactor chamber enabled feeding and biogas/digestate collection. The reactors were inoculated with anaerobic sludge from the sludge digestion chambers of a municipal wastewater treatment plant. Once a day, 0.133 mL of digestate was withdrawn from each reactor and the same volume of feedstock (with 9% VS concentration) was then added. One reactor was fed with Z. mays L. feedstock (R-ZM) and the second was fed with M. sacchariflorus feedstock (R-MS). To avoid sudden pH drops shortly after feedstock was added and to prevent reactor acidification, NaOH was added to the feedstock at a concentration of 0.5 g L−1 of feedstock during the first 10 days of the experiment. Both digesters were operated at a hydraulic retention time (HRT) of 45 days and an organic loading rate (OLR) of 2.07 g L−1 day−1. The temperature was maintained at 39 °C.

Analytical Methods

TS and VS were measured according to standard methods for examination of water and wastewater [5]. Content of water soluble carbohydrates was measured by the anthrone method; that of raw protein (RP), by total organic nitrogen (the difference between total Kjeldahl nitrogen and total ammonium nitrogen) and multiplying by a factor of 6.25. Content of raw lipids (RL) was found by Soxhlet extraction using petroleum ether as a solvent; that of raw fibre (RF), according to EN ISO 6865:2000 [6], and that of neutral detergent fiber, according to EN ISO 16472:2006 [7].

Nitrogen-free extract (NfE) was measured by subtracting the amounts of RP, RL and RF from VS [8]. The Van Soest procedure was used to measure the content of cellulose (Cel), hemicellulose (Hem) and lignin. The Van Soest extension introduces three fractions: neutral detergent fibre (NDF), acid detergent fibre (ADF) and acid detergent lignin (ADL). NDF is mainly composed of cellulose, hemicellulose and lignin; the main fractions in ADF are cellulose, lignin and ash from cell walls; whereas ADL includes lignin and acid insoluble ash.

Starch content was estimated by taking the difference between total carbohydrates (as the sum of RF and NfE) and NDF; the content of hemicellulose was calculated as the difference between NDF and ADF; whereas cellulose content, as the difference between ADF and ADL [8]. Content of non-fibrous carbohydrates (NFC) was calculated by summing the contents of starch and water-soluble carbohydrates.

The C, H and N contents of the biomass were measured using a FLASH 2000 (Thermo Scientific). In the feedstocks, TS and VS were measured. In the liquid phase of the feedstocks, the following variables were measured: soluble chemical oxygen demand (CODs), ammonium and phosphates concentration (according to standard methods; [5]), pH (immediately after sampling using a pH meter Hanna HI 221, USA), total alkalinity (by titration to pH 4.3 with normalized 0.1 M HCl using a Schott titroline system), and VFAs (using a gas chromatograph GC, Varian 3800, Australia) according to Gilroyed et al. [9]. The total concentration of VFAs was the sum of the concentrations of the individually analyzed VFAs. CODs, pH, ammonium, alkalinity and VFAs were measured in filtered supernatant samples previously centrifuged at 8693×g for 10 min. The biogas volume was measured using a gas meter (Aalborg, model XFM 17). Methane and carbon dioxide content were measured in samples collected in Tedlar bags by a GA2000 + automatic analyzer (Geotechnical Instruments). All measurements were performed daily.

Molecular Analysis of the Structure of the Archaeal Community

DNA Extraction

Genomic DNA was extracted from biomass samples collected from both reactors at the beginning of the experiments (day 0), then on days 5, 8, 12, 15, 22, 33, 36, 50, 57 and 68, and at the end of the experiments (day 75). DNA was extracted from the samples by the following steps: 0.1 g of biomass was washed in sodium phosphate buffer (0.1 M; pH 8.0) and pelleted by centrifugation. After rejection of the supernatant, 1 mL of the extraction buffer (100 mM Tris–HCl; 100 mM EDTA; 1.5 M NaCl; pH 8) and 0.3 g of glass beads (Ø 0.25–0.5 mm; Carl Roth, Karlsruhe, Germany) were added. Samples were then shaken for 20 min in a bead-beating device (Uniequip, Germany). Next, samples were incubated for one additional hour at 65 °C with 0.2 mL of the 10% SDS solution to improve lysis of the Archaeal cells. After centrifugation at 16,100×g for 10 min, the DNA solution was purified in silica washing columns (A&A Biotechnology, Gdynia, Poland). The purified DNA was suspended in 50 µL of deionized, DNAase free water and stored at 20 °C.

Polymerase Chain Reaction

Archaeal 16S rRNA gene fragments were amplified with these primers: 0357F with an additional GC clamp, and 0505R, as described by Watanabe et al. [10]. PCR was performed in a Gene-Amp® PCR System 9700 (Applied Biosystems). 50 ng of extracted DNA were used for the PCR mixture, which contained 0.5 µM of each primer, 100 µM of deoxynucleoside triphosphate (Promega, Wisconsin, U.S.A.), 0.6 U of Hypernova DNA polymerase (DNA-Gdansk, Poland), 3 µL of reaction buffer (100 mM Tris–HCl, 500 mM KCl, 1.5% Triton X-100), 1.5 mM MgCl2 and sterile water, for a final volume of 30 µL. The temperature program for DNA amplification was as follows: 94 °C for 5 min; 35 cycles of denaturation at 94 °C for 45 s, annealing at 58 °C for 45 s, extension at 72 °C for 45 s, and a single final elongation at 72 °C for 10 min. The PCR products were checked via electrophoresis on 1.0% agarose gels stained with ethidium bromide. Each sample was amplified in triplicate.

Denaturing Gradient Gel Electrophoresis

Polymerase chain reaction products with a GC clamp were resolved on a 6% polyacrylamide gel (37.5:1 acrylamide:bisacrylamide) with a gradient ranging from 30 to 60% of denaturating agent (formamide and urea). Electrophoresis was performed for 15 h at 60 V in 1× TAE buffer (2 M Tris base, 2 M acetic acid, 0.05 M EDTA) using a Dcode system (Bio-Rad Laboratories Inc., Hercules, Canada). The temperature was set at 60 °C. Bands were visualized by staining with 1:10,000 SybrGold (Invitrogen) for 20 min followed by UV trans-illumination. Images were recorded and analyzed with KODAK 1D 3.6 Image Analysis Software (Eastman Kodak Company, USA). Based on the band intensity on the gel tracks, as indicated by the peak heights of the densitometric curves, a Shannon index of general diversity H [11] was calculated. The H value was calculated using the following equation: H = − ∑(ni/N) × log(ni/N); where ni is the height of the peak and N the sum of all peak heights of the densitometric curve.

DNA Sequencing

The 16S rRNA gene bands showing the greatest intensity were excised from the denaturing gel, transferred into 50 µL sterile water and frozen at − 20 °C for 24 h. The samples were thawed at room temperature, and the gel fragments were rubbed mechanically. Eluted PCR products were then re-amplified using the same primer set and purified with a Clean-up kit (A&A Biotechnology, Gdynia, Poland). Sequencing reactions were performed by Macrogen Europe (Amsterdam, Netherlands) with ABI3730XL (PE Applied Biosystems, Foster City, CA, USA). All reactions were run following the manufacturers’ protocols. The obtained sequences were aligned using the ClustalW program and analyzed with the BLAST program.

Statistical Analysis

To test the null hypothesis of no differences in concentrations of individual VFAs between R-ZM and R-MS, the Tukey HSD test was used (STATISTICA 9.0, StatSoft Inc.), with p < 0.05 set as the criteria for significance. Pearson product-moment correlation coefficients were computed to assess the associations between the characteristics of the digestate and the intensity of DNA bands, the biogas production rate and the CH4 content in the biogas (using STATISTICA 9.0).

Calculations

The efficiency of VS removal was calculated using the following equation Eq. (1):

$${\eta _{{\text{VS}}}}=\frac{{\left( {{\text{V}}{{\text{S}}_{{\text{feedstock}}}} - {\text{V}}{{\text{S}}_{{\text{digestate}}}}} \right) \times 100}}{{{\text{V}}{{\text{S}}_{{\text{feedstock}}}}}}$$
(1)

where ηVS is the efficiency of VS removal (%), and VSfeedstock and VSdigestate are the concentrations of VS (g kg−1) in the feedstock and the digestate, respectively.

The efficiency of removal of individual components of both silage feedstocks was calculated with Eq. (2):

$${\eta _i}=\frac{{\left( {{C_{i{\text{,feedstock}}}} - {C_{i{\text{,digestate}}}}} \right) \times 100}}{{{C_{i{\text{,feedstock}}}}}}$$
(2)

where ηi is removal efficiency of individual components in the feedstock (RP, RL, NFC, Cel, Hem and ADL) (%), Ci,feedstock and Ci,digestate are the concentrations (g kg−1) of individual components in the feedstock and the digestate, respectively.

The specific biogas production rate was calculated using Eq. (3):

$${r_{\text{B}}}=\frac{{{\nu _{\text{B}}}}}{V}$$
(3)

where rB is the specific biogas production rate (L L−1 day−1), νB is the daily biogas production (L day−1) and V is the working volume of the digester (L).

The observed methane yield was calculated using Eq. (4):

$${Y_{\text{M}}}=\frac{{{\nu _{\text{B}}} \times {C_{\text{M}}}}}{{V{S_{{\text{feedstock}}}} \times {\rho _{{\text{feedstock}}}} \times {V_{{\text{feedstock}}}} \times 100}}$$
(4)

where YM is observed methane yield (L g−1 VS), νB is the daily biogas production (L day−1), CM is the concentration of methane in biogas (%), VSfeedstock is the concentration of VS in the feedstock (g kg−1), ρfeedstock is the density of the feedstock (assumed to be ρfeedstock ≈ 1.0 kg L−1) and Vfeedstock is the volume of the feedstock added to the digester each day (L day−1).

The theoretical yield of methane (YM, theor) for silage feedstock was estimated using the stoichiometric O’Rourke’s equation, as described in detail in Klimiuk et al. [12]. Because only part of the VS are converted to methane, due to the presence of a non-biodegradable fraction in the feedstock, the amount of methane which theoretically could be produced in these experimental conditions (YM,VSrem) was calculated (Eq. 5) with the assumption that the efficiency of volatile solids removal (ηVS) corresponded to the concentration of the biodegradable fraction in the feedstock.

$${Y_{{\text{M,V}}{{\text{S}}_{{\text{rem}}}}}}=\frac{{{Y_{{\text{M,theor}}}} \times {\eta _{{\text{VS}}}}}}{{100}}$$
(5)

where YM, VSrem is theoretical methane yield estimated based on VS removal (L g−1 VS), and ηVS is the efficiency of VS removal (%).

Results and Discussion

Digestate Performance and Organics Removal Efficiency

The results of selected indicators of reactor performance are shown in Table 1. To maintain appropriate alkalinity and pH in R-ZM and R-MS, NaOH was added.

Table 1 Selected indicators characterizing anaerobic digestion performance

The concentration of ammonium nitrogen in silages was low, but both substrates contained proteins that underwent decomposition, resulting in ammonium nitrogen release to digestate. In R-ZM, the concentration of ammonium nitrogen was 117 mg L−1, whereas in R-MS, it was 120 mg L−1, both within the optimal range of 50–200 mg L−1 [13]. In R-ZM, the average concentration of VFAs was 2150 mg L−1, while in R-MS it was 3000 mg L−1.

Although the OLR was similar in both series, the organics (expressed as VS) were removed with higher efficiency during the digestion of Z. mays L. silage than during that of M. sacchariflorus. This was due to different loads of individual fibrous and non-fibrous components (Fig. 1). Z. mays L. silage contained a higher concentration of non-fibrous carbohydrates than M. sacchariflorus. In contrast, M. sacchariflorus had about two times more cellulose and four times more lignin.

Fig. 1
figure 1

Composition of the feedstock and percent removal of individual components with aZ. mays (R-ZM) L.; bM. sacchariflorus (R-MS). Error bars represent standard deviations. RP raw protein, RL raw lipids, NFC non-fibrous carbohydrates, Cel cellulose, Hem hemicellulose, ADL acid detergent lignin

Non-fibrous materials were removed more efficiently in R-ZM (RP, 55.6%; RL, 92.1%; NFC, 90.2%) than in R-MS (RP, 51.5%; RL, 48.1%; NFC, 47.8%). As for the fibrous materials, hemicellulose was removed less efficiently in R-ZM (68.2%) than in R-MS (74.1%). In contrast, cellulose was removed more efficiently in R-ZM (65.4%) than R-MS (58.2%). The removal efficiency probably differed because the higher content of lignin in M. sacchariflorus reduced cellulose hydrolysis by two mechanisms: providing a physical barrier that impeded or prevented cellulases from accessing the cellulose, and irreversibly adsorbing and binding these enzymes, thus preventing their action on cellulose. Mussatto et al. [14] have shown that the cellulose conversion ratio (defined as a sum of glucose yield and cellobiose yield) is four times higher after removing lignin and hemicellulose from brewer’s spent grain than before removing these substances.

Biogas/Methane Production

The specific biogas production rate (rB) and the percent concentration of methane in biogas are presented in Fig. 2.

Fig. 2
figure 2

Specific biogas production rate (rB), and CH4 and CO2 content in biogas during digestion of a and cZ. mays L. (R-ZM); b and dM. sacchariflorus (R-MS)

The specific biogas production rate of R-ZM increased over time, whereas that of R-MS was stable. The specific biogas production rate of R-MS was half that of R-ZM, indicating that M. sacchariflorus is a less efficient substrate. Moreover, during digestion of M. sacchariflorus, CO2 concentration increased, whereas that of CH4 decreased. The theoretical and calculated methane yields are presented in Table 2.

Table 2 Methane yield during anaerobic digestion of Z. mays L. (R-ZM) and M. sacchariflorus (R-MS)

It was found that the ratio of YM/YM,VSrem in R-ZM was 0.94, whereas in R-MS it was only 0.69. This means that, with M. sacchariflorus as a substrate, degradation of organic compounds and methane production were uncoupled. This uncoupling may have occurred because the methanogens were inhibited by components of M. sacchariflorus biomass, or by intermediate metabolites formed during its digestion, or both.

The M. sacchariflorus biomass in the present study contained a high concentration of lignin, which based on the measured difference between lignin concentration in the feedstock and that in the digestate, was removed during digestion with an efficiency of 17.7% (Fig. 1b). Although this is not a direct measurement of intermediates of lignin digestion, recent research has shown that several species of bacteria can degrade lignin under anaerobic conditions [15, 16]. Kato et al. [17] demonstrated that lignin-derived aromatics are degraded under methanogenic conditions due to cooperation among acetogens, ring-cleaving fermenters/syntrophs and acetoclastic/hydrogenotrophic methanogens. They showed that acetogenic bacteria enriched from rice paddy field soil can utilize the O-methyl groups on these aromatics, but they are unable to degrade the aromatic ring structure, which can be degraded by some sulfate-reducing bacteria and fermentative bacteria. Molinuevo-Salces et al. [18] found that anaerobic sludge from the anaerobic digester of a municipal wastewater treatment plant degraded up to 80% of lignin under mesophilic conditions when a mixture of swine manure and vegetable waste (50:50%) was used as feedstock. They also reported that the efficiencies of lignin removal from vegetable waste as sole feedstock and from this waste in a mixture with poultry litter (50:50%) were 50% and 35%, respectively, after 80 days of batch experiments. More recently, Candia-Garcia et al. [19] reported that 32.62% of lignin was removed from rice straw during 60 days of treatment at room temperature in a batch reactor inoculated with rumen content. These efficiencies of lignin removal in these reports are higher than that observed in the present study with M. sacchariflorus as a substrate. This difference could be due to the shorter HRT that was used here (45 day) and differences in the recalcitrance of lignin and other lignocellulosic feedstocks. For example, when Triolo et al. [20] studied the influence on biochemical methane potential of the lignin in energy crops and that in animal manure, both alone and in combination, they found that the lignin in energy crops is degraded to a lesser extent than that in animal manure.

It is well documented that the presence of a lignin fraction, especially a low molecular fraction of this substance, induces methanogenic toxicity and inhibits bacterial growth [21]. Rodriguez-Chiang et al. [22] evaluated the methane potential and biodegradability of acetate and lignin-rich pulp mill effluent at different ratios using mesophilic anaerobic sludge as inoculum. They found a negative linear correlation between lignin content and methane yield. The presence of lignin resulted in a 41% decrease in methane yield when the lignin fraction was increased by a factor of 3. A negative linear correlation between lignin content and biodegradability was also observed. The inhibiting effect of lignin derivatives on methanogenesis of lignocellulosic biomass was recently reported by Schroyen et al. [23], who found that phenolic compounds released from Miscanthus severely inhibited AD, thus lowering biogas production. Thus, the fact that, in the present study, VS removal efficiency was 1.5-times lower in R-MS than in R-ZM, and the observed methane yield was 33% lower than that which theoretically could be produced in these experimental conditions could be due to the inhibitory effect of M. sacchariflorus lignin.

VFAs Formation

The profiles of VFA concentrations in R-ZM and R-MS are shown in Fig. 3. With each silage, the concentrations of VFAs during AD were different. In R-ZM the concentration of propionate (Pro) was highest. The ratio of propionate to acetate (Ac) decreased with time, averaging 2.23. The concentration of valerate (Va) was about two times lower than that of isovalerate (i-Va). Isobutyrate (i-Bu) was found sporadically, and an increase in its concentration was accompanied by a decrease in the concentration of propionate. From highest to lowest, the average concentrations of individual VFAs were as follows: HPro > Hi−Va > (HAc ≅ HVa) > Hi−Bu; the differences were statistically significant (p < 0.05). The fact that these VFAs were found in these proportions is probably due to the rate constants for their degradation. VFAs are degraded according to first order kinetics, as reported by Wang et al. [24] after investigating AD of pretreated (ultrasonic disintegration, heating, and freezing) and untreated waste activated-sludge. They found that propionate is degraded more slowly than acetate (kHPro = 0.02878 h−1 and kHAc = 0.06119 h−1, respectively). Degradation of isobutyrate is slower than that of butyrate (Bu) (0.04385 h−1 and 0.16098 h−1, respectively); and isovalerate degradation is slower than that of valerate (0.02169 h−1 and 0.04504 h−1, respectively). These differences in rate constants might be due to the different structures of the isoforms.

Fig. 3
figure 3

Profiles of VFAs during digestion of aZ. mays L (R-ZM).; bM. sacchariflorus (R-MS); Ac acetate, Pro propionate, Bu butyrate, i-Bu isobutyrate, Va valerate, i-Va isovalerate

In R-MS, the concentration of propionate was also higher than that of acetate, but the ratio of propionate to acetate (1.26) was lower than in R-ZM. As with acetate, the isovalerate to acetate ratio was also 1.26 with in R-MS. The concentration of isovalerate did not differ significantly from the concentration of propionate. As the concentration of valerate increased or decreased, so did the concentrations of acetate and propionate. The amount of isobutyrate remained stable over time, averaging 368 mg L−1. The average concentrations of the VFAs were in the following order: (HPro ≅ Hi−Va) > HAc > (HVa ≅ Hi−Bu). In R-MS, the concentrations of isobutyrate and isovalerate were higher than those of butyrate and valerate. Moreover, concentrations of the isoforms of these VFAs were also higher than those in R-ZM (3.08- and 1.7-times higher, respectively). These AD intermediates have been shown to build up when methanogenesis is inhibited by coenzyme M during AD of ferulic acid (an intermediate of lignin degradation) [25]. Isovaleric and isobutyric acid are considered to be the most sensitive indicators of digestion imbalance [4, 26], and their occurrence in R-MS could have inhibited methanogenic activity. This could also explain why the observed methane yield was 31% lower than the methane yield that theoretically could be produced in these experimental conditions.

In both digesters, propionate predominated during the entire experiment. Degradation of propionate is known to be the most thermodynamically unfavorable step in the AD system. It is degraded under a relatively low partial pressure of hydrogen because only under these conditions is the process thermodynamically favorable [26]. Thus, the high concentration of propionate in the digestates from R-ZM and R-MS indicted that propionate degradation was inhibited by a high partial pressure of hydrogen, and also indicated that hydrogenotrophic methanogens were inhibited [27]. This result is consistent with the lower diversity of hydrogenotrophic methanogens in these two digesters (presented in the next section).

The Community Structure of Methanogens in Anaerobic Sludge as Determined by DGGE

In this study, the diversity and structure of the Archaeal community was determined by PCR-DGGE followed by phylogenetic analysis. DGGE fingerprints of samples from R-ZM and R-MS are shown in Fig. 4. In both series, 24 unique bands were detected by image processing software using the peak heights of the densitometric curves. In R-ZM, the values of the Shannon diversity index increased up to the fifteenth day of the experiment (Fig. 5). After this time, index values oscillated around an average value of 2.82 ± 0.22. In R-MS, in contrast, the Shannon index decreased from 2.88 to 2.31 between the twelfth and the last day of the experiment. Our study shows that the type of substrate affects the diversity of methanogens. The influence of substrate on methanogens diversity was previously shown by Dabrowska et al. [28]. In their work, the addition of glycerol and pig slurry to plant silage caused differences in methanogen community structures.

Fig. 4
figure 4

DNA fingerprints obtained by PCR-DGGE of Archaeal 16S rDNA showing Archaeal community changes during anaerobic digestion of: aZ. mays L. (R-ZM), bM. sacchariflorus (R-MS). The number below each lane shows the sampling day. The most intense bands are indicated by numbers, and their taxa were identified (details in Fig. 6.)

Fig. 5
figure 5

Changes in Archaeal diversity measured using Shannon’s index. Error bars show the standard deviation

Fig. 6
figure 6

Neighbor-joining tree showing the phylogenetic affiliation of DGGE band sequences

In order to identify the taxa to which the methanogens belonged, the bands with the highest intensity were excised from the DGGE gels and sequenced. The phylogenetic tree of the Archaea, based on 16S rRNA gene sequences, is presented in Fig. 6.

With both feedstocks, acetoclastic and hydrogenotrophic methanogens were present. Five bands contained sequences identifying them as belong to the Methanosaetaceae family, which has an affinity for acetate (D-03, D-13, D-06, D-05 and D-15). In addition, band D-02 included a sequence that matched that of Methanosarcina thermophila (JQ346758). Bands D-04 and D-01 were grouped in Methanomicrobiales, which are hydrogenotrophs. Also, DNA sequences derived from the D-10 and D-08 bands formed a distinct clade separated from Methanobacteriales and Methanococcales, and showed close affiliation with a large phylum of Euryarchaeota known as ArcI; whereas band D-08 matched an uncultured ArcI archaeon (CU917025). Microorganisms from the ArcI clade have been recognized as methanogenic, most likely utilizing either H2/CO2 or formate as a substrate for methanogenesis [29]. They are commonly found in wastewater sludge, sediment and other methanogenic environments.

Marquez et al. [30] studied the diversity of methanogen communities in three single batch digesters operated with marine biomass composed of 59% sea grasses and 41% seaweeds during 40 days under thalassic ambient conditions. The reactors were inoculated with cow manure, sea wrack-associated microflora or marine sediment. In all three digesters, hydrogenotrophic methanogens belonging to Methanobacteriales and Methanomicrobiales dominated in the methanogenic populations. Moreover, in the digester inoculated with marine sediment, mixotrophic Methanosarcinaceae co-dominated. Nettman et al. [31] analyzed methanogenic communities in six full-scale biogas plants operated under mesophilic conditions. Hydrogenotrophic Methanomicrobiales predominated in a digester fed with a mixture of 82% maize silage, 12% barley grain and 6% water, working at an HRT of 107 days and an OLR of 3.4 kg organic dry substance m−3 day−1. In contrast, Methanosaetaceae predominated in a digester fed with a mixture of cattle manure (76%), maize silage (13%), grass silage (5%), cattle dung (4%) and grain (2%), working at an HRT of 46.9 day and an OLR of 2.5 kg organic dry substance m−3 day−1. Both genera belong to the order Methanosarcinales, which are acetoclastic methanogens that have been reported to be responsible for approximately 70% of the methane produced in biogas reactors [17]. Members of the genus Methanosaeta have been reported to dominate in reactors with low levels of NH3 and VFAs [4], whereas Methanosarcina have high metabolic versatility and the ability to use acetate, hydrogen, formate, secondary alcohols and methyl compounds as energy sources [32].

In our research, band D-05, with a DNA sequence affiliated with Methanosaeta concilii (AB679168), was found in both R-ZM and R-MS throughout the experiments. This was probably due to the low ammonium nitrogen concentration, which did not exceed 166 mg NH4+ L−1 in both digesters. M. concilii was found to be the most ammonium-sensitive methanogen among a group of pure cultures; it was completely inhibited at a concentration of 560 mg (total) NH4+ L− 1 at a suboptimal pH of 7.0. At consistently higher ammonium concentrations, abundant M. concilii are replaced with other methanogens as the microbial community adapts during an operational period [33].

Correlations Between the Concentrations of Individual VFAs and the Abundance of Methanogens

It was assumed that shifts in Archaeal community structure were related to changes in the concentration of individual VFAs because HRT, OLR, pH, alkalinity and ammonium were held constant throughout the experiments. The intensity of the DGGE bands varied throughout the experiment (Fig. 4), and because PCR-DGGE is a ‘semi-quantitative’ approach, the differences over time and between groups allowed us to estimate the relative abundance of the Archaea. To evaluate the associations between their abundance and the concentrations of VFAs, Pearson product-moment correlation coefficients were calculated.

There was a weak negative correlation between the intensity of band D-08 (hydrogenotrophic methanogens affiliated with uncultured ArcI) and acetate concentration (r = − 0.6823, p = < 0.05) (Table 3) in R-ZM.

Table 3 Pearson correlations between VFA concentrations and DGGE band intensities during anaerobic digestion of Z. mays L. (R-ZM)

Surprisingly, there was a weak positive correlation (r = 0.6486, p < 0.05) between the acetate concentration and the intensity of band D-01, whose sequence matched Methanolinea mesophila, belonging to the hydrogenotrophic methanogens. This can happen when hydrogen-consuming methanogens participate in methane production pathways as partners of acetate-oxidizing bacteria. Syntrophic oxidation of acetate has been described by Zinder and Koch [34] as consisting of two reactions: first, oxidation of acetate to hydrogen and carbon dioxide, and second, conversion of hydrogen and carbon dioxide to methane by hydrogenotrophic methanogens. Thus, for successful acetate degradation by this process, the availability of hydrogen is most important. Hydrogen-consuming methanogens belonging to Methanoculleus sp. have been found in all investigations of organisms that participate in mesophilic syntrophic oxidation of acetate [35].

The intensity of band D-01, affiliated with hydrogenotrophic M. mesophila, and the concentration of propionate also correlated positively (r = 0.6848, p < 0.05) (Table 3). Propionate is degraded by acetogenotrophs to produce 1 mol of acetate, 1 mol of carbon dioxide and 3 moles of hydrogen.

Thus, it is not surprising that syntrophic interactions between hydrogen-utilizing M. mesophila and acetogenotrophs have been shown in the literature. Sakai et al. [36] isolated M. mesophila from an anaerobic, propionate-degradation enrichment culture that was originally established from a soil sample from a rice field in Taiwan.

The concentration of valerate correlated strongly with the intensity of band D-10 (r = 0.7708, p < 0.01) and weakly with band D-04 (r = 0.6912, p < 0.05), both of which are associated with hydrogenotrophic methanogens. This is probably because valerate is most often degraded via β-oxidation to 1 mol of acetate, 1 mol of propionate and 2 mol of hydrogen [37]. In contrast to the hydrogenotrophic methanogens, bands D-06, D-05 and D-03, with DNA sequences affiliated with acetoclastic methanogens, correlated strongly with the concentration of isovalerate. This “interspecies acetate transfer” may be another type of syntrophic cooperation between anaerobic bacteria, in which acetate has to be kept at minimum concentrations in the presence of hydrogen. The degradation product of isovalerate is mainly acetate. The reaction occurs with CO2 as a co-substrate, producing 3 mol of acetate and 1 mol of H2 [37].

The DGGE analysis of samples from R-MS suggested that the methanogenic community in these digesters is less diverse than that in R-ZM. (Fig. 4b). Two bands that showed high intensity in R-ZM, band D-15 (affiliated with acetoclastic Methanosaetaceae) and D-01 (affiliated with Methanoregulaceae), were not present in R-MS.

In R-MS, correlations between Methanosarcina abundance and the concentrations of all measured VFAs were statistically significant (Table 4). Franke-Whittle et al. [4] made the same observation when investigating the influence of VFA concentrations on methanogenic communities in mesophilic conditions. In our study, band D-05, closely related to Methanosaeta sp., also had a significant positive correlation with the concentrations of acetate, isobutyrate and isovalerate (Table 4). Roest et al. [38] reported that Methanosarcina and Methanoseta species predominate during anaerobic treatment of paper and pulp mill wastewater.

Table 4 Pearson correlations between VFA concentrations and DGGE band intensities during anaerobic digestion of M. sacchariflorus (R-MS)

In the present study, the sum of the concentrations of VFAs was 2150 mg L−1 (27.2 mM) in R-ZM and 3000 mg L−1 (37.5 mM) in R-MS. Ahring et al. [2] found that biogas production from manure was not inhibited by concentrations of individual VFAs up to 50 mM. Thus, any variation in biogas production in the present study was probably not due to excessive concentrations of VFAs, as also indicated by the positive correlations between the intensity of DGGE bands and the concentrations of individual VFAs.

Conclusions

The results of this study indicate that the content of lignocellulosic components in biomass influences the efficiency of biogas synthesis. Organics were removed with higher efficiency during the digestion of Z. mays silage, due to the higher content of non-fibrous carbohydrates in this substrate. In addition, non-fibrous materials were removed more efficiently from Z. mays than from M. sacchariflorus silages. The concentrations of individual VFAs differed slightly between the digesters, but propionate predominated in both. In both digesters, independently of the substrate used, acetoclastic Methanosaeta sp. and Methanosarcina sp. predominated, and hydrogenotrophic methanogens were also detected. The abundance of hydrogenotrophic methanogens strongly correlated with valerate concentration, whereas the abundance of acetoclastic methanogens correlated with isovalerate concentration. To sum up, this study indicates that the presence in a substrate of lignin, or metabolites of lignin degradation, or both, uncouples organics removal and methane production, and influences the structure of the microbial communities responsible for anaerobic digestion.