- Open Access
Comparative analysis of deep sequenced methanogenic communities: identification of microorganisms responsible for methane production
© The Author(s) 2018
- Received: 16 October 2018
- Accepted: 11 December 2018
- Published: 20 December 2018
Although interactions between microorganisms involved in biogas production are largely uncharted, it is commonly accepted that methanogenic Archaea are essential for the process. Methanogens thrive in various environments, but the most extensively studied communities come from biogas plants. In this study, we employed a metagenomic analysis of deeply sequenced methanogenic communities, which allowed for comparison of taxonomic and functional diversity as well as identification of microorganisms directly involved in various stages of methanogenesis pathways.
A comprehensive metagenomic approach was used to compare seven environmental communities, originating from an agricultural biogas plant, cattle-associated samples, a lowland bog, sewage sludge from a wastewater treatment plant and sediments from an ancient gold mine. In addition to the native consortia, two laboratory communities cultivated on maize silage as the sole substrate were also analyzed. Results showed that all anaerobic communities harbored genes of all known methanogenesis pathways, but their abundance varied greatly between environments and that genes were encoded by different methanogens. Identification of microorganisms directly involved in different stages of methane production revealed that hydrogenotrophic methanogens, such as Methanoculleus, Methanobacterium, Methanobrevibacter, Methanocorpusculum or Methanoregula, predominated in most native communities, whereas acetoclastic Methanosaeta seemed to be the key methanogen in the wastewater treatment plant. Furthermore, in many environments, the methylotrophic pathway carried out by representatives of Methanomassiliicoccales, such as Candidatus Methanomethylophilus and Candidatus Methanoplasma, seemed to play an important role in methane production. In contrast, in stable laboratory reactors substrate versatile Methanosarcina predominated.
The metagenomic approach presented in this study allowed for deep exploration and comparison of nine environments in which methane production occurs. Different abundance of methanogenesis-related functions was observed and the functions were analyzed in the phylogenetic context in order to identify microbes directly involved in methane production. In addition, a comparison of two metagenomic analytical tools, MG-RAST and MetAnnotate, revealed that combination of both allows for a precise characterization of methanogenic communities.
- Methanogenesis pathways
- Environmental communities
- Deep sequencing
Biogas is one of the most promising solutions for energy production associated with degradation of various types of agricultural and industrial wastes, such as food waste, animal manure, crops and wastewater sludge. Under natural conditions, methane emission occurs in diverse anaerobic environments, such as animal digestive tracts, peatlands, anaerobic digesters, wetlands, rice field soils, marine sediments and hydrothermal habitats .
Many microorganisms responsible for the methanogenesis process have been identified and described so far. Most known biological producers of methane are represented within the Euryarcheota phylum: orders Methanobacteriales, Methanocellales, Methanococcales, Methanomicrobiales, Methanopyrales, Methanosarcinales and Methanomassiliicoccales [3, 4]. Most of methanogens carry out reduction of carbon dioxide connected with the consumption of hydrogen. However, within the orders of Methanosarcinales, Methanobacteriales and Methanomassiliicoccales there are microorganisms with the ability to convert methylated compounds into methane, while members of Methanosarcinales are capable of utilizing also acetate. Till now, only Methanosarcina of Methanosarcinales has been considered capable of carrying out all three methanogenesis pathways [4–7].
Development of metagenome sequencing allowed much better understanding of microorganisms’ role in anaerobic digestion and methane production by the analysis of both taxonomic and functional diversity [8–10]. In recent past, methanogenesis potential was often concluded from the abundance of a given methanogen only. However, genome-centric studies showed that conclusions drawn only from a taxonomic classification can lead to an underestimation of the true methanogenesis potential. For example, Rotaru and colleagues revealed that Methanosaeta, a genus assumed to be strictly acetoclastic, had a complete set of genes necessary for the reduction of carbon dioxide to methane via the hydrogenotrophic pathway .
Recent studies have suggested that we are only beginning to understand the diversity of methanogens and the methanogenesis process itself. The best studied methane-producing communities come from biogas reactors (e.g. [8–10]), but many waste materials used as a supplement in biogas systems could provide novel microorganisms adapted to the production of methane and the degradation of organic matter [12–15]. It seems right to assume that environmental communities are much more diverse but poorly characterized. This view has been reflected by the latest description of novel hydrogen-dependent methylotrophy in the Methanomassiliicoccales order  or by an identification of distant homologues of methanogenesis genes in Bathyarchaeota, which are still not considered methanogens . Moreover there are recent reports demonstrating biogenic methane production by Cyanobacteria within the Bacteria domain which challenge the paradigm that methanogenesis is exclusive to Archaea .
The aim of this study was to characterize and compare the content of methanogenesis-related genes for very different anaerobic microbial communities which operates on diverse substrates and coming from methanogenic environments, such as agricultural biogas reactors (maize silage), a wastewater treatment plant (industrial and municipal wastes), cattle-associated habitats (grass and grains), a lowland bog (peat moss) and an ancient gold mine (decaying wooden rafters and other abandoned organic materials). With the use of widely available metagenomic tools, we performed an in-depth analysis of the methanogenesis process. Special focus was given to the identification of microorganisms directly involved in methane production in native communities by the phylogenetic classification of methanogenesis-related genes, especially for methyl-CoM reductase (mcr) that participates in the final step of methane release, as well as for formylmethanofuran dehydrogenase (fmd); CO dehydrogenase/acetyl-CoA synthase (cdh), methanol-specific methyltransferase (mta) and methylamine-specific methyltransferases (mtm, mtb, mtt) as they are responsible for triggering carbon dioxide (CO2), acetate (CH3COOH), methanol (CH3OH) and methylamines ((CH3)xNH)) utilization into methane, respectively [2, 4, 5, 8]. In addition to the deep-sequenced metagenome analysis, the methanogenesis potential was also assessed by a simple and commonly used method of semi-quantitative marker gene amplification .
This methodology was also applied to explore the taxonomic and functional community change after laboratory cultivation. The first laboratory community, originated from a fermenter tank of a biogas plant, was considered a model sample with an ability to produce biogas with methane content above 50% . The second laboratory community came from raw sewage sludge from a wastewater treatment plant, as the sewage sludge is often used in biogas plants as a co-substrate source. In this type of samples, residual methanogenesis potential is usually low and methane content does not exceed 5% .
General description of sequenced metagenomes and bioinformatic strategy
A metagenomic approach was used to characterize different anaerobic communities, and to assess their ability for methane production. Among the native communities, two originated from full-scale biogas reactors. Samples were collected from an agricultural biogas fermenter (ABF) and an agricultural biogas hydrolyzer (ABH) and were treated as a reference source of an efficient methanogenic community residing in a two-stage biogas plant. Aside from the extensively studied biogas reactor communities, we examined also more natural communities from environments where methane emission is detected. These other communities came from cattle manure (CM), cattle slurry (CS) and sewage sludge from a wastewater treatment plant (WTP), as biogas plants are often supplemented with these waste materials. The last two environmental samples originated from a lowland bog (LB) and an ancient gold mine (GM), as these habitats are a reservoir of diverse uncultivated microorganisms potentially beneficial to the methanogenesis process [21–23]. For comparative analyses of microbial community structural and functional change following a cultivation process, two laboratory consortia (ABF_TS and WTP_TS) were also implemented in the study. They were selected in batch reactors, and then stabilized in two-stage (TS) bioreactors.
General analyses of community taxonomic and functional profiles were carried out through identification of similar sequences in reference databases with the use of the MG-RAST pipeline . To obtain a more detailed insight into the methanogenesis process and the microbes involved, we used MetAnnotate , in which a HMM-based search is implemented and taxonomy is determined by the best hit or phylogenetic tree placement of obtained hits. In this study, we compared results from both approaches.
Among the investigated metagenomes, samples coming from agricultural biogas plant (ABF, ABH) were considered model consortia that harbor microorganisms essential for all steps of high-performance biogas production. The most specific feature of the agricultural biogas plant fermenter sample (ABF) was the high abundance of methanogenic archaeons, i.e. Methanoculleus (4.8%) and Methanosarcina (2.3%), whereas in other native metagenomes their abundance was significantly lower, close to 0.7% (ABH and LB samples) or even not exceeding 0.1% (CM, CS, GM and WTP samples (see Additional file 1: Table S3). Furthermore, in the agricultural biogas plant samples, we detected 2.6% (ABF) and 3.5% (ABH) of metagenomic sequences classified as Candidatus Cloacamonas, while in other native metagenomes sequence counts of this bacterium accounted for a maximum of 0.2%. Moreover, the ABH sample was highly enriched in polysaccharide-degrading Prevotella (28.5%), Lactobacillus (4.1%) and numerous genera belonging to Veillonellaceae family (2.2% cumulatively).
Interestingly, in the cattle-associated samples we observed that Pseudomonas (7.5%; 6.7%; CM and CS, respectively), Acinetobacter (1.5%; 4.2%) and Bacillus (2.4%; 2.2%) were more abundant than in bioreactor samples. In addition, CM was richer in members of Acholeplasma (2.4%). In the case of the community originated from a wastewater treatment plant (WTP), we observed a high number of sequences assigned to Acinetobacter (9.6%), Arcobacter (6.1%), Aeromonas (4.6%) and Acidovorax (3.3%). The last two environmental samples analyzed, originating from a gold mine (GM) and a lowland bog (LB), were the most diverse and had few common genera with each other and with the remaining metagenomes (Fig. 3; Additional file 1: Fig. S1). The most unique feature of GM sample was a very high content of sequences assigned to methanotrophic Methylobacter (11.3%), Methylococcus (2.5%) and Methylotenera (2.3%). As a comparison, in other metagenomes these bacteria accounted for a maximum of 0.5%. In the case of LB community, the highest abundance was observed for Geobacter (3.3%) and Candidatus Solibacter (2.4%) (Fig. 3b). Importantly, LB community was enriched in numerous, low abundance environmental Archaea, including both methanogenic Archaea (3.9%) and other Archaea (4.7%) without a confirmed methane production ability (see Tax MG-RAST, Additional file 2).
Both laboratory consortia, originated from an agricultural biogas fermenter (ABF_TS) and sewage sludge from a wastewater treatment plant (WTP_TS), were dominated by Bacteria (93.0% in ABF_TS; 94.6% in WTP_TS), followed by Archaea (6.1% in ABF_TS, 4.9% in WTP_TS). Taxonomic annotations of the ABF_TS sample revealed that despite laboratory cultivation, the structure of microorganisms has been largely preserved (Fig. 3b; Additional file 1: Fig. S1). Nevertheless, we observed at least 1% change of the relative abundance of the following dominant genera: Bacteroides (from 10.9 to 14.1%), Clostridium (from 6.0 to 8.5%) and Parabacteroides (from 2.2 to 3.3%), Pseudomonas (from 0.3 to 1.3%), Prevotella (from 4.4 to 2.7%) and Methanoculleus (from 4.8 to 0.5%) (Fig. 3b; Additional file 1: Table S3). In the case of the consortium selected from sewage sludge of a wastewater treatment plant (WTP_TS), a substantial change of the community structure was observed compared to the native WTP sample (Fig. 3b; Additional file 1: Fig. S1). In general, abundance of Proteobacteria representatives diminished from 56.7 to 11.7%, while abundance of Bacteroidetes and Firmicutes increased from 26.8 and 10.0% in inocula to 34.6 and 33.8% in the laboratory consortium, respectively. Moreover, Archaea abundance increased from 0.5 to 4.9% of the total microbial structure. It seems that the most important change of laboratory consortium in comparison to its inoculum was the increase of fermentative Bacteroides (from 11.5 to 19.5%), Clostridium (from 2.1 to 10.7%), Parabacteroides (from 2.0 to 3.3%), Eubacterium (from 0.8 to 2.1%), Candidatus Cloacamonas (from 0.1 to 3.1%), Ruminococcus (from 0.4 to 1.6%), Syntrophomonas (from 0.1 to 1.6%) and methanogenic Methanosarcina (from 0.1% to 3.1%) (see Additional file 1: Table S3). Simultaneously abundance of other genera such as Acinetobacter, Arcobacter, Aeromonas, Acidovorax, Streptococcus, Flavobacterium, Thauera and Tolumonas decreased from 9.6% in native community (WTP) to max. level of 0.6% in laboratory community (WTP_TS) (Fig. 3b; Additional file 1: Table S3).
To explore the metabolic potential of the studied communities, we performed a detailed analysis of metagenomic sequences annotated against SEED subsystems within the MG-RAST pipeline. We detected on average 5590 functional categories with at least 0.001% annotated reads (see Additional file 1: Table S2 and Fun MG-RAST, Additional file 2). Analysis of commonly used diversity and evenness indices showed that all metagenomes were quite similar in terms of diversity. The most functionally diverse and even was WTP community (Shannon–Wiener index—7.560, Pielou index 0.834) while the least one was ABF metagenome (Shannon–Wiener index—7.196, Pielou index 0.816) (see Additional file 1: Table S5). Furthermore, similarly to RefSeq Bray–Curtis distances calculation, some samples located near each other e.g. ABF and ABF_TS, while GM and LB samples were the most different from majority of the analyzed metagenomes (see Additional file 1: Fig. S2).
Considering the percentage of methanogenesis-related annotations, ABF community followed by LB, ABF_TS and WTP_TS and ABH consortia had the highest potential for methane production via different pathways. Communities from an agricultural biogas plant and laboratory reactors (ABF, ABH, ABF_TS, WTP_TS) had similar sequence profiles of the selected genes. However, the functional profile of lowland bog (LB) metagenome differed significantly (Fig. 4). In the LB sample, we observed an overrepresentation of acs, cdh and hdr genes compared to other methanogenesis genes, including the key mcr genes. An even greater overrepresentation of genes involved in acetate utilization was observed for cattle-associated (CM, CS) and wastewater treatment plant (WTP) samples. However, in these samples most of the selected sequences were classified to ack and pta genes but few to cdh genes that directly mediate the transfer of carbon to methane pathway (Fig. 4). In the case of the native consortium from a gold mine (GM), we observed a relatively high number of sequences of fmd, ftr, mch, hmd, ack, pta, acs genes involved in the utilization of carbon dioxide and acetate, and very few sequences of mcr genes encoding the key methyl-CoM reductase enzyme.
The comparison of laboratory consortia and their native counterparts indicated that during laboratory cultivation the cumulative abundance of methanogenesis genes increased for WTP_TS and decreased for ABF_TS. In the case of WTP_TS, the increase concerned all genes except ack, pta and acs, while for ABF_TS the reduction concerned all studied genes.
Methanogenesis genes-specific phylogenetic characterization
A sequence-based analysis of metagenomic sequences with the use of the MG-RAST pipeline gives an opportunity to explore metabolic potential of complex communities. Nevertheless, metagenomic data enables also a simultaneous identification of function and microorganisms responsible for specific processes. In the following part of our study, we present detailed results of phylogenetic placement of methanogenesis-related sequences. For this type of an analysis, the HMM search and the taxonomic classification approach implemented in MetAnnotate  were used. We focused on genes encoding methyl CoM reductase (mcr) responsible for the final release of methane, as well as on sequences of genes encoding enzymes which are considered crucial in the utilization of various substrates, such as: CO2—formylmethanofuran dehydrogenase (fmd); CH3COOH—CO dehydrogenase/acetyl-CoA synthase (cdh); CH3OH—methanol methyltransferase (mta); (CH3)xNH—methylamine methyltransferase (mtm, mtb, mtt). A summary of key microorganisms involved in specific steps of the methanogenesis superpathway for each analyzed sample is presented in Additional file 1: Fig. S3A–I.
Apart from the identification of microorganisms responsible for the last step of methanogenesis, we aimed for identifying those involved in carbon incorporation from different substrates, such as CO2, CH3COOH, CH3OH and (CH3)xNH. An analysis of formylmethanofuran dehydrogenase subunits ACE (fmdACE), indicators of hydrogenotrophic pathway, showed that among the native communities most hits were detected for ABF consortium, followed by GM, LB, ABH, ABF_TS, WTP_TS, CM, CS and WTP communities (Fig. 5b, bar graph panel). There was an up to 19-fold difference in the fmdACE abundance between the ABF and WTP samples. In the ABF sample, hydrogenotrophic Methanoculleus (66%) was the organism predominantly assigned to fmd genes. In the WTP sample, a majority of fmd sequences mapped to Thauera (21%) and to numerous taxa of Archaea and Bacteria, which did not exceed 6% (Fig. 5b). Following the laboratory cultivation of the aforementioned communities, Methanosarcina predominated in ABF_TS (28%) and WTP_TS (54%) fmd sequences. In comparison, for LB metagenome many of the fmd sequences mapped to various low-abundant Archaea and Bacteria with the highest counts for Desulfobacter (7%) and Methylobacterium (7%). Interestingly, in the case of GM metagenome, most fmd sequences were classified to methanotrophic bacteria, such as Methylobacter, Methyloglobulus, Methylovulum (12% each), but few were classified to methanogenic Archaea (Fig. 5b). In contrast, agricultural biogas hydrolyzer (ABH), cattle manure (CM) and cattle slurry (CS) were dominated mainly by hydrogenotrophic methanogens, such as Methanoculleus (35%) and Methanobacterium (32%) in ABH sample, Methanocorpusculum (23%) in CM sample and Methanobrevibacter (22%) in CS sample. Furthermore, both CM and CS had many hits (~ 10%) for Methanospirillum, Methanoregula, Anaerosalibacter and Blautia (Fig. 5b).
In order to identify microorganisms involved in the acetoclastic pathway, we analyzed phylogenetic assignments of the D subunit of CO dehydrogenase/acetyl-CoA synthase (cdhD), as it is directly involved in the transmission of a methyl group from acetate during acetoclastic methanogenesis [4, 5]. The abundance of cdhD was highest for ABH, followed by LB, ABF, ABF_TS, WTP_TS, CM, CS samples, while in GM and WTP metagenomes it was scarce (Fig. 5c, bar graph panel). Identification of phylogenetic matches of cdhD indicated that in most of the metagenomes acetate utilization was mediated by two to four dominant genera with percentage above 10%. In the case of laboratory communities, cdhD sequences were mostly assigned to one methanogen, namely Methanosarcina, which accounted for 65% (ABF_TS) and 94% (WTP_TS). In the case of the other two samples, approximately half of the annotations were classified to one methanogen, namely Methanobacterium (53%) in ABH and Methanoregula (49%) in LB (Fig. 5c). For other metagenomes, cdhD sequences were assigned to few dominant genera. These were Methanosarcina (32%), Methanoculleus (24%), Methanosaeta (17%), and Methanobacterium (12%) in the ABF metagenome; Methanoregula (24%) and Desulfocapsa (15%) in GM; Methanosaeta (31%) and Methanoregula (23%) in WTP and Methanosarcina (24%, 32%), Blautia (25%, 15%), Treponema (13%, 14%) in CM and CS samples, respectively (Fig. 5c). Furthermore, in the CM, CS, GM, LB, WTP samples, cdhD was often encoded by various Archaea and Bacteria.
Amplicon-based analysis of methanogenesis-related genes
In order to test if methanogenesis potential of deeply sequenced communities could be grasped by widely used and simple amplicon-based approach, we performed a semi-quantitative amplification of marker genes designed to monitor the methanogenesis potential of environmental samples . PCR-amplified genes included alpha, beta and gamma subunits of methyl-CoM reductase (mcrABG), methanol-specific methyltransferase (mtaB) and methylamine-specific methyltransferase (mtbA).
Metagenome analysis is a very useful approach for a comprehensive description of complex microbial communities. With various tools, a different level of insight into community performance may be gained. Here, we applied deep shotgun metagenomic sequencing coupled with two commonly used metagenomic analytical tools, MG-RAST and MetAnnotate, to describe and compare nine samples coming from the environments where methane production was detected. A more general description of taxonomic and functional structure, achievable with MG-RAST, was complemented with detailed analyses, including phylogenetic placement of methanogenesis-related genes, using MetAnnotate. Our approach enabled us also to compare applicability of both tools for metagenomic analysis of methanogenic environments.
MG-RAST-based diversity analysis, presented in this study, demonstrated that organic matter degraders, such as Bacteroides, Clostridium, Parabacteroides, Prevotella, Candidatus Cloacamonas, Candidatus Solibacter or Geobacter, predominated in native consortia. These microorganisms were previously detected in other methanogenic communities [26–35]. Both our results and literature data [36–39] suggest that high abundance of Bacteroides could increase overall carbohydrate hydrolytic performance. Despite an important role of organic matter degraders, the key players in biogas systems are methanogenic Archaea, which often constitute only a small fraction of a microbial community. Among the communities analyzed in this study, ABF was the most enriched in Archaea, followed by LB, ABF_TS, WTP_TS, ABH, CS, CM, WTP and GM (see Additional file 1: Table S3). Archaea abundance corresponded well with the proportions of functional annotations related to methanogenesis in different samples (Fig. 4). However, the GM metagenome had a relatively high abundance of genes of the hydrogenotrophic pathway despite a low abundance of Archaea.
Metagenomic analysis with the MG-RAST pipeline offered an insight into the Archaea community structure and the abundance of genes involved in methanogenesis. The identity of microbes involved in pathways cannot be, however, easily inferred without the use of additional tools. With such general approach (separate taxonomic and functional analyses), it is difficult to determine interactions between microorganisms involved in a given pathway, partially due to the dispersion of gene abundance in various low abundant microorganisms. In contrast, the second tool used in this study, MetAnnotate, allows a more detailed view and the identification of a specific function with the simultaneous assignment to a taxonomic group. In most cases in this study, results obtained with both tools were consistent. Microorganisms that were among the most abundant archaeons identified by MG-RAST dominated also among microorganisms with methanogenesis-related functions assigned by MetAnnotate. This result proved that some conclusions on the functional structure of a community can be drawn from a phylogenetic structure only. In the analysis of samples coming from a mesophilic agricultural biogas plant, the model sample of methanogenic community used in this study, both tools indicated that hydrogenotrophic Methanoculleus and substrate-versatile Methanosarcina were the main methanogens involved in methane production in ABF samples (see Additional file 1: Fig. S3A). In contrast, for ABH, the results obtained with different tools were less consistent. The diversity analysis by MG-RAST revealed that Methanoculleus, Methanothermobacter and Methanosarcina are the major Archaea in the analyzed sample (Fig. 3c), suggesting that they could be the major contributors to methane production. However, MetAnnotate showed that most of the methanogenesis-related sequences mapped to Methanobacterium and Methanoculleus (see Additional file 1: Fig. S3C). However, this difference in the annotations probably resulted from a database bias, as MG-RAST did not assign any sequence in any sample to Methanobacterium but rather to other genera within the Methanobacteriaceae family (Additional file 1: Table S1). As Methanoculleus, Methanosarcina and Methanobacterium, and not Methanothermobacter, are often identified in mesophilic biogas reactors [1, 9, 40–42], the results obtained with MG-RAST should be processed with caution.
MetAnnotate-based phylogenetic classification of functionally annotated sequences proved to be more informative also for communities from less studied environments. In the case of samples from cattle manure (CM) and cattle slurry (CS), both MG-RAST and MetAnnotate indicated that hydrogenotrophic Methanobrevibacter and Methanocorpusculum were the dominant methanogens. However, MetAnnotate showed also that microorganisms from the seventh order of methanogens, such as Candidatus Methanomethylophilus and Candidatus Methanoplasma, substantially contributed to the methylotrophic pathway and to the final release of methane, as they encode mta, mtm, mtb, mtt, mcr genes (see Additional file 1: Fig. S3D, E). Whereas the predominance of Methanobrevibacter in cattle digestive tract is supported by the results of other studies, Methanocorpusculum were rarely detected in rumen [15, 43–45]. Furthermore, the identification of Candidatus Methanomethylophilus and Candidatus Methanoplasma in cattle-derived samples in this study is in agreement with the recent detection of Methanomassiliicoccales in wetlands and gastrointestinal tracts of various animals . The identification of Methanomassiliicoccales representatives in various environments suggests that the methylotrophic pathway may be more important than previously anticipated and can be carried out in a hydrogen-dependent manner .
In the case of very diverse environmental communities, such as samples from a gold mine (GM) or a lowland bog (LB), results derived from the MG-RAST analysis showed that the abundance of methanogenesis genes differs, but it is difficult to infer which microorganisms are involved in methane production, because many microbes were present in an even and low abundance. As a comparison, MetAnnotate offers a direct assignment of microbes to a specific function. Methanogenesis-related sequences of GM and LB communities mapped to multiple genera of both Archaea and Bacteria, suggesting that methane production in GM and LB samples could be negatively affected by substrate competition or the need for an interspecies intermediate transmission. It seems that hydrogenotrophic Methanoregula and Methanoculleus and representatives of the seventh order of methanogens were the major contributors to methane production in the lowland bog (see Additional file 1: Fig. S3G) and in the gold mine (see Additional file 1: Fig. S3F), respectively. However, the classification of many sequences to microorganisms that can use methane as a substrate, such as Archaeoglobus in LB sample and methanotrophic bacteria (e.g. Methylomonas) in GM sample, suggested that in these environments methane could be immediately utilized [47, 48]. Additionally, high abundance of hydrogenothropic-pathway genes may suggest that methylated compounds are rapidly utilized as a source of carbon, rather than used for methane production. This is in agreement with the previous study showing that methylothrophic bacteria can be responsible for carbon assimilation and cycling in a gold mine environment .
The last native community analyzed in this study was isolated from sewage sludge of a wastewater treatment plant. From the diversity analysis performed by the MG-RAST pipeline it could be concluded that Methanosarcina was one of the main methanogens in the WTP sample (Fig. 3c), however a domain profile search for methanogenesis-related genes with the use of MetAnnotate revealed that methane production was mediated mainly by Methanosaeta and Methanomethylovorans (see Additional file 1: Fig. S3H). This is in agreement with other studies [50, 51] and suggests the importance of acetoclastic and methylotrophic methanogenesis in industrial wastewater plants.
Aside from characterization of native communities, we compared the results for laboratory communities ABF_TS and WTP_TS. Both MG-RAST and MetAnnotate analyses suggested that after laboratory cultivation, a functional and taxonomic change occurred both for ABF_TS and WTP_TS. It seems that the selection of microbial community structure could be linked to the operating conditions of laboratory reactors, with the type of substrate as the presumably most important factor. Substrate impact on microbial diversity and dynamics were showed elsewhere . In our study, a drastic reorganization of the microbial structure occurred for WTP_TS community (see Additional file 1: Figs. S1, S2), probably due to organic material change from mixed protein:carbohydrates:lipids in natural environment of wastewater treatment plant to solely carbohydrates in laboratory bioreactors. In contrast, the operating conditions were similar for native ABF and laboratory ABF_TS, however some community profile change was also observed. It seems that, in both ABF_TS and WTP_TS samples, methane production was mediated mainly by Methanosarcina (see Additional file 1: Fig. S3B, I). However, ABF_TS retained the higher biodiversity of methanogens. While the dominance of Methanosarcina in laboratory reactors was previously demonstrated (e.g. [41, 42, 52]), our methodology enabled us to follow the changes by a direct, methanogenesis-specific taxonomic profiling.
Metagenomic approach proved to be effective in characterization of complex microbial communities, yet we also checked whether it is possible to gain insight into methanogenesis potential by a simple PCR-based amplification. The results showed that PCR amplification correctly predicts the presence of methanogenesis genes and thus the ability to produce methane by different methanogenic communities. However, it should be kept in mind that a PCR reaction is highly dependent on primer sequences and that—due to incompleteness of reference databases—not all sequences of environmental microbes could be detected by this method. Nevertheless, PCR amplification seems to be a convenient and simple method for a rapid, preliminary assessment of methanogenesis potential of communities coming particularly from well-studied habitats e.g. industrial reactors. For less studied and more diverse environments we suggest a more complex analysis of deeply sequenced metagenomes, as it offers the opportunity to identify even the low abundant microbes that contribute to methane emission.
It is important to note that due to the ongoing development of public databases and thus incompleteness of environmental microbial sequences there could still be microorganisms involved in methanogenesis process that are omitted in a metagenomic analysis. However, based on our results, metagenomic approach is the least biased method that could provide information about complex microbial communities, which are often very difficult to cultivate. As showed by Campanaro and colleagues , genome-centric metagenomic approach could shed even more light on syntrophic interactions of microorganisms. For a more comprehensive overview of environmental processes such as organic matter degradation and methanogenesis, other meta-omic approaches, such as metatranscriptomic and metaproteomic ones, should be employed alongside metagenomic studies [53–56].
The approach presented in this study allowed to explore in detail complex microbial communities coming from methane-producing environments. Communities predisposed to efficient methane production would be expected to contain a high abundance of genes of different steps of hydrogenotrophic, acetoclastic and methylotrophic pathways, which optimally are encoded by a few microbes. This view held true for engineered environments, such as industrial biogas reactors or laboratory cultures, but in most of the native environments—which have rarely been studied so far—we observed different levels of methanogenesis genes and their dispersion amongst various microorganisms. This was especially apparent for the lowland bog community and could suggest that the less described habitats are reservoirs of little-known microorganisms that contribute to methane cycle. Additionally, in such natural communities, it is more essential to remove just enough intermediate metabolites via methane pathways to keep the community functional than to cover the whole methanogenesis pathway. In contrast, engineered communities are specifically selected for high methane yield capacity.
Still a lot of work has to be done for the comprehensive characterization of methanogenic communities. A general analysis by the MG-RAST pipeline proved to be useful, however, for less described environments with many microorganisms present in low abundance, inferring microbial contributions with MG-RAST could be problematic. A deeper insight into microbial interactions could be obtained by searching for domain sequences with the MetAnnotate pipeline, which links specific metabolic functions to a specific microorganism. Moreover, this phylogenetic assignment of methanogenesis annotations seems to work much better for native communities. As one example, MetAnnotate-based analysis of the methylotrophic pathway suggested that it is carried out by the seventh order of methanogens and may play an important role in methane production in diverse environments.
Environmental sample collection and DNA isolation
Microbial consortium samples analyzed in this study were collected from environments specialized in anaerobic digestion and methane production, such as (I) fermenter (ABF) and (II) hydrolyzer (ABH) tank of an agricultural biogas plant in Miedzyrzec Podlaski, Poland; (III) cattle slurry (CS) and (IV) cattle manure (CM) from a farm in Mikanow, Poland; (V) bottom sediments of effluents from an ancient gold mine (GM) in Zloty Stok, Poland; (VI) peat from a lowland bog (LB) in Otwock, Poland; (VII) raw sewage sludge from a wastewater treatment plant (WTP) “Czajka” in Warsaw, Poland. Semi-liquid samples, such as sludge from the agricultural plant, the wastewater treatment plant and the cattle slurry were collected after removal of the material located in the vicinity of a drain valve by the release of at least 20 L. Likewise for more stable samples, such as cattle manure and peat, an upper surface was excluded and the samples were collected from the depth of approximately 30 cm (manure) and 90 cm (peat). In the case of sample from the gold mine, the material was collected as bottom sediments and surrounding liquids. Where possible (without undesirable aeration), samples were thoroughly mixed. Samples comprised of solids and liquids were maintained under native conditions for a maximum of 16 h prior to DNA extraction. When it was not possible to maintain native conditions, the samples were stored in dry ice.
Isolation of metagenomic DNA was performed according to Dziewit et al. . Briefly, 1 g of a sample biomass (crude sample of CM and LB, or pellet material after centrifugation of ABF, ABH, CS, GM, WTP samples) was resuspended in 2 mL of a lysis buffer (100 mM Tris–HCl (pH 8.0); 100 mM EDTA (pH 8.0); 100 mM Na2HPO4 (pH 8.0); 1.5 M NaCl; 1% (w/v) CTAB). Then the metagenomic DNA was extracted by a five-step bead-beating protocol, combined with freezing and thawing. The isolation of metagenomic DNA was performed in triplicate per sample. The final DNA purification from proteins, humic substances and other compounds was carried out using CsCl density gradient ultracentrifugation. The concentration and quality of the purified DNA were estimated using a NanoDrop 2000c spectrophotometer (NanoDrop Technologies) and by agarose gel electrophoresis.
The isolated whole community DNA (combined triplicates) was used as a template for amplification of methanogenesis markers as described by Dziewit and co-workers . The amplified genes included methyl-CoM reductase (subunits mcrA, mcrB, mcrG1) as well as subunits of methanol-specific methyltransferase (mtaB) and methylamine-specific methyltransferase (mtbA). Primers used were: MLf/MLr; LMCRB/RMCRB; LMCRG1/RMCRG1; LMTAB/RMTAB; LMTBA/RMTAB [18, 57]. Additionally, as a control of polymerase chain reaction (PCR) and purity of each metanogenomic DNA, bacterial and archaeal 16S rDNA fragments were amplified (primers: S-D-Arch-0349-a-S-17/S-D-Arch-0786-a-A-20 and S-D-Bact-0341-b-S-17/S-D-Bact-0785-a-A-21) . All PCR reactions were performed in a TProfessional Thermocycler (Biometra) with Phusion High-Fidelity DNA Polymerase (Thermo Scientific).
Laboratory reactors operation
The cultivation experiment was carried out in lab-scale bioreactors with a working volume of 800 mL, made of 1 L GL 45 glass bottles (Schott Duran, Germany) connected with Dreschel scrubbers and 1 L Tedlar gas bags (Sigma, Germany) as a biogas collector. Batch reactors were inoculated with 10 gvs L−1 of different methanogenic consortia and supplemented with 9.6 gvs L−1 of maize silage. The bioreactors were filled up with spring low-mineral water to the working volume of 800 mL and then pH was adjusted to 7.2 with sodium carbonate. Following our preliminary studies, anaerobic digestion was performed at 37 °C for 21 days for twelve passages. Passaging was carried out every 21 days, with 20% (160 mL) of working volume of the new bioreactors (passages from 2 to 12) coming from the previous passage as an inoculum. Furthermore, batch cultivation was divided into two steps based on substrate input: (I) 9.6 gvs L−1 of maize silage (passages 1–7) and (II) 28.8 gvs L−1 of maize silage (passages 8–12) as the grading of substrate concentrations is a frequently used method for the adaptation of microbial communities. In order to stabilize the best performing consortia, a semi-continuous cultivation in two-stage bioreactors was performed. For this purpose, biomass remains from passages 8–12 were subsampled and further cultivated in batch reactors in order to achieve a sufficient amount of consortia for inoculation of two-stage biogas reactors and to accelerate the start-up phase of the process. Two-stage bioreactor was constructed according to the Polish Patent no. PL197595 . The reactor was equipped with hydraulic agitation and operated in a quasi-continuous mode. The metagenomic DNA representing laboratory consortia was isolated after 30 days of cultivation in two-stage bioreactors. The isolation procedure was identical as described for environmental samples.
In all experiments, the bioreactors were fed with maize silage provided by a farm located in Mikanow, Poland. A bulk amount of maize silage was transported from Mikanow to the laboratory at room temperature, portioned into plastic bags, and stored at 4 °C.
To characterize the physico-chemical profiles of the studied environments, the following parameters were determined: methane content, volatile fatty acids (VFAs) content, total solids (TS) content, volatile solids (VS) content, chemical oxygen demand (COD) and pH. The TS and VS analyses were performed according to the American Public Health Association Standard Methods . The VFAs content and COD were determined using Nanocolor® kits (Macherey–Nagel, Germany). Methane content was analyzed by GC–MS gas chromatography (Agilent, USA).
Library preparation and sequencing
Metagenomic DNA isolated from environmental and laboratory communities (combined triplicates) was used for library preparation with an Illumina TruSeq DNA Sample Preparation Kit according to the manufacturer’s protocol. Purifications of DNA fragments were performed with Agencourt AMPure XP beads (Beckman Coulter). The libraries were analyzed by electrophoresis on 2% agarose gels (1× TAE buffer) with GelGreen staining, 2100 Bioanalyzer (Agilent) High-Sensitivity DNA Assay and KAPA Library Quantification Kit for Illumina. The libraries obtained were sequenced on the Illumina HiSeq 1500 platform (HiSeq Reagent Kit v2, 300 cycles) in a pair-end mode with a read length of 150 bp.
Sequenced data analysis
Metagenomic raw sequences were uploaded to Metagenomic Rapid Annotations using Subsystems Technology (MG-RAST) server . The metagenomes used in this work are available under the project accession mgp16315. Taxonomic profiles of consortia were created against the RefSeq database and the functional profiles were generated using the matches to the SEED Subsystems database with default parameters.
Phylogenetic classification of functional annotations relevant to the methanogenesis process was performed using MetAnnotate  based on HMM search and phylogenetic placement and best-hit approach with default parameters. The available HMM profiles of PFAM  and TIGRFAM  protein families were included in the analysis for the following enzymes: formylmethanofuran dehydrogenase (PF02663, TIGR03121, TIGR03122); formylmethanofuran-tetrahydromethanopterin formyltransferase (PF01913); methenyltetrahydromethanopterin cyclohydrolase (TIGR03120); H2-forming N5,N10-methylene-tetrahydromethanopterin dehydrogenase (PF03201); methylene-5,6,7,8-tetrahydromethanopterin dehydrogenase (PF01993); 5,10-methylenetetrahydromethanopterin reductase (TIGR03555); acetate kinase (TIGR00016); phosphate acetyltransferase (TIGR00651); acetyl-CoA synthetase (PF16177); CO dehydrogenase/acetyl-CoA synthase (PF03598, PF03599, TIGR00314, TIGR00315, TIGR00381); tetrahydromethanopterin S-methyltransferase (PF04208, PF05440, PF04211, PF04207, PF04206, PF09472, PF04210, PF02007, TIGR01111, TIGR04166, TIGR01148, TIGR01112, TIGR01113, TIGR02507, TIGR01114, TIGR00314, TIGR00315, TIGR00381), methanol methyltransferase (PF12176); methylamine methyltransferase (PF05369, PF06253, PF09505, TIGR02368, TIGR02369); methyl-CoM reductase (PF02249, PF02745, PF02241, PF02783, PF04609, PF02505, PF02240, TIGR03256, TIGR03257, TIGR03259, TIGR03260, TIGR03264); CoB-CoM heterodisulfide reductase (TIGR03288, TIGR03290).
Phylogenetic assignments of methanogenesis sequences are listed in Tax-Fun MetAnnotate, Additional file 2. For better readability, key methanogenesis annotations were combined (see Tax-Fun MetAnnotate Combined, Additional file 2) and presented as heatmaps created with the Statistical Analysis of Metabolic Profile (STAMP) software .
Taxonomic diversity (Shannon–Wiener index, H) was calculated using diversity function from vegan package in R . Eveness (J, Pielou index) was calculated using formula J = H/log(S), where H is Shannon–Wiener index and S is the species richness for given sample. Bray–Curtis distances were calculated using vegdist function from vegan package and processed using metaMDS function to produce multidimensional scaling (MDS) plots.
AP was involved in planning of the metagenomics approach and was responsible for isolation of metagenomic DNA, deep sequencing and bioinformatic analyses; MC was involved in the laboratory cultivation experiments and physico-chemical measurements; PSK participated in bioinformatic and statistical analyses; LD was involved in planning and supervision of the laboratory cultivation experiments. The manuscript was written by AP, consulted and corrected by LL, AS, LD, PSK and AD; LL and AS conceived and directed the study. Funding for this work was provided by LL and LD. All authors read and approved the final manuscript.
We acknowledge Joanna Karczewska-Golec for critical reading of this manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files. The metagenome sequences are available at MG-RAST server under the project accession mgp16315. Data obtained from MG-RAST and MetAnnotate analysis were collected in Additional file 2.
Consent for publication
Ethics approval and consent to participate
This study was funded by the Applied Research Program of the National Centre for Research and Development (Grant Number PBS1/A8/3/2012) and by the European Union’s European Regional Development Fund through the Innovative Economy Operational Program, 2007–2013 (Project Support Agreement POIG.01.01.02-14-054/09-00). Experiments were carried out with the use of CePT infrastructure financed by the European Union’s European Regional Development Fund through the Innovative Economy Operational Program, 2007–2013 (Project Support Agreement POIG.02.02.00-14-024/08-00).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Chaban B, Ng SY, Jarrell KF. Archaeal habitats-from the extreme to the ordinary. Can J Microbiol. 2006;52(2):73–116.PubMedGoogle Scholar
- Hedderich R, Whitman WB. Physiology and biochemistry of the methane-producing Archaea. In: Dworkin M, Falkow S, Rosenberg E, et al., editors. The prokaryotes, vol. 3. New York: Springer; 2006. p. 1050–79.Google Scholar
- Borrel G, O’Toole PW, Harris HM, Peyret P, Brugere JF, Gribaldo S. Phylogenomic data support a seventh order of methylotrophic methanogens and provide insights into the evolution of methanogenesis. Genome Biol Evol. 2013;5:1769–80.PubMedPubMed CentralGoogle Scholar
- Liu Y, Whitman WB. Metabolic, phylogenetic, and ecological diversity of the methanogenic archaea. Ann N Y Acad Sci. 2008;1125:171–89.PubMedGoogle Scholar
- Thauer RK, Kaster AK, Seedorf H, Buckel W, Hedderich R. Methanogenic archaea: ecologically relevant differences in energy conservation. Nat Rev Microbiol. 2008;6(8):579–91.PubMedGoogle Scholar
- Jablonski S, Rodowicz P, Lukaszewicz M. Methanogenic archaea database containing physiological and biochemical characteristics. Int J Syst Evol Microbiol. 2015;65(4):1360–8.PubMedGoogle Scholar
- Lang K, Schuldes J, Klingl A, Poehlein A, Daniel R, Brunea A. New mode of energy metabolism in the seventh order of methanogens as revealed by comparative genome analysis of “Candidatus Methanoplasma termitum”. Appl Environ Microbiol. 2015;81(4):1338–52.PubMedPubMed CentralGoogle Scholar
- Yang Y, Yu K, Xia Y, Lau FT, Tang DT, Fung WC, Fang HH, Zhang T. Metagenomic analysis of sludge from full-scale anaerobic digesters operated in municipal wastewater treatment plants. Appl Microbiol Biotechnol. 2014;98(12):5709–18.PubMedGoogle Scholar
- Stolze Y, Zakrzewski M, Maus I, Eikmeyer F, Jaenicke S, Rottmann N, Siebner C, Pühler A, Schlüter A. Comparative metagenomics of biogas-producing microbial communities from production-scale biogas plants operating under wet or dry fermentation conditions. Biotechnol Biofuels. 2015;8:14.PubMedPubMed CentralGoogle Scholar
- Luo G, Fotidis IA, Angelidaki I. Comparative analysis of taxonomic, functional, and metabolic patterns of microbiomes from 14 full-scale biogas reactors by metagenomic sequencing and radioisotopic analysis. Biotechnol Biofuels. 2016;9:51.PubMedPubMed CentralGoogle Scholar
- Rotaru AE, Shrestha PM, Liu F, Shrestha M, Shrestha D, Embree M, Zengler K, Wardman C, Nevin P, Lovley DR. A new model for electron flow during anaerobic digestion: direct interspecies electron transfer to Methanosaeta for the reduction of carbon dioxide to methane. Energy Environ Sci. 2014;7(1):408–15.Google Scholar
- Hess M, Sczyrba A, Egan R, Kim TW, Chokhawala H, Schroth G, Luo S, Clark DS, Chen F, Zhang T, Mackie RI, Pennacchio LA, Tringe SG, Visel A, Woyke T, Wang Z, Rubin EM. Metagenomic discovery of biomass-degrading genes and genomes from cow rumen. Science. 2011;331(6016):463–7.PubMedGoogle Scholar
- Borrel G, Parisot N, Harris HM, Peyretaillade E, Gaci N, Tottey W, Bardot O, Raymann K, Gribaldo S, Peyret P, O’Toole PW, Brugère JF. Comparative genomics highlights the unique biology of Methanomassiliicoccales, a thermoplasmatales-related seventh order of methanogenic archaea that encodes pyrrolysine. BMC Genomics. 2014;15:679.PubMedPubMed CentralGoogle Scholar
- Chaudhary PP, Gaci N, Borrel G, O’Toole PW, Brugère JF. Molecular methods for studying methanogens of the human gastrointestinal tract: current status and future directions. Appl Microbiol Biotechnol. 2015;99(14):5801–15.PubMedGoogle Scholar
- Wallace RJ, Rooke JA, McKain N, Duthie CA, Hyslop JJ, Ross DW, Waterhouse A, Watson M, Roehe R. The rumen microbial metagenome associated with high methane production in cattle. BMC Genomics. 2015;16:839.PubMedPubMed CentralGoogle Scholar
- Evans PN, Parks DH, Chadwick GL, Robbins SJ, Orphan VJ, Golding SD, Tyson GW. Methane metabolism in the archaeal phylum Bathyarchaeota revealed by genome-centric metagenomics. Science. 2015;350(6259):434–8.PubMedGoogle Scholar
- Bizic-Ionescu M, Klintzch T, Ionescu D, Hindiyeh MY, Guenthel M, Pastor AMM, Keppler F, Grossart HP. Widespread formation of methane by Cyanobacteria in aquatic and terrestrial environments. bioRxiv. 2018. https://doi.org/10.1101/398958 Google Scholar
- Dziewit L, Pyzik A, Romaniuk K, Sobczak A, Szczesny P, Lipinski L, Bartosik D, Drewniak L. Novel molecular markers for the detection of methanogens and phylogenetic analyses of methanogenic communities. Front Microbiol. 2015;6:694.PubMedPubMed CentralGoogle Scholar
- Weiland P. Biogas production: current state and perspectives. Appl Microbiol Biotechnol. 2010;85(4):849–60.PubMedGoogle Scholar
- Daelman MR, van Voorthuizen EM, van Dongen UG, Volcke EI, van Loosdrecht MC. Methane emission during municipal wastewater treatment. Water Res. 2012;46(11):3657–70.PubMedGoogle Scholar
- Takai K, Moser DP, DeFlaun M, Onstott TC, Fredrickson JK. Archaeal diversity in waters from deep South African gold mines. Appl Environ Microbiol. 2001;67(12):5750–60.PubMedPubMed CentralGoogle Scholar
- Tveit AT, Urich T, Svenning MM. Metatranscriptomic analysis of Arctic peat soil microbiota. Appl Environ Microbiol. 2014;80(18):5761–72.PubMedPubMed CentralGoogle Scholar
- Schmidt O, Hink L, Horn MA, Drake HL. Peat: home to novel syntrophic species that feed acetate- and hydrogen-scavenging methanogens. ISME J. 2016;10(8):1954–66.PubMedPubMed CentralGoogle Scholar
- Meyer F, Paarmann D, D’Souza M, Olson R, Glass EM, Kubal M, Paczian T, Rodriguez A, Stevens R, Wilke A, Wilkening J, Edwards RA. The metagenomics RAST server—a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinform. 2008;9:386.Google Scholar
- Petrenko P, Lobb B, Kurtz DA, Neufeld JD, Doxey AC. MetAnnotate: function-specific taxonomic profiling and comparison of metagenomes. BMC Biol. 2015;13:92.PubMedPubMed CentralGoogle Scholar
- Ward NL, Challacombe JF, Janssen PH, Henrissat B, Coutinho PM, Wu M, Xie G, Haft DH, Sait M, Badger J, Barabote RD, Bradley B, Brettin TS, Brinkac LM, Bruce D, Creasy T, Daugherty SC, Davidsen TM, DeBoy RT, Detter JC, Dodson RJ, Durkin AS, Ganapathy A, Gwinn-Giglio M, Han CS, Khouri H, Kiss H, Kothari SP, Madupu R, Nelson KE, Nelson WC, Paulsen I, Penn K, Ren Q, Rosovitz MJ, Selengut JD, Shrivastava S, Sullivan SA, Tapia R, Thompson LS, Watkins KL, Yang Q, Yu C, Zafar N, Zhou L, Kuske CR. Three genomes from the phylum Acidobacteria provide insight into the lifestyles of these microorganisms in soils. Appl Environ Microbiol. 2009;75(7):2046–56.PubMedPubMed CentralGoogle Scholar
- Purushe J, Fouts DE, Morrison M, White BA, Mackie RI, North American Consortium for Rumen Bacteria, Coutinho PM, Henrissat B, Nelson KE. Comparative genome analysis of Prevotella ruminicola and Prevotella bryantii: insights into their environmental niche. Microb Ecol. 2010;60(4):721–9.PubMedGoogle Scholar
- Morita M, Malvankar NS, Franks AE, Summers ZM, Giloteaux L, Rotaru AE, Rotaru C, Lovley DR. Potential for direct interspecies electron transfer in methanogenic wastewater digester aggregates. MBio. 2011;2(4):e00159-11.PubMedPubMed CentralGoogle Scholar
- Wirth R, Kovács 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.PubMedPubMed CentralGoogle Scholar
- Sundberg C, Al-Soud WA, Larsson M, Alm E, Yekta SS, Svensson BH, Sørensen SJ, Karlsson A. 454 pyrosequencing analyses of bacterial and archaeal richness in 21 full-scale biogas digesters. FEMS Microbiol Ecol. 2013;85(3):612–26.PubMedGoogle Scholar
- Solli L, Håvelsrud OE, Horn SJ, Rike AG. A metagenomic study of the microbial communities in four parallel biogas reactors. Biotechnol Biofuels. 2014;7(1):146.PubMedPubMed CentralGoogle Scholar
- Dai X, Tian Y, Li J, Luo Y, Liu D, Zheng H, Wang J, Dong Z, Hu S, Huang L. Metatranscriptomic analyses of plant cell wall polysaccharide degradation by microorganisms in the cow rumen. Appl Environ Microbiol. 2015;81(4):1375–86.PubMedPubMed CentralGoogle Scholar
- De Francisci D, Kougias PG, Treu L, Campanaro S, Angelidaki I. Microbial diversity and dynamicity of biogas reactors due to radical changes of feedstock composition. Bioresour Technol. 2015;176:56–64.PubMedGoogle Scholar
- Nathani NM, Kothari RK, Patel AK, Joshi CG. Functional characterization reveals novel putative coding sequences in Prevotella ruminicola genome extracted from rumen metagenomic studies. J Mol Microbiol Biotechnol. 2015;25(4):292–9.PubMedGoogle Scholar
- Poszytek K, Pyzik A, Sobczak A, Lipinski L, Sklodowska A, Drewniak L. The effect of the source of microorganisms on adaptation of hydrolytic consortia dedicated to anaerobic digestion of maize silage. Anaerobe. 2017;46:46–55.PubMedGoogle Scholar
- Chen S, Cheng H, Wyckoff KN, He Q. Linkages of Firmicutes and Bacteroidetes populations to methanogenic process performance. J Ind Microbiol Biotechnol. 2016;43(6):771–81.PubMedGoogle Scholar
- Schwalm ND 3rd, Groisman EA. Navigating the gut buffet: control of polysaccharide utilization in Bacteroides spp. Trends Microbiol. 2017;25(12):1005–15.PubMedGoogle Scholar
- Güllert S, Fischer MA, Turaev D, Noebauer B, Ilmberger N, Wemheuer B, Alawi M, Rattei T, Daniel R, Schmitz RA, Streit WR, Grundhoff A. Deep metagenome and metatranscriptome analyses of microbial communities affiliated with an industrial biogas fermenter, a cow rumen, and elephant feces reveal major differences in carbohydrate hydrolysis strategies. Biotechnol Biofuels. 2016;9(1):121.PubMedPubMed CentralGoogle Scholar
- Campanaro S, Treu L, Kougias PG, Luo G, Angelidaki I. Metagenomic binning reveals the functional roles of core abundant microorganisms in twelve full-scale biogas plants. Water Res. 2018;140:123–34.PubMedGoogle Scholar
- Nettmann E, Bergmann I, Mundt K, Linke B, Klocke M. Archaea diversity within a commercial biogas plant utilizing herbal biomass determined by 16S rDNA and mcrA analysis. J Appl Microbiol. 2008;105(6):1835–50.PubMedGoogle Scholar
- De Vrieze J, Hennebel T, Boon N, Verstraete W. Methanosarcina: the rediscovered methanogen for heavy duty biomethanation. Bioresour Technol. 2012;112:1–9.PubMedGoogle Scholar
- Goux X, Calusinska M, Fossépré M, Benizri E, Delfosse P. Start-up phase of an anaerobic full-scale farm reactor—appearance of mesophilic anaerobic conditions and establishment of the methanogenic microbial community. Bioresour Technol. 2016;212:217–26.PubMedGoogle Scholar
- Singh KM, Patel AK, Shah RK, Reddy B, Joshi CG. Potential functional gene diversity involved in methanogenesis and methanogenic community structure in Indian buffalo (Bubalus bubalis) rumen. J Appl Genet. 2015;56(3):411–26.PubMedGoogle Scholar
- St-Pierre B, Wright AD. Diversity of gut methanogens in herbivorous animals. Animal. 2013;7(Suppl 1):49–56PubMedGoogle Scholar
- Fernandes KA, Kittelmann S, Rogers CW, Gee EK, Bolwell CF, Bermingham EN, Thomas DG. Faecal microbiota of forage-fed horses in New Zealand and the population dynamics of microbial communities following dietary change. PLoS ONE. 2014;9(11):e112846.PubMedPubMed CentralGoogle Scholar
- Söllinger A, Schwab C, Weinmaier T, Loy A, Tveit AT, Schleper C, Urich T. Phylogenetic and genomic analysis of Methanomassiliicoccales in wetlands and animal intestinal tracts reveals clade-specific habitat preferences. FEMS Microbiol Ecol. 2016;92(1):fiv149.PubMedGoogle Scholar
- Timmers PH, Welte CU, Koehorst JJ, Plugge CM, Jetten MS, Stams AJ. Reverse methanogenesis and respiration in methanotrophic archaea. Archaea. 2017;2017:1654237.PubMedPubMed CentralGoogle Scholar
- Iguchi H, Yurimoto H, Sakai Y. Interactions of methylotrophs with plants and other heterotrophic bacteria. Microorganisms. 2015;3(2):137–51.PubMedPubMed CentralGoogle Scholar
- Drewniak L, Krawczyk PS, Mielnicki S, Adamska D, Sobczak A, Lipinski L, Burec-Drewniak W, Sklodowska A. Physiological and metagenomic analyses of microbial mats involved in self-purification of mine waters contaminated with heavy metals. Front Microbiol. 2016;7:1252.PubMedPubMed CentralGoogle Scholar
- Calderón K, González-Martínez A, Gómez-Silván C, Osorio F, Rodelas B, González-López J. Archaeal diversity in biofilm technologies applied to treat urban and industrial wastewater: recent advances and future prospects. Int J Mol Sci. 2013;14(9):18572–98.PubMedPubMed CentralGoogle Scholar
- Wu JH, Chuang HP, Hsu MH, Chen WY. Multiplexed relative abundance analysis of methanogens in anaerobic digestion systems using a hierarchical oligonucleotide primer extension approach. Appl Environ Microbiol. 2013;79(24):7598–609.PubMedPubMed CentralGoogle Scholar
- Wojcieszak M, Pyzik A, Poszytek K, Krawczyk PS, Sobczak A, Lipinski L, Roubinek O, Palige J, Sklodowska A, Drewniak L. Adaptation of methanogenic inocula to anaerobic digestion of maize silage. Front Microbiol. 2017;8:1881.PubMedPubMed CentralGoogle Scholar
- Bremges A, Maus I, Belmann P, Eikmeyer F, Winkler A, Albersmeier A, Pühler A, Schlüter A, Sczyrba A. Deeply sequenced metagenome and metatranscriptome of a biogas-producing microbial community from an agricultural production-scale biogas plant. Gigascience. 2015;4(1):33.PubMedPubMed CentralGoogle Scholar
- Maus I, Koeck DE, Cibis KG, Hahnke S, Kim YS, Langer T, Stolze Y, et al. Unraveling the microbiome of a thermophilic biogas plant by metagenome and metatranscriptome analysis complemented by characterization of bacterial and archaeal isolates. Biotechnol Biofuels. 2016;9(1):171.PubMedPubMed CentralGoogle Scholar
- Hassa J, Maus I, Off S, Pühler A, Scherer P, Klocke M, Schlüter A. Metagenome, metatranscriptome, and metaproteome approaches unraveled compositions and functional relationships of microbial communities residing in biogas plants. Appl Microbiol Biotechnol. 2018;102(12):5045–63.PubMedPubMed CentralGoogle Scholar
- Speda J, Jonsson BH, Carlsson U, Karlsson M. Metaproteomics-guided selection of targeted enzymes for bioprospecting of mixed microbial communities. Biotechnol Biofuels. 2017;10(1):128.PubMedPubMed CentralGoogle Scholar
- Luton PE, Wayne JM, Sharp RJ, Riley PW. The mcrA gene as an alternative to 16S rRNA in the phylogenetic analysis of methanogen populations in landfill. Microbiology. 2002;148(Pt 11):3521–30.PubMedGoogle Scholar
- Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, Glöckner FO. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013;41(1):e1.PubMedGoogle Scholar
- Krylowicz A, Chrzanowski K, Usidus J. Sposób i układ wytwarzania metanu i energii elektrycznej i cieplnej. Polish Patent No 19759. Warsaw: PL Patent and Trademark Office; 2001.Google Scholar
- American Public Health Association [APHA]. Standard methods for the examination of water and wastewater. 18th ed. Washington, DC: APHA; 1998.Google Scholar
- Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, Sonnhammer EL, Tate J, Punta M. Pfam: the protein families database. Nucleic Acids Res. 2014;42(Database issue):D222–30.PubMedGoogle Scholar
- Haft DH, Selengut JD, Richter RA, Harkins D, Basu MK, Beck E. TIGRFAMs and genome properties in 2013. Nucleic Acids Res. 2013;41(Database issue):D387–95.PubMedGoogle Scholar
- Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30(21):3123–4.PubMedPubMed CentralGoogle Scholar
- Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Szoecs E, Wagner H. 2018. vegan: community ecology package. R package version 2.5-3. https://CRAN.R-project.org/package=vegan.