Metagenomic Analysis of Anaerobic Microbial Communities Degrading Short-Chain Fatty Acids as Sole Carbon Sources

Analyzing microbial communities using metagenomes is a powerful approach to understand compositional structures and functional connections in anaerobic digestion (AD) microbiomes. Whereas short-read sequencing approaches based on the Illumina platform result in highly fragmented metagenomes, long-read sequencing leads to more contiguous assemblies. To evaluate the performance of a hybrid approach of these two sequencing approaches we compared the metagenome-assembled genomes (MAGs) resulting from five AD microbiome samples. The samples were taken from reactors fed with short-chain fatty acids at different feeding regimes (continuous and discontinuous) and organic loading rates (OLR). Methanothrix showed a high relative abundance at all feeding regimes but was strongly reduced in abundance at higher OLR, when Methanosarcina took over. The bacterial community composition differed strongly between reactors of different feeding regimes and OLRs. However, the functional potential was similar regardless of feeding regime and OLR. The hybrid sequencing approach using Nanopore long-reads and Illumina MiSeq reads improved assembly statistics, including an increase of the N50 value (on average from 32 to 1740 kbp) and an increased length of the longest contig (on average from 94 to 1898 kbp). The hybrid approach also led to a higher share of high-quality MAGs and generated five potentially circular genomes while none were generated using MiSeq-based contigs only. Finally, 27 hybrid MAGs were reconstructed of which 18 represent potentially new species-15 of them bacterial species. During pathway analysis, selected MAGs revealed similar gene patterns of butyrate degradation and might represent new butyrate-degrading bacteria. The demonstrated advantages of adding long reads to metagenomic analyses make the hybrid approach the preferable option when dealing with complex microbiomes.

Proteiniphilum and Methanothrix harundinacea became dominant acetate utilizers in a methanogenic reactor operated under strong ammonia stress

Microorganisms in anaerobic digestion (AD) are easily affected by ammonia, especially acetoclastic methanogens. Thus, in ammonia-suppressed AD systems, acetate degradation is reported to be carried out mainly by the cooperation of syntrophic acetate oxidizers and hydrogenotrophic methanogens. Previous studies have revealed ammonia inhibition on microbial flora by AD performance, but the effect mechanism of ammonia on microbial metabolism remains poorly understood. In this study, we constructed a mesophilic chemostat fed with acetate as the sole carbon source, gradually increased the total ammonia nitrogen (TAN) concentration from 1 g L^-1 to 6 g L^-1, and employed the 16S rRNA gene, metagenomics, and metatranscriptomics analysis to characterize the microbial community structure and metabolic behavior. The results showed that even at the TAN of 6 g L^-1 (pH 7.5), the methanogenesis kept normal, the biogas production was approximately 92% of that at TAN of 1 g L^-1 and the acetate degradation ratio reached 99%, suggesting the strong TAN tolerance of the microbial community enriched. 16S rRNA gene analysis suggested that the microbial community structure changed along with the TAN concentration. Methanothrix predominated in methanogens all the time, in which the dominant species was gradually replaced from M. soehngenii to M. harundinacea with the increased TAN. Dominant bacterial species also changed and Proteiniphilum showed a significant positive correlation with increased TAN. Meta-omics analysis showed that the absolute dominant microorganisms at TAN of 6 g L^-1 were M. harundinacea and Proteiniphilum, both of which highly expressed genes for anti-oxidative stress. M. harundinacea and the second dominant methanogen Methanosarcina highly expressed both acetate cleavage and CO₂ reduction pathways, suggesting the possibility that these two pathways contributed to methanogenesis together. Proteiniphilum and some other species in Firmicutes and Synergistetes were likely acetate oxidizers in the community as they highly expressed genes for syntrophic acetate oxidization, H₂ generation, and electron transfer. These results suggested that Proteiniphilum as well as M. harundinacea have strong ammonia tolerance and played critical roles in acetate degradation under ammonia-suppressed conditions. The achievements of the study would contribute to the regulation and management of the AD process.

