Use of the ADM1 to investigate the effects of acetoclastic methanogen population dynamics on mesophilic digester stability

Proposed protocol for rate-limiting step determination during anaerobic digestion of complex substrates

Anaerobic digestion is a complex process, involving various microorganism groups and, consequently, several reactions. An easy-to-use protocol for the rate-limiting step determination of the process is proposed. The hydrogen production, acetate production, and acetate consumption rates can be calculated, according to a structured algorithm. During the rate limiting step determination, several compounds (biopolymer and monomer representatives, as well as sodium acetate) were used, combined or not with the substrate, to draw the corresponding conclusions. Three substrates were tested, characterized by specific organic compound groups (carbohydrates, proteins, and fats). All three substrates followed the acetate-consuming pathway for the organic matter conversion to methane. In this study, the rate-limiting step for the pathway of acetate consumption was acetate production. Determining the rate-limiting step through the proposed protocol can point to the appropriate actions needed to boost methane production, like substrate pretreatment, using an acidogenic reactor, or checking for the presence of inhibitors.

Integrating Genome-Resolved Metagenomics with Trait-Based Process Modeling to Determine Biokinetics of Distinct Nitrifying Communities within Activated Sludge

Conventional bioprocess models for wastewater treatment are based on aggregated bulk biomass concentrations and do not incorporate microbial physiological diversity. Such a broad aggregation of microbial functional groups can fail to predict ecosystem dynamics when high levels of physiological diversity exist within trophic guilds. For instance, functional diversity among nitrite-oxidizing bacteria (NOB) can obfuscate engineering strategies for their out-selection in activated sludge (AS), which is desirable to promote energy-efficient nitrogen removal. Here, we hypothesized that different NOB populations within AS can have different physiological traits that drive process performance, which we tested by estimating biokinetic growth parameters using a combination of highly replicated respirometry, genome-resolved metagenomics, and process modeling. A lab-scale AS reactor subjected to a selective pressure for over 90 days experienced resilience of NOB activity. We recovered three coexisting Nitrospira population genomes belonging to two sublineages, which exhibited distinct growth strategies and underwent a compositional shift following the selective pressure. A trait-based process model calibrated at the NOB genus level better predicted nitrite accumulation than a conventional process model calibrated at the NOB guild level. This work demonstrates that trait-based modeling can be leveraged to improve our prediction, control, and design of functionally diverse microbiomes driving key environmental biotechnologies.

Determination of Microbial Maintenance in Acetogenesis and Methanogenesis by Experimental and Modeling Techniques

For biogas-producing continuous stirred tank reactors, an increase in dilution rate increases the methane production rate as long as substrate input can be converted fully. However, higher dilution rates necessitate higher specific microbial growth rates, which are assumed to have a strong impact on the apparent microbial biomass yield due to cellular maintenance. To test this, we operated two reactors at 37°C in parallel at dilution rates of 0.18 and 0.07 days-1 (hydraulic retention times of 5.5 and 14 days, doubling times of 3.9 and 9.9 days in steady state) with identical inoculum and a mixture of volatile fatty acids as sole carbon sources. We evaluated the performance of the Anaerobic Digestion Model No. 1 (ADM1), a thermodynamic black box approach (TBA), and dynamic flux balance analysis (dFBA), to describe the experimental observations. All models overestimated the impact of dilution rate on the apparent microbial biomass yield when using default parameter values. Based on our analysis, a maintenance coefficient value below 0.2 kJ per carbon mole of microbial biomass per hour should be used for the TBA, corresponding to 0.12 mmol ATP per gram dry weight per hour for dFBA, which strongly deviates from the value of 9.8 kJ Cmol h-1 that has been suggested to apply to all anaerobic microorganisms at 37°C. We hypothesized that a decrease in dilution rate might select taxa with minimized maintenance expenditure. However, no major differences in the dominating taxa between the reactors were observed based on amplicon sequencing of 16S rRNA genes and terminal restriction fragment length polymorphism analysis of mcrA genes. Surprisingly, Methanosaeta dominated over Methanosarcina even at a dilution rate of 0.18 days-1, which contradicts previous model expectations. Furthermore, only 23-49% of the bacterial reads could be assigned to known syntrophic fatty acid oxidizers, indicating that unknown members of this functional group remain to be discovered. In conclusion, microbial maintenance was found to be much lower for acetogenesis and methanogenesis than previously assumed, likely due to the exceptionally low growth rates in anaerobic digestion. This finding might also be relevant for other microbial systems operating at similarly low growth rates.

Ammonia Inhibition of Anaerobic Volatile Fatty Acid Degrading Microbial Communities

