Skip to main content
  • Research article
  • Open access
  • Published:

Comparative metagenomics of biogas-producing microbial communities from production-scale biogas plants operating under wet or dry fermentation conditions

Abstract

Background

Decomposition of biomass for biogas production can be practiced under wet and dry fermentation conditions. In contrast to the dry fermentation technology, wet fermentation is characterized by a high liquid content and a relatively low total solid content. In this study, the composition and functional potential of a biogas-producing microbial community in an agricultural biogas reactor operating under wet fermentation conditions was analyzed by a metagenomic approach applying 454-pyrosequencing. The obtained metagenomic dataset and corresponding 16S rRNA gene amplicon sequences were compared to the previously sequenced comparable metagenome from a dry fermentation process, meeting explicitly identical boundary conditions regarding sample and community DNA preparation, sequencing technology, processing of sequence reads and data analyses by bioinformatics tools.

Results

High-throughput metagenome sequencing of community DNA from the wet fermentation process applying the pyrosequencing approach resulted in 1,532,780 reads, with an average read length of 397 bp, accounting for approximately 594 million bases of sequence information in total. Taxonomic comparison of the communities from wet and dry fermentation revealed similar microbial profiles with Bacteria being the predominant superkingdom, while the superkingdom Archaea was less abundant. In both biogas plants, the bacterial phyla Firmicutes, Bacteroidetes, Spirochaetes and Proteobacteria were identified with descending frequencies. Within the archaeal superkingdom, the phylum Euryarchaeota was most abundant with the dominant class Methanomicrobia. Functional profiles of the communities revealed that environmental gene tags representing methanogenesis enzymes were present in both biogas plants in comparable frequencies. 16S rRNA gene amplicon high-throughput sequencing disclosed differences in the sub-communities comprising methanogenic Archaea between both processes. Fragment recruitments of metagenomic reads to the reference genome of the archaeon Methanoculleus bourgensis MS2T revealed that dominant methanogens within the dry fermentation process were highly related to the reference.

Conclusions

Although process parameters, substrates and technology differ between the wet and dry biogas fermentations analyzed in this study, community profiles are very similar at least at higher taxonomic ranks, illustrating that core community taxa perform key functions in biomass decomposition and methane synthesis. Regarding methanogenesis, Archaea highly related to the type strain M. bourgensis MS2T dominate the dry fermentation process, suggesting the adaptation of members belonging to this species to specific fermentation process parameters.

Background

Rising energy costs and considerations on long term environmental sustainability have placed renewable energy sources in the focus of debate. The development of renewable energy resources offers the chance to replace traditional fossil fuels and can help to reduce carbon dioxide emissions [1,2]. An economically attractive technology to generate bioenergy is the production of biogas that is a mixture of methane (CH4) and carbon dioxide (CO2) as the main components, with small amounts of hydrogen sulfide (H2S), nitrogen (N2), hydrogen (H2), ammonia (NH3) and carbon monoxide (CO) [3]. The most common and widespread utilization of biogas is the production of electricity and heat by its combustion in combined heat and power units.

The process of biogas production takes place under anaerobic conditions and involves microbial decomposition of organic matter, yielding methane as the main final product of underlying metabolic pathways. In Germany, mostly maize silage combined with liquid manure is utilized as the substrate for biogas production [1,4]. Complex consortia of microorganisms are responsible for biomass decomposition and biogas production involving the stages substrate hydrolysis, acidogenesis, acetogenesis and methanogenesis. However, most of these microbes, as well as their roles in biogas production, are currently unknown. Recently, the analysis of the structure, composition and activity of microbial communities in relation to input substrates and fermentation parameters in biogas plants have become the focus of research [5-7]. It is generally accepted that a better understanding of the composition and activity of the multifarious microbial community is crucial for further optimization of reactor performance and fermentation process technologies. Moreover, to increase the yield of biogas, a detailed insight into relevant microbial metabolic pathways involved in methane synthesis and syntrophy is necessary.

Previous studies analyzed the taxonomic structure and enzymatic potential of biogas communities residing in agricultural biogas reactors. A classical microbiological approach for the analysis of the communities’ taxonomy is the generation of 16S rRNA gene clone libraries, followed by Sanger sequencing of the 16S rRNA gene fragments [8-12]. Sequencing of 16S rRNA gene clone libraries is limited since coverage of the microbial complexity frequently is laborious, costly and time-consuming. Moreover, sequence information on community 16S rRNA marker genes does not provide direct insights into functions of microorganisms. To achieve deeper insights into community structure and function, metagenome analyses applying high-throughput sequencing technologies were carried out [13-17]. Elaborate bioinformatics methods and analysis platforms facilitated metagenome sequence data interpretation and comparison [13,14,16,18]. Another approach for comparison of metagenome datasets is fragment recruitment of metagenomic sequences related to selected genomes of reference species. This approach provides insights into the degree of relatedness of indigenous species within a given habitat to known reference species. Recently, fragment recruitment has been applied for marine and silage microbial communities [19,20].

In principle, decomposition of biomass for biogas production can be practiced under wet or dry fermentation conditions. Wet fermentation is characterized by a high liquid content and a low total solid content, which usually is below 10%. In contrast to this, in dry fermentation the total solids content is between 15 and 35%. Biogas plants operating under dry fermentation conditions apply mostly maize silage, green rye (and similar biomass), dung (cow dung, poultry dry excrement and so forth) or municipal solid wastes as substrates without any continuous supplementation of liquid manure, which consequently leads to a low liquid content [1,21]. To control the water content, recirculation of digestate or liquid is applied, which may have a great impact on the activity of the underlying community. Recirculation may influence the pH, salt and organic loads, which could inhibit the microbial activity. Dry fermentation proved to be a convenient technology for the fermentation of substrates possessing relatively high dry matter contents. According to the German Renewable Energy Law (EEG), a technology bonus was granted for dry fermentation biogas plants built before the year 2008. Previously, a dry fermentation process of a production-scale biogas plant was characterized at the metagenomic level [13,14,16]. On the other hand, wet fermentation utilizing maize silage and liquid manure from cattle or swine is performed in most mid-sized, agricultural biogas plants in Germany. In this study, the microbial community of an agricultural biogas plant performing wet fermentation was analyzed by applying a metagenomic approach. Obtained results were compared to taxonomic community profiles deduced from a dry fermentation biogas plant analyzed previously. The present study adopted exactly the same methodology for processing samples, preparing total community DNA and metagenome sequence data analysis as for metagenome analysis of the dry fermentation process in the study mentioned above. It is hypothesized that biogas-producing microbial communities comprise a ‘core’ microbiome and variable sub-communities that respond to specific conditions and process parameters prevailing in particular reactor environments. The objective of this study was to differentiate biogas communities from biogas plants performing wet and dry fermentation, with respect to their taxonomic profiles, and to deduce correlations between these profiles and process parameters collected for both fermentation types. Another aim of this study was to identify key species specifically adapted to one process or the other and their predicted functions, focusing on methanogenic species.

Results and discussion

Analyzed biogas production plants differ in substrate input and chemical parameters

To compare taxonomic and functional profiles of the biogas-producing microbial communities from production-scale biogas plants operating under dry or wet fermentation conditions, samples from the primary digesters of two agricultural biogas plants differing in these fermentation types were analyzed. The biogas plant operating under dry fermentation conditions (BGP_DF) was sampled previously [16], whereas samples from the biogas plant operating under wet fermentation conditions (BGP_WF) were taken in March 2011. One of the major differences between BGP_WF and BGP_DF is their dry matter content (BGP_DF: 14% ± 2% and BGP_WF: 9% ± 1%). BGP_DF was fed with high amounts of plant silages and low amounts of chicken manure [16], whereas the substrate composition of BGP_WF mainly consisted of maize silage and a relatively high amount of liquid pig manure (Table 1). Moreover, both biogas plants can be distinguished according to their process parameters, such as volatile organic acids, total inorganic carbon, acetic and propionic acid concentrations and ammonium contents (see Table 1). BGP_DF is characterized by higher acetic and propionic acid concentrations as compared to BGP_WF, suggesting that consumption of these compounds is limited in BGP_DF. The biogas and methane yields (698.2 l/kg and 350.5 l/kg organic dry matter (oDM), see Table 1) are in the normal range of production in mesophilic production-scale biogas plants.

Table 1 Characteristics of the studied biogas plants performing wet or dry fermentation technology

