Comparative analysis of taxonomic, functional, and metabolic patterns of microbiomes from 14 full-scale biogas reactors by metagenomic sequencing and radioisotopic analysis
Biotechnology for Biofuels volume 9, Article number: 51 (2016)
Biogas production is a very complex process due to the high complexity in diversity and interactions of the microorganisms mediating it, and only limited and diffuse knowledge exists about the variation of taxonomic and functional patterns of microbiomes across different biogas reactors, and their relationships with the metabolic patterns. The present study used metagenomic sequencing and radioisotopic analysis to assess the taxonomic, functional, and metabolic patterns of microbiomes from 14 full-scale biogas reactors operated under various conditions treating either sludge or manure.
The results from metagenomic analysis showed that the dominant methanogenic pathway revealed by radioisotopic analysis was not always correlated with the taxonomic and functional compositions. It was found by radioisotopic experiments that the aceticlastic methanogenic pathway was dominant, while metagenomics analysis showed higher relative abundance of hydrogenotrophic methanogens. Principal coordinates analysis showed the sludge-based samples were clearly distinct from the manure-based samples for both taxonomic and functional patterns, and canonical correspondence analysis showed that the both temperature and free ammonia were crucial environmental variables shaping the taxonomic and functional patterns. The study further the overall patterns of functional genes were strongly correlated with overall patterns of taxonomic composition across different biogas reactors.
The discrepancy between the metabolic patterns determined by metagenomic analysis and metabolic pathways determined by radioisotopic analysis was found. Besides, a clear correlation between taxonomic and functional patterns was demonstrated for biogas reactors, and also the environmental factors that shaping both taxonomic and functional genes patterns were identified.
Anaerobic digestion has been widely used for the treatment of organic wastes with simultaneous production of biogas. Biogas production from organic wastes includes four sequential metabolic steps (hydrolysis, fermentation, acetogenesis, and methanogenesis) and involves a number of different microorganisms (i.e., bacteria, archaea, fungi, and protozoa) . The role and interaction of these microorganisms is very complex [2, 3], and it is vital to understand the taxonomic and functional dynamics of anaerobic digestion microbiomes at different settings of the biomethanation process, and their correlation with the metabolic pathway, in order to improve the process performance.
Several culture-independent molecular methods based on 16S rRNA genes have been developed to investigate and characterize the microbiomes in biogas reactors [4–9]. The culture-independent methods include polymerase chain reaction (PCR)-denaturing gradient gel electrophoresis, PCR-terminal restriction fragment length polymorphism, PCR-cloning, and the recently developed PCR-high-throughput sequencing . Numerous studies on microbial composition with respect to physical, chemical, and biological characteristics of biogas reactors have been published [1, 3, 11–14]. It is now known that the microbial composition is influenced by environmental variables such as temperature, feedstock, biogas reactor configurations, et al. [11, 12, 15–18]. Moreover, it is known that not only aceticlastic methanogens but also, depending on the operational conditions (e.g., ammonia concentration, et al.), hydrogenotrophic methanogens play a significant role in methanogenesis [19–21]. In addition, biogas reactors, operating at constant conditions (feedstock, temperature, et al.), have demonstrated an unprecedented level of stability with a unique community structure .
Our understanding of the microbial composition in biogas reactors has been increased greatly with the establishment of culture-independent molecular methods . However, these molecular methods have several limitations, such as PCR bias , and lack of information about the functional genes of the microbiomes . The ongoing development of high-throughput molecular tools and bioinformatics allows sequencing of the bulk DNA instead of only 16S rRNA genes and thereby provides both taxonomic and functional information of microbiomes to an extent that was unimaginable even a few years ago . It should be noted that traditional microbiological methodologies (e.g., isolation and cultivation of pure strains) have to be employed in order to study the physiology, metabolism, et al. for new isolates derived from biogas reactors, which could not be accomplished by metagenomic sequencing. Therefore, the combination of the new molecular technologies with traditional microbiological methodologies is necessary for future studies .
Metagenomic sequencing has been performed on different environments (agricultural soil, acid mine biofilm, sea, et al. ), and the first metagenomic analysis of biogas reactors was reported in 2008 . Metagenomic studies on biogas reactors lacked an understanding of how functional genes encoded in their collective genomes, act across different biogas reactors, especially in correlation with different and/or changing environmental parameters (e.g., feedstock, temperature, process by-products such as ammonium, free ammonia nitrogen, or acids) [23, 26, 27]. In addition, previous studies estimated metabolic pathways (especially for methanogenesis) based on the corresponding functional genes by metagenomic analysis [26, 27]. Nevertheless, functional genes from metagenomic analysis only reflect the potential enzymes that could be synthesized by the microbes, and it is still not clear whether there is a direct correlation between metagenomic results and actual metabolic pathways taking place in the biogas reactors. Radioisotopic analysis and proteomic analysis have been widely used in the identification of the actual dominant methanogenic pathway in biogas reactors [9, 28].
Therefore, based on the above considerations, the aim of the present study was to conduct a detailed comparative analysis of both taxonomic and functional patterns of microbiomes from 14 different full-scale biogas reactors with a broad range of operational conditions by metagenomic sequencing and also to determine the predictability of actual metabolism by taxonomic and functional information. Specifically, the current study was designed to address the following questions: How is the relationship between the taxonomic and functional patterns of the microbiomes from biogas reactors and their correlation with the environmental variables? Is it possible to use the information on the taxonomic and functional compositions to predict the metabolic pathways (e.g., methanogenesis)?
Results and discussion
Overview of the metagenomics data
In total, around 400 million of 100 bp metagenomic reads for the 14 samples were obtained by paired-end sequencing (Additional file 1: Table S1), and the paired-end sequences were then joined to be joined reads with length around 170 bp (8,021,985–26,750,735 per sample). All the sequences were sub-sampled to 8,021,985 joined reads and then submitted to MG-RAST for further analysis. The percentages of identified 16S rRNA genes, used for the taxonomic assignments , were between 0.042 and 0.132 % for all the samples, which were consistent with previous studies analyzing microbes in biogas reactors or in wastewater treatment bioreactors [5, 29]. The ratios of sequences annotated based on SEED subsystems by MG-RAST were between 19 and 35 % of the total sequences for all the samples (Additional file 1: Table S1). The annotated sequences levels were similar to those reported in previous studies, where metagenomic sequencing was used to characterize microbiomes in soil [30, 31] and in other highly diverse microbial habitats .
In detail, six and eleven phyla (with relative abundance higher than 1 %) were identified in manure-based and sludge-based samples, respectively, demonstrating more diverse microbial communities for the sludge-based samples (Fig. 1). Firmicutes (54.8–75.8 %) and Bacteroidetes (3.5–20.2 %) were dominant in all the manure-based samples, while Proteobacteria (26.4–34.5 %) followed by Firmicutes (9.1–15.2 %) and Bacteroidetes (9.1–21.9 %) were dominant in all the sludge-based samples. The relative abundance of Firmicutes in manure-based samples (>50 %) was much higher than the sludge-based samples (<15 %). The dominance of Firmicutes in manure-based samples has also been reported in other studies [1, 10]. Clostridia, belonging to Firmicutes, were the most abundant class (>40 % of all the bacteria sequences) (Additional file 1: Table S2). For the manure-based samples, the relative abundances of Firmicutes and Bacteroidetes were more diverse for mesophilic samples compared with the thermophilic ones. The genus Dechloromonas , Syntrophorhabdus , and Syntrophus , which were chlorate-reducing, aromatic-degrading, and benzoate-degrading bacteria, were mainly found in sewage sludge-based samples, and it could be related to the corresponding pollutants contained in sewage sludge . The cellulolytic bacteria Halocella  was mainly found in manure-based samples, and it could be due to the high lignocellulose contains in manure . The presence of Halocella in biogas reactors, especially thermophilic reactors, was also reported previously .The observation of Syntrophaceticus mainly in samples MT2a, MT2b, MT3a, MT3b, and MT4 suggested syntrophic acetate-oxidizing occurred in the corresponding biogas reactors , which will be discussed later with the radioisotopic analysis data. In addition, it was found as expected that serial operated biogas reactors (MT2a and MT2b, MT3a, and MT3b) had similar microbial communities in both process steps (Additional file 1: Table S2). This is in accordance with our previous study where the microbial communities in the second reactor were found to be strongly correlated with the first reactor .
The relative abundance of sequences assigned to Archaea (i.e., methanogens) was low (from 1.2 to 11.4 %) for all samples (Additional file 1: Table S1), which was consistent with other studies . Nevertheless, it does not necessarily mean that archaea had lower metabolic activity, since they might have higher transcriptional activity compared to bacteria, as has been recently demonstrated by comparing taxonomic results from 16S rRNA genes, metagenomic, and metatranscriptomic analyses . Methanomicrobiales, Methanobacteriales, and Methanosarcinales were the three dominant orders among all samples tested (Fig. 1b). Moreover, it should be noted that the dominant genus of the order Methanosarcinales was different for manure-based (Methanosarcina) and sludge-based samples (Methanosaeta) (Table 2), which could be related to the low VFA concentration in sludge-based samples compared to manure-based samples. Methanosaeta has higher affinity for acetate and grow better with lower acetate concentrations compared to Methanosarcina [44, 45].
For manure-based samples from thermophilic biogas reactors, MT1 sample was dominated by Methanosarcinales (Mainly Methanosarcina (around 90 %) as seen in Table 2, mediating aceticlastic methanogenesis and in some cases also hydrogenotrophic methanogenesis. ), while MT2a, MT2b, MT3a, and MT3b samples were dominated by hydrogenotrophic methanogens Methanomicrobiales. Differences between the relative abundances of the archaea orders among sludge-based samples were also observed. Sample SM1 had higher percentage of hydrogenotrophic methanogens Methanomicrobiales, while all the other samples (SM2–SM5) had higher percentages of Methanosarcinales (Mainly Methanosaeta (40–60 %) as seen in Table 2, mediating both aceticlastic methanogenesis). Together, these observations suggest that there might be different dominant microorganisms mediating either hydrogenotrophic or aceticlastic methanogenesis for the samples from similar environments (e.g., manure-based samples from thermophilic biogas reactors and sludge-based samples from mesophilic biogas reactors).
The major functional categories for all the samples were those involved in metabolism of carbohydrates, clustering-based subsystems (containing such functions as proteosomes, ribosomes, and recombination-related clusters ), protein metabolism, amino acids, and derivatives (Additional file 1: Fig S1). Similar major functional genes were also found in the microbiomes of other ecosystems such as wastewater treatment system (activated sludge) , desert soil , and freshwater . The comparison of the functional genes of sludge-based and manure-based samples (Fig. 2) demonstrated that 24 out of the total 25 major functional categories were found to be significantly different (p < 0.05 based on ANOVA analysis). There were no significant differences between manure-based and sludge-based samples for the functional category “Regulation and cell signaling.” There were several difference (e.g., nitrite, poly-3-hydroxybutyrate/poly-3-hydroxyvalerate, aromatic compounds, et al.) in the composition of manure and sludge which promoted the functionality variances observed. For example, the genes assigned into nitrogen metabolism and phosphorus metabolism had significantly higher abundances in sludge-based samples than manure-based samples, which was due to these specific processes taking place in wastewater treatment plants. In detail, the sludge-based samples were derived from wastewater treatment process and inevitably nitrate (from the aeration tank)  and poly-3-hydroxybutyrate/poly-3-hydroxyvalerate (accumulated in the waste-activated sludge for phosphorus removal)  were entering the biogas reactors. Furthermore, the sludge-based samples had more genes related to the metabolism of aromatic compounds. This can be explained by the higher abundance of aromatic compounds and other organic contaminants in the sludge originating from the wastewater, compared to manure-based digesters, which would promote establishment of microbes involved with metabolism of these compounds , and it was also consistent with the observation of aromatic–degrading Syntrophorhabdus mainly in sewage sludge samples (Additional file 1: Table S2) .
In order to make a more detailed analysis on the functional genes, the joined reads were annotated to metabolic pathways based on KEGG database. Protein and carbohydrate were the dominant compounds in manure and sewage sludge [41, 52], and therefore, the key genes relating to protein and carbohydrate anaerobic degradation were analyzed as shown in Additional file 1: Fig S2. The relative abundance of genes relating to cellulose and hemicellulose degradation were all higher in manure-based samples compared to sewage sludge-based samples, which was consistent with the high fiber content (lignocellulose materials) in manure . For the genes relating to protein degradation, the relative abundances were similar for both sewage and manure-based samples. The genes relating to the metabolic pathways of methanogenesis [hydrogenotrophic (H2/CO2), acetoclastic (acetate), and methylotrophic (methanol)] were also analyzed, and the results are shown in Fig. 3 and Additional file 1: Table S3 [26, 53]. The dominant genes in hydrogenotrophic methanogenesis were formate dehydrogenase (EC: 126.96.36.199) and formylmethanofuran dehydrogenase (EC: 188.8.131.52), which were involved in the initial step of hydrogenotrophic pathway. The dominant genes in aceticlastic methanogenesis were acetyl-CoA synthetase (EC: 184.108.40.206) and acetyl-CoA decarbonylase/synthase complex (ACDS), and they were essential in the synthesis of acetyl-CoA from acetate. The most abundant genes were found to be related to hydrogenotrophic and aceticlastic methanogenic pathways in this study, and an interesting phenomenon was that hydrogenotrophic methanogens were absolutely dominant (the relative abundance higher than 95 % of archaea) in samples MM2, MT2a, MT2b, MT3a, and MT3b (Additional file 1: Fig S1). Additionally, genes encoding for aceticlastic methanogenesis is unique for this pathway and can therefore not be detected in hydrogenotrophic methanogens . In fact, acetyl-CoA synthetase and acetyl-CoA decarbonylase/synthase complex also exist in bacteria and are involved in other metabolic pathways (acetate oxidation to produce H2 and CO2, homoacetogenesis to produce acetate, et al.) . Thus, the claim that the abundance of different genes can be used to assess the methanogenic pathways in the biogas reactors [26, 27] seems not hold merit.
Methanogenic pathway determined by radioisotopic analysis
It is generally assumed that aceticlastic methanogenesis is the dominant pathway when 14CO2/14CH4<1, while the hydrogenotrophic methanogenesis is the main pathway when 14CO2/14CH4>1 . The results from radioisotopic analysis showed (Additional file 1: Fig S3) that all the sludge-based samples had 14CO2/14CH4<1, indicating the dominance of aceticlastic methanogenesis pathway. However, the dominant methanogenic pathway for all the manure-based samples except MT1 was hydrogenotrophic (14CO2/14CH4>1). The hydrogenotrophic pathway is coupled with syntrophic acetate oxidation that first converts acetate to H2/CO2 and then H2/CO2 is converted to methane by hydrogenotrophic methanogens. Syntrophic acetate oxidation is generally favored in the presence of inhibitors, particularly ammonium and volatile fatty acids , which may inhibit the aceticlastic methanogens. As shown in Table 1, all the manure-based samples had relatively higher ammonia and volatile fatty acids compared with sludge-based samples, which might be the reason for the dominance of hydrogenotrophic methanogenesis pathway in the manure-based samples. The dominance of syntrophic acetate-oxidizing Syntrophaceticus in MT2a, MT2b, MT3a, MT3b, and MT4 agreed well with the above results (Additional file 1: Table S2) . VFA and ammonia had synergetic effect on the dominant methanogenic pathways as revealed by Lu et al. . For MT1, the lower VFA and ammonia concentrations compared with the other manure-based samples might led to the dominance of aceticlastic methanogenesis pathway. It should be noted that Methanosarcina was dominant (90 %, Table 2) in MT1, while Methanoculleus was dominant in MT2a, MT2b, MT3a, MT3b, and MT4 (50–82 %, Table 2).
Methanosarcinales is the only order that could mediate methane production from acetate . Methanosaeta (strict acetoclastic methanogens) and Methanosarcina (aceticlastic or hydrogenotrophic methanogens) were found to be the two main genus found in the order Methanosarcinales (Table 2). The correlation between the percentage of Methanosarcinales in total archaea sequences and the values of 14CO2/14CH4 from radioisotopic analysis is shown in Additional file 1: Fig S4. It is obvious that when the ratio of 14CO2/14CH4 was higher than two, Methanosarcinales were absent indicating an absolute dominance of hydrogenotrophic methanogens. Therefore, a 14CO2/14CH4 ratio threshold for absolute dominance of hydrogenotrophic methanogens in full-scale biogas reactors can be proposed. When the 14CO2/14CH4 ratio was between one and two, the percentage of Methanosarcinales was between 0–43 % which still showed the dominance of hydrogenotrophic methanogens. However, when the values of 14CO2/14CH4 were lower than one, the percentage of Methanosarcinales was above 43 %, with the exception of sample SM1, for which it was 23 %. It was further found that strict aceticlastic methanogens Methanosaeta was the predominant genus of Methanosarcinales for SM1 (22.6 % as shown in Table 2), while Methanosarcina was not found. It means that SM1 had dominant aceticlastic methanogenic pathway but at the same time lower relative abundance of aceticlastic methanogens. The reason could be that the analysis of relative abundance of Methanosarcinales was based on DNA, and it was not directly correlated with the activity of the relevant enzymes and the expression level of the genes. On contrary, radioisotopic analysis is relatively straightforward method to detect methanogenic activity by measuring the produced 14CH4 and 14CO2 from labeled acetate. Therefore, it seems more reasonable to use radioisotopic analysis instead of the relative abundance of the corresponding microorganisms to determine the dominant methanogenic pathway. Nonetheless, when Methanosarcinales were absent (or at very low abundance (<1 %), the dominant hydrogenotrophic methanogenesis pathway can also be predicted by the taxonomic composition. In order to get a clear metabolic patterns of the microbiomes in biogas reactors, metatranscriptomic and metaproteomic analysis unveiling the expression level of the genes and relevant enzymes, respectively, is necessary in the future studies .
Variation of taxonomic and functional patterns and their correlations with environmental variables
Principal coordinate analysis (PCoA) based on Bray-Curtis distance was used to visualize the distances and variations between samples (Fig. 4). The comparison of the microbial communities among the 14 samples (Fig. 4a) indicated that feedstock and temperature played important roles in the shaping of microbial communities in biogas reactors. The sludge-based biogas reactors harbored communities, clustered clearly apart from the manure-based biogas reactors’ communities. It was consistent with different organic components in sludge (microbial cells) and manure (lignocellulose) [38, 52]. In addition, there were clear differences between the samples from thermophilic and mesophilic manure-based reactors, which was expected since different microorganisms adapted to different temperatures . However, mesophilic sample MM2 had a microbial community pattern close to those of thermophilic samples. This suggests that taxonomic patterns of mesophilic manure-based biogas reactors were more diverse compared with those of thermophilic manure-based and mesophilic sludge-based biogas reactors.
The comparison of the functional patterns among the 14 samples (Fig. 4b) showed that the sludge-based samples were clustered together and were well separated from manure-based samples. The functional pattern of sample MM2 (obtained from mesophilic manure-based biogas reactor) was closely clustered together with thermophilic manure-based biogas reactors. Further Procrustes analysis (Additional file 1: Fig S5) showed that generally functional patterns had significant correlation (p = 0.001) with taxonomic patterns by taking all the samples into consideration. The above results indicate that the functional patterns were identical to the taxonomic patterns. Therefore, one could predict the overall functional patterns from the taxonomic composition of the microbiomes, which was reported for the first time for biogas process, although similar conclusions were drawn from previous studies analyzing other ecosystems such as soil and gut microbiomes [60–62].
Canonical correspondence analysis (CCA) was used to determine the most significant environmental variables to shape the taxonomic and functional patterns (Fig. 5). Four significant environmental variables (Temperature, VFA, HRT, and free ammonia) were chosen from Table 1 on the basis of VIFs <10 . The CCA model explained 56.4 and 74.1 % of the total variance for taxonomic and functional patterns, respectively. The results showed that the selected environmental variables explained more for the variation of functional patterns compared to the variation of taxonomic patterns between samples. It should be noted that acetate accounted for more than 85 % of the total VFA as seen in Additional file 1: Table S4, and therefore the effect of VFA on taxonomic and functional patterns was mainly related with acetate. Besides, the distributions of samples in Fig. 5a, b were similar, which further suggests there is a correlation between taxonomic and functional patterns. Of the four selected environmental variables, temperature and free ammonia appears to be the most important environmental variables for both taxonomic and functional patterns. Although temperature and free ammonia were shown to affect the biogas production and microbial communities in biogas process [16, 20, 64], our study, for the first time, revealed that functional gene distributions were also determined by temperature and free ammonia.
Our study made the first step to reveal the variation of taxonomic and functional patterns of microbiomes across different biogas reactors. Moreover, the relationship between the metabolic patterns determined by metagenomic analysis and metabolic pathways determined by radioisotopic analysis was elucidated. Although the microbial composition and functional genes in the complex biogas process can be assessed with the information from metagenomic sequencing, it was not possible to identify, with certainty, the true dominant methanogenic pathways of the biogas process using only the metagenomic sequencing assessment. In addition, a clear correlation between taxonomic and functional genes patterns was demonstrated for biogas reactors. The sludge-based samples were clearly distinct from the manure-based samples for both taxonomic and functional patterns based on PCoA analysis, and temperature and free ammonia were identified to be the important environmental variables shaping both taxonomic and functional patterns.
As shown in Table 1, 14 samples were collected from 12 Danish full-scale biogas plants operated either with manure and sewage sludge as main feedstocks and at different temperatures, hydraulic retention times, etc. It should be noted that two of the biogas plants (Blåhoj and Lemvig biogas plants) have two biogas reactors running in series, and thus samples (samples names: MT2a, MT2b, MT3a, and MT3b) were collected from both steps of the series process. All the biogas plants had been running for more than 2 years under similar operational conditions. All the biogas plants reported normal operational conditions at the time of sampling, and no major changes had occurred prior to sampling, which ensured the representativeness of the samples. The samples for the microbial analysis were collected in sterile tubes (15 mL) and were frozen immediately in a cooler with dry ice. The samples for chemical analysis were collected in 0.5 L bottles and put in a cooler box with ice, while the samples for biological activities test were kept at ambient temperature. All biogas reactors had sampling points in the effluent lines close to the reactors ensuring good representative samples of the reactor biomass. The sampling valve was opened for 5 min before sample acquisition to flush the sampling valve and tube. After sampling, the samples were transported to the laboratory within 1 day. The parameters including VFA and ammonia were analyzed in our lab, while all the other parameters were obtained from the records of the specific biogas plant.
DNA extraction and metagenomic analysis
Total genomic DNA of each collected sample was extracted using QIAamp DNA Stool Mini Kit (QIAGEN, 51504) according to the manufacturer’s instructions. Libraries with insert size of 180 bp were constructed according to the manufacturer’s instructions (Illumina) for the samples. Sequencing was conducted using Illumina Hiseq 2000 platform by applying 101 bp paired-end strategy.
Sequence reads were initially filtered to remove those containing bases with quality score lower than 30 and containing one or more uncalled bases . All the pair-end reads of each dataset were joined to decrease the sequencing errors, and the reads that did not overlap were removed. The joined reads had an average length around 170 bp. To obtain a quantitative picture of the taxonomic and functional patterns, all the joined reads after sub-sampled to the same sequencing depth were uploaded to MG-RAST (Rapid Annotation using Subsystems Technology for Metagenomes) for downstream analyses with the project ID 6474 . For taxonomic analysis, the 16S rRNA gene sequences with hits were extracted from the results of BLAST analysis against the Ribosomal Database Project (RDP) database with an E-value cutoff <10−5. Only 16S rRNA fragments with at least 90 % assignment confidence were considered. The 16S rRNA gene sequences were then submitted to the RDP database (http://rdp.cme.msu.edu/) for classification with 50 % confidence. The 50 % confidence is recommended by RDP and used in many studies [6, 66]. For functional analysis, SEED subsystems and KEGG annotation were used to assign joined reads to different functional groups (SEED) or metabolic pathways (KEGG) in MG-RAST using the parameters of E-value cutoff 10−5 and minimum alignment length 50 bp, and only the annotated sequences were used for further analysis [30, 31]. Principal coordinates analysis (PCoA) and canonical correspondence analysis (CCA) were conducted by Canoco 5.0 to explore the taxonomic and functional relationships between the samples and also to identify the key environmental variables shaping the taxonomic and functional gene compositions. Procrustes analysis was performed by R (v.2.13.1; http://www.r-project.org/) with packages VEGAN. The relative abundance of 16S rRNA gene reads assigned to Genus by RDP and joined reads assigned to Subsystem Level three by MG-RAST were used for both PCoA and CCA analysis, which should yield a more conservative estimate of the distance between the samples compared to the reads assigned to species and individual genes. PCoA is the most commonly used dimensionality reduction techniques in microbial ecology to visualize the different patterns of various samples . CCA was used to explore the relationship between different patterns and environmental variables .
The methanogenic pathway of acetate degradation was determined by measuring the production of 14CH4 and 14CO2 from acetate labeled in the methyl group (C-2) . In the radioisotopic experiments, 118-mL glass batch reactors were used and 40 mL of inoculum was dispensed anaerobically under a N2/CO2 (80/20 %) headspace. For each of the 14 samples derived from the full-scale biogas reactors, the following batch reactors were included: (a) three reactors (n = 3) with inoculum only for estimation of residual methane production and (b) three reactors (n = 3) containing 47.58 ± 5.49 KBq L−1 [2-14C] sodium acetate (Amersham Pharmacia Biotech, England) for identification of the methanogenic pathway (calculation of 14CO2/14CH4 ratio). Afterward, the batch reactors were closed with butyl rubber stoppers, sealed with aluminum caps, and then incubated in a thermostatic incubator at 37 or 55 °C based on their original temperature of the biogas reactors (Table 1), until methane production ceased (around 1 month). Then, bottles containing labeled acetate were acidified (final pH = 0.95 ± 0.1) with 7.2 M HCl, resulting in the conversion of dissolved bicarbonate to CO2. The liquid and headspace of each bottle was sparged with approximately 2 L of O2, and the labeled 14CO2 was trapped with a carbon dioxide absorber for liquid scintillation counting (10 mL of CarbosorbR-E; Perkin-Elmer Company). Subsequently, the labeled 14CH4 was combusted to 14CO2 in a tube furnace above 800 °C, and the 14CO2 generated in the furnace was then trapped in 10 mL of CarbosorbR-E. For counting, the 10 mL of CarbosorbR-E were mixed with 10 mL Permaflour RE (Perkin-Elmer Company) scintillation fluid. All radioactivity measurements were performed using a liquid scintillation counter (Tri-Carb 1600; Perkin-Elmer Company).
Total ammonia was determined by the Kjeldahl method according to American Public Health Association’s Standard Methods . The concentrations of acetate, propionate, iso-butyrate, butyrate, iso-valerate, and valerate were determined by gas chromatograph (GC) (Hewlett Packard, HP5890 series II) equipped with a flame ionization detector and HP FFAP column (30 m × 0.53 mm × 1.0 μm). CH4 was analyzed by GC-TCD fitted with parallel column of 1.1 m × 3/16 "Molsieve 137 and 0.7 m × 1/4" chromosorb 108. Detailed information about the operational conditions of above GC was described previously . All analyses were made in triplicate, and the averages are presented along with the corresponding standard deviations (SD) calculated from the analyses. An analysis of variance (ANOVA) was used to test the significance of results, and p < 0.05 was considered to be statistically significant.
hydraulic retention time
volatile fatty acids
principal coordinates analysis
canonical correspondence analysis
Sundberg C, Al-Soud W, Larsson M, Alm E, Yekta S, Svensson B, et al. 454 pyrosequencing analyses of bacterial and archaeal richness in 21 full-scale biogas digesters. FEMS Microbiol Ecol. 2013;85:612–26.
Riviere D, Desvignes V, Pelletier E, Chaussonnerie S, Guermazi S, Weissenbach J, et al. Towards the definition of a core of microorganisms involved in anaerobic digestion of sludge. ISME J. 2009;3(6):700–14.
Werner JJ, Knights D, Garcia ML, Scalfone NB, Smith S, Yarasheski K, et al. Bacterial community structures are unique and resilient in full-scale bioenergy systems. Proc Nat Acad Sci USA. 2011;108(10):4158–63.
Delbès C, Moletta R, Godon J-J. Monitoring of activity dynamics of an anaerobic digester bacterial community using 16S rRNA polymerase chain reaction–single-strand conformation polymorphism analysis. Environ Microbiol. 2000;2(5):506–15.
Ye L, Zhang T, Wang TT, Fang ZW. Microbial structures, functions, and metabolic pathways in wastewater treatment bioreactors revealed using high-throughput sequencing. Environ Sci Technol. 2012;46(24):13244–52.
Luo G, Wang W, Angelidaki I. Anaerobic digestion for simultaneous sewage sludge treatment and CO biomethanation: process performance and microbial ecology. Environ Sci Technol. 2013;47(18):10685–93.
Sekiguchi Y, Kamagata Y, Harada H. Recent advances in methane fermentation technology. Curr Opin Biotechnol. 2001;12(3):277–82.
Nielsen HB, Uellendahl H, Ahring BK. Regulation and optimization of the biogas process: propionate as a key parameter. Biomass Bioenerg. 2007;31(11–12):820–30.
Heyer R, Kohrs F, Reichl U, Benndorf D. Metaproteomics of complex microbial communities in biogas plants. Microb Biotechnol. 2015;8(5):749–63.
Luo G, Angelidaki I. Analysis of bacterial communities and bacterial pathogens in a biogas plant by the combination of ethidium monoazide, PCR and Ion Torrent sequencing. Water Res. 2014;60:156–63.
Rincón B, Borja R, González JM, Portillo MC, Sáiz-Jiménez C. Influence of organic loading rate and hydraulic retention time on the performance, stability and microbial communities of one-stage anaerobic digestion of two-phase olive mill solid residue. Biochem Engin J. 2008;40(2):253–61.
Delmont TO, Prestat E, Keegan KP, Faubladier M, Robe P, Clark IM, et al. Structure, fluctuation and magnitude of a natural grassland soil metagenome. ISME J. 2012;6(9):1677–87.
Gonzalez A, Knight R. Advancing analytical algorithms and pipelines for billions of microbial sequences. Curr Opin Biotechnol. 2012;23(1):64–71.
Swindles GT, Reczuga M, Lamentowicz M, Raby CL, Turner TE, Charman DJ, et al. Ecology of testate amoebae in an amazonian peatland and development of a transfer function for palaeohydrological reconstruction. Microb Ecol. 2014;68(2):284–98.
Ariesyady HD, Ito T, Okabe S. Functional bacterial and archaeal community structures of major trophic groups in a full-scale anaerobic sludge digester. Water Res. 2007;41(7):1554–68.
De Vrieze J, Saunders AM, He Y, Fang J, Nielsen PH, Verstraete W, et al. Ammonia and temperature determine potential clustering in the anaerobic digestion microbiome. Water Res. 2015;75:312–23.
Regueiro L, Veiga P, Figueroa M, Lema JM, Carballa M. Influence of transitional states on the microbial ecology of anaerobic digesters treating solid wastes. Appl Microbiol Biotechnol. 2013;98(5):2015–27.
Regueiro L, Spirito CM, Usack JG, Hospodsky D, Werner JJ, Angenent LT. Comparing the inhibitory thresholds of dairy manure co-digesters after prolonged acclimation periods: part 2—correlations between microbiomes and environment. Water Res. 2015;87:458–66.
Karakashev D, Batstone DJ, Trably E, Angelidaki I. Acetate oxidation is the dominant methanogenic pathway from acetate in the absence of Methanosaetaceae. Appl Environ Microbiol. 2006;72(7):5138–41.
Fotidis IA, Wang H, Fiedel NR, Luo G, Karakashev DB, Angelidaki I. Bioaugmentation as a solution to increase methane production from an ammonia-rich substrate. Environ Sci Technol. 2014;48(13):7669–76.
Lu F, Hao LP, Guan DX, Qi YJ, Shao LM, He PJ. Synergetic stress of acids and ammonium on the shift in the methanogenic pathways during thermophilic anaerobic digestion of organics. Water Res. 2013;47(7):2297–306.
Becker S, Boger P, Oehlmann R, Ernst A. PCR bias in ecological analysis: a case study for quantitative Taq nuclease assays in analyses of microbial communities. Appl Environ Microbiol. 2000;66(11):4945–53.
Schluter A, Bekel T, Diaz NN, Dondrup M, Eichenlaub R, Gartemann KH, 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.
Scholz MB, Lo CC, Chain PSG. Next generation sequencing and bioinformatic bottlenecks: the current state of metagenomic data analysis. Curr Opin Biotechnol. 2012;23(1):9–15.
Tringe SG, Rubin EM. Metagenomics: dNA sequencing of environmental samples. Nat Rev Genet. 2005;6(11):805–14.
Li A, Chu YN, Wang XM, Ren LF, Yu J, Liu XL, et al. A pyrosequencing-based metagenomic study of methane-producing microbial community in solid-state biogas reactor. Biotechnol Biofuels. 2013;6:3.
Yang Y, Yu K, Xia Y, Lau FTK, Tang DTW, Fung WC, et al. Metagenomic analysis of sludge from full-scale anaerobic digesters operated in municipal wastewater treatment plants. Appl Microbiol Biotechnoly. 2014;98(12):5709–18.
Fotidis IA, Karakashev D, Angelidaki I. The dominant acetate degradation pathway/methanogenic composition in full-scale anaerobic digesters operating under different ammonia levels. Int J Environ Sci Technol. 2014;11(7):2087–94.
Jaenicke S, Ander C, Bekel T, Bisdorf R, Droge M, Gartemann KH, et al. Comparative and Joint Analysis of Two Metagenomic Datasets from a Biogas Fermenter Obtained by 454-Pyrosequencing. PLoS ONE. 2011;6(1):e14519.
Fierer N, Leff JW, Adams BJ, Nielsen UN, Bates ST, Lauber CL, et al. Cross-biome metagenomic analyses of soil microbial communities and their functional attributes. Proc Nat Acad Sci USA. 2012;109(52):21390–5.
Delmont TO, Prestat E, Keegan KP, Faubladier M, Robe P, Clark IM, et al. Structure, fluctuation and magnitude of a natural grassland soil metagenome. ISME J. 2012;6(9):1677–87.
Quaiser A, Zivanovic Y, Moreira D, Lopez-Garcia P. Comparative metagenomics of bathypelagic plankton and bottom sediment from the Sea of Marmara. ISME J. 2011;5(2):285–304.
Achenbach LA, Michaelidou U, Bruce RA, Fryman J, Coates JD. Dechloromonas agitata gen. nov., sp. nov. and Dechlorosoma suillum gen. nov., sp. nov., two novel environmentally dominant (per)chlorate-reducing bacteria and their phylogenetic position. Int J Syst Evol Microbiol. 2001;51(2):527–33.
Qiu Y-L, Hanada S, Ohashi A, Harada H, Kamagata Y, Sekiguchi Y. Syntrophorhabdus aromaticivorans gen. nov., sp. nov., the first cultured anaerobe capable of degrading phenol to acetate in obligate syntrophic associations with a hydrogenotrophic methanogen. Appl Environ Microbiol. 2008;74(7):2051–8.
Schöcke L, Schink B. Energetics of methanogenic benzoate degradation by Syntrophus gentianae in syntrophic coculture. Microbiology. 1997;143(7):2345–51.
Cai QY, Mo CH, Wu QT, Zeng QY, Katsoyiannis A. Occurrence of organic contaminants in sewage sludges from eleven wastewater treatment plants, China. Chemosphere. 2007;68(9):1751–62.
Simankova MV, Chernych NA, Osipov GA, Zavarzin GA. Halocella cellulolytic gen-nov, sp-nov, a new obligately anaerobic, halophilic, cellulolytic bacterium. Syst Appl Microbiol. 1993;16(3):385–9.
Bruni E, Jensen AP, Angelidaki I. Comparative study of mechanical, hydrothermal, chemical and enzymatic treatments of digested biofibers to improve biogas production. Bioresour Technol. 2010;101(22):8713–7.
Azman S, Khadem AF, van Lier JB, Zeeman G, Plugge CM. Presence and role of anaerobic hydrolytic microbes in conversion of lignocellulosic biomass for biogas production. Crit Rev Environ Sci Technol. 2015;45(23):2523–64.
Westerholm M, Roos S, Schnurer A. Syntrophaceticus schinkii gen nov, sp nov, an anaerobic, syntrophic acetate-oxidizing bacterium isolated from a mesophilic anaerobic filter. FEMS Microbiol Lett. 2010;309(1):100–4.
Boe K, Angelidaki I. Serial CSTR digester configuration for improving biogas production from manure. Water Res. 2009;43(1):166–72.
Solli L, Havelsrud O, Horn S, Rike A. A metagenomic study of the microbial communities in four parallel biogas reactors. Biotechnol Biofuels. 2014;7(1):146.
Zakrzewski M, Goesmann A, Jaenicke S, Junemann S, Eikmeyer F, 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.
Karakashev D, Batstone DJ, Angelidaki I. Influence of environmental conditions on methanogenic compositions in anaerobic biogas reactors. Appl Environ Microbiol. 2005;71(1):331–8.
De Vrieze J, Hennebel T, Boon N, Verstraete W. Methanosarcina: the rediscovered methanogen for heavy duty biomethanation. Bioresour Technol. 2012;112:1–9.
Ma J, Zhao B, Frear C, Zhao Q, Yu L, Li X, et al. Methanosarcina domination in anaerobic sequencing batch reactor at short hydraulic retention time. Bioresour Technol. 2013;137:41–50.
Ju F, Guo F, Ye L, Xia Y, Zhang T. Metagenomic analysis on seasonal microbial variations of activated sludge from a full-scale wastewater treatment plant over 4 years. Environ Microbiol Rep. 2014;6(1):80–9.
Pandit SN, Kolasa J, Cottenie K. Contrasts between habitat generalists and specialists: an empirical extension to the basic metacommunity framework. Ecology. 2009;90(8):2253–62.
Ahn Y-H. Sustainable nitrogen elimination biotechnologies: a review. Proc Biochem. 2006;41(8):1709–21.
Bashan LE, de-Bashan Y. Recent advances in removing phosphorus from wastewater and its future use as fertilizer (1997–2003). Water Res. 2004;38(19):4222–46.
Pérez S, Guillamón M, Barceló D. Quantitative analysis of polycyclic aromatic hydrocarbons in sewage sludge from wastewater treatment plants. J Chromatogr A. 2001;938(1–2):57–65.
Yuan HY, Chen YG, Zhang HX, Jiang S, Zhou Q, Gu GW. Improved bioproduction of short-chain fatty acids (SCFAs) from excess sludge under alkaline conditions. Environ Sci Technol. 2006;40(6):2025–9.
Whiticar MJ. Carbon and hydrogen isotope systematics of bacterial formation and oxidation of methane. Chem Geol. 1999;161(1–3):291–314.
Fricke WF, Seedorf H, Henne A, Krüer M, Liesegang H, Hedderich R, et al. The genome sequence of Methanosphaera stadtmanae reveals why this human intestinal archaeon is restricted to methanol and H2 for methane formation and ATP synthesis. J Bacteriol. 2006;188(2):642–58.
Dai Y-R, Reed DW, Millstein JH, Hartzell PL, Grahame DA, DeMoll E. Acetyl-CoA decarbonylase/synthase complex from Archaeoglobus fulgidus. Arch Microbiol. 1998;169(6):525–9.
Schnurer A, Zellner G, Svensson BH. Mesophilic syntrophic acetate oxidation during methane formation in biogas reactors. FEMS Microbiol Ecol. 1999;29(3):249–61.
Kendall MM, Boone DR. The order Methanosarcinales. The Prokaryotes. 2006: 244–56.
Simon C, Daniel R. Metagenomic analyses: past and future trends. Appl Environ Microbiol. 2011;77(4):1153–61.
Amani T, Nosrati M, Sreekrishnan TR. Anaerobic digestion from the viewpoint of microbiological, chemical, and operational aspects—a review. Environ Rev. 2010;18(NA):255–78.
Muegge BD, Kuczynski J, Knights D, Clemente JC, González A, Fontana L, et al. Diet drives convergence in gut microbiome functions across mammalian phylogeny and within humans. Science. 2011;332(6032):970–4.
Fierer N, Lauber CL, Ramirez KS, Zaneveld J, Bradford MA, Knight R. Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 2012;6(5):1007–17.
Gilbert JA, Field D, Swift P, Thomas S, Cummings D, Temperton B, et al. The taxonomic and functional diversity of microbes at a temperate coastal site: a ‘multi-omic’ study of seasonal and diel temporal variation. PLoS ONE. 2010;5(11):e15545.
Agler MT, Werner JJ, Iten LB, Dekker A, Cotta MA, Dien BS, et al. Shaping reactor microbiomes to produce the fuel precursor n-butyrate from pretreated cellulosic hydrolysates. Environ Sci Technol. 2012;46(18):10229–38.
Nielsen HB, Mladenovska Z, Westermann P, Ahring BK. Comparison of two-stage thermophilic (68 and #xB0;C/55 and #xB0;C) anaerobic digestion with one-stage thermophilic (55 & #xB0;C) digestion of cattle manure. Biotechnol Bioeng. 2004;86(3):291–300.
Meyer F, Paarmann D, D’Souza M, Olson R, Glass EM, Kubal M, et al. The metagenomics RAST server—a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics. 2008;9:386.
Ye L, Zhang T. Pathogenic bacteria in sewage treatment plants as revealed by 454 pyrosequencing. Environ Sci Technol. 2011;45(17):7173–9.
APHA. Standard methods for the examination of water and wastewater. 19th ed. New York: American Public Health Association; 1995.
Luo G, Talebnia F, Karakashev D, Xie L, Zhou Q, Angelidaki I. Enhanced bioenergy recovery from rapeseed plant in a biorefinery concept. Bioresour Technol. 2010;102(7):1310–3.
GL and IA designed the experiment. GL carried out the experiment. IF participated in the radioisotopic analysis. GL made the bioinformatics analysis and drafted the manuscript. All the authors read and approved the final manuscript.
This study was funded by the Yangfan project from Science and Technology Commission of Shanghai Municipality (14YF1400400), National Natural Science Foundation of China (51408133), and individual postdoctoral grant from The Danish Council for Independent Research (12-126632).
The authors declare that they have no competing interests.
Additional file 1. Table S1 Summary of the metagenomic sequencing data; Table S2 Taxonomic classification of the bacterial sequences (Only relative abundances of identified Class and Genus higher than 1 % were listed); Table S3 Relative abundance of each gene involved in methanogenesis pathways from metagenomic datasets of the 14 samples; Table S4 Individual VFAs in each sample; Figure S1 Relative abundances of major categories of functional genes in the shotgun metagenomes obtained from the 14 samples; Figure S2 Genes involved in carbohydrate and protein pathways from metagenomic datasets of the 14 samples; Figure S3 The ratio of 14CO2/14CH4 of the 14 samples; Figure S4 Correlation between the ratio of 14CO2/14CH4 and percentage of Methanosarcinales; Figure S5 Procrustes analyses of taxonomic and functional patterns.
About this article
Cite this article
Luo, G., Fotidis, I.A. & 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 9, 51 (2016). https://doi.org/10.1186/s13068-016-0465-6