- 1Department of Environmental Microbiology, UFZ–Helmholtz Center for Environmental Research, Leipzig, Germany
- 2Biochemical Conversion Department, DBFZ-Deutsches Biomasseforschungszentrum gGmbH, Leipzig, Germany
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.
Introduction
Biogas production is an important renewable energy source and organic waste treatment technology (Plugge, 2017). For nitrogen-rich organic waste, the accumulation of ammonia can become a major problem. Ammonia inhibition has been held responsible for heavy economic losses and even reactor failures (Rajagopal et al., 2013). Suggested solutions to ammonia inhibition are based on the direct removal of ammonia from the reactor, the prevention of high ammonia concentrations by dilution or co-digestion with nitrogen-poor substrates (e.g., maize silage), or by adaptation of the microbial community (Krakat et al., 2017). Bioaugmentation (Fotidis et al., 2014) or support media (Poirier et al., 2017) have been suggested to speed up this adaptation.
Several, partly contradicting theories have been presented on ammonia inhibition in anaerobic digestion. A major controversy is whether acetoclastic or hydrogenotrophic methanogens are more strongly inhibited, with experimental evidence for both cases (Krakat et al., 2017). Furthermore, the shift toward more syntrophic acetate oxidation (SAO) instead of acetoclastic methanogenesis at elevated ammonia concentrations has received much attention (Schnürer and Nordberg, 2008; Werner et al., 2014; Luo et al., 2016). Commonly, free ammonia is thought responsible for ammonia inhibition because it can diffuse into the cells (Rajagopal et al., 2013), but also the ammonia ion is thought to cause inhibition (Astals et al., 2018). The concentration of free ammonia depends on total ammonia nitrogen (TAN) concentration, pH, and temperature. Several underlying and partly connected mechanisms of free ammonia inhibition after diffusion into a cell have been put forward and summarized by Krakat et al. (2017): proton imbalance, change of the intracellular pH, increase in maintenance energy requirement, and inhibition of specific enzymatic reactions. Interestingly, these mechanisms have been solely discussed as explanations for ammonia inhibition in the context of methanogenesis. However, they might also apply to other functional groups, for example proton-reducing bacteria degrading propionic or butyric acid.
Propionic acid and in fewer cases and less in strength also butyric acid accumulations have been observed repeatedly in the context of ammonia inhibition (Li et al., 2017b; Yirong et al., 2017; Peng et al., 2018; Yang et al., 2018). Nevertheless, the inhibition of syntrophic propionic and butyric acid oxidizing bacteria is often neglected in mechanistic descriptions of ammonia inhibition. For example, in the Anaerobic Digestion Model No. 1 (ADM1), only the inhibition of acetoclastic methanogens is included (Batstone et al., 2002). In a more recent adaptation of ADM1, ammonia inhibition was implemented for syntrophic acetic acid oxidizing bacteria (Wett et al., 2014) but still not for syntrophic propionic and butyric acid oxidizers. Quantitative descriptions of the inhibition of volatile fatty acid (VFA) degradation are difficult in complex systems because these acids are simultaneously produced and consumed. While there are several studies on ammonia inhibition using acetic acid as sole carbon source (Steinhaus et al., 2007; Hao et al., 2015; Westerholm et al., 2017), there is only one study on the impact of ammonia inhibition on an open, methanogenic culture fermenting propionic acid as sole carbon source (Li et al., 2017a). In their study, methanogenesis from propionic acid was strongly inhibited, but it could not be concluded if propionic acid oxidizing bacteria are directly inhibited by ammonia or indirectly inhibited by accumulation of hydrogen via an inhibition of the hydrogenotrophic methanogens as suggested earlier (Wiegant and Zeeman, 1986). There are no studies on ammonia inhibition using butyric acid as sole carbon source.
The goal of our study was to compare ammonia inhibition of acetic acid degradation with that of propionic and butyric acid degradation. Therefore, we used a synthetic mixture of acetic, propionic, and butyric acid as substrate amended with micronutrients in three continuous stirred tank reactors (CSTRs), and followed VFA degradation efficiencies and microbial community dynamics over time. Deliberate ammonia inhibition was induced by increasing ammonium bicarbonate concentration in the substrate in two reactors while the third reactor remained unchanged as a control. Furthermore, we added HCl to one of the reactors to reduce the pH and thus the share of free ammonia, which has been shown previously to successfully alleviate ammonia inhibition (Strik et al., 2006). Microbial community compositions were analyzed over the course of the experiment using 16S rRNA gene amplicon sequencing for bacteria and terminal restriction fragment length polymorphism (T-RFLP) profiling of mcrA genes for methanogenic archaea. At the end of the experiment, the functional resilience of the microbial communities was studied by deliberate disturbances in the form of pulse feedings. Finally, ADM1 was amended by second populations for acetic and propionic acid degradation, respectively, to model the effect of microbial community changes on the VFA concentration.
Materials and Methods
Laboratory-Scale CSTR Experiments
Three CSTRs (Rctrl, RNH3, RNH3, HCl) with working volumes of 6 L were operated in parallel at 37°C at a hydraulic retention time (HRT) of 5.5 d. All CSTRs were fed continuously with a synthetic, liquid substrate containing a mixture of VFAs as the only carbon source (45% acetic, 10% propionic, and 45% butyric acid based on chemical oxygen demand, COD) with a total concentration of 37.2 gCOD L−1 in a mineral medium containing all necessary trace elements, macronutrients, and vitamins (see Supplementary Material A, Table S1). The CSTRs were inoculated from a lab-scale digester operated at a HRT of 5.5 d with the same synthetic substrate. The ammonia concentration in the substrate was 52 mM at the start for all CSTRs, and was increased to 277 mM for RNH3 and RNH3, HCl on day 21 until the end of the experiment (see Figure 1). For RNH3, HCl, 100 mM HCl were added starting on day 38. On day 55, 62, and 69, RNH3 and RNH3, HCl were fed not continuously but with the volume of one daily feeding at once as deliberate disturbances. On the day after each disturbance, continuous feeding was resumed as usual. On day 73, HCl addition to the medium of RNH3HCl was ceased until the end of the experiment on day 79. Biogas composition (CH4, CO2, O2, H2, and H2S content), biogas production rate, pH, VFA concentrations, total solids content (TS), and volatile solids content (VS) were determined as described previously (Mulat et al., 2016). The TS and VS content measurements were modified in a way that a reactor sample was not added directly to a crucible, but only the pellet after centrifuging 100 mL reactor content (10,000 × g, 10 min, 10°C) to increase the amount of biomass per analysis and washing the pellet with phosphate buffered saline (PBS) solution (140 mM NaCl, 10 mM Na2HPO4 × 2 H2O, 2.7 mM KCl, 1.8 mM KH2PO4, pH 7.4) to reduce the potential error of VFAs and ammonia in the reactor sample.
Figure 1. Experimental set-up. Three CSTRs (Rctrl, RNH3, RNH3, HCl) were fed with a synthetic medium containing VFAs, mineral medium and varying concentrations of total ammonia nitrogen (TAN) and HCl. Dashed vertical lines indicate the start of the increase in ammonia concentration in the influent of RNH3 and RNH3, HCl. Solid vertical lines indicate deliberate disturbances of RNH3 and RNH3, HCl.
The degradation efficiencies Ei (t) for each individual VFA i were calculated for each sampling date based on mass balances over 1 day considering the VFA input via the influent over the last day Si,in/HRT·Δt, the VFA output via the effluent over the last day ΔSi, out (t), the amount of VFAs in the reactor at the sampling date minus the amount at the previous day ΔSi(t) and, for acetic acid degradation only, the production of acetic acid from propionic and butyric acid oxidation over the last day ΔSi, p (t). VFA mass degraded during 1 day is then given by Si,in/HRT.Δt + ΔSi,p(t) − ΔSi,out(t) − ΔSi(t). Based on a constant reaction volume, the masses of each component can be calculated based on the respective concentrations in the influent and reactor digestate so that degradation efficiency can be expressed as:
with fac,pro and fac,bu being the stoichiometric conversion factors of propionate oxidation (1 mol acetate per mol propionate) and of butyrate oxidation (2 mol acetate per mol of butyrate), respectively. Between sampling dates, VFA concentrations were linearly interpolated. E values of 1 correspond to a complete degradation in steady state. E values can reach values higher than 1 when accumulated VFAs are degraded (ΔSi(t) < 0).
Microbial Community Analyses
DNA and RNA Extraction
Reactor samples were taken using 5 mL pipetting tips that were cut at the top to prevent the exclusion of aggregates of microorganisms and thereby a biased community analysis. Samples were centrifuged at −7°C, 15,000 × g for 2 min and the supernatant was discarded. The low centrifugation temperature and short duration did not lead to sample freezing and were chosen in order to prevent RNA degradation. Pellets were stored at −80°C prior to RNA extraction and −20°C prior to DNA extraction. NucleoSpin® Soil Kit (MACHEREY-NAGEL GmbH & Co. KG, Germany, buffer SL2, no enhancer) was used to extract DNA from the pellets. The quality and quantity of purified DNA were determined by agarose gel electrophoresis and by NanoDrop ND 1000 spectral photometer (Thermo Fisher Scientific, USA). DNA was stored at −20°C. The ZR Soil/Fecal RNA MicroPrepTM Kit (Zymo Research, USA) was used for RNA extraction. DNA was removed using the DNA-free™ DNA Removal Kit (InvitrogenTM, Thermo Fisher Scientific, USA). The quality and quantity of extracted RNA were determined by Qubit® Fluorometer 3.0 (Life technologies, USA, Oregon, Eugene) with the high sensitivity Qubit® RNA HS Assay Kit. RNA was converted to cDNA using the RevertAid H Minus First Strand cDNA Synthesis Kit (Thermo Scientific™, Thermo Fischer Scientific, USA). The resulting cDNA was stored at −20°C. The concentration of cDNA was quantified with a NanoDrop® ND-1000 spectral photometer (Thermo Fisher Scientific, USA).
Composition of the Methanogenic Community
T-RFLP was used to analyze the methanogenic community composition based on mcrA gene amplicons as described previously (Sträuber et al., 2016). The database of Bühligen et al. (2016) was used for the taxonomic assignment. PCR primers mlas (GGT GGT GTM GGD TTC ACM CAR TA) and mcrA-rev (CGT TCA TBG CGT AGT TVG GRT AGT) and PCR program were applied as described previously (Steinberg and Regan, 2008). BstNI (New England Biolabs) was used as restriction enzyme.
Composition of the Bacterial Community
Amplicon sequencing of 16S rRNA genes was used to analyze the bacterial community compositions. PCR amplification using the primers 341f (CCT ACG GGN GGC WGC AG) and 785r (GAC TAC HVG GGT ATC TAA KCC) (Klindworth et al., 2013) and amplicon sequencing with the MiSeq platform (V3, 2 × 300 bp, Illumina) were performed by LGC Genomics GmbH (Berlin, Germany). De-multiplexing and removal of barcodes (allowing 1 mismatch), adapter and primer sequences (allowing 3 mismatches) were performed by LGC Genomics. Forward and reverse reads were merged using the BBMerge 34.48 software (http://bbmap.sourceforge.net/). Merged reads were processed using the QIIME 1.9.1 Virtual Box release (Caporaso et al., 2010). Quality filtering was applied removing low quality reads (quality threshold lower than 20) and allowing no ambiguous base calls. Chimeric sequences were removed and operational taxonomic unit (OTU) clustering was performed using the usearch tool (Edgar, 2010). Taxonomic assignment was performed using the latest MiDAS taxonomy 2.1 (McIlroy et al., 2015) and the RDP classifier 2.2 with a confidence threshold of 0.8 (Wang et al., 2007). The OTU tables were rarefied to 33,013 sequences per sample. The rarefied OTU tables were filtered for bacterial OTUs because the applied primers only partially amplify archaeal 16S rRNA genes and the resulting methanogenic community composition is therefore potentially biased (Klindworth et al., 2013). Raw de-multiplexed sequence data was deposited at EMBL European Nucleotide Archive (ENA) under accession number PRJEB27940 (http://www.ebi.ac.uk/ena/data/view/ PRJEB27940). Relative abundances of 16S rRNA genes were converted to relative genome abundances using the average 16S rRNA gene copy number per genome taken from the rrnDB database (version 5.2).
ADM1 Model Extension
We modified the original ADM1 model structure (Batstone et al., 2002) in various aspects. Mass balances were closed by amending the ordinary differential equations for inorganic carbon and inorganic nitrogen. Biogas production rate was calculated considering the overpressure in the headspace. If not indicated otherwise, parameter values from the benchmark model by Rosen and Jeppsson (2006) were used.
The acetoclastic methanogens (Xac) were replaced by two competing populations: Methanosarcina (Xac,1, see Equation 2a) and Methanosaeta (Xac,2, see Equation 2b). Accordingly, the differential equation for acetic acid was adapted to Equation 2c. The original ADM1 inhibition functions (Batstone et al., 2002) were used with different parameter values for Xac,1 and Xac,2 concerning the ammonia inhibition function (Equation 2d, Table 1).
with Y being the biomass yield, km the maximum substrate uptake rate, KS the half saturation constant, kdec the decay rate, I the inhibition functions, S the substrate concentration, HRT the hydraulic retention time, ac the index for acetic acid, ac1 the index for Methanosarcina, ac2 the index for Methanosaeta, the acetic acid producing reactions (see Peterson Matrix, Supplementary Material B), Inh3,i the inhibition ammonia inhibition function, Snh3 the ammonia concentration and KI,nh3,i the empirical ammonia inhibition parameter with different values for Methanosarcina (i = ac1) and Methanosaeta (i = ac2).
The propionic acid degraders (Xpro) were replaced by two competing populations: an unknown taxon (Xpro,1) and Syntrophobacter (Xpro,2) following the same approach as described above for the acetoclastic methanogens (Supplementary Material A, Section Implementation of two competing propionic acid oxidizers in ADM1). Ammonia inhibition was introduced for both populations following the same inhibition equation that is used for acetic acid degraders in the original ADM1 structure with different parameter values for Xpro,1 and Xpro,2 (see Table 1). A detailed description of the extended model structure is included in the Supplementary Material (Peterson Matrix in Supplementary Material B). All differential and algebraic equations were implemented and solved in the software environment Matlab R2014b (Mathworks).
Cation and anion concentrations in the influent were calculated based on the composition of our synthetic media to 78.3 and 19.8 mM, respectively. Inorganic carbon concentrations were equal to inorganic nitrogen concentrations in the influent with 52 mM before and 277 mM after day 21 based on the ammonium bicarbonate concentration in the synthetic media. Initial concentrations for all state variables (except Xac,1 and Xpro,1) were approximated based on simulation results for steady state concentrations for the inorganic nitrogen concentration of 52 mM in the influent (ADM1 Simulation in Supplementary Material B).
Only parameter values for acetic acid and propionic acid degraders were changed compared to the benchmark model parameter values by Rosen and Jeppsson (2006). The initial concentrations and half-saturation constants of the ammonia tolerant populations (Xac1(t = 0), Xpro1(t = 0), KS,ac1, KS,pro1) and ammonia inhibition constants of all acetic acid and propionic acid degraders (KI,nh3,ac1, KI,nh3,ac2, KI,nh3,pro1, KI,nh3,pro2) were estimated using the solver fmincon (algorithm: interior point) in Matlab R2014b (Mathworks) minimizing the sum of squared error between experimental and simulation results, first for propionic acid concentration and second for acetic acid concentration.
Results
CSTR Performance
On day 21, after 3 HRTs under low ammonia concentrations of 52 mM, all three CSTRs (see Figure 2) had total VFA concentrations of maximal 1.5 g L−1 and pH values of maximal 7.35. Increasing the ammonia concentration in the influent to 277 mM in RNH3 and RNH3, HCl starting on day 21 was followed by an accumulation of VFAs and decline in methane production rates. The strongest accumulation was observed for acetic acid with up to 8.86 g L−1 (RNH3 on day 42). Propionic acid accumulated up to 1.63 g L−1 (RNH3 on day 49) and butyric acid showed the smallest accumulation with up to 0.52 g L−1 (RNH3 on day 42). On day 55, acetic, propionic, and butyric acid concentrations had declined to 0.85, 1.26, and 0.02 g L−1 for RNH3 as well as 0.69, 0.48, and 0.01 g L−1 for RNH3, HCl, respectively. This decline was accompanied by increases in pH values to 7.7 and 7.53 for RNH3 and RNH3, HCl, respectively. Starting on day 73, after HCl was omitted from the medium of RNH3, HCl, the pH value in RNH3, HCl rose to 7.7 toward the end of the experiment, equal to RNH3, while total VFA concentrations remained at about 1 gCOD L−1. Hydrogen partial pressures in the biogas were hardly affected by the increasing ammonia concentration and remained below 17 ppm in all CSTRs between day 21 and 55.
Figure 2. VFA concentration, methane production rate and pH for the three CSTRs Rctrl, RNH3, and RNH3, HCl. Dashed vertical lines indicate the start of the increase in ammonia concentration for RNH3 and RNH3, HCl. Solid vertical lines indicate deliberate disturbances of RNH3 and RNH3, HCl. Error bars represent one standard error of the mean.
On days 55, 62, and 69, RNH3 and RNH3, HCl were given the amount of one daily feeding within 20 min instead of continuous feeding. Both reactors recovered within 3 days to pre-disturbance acetic and butyric acid degradation efficiencies except for RNH3, HCl on day 55 when an acetic acid accumulation was observed followed by a recovery over about one HRT. In both reactors, propionic acid accumulation was observed after a disturbance followed by a recovery over about two HRTs.
During the start of the experiment, the degradation efficiencies of acetic and propionic acid in RNH3 and RNH3, HCl showed declines but recovered again until day 21. On day 21, the degradation efficiencies of all VFAs were higher than 86% in all reactors (Figure 3). After the increase in ammonia concentration on day 21, the lowest degradation efficiencies were 31% for propionic and 65% for acetic acid. Butyric acid degradation was almost not affected showing efficiencies of higher than 94% throughout the experiment.
Figure 3. VFA degradation efficiencies for the three reactors. Dashed lines mark the beginning of the ammonia inhibition and solid lines the deliberate disturbances.
After increasing the ammonia concentration, the formation of aggregates was observed, which reached a size of up to one mm in diameter toward the end of the experiment (see Supplementary Material A, Figure S1).
Methanogenic Community Dynamics
On day 21, all CSTRs showed a similar methanogenic community composition on DNA level, which was dominated by the acetoclastic Methanosaeta (47–56%) and the hydrogenotrophic Methanomicrobiaceae (25–30%). Methanosarcina, which can perform both acetoclastic and hydrogenotrophic methanogenesis, reached a relative abundance of up to 13% (Figure 4).
Figure 4. Methanogenic community analyses by T-RFLP profiling using BstNI as restriction enzyme (A) on DNA basis and (B) on mRNA (cDNA) basis. Note the different time scales in (A) and (B). Days labeled with an asterisk indicate samples taken 5 h after the deliberate disturbances. T-RFs with a length of 340–346 bp can be Methanocorpusculaceae, Methanomicrobiaceae, Methanoregulaceae, and/or Methanospirillacea. T-RFs with a length of 467–471 bp can be Methanobacteriales, Methanococcales, and/or Methanomassiliicoccales.
After increasing the ammonia concentration in the influent to 277 mM in RNH3 and RNH3, HCl starting on day 21, a strong increase in the relative abundance of Methanosarcina with a simultaneous decrease of Methanosaeta was detected for RNH3 and RNH3, HCl, while Methanosaeta continued to dominate over Methanosarcina in Rctrl (Figure 4A). Methanomicrobiaceae remained in all reactors throughout the experiments and reached a relative abundance of about 20% in RNH3 and RNH3, HCl at the end of the experiment.
On cDNA level, the shift from Methanosaeta to Methanosarcina as the dominant methanogen was even more pronounced (Figure 4B). For RNH3, Methanosarcina reached a relative abundance of 90% already on day 36 on cDNA level compared to day 49 on DNA level. Relative T-RF abundances were in general higher on cDNA level compared to DNA level for Methanosarcina in RNH3 and RNH3, HCl and for Methanosaeta in Rctrl. No major impacts of the deliberate disturbances on the relative T-RF abundances of mcrA transcripts were observed.
Bacterial Community Composition
On day 21, 96 to 110 OTUs were detected in the three CSTRs. All CSTRs were dominated by the OTUs Syntrophomonas (8–27% of bacterial reads), Syntrophobacter (13–16%), Thermovirga (15–22%), and Blvii28 wastewater-sludge group (6–38%).
After the increase in ammonia concentration, the community composition of RNH3 and RNH3, HCl showed strong changes. Syntrophobacter, Thermovirga and Blvii28 wastewater-sludge group decreased in relative abundance to < 0.1% of all bacterial reads in both RNH3 and RNH3, HCl toward the end of the experiment (Figure 5). Aminobacterium and Lutispora increased in both reactors from below 1.9% to 4.4–12.6% and from below 0.1% to 1.8–13.9% of all bacterial reads, respectively. The OTU Lineage 1 (Endomicrobia) strongly increased in relative abundance shortly after the ammonia inhibition in RNH3 and decreased again toward the end of the experiment. A similar behavior was observed for Parcubacteria in RNH3, HCl. Syntrophomonas remained one of the most abundant bacterial genera throughout the experiment in all CSTRs.
Figure 5. Bacterial community composition based on 16S rRNA gene amplicon sequencing. Dashed lines mark the beginning of the ammonia inhibition and single solid lines the deliberate disturbances. Only Syntrophomonas, Syntrophobacter, Pelotomaculum, and the seven most abundant OTUs (>12% relative abundance to total bacteria) are shown. Known butyric acid oxidizers are shown in blue, known propionic acid oxidizers in red, and other OTUs in gray.
Our primer pair also led to the amplification of 16S rRNA genes from methanogens despite the lower coverage of archaea compared to bacteria (Klindworth et al., 2013). The ratio of bacterial to total reads varied only slightly throughout the experiment. On day 21, bacterial reads made up 18% and 16% of total reads for RNH3 and RNH3, HCl, respectively, compared to 13 and 17%, respectively, on day 79.
ADM1 Simulation
Simulations using the default ADM1 implementation of ammonia inhibition did not resemble our observed VFA concentrations (Figure 6A, r acid = 0.05, r acid = 0.07). Only acetic acid accumulation was predicted, but not the acclimatization toward the end of the experiment. After implementing our model modifications, i.e., second populations for acetic and propionic acid degradation, respectively, the simulations represented the observed VFA concentrations better (Figure 6B, r acid = 0.98, r acid = 0.56). After the increase of ammonia concentration in the influent on day 21, the ammonia sensitive populations in the simulations (Xac2, Xpro2) became strongly inhibited and were washed out, while accumulations of 8.2 and 2.5 gCOD L−1 occurred for acetic and propionic acid, respectively. At the same time, ammonia tolerant populations (Xac1, Xpro1) increased in relative abundance (Figures 6C,D). On day 55, Xac1 reached a relative abundance to total methanogens of 0.72 and acetic acid accumulation was reduced to 0.5 gCOD L−1. On day 62, Xpro1 reached a relative abundance to total bacteria of 0.12 and propionic acid accumulation was reduced to 0.4 gCOD L−1.
Figure 6. Experimental observations vs. ADM1 simulations: VFA concentrations for (A) the original ADM1 structure and (B) our model including second competing populations for acetic and propionic acid degradation, respectively, as well as ammonia inhibition of propionic acid degradation. (C,D) Relative abundances of acetic acid degraders and propionic acid degraders, respectively. Dashed vertical lines mark the beginning of the ammonia inhibition and solid vertical lines the deliberate disturbances.
Discussion
Apart from a lower pH, the addition of HCl to the influent of RNH3, HCl starting on day 36 led only to minor differences in VFA concentrations, acclimatization time and microbial community composition compared to RNH3. Apparently, at the time of HCl addition, both reactors were already successfully acclimatizing to ammonia inhibition. Furthermore, both reactors were overall little affected by the disturbances. Therefore, both reactors are discussed in the following with the focus on their common response to ammonia inhibition. The VFA accumulations during the start-up of RNH3 and RNH3, HCl seemed to have no decisive effects on the overall experiments since both reactors recovered before the ammonia inhibition was induced.
Strong Ammonia Inhibition of Both Acetic and Propionic Acid Degradation
The inhibition of VFAs other than acetic acid are neglected in ADM1 (Batstone et al., 2002). However, looking at VFA degradation efficiencies (Figure 3), it became clear that propionic acid degradation was even more inhibited than acetic acid degradation. We also observed a stronger or similarly strong inhibition of propionic acid degradation compared to acetic acid degradation in several additional experiments we conducted (see Supplementary Material A). This suggests that ammonia inhibition of propionic acid degradation should receive more attention and be included in anaerobic digestion models.
Calculating VFA degradation efficiencies is only possible in synthetic systems like ours and not for more complex substrates such as manure, because the amounts of individual VFAs produced in acidogenesis from complex substrates are hard to quantify and thus, the VFA degradation efficiencies cannot be calculated. Therefore, the impact of ammonia inhibition on propionic acid degradation might have appeared weaker than it actually was in many studies because lower propionic acid concentrations than acetic acid concentrations were reached.
Concerning the mechanism of ammonia inhibition of propionic acid degradation, an indirect inhibition mechanism was suggested by Wiegant and Zeeman (1986) who argued that a strong ammonia inhibition of hydrogenotrophic methanogenesis can lead to accumulation of hydrogen, which inhibits propionic acid degradation by increasing the Gibbs energy change of catabolism to near or above zero (Figure 7A). However, this indirect route seemed not to play a major role in our experiment since we did not observe any hydrogen accumulation. Furthermore, butyric acid degradation was almost not inhibited in our experiment, which would have been the case if hydrogen accumulated. Therefore, we assumed direct ammonia inhibition of propionic acid degradation as the working hypothesis for our study. Based on our observations, acetic and propionic acid degradation were strongly inhibited in our experiment while butyric acid degradation and hydrogen conversion were hardly affected (Figure 7B).
Figure 7. Inhibition scheme of acetogenesis and methanogenesis following (A) Wiegant and Zeeman (1986) and (B) our working hypothesis. Hydrogen “inhibition” means that increased hydrogen partial pressures increase the Gibbs free energy change of catabolism of propionic and butyric acid oxidation. C2H4O2, C3H6O2, C4H8O2, CO2, H2, CH4, NH3, and NH are the molecular formulas of acetic, propionic, and butyric acid, as well as carbon dioxide, hydrogen, methane, free ammonia and ammonium cation, respectively.
Relationship of VFA Degradation Rates and Microbial Community Composition
After the increase in ammonia concentration started, it took about 5 HRTs until VFA concentrations reached values similar to the start of the experiment in RNH3 and RNH3, HCl. These increases and decreases in VFA concentrations were accompanied by strong shifts in the microbial community.
The propionic acid oxidizing Syntrophobacter was washed out of the reactors after ammonia inhibition started. Also in several other experiments we conducted at elevated ammonia concentrations, Syntrophobacter was washed out (see Supplementary Materials A,C). The absence of Syntrophobacter at high ammonia concentrations was also observed in a mesophilic reactor treating household waste at TAN concentrations of 386–414 mM (Westerholm et al., 2015), supporting our observation that Syntrophobacter is an ammonia sensitive taxon. The only other known syntrophic propionic acid degrading genus in our experiments was Pelotomaculum, which was not competitive with a relative abundance of about 0.1% independent of the ammonia concentration.
Since propionic acid oxidation continued after these two taxa reached negligible relative abundances, one or several yet unknown ammonia tolerant propionic acid oxidizing taxa must be among the detected bacteria. The OTUs that increased in relative abundance after the ammonia inhibition such as Aminobacterium and Lutispora are possible candidates; however, they could grow on butyric acid or decaying biomass as well. Other studies support the ammonia tolerance of these two taxa. Lutispora was found to correlate with the recovery of biogas production in an ammonia inhibited biogas reactor digesting wastewater treatment plant sludge (Chen et al., 2018). Aminobacterium increased in abundance after an increase in ammonia concentration in an anaerobic digester treating chicken manure and feathers with wood shavings (Belostotskiy et al., 2015).
The butyric acid degrading Syntrophomonas tolerated the increase in ammonia concentration in our experiment. A tolerance to TAN of up to 6.5 g L−1 was observed earlier (Lee et al., 2018). Poirier et al. (2017) even observed relative abundances to total bacteria of 22% for ammonia concentrations as high as 25 g L−1, supporting our observation that Syntrophomonas is ammonia tolerant.
The bacterial taxa Thermovirga and Blvii28 wastewater sludge group became abundant with up to 30% relative abundance to total bacteria each during the first 21 days. It is unlikely that they were mainly propionic acid degraders in our reactors because of their high abundance compared to the low share of propionic acid in the feed and it is unlikely that they were mainly acetic acid degraders because at the low free ammonia concentrations during the first 21 days in all reactors (<0.028 g L−1), acetoclastic methanogenesis commonly dominates over SAO (Luo et al., 2016). Therefore, they were most likely syntrophic butyric acid degraders. Since both taxa were washed out after the increase in ammonia concentration in our experiment as well as in an experiment by Lee et al. (2018), they likely are ammonia sensitive.
Concerning acetoclastic methanogens, Methanosaeta dominated all CSTRs at the start of the experiments when ammonia and acetic acid concentrations were low. Its advantage over Methanosarcina at low acetic acid concentrations can be explained by its higher substrate affinity and minimum concentration threshold for acetic acid (Jetten et al., 1992). After the increase in ammonia concentration, Methanosaeta was washed out during the ammonia inhibition and replaced by Methanosarcina. The sensitivity of Methanosaeta to ammonia has been shown in a pure culture study (Steinhaus et al., 2007). Methanosarcina is known to be more tolerant against ammonia inhibition than Methanosaeta (De Vrieze et al., 2012). This tolerance has been connected with their ability to form aggregates (De Vrieze et al., 2012). Aggregates with diameters up to 1 mm were also observed in our experiment. After Methanosarcina completely dominated the acetoclastic methanogens, acetic acid degradation efficiencies recovered up to almost 100% toward the end of the experiment, showing that the applied ammonia concentrations were not inhibitory for Methanosarcina.
The hydrogenotrophic methanogens of the family Methanomicrobiaceae were also tolerant to the increase in ammonia concentration and remained abundant until the end of the experiment. We observed this tolerance in several other experiments on ammonia inhibition that we conducted (see Supplementary Material A). Methanomicrobiaceae have been found previously to be tolerant to TAN concentrations of about 3 g L−1 in pure culture (Schnürer et al., 1999; Nettmann et al., 2010; Wang et al., 2015) and mixed culture studies (Angenent et al., 2002; Westerholm et al., 2011, 2012).
Theoretically, it is possible that all acetic acid was converted via SAO and hydrogenotrophic methanogenesis after day 21. In that case, Methanosarcina would be hydrogenotrophic instead of acetoclastic as assumed above. However, this appears to be unlikely. First, a switch to SAO should have led to an increase in relative abundance of bacteria compared to methanogens that we did not observe. Second, genera known for SAO were observed only in negligible abundances in our reactors. Five species capable of SAO have been cultured belonging to the genera Pseudothermotoga, Thermacetogenium, Clostridium, Syntrophaceticus, and Tepidanaerobacter (Müller et al., 2016). From them, only Clostridium and Syntrophaceticus were detected with maximum relative abundances of 0.2% to total bacteria. Third, anaerobic digesters dominated by SAO have been observed to be dominated by Methanoculleus or Methanobrevibacter instead of Methanosarcina (Luo et al., 2016). Finding acetoclastic Methanosarcina dominant in ammonia inhibited biogas plants is not unusual. Luo et al. (2016) observed the dominance of acetoclastic methanogenesis in industrial biogas plants with up to 0.46 g L−1 of free ammonia which is almost double the maximum free ammonia concentration of 0.25 g L−1 that we observed. Schnürer and Nordberg (2008) observed a dominance of acetoclastic methanogenesis over SAO at free ammonia nitrogen concentrations of up to about 0.3 g L−1 in lab-scale digesters. Therefore, we assumed in our ADM1 simulations that Methanosarcina replaces Methanosaeta as an acetoclastic methanogen. Still, it cannot be excluded that Methanosarcina additionally took up hydrogen and that some SAO occurred in our reactors.
Implementation of Taxon-Specific Ammonia Inhibition in ADM1
The relationship between microbial community dynamics and ammonia as well as VFA concentrations could be successfully reproduced with our extended ADM1 model. However, major changes to the original ADM1 structure were necessary. First of all, inhibition of propionic acid degradation was not implemented in the original ADM1 structure. However, we could show in our experiment that propionic acid degradation was even stronger inhibited than acetic acid degradation, and therefore, the addition of an ammonia inhibition term for propionic acid degradation was necessary. Secondly, simulations based on the original ADM1 structure did not (and in principle cannot) lead to a process recovery after the increase in VFA concentrations following ammonia inhibition because only one taxon of each functional group is part of the original structure. For example, only one taxon is capable of acetoclastic methanogenesis. As a consequence, our simulations using the original ADM1 structure could only lead to no VFA accumulation, to a constant elevated VFA concentration, or process breakdown. However, our experimental data clearly showed an increase in VFA concentration due to ammonia inhibition, followed by a decrease in VFA concentration toward the end of the experiment after the microbial community adapted to the new conditions. Therefore, the addition of at least a second taxon of acetic and propionic acid degraders was necessary in ADM1 to reflect the reactor performance dynamics as a consequence of the microbial community dynamics.
When fitting the simulations to the experimental results, we only changed the ammonia inhibition constants of Xac,1, Xac,2, Xpro,1, and Xpro2, and the initial concentration and half saturation constants of Xac,1 and Xpro,1 compared to the benchmark model parameter values by Rosen and Jeppsson (2006). The fitting resulted in higher substrate affinities (lower half saturation constants) for the ammonia sensitive populations (Xac,2 and Xpro,2) compared to the ammonia tolerant populations (Xac,1 and Xpro,1). This reflects the higher substrate affinities previously observed for Methanosaeta compared to Methanosarcina (Jetten et al., 1992; Straub et al., 2006). However, the advantage of Syntrophobacter over its competitors is not known. Other advantages than substrate affinity are possible, such as higher maximum substrate uptake rates, higher specific growth yields or a combination thereof. A high sensitivity for the parameters KI,nh3,ac, 2, KS,ac,1, KI,nh3,pro,2, KS,ac,2 and the initial concentrations of Xac,1 and Xpro,1 was found (see Sensitivity Analysis in Supplementary Material A). By contrast, increasing the values of the ammonia inhibition constants of the ammonia resistant taxa (KI,nh3,ac,1 and KI,nh3,pro,1), i.e., decreasing the ammonia sensitivity of these populations, did not change the simulation results. This indicates that higher ammonia concentrations are necessary to unambiguously determine these parameter values (see Supplementary Material A, Figure S6). Hydrogenotrophic methanogenesis and butyric acid oxidation were not changed in ADM1 because their performance was hardly impacted by the increase in ammonia concentration in our experiment.
While our model helps to illustrate the relationship between microbial community dynamics and VFA accumulation, extending it for the use as predictive model remains a challenge, in particular regarding parameter identification. In communities with competing taxa, the quantification of taxon-specific substrate uptake rates is a major obstacle. Fitting these rates based on gross consumption likely leads to non-unique solutions and a loss of generality and predictive power. A predictive model would require the physiological characterization of all relevant taxa in pure cultures or defined co-cultures on defined media. Isolating prokaryotes was not always successful in the past (Pelletier et al., 2008). However, the advent of metaomics techniques gives hope that suitable cultivation conditions can be inferred more easily in the future (Overmann et al., 2017). Success in cultivation would be rewarded with high benefits: A predictive multi-taxa ADM1 model could be an essential resource in gaining deterministic control over the microbial community composition, for example by changing process parameters and/or bioaugmentation, in order to increase the productivity of biogas plants suffering from ammonia and other inhibitions.
In conclusion, ammonia inhibition is a major challenge for the biogas process, resulting in economic losses and even process failures. Both our experimental and simulation results showed the importance of ammonia-sensitive taxa, such as Methanosaeta and Syntrophobacter, and ammonia-tolerant taxa, such as Methanomicrobiaceae and Syntrophomonas, for understanding the reactor performance as a result of microbial community dynamics of anaerobic digesters impacted by high ammonia concentrations.
Data Availability Statement
The datasets generated for this study can be found in the Supplementary Material B and the EMBL European Nucleotide Archive (ENA) under accession number PRJEB27940 (http://www.ebi.ac.uk/ena/data/view/ PRJEB27940).
Author Contributions
FB, HS, DP, SK, and FC designed the experiments. FB performed the experiments and simulations and analyzed the reactor and T-RFLP data. SW and FB developed the ADM1 model extension. DP performed the 16S rRNA gene sequencing data analysis. HH contributed to the discussion of the results. All authors contributed to the preparation of the manuscript. All authors read and approved the manuscript.
Funding
This work was funded by the German Federal Ministry of Education and Research (e:Bio project McBiogas, FKZ 031A317).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank Sophia Klink, Ute Lohse, Eric Witt, and Amit Nain for their technical assistance in the T-RFLP analyses, as well as Bärbel Haase, Martin Apelt, Peter Fischer, and Susann Hoffmann for their technical support in analytics of abiotic parameters.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.02921/full#supplementary-material
References
Angenent, L. T., Sung, S., and Raskin, L. (2002). Methanogenic population dynamics during startup of a full-scale anaerobic sequencing batch reactor treating swine waste. Water Res. 36, 4648–4654. doi: 10.1016/S0043-1354(02)00199-9
Astals, S., Peces, M., Batstone, D. J., Jensen, P. D., and Tait, S. (2018). Characterising and modelling free ammonia and ammonium inhibition in anaerobic systems. Water Res. 143, 127–135. doi: 10.1016/j.watres.2018.06.021
Batstone, D. J., Keller, J., Angelidaki, I., Kalyuzhnyi, S. V., Pavlostathis, S. G., Rozzi, A., et al. (2002). Anaerobic Digestion Model No.1. (ADM1). London: IWA Publishing.
Belostotskiy, D. E., Ziganshina, E. E., Siniagina, M., Boulygina, E. A., Miluykov, V. A., and Ziganshin, A. M. (2015). Impact of the substrate loading regime and phosphoric acid supplementation on performance of biogas reactors and microbial community dynamics during anaerobic digestion of chicken wastes. Bioresour. Technol. 193, 42–52. doi: 10.1016/j.biortech.2015.06.066
Bühligen, F., Lucas, R., Nikolausz, M., and Kleinsteuber, S. (2016). A T-RFLP database for the rapid profiling of methanogenic communities in anaerobic digesters. Anaerobe 39, 114–116. doi: 10.1016/j.anaerobe.2016.03.013
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Chen, S., He, J., Wang, H., Dong, B., Li, N., and Dai, X. (2018). Microbial responses and metabolic pathways reveal the recovery mechanism of an anaerobic digestion system subjected to progressive inhibition by ammonia. Chem. Eng. J. 350, 312–323. doi: 10.1016/j.cej.2018.05.168
De Vrieze, J., Hennebel, T., Boon, N., and Verstraete, W. (2012). Methanosarcina: The rediscovered methanogen for heavy duty biomethanation. Bioresour. Technol. 112, 1–9. doi: 10.1016/j.biortech.2012.02.079
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Fotidis, I., Wang, H., Fiedel, N. R., Luo, G., Karakashev, D. B., and Angelidaki, I. (2014). Bioaugmentation as a solution to increase methane production from an ammonia-rich substrate. Environ. Sci. Technol. 48, 7669–7676. doi: 10.1021/es5017075
Hao, L., Lü, F., Mazéas, L., Desmond-Le Quéméner, E., Madigou, C., Guenne, A., et al. (2015). Stable isotope probing of acetate fed anaerobic batch incubations shows a partial resistance of acetoclastic methanogenesis catalyzed by Methanosarcina to sudden increase of ammonia level. Water Res. 69, 90–99. doi: 10.1016/j.watres.2014.11.010
Jetten, M. S. M., Stams, A. J. M., and Zehnder, A. J. B. (1992). Methanogenesis from acetate: a comparison of the acetate metabolism in Methanothrix soehngenii and Methanosarcina spp. FEMS Microbiol. Lett. 88, 181–198. doi: 10.1111/j.1574-6968.1992.tb04987.x
Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast, C., Horn, M., et al. (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41, e1–e1. doi: 10.1093/nar/gks808
Krakat, N., Anjum, R., Dietz, D., and Demirel, B. (2017). Methods of ammonia removal in anaerobic digestion: a review. Water Sci. Technol. 76, 1925–1938. doi: 10.2166/wst.2017.406
Lee, J., Kim, E., Han, G., Tongco, J. V., Shin, S. G., and Hwang, S. (2018). Microbial communities underpinning mesophilic anaerobic digesters treating food wastewater or sewage sludge: a full-scale study. Bioresour. Technol. 259, 388–397. doi: 10.1016/j.biortech.2018.03.052
Li, Y., Zhang, Y., Kong, X., Li, L., Yuan, Z., Dong, R., et al. (2017a). Effects of ammonia on propionate degradation and microbial community in digesters using propionate as a sole carbon source. J. Chem. Technol. Biotechnol. 92, 2538–2545. doi: 10.1002/jctb.5260
Li, Y., Zhang, Y., Sun, Y., Wu, S., Kong, X., Yuan, Z., et al. (2017b). The performance efficiency of bioaugmentation to prevent anaerobic digestion failure from ammonia and propionate inhibition. Bioresour. Technol. 231, 94–100. doi: 10.1016/j.biortech.2017.01.068
Luo, G., Fotidis, I., and Angelidaki, I. (2016). 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. doi: 10.1186/s13068-016-0465-6
McIlroy, S. J., Saunders, A. M., Albertsen, M., Nierychlo, M., McIlroy, B., Hansen, A. A., et al. (2015). MiDAS: the field guide to the microbes of activated sludge. Database 2015:bav062. doi: 10.1093/database/bav062
Mulat, D. G., Jacobi, H. F., Feilberg, A., Adamsen, A. P. S., Richnow, H.-H., and Nikolausz, M. (2016). Changing feeding regimes to demonstrate flexible biogas production: effects on process performance, microbial community structure, and methanogenesis pathways. Appl. Environ. Microbiol. 82, 438–449. doi: 10.1128/AEM.02320-15
Müller, B., Sun, L., Westerholm, M., and Schnürer, A. (2016). Bacterial community composition and fhs profiles of low- and high-ammonia biogas digesters reveal novel syntrophic acetate-oxidising bacteria. Biotechnol. Biofuels 9:48. doi: 10.1186/s13068-016-0454-9
Nettmann, E., Bergmann, I., Pramschufer, S., Mundt, K., Plogsties, V., Herrmann, C., et al. (2010). Polyphasic analyses of methanogenic archaeal communities in agricultural biogas plants. Appl. Environ. Microbiol. 76, 2540–2548. doi: 10.1128/AEM.01423-09
Overmann, J., Abt, B., and Sikorski, J. (2017). Present and future of culturing bacteria. Annu. Rev. Microbiol. 71, 711–730. doi: 10.1146/annurev-micro-090816-093449
Pelletier, E., Kreimeyer, A., Bocs, S., Rouy, Z., Gyapay, G., Chouari, R., et al. (2008). “Candidatus Cloacamonas Acidaminovorans”: genome sequence reconstruction provides a first glimpse of a new bacterial division. J. Bacteriol. 190, 2572–2579. doi: 10.1128/JB.01248-07
Peng, X., Zhang, S., Li, L., Zhao, X., Ma, Y., and Shi, D. (2018). Long-term high-solids anaerobic digestion of food waste: effects of ammonia on process performance and microbial community. Bioresour. Technol. 262, 148–158. doi: 10.1016/j.biortech.2018.04.076
Poirier, S., Madigou, C., Bouchez, T., and Chapleur, O. (2017). Improving anaerobic digestion with support media: Mitigation of ammonia inhibition and effect on microbial communities. Bioresour. Technol. 235, 229–239. doi: 10.1016/j.biortech.2017.03.099
Rajagopal, R., Massé, D. I., and Singh, G. (2013). A critical review on inhibition of anaerobic digestion process by excess ammonia. Bioresour. Technol. 143, 632–641. doi: 10.1016/j.biortech.2013.06.030
Rosen, C., and Jeppsson, U. (2006). Aspects on ADM1 Implementation within the BSM2 Framework. Lund. Available online at: http://www.iea.lth.se/publications/Reports/LTH-IEA-7224.pdf
Schnürer, A., and Nordberg, A. (2008). Ammonia, a selective agent for methane production by syntrophic acetate oxidation at mesophilic temperature. Water Sci. Technol. 57, 735–740. doi: 10.2166/wst.2008.097
Schnürer, A., Zellner, G., and Svensson, B. H. (1999). Mesophilic syntrophic acetate oxidation during methane formation in biogas reactors. FEMS Microbiol. Ecol. 29, 249–261. doi: 10.1016/S0168-6496(99)00016-1
Steinberg, L. M., and Regan, J. M. (2008). Phylogenetic comparison of the methanogenic communities from an acidic, oligotrophic fen and an anaerobic digester treating municipal wastewater sludge. Appl. Environ. Microbiol. 74, 6663–6671. doi: 10.1128/AEM.00553-08
Steinhaus, B., Garcia, M. L., Shen, A. Q., and Angenent, L. T. (2007). A portable anaerobic microbioreactor reveals optimum growth conditions for the methanogen Methanosaeta concilii. Appl. Environ. Microbiol. 73, 1653–1658. doi: 10.1128/AEM.01827-06
Straub, A. J., Conklin, A. S. Q., Ferguson, J. F., and Stensel, H. D. (2006). Use of the ADM1 to investigate the effects of acetoclastic methanogen population dynamics on mesophilic digester stability. Water Sci. Technol. 54, 59–66. doi: 10.2166/wst.2006.526
Sträuber, H., Lucas, R., and Kleinsteuber, S. (2016). Metabolic and microbial community dynamics during the anaerobic digestion of maize silage in a two-phase process. Appl. Microbiol. Biotechnol. 100, 479–491. doi: 10.1007/s00253-015-6996-0
Strik, D. P. B. T. B., Domnanovich, A. M., and Holubar, P. (2006). A pH-based control of ammonia in biogas during anaerobic digestion of artificial pig manure and maize silage. Process Biochem. 41, 1235–1238. doi: 10.1016/j.procbio.2005.12.008
Wang, H., Fotidis, I. a, and Angelidaki, I. (2015). Ammonia effect on hydrogenotrophic methanogens and syntrophic acetate oxidizing bacteria. FEMS Microbiol. Ecol. 91:fiv130. doi: 10.1093/femsec/fiv130
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Werner, J. J., Garcia, M. L., Perkins, S. D., Yarasheski, K. E., Smith, S. R., Muegge, B. D., et al. (2014). Microbial community dynamics and stability during an ammonia-induced shift to syntrophic acetate oxidation. Appl. Environ. Microbiol. 80, 3375–3383. doi: 10.1128/AEM.00166-14
Westerholm, M., Dolfing, J., Sherry, A., Gray, N. D., Head, I. M., and Schnürer, A. (2011). Quantification of syntrophic acetate-oxidizing microbial communities in biogas processes. Environ. Microbiol. Rep. 3, 500–505. doi: 10.1111/j.1758-2229.2011.00249.x
Westerholm, M., Levén, L., and Schnürer, A. (2012). Bioaugmentation of syntrophic acetate-oxidizing culture in biogas reactors exposed to increasing levels of ammonia. Appl. Environ. Microbiol. 78, 7619–7625. doi: 10.1128/AEM.01637-12
Westerholm, M., Müller, B., Isaksson, S., and Schnürer, A. (2015). Trace element and temperature effects on microbial communities and links to biogas digester performance at high ammonia levels. Biotechnol. Biofuels 8, 1–19. doi: 10.1186/s13068-015-0328-6
Westerholm, M., Müller, B., Singh, A., Karlsson Lindsjö, O., and Schnürer, A. (2017). Detection of novel syntrophic acetate-oxidizing bacteria from biogas processes by continuous acetate enrichment approaches. Microb. Biotechnol. doi: 10.1111/1751-7915.13035
Wett, B., Takács, I., Batstone, D., Wilson, C., and Murthy, S. (2014). Anaerobic model for high-solids or high-temperature digestion - Additional pathway of acetate oxidation. Water Sci. Technol. 69, 1634–1640. doi: 10.2166/wst.2014.047
Wiegant, W. M., and Zeeman, G. (1986). The mechanism of ammonia inhibition in the thermophilic digestion of livestock wastes. Agric. Wastes 16, 243–253. doi: 10.1016/0141-4607(86)90056-9
Yang, Z., Wang, W., He, Y., Zhang, R., and Liu, G. (2018). Effect of ammonia on methane production, methanogenesis pathway, microbial community and reactor performance under mesophilic and thermophilic conditions. Renew. Energy 125, 915–925. doi: 10.1016/j.renene.2018.03.032
Keywords: anaerobic digestion, biogas, methanogens, syntrophic propionate oxidation, 16S rRNA amplicon sequencing, T-RFLP, ADM1, chemostat
Citation: Bonk F, Popp D, Weinrich S, Sträuber H, Kleinsteuber S, Harms H and Centler F (2018) Ammonia Inhibition of Anaerobic Volatile Fatty Acid Degrading Microbial Communities. Front. Microbiol. 9:2921. doi: 10.3389/fmicb.2018.02921
Received: 27 July 2018; Accepted: 14 November 2018;
Published: 30 November 2018.
Edited by:
Pankaj Kumar Arora, Babasaheb Bhimrao Ambedkar University, IndiaReviewed by:
Stefan Junne, Technische Universität Berlin, GermanyIn-Hyun Nam, Korea Institute of Geoscience and Mineral Resources, South Korea
Copyright © 2018 Bonk, Popp, Weinrich, Sträuber, Kleinsteuber, Harms and Centler. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Florian Centler, Zmxvcmlhbi5jZW50bGVyQHVmei5kZQ==