Ammonia inhibition is an important reason for reactor failures and economic losses in anaerobic digestion. Its impact on acetic acid degradation is well-studied, while its effect on propionic and butyric acid degradation has received little attention and is consequently not considered in the Anaerobic Digestion Model No. 1 (ADM1). To compare ammonia inhibition of the degradation of these three volatile fatty acids (VFAs), we fed a mixture of them as sole carbon source to three continuous stirred tank reactors (CSTRs) and increased ammonium bicarbonate concentrations in the influent from 52 to 277 mM. The use of this synthetic substrate allowed for the determination of degradation efficiencies for the individual acids. While butyric acid degradation was hardly affected by the increase of ammonia concentration, propionic acid degradation turned out to be even more inhibited than acetic acid degradation with degradation efficiencies dropping to 31 and 65% for propionic and acetic acid, respectively. The inhibited reactors acclimatized and approximated pre-disturbance degradation efficiencies toward the end of the experiment, which was accompanied by strong microbial community shifts, as observed by amplicon sequencing of 16S rRNA genes and terminal restriction fragment length polymorphism (T-RFLP) of mcrA genes. The acetoclastic methanogen Methanosaeta was completely replaced by Methanosarcina. The propionic acid degrading genus Syntrophobacter was replaced by yet unknown propionic acid degraders. The butyric acid degrading genus Syntrophomonas and hydrogenotrophic Methanomicrobiaceae were hardly affected. We hypothesized that the ammonia sensitivity of the initially dominating taxa Methanosaeta and Syntrophobacter led to a stronger inhibition of the acetic and propionic acid degradation compared to butyric acid degradation and hydrogenotrophic methanogenesis, which were facilitated by the ammonia tolerant taxa Syntrophomonas and Methanomicrobiaceae. We implemented this hypothesis into a multi-taxa extension of ADM1, which was able to simulate the dynamics of both microbial community composition and VFA concentration in the experiment. It is thus plausible that the effect of ammonia on VFA degradation strongly depends on the ammonia sensitivity of the dominating taxa, for syntrophic propionate degraders as much as for acetoclastic methanogens.

Intermittent fasting for microbes: how discontinuous feeding increases functional stability in anaerobic digestion

BACKGROUND

Demand-driven biogas production could play an important role for future sustainable energy supply. However, feeding a biogas reactor according to energy demand may lead to organic overloading and, thus, to process failures. To minimize this risk, digesters need to be actively steered towards containing more robust microbial communities. This study focuses on acetogenesis and methanogenesis as crucial process steps for avoiding acidification. We fed lab-scale anaerobic digesters with volatile fatty acids under various feeding regimes and disturbances. The resulting microbial communities were analyzed on DNA and RNA level by terminal restriction fragment length polymorphism of the mcrA gene, 16S rRNA gene amplicon sequencing, and a [2-13C]-acetate assay. A modified Anaerobic Digestion Model 1 (ADM1) that distinguishes between the acetoclastic methanogens Methanosaeta and Methanosarcina was developed and fitted using experimental abiotic and biotic process parameters.

RESULTS

Discontinuous feeding led to more functional resilience than continuous feeding, without loss in process efficiency. This was attributed to a different microbial community composition. Methanosaeta dominated the continuously fed reactors, while its competitor Methanosarcina was washed out. With discontinuous feeding, however, the fluctuating acetic acid concentrations provided niches to grow and co-exist for both organisms as shown by transcription analysis of the mcrA gene. Our model confirmed the higher functional resilience due to the higher abundance of Methanosarcina based on its higher substrate uptake rate and higher resistance to low pH values. Finally, we applied our model to maize silage as a more complex and practically relevant substrate and showed that our model is likely transferable to the complete AD process.

CONCLUSIONS

The composition of the microbial community determined the AD functional resilience against organic overloading in our experiments. In particular, communities with higher share of Methanosarcina showed higher process stability. The share of these microorganisms can be purposefully increased by discontinuous feeding. A model was developed that enables derivation of the necessary feeding regime for a more robust community with higher share of Methanosarcina.

Population dynamics at digester overload conditions

Two different case studies concerning potential overload situations of anaerobic digesters were investigated and mathematically modelled by means of the Anaerobic Digestion Model No. 1 (ADM1). The first scenario included a digester failure at a municipal WWTP which occurred during revision works of the upstream digester within a two-step digestion system when the sludge was directly by-passed to the 2nd-step reactor. Secondly, the non-occurrence of a highly expected upset situation in a lab-scale digester fed with cattle manure was investigated. ADM1 was utilized to derive indicators which were used to investigate the relationship between digester stability and biomass population dynamics. Conventional design parameters such as the organic loading rate appeared unsuitable for process description under dynamic conditions. Indicators reflecting the biokinetic state (e.g. F(net)/M(net) or the VFA/alkalinity ratio) are more adequate for the assessment of the stability of reactors in transient situations.

Monofermentation of grass silage under mesophilic conditions: measurements and mathematical modeling with ADM 1

In this paper experimental data from grass fermentation and simulation results with the Anaerobic Digestion Model (ADM) No. 1 are described. Two laboratory reactors were operated under mesophilic conditions with volumetric loading rates in between 0.3 and 2.5 kg(VS)/(m(3) x d). Two different kinds of grass silage were used as substrates, resulting in an average specific biogas production of 600 L/kg(VS). The ADM 1 was calibrated both manually and with the help of a Genetic Algorithm in Matlab/Simulink. Results from calibration indicate that the NH3 inhibition constant used to model the inhibition of acetate uptake is three to five times higher compared with digested activated sludge. The hydrogen inhibition constants applied for propionate and valerate/butyrate uptake are around two orders of magnitude lower than for sludge digestion.