- Open Access
Metabolic dependencies govern microbial syntrophies during methanogenesis in an anaerobic digestion ecosystem
Microbiome volume 8, Article number: 22 (2020)
Methanogenesis, a biological process mediated by complex microbial communities, has attracted great attention due to its contribution to global warming and potential in biotechnological applications. The current study unveiled the core microbial methanogenic metabolisms in anaerobic vessel ecosystems by applying combined genome-centric metagenomics and metatranscriptomics. Here, we demonstrate that an enriched natural system, fueled only with acetate, could support a bacteria-dominated microbiota employing a multi-trophic methanogenic process. Moreover, significant changes, in terms of microbial structure and function, were recorded after the system was supplemented with additional H2. Methanosarcina thermophila, the predominant methanogen prior to H2 addition, simultaneously performed acetoclastic, hydrogenotrophic, and methylotrophic methanogenesis. The methanogenic pattern changed after the addition of H2, which immediately stimulated Methanomicrobia-activity and was followed by a slow enrichment of Methanobacteria members. Interestingly, the essential genes involved in the Wood-Ljungdahl pathway were not expressed in bacterial members. The high expression of a glycine cleavage system indicated the activation of alternative metabolic pathways for acetate metabolism, which were reconstructed in the most abundant bacterial genomes. Moreover, as evidenced by predicted auxotrophies, we propose that specific microbes of the community were forming symbiotic relationships, thus reducing the biosynthetic burden of individual members. These results provide new information that will facilitate future microbial ecology studies of interspecies competition and symbiosis in methanogenic niches.
Microbial methanogenic metabolism is considered as one of the oldest bio-activities on earth and draws great attention because of its global warming potential  , which is 28 times higher than carbon dioxide (CO2) on a 100-year horizon . In a natural ecosystem, around one billion tons of methane (CH4) is formed through microbial activity as an intermediate step of the global carbon cycle . Nevertheless, an enhanced and well-controlled methanogenic process has implications for energy generation  due to its high calorific value. Microbial methanation was extensively employed in vessel ecosystems, i.e. biogas reactors, to attain large-scale production as a sustainable energy source. It is postulated that methanogenesis is performed mainly through acetoclastic, hydrogenotrophic, and secondary through methylotrophic pathways in oxygen-depleted environments. The known methanogenic members belong mainly to phylum Euryarchaeota, with few exceptions, which were recently assigned to candidate phyla “Bathyarchaeota”  and “Verstraetearchaeota” . All methanogens are physiologically specialized and able to scavenge the electrons from hydrogen (H2), formate, methanol, and acetate, having CH4 as the final product. Archaeal growth and activity can create ecological niches for the oxidizing (H2 producing) bacteria, and form syntrophic relations in a complex community.
In the past years, genome-centric metagenomics was extensively used to describe complex syntrophic microbial communities, and successfully revealed essential knowledge regarding the microbial functions of the keystone species mainly based on their gene profiles [7, 8] . The majority of studies regarding methanogenic process were focused on specific microbes contributing to the degradation of recalcitrant substrates  or the involvement of rare taxon in the methanogenic process [10, 11], while few attempts have been made to underlie the basic mechanisms of microbial
community assembly and function [7, 12, 13]. In natural ecosystems, the holistic untangling of the intricate methanogenic process was hampered by the inextricable influence of numerous environmental variables occurring simultaneously. Moreover, the in-situ activity of the individual members in microbial communities and the ecological relationships existing among microbes were extremely difficult to elucidate during the digestion of complexed substrates. Thus, simplified model systems are required to unveil the fundamental metabolic insights into methanogenic activities. A previous study dissected the complex methanogenic consortium into tractable model sections by substrates specification in continuous reactor operation and successfully assigned putative functional roles to the de-novo reconstructed genomes . However, a crucial limitation of studies based solely on metagenomic surveys is the lack of direct evidence for the activity of individual microbes. Therefore, other –omics approaches and advanced molecular tools, such as transcriptomics, proteomics, metabolomics, and stable isotope labelling were gradually introduced to analyze the microbial activity during the methanogenic process [14,15,16,17].
The current study is dedicated to unveil the core microbial methanogenic metabolisms with combined genome-centric metagenomic and metatranscriptomic strategies. The methanogenic metabolism was favoured in microcosms where the microbial communities were simplified by providing a chemically-defined substrate (acetate). The study revealed the in situ activity of methanogens in syntrophic microbial communities and their affinity to H2 provision. Moreover, this work also provided mechanistic understandings of the bacterial functionalities both for acetate oxidation and revealed important auxotrophic dependencies, as well as community structure maintenance during methanogenesis.
Materials and methods
The triplicate lab-scale biogas continuous stirred-tank reactors (working volume 1.8L) were inoculated with digestate from full-scale thermophilic biogas plant (Snertinge, Denmark). The plant was fed with 70–90% animal manure and 10–30% food industrial organic waste; therefore, the inoculum provided the microbial community to adopt to heterogeneous substrate degradation. During the experiment, the reactors were fed with synthetic medium, in which only acetic acid was supplied as an organic carbon source. Other nutrients were provided by basal anaerobic medium . The reactors operated under thermophilic condition (55 °C) and the operational parameters were chosen according to empirical experiences of highly efficient thermophilic biogas reactors, i.e. the organic loading rate was 1g acetic acid/day. L-reactor and the hydraulic retention time was 15 days. The reactors were fed four times per day with peristaltic pumps to achieve the desired organic loading rate and HRT. Once the reactors reached the steady state, H2 gas was supplemented to each reactor with two stainless steel diffusers (pore size 2 μm) at the rate of 1 mL/min. To ensure efficient H2 utilization, the gas phase of the reactors was constantly recirculated into liquid phase with peristaltic pumps. Throughout the experiment, biogas production was recorded with water-replacement gas metres; biogas composition was measured using a gas chromatograph (Mikrolab, Aarhus A/S, Denmark), equipped with a thermal conductivity detector (TCD). The volatile fatty acids and ethanol were measured with a gas chromatograph (Shimadzu GC-2010 AF, Kyoto, Japan), equipped with a flame ionization detector (FID) . Biomass formation was estimated through volatile suspended solids measured according to the Standard Methods for the Examination of Water and Wastewater . All determinations and measurements were done in triplicate samples.
Sample collection and sequencing
The liquid samples were acquired from the triplicate reactors before, 18 hours after and 36 days after H2 addition (Sample point 1, 2, and 3, respectively). For all the samples, the genomic DNA was extracted with PowerSoil® DNA Isolation Kit and the total RNA was extracted PowerMicrobiome® RNA Isolation Kit (Mo Bio Laboratories, Inc., Carlsbad, USA). All the extractions were performed with additional phenol cleaning steps in order to improve the quality of the extractives. The ribosome RNA was removed from total RNA samples with Ribo-Zero® rRNA Removal Kit (Bacteria) (Illumina, San Diego, USA). The DNA and RNA samples were sent to Ramaciotti Centre for Genomics (UNSW, Sydney, Australia) for cDNA construction, library preparation, and sequencing (Illumina NextSeq).
The DNA sequences from 9 samples were filtered with Trimmomatic software , co-assembled with metaSPAdes , and atomically binned with MetaBAT . The quality of the metagenome-assembled genomes (MAGs) were examined with CheckM  and evaluated with a MAG quality standard developed by Genomic Standards Consortium . The average nucleotide identity analysis (ANI) was performed against all the genomes that were deposited in NCBI Reference Sequence Database . The genomes hits with ANI higher than 97% were used to classify the MAGs at the species level [27, 28]. The putative taxonomy classification of the unclassified MAGs were further assessed based on ubiquitous proteins with PhyloPhlAn . The genes of the co-assembled metagenome were predicted and annotated with Integrated Microbial Genomes & Microbiomes (IMG) . For more comprehensive methanogenic pathway reconstruction, all archaeal MAGs were resubmitted to IMG as genomes assembled from metagenome for gene prediction and annotation.
The microbial communities were profiled through reads recruitment from the sequencing samples. The average coverage of MAGs in each metagenome sample was calculated based on the number of reads aligned by Bowtie2  and the detailed procedure were described by Campanaro et al. (2016) . The relative abundance of the MAGs in a community was determined with the average coverage of the MAG in one metagenome sequencing sample, according to:
The sequenced transcripts were aligned to assembled metagenomes with Bowtie2 and quantified with HTSeq-count [31, 33]. Therefore, instead of de-novo assembling the RNA sequences, the metatranscriptomes inherited the annotation from the corresponding metagenomes. The expression level of genes was evaluated by reads per reads per kilobase of exon model per million mapped reads (RPKM) .For comparison purpose, the RPKM numbers were normalized considering the expression level of methyl-coenzyme M reductase gene (subunit alpha) and CH4 production rates of the reactors during the time that each sample was collected. Moreover, the metatranscriptomes were dissected according to the binning results in order to generate individual expression profiles for MAGs. The overall activity of a MAG was evaluated by the average gene RPKM within the genome. The relative activity of MAGs in a community was measured according to a similar formula as relative abundance:
The comparison between the relative abundance and activity of a MAG suggested its activity level. More specifically, a low abundance/ activity ratio represented an active member, who undertook many microbial metabolisms with few numbers of cells. The overall microbial community composition was visualized by Anvi’o .
The gene expression profiles derived from genome-dissected metatranscriptome were also used to indicate the functional role of each MAG during the methanogenic process. Genes were categorized based on KEGG modules and the average RPKM of all the genes was calculated for each module. The differential expression of each gene in coordination to the H2 addition was examined using edgeR package . Other statistical tests (Student’s t-tests and correlation tests) were performed with Excel.
Specific metabolisms, i.e. methanogenesis and acetate uptake, were tentatively distributed to individual MAGs based on the expression level of signature genes. For instance, the methanogenic activity of individual archaeal MAGs was determined based on the expression level of MAG-specific mcrA comparing to the overall expression of all mcrA . Moreover, MAG-specific acetate kinase (ack) as well as acyl-CoA synthetase (acs) were used to correlate the acetate metabolism to individual members among the microbial community.
The raw sequence data were deposited on sequence read archive with accession no PRJNA525781, The biosample metadata were deposited in Genomes OnLine Database (GOLD) as study Gs0128993.The metagenome annotation was deposited as analysis project Ga0214976. The annotation of methanogen MAGs was deposited as analysis projects Ga0214977, Ga0214981, Ga0214989, Ga0214990 and Ga0214991.
Methanogenic microcosms enriched by acetate and H2
The tractable low-complexity methanogenic microbial communities were obtained from triplicate lab-scale continuous biogas reactors operated under thermophilic conditions providing acetate as the only organic carbon source. After establishment of stable conditions, external H2 gas was injected into all reactors with stainless steel diffusers to assess the microbes’ affinity to H2 partial pressure. During the entire experimental operation, the pH in each reactor was self-stabilized within the optimal range of methanogenesis (7-7.5). The triplicate reactors performed consistently during the two steady operational conditions (Sample Points 1 and 3). Nevertheless, a significant discrepancy among triplicate reactors was observed during the transition before and after H2 addition (Sample Point 2), which was mainly attributed to the instability of the microbial community adaptation process. CH4, along with inorganic carbon compounds including CO2, bicarbonate (HCO3-) and carbonate (CO32-), were the main digestion products. The digestion profiles are described as mol of carbon contained in each products (Fig. 1). In addition, approximately 4% of carbon (mol of carbon in biomass / mol as carbon in acetate) was used to build microbial biomass. The methanation process was extremely efficient as, less than 0.5% of the carbon (mol of carbon in acetate / mol as carbon in acetate) left the system as undigested acetate. All of the injected H2 was consumed and the CH4 yield significantly increased from 299.8±4.4 mL/g acetate to 409.3±14.6 mL/g acetate (Fig. 1, Additional file 1).
Metagenome-assembled genome reconstruction and taxonomy assignment
Samples for shotgun sequencing were collected from triplicate reactors at three time points:1) before H2 addition, 2) 18 hours after H2 addition, and 3) 36 days after H2 addition. Point 1 and Point 3 were chosen during the operational steady states where the CH4 production rate of each reactor varied less than 10% for 10 consecutive days. Sequences from all samples were co-assembled and automatically binned in order to extract metagenome assembled genomes (MAGs). In total, 79 MAGs were extracted with total coverage ranging from 5× to 9595× (Fig. 2, Additional File 1). According to the completeness and contamination values assessed by CheckM, and the quality standards developed by Genomic Standards Consortium , 36 MAGs were assigned to the “high-quality” group, 25 MAGs to the “medium quality” group, and 18 MAGs to the “low quality” group. It is noteworthy that over 95% of total shotgun sequences including DNA and RNA in all samples, could be aligned to the medium-high quality MAGs, suggesting that the majority of the microbial diversity has been recovered with the assembly and binning process.
Five nearly complete archaeal MAGs were present in the community (> 98% complete) , four out of which were characterized at species level as Methanothermobacter thermautotrophicus [38, 39] DTU592, Methanosarcina thermophila  DTU593, Methanoculleus thermophilus  DTU608 and Methanobacterium sp. MB 1 DTU624 (Additional File 2). The unknown archaeal MAG (unclassified Methanomicrobia DTU639), which was also previously found in biogas reactors , could be assigned to a member of class Methanomicrobia based on tentative phylogenetic classification. In contrast to archaeal MAGs, 51 out of 56 bacterial MAGs could only be classified at the family or higher taxonomic level. According to the relative abundance, more than 90% of the microbial community could be represented by 18 most abundant MAGs, 15 of which belonged to domain Bacteria. Among bacteria, 11 MAGs were assigned to Firmicutes (4 MAGs), Bacteroidetes (2 MAGs) Synergistetes (2 MAGs), Proteobacteria (2 MAGs), and Thermotogae (1 MAG). The remaining four MAGs were unclassified Bacteria spp. (Additional File 3).
Microbial community composition and transcriptional activity
The microbial community composition and the transcriptional activity profiles were determined using the average genome coverage of each MAG and the average gene expression level (reads per kb per million mapped reads, RPKM) of all protein-coding genes in each MAG (Fig. 3, Additional File 4 and 5 ). Interestingly, a robust bacterial activity was observed in the reactor, although the present acetoclastic methanogens (M. thermophila DTU593) could theoretically undertake the majority of the acetate methanation process. In fact, the five methanogens constituted only a small part of the total microbial community, which is 19–37% of the abundance and 7–27% of the activity. It was surprising that methanogens constituted the minority both in respect to relative abundance and activity since only methanogenic substrates (acetate and H2-CO2) were fuelling the process. Before H2 addition, the most abundant MAG (DTU593) among the entire microbial community was classified as Methanosarcina thermophila, accounting for 17% of the total community (Additional File 4). However, its activity was calculated as 5.3% among the entire microbial community (Additional File 5). The relatively high RNA/DNA ratio indicated a high cellular protein synthesis potential of M. thermophila DTU593 , suggesting a possible high growth rate under this condition [44, 45]. In contrast, unclassified Bacteria sp. (DTU645) and unclassified Synergistaceae sp. (DTU638), which were the second and third most abundant MAGs (accounting for 12% of the community each), were responsible for 19% and 18% of the activity respectively (Additional File 4 and 5). After H2 addition, the microbial abundance (based on genome coverage) and transcriptional activity (based on average gene RPKM) profiles changed significantly as a result of community adaptation .The overall archaeal activity changes correlated with the CH4 production rate of the reactor in steady state (R2=0.84), whereas the correlation between the overall archaeal relative abundance and steady state CH4 production was lower (R2=0.53) (Additional File 5). In addition to methanogenic archaea, the supplemented H2 also reshaped the bacterial community. The most significant change was the increase of Coprothermobacter proteolyticus DTU632, which became the most abundant MAG, accounting for 19% of the total community. Interestingly, C. proteolyticus DTU632 only contributed to 6.8% of activity, which was lower than the hydrogenotrophic methanogen M. thermophilus DTU608 (18.6%) and unclassified Bacteria sp. DTU645 (7.2%) (Additional files 4 and 5).
Metabolism of the methanogens
Acetate, the only organic carbon source supplied to the reactors, was taken up by the microbes through two pathways: inversed phosphotransacetylase-acetate kinase pathway (PTA-ACKA) and AMP-forming acyl-CoA synthetase pathway (AMP-ACS). Thus, acetate utilization by individual microbes could be estimated according to the expression level of MAG-specific acyl-CoA synthetase (acs) and acetate kinase (ack) genes (Additional File 7). The main archaeal acetate consumer was M. thermophila DTU593, which consumed less than 50% of the acetate supplied to the reactors through inversed PTA-ACKA before H2 addition. After H2 addition, the expression level of M. thermophila-specific ack was decreased significantly, while expression of bacterial ack was increased. Considering that CH4 was produced by the five archaea, the methanogenic activity was tentatively distributed among them (expressed as %) based on the expression level of MAG-specific mcrA (Fig. 2, Additional File 8). Before H2 addition, the methanogenic activity was highest in M. thermophila DTU593 (86%) and M. thermophilus DTU608 (11%). After reaching the steady state, external H2 gas was supplied in order to trigger a metabolic shift towards hydrogenotrophic methanogenesis. The amount of H2 injected into the reactor was chosen stoichiometrically to reduce half of the CO2 that was produced from acetate during the methanogenic process. The addition of external H2 gas changed the methanogenic activity of archaeal MAGs. Specifically, the activity of M. thermophilus DTU608 was significantly enhanced and it became the main methanogenic contributor in the microbial community shortly after the H2 addition (98% of methanogenic activity). By the end of the experiment, the reactor stabilized in a new operational steady state, during which 56 ± 0.4% of acetate-carbon was converted to CH4 (16% more compared with previous states). After long-term adaptation to the H2 addition, although M. thermophilus DTU608 maintained its dominance (71%), a small but significant fraction of methanogenic activity was taken over by M. thermophila DTU593 (9%) and other hydrogenotrophic methanogens (15%).
Pathways related to methanogenesis and relevant energy conservation systems were reconstructed in all archaeal MAGs (Fig. 4). The expression levels of those genes (normalized according to the expression level of mcrA gene and CH4 production rate) were examined before, shortly after and long after H2 addition (Additional File 8). M. thermophila DTU593 expressed all the genes involved in hydrogenotrophic, methylotrophic and acetoclastic methanogenesis, indicating its multi-trophic functional role in anaerobic digestion. In particular, the expression of methylamine/methanol-specific coenzyme M methyltransferase genes (mta, mtb) suggested a considerable contribution of methylotrophic pathways (Additional File 8). For energy conservation, M. thermophila DTU593 obtained the electron from intermediate H2 through methanophenazine-reducing hydrogenase (Vho), coenzyme F420-reducing hydrogenase (Frh), and Escherichia coli hydrogenase 3 (Ech). The electrons provided by Ech were transferred to ferredoxin and used for CO2 reduction. The electrons carried by F420H2 were not only used for CHO-H4MPT reduction in hydrogenotrophic methanogenesis but also transferred to methanophenazine through F420H2 dehydrogenase (Fpo). Finally, methanophenazine reduced by Fpo and Vho transferred the electrons to CoM-S-S-CoB through the membrane-bound heterodisulfide reductase (hdr DE).
In contrast, the methanogenic activity of M. thermophilus DTU608 was restricted to hydrogenotrophic pathways (Fig. 4). M. thermophilus DTU608 lacked cytochromes but possessed the electron bifurcation system, which allows coupling CO2 reduction and CoM-S-S-CoB reduction with Mvh-Hdr complex oxidation. In addition, M. thermautotrophicus DTU592 and Methanobacterium sp. gradually increased their relative abundance and activity only after long term operation. These two ‘slowly emerged’ archaeal MAGs encoded core hydrogenotrophic methanogenesis pathways similar to M. thermophilus DTU608; however, they possessed different hydrogenase complexes for H2 uptake. An important difference is that both M. thermautotrophicus DTU592 and Methanobacterium sp. DTU624 used energy converting hydrogenase B (Ehb) instead of Ech, which was present in M. thermophilus DTU608 (Fig. 4, Additional file 9).
Ehb was proven to be related to autotrophic CO2 assimilation, which could confer an advantage to Methanobacteriaceae spp. by increasing its relative abundance in the microbial community during long-term organic carbon starvation . Moreover, both M. thermautotrophicus DTU592 and Methanobacterium sp. DTU624 significantly upregulated carbon monoxide dehydrogenase (coo) and acetyl-CoA decarbonylase/synthase (Cdh) genes, supporting carbon fixation activity for autotrophic growth (Fig. 4, Additional File 8), while M. thermophilus DTU608 relied on external acetate (heterotrophic) as indicated by the expression of NDP forming acyl-CoA synthetase genes.
Metabolism of the bacteria
More than 50% of the acetate, which was not taken up by archaea, was metabolized by bacterial members in the community (Additional File 7). According to the transcriptional activity of formyltetrahydrofolate synthetase gene (fhs), about half of the bacteria community (31 out of 79) were indicated to have a syntrophic lifestyle in association with methanogenic Archaea [47, 48]. Interestingly, the genes encoding the enzyme to directly break the bonds between the carbonyl and methyl branch in the acetate (acetyl-CoA decarbonylase , cdh) were expressed at an extremely low level (not significant according to edgeR normalization) in all bacterial MAGs (Additional File 10). This result indicated that these bacteria may have acetate utilization pathways other than the conventional Wood-Ljungdahl (WL) pathway, similar to the alternative pathway previously proposed in Thermotogae spp. . High expression levels of ack and glycine decarboxylase genes (gcvP) were found in MAGs belonging to diverse taxa, indicating that the glycine cleavage system proposed for the Thermotogae phylum might be more widely used for bacterial acetate utilization (Additional File 10). Bacteria adopted versatile strategies to transform acetate to glycine, to further metabolize the intermediates released from the glycine cleavage system, as well as to conserve energy. The detailed acetate degradation pathways were proposed based on highly expressed genes (50% quantile among all the genes expressed in the genome) in the most abundant acetate consuming bacterial MAGs (Fig. 5, Additional File 11, 12 and 13).
The Thermotogae-like acetate utilization pathway was reconstructed in unclassified Bacteria sp. DTU645, which phylogenetically clustered with Firmicutes (Fig. 1). The high expression level of glycine reductase gene (grd) in unclassified Bacteria sp. DTU645 suggested an alternative path for acetate to enter the glycine cleavage system, where acetyl phosphate was directly converted to glycine through reversed glycine reduction (Additional File 11). The glycine reductase gene was found highly expressed in many other bacterial acetate utilizers, e.g. unclassified Synergistaceae sp. DTU638 and C. proteolyticus DTU632 (Additional File 12 and 13). The glycine cleavage system catalysed the decarboxylation of glycine and released methylene-tetrahydrofolate, NH3 and CO2. The methylene-tetrahydrofolate could be further oxidized through a partial reversed WL pathway in unclassified Synergistaceae sp. DTU638 and unclassified Bacteria sp. DTU645, having CO2 as the final product (Additional File 11 and 12). Interestingly, the gene set mediating methylene-tetrahydrofolate oxidation was completely absent in C. proteolyticus DTU632, whose genome was 100% complete according to CheckM (Additional File 2). C. proteolyticus was previously characterized as a proteolytic bacterium that produces acetate, CO2 and H2 as main fermentative products . However, its high relative abundance and activity in this study indicated its involvement in the acetate metabolism with additional H2 supplements. Considering C. proteolyticus is known to metabolize amino acids, a Stickland-like reaction was tentatively reconstructed in strain DTU632 after considering the highly expressed genes (Fig. 5, Additional file 13). Specifically, the methylene-tetrahydrofolate released from the glycine cleavage system was combined with another glycine to create serine, and eventually entered the pathways for amino acid metabolism. In all the proposed pathways, the electrons were balanced from acetate-uptake to glycine decarboxylation, and additional electron disposal was required for further oxidation of methylene-tetrahydrofolate. For unclassified Synergistaceae sp. DTU638, the electron was disposed of as H2 as suggested by the high expression of membrane-bound hydrogenase and Fe-S-cluster-containing hydrogenase (Additional File 12), which explained why H2 inhibited the activity of unclassified Synergistaceae sp. DTU638. As a consequence of external H2 addition, the acetate uptake activity was taken over by unclassified Bacteria spp. Unlike unclassified Synergistaceae sp. DTU638, the formate dehydrogenase operon of unclassified Bacteria sp. DTU645 did not contain a hydrogenase gene (Additional File 11), suggesting the electrons could be disposed of in other forms than H2. This observation explained the increment in relative abundance of DTU645 and other bacterial members after H2 addition.
Overall metabolism of microbial community
In order to maintain the methanogenic activity of the microbial community, a syntrophic behaviour is needed to synthesize numerous metabolites. The holistic microbial community activity could be evaluated by the average RPKM of genes in each KEGG module. An overall shift of the microbial activity was observed in the majority of the KEGG modules after H2 addition. Specifically, the expression level of the KEGG modules related to methanogenesis, including both reactions directly linked to CH4 formation and biosynthesis of cofactors (F420) increased roughly 1.5-fold after H2 addition (Additional File 14). Moreover, H2 also enhanced the activity of the glyoxylate cycle and the biosynthesis of lipids and specific amino acids (Additional File 14).
Although both abundance and activity of individual MAGs changed significantly in the different H2 adaptation stages, ubiquitous metabolic pathways were found to be essential for maintaining the complex microbial community. Specifically, the results showed that the core metabolisms carried out by the dominant bacteria community before H2 addition could be maintained by other members proliferating after H2 addition (Additional File 14). These metabolic steps were catalysed by proteins encoded by constitutively expressed genes that maintain basic cellular function, such as biosynthesis, energy conservation, repair, and regulatory systems.
The investigation of each individual MAGs’ expression profile showed that the biosynthesis of common cofactors such as coenzyme A, NAD and riboflavin were evenly expressed in the dominant microbes, whereas the biosynthesis of several energy-expensive amino acids  and cofactors (such as biotin) were absent from some MAGs (Fig. 6). Specific metabolic traits are suggested for individual microbes based on their gene expression profile. For example, C. proteolyticus DTU632 lacked efficient pathways for electron disposal and energy-efficient acetate catabolism but showed high activity of reductive citric acid cycle. Therefore, C. proteolyticus might grow as energy-expensive amino acid auxotrophs to reduce the biosynthetic burden. The high expression of many amino acid transport systems indicated that the growth of C. proteolyticus DTU632 was supported by external amino acid uptake, such as tryptophan, tyrosine, and cysteine, (Additional File 13). The expression profiles of unclassified Bacteria sp. DTU628, unclassified Clostridiales sp. DTU630, and unclassified Rhodocyclaceae sp. DTU583 implied that these bacteria could synthesize relevant amino acids during H2 addition (Fig. 6). Another interesting observation was that the genes involved in the biosynthesis of biotin were only found in unclassified Clostridiaceae sp. DTU570, unclassified Gammaproteobacteria sp. DTU594, and unclassified Clostridiales sp. DTU630. It was previously demonstrated that the growth of some methanogens required an external supply of biotin . Considering the consistent expression of genes encoding biotin-specific transporters in the methanogens, biotin auxotrophy might have forced methanogens to scavenge metabolic products for methanogenesis, thereby leading to syntrophic behaviour between bacteria and archaea.
The combination of genome-centric metagenomics and metatranscriptomics successfully revealed individual functional roles of microbial members in methanogenic microcosms. The results assigned a multi-trophic role to Methanosarcina thermophila, suggesting its ability to perform simultaneous methanogenesis from acetate, CO2 and methanol/ methylamine. Although the use of cytochromes in M. thermophila would impose thermodynamic limitations to compete for H2 during low H2 partial pressure , Fpo-Hdr mediated heterodisulfide reduction promoted the activity of Frh, leading to the activation of the hydrogenotrophic pathway. Therefore, the H2 produced as an intermediate during anaerobic digestion not only promote the growth of hydrogenotrophic methanogens but also provide a favourable ecological niche for M. thermophila. The complex association between acetoclastic methanogens and acetate-oxidising bacteria could be one cause of functional redundancy in microbial communities involved in biomethanation. In this experiment, although M. thermophila had the metabolic potential to perform methanogenesis through a mainly acetoclastic pathway, a bacteria-dominated microbial consortium was formed, which resulted in a multi-trophic methanogenesis strategy. The results also underlined the importance of methanol/methylamine methanogenic pathways, although significant methanol concentrations were not detected during the process. In fact, the methanol/ methylamine-specific methanogens were previously identified as pivotal members in many other biogas reactors fed with manure . From this result, we believe that the maintenance of the relevant metabolites (such as methanol/methylamine) at low concentration in an efficient anaerobic digestion system. The addition of external H2 greatly enhanced the activity and the relative abundance of hydrogenotrophic methanogens, including M. thermophilus, whose activity was inherent in the microbial community before H2 addition and Methanobacteriaceae spp, which was nonexistent before H2 addition but significantly increased later after a long period of adaptation to external H2. The stimulation of Methanomicrobia members was in accordance with previous research, where anaerobic digestion systems were exposed to different H2 partial pressures [53,54,55]. For instance, a study on biogas biological upgrading systems  concluded that the microbial community would turn over from a “Methanoculleus-dominated” microbial community to a “Methanothermobacter-dominated” community after a 2-year stable operation with external supplemented H2. Several hypotheses were proposed to explain the methanogens differing affinities to H2 concentration; these hypotheses were based on gene expression regulation, or considered energy conservation strategies and syntrophic associations with bacterial partners . This work suggests that the competition among the different hydrogenotrophic methanogens can be explained by a bargain between methanogenic activity and autotrophic growth. The hydrogenase genes encoded by Methanobacteriaceae spp. (ehb) might especially support growth with external H2 and promote growth during long-term H2 adaptation and limited carbon sources.
The bacterial metabolic pathways were essential for their contribution to acetate oxidation, as well as for their role in maintaining the microbial community structure. Bacterial acetate oxidation under anaerobic conditions is postulated to be performed through the reductive WL pathway, which was annotated in the known syntrophic oxidizer, Syntrophaceticus schinkii . However, it was also recently found that many genomes from acetate utilizers, including both MAGs and isolates, possessed only a subset of WL pathway genes [49, 57]. The results of this study showed extensive use of a glycine cleavage system by many members of the community to circumvent the direct break of the carbonyl and methyl carbon bonds of acetate. The glycine cleavage system could be used in the previously proposed manner, where it was combined with a partial WL pathway to produce CO2/H2 and support the syntrophic activity with hydrogenotrophic methanogens . Moreover, in the present study, a completely new Stickland-like path was proposed for C. proteolyticus DTU632. Unlike the conventional Stickland reaction where the amino acids were provided to the microbes as a carbon source, in this newly proposed pathway, acetate was converted to glycine which served as both an electron donor and acceptor for further metabolism. The oxidation of carbonyl groups was performed through the glycine cleavage system and the methyl carbons were used for amino acids biosynthesis as previously suggested in organohalide-respiring Dehalococcoides mccartyi . C. proteolyticus’s capability to utilize acetate was not revealed in studies performed on pure cultures  and this metabolic trait might only be activated under specific conditions. The current experiment imposed a selective pressure on C. proteolyticus, where acetate was the sole organic source, external H2 was supplemented, and microbial partners were present to form syntrophic relationships. This finding encourages future studies to explore metabolic potential in diverse environments and to prove that the functional roles of individual members of a microbial community could go beyond the physiological characterization of the corresponding isolates. Lastly, some crucial transcriptomic activities, such as biosynthesize of amino acids and cofactors, were absent in the most abundant MAGs, which indicated potential exchange of carbon sources, amino acids, and cofactors among bacterial and archaeal members. These results underlined the importance of auxotrophy in the microbial communities, which was previously proposed to reduce biosynthetic burden [60, 61]. This finding may be considered one of the most important reasons for maintaining/forming a complex microbial community even during growth on simple substrates (e.g. acetate). Auxotrophy could also provide explanations for previous observations, as for example, the unexpected proteolytic activity of C. proteolyticus , which was observed during cellulose degradation (without protein as substrate), and required an external source of energy-expensive amino acids.
The combined genome-centric metagenomics and metatranscriptomics strategy used in the present work was extremely informative to characterize unknown microbial communities and elucidate the metabolic activity of individual microbial species. Especially, the distribution of metabolic activities based on genome-dissected metatranscriptomes directly revealed the contribution of individual MAGs to the global activity of the microbiome. The novel microbial insights illustrated in the current study expanded the current knowledge regarding metabolisms in methanogenic systems and the results obtained can open new horizons for future microbial ecology studies of interspecies competition or symbiosis.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the sequence read archive (SRA, https://www.ncbi.nlm.nih.gov/sra) and Genomes OnLine Database (GOLD, https://gold.jgi.doe.gov/) and Integrated Microbial Genomes & Microbiomes (IMG, https://img.jgi.doe.gov/). Additional data are all provided as Supplementary Datasets in Additional files.
- ack :
Acetate kinase gene
- acs :
Acyl-CoA synthetase gene
AMP-forming acyl-CoA synthetase pathway
- Cdh :
- coo :
Carbon monoxide dehydrogenase gene
- Ech :
Energy converting hydrogenase
- Ehb :
Energy-converting hydrogenase B
- Fpo :
- Frh :
Coenzyme F420-reducing hydrogenase
- gcv :
Glycine decarboxylase genes
- grd :
Glycine reductase gene
- hdr :
Membrane-bound heterodisulfide reductase
Metagenome assembled genome
- mcrA :
Methyl coenzyme M reductase gene
- mta, mtb :
Coenzyme M methyltransferase genes
Phosphotransacetylase-acetate kinase pathway
Reads per kilobase of exon model per million mapped reads
- Vho :
Sorokin DY, Makarova KS, Abbas B, Ferrer M, Golyshin PN, Galinski EA, et al. Discovery of extremely halophilic, methyl-reducing euryarchaea provides insights into the evolutionary origin of methanogenesis. Nat Microbiol. Nature Publishing Group. 2017;2:17081.
Myhre G, Shindell D, Bréon F-M, Collins W, Fuglestvedt J, Huang J, et al. Anthropogenic and natural radiative forcing. Clim Chang. 2013;423:658–740.
Thauer RK, Shima S. Biogeochemistry: Methane and microbes. Nature. Nature Publishing Group; 2006;440:878.
Angelidaki I, Karakashev D, Batstone DJ, Plugge CM, Stams AJM. Biomethanation and its potential. Methods Enzymol. Elsevier; 2011. p. 327–51.
Evans PN, Parks DH, Chadwick GL, Robbins SJ, Orphan VJ, Golding SD, et al. Methane metabolism in the archaeal phylum Bathyarchaeota revealed by genome-centric metagenomics. Science (80- ). American Association for the Advancement of Science; 2015;350:434–438.
Vanwonterghem I, Evans PN, Parks DH, Jensen PD, Woodcroft BJ, Hugenholtz P, et al. Methylotrophic methanogenesis discovered in the archaeal phylum Verstraetearchaeota. Nat Microbiol. Nature Publishing Group; 2016;1:16170.
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. Elsevier. 2018;140:123–34.
Kougias P, Campanaro S, Treu L, Tsapekos P, Armani A, Angelidaki I. Genome-centric metagenomics revealed the spatial distribution and the diverse metabolic functions of lignocellulose degrading uncultured bacteria. bioRxiv. Cold Spring Harbor Laboratory; 2018;328989.
Lykidis A, Chen C-L, Tringe SG, McHardy AC, Copeland A, Kyrpides NC, et al. Multiple syntrophic interactions in a terephthalate-degrading methanogenic consortium. ISME J. Nature Publishing Group; 2011;5:122.
Kirkegaard RH, Dueholm MS, McIlroy SJ, Nierychlo M, Karst SM, Albertsen M, et al. Genomic insights into members of the candidate phylum Hyd24-12 common in mesophilic anaerobic digesters. ISME J: Nature Publishing Group; 2016.
Stolze Y, Bremges A, Rumming M, Henke C, Maus I, Pühler A, et al. Identification and genome reconstruction of abundant distinct taxa in microbiomes from one thermophilic and three mesophilic production-scale biogas plants. Biotechnol Biofuels. BioMed Central; 2016;9:156.
Zhu X, Campanaro S, Treu L, Kougias PG, Angelidaki I. Novel ecological insights and functional roles during anaerobic digestion of saccharides unveiled by genome-centric metagenomics. Water Res. 2019;151.
Narihiro T, Nobu MK, Kim N, Kamagata Y, Liu W. The nexus of syntrophy-associated microbiota in anaerobic digestion revealed by long-term enrichment and community survey. Environ Microbiol. Wiley Online Library. 2015;17:1707–20.
Fontana A, Kougias PG, Treu L, Kovalovszki A, Valle G, Cappa F, et al. Microbial activity response to hydrogen injection in thermophilic anaerobic digesters revealed by genome-centric metatranscriptomics. Microbiome. BioMed Central; 2018;6:194.
Lü F, Bize A, Guillot A, Monnet V, Madigou C, Chapleur O, et al. Metaproteomics of cellulose methanisation under thermophilic conditions reveals a surprisingly high proteolytic activity. ISME J. Nature Publishing Group. 2014;8:88–102.
Kindt A, Liebisch G, Clavel T, Haller D, Hörmannsperger G, Yoon H, et al. The gut microbiota promotes hepatic fatty acid desaturation and elongation in mice. Nat Commun. Nature Publishing Group; 2018;9:3760.
Singer E, Wagner M, Woyke T. Capturing the genetic makeup of the active microbiome in situ. ISME J. Nature Publishing Group; 2017;11:1949.
Angelidaki I, Petersen SP, Ahring BK. Effects of Lipids on Thermophilic Anaerobic-Digestion and Reduction of Lipid Inhibition Upon Addition of Bentonite. Appl Microbiol Biotechnol. 1990;33:469–72.
Kougias PG, Boe K, Einarsdottir ES, Angelidaki I. Counteracting foaming caused by lipids or proteins in biogas reactors using rapeseed oil or oleic acid as antifoaming agents. Water Res. Elsevier. 2015;79:119–27.
APHA. Standard Methods for the Examination of Water and Wastewater: Stand. Methods. American Public Health Association; 2005.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;btu170.
Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. Cold Spring Harbor Lab. 2017;27:824–34.
Kang DD, Froula J, Egan R, Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ. PeerJ Inc.; 2015;3:e1165.
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. Cold Spring Harbor Lab. 2015;25:1043–55.
Bowers RM, Kyrpides NC, Stepanauskas R, Harmon-Smith M, Doud D, Reddy TBK, et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol. 2017;35:725.
O’Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. Oxford University Press. 2015;44:D733–45.
Rinke C, Schwientek P, Sczyrba A, Ivanova NN, Anderson IJ, Cheng J-F, et al. Insights into the phylogeny and coding potential of microbial dark matter. Nature. Nature Publishing Group; 2013;499:431.
Konstantinidis KT, Ramette A, Tiedje JM. The bacterial species definition in the genomic era. Philos Trans R Soc B Biol Sci. The Royal Society. 2006;361:1929–40.
Segata N, Börnigen D, Morgan XC, Huttenhower C. PhyloPhlAn is a new method for improved phylogenetic and taxonomic placement of microbes. Nat Commun. NIH Public Access; 2013;4:2304.
Markowitz VM, Mavromatis K, Ivanova NN, Chen I-MA, Chu K, Kyrpides NC. IMG ER: a system for microbial genome annotation expert review and curation. Bioinformatics. Oxford University Press. 2009;25:2271–8.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. Nature Research. 2012;9:357–9.
Campanaro S, Treu L, Kougias PG, Francisci D, Valle G, Angelidaki I. Metagenomic analysis and functional characterization of the biogas microbiome using high throughput shotgun sequencing and a novel binning strategy. Biotechnol Biofuels. BioMed Central; 2016;9:1.
Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. Oxford University Press. 2015;31:166–9.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. Nature Publishing Group; 2008;5:621.
Eren AM, Esen ÖC, Quince C, Vineis JH, Morrison HG, Sogin ML, et al. Anvi’o: an advanced analysis and visualization platform for ‘omics data. PeerJ. PeerJ Inc.; 2015;3:e1319.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. Oxford University Press. 2010;26:139–40.
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 landfillb. Microbiology. Microbiology Society. 2002;148:3521–30.
Pennings JLA, Keltjens JT, Vogels GD. Isolation and characterization of Methanobacterium thermoautotrophicum ΔH mutants unable to grow under hydrogen-deprived conditions. J Bacteriol. Am Soc Microbiol. 1998;180:2676–81.
Wasserfallen A, Nölling J, Pfister P, Reeve J, De Macario EC. Phylogenetic analysis of 18 thermophilic Methanobacterium isolates supports the proposals to create a new genus, Methanothermobacter gen. nov., and to reclassify several isolates in three species, Methanothermobacter thermautotrophicus comb. nov., Methano. Int J Syst Evol Microbiol. Microbiology Society. 2000;50:43–53.
Zinder SH, Sowers KR, Ferry JG. Methanosarcina thermophila sp. nov., a thermophilic, acetotrophic, methane-producing bacterium. Int J Syst Evol Microbiol. Microbiology Society. 1985;35:522–3.
RIVARD CJ, SMITH PH. Isolation and characterization of a thermophilic marine methanogenic bacterium, Methanogenium thermophilicum sp. nov. Int J Syst Evol Microbiol. Microbiology Society. 1982;32:430–6.
Maus I, Wibberg D, Stantscheff R, Cibis K, Eikmeyer F-G, König H, et al. Complete genome sequence of the hydrogenotrophic archaeon Methanobacterium sp. Mb1 isolated from a production-scale biogas plant. J Biotechnol. Elsevier. 2013;168:734–6.
Denef VJ, Fujimoto M, Berry MA, Schmidt ML. Seasonal succession leads to habitat-dependent differentiation in ribosomal RNA: DNA ratios among freshwater lake bacteria. Front Microbiol. Frontiers. 2016;7:606.
Prosser JI. Dispersing misconceptions and identifying opportunities for the use of’omics’ in soil microbial ecology. Nat Rev Microbiol. Nature Publishing Group; 2015;13:439.
Dortch Q, Roberts TL, Clayton JR Jr, Ahmed SI. RNA/DNA ratios and DNA concentrations as indicators of growth rate and biomass in planktonic marine organisms. Mar Ecol Prog Ser Oldend. 1983;13:61–71.
Porat I, Kim W, Hendrickson EL, Xia Q, Zhang Y, Wang T, et al. Disruption of the operon encoding Ehb hydrogenase limits anabolic CO2 assimilation in the archaeon Methanococcus maripaludis. J Bacteriol. Am Soc Microbiol. 2006;188:1373–80.
Müller B, Sun L, Westerholm M, Schnürer A. Bacterial community composition and fhs profiles of low-and high-ammonia biogas digesters reveal novel syntrophic acetate-oxidising bacteria. Biotechnol Biofuels. BioMed Central; 2016;9:48.
Mosbæk F, Kjeldal H, Mulat DG, Albertsen M, Ward AJ, Feilberg A, et al. Identification of syntrophic acetate-oxidizing bacteria in anaerobic digesters by combined protein-based stable isotope probing and metagenomics. ISME J. Nature Publishing Group. 2016;10:2405–18.
Nobu MK, Narihiro T, Rinke C, Kamagata Y, Tringe SG, Woyke T, et al. Microbial dark matter ecogenomics reveals complex synergistic networks in a methanogenic bioreactor. ISME J. Nature Publishing Group. 2015;9:1710–22.
Akashi H, Gojobori T. Metabolic efficiency and amino acid composition in the proteomes of Escherichia coli and Bacillus subtilis. Proc Natl Acad Sci. National Acad Sciences. 2002;99:3695–700.
Widdel F. Growth of methanogenic bacteria in pure culture with 2-propanol and other alcohols as hydrogen donors. Appl Environ Microbiol. Am Soc Microbiol. 1986;51:1056–62.
Thauer RK, Kaster A-K, Seedorf H, Buckel W, Hedderich R. Methanogenic archaea: ecologically relevant differences in energy conservation. Nat Rev Microbiol. Nature Publishing Group. 2008;6:579–91.
Hori T, Haruta S, Ueno Y, Ishii M, Igarashi Y. Dynamic transition of a methanogenic population in response to the concentration of volatile fatty acids in a thermophilic anaerobic digester. Appl Environ Microbiol. 2006;72:1623–30.
Treu L, Kougias PG, de Diego-Díaz B, Campanaro S, Bassani I, Fernández-Rodríguez J, et al. Two-year microbial adaptation during hydrogen-mediated biogas upgrading process in a serial reactor configuration. Bioresour Technol: Elsevier; 2018.
Treu L, Campanaro S, Kougias PG, Sartori C, Bassani I, Angelidaki I. Hydrogen-fueled microbial pathways in biogas upgrading systems revealed by genome-centric metagenomics. Front Microbiol. Frontiers Media SA. 2018;9.
Westerholm M, Roos S, Schnürer A. Syntrophaceticus schinkii gen. nov., sp. nov., an anaerobic, syntrophic acetate-oxidizing bacterium isolated from a mesophilic anaerobic filter. FEMS Microbiol Lett. Blackwell Publishing Ltd Oxford, UK; 2010;309:100–4.
Manzoor S, Schnürer A, Bongcam-Rudloff E, Müller B. Genome-Guided Analysis of Clostridium ultunense and Comparative Genomics Reveal Different Strategies for Acetate Oxidation and Energy Conservation in Syntrophic Acetate-Oxidising Bacteria. Genes (Basel). Multidisciplinary Digital Publishing Institute; 2018;9:225.
Zhuang W-Q, Yi S, Bill M, Brisson VL, Feng X, Men Y, et al. Incomplete Wood–Ljungdahl pathway facilitates one-carbon metabolism in organohalide-respiring Dehalococcoides mccartyi. Proc Natl Acad Sci. National Acad Sciences. 2014;111:6419–24.
Kersters I, Maestrojuan GM, Torck U, Vancanneyt M, Kersters K, Verstraete W. Isolation of Coprothermobacter proteolyticus from an anaerobic digest and further characterization of the species. Syst Appl Microbiol. Elsevier. 1994;17:289–95.
Lawson CE, Wu S, Bhattacharjee AS, Hamilton JJ, McMahon KD, Goel R, et al. Metabolic network analysis reveals microbial community interactions in anammox granules. Nat Commun. Nature Publishing Group; 2017;8:15416.
Mee MT, Collins JJ, Church GM, Wang HH. Syntrophic exchange in synthetic microbial communities. Proc Natl Acad Sci. National Acad Sciences. 2014;111:E2149–56.
We thank Hector Garcia and Hector Diaz for technical assistance. Sequencing was performed at the Ramaciotti Centre for Genomics (Sydney, Australia).
This work was supported by the Innovation Fund under the project “SYMBIO—Integration of biomass and wind power for biogas enhancement and upgrading via hydrogen assisted anaerobic digestion,” contract 12-132654.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhu, X., Campanaro, S., Treu, L. et al. Metabolic dependencies govern microbial syntrophies during methanogenesis in an anaerobic digestion ecosystem. Microbiome 8, 22 (2020). https://doi.org/10.1186/s40168-019-0780-9
- Anaerobic digestion
- Microbial community
- Syntrophic acetate oxidation
- Glycine cleavage
- Methanogenic pathways