Methanogen Community Dynamics and Methanogenic Function Response to Solid Waste Decomposition

Methane production during solid waste decomposition is a typical methanogen-mediated and enzyme-catalyzed anaerobic digestion (AD). Methanogen community dynamics and metabolic diversity during the decomposition are not known. In this study, we investigated methanogen community dynamics and methanogenic pathways during solid waste decomposition in a bioreactor using high-throughput Illumina MiSeq sequencing and phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt), respectively. We also related the methanogen community differences with solid waste and leachate physiochemical parameters. Results showed that the percentage of biodegradable matter (BDM) in solid waste decreased from 55 ± 5% in aerobic phase (AP) to 30 ± 2% in anaerobic acid phase (ACP), and to 13 ± 11% in methanogenic phase (MP), resulting in 76% BDM consumption by microbes. Methanogen community structure varied in AP, ACP, and MP, showing that Methanomicrobiales and Methanosarcinales were dominant in AP and MP and archaea E2 was abundant in ACP. Each phase had abundant core methanogen orders, genera, and species with significant difference (p < 0.05). Redundancy analysis showed that biochemical oxygen demand (BOD₅) and nitrate significantly influenced methanogen community composition, suggesting that methanogen community structure is nutrient-dependent. Two methanogenic pathways including acetoclastic and hydrogenotrophic pathways with associated functional genes differed at three phases. ACP had the lowest abundance of these genes, indicating that methanogenesis was inhibited in acidogenesis. Abundant hydrogenotrophic and acetoclastic methanogenesis functional genes in MP and AP are in response to the abundance of Methanomicrobiales and Methanosarcinales. The findings provide previously unidentified insight into the mechanism of methanogen community structure and function during solid waste bioconversion for methane.

Metabolic Adaptation of Methanogens in Anaerobic Digesters Upon Trace Element Limitation

Anaerobic digestion (AD) is a complex multi-stage process relying on the activity of highly diverse microbial communities including hydrolytic, acidogenic and syntrophic acetogenic bacteria as well as methanogenic archaea. The lower diversity of methanogenic archaea compared to the bacterial groups involved in AD and the corresponding lack of functional redundancy cause a stronger susceptibility of methanogenesis to unfavorable process conditions such as trace element (TE) deprivation, thus controlling the stability of the overall process. Here, we investigated the effects of a slowly increasing TE deficit on the methanogenic community function in a semi-continuous biogas process. The aim of the study was to understand how methanogens in digester communities cope with TE limitation and sustain their growth and metabolic activity. Two lab-scale biogas reactors fed with distillers grains and supplemented with TEs were operated in parallel for 76 weeks before one of the reactors was subjected to TE deprivation, leading to a decline of cobalt and molybdenum concentrations from 0.9 to 0.2 mg/L, nickel concentrations from 2.9 to 0.8 mg/L, manganese concentrations from 38 to 18 mg/L, and tungsten concentrations from 1.4 to 0.2 mg/L. Amplicon sequencing of mcrA genes revealed Methanosarcina (72%) and Methanoculleus (23%) as the predominant methanogens in the undisturbed reactors. With increasing TE limitation, the relative abundance of Methanosarcina dropped to 67% and a slight decrease of acetoclastic methanogenic activity was observed in batch tests with ¹³C-methyl-labeled acetate, suggesting a shift toward syntrophic acetate oxidation coupled to hydrogenotrophic methanogenesis. Metaproteome analysis revealed abundance shifts of the enzymes involved in methanogenic pathways. Proteins involved in methylotrophic and acetoclastic methanogenesis decreased in abundance while formylmethanofuran dehydrogenase from Methanosarcinaceae increased, confirming our hypothesis of a shift from acetoclastic to hydrogenotrophic methanogenesis by Methanosarcina. Both Methanosarcina and Methanoculleus increased the abundance of N5-methyltetrahydromethanopterin-coenzyme M methyltransferase and methyl-coenzyme M reductase. However, these efforts to preserve the ion motive force for energy conservation were seemingly more successful in Methanoculleus. We conclude that both methanogenic genera use different strategies to stabilize their energy balance under TE limitation. Methanosarcina switched from TE expensive pathways (methylotrophic and acetoclastic methanogenesis) to hydrogenotrophic methanogenesis. Methanoculleus showed a higher robustness and was favored over the more fastidious Methanosarcina, thus stabilizing reactor performance under TE limitation.

