- Open Access
Lessons learned from the microbial ecology resulting from different inoculation strategies for biogas production from waste products of the bioethanol/sugar industry
Biotechnology for Biofuelsvolume 9, Article number: 144 (2016)
During strategic planning of a biogas plant, the local availability of resources for start-up and operation should be taken into consideration for a cost-efficient process. Because most bioethanol/sugar industries in Brazil are located in remote areas, the use of fresh cattle manure from local farms could be a solution for the inoculation of the biogas process. This study investigated the diversity and dynamics of bacterial and archaeal communities and the performance of biogas reactors inoculated with manure and a mixed inoculum from different biogas reactors as for a controlled start-up until steady state.
Laboratory-scale biogas reactors were fed semi-continuously with sugarcane filter cake alone (mono-digestion) or together with bagasse (co-digestion). At the initial start-up, the reactors inoculated with the mixed inoculum displayed a less diverse taxonomic composition, but with higher presence of significant abundances compared to reactors inoculated with manure. However, in the final steady state, the communities of the differently inoculated reactors were very similarly characterized by predominance of the methanogenic genera Methanosarcina and Methanobacterium, the bacterial families Bacteroidaceae, Prevotellaceae and Porphyromonadaceae (phylum Bacteroidetes) and Synergistaceae (phylum Synergistetes). In the mono-digestion reactors, the methanogenic communities varied greater than in the co-digestion reactors independently of the inoculation strategy.
The microbial communities involved in the biogas production from waste products of the Brazilian bioethanol/sugar industry were relatively similar and stable at the reactor’s steady phase independently of the inoculum source (manure or mixed inoculum). Therefore, the locally available manure can be used as inoculum for start-up of the biogas process, since it also contains the microbial resources needed. The strong fluctuation of methanogenic communities in mono-digestion reactors indicates higher risk of process instability than in co-digestion reactors.
The Brazilian bioethanol/sugar industry has been previously reported in our studies to have a big potential to improve the local bioeconomy while reducing greenhouse gas emission by applying biogas technology to the treatment of its waste products [1–3]. However, to make the biogas technology profitable and reliable for the Brazilian bioethanol/sugar industry, the reactor design and start-up requires strategic considerations.
In regions where biogas plants are widespread and well developed, a biogas reactor can be started with inocula from already established processes. Contrary to this scenario, in Brazil there are only very few plants applying the anaerobic digestion (AD) process to treat waste at large scale and these plants are spread across the country and separated by long distances. Therefore, we have addressed the possibility of reactor inoculation with fresh cattle manure (FCM) as a locally available, potential inoculum. For comparison in terms of microbial robustness, we have also prepared an engineered mixed inoculum (referred to as MIX) originating from the digestate of different biogas reactors fed mainly with energy crops and agricultural wastes such as maize silage, thin stillage, straw and chicken manure. The first results regarding biogas production and feasibility of these two different inoculation strategies were recently reported by Janke et al. . Nevertheless, the microbiological background of this experiment remained to be investigated.
In the biogas process, a complex metabolic network of microorganisms is responsible for organic matter degradation, which proceeds in four steps: hydrolysis, acidogenesis, acetogenesis and methanogenesis. While bacteria are involved in the first three steps, methanogenic archaea are responsible for the last step. Methanogens are very sensitive to process changes due to the relative lack of functional redundancy and low diversity . Besides that, the methanogens are very important for the AD process stability, because they are directly involved in the removal of fermentation product acetate or makes the syntrophic oxidation of acetate and other fermentation products thermodynamically feasible by keeping hydrogen partial pressure low . Therefore, the methanogenic community requires particular regard for the development of an efficient and robust AD process.
The methane formation by the methanogens is carried out either via direct acetate conversion (aceticlastic methanogenesis) or reduction of CO2 with H2 (hydrogenotrophic methanogenesis). The determination of different methanogenic pathways has been described to have crucial implications for the design and operation of biogas reactors [5, 7–10], since the aceticlastic and hydrogenotrophic methanogens may have different growth rates depending on the reactor’s conditions . However, the factors controlling the balance of methanogenic pathways is still not clear and seems to depend on the substrate used and process conditions such as organic loading rate, reactor type and temperature [12–14].
Knowledge of microbial adaptation to environmental conditions is also very important to understand the complex interplay between bacteria and methanogens [15, 16], particularly during acclimatization leading to successful and efficient reactor operation . Moreover, correlations of methanogenic community data with process parameters have been reported to contain decisive information about shifting pathway dominance during biogas production [9, 18, 19].
In a comprehensive study, we investigated the impact of different inoculation strategies and digestion setups (mono-digestion of filter cake and its co-digestion with bagasse) on the microbial community composition and dynamics along the operation of six mesophilic laboratory-scale continuously stirred tank reactors. In addition, the bacterial and methanogenic communities were assessed to identify key microorganisms of indicator value for the AD of the waste products from the bioethanol/sugar industry. Moreover, this study evaluated the links between microbial composition and reactor performance as a basis for an efficient future microbial resource management. The microbial communities were analyzed by DNA-based fingerprinting techniques and next-generation sequencing. The correlation between the microbial communities and reactor parameters was established by multivariate data analysis. In addition, the dominant methanogenic pathway was determined by the stable isotope fingerprinting of the produced biogas, following approaches applied to various anaerobic systems  and biogas reactors with several substrates [7, 21, 22].
Six identical laboratory-scale continuously stirred tank reactors with working volumes of 3 L under mesophilic conditions (38 ± 1 °C) were established to operate three experiments in duplicate: R3.3 and R3.4 were inoculated with MIX and fed with a single substrate (mono-digestion of filter cake); R3.5 and R3.6 obtained the same inoculum, but the filter cake was co-digested with bagasse; R3.7 and R3.8 were inoculated with FCM and performed co-digestion similarly to R3.5 and R.3.6. Figure 1 shows the major process parameters for each reactor along the whole experiment. The analytical methods and detailed description of the start-up and performance of the reactors are reported in our previous study . Briefly, the experiments of 137 days duration were divided into two phases: start-up (until day 69) and steady state (from day 70 on), as shown in Table 1. In our experiments, the day 0 does not represent the inoculation with MIX and FCM, but the end of an acclimatization time of around 24 h. On day 0, the feeding of the reactors was started. During the start-up phase, the organic loading rate (OLR) varied from 1.0 (only R3.7 and R3.8) or 2.0 (the other four reactors) to 2.5 gvs L−1 days−1 (all reactors), whereas during the steady state the OLR was constant at 3.0 gvs L−1 days−1. The hydraulic retention times (HRT) at the final stage of the steady state were 23 and 28 days for mono- and co-digestion, respectively.
Microbial community analysis
Duplicate digestate samples from each reactor were taken for molecular analysis on specific days and stored at −20 °C until further analysis. Simultaneously, biogas samples were taken from the reactors’ headspace for stable isotope analysis. The total genomic DNA of the bacterial and methanogenic communities was extracted with the ‘NucleoSpin Soil’ kit (Macherey–Nagel) as recommended by the supplier. The buffers SL2 and SX were used.
PCR amplifications for terminal restriction fragment length polymorphism (T-RFLP) screening of the methanogenic community were targeting the mcrA genes using the forward primer mlas and the reverse primer mcrA-rev and following the PCR protocol of Steinberg and Regan . T-RFLP analysis of purified PCR products was conducted after digestion with the restriction enzyme BstNI using the fragment size standard GeneScan-500 ROX (Applied Biosystems GmbH, Weiterstadt, Germany). T-RFLP electropherograms were processed as described by Lucas et al. . During statistical analysis in R, signals with low peak areas were removed using a cutoff of 12 times the standard deviation of the data sets. The reproducibility of the T-RFLP was validated by comparing the results with duplicate samples from each reactor at a particular sampling day (Additional file 1: Figure S1). The mcrA-derived T-RFs were assigned taxonomically using cloned mcrA amplicons database from anaerobic digester sample analyses performed in our laboratory [21, 24–27].
The bacterial community analysis was performed only for the four co-digestion reactors inoculated with either FCM or MIX. Samples from two times (days 0 and 44) in the start-up phase and from one time (day 113) of the steady state were processed based on the 16S ribosomal RNA genes and further analyzed on the 454-pyrosequencing platform GS Junior (Roche) as described by Ziganshin et al. . The variable regions V1–V3 of the bacterial 16S rRNA gene fragments were amplified with the primers Bac27F (5′-AGAGTTTGATCMTGGCTCAG-3′) and Bac519R (5′-GWATTACCGCGGCKGCTG-3′) using the Phire Hot Start II DNA Polymerase (Thermo Scientific). The raw sequence data were assessed with the QIIME 1.8.0 Virtual Box release . Further data processing was performed according to Lucas et al.  and Sun et al. . In summary, the dataset was firstly quality filtered by excluding sequences that were shorter than 150 and longer than 590 bp in lengths, comprised an average quality score below 25, held 50 bp at the end section below the quality score threshold of 25, comprised ambiguous bases, held a homopolymer run with more than 6 bp, or did not comprise any primer or barcode sequence. The USEARCH pipeline was applied on the sequences for further quality filtering based on non-chimeric sequences and for clustering into operational taxonomic units (OTUs) consisting of 97 % identity threshold . The taxonomic classification based on representative sequences was performed using the Greengenes core set (gg_13_8)  and the Ribosomal Database Project classifier 2.2 . For the taxonomic alignment, the Infernal algorithm with default setting was used . Finally, the summarized OTU tables were constructed according to their taxonomy and abundance. Further, the visualization of the OTU tables was processed via the spreadsheet program. De-multiplexed sequences of the 12 samples were deposited under the EMBL-EBI accession number PRJEB12073 (http://www.ebi.ac.uk/ena/data/view/PRJEB12073). The ecological data analyses leading to chao1, Shannon and Simpson indices [34, 35] and rarefaction curve were also performed with the QIIME software based on alpha diversity. Due to the differences in sequencing library size between the samples, we have used QIIME further to subsample (rarefy) the libraries down to 8000 sequences per sample for comparative diversity analyses and to calculate the beta diversity (pairwise sample dissimilarity).
The ordination of the dissimilarity matrices achieved by non-metric multidimensional scaling (NMDS) was processed as reported by Lucas et al. . Shortly, the variability of the microbial communities was evaluated by the Bray–Curtis dissimilarity index based on the presence and relative abundance. Thus, e.g., highly similar community composition is indicated by tiny distances. The correlation of reactor parameters and microbial communities based on the relative abundance was analyzed with the ‘envfit’ function and its significance was tested by a Monte Carlo test with 999 permutations. The significance threshold was set to a maxima of 0.001 and 0.05 for the methanogenic and the bacterial communities, respectively.
Stable isotope fingerprinting
The carbon and hydrogen stable isotope compositions of CH4 and CO2 from each reactor were measured in triplicate biogas samples collected in 20-mL gas-tight pre-evacuated vials. For analysis, an isotope ratio mass spectrometry system (Finnigan MAT 253, Thermofinnigan Bremen) coupled to a gas chromatograph (GC) (HP 6890 Series, Agilent Technology, USA) via a combustion device and a pyrolysis unit (with a water-removal assembly) was used for carbon and hydrogen measurements, respectively. Fifty-microliter biogas sample was injected into the inlet tube of the GC instrument equipped with a CP-Porabond Q column (50 m × 0.32 mm ID Varian, USA) held at a constant temperature of 40 °C. Helium was used as a carrying gas at the split ratio of 1:50 for carbon and 1:5 for hydrogen analysis. The isotope ratios of all samples are given in delta notation (δ13C and δ2H) in per mil (‰) units according to the standards VPDB (Vienna Pee Dee Belemnite) for carbon and VSMOW (Vienna Standard Mean Ocean Water) for hydrogen.
Results and discussion
Bacterial community succession
The microbial profiles of anaerobic digesters have been reported to be very specific for each type of reactor and substrate feeding [12, 28, 36]. Thus, to investigate the shaping forces of novel substrates (filter cake and bagasse) on the inocula (FCM and MIX), the bacterial community in the co-digestion reactors (R3.5, R3.6, R3.7 and R3.8) was assessed by amplicon pyrosequencing at three sampling times, i.e., days 0 and 44 in the start-up phase and day 113 during steady state. The bacterial community succession of each reactor is shown on phylum level in Fig. 2. The parallel reactors R3.5 and R3.6, and R3.7 and R3.8 had very similar bacterial profiles. The community similarity between the reactor samples based on the beta diversity are shown in Additional file 1: Table S1 and Figure S3. All OTUs obtained from the 12 samples are presented in Additional file 2: Table S2, and the respective rarefaction curves are shown in the Additional file 1: Figure S2.
We obtained on average around 12,000 high-quality sequence reads per reactor, varying from circa 9000 to 13,000. Along the operation of the co-digestion reactors during start-up and steady state, the taxonomic composition of our analyzed reactors was very diverse with a total of 1137 OTUs, but only 18 core OTUs (1.6 %), which were found in all reactors at all sampling times (Fig. 3), accounting for 18 % of all sequence reads. Thus, the variation of bacterial communities was very high with only a few microorganisms being significantly abundant in all reactors throughout the experiment. Examples were the families Porphyromonadaceae and Synergistaceae accounting for 8 and 4 % of all sequence reads, respectively.
Figure 3 also shows the variation of the numbers of core OTUs and phyla along the experiment. Day 0 samples contained the highest numbers of unique OTUs shared either only by the MIX parallels R3.5 and R3.6 (242, 21 % of all OTUs) or by the FCM parallels R3.7 and R3.8 (313, 28 % of all OTUs). The parallel FCM reactors R3.7 and R3.8 were more diverse in terms of OTUs, but only few microorganisms of high relative abundance were detected. On the other hand, parallel MIX reactors R3.5 and R3.6 presented lower diversity comprising microorganisms with high relative abundances (Fig. 2). On day 44, the number of total shared OTUs and phyla for all reactors was significantly increased, whereas the amount of unique OTUs per parallel reactors was drastically decreased (57 for R3.5 and R3.6 and 90 for R3.7 and R3.8). On day 113, the number of total OTUs and phyla shared by all reactors was decreased to 110 and 16, respectively. Moreover, the number of unique OTUs shared by parallel reactors was much lower.
These variations indicate that despite different inoculation, the bacterial communities were already very similar after 44 days of reactor operation. Nevertheless, the consolidation and stabilization of the bacterial communities proceeded until the steady state. Figure 3 shows that there was uniqueness in terms of OTUs and phyla of individual reactors along the experiment, although R3.5 and R3.6 as well as R3.7 and R3.8 were operated in parallel under the same conditions.
The estimated richness of the bacterial community in all samples is shown in Table 2. The number of OTUs decreased slightly during start-up, whereas at steady state a drastic drop was observed. The same trend was reflected by the OTU richness estimator chao 1 and the Shannon index. The Simpson diversity was stable between the two samplings during start-up, while it decreased during steady state. In MIX reactors R3.5 and R3.6, the Simpson index dropped more than in the FCM reactors R3.7 and R3.8, indicating higher community evenness in the former. Generally, the bacterial species richness was lower during steady state with more pronounced predominance of some microorganisms.
The phyla Firmicutes and Bacteroidetes dominated the bacterial communities in all analyzed reactors along the entire operation time (Fig. 2). Both together comprised 72 % of all sequence reads. The presence of most OTUs affiliated with these two phyla known to utilize carbohydrates has frequently been reported in studies about AD of different substrates, mainly maize silage and manure, in laboratory- and full-scale biogas reactors [12, 24, 28, 36, 37].
Although the succession toward steady state was not favorable for Firmicutes, the phylum was the most diverse taxonomic group in the reactors with 620 OTUs (54 %), represented mainly by the order Clostridiales (485 OTUs), in which the families Ruminococcaceae (153 OTUs), Clostridiaceae (53 OTUs) and Lachnospiraceae (41 OTUs) were predominant. Following the succession of these families during our experiment, Ruminococcaceae was the most constant family, without much variation in its abundance (around 5 % in all reactor samples). Clostridiaceae had also a constant abundance around 0.8 % in the reactors inoculated with MIX. However, the same family was the most variable one in the reactors inoculated with FCM, presenting 21.8 and 1.3 % relative abundances at the beginning (day 0) and end of the experiment (day 113), respectively. The family Lachnospiraceae was relatively constant at abundances around 1.8 %, except for the reactors at the beginning of the experiment with MIX. Hence, our results with sugarcane waste products were in agreement with former studies that also investigated the bacterial communities in the AD of plant-based biomass [37, 38], especially in terms of the prevalence of Clostridiales. Members of this order are equipped with the cellulosome, a multienzyme complex, which enables them to efficiently hydrolyze recalcitrant cellulosic and hemicellulosic structures in the plant cell wall . Our feeding substrates from the sugarcane industry contain high percentages of cellulose and hemicellulose, both together representing 75 and 55 % of the total carbohydrate and lignin content of bagasse and filter cake, respectively . Sequences affiliated to the cellulosome-producing bacterium Ruminococcus flavefaciens were detected in our reactors. These type of sequences related to the anaerobic cellulolytic rumen bacterium contributed significant proportions of sequence reads in our experiment (0.79 % of all), especially at the end at day 113 (0.34 %), which suggests that this phylotype was one of the specialists degrading bagasse and filter cake, since both inocula sustained its presence. Members of the genus Ruminococcus are also known to produce hydrogen and acetate, thus supplying the hydrogenotrophic methanogenic genus Methanobacterium .
Bacteroidetes as the second most diverse phylum with 128 OTUs showed a gradual increase in relative abundance toward steady state. Within this phylum, OTUs affiliated with the genera Bacteroides (4 OTUs, ca. 9.4 % of all sequence reads) and Prevotella (4 OTUs, ca. 7.2 % of all reads) and the family Porphyromonadaceae (39 OTUs, ca. 9.8 % of all reads) were predominant. Since most of these sequence reads were detected at the end of the experiment, we assume that these microorganisms are crucial for the hydrolysis and fermentation of the filter cake and bagasse. The genus Bacteroides represented by the species B. cellulosolvens  is also known for its cellulosome. Within the genus Prevotella, there are some species notably involved in the degradation of hemicellulose, e.g., Prevotella paludivivens .
On day 0, the phyla Firmicutes and Bacteroidetes together represented about 85 and 62 % abundance in the reactors inoculated with FCM and MIX, respectively (Fig. 2). In R3.5 and R3.6 (MIX), other phyla such as Actinobacteria, Fibrobacteres, Spirochaetes, Synergistetes, Tenericutes, Thermotogae and the candidate WWE1 were also presented in significant relative abundances. Most of these phyla are involved in the degradation of lignocellulose-rich substrates [14, 28, 36]. However, the ecophysiological role of some phyla such as WWE1 is still unclear .
On day 44, the bacterial communities were already very similar in all reactors independent of the inoculation source. The relative abundance of the phylum Firmicutes decreased, whereas that of the Bacteroidetes increased. While the relative abundance of Synergistetes was kept constant in R3.5 and R3.6 (MIX), it increased slightly in the other co-digestion reactors (R3.7 and R3.8). The third most abundant phylum at this time was the candidate WWE1 (12 %), which had been found in many other studies dealing with mesophilic AD in the frame of wastewater treatment [43, 44]. Moreover, this candidate phylum was recently reported in large-scale continuously stirred tank reactors digesting the cellulose-rich substrate maize silage .
The steady state on day 113 was characterized by a less diverse community with three major phyla, Bacteroidetes, Firmicutes and Synergistetes, comprising nearly 94 % of the entire community. The relative abundance of the phylum Bacteroidetes of approximately 66 % in all reactors was even greater than on day 44. On the contrary, the abundance of Firmicutes dropped further to about 16 %. Synergistetes were present in all reactors along the experiment and their increasing abundance (to ca. 11 %) demonstrated their participation in the degradation of the lignocellulosic bagasse and filter cake. Within this phylum, the class Synergistia has reported to degrade fiber-rich feedstock . Whereas all reactors had a very similar overall community structure, the phylum Chloroflexi (specifically, class Anaerolineae) was significantly abundant (3 %) only in R3.5 and R3.6 (MIX). This phylum had also been found in similar proportion in a biogas plant co-digesting maize silage, green rye and liquid manure .
The plot in Fig. 4 shows an NMDS analysis of bacterial 16S rRNA genes. The bacterial community compositions notably converged toward steady state (day 113). It is also visible that duplicate FCM reactors presented more diverse bacterial communities than MIX reactors. The correlation between the bacterial communities and the reactor parameters is indicated by the most significant vectors represented as arrows in the NMDS plot. The abundance of the phylum Synergistetes was strongly correlated with increasing OLR, whereas the abundance of the candidate phylum TM7 grew with the decrease of NH4–N. The more enriched values for δ13C and δ2H of methane were inversely correlated with NH4–N and positively with the bacterial community presented after the second sampling point (day 44). Furthermore, the decrease of pH influenced the bacterial communities primarily in the start-up phase. As already mentioned, the presence of phylum Bacteroidetes was strongly correlated with the steady state (at day 113).
Dynamics of the methanogenic community
The phylogenetic composition and dynamics of the methanogenic community were analyzed in all six reactors by mcrA-based T-RFLP fingerprinting. A heatmap presenting the relative abundance of the methanogenic communities along the whole experiment is shown in Fig. 5. Duplicate reactors (R3.3 and R3.4, R3.5 and R3.6, R3.7 and R3.8) had very similar methanogenic community compositions at each sampling time (Additional file 1: Figure S4). However, some non-abundant microorganisms, for examples Methanomassiliicoccaceae at the initial phase and Methanomassiliicoccus and Methanoregulaceae at the steady phase were occasionally detected in one of the duplicate reactors. Except during start-up, microorganisms affiliated with the genus Methanosarcina were the most abundant methanogens across all analyzed samples, reaching as much as 90 % T-RF abundance at day 51 in the final stage of the start-up phase.
The first sampling on day 0 was characterized by different methanogens in the MIX and FCM inocula. MIX-inoculated reactors were dominated by the genera Methanoculleus (ca. 50 %) and Methanosaeta (ca. 30 %) and less by Methanomassiliicoccus (ca. 8 %). Methanosarcina and Methanomassiliicoccaceae were not detected simultaneously in all four reactors at day 0. In the FCM-inoculated reactors, the most abundant genera were Methanospirillum (ca. 50 %), Methanosarcina (ca. 30 %) and Methanobacterium (ca. 20 %). In case of both inoculation strategies, the predominant methanogens were affiliated to either strictly hydrogenotrophic taxa or to the metabolically flexible Methanosarcina genus. Nevertheless, strict aceticlastic methanogens were also present at usually minor abundances.
The addition of bagasse as a co-substrate in MIX-inoculated reactors had an influence on the methanogenic community. The T-RF profiles of the mono- and co-digestion on day 28 were already significantly distinct. On the other hand, the methanogenic community structures of the co-digestion reactors (MIX and FCM) were similar after sampling day 21. In this latter case, the dynamics of the microbial community also followed the same trend toward the predominance of the Methanosarcina.
The genus Methanosarcina has been described to out-compete Methanosaeta at elevated acetate levels, due to its faster growth kinetics [45, 46]. However, recently Methanosaeta was reported to dominate over Methanosarcina even at high acetate concentration in the AD of dairy waste with poultry waste as co-substrate . In accordance with this recent finding, the high abundance of the strictly aceticlastic Methanosaeta was stable during volatile organic acids (VOA) accumulation (1273 mg L−1 of acetate and 55 mg L−1 of propionate) at day 35. This corroborates our previous study , in which Methanosaeta was also detected at high VOA concentrations.
In the final part of the start-up phase (days 42–69), Methanosarcina was the most abundant methanogen reaching the highest T-RF relative abundance in all reactors. The presence of Methanosarcina in AD is regarded to be advantageous. As discussed by De Vrieze et al. , it is very robust against many stress conditions, e.g., tolerate process overloading, sudden pH changes and high levels of ammonium and salt. In addition, Methanosarcina has high growth rates, can use both the aceticlastic and the hydrogenotrophic pathway and utilize methanol to produce methane.
In the late start-up stage, Methanoculleus and Methanospirillum, which had been dominant in the MIX and FCM inocula, respectively, were no longer abundant. The taxon represented by T-RF 180 in the mono-digestion (MIX) and Methanoregulaceae in the co-digestion (MIX) became abundant until steady state.
The relative abundance of Methanosarcina decreased at steady state to around 30 and 50 % in the mono- and co-digestion reactors, respectively. In the mono-digestion, the taxon represented by T-RF 180 and Methanobacterium became more abundant, while Methanoregulaceae appeared first in the process. In the four co-digestion reactors, the methanogenic communities were very similar during steady state. Generally, the abundances of Methanomassiliicoccus and Methanobacterium increased, whereas Methanoregulaceae was not abundant.
Figure 6 shows the results of a multivariate statistical analysis based on the correlation between T-RLFP data and reactor parameters such as biogas yield and composition (CH4 and CO2), volatile fatty acids (acetate, propionate and butyrate), pH, total VOA, volatile organic acids per total inorganic carbonate buffer (VOA/TIC), HRT, ammonium–nitrogen (NH4–N) and stable isotope data (δ13CCH4, δ13CCO2 and δ2HCH4). The effect of the inoculation strategy (MIX vs. FCM) using the co-digestion setup and the effect of the substrate (mono vs. co-digestion) applying MIX inoculum were analyzed separately for more comprehensible results.
Figure 6a compares the different inoculation strategies (MIX and FCM) in the co-digestion setup. The main difference between them was observed at the beginning of the start-up phase. After day 21, both setups, independent of the inoculum, had attained similar methanogenic community patterns with minor temporal compositional shifts. Methanoculleus and Methanosaeta were strongly correlated with the MIX inoculum at the beginning of the experiment, whereas Methanosarcina was predominantly present in the second part of the start-up phase. The methanogenic community dynamics was most significantly correlated (p < 0.001) with the pH, which gradually decreased during the start-up phase and kept constant values during steady state. On day 128, the T-RFLP profiles clustered together, indicating very high similarity of the methanogenic community structures, independent of the inoculum source.
Mono- and co-digestion with the same inoculum (MIX) were compared in an independent analysis (Fig. 6b). Although the substrate regimes differed, methanogenic community patterns were similar. Bigger changes in the structure of the methanogenic community were observed during mono-digestion, visible as a more extended dashed hull, while the samples from the co-digestion reactor clustered more closely together. The reactor parameters with most significant correlation (p < 0.001) to the community structures of the investigated samples were OLR, δ13CCH4, NH4–N and pH. Whereas OLR and δ13CCH4 values rose along the reactor operation, NH4–N concentration and pH value decreased gradually, mainly in the start-up phase.
Stable isotope characterization of the produced biogas
Isotopic variations as a function of inoculation strategies (MIX and FCM) and substrate regimes are shown in Fig. 7. Samples for δ13CCH4, δ13CCO2 and δ2HCH4 measurements were taken together with samples for the analyses of microbial profiles and reactor parameters. Substantial isotopic changes were mainly observed in the initial start-up phase (days 0–41). Thereafter, the isotope values followed very similar trends independent of the inoculum. Initially, the isotope values were strongly depleted, particularly in reactors inoculated with FCM (−64 ‰ for δ13CCH4, −6 ‰ for δ13CCO2 and −352 ‰ for δ2HCH4).
Figure 7a shows that the δ13CCH4 values during the initial start-up increased until day 44. As an exception on day 9 in the FCM-inoculated reactors, the δ13CCH4 values changed drastically to −25 ‰, but from day 21 onward the δ13CCH4 values followed the same trend as in the other reactors inoculated with MIX. A similar isotope signature for δ13CCH4 of about −32 ‰ on day 44 had also been observed in our previous study of co-digestion of filter cake and bagasse at gradual OLR increase  and in some other studies with C4 plant biomass, namely maize silage [7, 21]. Afterward, all reactors displayed constant isotope signatures also indicating stable methanogenic process.
The effect of different inocula and substrates on the δ13CCO2 composition is clearly observed in Fig. 7b. The reactors inoculated with MIX had very similar isotope values on days 0 and 9. But after that, the isotopic pattern diverged between mono- and co-digestion. Enriched δ13CCO2 values were found for mono-digestion (up to about 15 ‰). The reactors inoculated with FCM had much more depleted δ13CCO2 values on days 0 and 9. On day 79 during steady state, the δ13CCO2 values for all co-digestion reactors were very similar (around 11 ‰).
The δ2HCH4 values shown in Fig. 7c presented a similar tendency among the reactors at specific sampling times. After start feeding the freshly inoculated reactors (from day 9 on), the δ2HCH4 values indicated enrichment of methane in deuterium. Similarly to the dynamics of δ13CCO2 values, the influence of the co-digestion was also visible from the variation of the δ2HCH4 isotopic signatures, mainly at the last samplings. In co-digestion reactors, lower δ2HCH4 values were observed.
The δ13C composition of filter cake and bagasse were −14.30 and −13.64 ‰, respectively . The isotopic signature of the substrates has also influenced the trends observed in Fig. 7, in which the δ13C values became more enriched. The daily added tap water affects as well the δ2HCH4 isotopic signatures. Nikolausz et al.  reported that the δ2H of the tap water in Leipzig corresponds to about −64.7 ± 0.8 ‰. Therefore, more enriched δ2HCH4 values were found after the addition of water to the substrates.
The dominance of a specific methanogenic pathway along the reactor operation was assessed using the apparent fractionation factor (αC) (Fig. 8a), which was calculated based on the isotopic composition (δ13C) of methane and carbon dioxide. As previously reported [20, 49, 50], a predominance of hydrogenotrophic methanogenesis is indicated by αC >1.065, whereas αC <1.025 indicates the predominance of aceticlastic methanogenesis. Intermediate αC values mostly found along the experiment showed that both methanogenic pathways were involved in methane production. This is also in agreement with the methanogenic community profiles composed of hydrogenotrophic and aceticlastic methanogens (Fig. 5).
Only on days 0 and 9 in co-digestion reactors inoculated with FCM, the predominance of one methanogenic pathway was observed. On day 0, the αC value indicated the predominance of hydrogenotrophic methanogenesis in accordance with a 68 % predominance of the strictly hydrogenotrophic genera Methanospirillum (ca. 50 %) and Methanobacterium (18 %). In addition, the versatile genus Methanosarcina, which contributed 32 % to the community profile, most probably produced methane via the hydrogenotrophic pathway during this period. On day 9, the methanogenic pathway was shifted drastically to the dominance of aceticlastic methanogenesis. The αC value was slightly lower than 1.025 and the predominance of the versatile Methanosarcina (ca. 75 %) suggested that members of this genus became not just numerically abundant, but switched to methane production from acetate. On day 21, again a sudden change was observed from the αC dynamics. At this time, Methanosarcina was still the most abundant taxon, which, however, probably used both methanogenic pathways.
The methanogenic pathway dynamics during the experiment is possibly also observed in plots combining δ13C and δ2H isotopic variations of methane during start-up (Fig. 8b) and steady state (Fig. 8c). Nikolausz et al.  had used this combination of δ13CCH4 and δ2HCH4 to identify the dominant methanogenic pathway in laboratory-scale biogas reactors. Accordingly, the suggested predominance of the hydrogenotrophic methanogenesis at day 0 (marked in the graphic with the hull) for FCM-inoculated reactors is confirmed as suggested by Sugimoto and Wada . Although their classification did not correspond to the other sampling points in terms of predominant methanogenic pathway defined by αC and the molecular biological results, the differences between the inocula and the substrate were still observed in the combined δ13CCH4 and δ2HCH4 graphics (hulls). A general trend of less negative δ2HCH4 values observed in case of later samples in all reactors is probably due to the influence of the isotopic composition of the supplemented process water.
Our results confirmed that FCM is a reliable and efficient inoculum for the co-digestion of filter cake and bagasse, since very similar methane and biogas yields were obtained under the steady-state phase, independent of the inoculation strategy. The bacterial and methanogenic communities were also very similar at the end of the experiment, regardless of whether the reactors had been inoculated with FCM or MIX. Bacterial composition and succession showed that the major phyla involved in the anaerobic degradation of the waste products from the Brazilian bioethanol/sugar industry were Bacteroidetes, Firmicutes and Synergistetes. The co-digestion of filter cake and bagasse in both inocula setups led to the development of polysaccharide-degrading specialists affiliated to the genera Prevotella and Bacteroides (both comprised around 50 % relative abundance of the bacterial community at steady phase). Methanogenic communities varied mainly in the first 3 weeks of operation for co-digestion reactors and in the first 5 weeks for mono-digestion. The key methanogens were affiliated with Methanosarcina, Methanobacterium, Methanoregulaceae and Methanomassiliicoccus in co-digestion reactors, as opposed to Methanosarcina, Methanobacterium, Methanoregulaceae and a non-identified taxon in mono-digestion reactors. The most important reactor parameter correlating with the methanogenic community structure in both digestion setups was the pH. Stable isotope fingerprinting showed clearly, especially in case of carbon dioxide, the influence of different inoculum strategies and substrate feeding on the methanogenic activity. Based on the isotope analysis, in agreement with the molecular approach, both aceticlastic and hydrogenotrophic methanogenesis pathways contributed importantly to methane production throughout the experiment.
fresh cattle manure
hydraulic retention time
non-metric multidimensional scaling
engineered mixed inoculum
organic loading rate
operational taxonomic units
terminal restriction fragment length polymorphism
volatile organic acid
volatile organic acids per total inorganic carbonate buffer
Leite AF, Janke L, Harms H, Zang JW, Fonseca-Zang WA, Stinner W, Nikolausz M. Assessment of the variations in characteristics and methane potential of major waste products from the Brazilian bioethanol industry along an operating season. Energy Fuels. 2015;29(7):4022–9.
Janke L, Leite A, Nikolausz M, Schmidt T, Liebetrau J, Nelles M, Stinner W. Biogas production from sugarcane waste: assessment on kinetic challenges for process designing. Int J Mol Sci. 2015;16:20685–703.
Janke L, Leite A, Wedwitschka H, Schmidt T, Nikolausz M, Stinner W. Biomethane Production Integrated to the Brazilian Sugarcane Industry: The Case Study of São Paulo State. In: Proceedings of the 22nd European biomass conference and exhibition. 2014. p. 1295–9.
Janke L, Leite AF, Nikolausz M, Radetski CM, Nelles M, Stinner W. Comparison of start-up strategies and process performance during semi-continuous anaerobic digestion of sugarcane filter cake co-digested with bagasse. Waste Manag. 2016;48:199–208.
Demirel B. Major pathway of methane formation from energy crops in agricultural biogas digesters. Crit Rev Environ Sci Technol. 2014;44:199–222.
Oh ST, Martin AD. Long chain fatty acids degradation in anaerobic digester: thermodynamic equilibrium consideration. Process Biochem. 2010;45:335–45.
Nikolausz M, Walter RF, Strauber H, Liebetrau J, Schmidt T, Kleinsteuber S, Bratfisch F, Gunther U, Richnow HH. Evaluation of stable isotope fingerprinting techniques for the assessment of the predominant methanogenic pathways in anaerobic digesters. Appl Microbiol Biotechnol. 2013;97:2251–62.
Gehring T, Klang J, Niedermayr A, Berzio S, Immenhauser A, Klocke M, Wichern M, Lubken M. Determination of methanogenic pathways through carbon isotope (delta13C) analysis for the two-stage anaerobic digestion of high-solids substrates. Environ Sci Technol. 2015;49:4705–14.
Lv Z, Hu M, Harms H, Richnow HH, Liebetrau J, Nikolausz M. Stable isotope composition of biogas allows early warning of complete process failure as a result of ammonia inhibition in anaerobic digesters. Bioresour Technol. 2014;167:251–9.
Ferry JG. Fundamentals of methanogenic pathways that are key to the biomethanation of complex biomass. Curr Opin Biotechnol. 2011;22:351–7.
Zinder S. Physiological ecology of Methanogens. In: Ferry J, editor. Methanogenesis. Chapman & hall microbiology series. Berlin: Springer; 1993. p. 128–206.
Abendroth C, Vilanova C, Gunther T, Luschnig O, Porcar M. Eubacteria and archaea communities in seven mesophile anaerobic digester plants in Germany. Biotechnol Biofuels. 2015;8:87.
Nettmann E, Bergmann I, Pramschufer S, Mundt K, Plogsties V, Herrmann C, Klocke M. Polyphasic analyses of methanogenic archaeal communities in agricultural biogas plants. Appl Environ Microbiol. 2010;76:2540–8.
Sun L, Pope PB, Eijsink VG, Schnurer A. Characterization of microbial community structure during continuous anaerobic digestion of straw and cow manure. Microb Biotechnol. 2015;8:815–27.
Sun L, Müller B, Westerholm M, Schnürer A. Syntrophic acetate oxidation in industrial CSTR biogas digesters. J Biotechnol. 2014;171:39–44.
Batstone D, Picioreanu C, Van Loosdrecht M. Multidimensional modelling to investigate interspecies hydrogen transfer in anaerobic biofilms. Water Res. 2006;40:3099–108.
Westerholm M, Muller B, Isaksson S, Schnurer A. Trace element and temperature effects on microbial communities and links to biogas digester performance at high ammonia levels. Biotechnol Biofuels. 2015;8:154.
Polag D, Krapf LC, Heuwinkel H, Laukenmann S, Lelieveld J, Keppler F. Stable carbon isotopes of methane for real-time process monitoring in anaerobic digesters. Eng Life Sci. 2014;14:153–60.
Karakashev D, Batstone DJ, Angelidaki I. Influence of environmental conditions on methanogenic compositions in anaerobic biogas reactors. Appl Environ Microbiol. 2005;71:331–8.
Conrad R. Quantification of methanogenic pathways using stable carbon isotopic signatures: a review and a proposal. Org Geochem. 2005;36:739–52.
Lv Z, Leite AF, Harms H, Richnow HH, Liebetrau J, Nikolausz M. Influences of the substrate feeding regime on methanogenic activity in biogas reactors approached by molecular and stable isotope methods. Anaerobe. 2014;29:91–9.
Polag D, May T, Muller L, Konig H, Jacobi F, Laukenmann S, Keppler F. Online monitoring of stable carbon isotopes of methane in anaerobic digestion as a new tool for early warning of process instability. Bioresour Technol. 2015;197:161–70.
Steinberg LM, Regan JM. Phylogenetic comparison of the methanogenic communities from an acidic, oligotrophic fen and an anaerobic digester treating municipal wastewater sludge. Appl Environ Microbiol. 2008;74:6663–71.
Lucas R, Kuchenbuch A, Fetzer I, Harms H, Kleinsteuber S. Long-term monitoring reveals stable and remarkably similar microbial communities in parallel full-scale biogas reactors digesting energy crops. FEMS Microbiol Ecol. 2015. doi:10.1093/femsec/fiv004.
Popp D, Schrader S, Kleinsteuber S, Harms H, Strauber H. Biogas production from coumarin-rich plants-inhibition by coumarin and recovery by adaptation of the bacterial community. FEMS Microbiol Ecol. 2015;91:fiv103.
Leite A, Janke L, Lv Z, Harms H, Richnow H-H, Nikolausz M. Improved monitoring of semi-continuous anaerobic digestion of sugarcane waste: effects of increasing organic loading rate on methanogenic community dynamics. Int J Mol Sci. 2015;16:23210.
Buhligen F, Lucas R, Nikolausz M, Kleinsteuber S. A T-RFLP database for the rapid profiling of methanogenic communities in anaerobic digesters. Anaerobe. 2016;39:114–6.
Ziganshin AM, Liebetrau J, Proter J, Kleinsteuber S. Microbial community structure and dynamics during anaerobic digestion of various agricultural waste materials. Appl Microbiol Biotechnol. 2013;97:5161–74.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.
McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, Andersen GL, Knight R, Hugenholtz P. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6:610–8.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.
Nawrocki EP, Kolbe DL, Eddy SR. Infernal 1.0: inference of RNA alignments. Bioinformatics. 2009;25:1335–7.
Magurran AE. Measuring biological diversity. Hoboken: Wiley; 2004.
Hill TCJ, Walsh KA, Harris JA, Moffett BF. Using ecological diversity measures with bacterial communities. FEMS Microbiol Ecol. 2003;43:1–11.
Strauber H, Lucas R, Kleinsteuber S. Metabolic and microbial community dynamics during the anaerobic digestion of maize silage in a two-phase process. Appl Microbiol Biotechnol. 2016;100(1):479–91.
Krober M, Bekel T, Diaz NN, Goesmann A, Jaenicke S, Krause L, Miller D, Runte KJ, Viehover P, Puhler A, Schluter A. Phylogenetic characterization of a biogas plant microbial community integrating clone library 16S-rDNA sequences and metagenome sequence data obtained by 454-pyrosequencing. J Biotechnol. 2009;142:38–49.
Klocke M, Mahnert P, Mundt K, Souidi K, Linke B. Microbial community analysis of a biogas-producing completely stirred tank reactor fed continuously with fodder beet silage as mono-substrate. Syst Appl Microbiol. 2007;30:139–51.
Bayer EA, Belaich JP, Shoham Y, Lamed R. The cellulosomes: multienzyme machines for degradation of plant cell wall polysaccharides. Annu Rev Microbiol. 2004;58:521–54.
Latham M, Wolin M. Fermentation of cellulose by Ruminococcus flavefaciens in the presence and absence of Methanobacterium ruminantium. Appl Environ Microbiol. 1977;34:297–301.
Schwarz WH. The cellulosome and cellulose degradation by anaerobic bacteria. Appl Microbiol Biotechnol. 2001;56:634–49.
Ueki A, Akasaka H, Satoh A, Suzuki D, Ueki K. Prevotella paludivivens sp. nov., a novel strictly anaerobic, gram-negative, hemicellulose-decomposing bacterium isolated from plant residue and rice roots in irrigated rice-field soil. Int J Syst Evol Microbiol. 2007;57:1803–9.
Chouari R, Le Paslier D, Dauga C, Daegelen P, Weissenbach J, Sghir A. Novel major bacterial candidate division within a municipal anaerobic sludge digester. Appl Environ Microbiol. 2005;71:2145–53.
Pelletier E, Kreimeyer A, Bocs S, Rouy Z, Gyapay G, Chouari R, Riviere D, Ganesan A, Daegelen P, Sghir A, et al. “Candidatus Cloacamonas acidaminovorans”: genome sequence reconstruction provides a first glimpse of a new bacterial division. J Bacteriol. 2008;190:2572–9.
Blume F, Bergmann I, Nettmann E, Schelle H, Rehde G, Mundt K, Klocke M. Methanogenic population dynamics during semi-continuous biogas fermentation and acidification by overloading. J Appl Microbiol. 2010;109:441–50.
Griffin ME, McMahon KD, Mackie RI, Raskin L. Methanogenic population dynamics during start-up of anaerobic digesters treating municipal solid waste and biosolids. Biotechnol Bioeng. 1998;57:342–55.
Chen S, He Q. Persistence of Methanosaeta populations in anaerobic digestion during process instability. J Ind Microbiol Biotechnol. 2015;42(8):1129–37.
De Vrieze J, Hennebel T, Boon N, Verstraete W. Methanosarcina: the rediscovered methanogen for heavy duty biomethanation. Bioresour Technol. 2012;112:1–9.
Whiticar MJ, Faber E, Schoell M. Biogenic methane formation in marine and freshwater environments: CO2 reduction vs. acetate fermentation—isotope evidence. Geochim Cosmochim Acta. 1986;50:693–709.
Galand PE, Yrjälä K, Conrad R. Stable carbon isotope fractionation during methanogenesis in three boreal peatland ecosystems. Biogeosciences. 2010;7:3893–900.
Sugimoto A, Wada E. Hydrogen isotopic composition of bacterial methane: CO2/H2 reduction and acetate fermentation. Geochim Cosmochim Acta. 1995;59:1329–37.
Oliveros JC. VENNY. An interactive tool for comparing lists with Venn diagrams. 2007. http://bioinfogp.cnb.csic.es/tools/venny/index.html. Accessed Nov 2015.
AFL designed the experiment, performed the molecular and isotope analyses, evaluated the data and prepared the manuscript. LJ designed the experiment, contributed to the data interpretation and reviewed the manuscript. HH supervised the study and thoroughly revised the manuscript. HHR contributed to the evaluation and interpretation of the isotope data and revised the manuscript. MN designed the experiment, contributed to the data interpretation, supervised the study and thoroughly revised the manuscript. All authors read and approved the final manuscript.
We thank the funding agency of the Brazilian scientific mobility program, Science without Borders (Pt.: Ciência sem Fronteiras), for the scholarship Grants of Athaydes Francisco Leite (202024/2012-1) and Leandro Janke (237938/2012-0). We gratefully acknowledge our collaboration partners from the Department Biochemical Conversion of the Deutsches Biomasseforschungszentrum (DBFZ) for the contribution to the analytical measurements, especially Bärbel Haase for the technical support. We very much appreciate the endorsement of Dr. Sabine Kleinsteuber during pyrosequencing analysis. We also would like to thank Ute Lohse for the assistance with the molecular analyses and Ursula Günther for the support with the isotope measurements. Furthermore, we thank Mustafa Kaya and Ellen Guimarães for their contribution to the T-RFLP analysis and to the isotope data evaluation, respectively. We are grateful to the i-NoPa project “Sustainable bioeconomy in Brazil: Bioenergy from biogas using various types of waste substrates from the Brazilian bioethanol industry”, which partially financed the research.
The authors declare that they have no competing interests.
Availability of supporting data
De-multiplexed sequences from the pyrosequencing analysis were deposited under the EMBL-EBI accession number PRJEB12073 (http://www.ebi.ac.uk/ena/data/view/PRJEB12073).
Further supporting data related to this study were included as additional supporting materials.