In a recent study, a metagenome approach was carried out to study the taxonomic composition and functional potential of the microbial community in the biogas plant BGP_DF [13,16]. In total, 1,347,644 sequencing reads were generated with an average read length of 367.7 bases providing approximately 496 million bases sequence information (Table 2). Most biogas plants in Germany practice wet fermentation utilizing liquid manure and maize silage for the production of methane. To obtain insights into the microbial community composition of this process, a metagenome sequencing approach for the biogas plant applying wet fermentation (BGP_WF) was carried out. The same sample preparation, DNA-extraction method and sequencing technique were applied as previously described for BGP_DF to ensure comparability of the metagenome datasets. Sequencing of the samples originating from BGP_WF resulted in 1,532,780 sequencing reads, with an average read length of 387.3 bases, accounting for approximately 594 million bases sequence information (Table 2). To include only high quality sequences, both datasets were filtered for GC (G: Guanine, C: Cytosine bias and duplicates as previously described [22]. After this filtering step, 1,019,333 sequences from the BGP_DF and 1,097,549 sequences from the BGP_WF remained and were used for downstream taxonomic and functional analyses (see Table 2).

Table 2 Metagenome sequencing statistics of the DNA samples from the wet (BGP_WF) and dry (BGP_DF) fermentation biogas plants

Comparative analyses of taxonomic profiles obtained from wet and dry fermentation communities revealed high similarities

The community structures in the biogas plants operating under wet (BGP_WF) or dry fermentation (BGP_DF) conditions were studied using CARMA3 [23] and MetaSAMS [24]. The software CARMA3 was applied to calculate the taxonomic (microbial composition based on phylogenetic analyses) and functional profile (predicted Pfam protein families based on similarity searches) of the microbial community. CARMA3 is implemented in MetaSAMS (Metagenome Sequence Analysis and Management System), a software suited for the analysis of metagenome datasets. For the interpretation only taxonomic assignments with an E-value threshold of 10−5 were considered. In total, 711,293 sequences of the BGP_WF were assigned to a superkingdom representing 64.8% of the total number of analyzed sequences. In the BGP_DF dataset, 604,243 sequences (59.3% of the total sequences) were similar to known reference sequences at the rank superkingdom. Rarefaction analyses on the mean taxonomic richness on the taxonomic family rank showed a saturation from approximately 800,000 reads, which indicates that the majority of the microbial community has been captured, while the rarefaction analysis on the rank genus was nearly saturated (Additional file 1: Figure S1A and D).

At higher taxonomic ranks, the community compositions in BGP_WF and BGP_DF are very similar (Figure 1). Both communities are mainly composed of bacterial (59% in BGP_WF and 52.5% in BGP_DF of the total analyzed reads) and archaeal microorganisms (5.5% in BGP_WF and 6.6% in BGP_DF) (Figure 1). The bacterial superkingdom mainly comprises the phyla Firmicutes (32.5% in BGP_WF and 25.4% in BGP_DF), Bacteroidetes (8.4% in BGP_WF and 6.2% in BGP_DF), Spirochaetes (1.7% in BGP_WF and 0.5% in BGP_DF) and Proteobacteria (1.7% in BGP_WF and 1.2% in BGP_DF). These phyla were also common in other biogas plant microbiota [9,17,18]. Bacterial groups belonging to the taxa Firmicutes, Bacteroidetes and Spirochaetes are assumed to be involved in cellulolytic degradation, proteolysis, acidogenesis and homoacetogenesis [10]. Among the archaeal community, Euryarchaeota (5.2% in BGP_WF and 6.2% in BGP_DF of all analyzed sequences) is the most abundant phylum with Methanomicrobia as the largest class (Figure 1). Comparing both profiles on the phylum level, minor differences were observed in the relative abundances of Firmicutes. However, there are some noteworthy changes in the relative abundances of taxonomic groups between wet and dry fermentation digesters on the class level, in particular of classes belonging to the phylum Firmicutes.

Figure 1
figure 1

Taxonomic composition of microbial communities from the wet and dry fermentation processes. Taxonomic structure of the microbial community associated with a (A) dry (BGP_DF) or (B) wet (BGP_WF) fermentation process: The taxonomic composition of the microbial communities residing in the biogas plants analyzed was determined using CARMA3 and visualized by Krona plots. Only the ranks from superkingdom to family were considered.

Among the phylum Firmicutes in BGP_WF, Clostridia (19.5% of all analyzed sequences in BGP_WF) forms the largest class, followed by Bacilli (1.8%) and Erysipelotrichi (1.2%) (Figure 1). Likewise, Clostridia and Bacilli are also the most abundant classes of Firmicutes in the microbiome of the BGP_DF, with 14.3% and 1.4% of the total number of analyzed reads, respectively, whereas Erysipelotrichi is barely present (0.1%) in this digester. While Clostridia and Bacilli species are well described in the anaerobic digestion process in biogas plants, the evidence for Erysipelotrichi species is sparse [248]. Little is known about the family Erysipelotrichaceae. However, it was also identified in the microbial community of the gut [25]. An increase in members of this family was associated with a diet high in fat, increased body weight and decreased fecal short-chain fatty acid concentrations in mice [25]. Concurrently, occurrence of Erysipelotrichi members correlates with lower short-chain fatty acid concentrations in BGP_WF. Whether this observation is really due to the metabolic capabilities of this group of microorganisms remains to be determined.

In both biogas plants, the most abundant families belonging to the class Clostridia are Clostridiaceae (3.8% in BGP_WF and 1.8% in BGP_DF of all analyzed sequences), Ruminococcaceae (0.5% in BGP_WF and 0.3% in BGP_DF) and Lachnospiraceae (0.3% in BGP_WF and 0.1% in BGP_DF) (Figures 1 and 2). Clostridium, the prevalent genus within Clostridiaceae, seems to belong to the core set of organisms, as it is dominant in both biogas plants studied (Figure 1). Species of this genus, such as Clostridium thermocellum [26] and Clostridium clariflavum [27], produce cellulosomes, an extracellular multi-enzyme complex which is important for the degradation of complex carbohydrates such as cellulose. Indeed, environmental gene tags (EGTs) classified to the genus Clostridium encode enzymes relevant in the hydrolysis process of glycoside bonds (PF00150 and PF00759). In both biogas plants, Alkaliphilus is the second largest genus within the family Clostridiaceae. It has also been detected in high amounts in a biogas plant fed with plant biomass and pig manure slurry [17]. The species Alkaliphilus peptidofermentans, isolated from a soda lake, is described to ferment peptides to acetate and formate [28]. The functional profile of EGTs assigned to the genus Alkaliphilus in the wet fermentation process includes various peptidase families (PF00768 and PF05343).

Figure 2
figure 2

Scatterplot of sequence counts assigned to the taxonomic rank family for microbial communities from the wet and dry fermentation processes. For each taxonomic family, sequence abundances, normalized on the smallest dataset, with pseudocounts (for an explanation see Methods section) are plotted as a logarithm to the base 10. Red symbols indicate archaeal and green symbols indicate bacterial families. Only families for which the sum of the relative read abundances measured was at least 0.01% were considered.

Members of the genus Ruminococcus belonging to the family Ruminococcaceae (order Clostridiales) are cellulolytic organisms, which occur in the rumen community [29] and in biogas plants [17]. They play an important role in the digestion of plant cell wall material and produce acetate [29,30]. Among the functional profile of EGTs assigned to Ruminococcus, protein families for glycosyl hydrolase family 9, 26 and 48 (PF00759, PF02156 and PF04616) representing cellulose-, cellobiose- and hemicellulose-degrading enzymes were observed. Furthermore, gene fragments for dockerin and cohesin (PF00404 and PF00963), the scaffoldin units of cellulosomes, were detected in the functional profile of the wet fermentation process. Recently, the genome of Acetivibrio cellulolyticus of the family Ruminococcus was sequenced [31]. The genome carries genes for a complex cellulosome system including endoglucanases and cellobiohydrolases. Glycosyl family 8 and 9, both representing these enzymes, were detected among the EGTs assigned to the genus Acetivibrio.

The family Spirochaetaceae dominates the phylum Spirochaetes with 68% of the reads assigned to this phylum in BGP_DF (Figure 1). It is the third largest family, with 1.2% of all reads in the wet fermentation sample, whereas it belongs to the minor groups in the dry fermentation process (0.2%). Microorganisms related to Treponema have already been described in a mesophilic biogas digester treating pig manure using 16S rRNA clone libraries [10]. Genomes of the Treponema species encode proteolytic enzymes [32] and glycoside hydrolases [33]. EGTs classified to the genus Treponema in BGP_WF were also assigned to carbohydrate phosphorylase (PF00343), alpha amylase (PF00128) and 4-alpha-glucanotransferase (PF02446), which participate in starch and sucrose metabolism, carbohydrate phosphorylase (PF00343) and alphaamylase (PF00128). The family Succinivibrionaceae belonging to the Gammaproteobacteria is only predicted in the biogas plant operating wet fermentation (0.04% of all analyzed sequenced reads) and was sparsely detected in the dry fermentation process (Figures 1 and 2). This family is noted to use glucose and other carbohydrates as an energy source and to produce succinate and acetate [34]. Likewise, the family Fibrobacteraceae is more frequently present in the wet fermentation process. The known species Fibrobacter succinogenes is described to play a key role in the degradation process of cellulose [35].

Prevotellaceae belonging to the phylum Bacteroidetes is more abundant in the BGP_WF, with 0.18% of all analyzed sequences in the corresponding dataset as compared to 0.8% of all sequences in the BGP_DF (Figures 1 and 2). Species of the genus Prevotella are described to decompose hemicellulose, starch and pectin [36,37]. Various Pfam families were discovered in the wet fermentation digester that were assigned to Prevotella and are predicted to be involved in the degradation of hemicellulose (such as glycosyl hydrolase family 3 (PF00933, PF01915), glycosyl hydrolase 92 (PF07971) and alpha-L-arabinofuranosidase (PF06964)) or starch (such as glycosyl hydrolase family 31 (PF01055)).

Synergistaceae of the phylum Synergistetes belongs to the core-set of the families represented in both biogas plants (0.4% in BGP_WF and 0.1% in BGP_DF of all analyzed reads) (Figures 1 and 2). Anaerobaculum and Aminobacterium are genera within the phylum Synergistetes, which are both predicted in the biogas plants. Recently, Anaerobaculum was detected in anaerobic digestion of slaughterhouse waste mixture [38] and thermophilic sludge [11]. Species of Anaerobaculum are associated with the fermentation of peptides and produce short-chain fatty acids. Likewise, species of the genus Aminobacterium were isolated from anaerobic sludge and are described to ferment a range of amino acids to acetate, propionate and hydrogen [39,40].

Community profiles of the dry and wet fermentation process differ at lower taxonomic ranks

There are also a number of taxa that are slightly increased in the biogas plant operating the dry fermentation technology (Figures 1 and 2). For example, the two biogas plants differ in the proportion of sequences belonging to the family Acholeplasmataceae, of the phylum Tenericutes. The species Acholeplasma laidlawii was isolated from wastewater [15] and also has been identified in other biogas plants [1]. As a source for carbon, A. laidlawii utilizes glucose, fructose and galactose [15]. The genome of A. laidlawii harbors genes for enzymes that degrade starch, amino sugars and other sugars. In the functional profile of EGTs assigned to this family are glycoproteases (PF00814) and peptidases (PF01546) (data not shown).

In the dry fermentation process, Candidatus Cloacamonas is more prevalently present (Figure 2). The species Candidatus Cloacamonas acidaminovorans was recently detected to be highly abundant in other anaerobic digesters [17,41]. Previously, the genus was also identified in a 16S rRNA clone library of the same biogas plant [42]. Proteome analysis indicated that C. acidaminovorans might attain energy from sugars in the Embden-Meyerhof pathway and from the fermentation of amino acids, and thereby produces hydrogen and carbon dioxide [41].

Overall, the taxonomic profiles are similar for the wet and dry anaerobic digestion, especially at higher taxonomic ranks. In the BGP_WF, the families Erysipelotrichaceae, Fibrobacteraceae, Succinivibrionaceae and Clostridiaceae were found to be more abundant, whereas more sequences were assigned to Acholeplasmataceae and Candidatus Cloacamonas in the BGP_DF.

Hydrogenotrophic methanogens are dominant in the dry and wet fermentation processes

A detailed taxonomic analysis was performed for the archaeal sequences of both fermentation processes. In the wet fermentation process, less EGTs were assigned to Methanomicrobiaceae (3.4%) as compared to the dry fermentation digester (4.4%) (Figures 1 and 2). In both fermentation processes, Methanomicrobiaceae is the most abundant methanogenic family with Methanoculleus being the prevalent genus. Methanoculleus species conduct the hydrogenotrophic methanogenesis pathway, synthesizing methane from carbon dioxide and hydrogen. Further identified families that are capable of performing methanogenesis are Methanosarcinaceae, Methanosaetaceae and Methanobacteriaceae [43]. These families are only detected in low abundance in the microbial communities. In the BGP_WF, the acetoclastic methanogen Methanosaeta was measured in slightly higher frequencies as compared to the dry fermentation process (Figure 2). This genus was noted to be more dominant in biogas plants with low acetate concentrations [44], as it has a high affinity to acetate. The observation is supported by this study as the abundance of Methanosaeta correlates with the lower concentration of acetic acid in BGP_WF compared to BGP_DF (see Table 1).

In addition to the taxonomic comparison of the biogas plants, enzymes involved in the acetoclastic and hydrogenotrophic methanogenesis were categorized according to Pfam families and were searched in the functionally characterized metagenomes obtained from BGP_WF and BGP_DF. Reads assigned to the selected Pfam families and classified to the superkingdom Archaea were used for the subsequent analysis (see Figure 3). No EGTs were assigned to acetate kinase and phosphotransacetylase, which are key enzymes in the initial step of acetoclastic methanogenesis in Methanosarcina [45]. However, EGTs for all enzymes relevant in the hydrogenotrophic methanogenesis pathway were identified with a slightly higher amount in the dry fermenter (data not shown). Rarefaction analyses on the mean EGT richness showed a saturation at approximately 800,000 reads, which indicates that the majority of the microbial community has been captured (see Additional file 1: Figure S1B). Moreover, we emphasized the EGTs that are unique in each of the biogas plants’ metagenomes, which means that they are absent in the one plant while they have an abundance of at least five in the other plant (normalized on the smallest dataset). This analysis showed that there are more unique EGTs in the wet fermentation process metagenome. However, none of the identified unique EGTs can be correlated to functions directly associated with the fermentation processes leading to methane production (Additional file 2: Table S1).

Figure 3
figure 3

Comparative analysis of archaeal environmental gene tags (EGTs) representing methanogenesis-related enzymes in metagenome datasets for the wet or dry fermentation process. Comparison of the relative abundances of EGTs representing enzymes involved in the acetyl-CoA pathway, methanogenesis and Wood-Ljungdahl pathway. Only EGTs assigned to Archaea were considered and relative abundances were normalized based on the smallest dataset.

Differentiation of the methanogenic sub-community within the wet and dry fermentation process

To analyze the archaeal sub-community at a higher resolution, high-throughput 16S rRNA gene amplicon sequencing was carried out. The sequencing procedure for BGP_WF yielded 83,719 sequencing reads. A previous rRNA gene amplicon sequencing run of the sample obtained from BGP_DF resulted in 170,941 sequencing reads [46]. In total, 37 operational taxonomic units (OTUs) were of archaeal origin, representing 2,753 sequences of the BGP_WF dataset (8%). After manual inspection for chimeric sequences and excluding those sequences with a length below 100 bp, 15 OTUs remained comprising 2,691 reads. In the dry fermentation process, 22 OTUs representing 1,118 sequences (2%) were assigned to Archaea. After removing manually detected chimeric and short (<100 bp) sequences, nine OTUs remained representing 1,095 sequences. The number of sequences included in the OTUs was normalized to the number of the smallest filtered 16S rRNA gene amplicon dataset. Phylogenetic analysis of representative OTU sequences that were assigned to the superkingdom Archaea was carried out using FastTree. Rarefaction analyses on the mean OTU richness showed saturation at approximately 30,000 reads, which indicates that the majority of the microbial community has been captured (see Additional file 1: Figure S1C).

The largest archaeal OTUs of the wet fermentation and dry fermentation plant are similar to Methanoculleus bourgensis with 2,184 and 680 sequences, respectively (Figure 4). The presence of a phylogenetic cluster related to M. bourgensis is in agreement with previous studies based on 16S rRNA clone library analyses obtained from the biogas plant BGP_DF [42]. The representative sequences of the second largest OTU of each biogas plant cluster together. The representative amplicon sequences comprise 359 and 51, respectively, and are located in a large cluster formed by known Methanoculleus species. However, no sequence of a described Methanoculleus species could be identified in close proximity of these representative sequences, indicating a so far unknown archaeal species related to Methanoculleus (unknown Methanoculleus cluster II).

Figure 4
figure 4

Archaeal phylogenetic tree of representative operational taxonomic units (OTUs) sequences from the wet and dry fermentation processes. Representative OTU sequences and corresponding counts assigned from the wet (BGP_WF, red) or dry fermentation (BGP_DF, green) process are labeled with an OTU identifier in brackets, and the nodes are labeled with the respective boot strap values. Representative sequences were assigned to one known Methanobacteriaceae and three unknown Methanomicrobiaceae clusters (I to III), and to one known and two unknown Methanoculleus bourgensis clusters (I and II), respectively. Sequence counts were normalized according to the smallest dataset.

A further unknown Methanoculleus species is highly represented in the wet fermentation plant (unknown Methanoculleus cluster I). The corresponding OTU contains 98 sequence reads. A representative sequence of the dry fermentation biogas plant clusters close to the latter OTU, comprising only four reads. A so far unknown Methanoculleus species has also been described in the same biogas plant based on 16S rRNA library clones [42].

Moreover, representative sequences are located outside the phylogenetic cluster formed by the described Methanoculleus species. Corresponding archaeal species also belong to the class Methanomicrobia (cluster I and cluster II). The phylogenetic tree distantly affiliated a sequence from the dry fermentation plant to the recently identified species Methanomassiliicoccus luminyensis B10 from a human gut microbiome sample [47]. The representative sequence also clusters close to a so far unknown archaeal clone [48], which originates from an agricultural biogas plant supplied with water, maize silage and barley grains [49]. The distribution of similar sequences in various habitats suggests a wide occurrence of species related to M. luminyensis.

Besides Methanomicrobia, Methanobacteria were identified in the biogas plants based on phylogenetic analysis. Corresponding species are related to Methanobrevibacter and Methanobacterium. Both were identified in the 16S rRNA gene amplicon dataset, as well as in the whole metagenome approach, in the two biogas plants with minor frequencies. Some representative sequences (BGP_DF_8, BGP_WF_7 and BGP_WF_14) form a phylogenetic cluster with the 16S rRNA gene sequence of Methanobacterium kanagiense [50], a hydrogenotrophic archaeon isolated from an anaerobic propionate-oxidizing soil sample. Another representative sequence (BGP_WF_12) is 96% identical to the 16S rRNA gene sequence of Methanobrevibacter olleyae [51]. This species is present in sheep and bovine rumen and was described to use hydrogen and formate for methane formation. Further representative sequences (BGP_DF_9 and BGP_WF_15) are associated with Methanobrevibacter woesei and feature an identity of 97 to 98% to this species. M. woesei was isolated from goose gut and is a hydrogenotrophic archaeon that mainly uses hydrogen and carbon dioxide for methane production.

Finally, no representative sequences related to acetoclastic methanogens were identified in the biogas plants, showing a dominance of hydrogenotrophic methanogens in the biogas production process of these plants.

Overall, the constructed archaeal phylogenetic tree illustrates differences in the composition of the methanogenic sub-communities of both biogas plants. Besides the dominating M. bourgensis further Methanoculleus species are also present in the dry fermentation process (unknown Methanoculleus cluster I and II). In BGP_WF other methanogens comprising species related to unknown Methanomicrobiaceae species (cluster II and cluster III (see Figure 4) are prominent.

Fragment recruitments revealed that methanogens of the dry fermentation plant are closely related to the type strain Methanoculleus bourgensis MS2T at the genomic level

To analyze the degree of relatedness of biogas-producing community members to completely sequenced reference microorganisms, and to differentiate the metagenome datasets of both biogas plants referring to this, fragment recruitments were conducted as described previously [19]. Another objective of this analysis was to determine similarities of methanogens of both biogas plants to the genome of the type strain M. bourgensis MS2T [43]. For this purpose, metagenome sequence reads were searched for matches to completely sequenced microbial genomes stored in the National Center for Biotechnology Information (NCBI) database applying the BLASTn algorithm [Basic Local Alignment Search Tool on nucleotide level]. Counts of metagenome reads that match to reference genomes with more than 90% sequence similarity are listed for the top 20 genomes in Table 3. Reference organisms appearing in this analysis taxonomically represent the orders Methanomicrobiales, Clostridiales, Lactobacillales, Thermoanaerobacterales and Synergistales and could be classified to belong to one of the following functional groups: cellulolytic organisms, secondary fermenters (acidogenic), syntropic organisms (acetogenic) and methanogens (see Table 4). Also, for BGP_DF, the orders Thermotogales and Bacillales were identified among the top 20 recruitments.

Table 3 List of 20 reference genomes showing the highest similarities to the metagenome datasets from biogas plants operating under wet or dry fermentation conditions as analyzed by fragment recruitment analysis
Table 4 Species identified during a fragment recruitment analysis using metagenome sequence data of wet and dry fermentation biogas plants

Additionally, BLASTn results were graphically evaluated in fragment recruitment plots, in which the degree of similarity of a single hit is plotted against the position of this hit on the reference genome sequence. This was done for those three microorganisms showing the highest sequence similarities to matching metagenome sequence reads (see Table 3), namely M. bourgensis MS2T, C. clariflavum DSM 19732 and C. thermocellum ATCC 27505 (see Figure 5). Hit distribution plots show the numbers of hits with a similar degree of similarity to the reference genome sequence. Fragment recruitments revealed that methanogens within fermentation samples of both biogas plants are related to the reference type strain M. bourgensis MS2T. It was observed that dominant methanogens of BGP_DF are more closely related to the reference genome as compared to those of BGP_WF (see Figure 5A). It was supposed that M. bourgensis may be very well adapted to adverse conditions prevailing in BGP_DF, such as high ammonium concentrations and osmolarity [12,77,78], due to the fed substrate, which was chicken dry excrement.

Figure 5
figure 5

Fragment recruitment plots and hit distribution for metagenome sequence reads from the wet and dry fermentation processes on three reference genomes. Visualization of BLASTn analyses of metagenomic sequence reads to the reference genomes of Methanoculleus bourgensis MS2T (A), Clostridium clariflavum DSM 19732 (B) and Clostridium thermocellum ATCC 27405 (C). Within the fragment recruitment plot (I) the sequence identities (>55%) between each hit of a metagenomic sequence read and the chromosomal reference sequence are plotted against the position of the alignment. In the hit distribution plot (II) the normalized numbers of reads featuring hits to the reference genome were plotted in intervals of 1% for metagenome reads showing 55 to 100% sequence identity to the reference.

Surprisingly, M. bourgensis MS2T represents the only organism that can be regarded as reasonable reference for members of the biogas-producing communities analyzed in this study. All other microorganisms listed in Table 3 feature a lesser degree of conservation to biogas community members and hence are more distantly related to them (exemplarily see Figure 5B,C). Moreover, differentiation of the microbial communities of both biogas plants is not apparent from these fragment recruitment analyses, since results obtained for both metagenome datasets uncovered nearly the same set of prominent reference organisms, with only slight differences in rankings (see Tables 3 and 4). This result confirms deduced taxonomic profiles for both biogas plants at higher taxonomic ranks which did not uncover pronounced differences between both microbial communities (see previous sections). However, it was observed that the diversity within the genus Streptococcus is greater for BGP_WF which was fed with swine manure. It is known that particular Streptococcus species such as Streptococcus suis and Streptococcus pasteurianus belong to the common microbiome of the pig’s intestinal tract [52,79,80]. Moreover, it should be noted here that Desulfotomaculum carboxydivorans was identified at rank nine for BGP_WF and at rank eight for BGP_DF in fragment recruitments. This bacterium is able to catabolize low-molecular weight compounds (pyruvate, lactate, glucose etc.) and reduces sulfate as terminal electron acceptor to sulfide (hydrogen sulfide), which is deleterious in biogas since it leads to corrosion of combustion units and in pipelines.

In summary, fragment recruitments revealed that representation of biogas community members in databases is currently insufficient. The only exception is M. bourgensis MS2T which is very similar to dominant members of the methanogenic sub-community, at least in BGP_DF.

Conclusions

The methane-producing community of an agricultural biogas plant applying wet fermentation technology was characterized by means of a metagenomic approach and taxonomically compared with the microbial community of a biogas plant applying dry fermentation that was previously characterized [16]. This is the first study in which metagenomes of different biogas plants are compared using the same DNA preparation technique, sequencing technology and bioinformatics methods. Both communities’ metagenomes were sequenced using 454-pyrosequencing, and the sequence data were analyzed based on their taxonomic composition and functional profile, with focus on the methanogenic sub-community and genes involved in methanogenesis, using the same bioinformatics tools and pipelines. Despite differences in the process parameters, such as acetic acid concentration and pH value and the biogas plant’s substrate compositions, the microbial communities of BGP_WF and BGP_DF are similar in their composition on higher taxonomic ranks. Only a minor number of taxa on lower ranks differed between BGP_WF and BGP_DF. Accordingly, the majority of the taxa belonged to the core set of microorganisms residing in both biogas plants. Even though in both fermenters, the family of Methanomicrobiaceae is most dominant and Methanoculleus is the prevalent genus within this family, the composition of this genus differs between the different plants. As fragment recruitments revealed, within the genus Methanoculleus, M. bourgensis was the most dominant species in dry fermentation while this was not the case for the wet fermentation process.

The presented results clearly indicate that the hypothesis of a core community being present in biogas fermenters, operating either under wet or dry fermentation conditions, can be confirmed. This is in accordance with several studies examining the microbial taxonomic composition in biogas fermenters. Identified differences in taxonomic profiles can be associated with different process parameters of the fermenters, such as fed substrate, pH, acid and ammonium concentrations, as especially methanogens are sensitive regarding changes in acetate and propionate concentrations.

To further verify the hypothesis of a core microbiome and an adaptation of specific species to particular process parameters, further comparative studies with more biogas fermenters are required. Analyzing microbiota in biogas plants with altered physiochemical characteristics (temperature, pH and concentration of relevant metabolites) and substrate supply would aid in identifying the core species essential for the anaerobic digestion process, which in turn would provide information to control the biogas production process and prevent unstable conditions. Moreover, the metatranscriptome of the microbial biogas communities should be studied and compared to corresponding metagenomes to enable characterization not only of the taxonomic composition, but also to deduce the actual metabolic activity within the biogas fermenters.

Methods

Total community DNA isolation, purification and sequencing

The biogas plant featuring the mesophilic continuous dry fermentation technology (hereafter noted as BGP_DF) was designed for a capacity of 530 kWel (combined heat and power (CHP)) and a daily input of maize silage (63%), green rye (35%) and chicken manure (2%), divided into 24 feedings per day. The process comprises two digesters; the primary digester (BIOGAS NORD GmbH, Bielefeld, Germany) (1,557 m3 net volume, height of 6 m, diameter of 19 m) has an organic load of 4.8 kg oDM m−3 d−1, a theoretical retention time of 59 days and a temperature of 40°C. At the end, the digestate is stored in a closed non-heated final storage reactor (BIOGAS NORD GmbH, Bielefeld, Germany) (2,987 m3 net volume, height of 6 m, diameter of 26 m). The biogas and methane yields at the time of sampling were at 698.2 and 350.5 l/kg oDM, respectively.

The biogas plant applying the mesophilic wet fermentation technology (hereafter BGP_WF) was designed for a capacity of 537 kWel CHP. The daily input of maize silage (approximately 72%) and liquid pig manure (approximately 28%) was divided into 24 feedings per day. The biogas plant is composed of two digesters and the storage tank (BIOGAS NORD GmbH, Bielefeld, Germany). The digester (2,041 m3 net volume, height of 6.4 m, diameter of 21 m) has an organic load of 4.0 kg oDM m−3 d−1, a theoretical retention time of 55 days and a temperature of 40°C. The digestate is stored in a closed non-heated final storage tank (4,742 m3 net volume, height of 6 m, diameter of 32 m). The biogas and methane yields at the time of sampling were at 810.5 and 417.8 l/kg oDM, respectively.

Samples were taken from the primary digester of BGP_WF and total community DNA was extracted in triplicates applying the same procedure as described previously for the biogas plant featuring dry fermentation (BGP_DF). The triplicates were pooled prior to sequencing applying high-throughput sequencing using the Genome Sequencer FLX system (Roche Diagnostics GmbH, Mannheim, Germany) [22]. The BGP_WF metagenomic DNA was sequenced omitting the additive for high-GC DNA in order to ensure comparability between the metagenomes of the BGP_DF and BGP_WF. Subsequently, both sequence datasets were processed to remove emulsion-PCR duplicates and sequences affected by GC-biases in the course of sequencing, as described previously [22].

Taxonomic and functional analysis of metagenome sequence data

The filtered metagenome sequences of BGP_DF and BGP_WF were imported into MetaSAMS [81] and taxonomically characterized using CARMA3 [23]. CARMA3 computes taxonomic and functional assignments for EGTs derived from a microbial community. The filtered metagenome sequences of BGP_DF and BGP_WF were imported into MetaSAMS to apply CARMA3, which is implemented in MetaSAMS. Moreover, MetaSAMS enables exploration, analysis, management and visualization of calculated observations for metagenome sequences. Using MetaSAMS, rarefaction curves were calculated based on the levels ‘family’ and ‘genus’. Rarefaction analysis addresses the assessment of ‘taxon’ richness from different sub-sample sizes regarding metagenome sequence reads. In a rarefaction curve, the number of assigned taxa (on a specified level) is plotted as a function of the number of sequences within a selected sub-sample. The functional profile was calculated by searching for Pfam families matching with an E-value threshold of 10−5. Subsequently, the functional profile for selected taxonomic units was evaluated regarding their functional roles in the anaerobic digestion process.

To determine whether the relative abundance of specific taxa changed in the biogas fermenters, abundances were visualized in a scatter plot. First, the absolute counts of reads assigned to a specific family were normalized according to the smallest dataset. Then, the logarithm was used to decrease the influence of more dominant families. To allow the logarithm of zero, a pseudo-count of one sequence was added for each family prior to the logarithm. Moreover, differences in the methanogenesis step were studied. Briefly, EGTs which were taxonomically characterized as originating from archaeal species and functionally assigned to genes involved in the acetyl-CoA, methanogenesis or Wood-Ljungdahl pathway were counted and normalized according to the smallest sample [19]. The metagenomic sequence data of BGP_WF can be found at the European Molecular Biology Laboratory - European Bioinformatics Institute (EMBL-EBI) database under the accession number [EMBL: PRJEB5813].

16S rDNA amplicon generation, sequencing and analysis

16S rRNA gene amplicons were generated and sequenced as described recently [24], applying high-throughput sequencing using the Genome Sequencer FLX system (Roche Diagnostics GmbH, Mannheim, Germany). First, a PCR was performed to amplify a region covering V3 and V4 using universal 16S rDNA primers. Next, gel electrophoresis and gel extraction were applied to obtain only amplicons with the correct length. The PCR was repeated in order to attach barcode tags as well as adaptors to the amplicons. Finally, the amplicons were sequenced on a 454 Genome Sequencer (GS) device using FLX Titanium chemistry.

16S rRNA amplicon sequences of the BGP_DF and BGP_WF were simultaneously processed using the QIIME (Quantitative Insights Into Microbial Ecology) pipeline [82]. First, barcode and primer sequences were removed allowing 0 and 2 mismatches, respectively, and sequences with ambiguous bases were discarded. The option ‘truncate_only’ was used meaning that reverse primer sequences will be trimmed only if they are identified at the end of the amplicon sequences. To obtain high quality sequences for phylogenetic analyses, strict quality processing was carried out (window size of 25 bases, average quality score 25). Subsequently, the software package USEARCH version 6.0 was applied for denoising, chimera detection (de novo mode) and clustering into OTUs based on a 97% sequence identity [83,84]. Afterwards, representative sequences were selected for each cluster and assigned to taxonomic groups using the RDP Classifier 2.5 [85]. Only assignments with a confidence value of at least 0.8 were considered. Rare-faction curves based on OTUs were calculated to determine the coverage of the microbial community by the sequenced metagenome reads. For phylogenetic analysis, representative sequences assigned to Archaea with a confidence value of at least 0.8 were selected and aligned with the Infernal 1.1 software [86] using the Archaea SSU rRNA model (RF01959) from Rfam [48]. Finally, the alignment was used as a basis for tree reconstruction using FastTree [87]. The tree was rooted with the sequence of the Crenarchaeota Sulfolobus acidocaldarius covering the V3-V4 region.

Fragment recruitments

Fragment recruitments were performed as described previously [19,20] to compare the relatedness of metagenomic reads to the genomes of reference microorganisms. BLASTn analyses of metagenomic reads against a database containing all genomes of completely sequenced microorganisms were accomplished. Identity values of hits were computed by dividing the number of identical bases by the sequence length of the read. Hits with an identity value below 55% were discarded. Results of the BLASTn analysis were then visualized by plotting the calculated identity of each sequence read against the alignment position on the reference sequence. Moreover, a histogram was generated to display the distribution of hit identities. Numbers of hits displayed in fragment recruitments and in histograms were normalized based on the smallest sample size.

Abbreviations

BGP_DF:

Biogas plant operating under dry fermentation conditions

BGP_WF:

Biogas plant operating under wet fermentation conditions

EGT:

Environmental gene tag

oDM:

Organic dry matter

OTU:

Operational taxonomic unit

References

  1. Weiland P. Biogas production: current state and perspectives. Appl Microbiol Biotechnol. 2010;85:849–60.

    Article  CAS  Google Scholar 

  2. Yadvika S, Sreekrishnan TR, Kohli S, Rana V. Enhancement of biogas production from solid substrates using different techniques - a review. Bioresour Technol. 2004;95:1–10.

    Article  CAS  Google Scholar 

  3. Ohmiya K, Sakka K, Kimura T. Anaerobic bacterial degradation for the effective utilization of biomass. Biotechnol Bioprocess Eng. 2005;10:482–93.

    Article  CAS  Google Scholar 

  4. Braun R. Biogas from Energy Crop Digestion. In: IEA Task 37 Brochure. Paris, France: International Energy Agency; 2009.

    Google Scholar 

  5. Ács N, Kovács E, Wirth R, Bagi Z, Strang O, Herbel Z, et al. Changes in the Archaea microbial community when the biogas fermenters are fed with protein-rich substrates. Bioresour Technol. 2012;131:121–7.

    Article  Google Scholar 

  6. Kampmann K, Ratering S, Kramer I, Schmidt M, Zerr W, Schnell S. Unexpected stability of Bacteroidetes and Firmicutes communities in laboratory biogas reactors fed with different defined substrates. Appl Environ Microbiol. 2012;78(7):2106–19.

    Article  CAS  Google Scholar 

  7. Wei H, Tucker MP, Baker JO, Harris M, Luo Y, Xu Q, et al. Tracking dynamics of plant biomass composting by changes in substrate structure, microbial community, and enzyme activity. Biotechnol Biofuels. 2012;5:20.

    Article  CAS  Google Scholar 

  8. Bengelsdorf FR, Gerischer U, Langer S, Zak M, Kazda M. Stability of a biogas-producing bacterial, archaeal and fungal community degrading food residues. FEMS Microbiol Ecol. 2012;84(1):201–12.

    Article  Google Scholar 

  9. Klocke M, Mähnert 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.

    Article  CAS  Google Scholar 

  10. Liu FH, Wang SB, Zhang JS, Zhang J, Yan X, Zhou HK, et al. The structure of the bacterial and archaeal community in a biogas digester as revealed by denaturing gradient gel electrophoresis and 16S rDNA sequencing analysis. J Appl Microbiol. 2009;106(3):952–66.

    Article  CAS  Google Scholar 

  11. Weiss A, Jerome V, Freitag R, Mayer HK. Diversity of the resident microbiota in a thermophilic municipal biogas plant. Appl Microbiol Biotechnol. 2008;81(1):163–73.

    Article  CAS  Google Scholar 

  12. Weiss A, Jérôme V, Burghardt D, Likke L, Peiffer S, Hofstetter EM, et al. Investigation of factors influencing biogas production in a large-scale thermophilic municipal biogas plant. Appl Microbiol Biotechnol. 2009;84:987–1001.

    Article  CAS  Google Scholar 

  13. Jaenicke S, Ander C, Bekel T, Bisdorf R, Dröge M, Gartemann K-H, et al. Comparative and joint analysis of two metagenomic datasets from a biogas fermenter obtained by 454-pyrosequencing. PLoS One. 2011;6:e14519.

    Article  CAS  Google Scholar 

  14. Krause L, Diaz NN, Edwards RA, Gartemann K-H, Krömeke H, Neuweger H, et al. Taxonomic composition and gene content of a methane-producing microbial community isolated from a biogas reactor. J Biotechnol. 2008;136(1–2):91–101.

    Article  CAS  Google Scholar 

  15. Li A, Chu Y, Wang X, Ren L, Yu J, Liu X, et al. A pyrosequencing-based metagenomic study of methane-producing microbial community in solid-state biogas reactor. Biotechnol Biofuels. 2013;6:3.

    Article  CAS  Google Scholar 

  16. Schlüter A, Bekel T, Diaz NN, Dondrup M, Eichenlaub R, Gartemann K-H, et al. The metagenome of a biogas-producing microbial community of a production-scale biogas plant fermenter analysed by the 454-pyrosequencing technology. J Biotechnol. 2008;136(1–2):77–90.

    Article  Google Scholar 

  17. Wirth R, Kovacs E, Maróti G, Bagi Z, Rákhely G, Kovács KL. Characterization of a biogas-producing microbial community by short-read next generation DNA sequencing. Biotechnol Biofuels. 2012;5:41.

    Article  CAS  Google Scholar 

  18. Rademacher A, Zakrzewski M, Schlüter A, Schönberg M, Szczepanowski R, Goesmann A, et al. Characterization of microbial biofilms in a thermophilic biogas system by high-throughput metagenome sequencing. FEMS Microbiol Ecol. 2012;79(3):785–99.

    Article  CAS  Google Scholar 

  19. Eikmeyer F, Köfinger P, Poschenel A, Jünemann S, Zakrzewski M, Heinl S, et al. Metagenome analyses reveal the influence of the inoculant Lactobacillus buchneri CD034 on the microbial community involved in grass silaging. J Biotechnol. 2013;167(3):334–43.

    Article  CAS  Google Scholar 

  20. Rusch DB, Halpern AL, Sutton G, Heidelberg KB, Williamson S, Yooseph S, et al. The Sorcerer II global ocean sampling expedition: northwest Atlantic through eastern tropical Pacific. PLoS Biol. 2007;5(3):e77.

    Article  Google Scholar 

  21. Weiland P. Biomass digestion in agriculture: a successful pathway for the energy production and waste treatment in Germany. Eng Life Sci. 2006;6(3):302–9.

    Article  CAS  Google Scholar 

  22. Jans C, Gerber A, Bugnard J, Njage PM, Lacroix C, Meile L. Novel Streptococcus infantarius subsp. infantarius variants harboring lactose metabolism genes homologous to Streptococcus thermophilus. Food Microbiol. 2012;31(1):33–42.

    Article  CAS  Google Scholar 

  23. Gerlach W, Stoye J. Taxonomic classification of metagenomic shotgun sequences with CARMA3. Nucleic Acids Res. 2011;39:e91.

    Article  CAS  Google Scholar 

  24. Zakrzewski M, Bekel T, Ander C, Pühler A, Rupp O, Stoye J, et al. MetaSAMS-A novel software platform for taxonomic classification, functional annotation and comparative analysis of metagenome datasets. J Biotechnol. 2012;167(2):156–65.

    Article  Google Scholar 

  25. Fleissner CK, Huebel N, Abd El-Bary MM, Loh G, Klaus S, Blaut M. Absence of intestinal microbiota does not protect mice from diet-induced obesity. Br J Nutr. 2010;104(6):919–29.

    Article  CAS  Google Scholar 

  26. Brown SD, Lamed R, Morag E, Borovok I, Shoham Y, Klingeman DM, et al. Draft genome sequences for Clostridium thermocellum wild-type strain YS and derived cellulose adhesion-defective mutant strain AD2. J Bacteriol. 2012;194(12):3290–1.

    Article  CAS  Google Scholar 

  27. Izquierdo JA, Goodwin L, Davenport KW, Teshima H, Bruce D, Detter C, et al. Complete genome sequence of Clostridium clariflavum DSM 19732. Stand Genomic Sci. 2012;6(1):104–15.

    Article  CAS  Google Scholar 

  28. Zhilina TN, Zavarzina DG, Kolganova TV, Lysenko AM, Turova TP. New alkaliphilic peptide fermenting and Fe(III) regenerating Bacteria Alkaliphilus peptidofermentans sp.nov. from soda lake. Mikrobiologiia. 2012;78(4):496–505.

    Google Scholar 

  29. Berg Miller ME, Antonopoulos DA, Rincon MT, Band M, Bari A, Akraiko T, et al. Diversity and strain specificity of plant cell wall degrading enzymes revealed by the draft genome of Ruminococcus flavefaciens FD-1. PLoS One. 2009;4(8):e6650.

    Article  Google Scholar 

  30. Miller TL, Wolin MJ. Bioconversion of cellulose to acetate with pure cultures of Ruminococcus albus and a hydrogen-using acetogen. Appl Environ Microbiol. 2012;61(11):3832–5.

    Google Scholar 

  31. Dassa B, Borovok I, Lamed R, Henrissat B, Coutinho P, Hemme CL, et al. Genome-wide analysis of Acetivibrio cellulolyticus provides a blueprint of an elaborate cellulosome system. BMC Genomics. 2012;13:210.

    Article  CAS  Google Scholar 

  32. Hu P, Yang M, Zhang A, Wu J, Chen B, Hua Y, et al. Complete genome sequence of Streptococcus suis serotype 14 strain JS14. J Bacteriol. 2012;193(9):2375–6.

    Article  Google Scholar 

  33. Kasperowicz A, Stan-Glasek K, Guczynska W, Piknová M, Pristas P, Nigutová K, et al. Fructanolytic and saccharolytic enzymes of the rumen bacterium Pseudobutyrivibrio ruminis strain 3–preliminary study. Folia Microbiol. 2010;55(4):329–31.

    Article  CAS  Google Scholar 

  34. Hippe H, Hagelstein A, Kramer I, Swiderski J, Stackebrandt E. Phylogenetic analysis of Formivibrio citricus, Propionivibrio dicarboxylicus, Anaerobiospirillum thomasii. Succinimonas amylolytica and Succinivibrio dextrinosolvens and proposal of Succinivibrionaceae fam. nov. Int J Syst Bacteriol. 1999;49(Pt 2):779–82.

    Article  Google Scholar 

  35. Suen G, Weimer PJ, Stevenson DM, Aylward FO, Boyum J, Deneke J, et al. The complete genome sequence of Fibrobacter succinogenes S85 reveals a cellulolytic and metabolic specialist. PLoS One. 2011;6(4):e18814.

    Article  CAS  Google Scholar 

  36. Deublein D, Steinhauser A. Biogas from Waste And Renewable Resources: An Introduction. Hoboken, New Jersey, USA: John Wiley & Sons; 2010.

    Book  Google Scholar 

  37. Purushe J, Fouts DE, Morrison M, White BA, Mackie RI, North American Consortium for Rumen Bacteria, et al. Comparative genome analysis of Prevotella ruminicola and Prevotella bryantii: insights into their environmental niche. Microb Ecol. 2010;60(4):721–9.

    Article  Google Scholar 

  38. Palatsi J, Viñas M, Guivernau M, Fernandez B, Flotats X. Anaerobic digestion of slaughterhouse waste: main process limitations and microbial community interactions. Bioresour Technol. 2011;102(3):2219–27.

    Article  CAS  Google Scholar 

  39. Baena S, Fardeau ML, Labat M, Ollivier B, Thomas P, Garcia JL, et al. Aminobacterium colombiensegen nov. sp. nov., an amino acid-degrading anaerobe isolated from anaerobic sludge. Anaerobe. 1998;4(5):241–50.

    Article  CAS  Google Scholar 

  40. Chertkov O, Sikorski J, Brambilla E, Lapidus A, Copeland A, Glavina Del Rio T, et al. Complete genome sequence of Aminobacterium colombiense type strain (ALA-1). Stand Genomic Sci. 2010;2(3):280–9.

    Article  Google Scholar 

  41. Pelletier E, Kreimeyer A, Bocs S, Rouy Z, Gyapay G, Chouari R, et al. “Candidatus Cloacamonas acidaminovorans”: genome sequence reconstruction provides a first glimpse of a new bacterial division. J Bacteriol. 2008;190(7):2572–9.

    Article  CAS  Google Scholar 

  42. Kröber M, Bekel T, Diaz NN, Goesmann A, Jaenicke S, Krause L, et al. 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(1):38–49.

    Article  Google Scholar 

  43. Maus I, Wibberg D, Stantscheff R, Eikmeyer FG, Seffner A, Boelter J, et al. Complete genome sequence of the hydrogenotrophic, methanogenic archaeon Methanoculleus bourgensis strain MS2(T), isolated from a sewage sludge digester. J Bacteriol. 2012;194(19):5487–8.

    Article  CAS  Google Scholar 

  44. Walter A, Knapp BA, Farbmacher T, Ebner C, Insam H, Franke-Whittle IH. Searching for links in the biotic characteristics and abiotic parameters of nine different biogas plants. Microb Biotechnol. 2012;5(6):717–30.

    Article  Google Scholar 

  45. Smith KS, Ingram-Smith C. Methanosaeta, the forgotten methanogen? Trends Microbiol. 2007;15(4):150–5.

    Article  CAS  Google Scholar 

  46. Zakrzewski M, Goesmann A, Jaenicke S, Jünemann S, Eikmeyer FG, Szczepanowski R, et al. Profiling of the metabolically active community from a production-scale biogas plant by means of high-throughput metatranscriptome sequencing. J Biotechnol. 2012;158(4):248–58.

    Article  CAS  Google Scholar 

  47. Dridi B, Fardeau ML, Ollivier B, Raoult D, Drancourt M. Methanomassiliicoccus luminyensis gen. nov., sp. nov., a methanogenic archaeon isolated from human faeces. Int J Syst Evol Microbiol. 2012;62(Pt 8):1902–7.

    Article  CAS  Google Scholar 

  48. Griffiths-Jones S, Bateman A, Marshall M, Khanna A, Eddy SR. Rfam: an RNA family database. Nucleic Acids Res. 2003;31(1):439–41.

    Article  CAS  Google Scholar 

  49. Nettmann E, Bergmann I, Pramschüfer S, Mundt K, Plogsties V, Herrmann C, et al. Polyphasic analyses of methanogenic archaeal communities in agricultural biogas plants. Appl Environ Microbiol. 2010;76(8):2540–8.

    Article  CAS  Google Scholar 

  50. Kitamura K, Fujita T, Akada S, Tonouchi A. Methanobacterium kanagiense sp. nov., a hydrogenotrophic methanogen, isolated from rice-field soil. Int J Syst Evol Microbiol. 2011;61(Pt 6):1246–52.

    Article  Google Scholar 

  51. Rea S, Bowman JP, Popovski S, Pimm C, Wright AD. Methanobrevibacter millerae sp. nov. and Methanobrevibacter olleyae sp. nov., methanogens from the ovine and bovine rumen that can utilize formate for growth. Int J Syst Evol Microbiol. 2007;57(Pt 3):450–6.

    Article  CAS  Google Scholar 

  52. Hu P, Yang M, Zhang A, Wu J, Chen B, Hua Y, et al. Complete genome sequence of Streptococcus suis serotype 3 strain ST3. J Bacteriol. 2011;193(13):3428–9.

    Article  CAS  Google Scholar 

  53. Feinberg L, Foden J, Barrett T, Davenport KW, Bruce D, Detter C, et al. Complete genome sequence of the cellulolytic thermophile Clostridium thermocellum DSM1313. J Bacteriol. 2011;193(11):2906–7.

    Article  CAS  Google Scholar 

  54. Rusniok C, Couvé E, Da Cunha V, El Gana R, Zidane N, Bouchier C, et al. Genome sequence of Streptococcus gallolyticus: insights into its adaptation to the bovine rumen and its ability to cause endocarditis. J Bacteriol. 2010;192(8):2266–76.

    Article  CAS  Google Scholar 

  55. Heinl S, Wibberg D, Eikmeyer F, Szczepanowski R, Blom J, Linke B, et al. Insights into the completely annotated genome of Lactobacillus buchneri CD034, a strain isolated from stable grass silage. J Biotechnol. 2012;161(2):153–66.

    Article  CAS  Google Scholar 

  56. Lazarev VN, Levitskii SA, Basovskii YI, Chukin MM, Akopian TA, Vereshchagin VV, et al. Complete genome and proteome of Acholeplasma laidlawii. J Bacteriol. 2011;193(18):4943–53.

    Article  CAS  Google Scholar 

  57. Sikorski J, Teshima H, Nolan M, Lucas S, Hammon N, Deshpande S, et al. Complete genome sequence of Mahella australiensis type strain (50–1 BON). Stand Genomic Sci. 2011;4(3):331–41.

    Article  CAS  Google Scholar 

  58. Anderson I, Ulrich LE, Lupa B, Susanti D, Porat I, Hooper SD, et al. Genomic characterization of Methanomicrobiales reveals three classes of methanogens. PLoS One. 2009;4(6):e5797.

    Article  Google Scholar 

  59. Anderson IJ, Sieprawska-Lupa M, Lapidus A, Nolan M, Copeland A, Glavina Del Rio T, et al. Complete genome sequence of Methanoculleus marisnigri Romesser et al. 1981 type strain JR1. Stand Genomic Sci. 2009;1(2):189–96.

    Article  Google Scholar 

  60. Parshina SN, Sipma J, Nakashimada Y, Henstra AM, Smidt H, Lysenko AM, et al. Desulfotomaculum carboxydivorans sp. nov., a novel sulfate-reducing bacterium capable of growth at 100% CO. Int J Syst Evol Microbiol. 2005;55(Pt 5):2159–65.

    Article  CAS  Google Scholar 

  61. He M, Sebaihia M, Lawley TD, Stabler RA, Dawson LF, Martin MJ, et al. Evolutionary dynamics of Clostridium difficile over short and long time scales. Proc Natl Acad Sci U S A. 2010;107(16):7527–32.

    Article  CAS  Google Scholar 

  62. Chovatia M, Sikorski J, Schröder M, Lapidus A, Nolan M, Tice H, et al. Complete genome sequence of Thermanaerovibrio acidaminovorans type strain (Su883). Stand Genomic Sci. 2009;1(3):254–61.

    Article  Google Scholar 

  63. Sieber JR, Sims DR, Han C, Kim E, Lykidis A, Lapidus AL, et al. The genome of Syntrophomonas wolfei: new insights into syntrophic metabolism and biohydrogen production. Environ Microbiol. 2010;12(8):2289–301.

    CAS  Google Scholar 

  64. Seedorf H, Fricke WF, Veith B, Brüggemann H, Liesegang H, Strittmatter A, et al. The genome of Clostridium kluyveri, a strict anaerobe with unique metabolic features. Proc Natl Acad Sci U S A. 2008;105(6):2128–33.

    Article  CAS  Google Scholar 

  65. Hemme CL, Mouttaki H, Lee YJ, Zhang G, Goodwin L, Lucas S, et al. Sequencing of multiple clostridial genomes related to biomass conversion and biofuel production. J Bacteriol. 2010;192(24):6494–6.

    Article  CAS  Google Scholar 

  66. Hinse D, Vollmer T, Rückert C, Blom J, Kalinowski J, Knabbe C, et al. Complete genome and comparative analysis of Streptococcus gallolyticus subsp. gallolyticus, an emerging pathogen of infective endocarditis. BMC Genomics. 2011;12:400.

    Article  CAS  Google Scholar 

  67. Schlegel L, Grimont F, Ageron E, Grimont PA, Bouvet A. Reappraisal of the taxonomy of the Streptococcus bovis/Streptococcus equinus complex and related species: description of Streptococcus gallolyticus subsp. gallolyticus subsp. nov., S. gallolyticus subsp. macedonicus subsp. nov. and S. gallolyticus subsp. pasteurianus subsp. nov. Int J Syst Evol Microbiol. 2003;53(Pt 3):631–45.

    Article  CAS  Google Scholar 

  68. Papadimitriou K, Ferreira S, Papandreou NC, Mavrogonatou E, Supply P, Pot B, et al. Complete genome sequence of the dairy isolate Streptococcus macedonicus ACA-DC 198. J Bacteriol. 2012;194(7):1838–9.

    Article  CAS  Google Scholar 

  69. Feng X, Xu Y, Chen Y, Tang YJ. MicrobesFlux: a web platform for drafting metabolic models from the KEGG database. BMC Syst Biol. 2012;6:94.

    Article  Google Scholar 

  70. Lee YE, Jain MK, Lee C, Zeikus JG. Taxonomic distinction of saccharolytic thermophilic anaerobes: description of Thermoanaerobacterium xylanolyticum gen. nov., sp. nov., and Thermoanaerobacterium saccharolyticum gen. nov., sp. nov.; reclassification of Thermoanaerobium brockii, Clostridium thermosulfurogenes, and Clostridium thermohydrosulfuricum E100-69 as Thermoanaerobacter brockii comb. nov., Thermoanaerobacterium thermosulfurigenes comb. nov., and Thermoanaerobacter thermohydrosulfuricus comb. nov., respectively; and transfer of Clostridium thermohydrosulfuricum 39E to Thermoanaerobacter ethanolicus. Int J Syst Bacteriol. 1993;43:41–51.

    Article  Google Scholar 

  71. Kosaka T, Uchiyama T, Ishii S, Enoki M, Imachi H, Kamagata Y, et al. Reconstruction and regulation of the central catabolic pathway in the thermophilic propionate-oxidizing syntroph Pelotomaculum thermopropionicum. J Bacteriol. 2006;188(1):202–10.

    Article  CAS  Google Scholar 

  72. Kosaka T, Kato S, Shimoyama T, Ishii S, Abe T, Watanabe K. The genome of Pelotomaculum thermopropionicum reveals niche-associated evolution in anaerobic microbiota. Genome Res. 2008;18(3):442–8.

    Article  CAS  Google Scholar 

  73. Tamaru Y, Miyake H, Kuroda K, Nakanishi A, Matsushima C, Doi RH, et al. Comparison of the mesophilic cellulosome-producing Clostridium cellulovorans genome with other cellulosome-related clostridial genomes. Microb Biotechnol. 2011;4(1):64–73.

    Article  CAS  Google Scholar 

  74. Lien T, Madsen M, Rainey FA, Birkeland NK. Petrotoga mobilis sp. nov., from a North Sea oil-production well. Int J Syst Bacteriol. 1998;48(Pt 3):1007–13.

    Article  CAS  Google Scholar 

  75. Su F, Tao F, Tang H, Xu P. Genome sequence of the thermophile Bacillus coagulans Hammer, the type strain of the species. J Bacteriol. 2012;194(22):6294–5.

    Article  CAS  Google Scholar 

  76. Djao OD, Zhang X, Lucas S, Lapidus A, Del Rio TG, Nolan M, et al. Complete genome sequence of Syntrophothermus lipocalidus type strain (TGB-C1). Stand Genomic Sci. 2010;3(3):268–75.

    Article  Google Scholar 

  77. Goberna M, Insam H, Franke-Whittle IH. Effect of biowaste sludge maturation on the diversity of thermophilic Bacteria and Archaea in an anaerobic reactor. Appl Environ Microbiol. 2009;75(8):2566–72.

    Article  CAS  Google Scholar 

  78. Schnürer J, Olsson J, Börjesson T. Fungal volatiles as indicators of food and feeds spoilage. Fungal Genet Biol. 1999;27(2–3):209–17.

    Article  Google Scholar 

  79. Lin IH, Liu TT, Teng YT, Wu HL, Liu YM, Wu KM, et al. Sequencing and comparative genome analysis of two pathogenic Streptococcus gallolyticus subspecies: genome plasticity, adaptation and virulence. PLoS One. 2011;6(5):e20519.

    Article  CAS  Google Scholar 

  80. Zheng X, Zheng H, Lan R, Ye C, Wang Y, Zhang J, et al. Identification of genes and genomic islands correlated with high pathogenicity in Streptococcus suis using whole genome tiling microarrays. PLoS One. 2011;6(3):e17987.

    Article  CAS  Google Scholar 

  81. Wibberg D, Blom J, Jaenicke S, Kollin F, Rupp O, Scharf B, et al. Complete genome sequencing of Agrobacterium sp. H13-3, the former Rhizobium lupini H13-3, reveals a tripartite genome consisting of a circular and a linear chromosome and an accessory plasmid but lacking a tumor-inducing Ti-plasmid. J Biotechnol. 2011;155(1):50–62.

    Article  CAS  Google Scholar 

  82. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.

    Article  CAS  Google Scholar 

  83. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.

    Article  CAS  Google Scholar 

  84. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2010;27(16):2194–200.

    Article  Google Scholar 

  85. 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(16):5261–7.

    Article  CAS  Google Scholar 

  86. Nawrocki EP, Kolbe DL, Eddy SR. Infernal 1.0: inference of RNA alignments. Bioinformatics. 2009;25(10):1335–7.

    Article  CAS  Google Scholar 

  87. Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26(7):1641–50.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

This research was supported by the Bundenministerium für Ernährung, Landwirtschaft und Verbraucherschutz (BMEL) and Fachagentur für nachwachsenende Rohstoffe (FNR) projects, grant numbers 22006712 and 22028711. FE and IM acknowledge the receipt of a scholarship from the CLIB Graduate Industrial Biotechnology. The Bioinformatics Resource Facility (BRF) is thanked for their support on running the bioinformatics software at the Center for Biotechnology (CeBiTec), Bielefeld, Germany. Many thanks to the sequencing department at the Center of Biotechnology (CeBiTec).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Andreas Schlüter.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

YS participated in functional sequence data analysis and coordinated and drafted the manuscript. MZ participated in sequence data processing, carried out the taxonomic analysis and participated in the functional analysis of the sequencing data. IM participated in manuscript preparation. FE did the laboratory work and participated in fragment recruitment analysis and presentation. SJ participated in sequence data processing and analysis. NR provided biogas plant and metadata information. AP supervised the study and contributed to the draft of the manuscript. AS conceived the study, participated in its design and the fragment recruitment analysis. All authors read and approved the final manuscript.

Yvonne Stolze and Martha Zakrzewski contributed equally to this work.

Additional files

Additional file 1: Figure S1.

Rarefaction analyses of sequenced metagenomes and 16S rRNA gene amplicons originating from dry (BGP_DF) and wet fermentation biogas plants (BGP_WF) microbial communities. Rarefaction analysis plots on (A) taxonomic mean richness at the family rank derived from metagenome data, (B) environmental gene tags (EGT) derived from metagenome data and (C) operational taxonomic unit (OTU) richness derived from 16S rRNA gene amplicons in correlation with the sampled reads.

Additional file 2: Table S1.

Unique environmental gene tags (EGTs) encoding different proteins found in metagenome datasets for the wet and dry fermentation processes. Only EGTs featuring an abundance of 0 in the one and at least five in the other metagenome dataset were taken into account. Relative abundances were normalized based on the smallest dataset.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Stolze, Y., Zakrzewski, M., Maus, I. et al. Comparative metagenomics of biogas-producing microbial communities from production-scale biogas plants operating under wet or dry fermentation conditions. Biotechnol Biofuels 8, 14 (2015). https://doi.org/10.1186/s13068-014-0193-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-014-0193-8

Keywords