Metabolic Adaptation of Methanogens in Anaerobic Digesters Upon Trace Element Limitation

Anaerobic digestion (AD) is a complex multi-stage process relying on the activity of highly diverse microbial communities including hydrolytic, acidogenic and syntrophic acetogenic bacteria as well as methanogenic archaea. The lower diversity of methanogenic archaea compared to the bacterial groups involved in AD and the corresponding lack of functional redundancy cause a stronger susceptibility of methanogenesis to unfavorable process conditions such as trace element (TE) deprivation, thus controlling the stability of the overall process. Here, we investigated the effects of a slowly increasing TE deficit on the methanogenic community function in a semi-continuous biogas process. The aim of the study was to understand how methanogens in digester communities cope with TE limitation and sustain their growth and metabolic activity. Two lab-scale biogas reactors fed with distillers grains and supplemented with TEs were operated in parallel for 76 weeks before one of the reactors was subjected to TE deprivation, leading to a decline of cobalt and molybdenum concentrations from 0.9 to 0.2 mg/L, nickel concentrations from 2.9 to 0.8 mg/L, manganese concentrations from 38 to 18 mg/L, and tungsten concentrations from 1.4 to 0.2 mg/L. Amplicon sequencing of mcrA genes revealed Methanosarcina (72%) and Methanoculleus (23%) as the predominant methanogens in the undisturbed reactors. With increasing TE limitation, the relative abundance of Methanosarcina dropped to 67% and a slight decrease of acetoclastic methanogenic activity was observed in batch tests with ¹³C-methyl-labeled acetate, suggesting a shift toward syntrophic acetate oxidation coupled to hydrogenotrophic methanogenesis. Metaproteome analysis revealed abundance shifts of the enzymes involved in methanogenic pathways. Proteins involved in methylotrophic and acetoclastic methanogenesis decreased in abundance while formylmethanofuran dehydrogenase from Methanosarcinaceae increased, confirming our hypothesis of a shift from acetoclastic to hydrogenotrophic methanogenesis by Methanosarcina. Both Methanosarcina and Methanoculleus increased the abundance of N5-methyltetrahydromethanopterin-coenzyme M methyltransferase and methyl-coenzyme M reductase. However, these efforts to preserve the ion motive force for energy conservation were seemingly more successful in Methanoculleus. We conclude that both methanogenic genera use different strategies to stabilize their energy balance under TE limitation. Methanosarcina switched from TE expensive pathways (methylotrophic and acetoclastic methanogenesis) to hydrogenotrophic methanogenesis. Methanoculleus showed a higher robustness and was favored over the more fastidious Methanosarcina, thus stabilizing reactor performance under TE limitation.

Improved Monitoring of Semi-Continuous Anaerobic Digestion of Sugarcane Waste: Effects of Increasing Organic Loading Rate on Methanogenic Community Dynamics

The anaerobic digestion of filter cake and its co-digestion with bagasse, and the effect of gradual increase of the organic loading rate (OLR) from start-up to overload were investigated. Understanding the influence of environmental and technical parameters on the development of particular methanogenic pathway in the biogas process was an important aim for the prediction and prevention of process failure. The rapid accumulation of volatile organic acids at high OLR of 3.0 to 4.0 gvs·L⁻¹·day⁻¹ indicated strong process inhibition. Methanogenic community dynamics of the reactors was monitored by stable isotope composition of biogas and molecular biological analysis. A potential shift toward the aceticlastic methanogenesis was observed along with the OLR increase under stable reactor operating conditions. Reactor overloading and process failure were indicated by the tendency to return to a predominance of hydrogenotrophic methanogenesis with rising abundances of the orders Methanobacteriales and Methanomicrobiales and drop of the genus Methanosarcina abundance.