Next Article in Journal
Pilot Scale Evaluation of Wild Saccharomyces cerevisiae Strains in Aglianico
Next Article in Special Issue
Biomethane Potential in Anaerobic Biodegradation of Commercial Bioplastic Materials
Previous Article in Journal
A Review Unveiling Various Machine Learning Algorithms Adopted for Biohydrogen Productions from Microalgae
Previous Article in Special Issue
Anaerobic Co-Digestion of Agricultural Residues Produced in Southern and Northern Greece
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Water Potential on Anaerobic Methane Production and a Microbial Consortium

1
Department of Plant Life & Environmental Science, Hankyong National University, Anseong 17579, Gyeonggi, Republic of Korea
2
Biogas Research Center, Hankyong National University, Anseong 17579, Gyeonggi, Republic of Korea
3
Plant Engineering Center, Institute for Advanced Engineering, Yongin 17528, Gyeonggi, Republic of Korea
*
Author to whom correspondence should be addressed.
Fermentation 2023, 9(3), 244; https://doi.org/10.3390/fermentation9030244
Submission received: 10 February 2023 / Revised: 28 February 2023 / Accepted: 28 February 2023 / Published: 3 March 2023

Abstract

:
This study probed the effect of the water potential (Ψ) on anaerobic methane production and a microbial consortium. The Ψ level of the investigated anaerobic digester (n = 20) was in the range from −0.10 to −2.09 MPa with a mean value of −1.23 MPa, and the Ψ level of the anaerobic digester was significantly correlated with the SCOD, TKN, NH4+-N, alkalinity, salinity (SPS), NH4+, Na+, K+, Cl, NO3, and PO43− (p < 0.001). The maximum methane production rate (Rm) of the Control (−0.40 MPa) was 8.11 mL day−1 and decreased to 1.70 mL day−1 at −3.91 MPa (K5), and the lag growth phase time (λ) was delayed to 35.96 and 25.34 days at −2.85 MPa (K4) and −3.91 MPa (K5), respectively. The ultimate methane potential (Bu) was 0.264 Nm3 kg−1-VSadded for the Control, and when Ψ was adjusted, Bu increased to 0.278 Nm3 kg−1-VSadded at −1.49 MPa (K3) but decreased to 0.203 and 0.172 Nm3 kg−1-VSadded at −2.85 MPa (K4) and −3.91 MPa (K5), respectively. Therefore, the methane yield was inhibited due to the decrease in Ψ, and the methane yield is predicted to be inhibited from about −1.65 MPa. In the genus-level taxonomic classification of the microbial community, the relative abundance of Methanosarcina decreased significantly to 36.76% at −3.91 MPa (K5) compared to 58.15% for the Control; however, the relative abundance of Methanoculleus significantly increased to 35.16% at −3.91 MPa (K5) compared to 14.85% for the Control.

1. Introduction

Anaerobic digestion technology is a biological conversion technology to produce methane (CH4) from organic matter using a sequential microbial reaction. The operating factors of a typical anaerobic digestion process are those factors known to affect the activity of methanogens, such as the pH, alkalinity, temperature, total volatile fatty acids (TVFAs), and hydraulic retention time (HRT) [1,2,3]. In order to maximize the efficiency of methane production of an anaerobic digester, it is important to maintain the conditions such that they optimally support the activity of the methanogens, which is inhibited outside of the optimal range of operating factors of the digester. Kroeker et al. [4] classified the factors that inhibit the activity of methanogenic bacteria in an anaerobic digestion process into two types: toxicity and inhibition. Toxicity is represented by the complete cessation of methanogenic activity, and inhibition occurs as a result of reducing the rate and extent of methanogenesis. In particular, ammonia (NH3), hydrogen sulfide (H2S), cations (Na+, K+, Ca2+, Mg2+), long chain fatty acids (LCFA), heavy metals, etc., have been reported as inhibitors of the activity of methanogens [5].
In general, organic waste from resources such as livestock manure, food waste, and sewage sludge is used in the anaerobic digestion process [6,7]. In Korea, since the methane production efficiency is low in the case of anaerobic digestion of livestock manure and sewage sludge, a co-digestion method, in which these two types of waste are mixed with food waste with relatively high methane production efficiency, is used [8]. However, the food waste generated in Korea contains a high concentration of salt of approximately 3% (based on dry matter content) [9], which causes an operational problem for the anaerobic digestion of food waste owing to salt stress. The inhibitory effect of salt on anaerobic digestion is known to be due to cations such as sodium ion (Na+), potassium ion (K+), calcium ion (Ca2+), and magnesium ion (Mg2+) [10]. In previous studies related to cation inhibition, the inhibition of anaerobic digestion efficiency was reported to occur in the concentration range of Na+ from 4100 to 11,000 mg L−1, K+ from 3000 to 28,000 mg L−1, Ca2+ from 5000 to 8000 mg L−1, and Mg2+ from 729 to 3000 mg L−1 [9,11,12,13,14,15,16,17,18]. As such, the cation inhibitory effect on anaerobic digestion is diverse in that it occurs across a wide concentration range of each ionic species, making it difficult to manage the inhibitory concentration for each type of cation when operating an anaerobic digester. Therefore, commercial anaerobic digestion reactors are operated without taking into consideration the inhibition of the anaerobic digester by the presence of excess cations.
The water potential (Ψ) is an indicator of the state of the potential energy of water and can be used to evaluate the availability of water for microorganisms. In an anaerobic digester, the Ψ is mainly composed of the osmotic potential (Ψo), which represents the colligative property of the solution affected by the concentration of salt, organic acid, etc. In general, high concentrations of salt are known to lower Ψ and induce osmotic stress by cell dehydration, which results in the inhibition of microbial activity and metabolism [19,20,21,22,23]. Generally, the cytoplasm of microorganisms has a relatively low Ψ due to the concentration of various solutes [23,24], but severe osmotic stress induces plasmolysis [25,26] and cytolysis [27] of microorganisms. Excessive salt in the anaerobic digester may therefore severely inhibit the activity and metabolism of anaerobic microorganisms because of the low Ψ. However, studies on the effect of salt on the anaerobic digestion efficiency have revealed a range of inhibitory concentrations for each ionic species, yet the dependence of the anaerobic digestion efficiency on Ψ, which represents a comprehensive indicator of salt stress, has not been reported. This motivated our study, in which we attempted to understand the effect of Ψ on anaerobic methane production by the microbial reaction of methanogens. The Ψ level of 20 anaerobic digesters was evaluated, and the effect of the physicochemical factors of these anaerobic digesters on the Ψ was analyzed. In addition, a batch-type anaerobic reactor was operated by adjusting the Ψ using potassium chloride (KCl). This enabled us to analyze the changes in the characteristics of a community of anaerobic microorganisms according to the changes in Ψ using next-generation sequencing (NGS) technology.

2. Materials and Methods

2.1. Materials

The physicochemical properties and the value of Ψ were determined for anaerobic digestates collected from 20 anaerobic digestion facilities in operation in Korea (Table 1). The facilities we selected for our investigation were those in which livestock manure, sewage sludge, and food waste are anaerobically digested or co-digested, and these facilities were operated under mesophilic conditions.

2.2. Water Potential (Ψ) Analysis

Ψ is defined as the difference in the chemical potential of water (J mol−1) per unit volume (J m−3) between a given water sample and pure free water at the same temperature. This is expressed as Equation (1). Rearrangement of the units in Equation 1 shows that Ψ is expressed in terms of the pressure (MPa), which is frequently used as the unit in Ψ measurements [28].
Ψ = μ W μ ° W V W
where Ψ is the water potential (MPa), μ W is the chemical potential of the water, μ W is the chemical potential of the free water at the same temperature, and V W is the partial molar volume of water in the system. The Ψ was measured using a dew point hygrometer (WP4C, METER Group, Inc., Pullman, MA, USA) based on the principle of the chilled mirror dew point. When the air in the chamber containing the sample flows over the cooling mirror, dew forms on the surface of the mirror. The Ψ of a liquid sample is found by relating the sample Ψ reading to the saturation vapor pressure of air in equilibrium with the sample. The relationship between the sample Ψ and the saturation vapor pressure of the air is determined using Equation (2) [29].
Ψ = R T M w l n e s T d e s T s
where Ψ is the water potential (MPa), e s T d is the saturation vapor pressure of the air at dew point temperature, e s T s is the saturation vapor pressure at sample temperature, R is the gas constant, 8.31 J mol−1 K−1, T is the temperature of the sample in Kelvin, and M w is the molecular mass of water.

2.3. Methane Production Potential

2.3.1. Theoretical Methane Potential (Bth)

The theoretical methane potential (Bth) was calculated stoichiometrically using Boyle’s equation based on the elemental analysis results of the samples (Equations (3) and (4)).
C a H b O c N d S e + a b 4 c 2 + 3 d 4 + e 2 H 2 O a 2 + b 8 c 4 3 d 8 e 4 C H 4 + a 2 b 8 + c 4 + 3 d 8 + e 4 C O 2 + d N H 3 + e H 2 S
B t h N m 3   k g 1 V S a d d e d = 22.4 ×   4 a + b 2 c 3 d 2 e / 8 ) 12 a + b + 16 c + 14 d + 32 e  

2.3.2. Water Potential (Ψ) Inhibition Assay

The ultimate methane potential (Bu) was assessed by the biochemical methane potential assay [30]. To assess the Bu in the Ψ inhibition assay, a batch-type anaerobic reactor was operated under mesophilic conditions (38 °C). The anaerobic inoculum was collected from a farm-scale anaerobic digester located in Icheon City, Republic of Korea. The chemical properties of the inoculum are provided in Table 2.
The batch anaerobic reactors for the Ψ inhibition assay were prepared as in Table 3. The inoculum for the Ψ inhibition assay was kept under mesophilic conditions at 38 °C for two weeks to remove any remaining biodegradable fraction. The substrate to inoculum ratio in all anaerobic batch reactors was equal to 0.5 (g-VSsubstrate/g-VSinoculum). The Ψ of each batch anaerobic reactor was adjusted to −0.57, −0.89, −1.49, −2.85, and −3.91 MPa using the KCl solution for the Ψ inhibition assay. The working volume for anaerobic batch fermentation was 200 mL of a 250 mL serum bottle. The headspace of the serum bottle was filled with N2 gas and sealed with a butyl rubber stopper. The anaerobic batch reactors for each sample and blank were incubated for up to 78 days in the convection incubator and manually mixed each day during the fermentation period. Then, the anaerobic batch reactions for each sample and blank were performed in triplicate. The biochemical methane potential was calculated based on the volatile solid (VS) content. The biochemical methane potentials of the samples were corrected using the blank value and calibrated under standard temperature and pressure (STP) conditions (0 °C, 1 atm). The modified Gompertz model (Equation (5)) was employed to interpret the progress of cumulative methane production. This enabled the cumulative methane production data to be optimized using the following equations [31].
M = P × exp e x p R m · e P λ t + 1
where M is the cumulative methane production (mL), t is the anaerobic fermentation time (days), P is the final methane production (mL), e is the exp (1), R m is the maximum methane production rate (mL day−1), and λ represents the lag growth phase time (days). The cumulative methane production curves in the Ψ inhibition assay were optimized with SigmaPlot (SigmaPlot Version 12.5, Systat Software Inc., Chicago, IL, USA) using the modified Gompertz model. The inhibition degree in the Ψ inhibition assay was used to describe the adverse impact of water stress on the methane yield by Equation (6). Herein, I D is the inhibition degree (%), B u C is the ultimate methane potential (Nm3 kg−1-VSadded) in the Control, and B u T is the ultimate methane potential (Nm3 kg−1-VSadded) after the treatment.
I D % = 1 B u C B u T

2.4. Analysis of the Microbial Consortium

DNA was extracted using a DNA kit (GD141-050 Gram positive, Biofact, Seoul, Republic of Korea) according to the manufacturer’s instructions. The extracted DNA was quantified using Quant-IT PicoGreen (Invitrogen, Waltham, MA, USA).

DNA Extraction and Quantification

The sequencing libraries were prepared according to the Illumina 16S Metagenomic Sequencing Library protocols to amplify the V3 and V4 regions. The input gDNA 2 ng was PCR amplified with 5× reaction buffer, 1 mM of dNTP mix, 500 nM each of the universal F/R PCR primer, and Herculase II fusion DNA polymerase (Agilent Technologies, Santa Clara, CA, USA). The cycling conditions for the 1st PCR were 3 min at 95 °C for heat activation, and 25 cycles of 30 s at 95 °C, 30 s at 55 °C, and 30 s at 72 °C, followed by a 5 min final extension at 72 °C. The universal primer pair with Illumina adapter overhang sequences used for the first amplifications were as follows: V3-F: 5’-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG-3’, V4-R: 5’-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACHVGGTATCTAATCC-3’. The 1st PCR product was purified with AMPure beads (Agencourt Bioscience, Beverly, MA, USA). Following purification, 2 µL of the 1st PCR product was PCR amplified for final library construction containing the index using the NexteraXT Indexed Primer. The cycling conditions for the 2nd PCR were the same as for the 1st PCR except that 10 cycles were used. The PCR product was purified with AMPure beads. The final purified product was then quantified using qPCR according to the qPCR Quantification Protocol Guide (KAPA Library Quantification kits for Illumina Sequencing platforms) and qualified using the TapeStation D1000 ScreenTape (Agilent Technologies, Waldbronn, Germany). The paired-end (2 × 300 bp) sequencing was performed by the Macrogen online sequencing system using the MiSeq™ platform (Illumina, San Diego, CA, USA).

2.5. Chemical Analysis

The physicochemical properties of the collected anaerobic digestates were determined and the contents of the reactors in the Ψ inhibition assay were analyzed based on standard methods [29]. The total solid (TS), volatile solid (VS), fixed solid (FS), total suspended solid (TSS), volatile suspended solid (VSS), fixed suspended solid (FSS), total chemical oxygen demand (TCOD), soluble chemical oxygen demand (SCOD), total Kjeldahl nitrogen (TKN), ammonium nitrogen (NH4+-N), and total volatile fatty acids (TVFAs) were analyzed to determine the organic composition. In addition, the pH, alkalinity, salinity (SPS), cations (Na+, K+, Mg2+, Ca2+), and anions (NO3, Cl, SO42−, PO43−) were analyzed to determine the chemical properties. More specifically, the VFAs were measured using a gas chromatograph (GC2010, Shimadzu Scientific Instruments, Inc., Columbia, MD, USA) equipped with a flame ionization detector with an automatic sampler. The concentrations of the soluble cations (Na+, K+, Mg2+, Ca2+) were determined by analyzing the filtrate from which TSS was removed using inductively coupled plasma optical emission spectrometry (ICP-OES) (Avio 550 Max, PerkinElmer, USA). The concentrations of the soluble anions (NO3, Cl, SO42−, PO43−) in the filtrate from which TSS was removed were analyzed by ion chromatography (IC) (833 Basic IC plus, Metrohm, Switzerland). The salinity was defined as the practical salinity (SPS), which was calculated from the electrical conductivity measured by an electrical conductivity meter (CM-42X, TOADKK, Japan) [32]. In the Ψ inhibition assay, the total gas production was measured daily for the first seven days and then every two or three days. The gas that was produced displaced an acidified brine solution in a burette and the volume of the displaced solution was recorded after correcting for atmospheric pressure [33]. The CH4 and CO2 concentrations in the gas samples were determined using a gas chromatograph (Clarus 680, PerkinElmer, Inc., Waltham, MA, USA) equipped with a thermal conductivity detector and a HayeSepQ packed column (CRS, Inc., Louisville, KY, USA). The column was operated with helium as a carrier gas at a flow rate of 5 mL/min. The temperatures of the injector, oven, and detector were set to 150, 90, and 150 °C, respectively [34]. Total ammonium nitrogen (TAN) has two forms, NH4+-N and NH3-N, which were calculated according to Equation (7) [35], where T is the absolute temperature (K).
F r e e   N H 3 N = T A N × 1 + 10 p H 10 0.09018 + 2729.92 T 1

2.6. Statistical Analysis

In this study, to examine the physicochemical parameters that affect the Ψ of the anaerobic digester, the Pearson correlation coefficients between Ψ and the physicochemical parameters (pH, WC, TS, VS, FS, TSS, VSS, FSS, TCOD, SCOD, TKN, NH4+-N, alkalinity, TVFAs, SPS, Na+, K+, Mg2+, Ca2+, NO3, Cl, SO42−, PO43−) were analyzed by R studio (R Version 4.1.3, R Foundation for Statistical Computing, Vienna, Austria). The tables in this article present the mean values and standard deviations of the data obtained from the experiments. The statistical analysis of the results of this experiment was conducted using the general linear model (GLM) procedure of the SAS program package (SAS ver. 9.4, SAS Instrument Inc., Cary, NC, USA), and the significant difference (p < 0.05) of the mean between treatments was tested using Duncan’s multiple range test.

3. Results and Discussion

3.1. Correlation Factors with Ψ

The physicochemical properties of 20 anaerobic digesters operating in Korea were investigated, and the results are reported in Table 4 and Table 5. The pH of the anaerobic digester had mean, maximum, and minimum values of 7.86, 8.78, and 6.89, respectively. In some facilities, it was necessary to adjust the pH of the digester to lower the pH, but the mean value was within the normal operating range. In addition, the TKN (4.91 g L−1) and alkalinity (16.82 g L−1) were in the optimum operating range of an anaerobic digester. The Ψ of the anaerobic digester had mean, maximum, and minimum values of −1.23 MPa, −0.10 MPa, and −2.09 MPa, respectively.
The value of Ψ of an anaerobic digester is mainly affected by ionic compounds (cations and anions), and as the concentration of ionic compounds increases, the Ψ decreases [36,37]. In addition, organic acids, amino acids, polypeptides, and proteins with acidic dissociation functional groups also cause Ψ to decrease in solution. However, among the ionic compounds, multivalent ionic species such as Ca2+, Mg2+, SO42−, and PO43− often exist as precipitated substances while forming complexes with each other, and these precipitated substances do not affect Ψ [38,39,40,41]. In this study, the NH4+ concentration of the anaerobic digestate was 3694 mg/L, which was high compared to other cationic compounds, and the concentrations of Na+ and K+ were 1789 and 2263 mg L−1, respectively. In this study, the types of ionic compounds that affect Ψ were identified by measuring the ionic compounds present in soluble forms after removing TSS. Measurement of the concentration after the exclusion of ionic compounds that are present in the anaerobic digester but are inactive due to adsorption and precipitation revealed the concentration of Cl to be the highest among those of the ionic compounds in the anaerobic digester. The physicochemical factors affecting the Ψ of the anaerobic digester were evaluated by calculating the Pearson correlation coefficient, and Ψ was correlated with the SCOD, TKN, NH4+-N, and alkalinity at a significance level of 99.9% (p < 0.001). The correlation between Ψ and the solid content was significant at TS (p < 0.05) and VS (p < 0.01). Ψ was correlated with SPS, NH4+, Na+, K+, Cl, NO3, and PO43− at a significance level of 99.9% (p < 0.001), and with Mg2+ and SO42− at a significance level of 99.0% (p < 0.01).
The correlation between the Ψ of the anaerobic digester and the Ca2+ content was evaluated as very low. As for the correlation between Ψ and ionic compounds, the divalent cations were reported to affect Ψ [10,37], but in this study, the effect of polyvalent cations on Ψ was evaluated as low. In this study, by evaluating the correlation between Ψ and the soluble forms of the ionic compounds in the anaerobic digestate, the small effect on Ψ of divalent cations such as Ca2+ is considered to be the consequence of the relatively high degree of adsorption and precipitation of these cations. The SCOD is an indicator of soluble organic substances such as sugars and organic acids present in the anaerobic digester [42] and was significantly correlated with Ψ of the anaerobic digester. TKN is an indicator of the organic (proteins and amino acids) and inorganic nitrogen (NH4+-N, NH3-N) content of the anaerobic digester, and has a significant correlation with Ψ [43]. In particular, NH4+-N is a weak acid (pKa = 9.245), and at pH lower than pH 9.245, the NH4+ ionic form is dominant, and above pH 9.245, the NH3 gas form is dominant [44]. Therefore, TAN and NH4+-N of the anaerobic digester were strongly correlated (p < 0.001) with Ψ. The alkalinity is an indicator of the amount of H+ buffer substances such as HCO3, CO32−, and OH, and had a strong correlation (p < 0.001) with Ψ along with the pH (p < 0.01) [45].

3.2. Water Potential (Ψ) Inhibition Assay

In order to understand the effect of Ψ on methane production in anaerobic digestion, the Ψ of the anaerobic batch reactors was adjusted to −0.40 (Control), −0.57 (K1), −0.89 (K2), −1.49 (K3), −2.85 (K4), and −3.91 MPa (K5) and they were operated for 78 days. During this period, the production and consumption of TVFAs in the anaerobic reactor were delayed as Ψ decreased (Table 6). The anaerobic digestion of organic matter proceeds via microbial chemical reactions in the stages of hydrolysis, acidogenesis, acetogenesis, and methanogenesis [46], and the results of this study showed that low Ψ levels simultaneously had inhibitory effects on the methanogens as well as the microorganisms involved in acetogenesis. However, during the operation of the batch reactors, the Ψ of the reactors was maintained as stable without significant fluctuations. These results showed that the effect of the variation in Ψ according to the change in the concentration of the TVFAs was not significant during the operation of the anaerobic digester.
Figure 1 shows the cumulative methane yield curves in the Ψ inhibition assay, and Table 5 lists the parameters and organic matter degradation rate (VSr) derived by optimization of the modified Gompertz model. The value of Bu was 0.264 Nm3 kg−1-VSadded for the Control (−0.40 MPa). When Ψ was adjusted, Bu increased to 0.278 Nm3 kg−1-VSadded at −1.49 MPa (K3), but decreased to 0.203 and 0.172 Nm3 kg−1-VSadded at −2.85 MPa (K4) and −3.91 MPa (K5), respectively. The effect on methane production according to the decrease in Ψ was clearly reflected in the maximum methane production rate (Rm) and lag growth phase time (λ); that is, Rm for the Control (−0.40 MPa) was 8.11 mL day−1 and decreased to 1.70 mL day−1 at K5 (−3.91 MPa), and λ was delayed to 35.96 and 25.34 days at K4 (−2.85 MPa) and K5 (−3.91 MPa), respectively (Table 7). Therefore, the level of Ψ that inhibits methane production in the anaerobic digestion of glucose was expected to be between −1.49 MPa (K+ 2.22%) and −2.85 MPa (K+ 3.07%) (Figure 2). In batch-type anaerobic digestion using cow manure, anaerobic digestion was reportedly inhibited by 50% at a K+ concentration of 2.80% (w/v) [16], and a similar K+ inhibitory concentration range was determined in this study. However, at very low Ψ levels of −2.85 MPa or less, the methane yield and methane production rate rapidly decreased, whereupon the activity of anaerobic microorganisms gradually recovered when sufficient acclimatization time elapsed.
In general, methanogens are known to have survived even in high concentrations of salt [19,47], and it has been reported that, by accumulating osmotic substances, the intracellular osmotic potential is lowered and adapted to the low Ψ level of the surroundings [48,49,50]. In general, in a low Ψ environment induced by cations such as K+, plasmolysis occurs as water flows out of the cells. However, conversion of the environment of these cells to a hypotonic solution would enable the cells to recover to their normal condition [25,26], although the inhibition of Ψ by substances with a high molecular mass such as polyethylene glycol (PEG) results in cytorrhysis because of a strong decrease in intracellular pressure [26,51,52]. Therefore, the fact that gradual methane production proceeded even at the lowest Ψ level (−3.91 MPa) in this study indicates that moisture with a Ψ of −3.91 MPa did not have a fatal effect such as the collapse of the cell wall of anaerobic microorganisms. This was evaluated as the Ψ level at which cell activity could be restored by acclimatization. Figure 3 shows the degree of inhibition of methane production according to the decrease in Ψ and the increase in the K+ concentration, respectively. Analysis of the Ψ inhibition level using linear regression in the Ψ range in which the methane yield was inhibited owing to the decrease in Ψ enabled us to predict that the methane yield was inhibited from about −1.65 MPa. At this time, the K+ concentration at which the methane yield was inhibited was found to be 11.69 g L−1. Several researchers conducted anaerobic digestion toxicity studies by K+ during anaerobic digestion (Table 8) and reported the concentration of K+ capable of inhibiting anaerobic digestion according to the feeding materials and operating temperatures in various reactors. In this study, the concentration of K+ at which anaerobic digestion is inhibited (11.69 g L−1) was within the inhibition range reported by previous researchers. Generally, the inhibition effect of K+ is rarely referenced in the literature. Low concentrations of potassium (less than 400 mg L−1) were reported to enhance the performance of both the thermophilic and mesophilic ranges, whereas there was an inhibitory effect in the thermophilic temperature range at higher K+ concentrations. However, the concentration of K+ at which anaerobic digestion is inhibited varies from 3000 to 28,000 mg L−1 depending on the type of feeding material and the temperature [5,11,14,15,16], and the mechanism of metabolic inhibition of anaerobic microorganisms by K+ has not been reported in detail. This suggests that the inhibitory effect on anaerobic digestion by K+ may be an indirect inhibitory effect due to the decrease in Ψ that accompanies high concentrations of K+ rather than a direct toxic effect of K+.

3.3. Microbial Community

The effect of Ψ on changes in the characteristics of a community of microorganisms in anaerobic digestion was investigated by performing the Ψ inhibition assay in the Ψ range of −0.40 to −3.91 MPa. Then, after the batch reactor had been in operation for 78 days, the characteristics of a community of microorganisms were analyzed using the NGS method. During anaerobic digestion, microbial methane is mainly produced by the methanogenic archaea, which could produce methane primarily from CO2/H2, methyl groups, or acetate under anaerobic conditions [53]. Methanogenic archaea are slow-growing microorganisms, and methane is produced by keeping the syntrophic mechanism with variable microorganisms concerned in the acetogenesis [54]. Methanogenic archaea are assumed to belong to the Euryarchaeota phylum. In the Euryarchaeota phylum, the methanogenic archaea are classified into eight orders, namely Methanobacteriales, Methanococcales, Methanomicrobiales, Methanosarcinales, Methanopyrales, Methanocellales, Methanomassiliicoccales, and Methanonatronarchaeales [54,55]. In this study, the four orders of Methanobacteriales, Methanomicrobiales, Methanomassiliicoccales, and Methanosarcinales constituted the community of methanogens (Figure 4), and the community of methanogens accounted for 73.47 to 80.14% of the total microbial community.
In particular, Methanosarcinales had the largest distribution in the methanogen community, with 58.34 in the Control (−0.40 MPa), 56.56 in K1 (−0.57 MPa), 58.88 in K2 (−0.89 MPa), 52.38 in K3 (−1.49 MPa), 55.80 in K4 (−2.85 MPa), and 37.19% K5 (−3.91 MPa), with a rapid decrease in the cluster distribution at the Ψ level of K5 (−3.91 MPa). Although the distribution of Methanomicrobiales was 14.85% in the Control (−0.40 MPa), the community distribution increased from 14.84% in K1 (−0.57 MPa) to 35.19% in K5 (−3.91 MPa) as Ψ decreased. The proportion of Methanomassiliicoccales in the methanogen community was as low as 1.27% in the Control (−0.40 MPa) and decreased from 1.21% in K1 (−0.57 MPa) to 0.92% in K5 (−3.91 MPa) according to the decrease in Ψ. In addition, the proportion of Methanobacteriales among methanogens was as low as 0.59% in the Control (−0.40 MPa). However, it showed community characteristics that increased from 0.86% in K1 (−0.57 MPa) to 1.33% in K5 (−3.91 MPa) according to the decrease in Ψ. The methanogenesis process consists of four main pathways: CO2 reduction, acetoclastic, methylotrophic, and methyl reduction pathways [55]. The Methanosarcinales order is capable of undergoing methanogenesis via all four methanogenic pathways. Generally, the Methanomicrobiales order can reduce CO2 to methane with electrons derived from the oxidation of H2 via the CO2 reduction pathway [56]. Contrary to this, the species in the Methanobacteriales and Methanomassiliicoccales orders proceed along a methyl reduction pathway by reducing methyl compounds to methane with electrons derived from H2 [55,57]. In the genus-level taxonomic classification of the microbial community, Methanosarcina and Methanothrix were found to be methanogens belonging to the Methanosarcinales order, and the community proportion of Methanosarcina, reported as an acetoclastic methanogen [58], decreased significantly to 36.76% in K5 (−3.91 MPa) in comparison to the 58.15% in the Control (−0.40 MPa) (Figure 5). However, in the Methanomicrobiales order, Methanoculleus, which has been reported as a hydrogenotrophic methanogen [58], formed the majority of the colony, and it significantly increased to 35.16% in K5 (−3.91 MPa) compared to 14.85% in the Control (−0.40 MPa).
Therefore, in this study, at a low Ψ level of −3.91 MPa, the population of acetoclastic methanogens (Methanosarcina) decreased and the population of hydrogenotrophic methanogens (Methanoculleus) increased, and these results suggested that hydrogenotrophic methanogens are advantageous for acclimation under low Ψ conditions. In general, hydrogenotrophic methanogens are syntrophic partners that are associated with microorganisms in the acetogenesis. This led to the speculation that changes in the community characteristics of the methanogen according to the decrease in Ψ are related to the decrease in the activity of microorganisms at the stages of acidogenesis and acetogenesis. In this study, the production of TVFAs was increasingly delayed during the operation of the anaerobic reactor as the Ψ level decreased (Table 6), and these results showed that the microorganisms in the acidogenesis and acetogenesis stages were also inhibited by the low Ψ. In particular, syntrophism occurs as a result of indirect interspecies electron transfer (IIET), in which hydrogen or formate, a product of acetogenesis, acts as an electron carrier between methanogens [59]. This means that, even though the adaptation mechanism of hydrogenotrophic methanogens to low Ψ is important to increase the community of hydrogenotrophic methanogens at a low Ψ, the adaptation mechanism of acetogenic microorganisms related to IIET should also be considered. Microorganisms belonging to the firmiutes phylum are acidogenic and acetogenic and are involved in the production of VFAs and acetate [60,61]. In particular, the Eubacteriales and thermoanaerobacterales orders are involved in acetogenesis, producing acetate, H2, and CO2. In this study, Clostridium belonging to the Eubacteriales order and Syntrophaceticus belonging to the thermoanaerobacterales order maintained or increased their community distribution at a low Ψ (Figure 5). In addition, Clostridium maintained a range of 1.05 to 1.64% in the Ψ range from −0.40 MPa to −3.91 MPa, and Syntrophaceticus had a population of 0.00% in the Control (−0.40 MPa) but it increased to 2.09% in K5 (−3.91 MPa).

4. Conclusions

This study was carried out to understand the effect of Ψ on anaerobic methane production and a microbial consortium. The Ψ of the investigated anaerobic digester (n = 20) ranged from −0.10 to −2.09 MPa with a mean value of −1.23 MPa, and the Ψ of the anaerobic digester was significantly correlated with the SCOD, TKN, NH4+-N, alkalinity, SPS, NH4+, Na+, K+, Cl, NO3, and PO43− (p < 0.001). In the Ψ inhibition assay, the Rm of the Control (−0.40 MPa) was 8.11 mL day−1 and decreased to 1.70 mL day−1 at −3.91 MPa (K5), and the λ was delayed to 35.96 and 25.34 days at −2.85 MPa (K4) and −3.91 MPa (K5), respectively. The Bu was 0.264 Nm3 kg−1-VSadded for the Control, and Bu increased to 0.278 Nm3 kg−1-VSadded at −1.49 MPa (K3) but decreased to 0.203 and 0.172 Nm3 kg−1-VSadded at −2.85 MPa (K4) and −3.91 MPa (K5), respectively. Therefore, the methane yield was inhibited due to the decrease in the Ψ, with inhibition of the methane yield predicted to occur from about −1.65 MPa, which corresponded to the K+ concentration of 11.69 g L−1. In the genus-level taxonomic classification of the microbial community, the relative abundance of Methanosarcina (acetoclastic methanogens) decreased significantly to 36.76% (58.15% for the Control) at the lowest Ψ level; however, the relative abundance of Methanoculleus (hydrogenotrophic methanogens) significantly increased to 35.16% (14.85% for the Control) at this Ψ level. In this study, not only methane-producing bacteria but also acid- and acetic acid-producing microorganisms were expected to be inhibited by low Ψ. Existing studies on salt inhibition in an anaerobic digester mainly focused on different types of cations (Na+, K+, Ca2+, Mg2+) and examined the phenomenon of inhibition of the anaerobic digestion efficiency by individual cations. Ψ is a colligative property that depends on the ratio of the number of solute particles to the number of solvent particles in a solution. In addition, Ψ is not only affected by the salt content but also by various micro- and macromolecules with TVFAs, amino acids, and proteins that contain electronic functional groups present in the anaerobic digester. However, few studies on the effect of Ψ on anaerobic digestion efficiency have been reported thus far, and future research on the inhibition of anaerobic digestion efficiency by cations would require in-depth research in relation to Ψ, which represents the colligative properties of the solution. These studies are anticipated to make it possible to use Ψ as a new operating factor for salt management in the operational management of an anaerobic digester.

Author Contributions

Conceptualization, Y.-M.Y. and J.Y.; methodology, J.Y.; software, J.-H.L. and E.S.; validation, C.-G.K., E.S. and Y.-M.Y.; formal analysis, J.Y.; investigation, J.Y.; resources, Y.-M.Y.; data curation, Y.-M.Y.; writing—original draft preparation, J.Y. and C.-G.K.; writing—review and editing, Y.-M.Y.; visualization, J.-H.L.; supervision, Y.-M.Y.; project administration, Y.-M.Y.; funding acquisition, Y.-M.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by MOLIT and KAIA’s development of combined plant compact technology (Project No. RS-2020-KA157945).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

This study was supported through MOLIT and KAIA’s development of combined plant compact technology (Project No. RS-2020-KA157945).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Björnsson, L.; Murto, M.; Jantsch, T.G.; Mattiasson, B. Evaluation of new methods for the monitoring of alkalinity, dissolved hydrogen and the microbial community in anaerobic digestion. Water Res. 2001, 35, 2833–2840. [Google Scholar] [CrossRef] [PubMed]
  2. Beuvink, J.; Spoelstra, S.; Hogendorp, R. An automated method for measuring time-course of gas production of feedstuffs incubated with buffered rumen fluid. Neth. J. Agric. Sci. 1992, 40, 401–407. [Google Scholar] [CrossRef]
  3. Ahring, B.K.; Sandberg, M.; Angelidaki, I. Volatile fatty acids as indicators of process imbalance in anaerobic digestors. Appl. Microbiol. Biotechnol. 1995, 43, 559–565. [Google Scholar] [CrossRef]
  4. Kroeker, E.; Schulte, D.; Sparling, A.; Lapp, H. Anaerobic treatment process stability. J. (Water Pollut. Control Fed.) 1979, 51, 718–727. [Google Scholar]
  5. Chen, Y.; Cheng, J.J.; Creamer, K.S. Inhibition of anaerobic digestion process: A review. Bioresour. Technol. 2008, 99, 4044–4064. [Google Scholar] [CrossRef]
  6. Oh, S.-E.; Lee, M.-K.; Kim, D.-H. Continuous mesophilic-dry anaerobic digestion of organic solid waste. J. Korean Soc. Environ. Eng. 2009, 31, 341–345. [Google Scholar]
  7. Oh, K.-S.; Hwang, J.-K.; Song, Y.-J.; Kim, M.-J.; Park, J.-G.; Pak, D.-W. Integrated digestion of thermal solubilized sewage sludge to improve anaerobic digestion efficiency of organic waste. J. Korean Soc. Water Environ. 2022, 38, 95–102. [Google Scholar] [CrossRef]
  8. Lee, S.; Yoon, Y.-S.; Kang, J.-G.; Kim, K.-H.; Shin, S.K. Anaerobic co-digestion characteristics of food waste leachate and sewage sludge. J. Korea Org. Resour. Recycl. Assoc. 2016, 24, 21–29. [Google Scholar] [CrossRef] [Green Version]
  9. Kim, D.-H.; Kim, S.-H.; Shin, H.-S. Sodium inhibition of fermentative hydrogen production. Int. J. Hydrogen Energy 2009, 34, 3295–3304. [Google Scholar] [CrossRef]
  10. McCarty, P.L.; McKinney, R.E. Salt toxicity in anaerobic digestion. J. (Water Pollut. Control Fed.) 1961, 33, 399–415. [Google Scholar]
  11. Zha, X.; Tsapekos, P.; Alvarado-Morales, M.; Lu, X.; Angelidaki, I. Potassium inhibition during sludge and biopulp co-digestion; experimental and model-based approaches. Waste Manag. 2020, 113, 304–311. [Google Scholar] [CrossRef]
  12. Schmidt, J.; Ahring, B. Effects of magnesium on thermophilic acetate-degrading granules in upflow anaerobic sludge blanket (UASB) reactors. Enzym. Microb. Technol. 1993, 15, 304–310. [Google Scholar] [CrossRef]
  13. Romero-Güiza, M.; Mata-Alvarez, J.; Chimenos, J.; Astals, S. The effect of magnesium as activator and inhibitor of anaerobic digestion. Waste Manag. 2016, 56, 137–142. [Google Scholar] [CrossRef] [PubMed]
  14. McCarty, P.L. Anaerobic waste treatment fundamentals. Public Works 1964, 95, 107–112. [Google Scholar]
  15. Jard, G.; Jackowiak, D.; Carrère, H.; Delgenès, J.-P.; Torrijos, M.; Steyer, J.-P.; Dumas, C. Batch and semi-continuous anaerobic digestion of Palmaria palmata: Comparison with Saccharina latissima and inhibition studies. Chem. Eng. J. 2012, 209, 513–519. [Google Scholar] [CrossRef]
  16. Fang, C.; Boe, K.; Angelidaki, I. Anaerobic co-digestion of desugared molasses with cow manure; focusing on sodium and potassium inhibition. Bioresour. Technol. 2011, 102, 1005–1011. [Google Scholar] [CrossRef] [PubMed]
  17. Dang, Y.; Zhang, R.; Wu, S.; Liu, Z.; Qiu, B.; Fang, Y.; Sun, D. Calcium effect on anaerobic biological treatment of fresh leachate with extreme high calcium concentration. Int. Biodeterior. Biodegrad. 2014, 95, 76–83. [Google Scholar] [CrossRef]
  18. Chen, Y.; Cheng, J.J. Effect of potassium inhibition on the thermophilic anaerobic digestion of swine waste. Water Environ. Res. 2007, 79, 667–674. [Google Scholar] [CrossRef]
  19. Zhao, W.; Wang, Y.; Liu, S.; Pan, M.; Yang, J.; Chen, S. Denitrification activities and N2O production under salt stress with varying COD/N ratios and terminal electron acceptors. Chem. Eng. J. 2013, 215, 252–260. [Google Scholar] [CrossRef]
  20. Wang, R.; Zheng, P.; Ding, A.-q.; Zhang, M.; Ghulam, A.; Yang, C.; Zhao, H.-P. Effects of inorganic salts on denitrifying granular sludge: The acute toxicity and working mechanisms. Bioresour. Technol. 2016, 204, 65–70. [Google Scholar] [CrossRef]
  21. Phae, C.-G.; Chu, Y.-S.; Park, J.-S. Investigation of effect on composting process and plant growth of salt concentration in food waste. J. Korea Org. Resour. Recycl. Assoc. 2002, 10, 103–111. [Google Scholar]
  22. Lupin, H.; Boeri, R.; Moschiar, S. Water activity and salt content relationship in moist salted fish products. Int. J. Food Sci. Technol. 1981, 16, 31–38. [Google Scholar] [CrossRef]
  23. Bremer, E.; Krämer, R. Responses of microorganisms to osmotic stress. Annu. Rev. Microbiol. 2019, 73, 313–334. [Google Scholar] [CrossRef] [PubMed]
  24. Wood, J.M. Bacterial osmoregulation: A paradigm for the study of cellular homeostasis. Annu. Rev. Microbiol. 2011, 65, 215–238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Csonka, L.N. Physiological and genetic responses of bacteria to osmotic stress. Microbiol. Rev. 1989, 53, 121–147. [Google Scholar] [CrossRef] [PubMed]
  26. Choi, D.-W.; Shin, H.-H.; Kim, J.-G. A study of dewatering phenomena of potato slice cytorrhysed by high molecules. Korean J. Food Nutr. 2006, 19, 358–365. [Google Scholar]
  27. Money, N.P. Measurement of pore size in the hyphal cell wall of Achlya bisexualis. Exp. Mycol. 1990, 14, 234–242. [Google Scholar] [CrossRef]
  28. Griffin, D. Water and microbial stress. Adv. Microb. Ecol. 1981, 5, e136. [Google Scholar]
  29. Rawlins, S.L.; Campbell, G.S. Water potential: Thermocouple psychrometry. Methods Soil Anal. Part 1 Phys. Mineral. Methods 1986, 5, 597–618. [Google Scholar] [CrossRef]
  30. Angelidaki, I.; Alves, M.; Bolzonella, D.; Borzacconi, L.; Campos, J.L.; Guwy, A.J.; Kalyuzhnyi, S.; Jenicek, P.; van Lier, J.B. Defining the biomethane potential (BMP) of solid organic wastes and energy crops: A proposed protocol for batch assays. Water Sci. Technol. 2009, 59, 927–934. [Google Scholar] [CrossRef] [Green Version]
  31. Lay, J.-J.; Li, Y.-Y.; Noike, T. Mathematical model for methane production from landfill bioreactor. J. Environ. Eng. 1998, 124, 730–736. [Google Scholar] [CrossRef]
  32. Unesco. Background Papers and Supporting Data on the Practical Salinity Scale 1978-Unesco; Joint Panel on Oceanographic Tables and Standards: Woods Hole, MA, USA, 1981. [Google Scholar]
  33. Oh, S.-Y.; Yoon, Y.-M. Energy recovery efficiency of poultry slaughterhouse sludge cake by hydrothermal carbonization. Energies 2017, 10, 1876. [Google Scholar] [CrossRef] [Green Version]
  34. Sørensen, A.H.; Winther-Nielsen, M.; Ahring, B.K. Kinetics of lactate, acetate and propionate in unadapted and lactate-adapted thermophilic, anaerobic sewage sludge: The influence of sludge adaptation for start-up of thermophilic UASB-reactors. Appl. Microbiol. Biotechnol. 1991, 34, 823–827. [Google Scholar] [CrossRef]
  35. Hansen, K.H.; Angelidaki, I.; Ahring, B.K. Anaerobic digestion of swine manure: Inhibition by ammonia. Water Res. 1998, 32, 5–12. [Google Scholar] [CrossRef]
  36. Trejo-Téllez, L.I.; Gómez-Merino, F.C. Nutrient solutions for hydroponic systems. In Hydroponics: A Standard Methodology for Plant Biological Research; IntechOpen: London, UK, 2012; pp. 1–22. [Google Scholar]
  37. Sterne, R.; Zentmyer, G.; Bingham, F. The effect of osmotic potential and specific ions on growth of Phytophthora cinnamomi. Phytopathology 1976, 66, 1398–1402. [Google Scholar] [CrossRef]
  38. Oren, A. Bioenergetic aspects of halophilism. Microbiol. Mol. Biol. Rev. 1999, 63, 334–348. [Google Scholar] [CrossRef] [Green Version]
  39. Fougere, F.; Le Rudulier, D.; Streeter, J.G. Effects of salt stress on amino acid, organic acid, and carbohydrate composition of roots, bacteroids, and cytosol of alfalfa (Medicago sativa L.). Plant Physiol. 1991, 96, 1228–1236. [Google Scholar] [CrossRef] [Green Version]
  40. Brown, A. Microbial water stress. Bacteriol. Rev. 1976, 40, 803–846. [Google Scholar] [CrossRef]
  41. Killham, K.; Firestone, M. Salt stress control of intracellular solutes in streptomycetes indigenous to saline soils. Appl. Environ. Microbiol. 1984, 47, 301–306. [Google Scholar] [CrossRef] [Green Version]
  42. Song, C.-W.; Kim, N.-C.; Ryu, J.-K.; Kim, J.-M. A study on the fuelization of livestock sludge using thermal hydrolysis. J. Korea Org. Resour. Recycl. Assoc. 2015, 23, 51–59. [Google Scholar] [CrossRef] [Green Version]
  43. Gunde-Cimerman, N.; Plemenitaš, A.; Oren, A. Strategies of adaptation of microorganisms of the three domains of life to high salt concentrations. FEMS Microbiol. Rev. 2018, 42, 353–375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Cardwell, T.; Cattrall, R. Determination of ammonia in waste waters by a differential pH method using flow injection potentiometry and a nonactin-based sensor. Analyst 1997, 122, 89–93. [Google Scholar] [CrossRef]
  45. Deublein, D.; Steinhauser, A. Biogas from Waste and Renewable Resources: An Introduction; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  46. Lohani, S.P.; Havukainen, J. Anaerobic digestion: Factors affecting anaerobic digestion process. In Waste Bioremediation; Springer: Singapore, 2018; pp. 343–359. [Google Scholar] [CrossRef]
  47. Migliore, G.; Alisi, C.; Sprocati, A.; Massi, E.; Ciccoli, R.; Lenzi, M.; Wang, A.; Cremisini, C. Anaerobic digestion of macroalgal biomass and sediments sourced from the Orbetello lagoon, Italy. Biomass Bioenergy 2012, 42, 69–77. [Google Scholar] [CrossRef]
  48. Sowers, K.R.; Gunsalus, R.P. Halotolerance in Methanosarcina spp.: Role of N (sup (epsilon))-acetyl-(beta)-lysine, (alpha)-glutamate, glycine betaine, and K (sup+) as compatible solutes for osmotic adaptation. Appl. Environ. Microbiol. 1995, 61, 4382–4388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. McGenity, T.J.; Sorokin, D.Y. Methanogens and methanogenesis in hypersaline environments. In Biogenesis of Hydrocarbons; Handbook of Hydrocarbon and Lipid Microbiology; Springer: Cham, Switzerland, 2010; pp. 665–680. [Google Scholar] [CrossRef]
  50. Lai, M.-C.; Sowers, K.R.; Robertson, D.E.; Roberts, M.F.; Gunsalus, R.P. Distribution of compatible solutes in the halophilic methanogenic archaebacteria. J. Bacteriol. 1991, 173, 5352–5358. [Google Scholar] [CrossRef] [Green Version]
  51. Harris, R. Effect of water potential on microbial growth and activity. Water Potential Relat. Soil Microbiol. 1981, 9, 23–95. [Google Scholar] [CrossRef]
  52. Carpita, N.; Sabularse, D.; Montezinos, D.; Delmer, D.P. Determination of the pore size of cell walls of living plant cells. Science 1979, 205, 1144–1147. [Google Scholar] [CrossRef]
  53. Wen, X.; Yang, S.; Horn, F.; Winkel, M.; Wagner, D.; Liebner, S. Global biogeographic analysis of methanogenic archaea identifies community-shaping environmental factors of natural environments. Front. Microbiol. 2017, 8, 1339. [Google Scholar] [CrossRef] [Green Version]
  54. Evans, P.N.; Boyd, J.A.; Leu, A.O.; Woodcroft, B.J.; Parks, D.H.; Hugenholtz, P.; Tyson, G.W. An evolving view of methane metabolism in the Archaea. Nat. Rev. Microbiol. 2019, 17, 219–232. [Google Scholar] [CrossRef]
  55. Mand, T.D.; Metcalf, W.W. Energy conservation and hydrogenase function in methanogenic archaea, in particular the genus Methanosarcina. Microbiol. Mol. Biol. Rev. 2019, 83, e00020-19. [Google Scholar] [CrossRef]
  56. Adams, M.W.; Maier, R.J.; Wiegel, J. Incredible anaerobes from physiology to genomics to fuels. Ann. N. Y. Acad. Sci. 2008, 1125. [Google Scholar] [CrossRef]
  57. Sorokin, D.Y.; Merkel, A.Y.; Abbas, B.; Makarova, K.S.; Rijpstra, W.I.C.; Koenen, M.; Damsté, J.S.S.; Galinski, E.A.; Koonin, E.V.; Van Loosdrecht, M.C. Methanonatronarchaeum thermophilum gen. nov., sp. nov. and Candidatus Methanohalarchaeum thermophilum, extremely halo (natrono) philic methyl-reducing methanogens from hypersaline lakes comprising a new euryarchaeal class Methanonatronarchaeia classis nov. Int. J. Syst. Evol. Microbiol. 2018, 68, 2199. [Google Scholar] [CrossRef] [PubMed]
  58. Xu, R.-z.; Fang, S.; Zhang, L.; Huang, W.; Shao, Q.; Fang, F.; Feng, Q.; Cao, J.; Luo, J. Distribution patterns of functional microbial community in anaerobic digesters under different operational circumstances: A review. Bioresour. Technol. 2021, 341, 125823. [Google Scholar] [CrossRef] [PubMed]
  59. Baek, G.; Kim, J.; Kim, J.; Lee, C. Role and potential of direct interspecies electron transfer in anaerobic digestion. Energies 2018, 11, 107. [Google Scholar] [CrossRef] [Green Version]
  60. Aydin, S. Enhanced biodegradation of antibiotic combinations via the sequential treatment of the sludge resulting from pharmaceutical wastewater treatment using white-rot fungi Trametes versicolor and Bjerkandera adusta. Appl. Microbiol. Biotechnol. 2016, 100, 6491–6499. [Google Scholar] [CrossRef] [PubMed]
  61. Alenzi, A.; Hunter, C.; Spencer, J.; Roberts, J.; Craft, J.; Pahl, O.; Escudero, A. Pharmaceuticals effect and removal, at environmentally relevant concentrations, from sewage sludge during anaerobic digestion. Bioresour. Technol. 2021, 319, 124102. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Cumulative methane yield curves in the water potential (Ψ) inhibition assay.
Figure 1. Cumulative methane yield curves in the water potential (Ψ) inhibition assay.
Fermentation 09 00244 g001
Figure 2. Effect of the water potential (Ψ) on the methane yield (vertical bars mean the standard deviations).
Figure 2. Effect of the water potential (Ψ) on the methane yield (vertical bars mean the standard deviations).
Fermentation 09 00244 g002
Figure 3. Effects of the water potential (Ψ) (a) and K+ concentration (b) on the degree of inhibition of methane production from food wastewater.
Figure 3. Effects of the water potential (Ψ) (a) and K+ concentration (b) on the degree of inhibition of methane production from food wastewater.
Fermentation 09 00244 g003
Figure 4. Order-level taxonomic classification of the microbial community in the water potential (Ψ) inhibition assay.
Figure 4. Order-level taxonomic classification of the microbial community in the water potential (Ψ) inhibition assay.
Fermentation 09 00244 g004
Figure 5. Genus-level taxonomic classification of the microbial community in the water potential (Ψ) inhibition assay.
Figure 5. Genus-level taxonomic classification of the microbial community in the water potential (Ψ) inhibition assay.
Fermentation 09 00244 g005
Table 1. Operating status of the facilities investigated in this study.
Table 1. Operating status of the facilities investigated in this study.
No.Input MaterialsCapacityDigester TemperatureHRT 1Location
(m3 Day−1)(°C)(Days)
1Sewage sludge286337.515.5Gangnam, Seoul
2Pig slurry, food wastewater20036.825.0Goyang, Gyeonggi
3Pig slurry, food wastewater, livestock byproducts15036.030.0Nonsan, Chungnam
4Sewage sludge, food wastewater200037.525.0Dongnae, Busan
5Pig slurry, food waste, sewage sludge32037.035.0Seosan, Chungnam
6Sewage sludge270037.530.0Seongnam, Gyeonggi
7Sewage sludge, food wastewater14034.520.0Ansan, Gyeonggi
8Pig slurry, food wastewater, sewage sludge38436.015.0Anyang, Gyeonggi
9Pig slurry, food waste13036.020.0Yangsan, Gyeongnam
10Pig slurry, food wastewater10040.035.0Yeoncheon, Gyeonggi
11Food waste, slaughterhouse byproducts22040.055.0Wonju, Gangwon
12Pig slurry, food wastewater10038.040.0Icheon, Gyeonggi
13Food wastewater65039.043.0Seo, Incheon
14Food wastewater30036.025.0Jeonju, Jeonbuk
15Pig slurry, food wastewater12538.011.0Seogwipo, Jeju
16Food wastewater15039.526.0Cheongju, Chungbuk
17Sewage sludge, food wastewater52538.030.7Chuncheon, Gangwon
18Food waste8038.019.0Chungju, Chungbuk
19Pig slurry, food wastewater10037.035.0Hongcheon, Gangwon
20Food wastewater12039.530.0Hwaseong, Gyeonggi
1 Hydraulic retention time.
Table 2. Chemical composition of inoculum.
Table 2. Chemical composition of inoculum.
ParameterspHTS 1VS 2TCOD 3SCOD 4TKN 5NH4+-N 6AlkalinityTVFAs 7
(-)(mg L−1)(mg L−1)(mg L−1)(mg L−1)(mg L−1)(mg L−1)(mg L−1 as CaCO3)(mg L−1 as Acetate)
Inoculum7.9630,04414,76721,0506012220113146389113
1 Total solid, 2 Volatile solid, 3 Total chemical oxygen demand, 4 Soluble chemical oxygen demand, 5 Total Kjeldahl nitrogen, 6 Ammonium nitrogen, 7 Total volatile fatty acids.
Table 3. Experimental conditions for the water potential (Ψ) inhibition assay.
Table 3. Experimental conditions for the water potential (Ψ) inhibition assay.
ParametersBlankControlTreatments
K1K2K3K4K5
Inoculum (mL)70707070707070
Substrate 1 (g)00.520.520.520.520.520.52
I/S ratio 2 (-)-0.50.50.50.50.50.5
WP 3
adjusting solution
KCl conc. (%)--0.460.981.954.235.85
Volume (mL)--130130130130130
DW 4 volume (mL)130130-----
Operation volume (mL)200200200200200200200
Ψ 5 (MPa)−0.36 (0.03) 6−0.40 (0.03)−0.57 (0.03)−0.89 (0.03)−1.49 (0.04)−2.85 (0.01)−3.91 (0.02)
K+ conc. 7 (%)N.A. 8N.A.0.240.511.022.223.07
1 D-Glucose used as substrate, 2 Substrate to inoculum ratio (g-VSsubstrate/g-VSinoculum), 3 Water potential,4 Distilled water. 5 Water potential of batch anaerobic reactor, 6 Standard deviation, 7 K+ concentration, 8 Not analyzed.
Table 4. Physicochemical properties of anaerobic digestates (n = 20) investigated in this study.
Table 4. Physicochemical properties of anaerobic digestates (n = 20) investigated in this study.
ParametersMeanMedianMaximumMinimumr 15
Ψ 1 (MPa)−1.23−1.25−0.10−2.09-
pH (-)7.867.838.786.89−0.50 **
WC 2 (wt. %, w.b.)96.2096.3098.5492.130.42 *
TCOD 3 (g L−1)30.9532.3452.4314.30−0.45 *
SCOD 4 (g L−1)9.339.6724.630.00−0.68 ***
TKN 5 (g L−1)4.915.666.951.05−0.87 ***
TAN 6 (g L−1)3.273.305.190.30−0.91 ***
NH4+-N 7 (g L−1)2.872.874.570.30-
Alkalinity (g L−1 as CaCO3)16.8216.5741.601.45−0.61 ***
TVFAs 8(cmolc L−1 as acetate)1.300.5010.340.19−0.12
(mg L−1 as acetate)7783026210113-
Solid
content
TS 9 (g L−1)37.9637.0378.7214.60−0.42 *
VS 10 (g L−1)20.0020.1528.7110.62−0.50 **
FS 11 (g L−1)17.9616.7853.633.98−0.32
TSS 12 (g L−1)31.0427.8564.9015.43−0.06
VSS 13 (g L−1)20.5119.8533.809.63−0.28
FSS 14 (g L−1)15.1810.5338.653.95−0.33
1 Water potential, 2 Water content, 3 Total chemical oxygen demand, 4 Soluble chemical oxygen demand, 5 Total Kjeldahl nitrogen, 6 Total ammonium nitrogen (NH4+-N + NH3-N), 7 Ammonium nitrogen, 8 Total volatile fatty acids, 9 Total solid, 10 Volatile solid, 11 Fixed solid, 12 Total suspended solid, 13 Volatile suspended solid, 14 Fixed suspended solid, 15 Pearson correlation coefficient between the water potential and other parameters (* p < 0.05, ** p < 0.01, *** p < 0.001).
Table 5. Chemical properties of anaerobic digestates (n = 20) investigated in this study.
Table 5. Chemical properties of anaerobic digestates (n = 20) investigated in this study.
ParametersMeanMedianMaximumMinimumr 3
Ψ 1 (MPa)−1.23−1.25−0.10−2.09-
SPS 2 (p.s.u.)15.4916.0126.521.27−0.90 ***
NH4+(cmolc L−1)20.4820.4732.622.16−0.89 ***
(mg L−1)369436935884389-
Na+(cmolc L−1)7.787.3427.920.52−0.51 ***
(mg L−1)178916886419119-
K+(cmolc L−1)5.795.6322.510.16−0.68 ***
(mg L−1)22632202880264-
Mg2+(cmolc L−1)0.380.370.810.03−0.47 **
(mg L−1)46.5844.9798.753.78-
Ca2+(cmolc L−1)1.170.7210.260.20−0.07
(mg L−1)234144205639-
Cl(cmolc L−1)10.949.6832.060.37−0.71 ***
(mg L−1)3879343211,366132-
NO3(cmolc L−1)0.030.000.080.00−0.60 ***
(mg L−1)17.900.0047.220.00-
SO42−(cmolc L−1)0.120.041.140.00−0.49 **
(mg L−1)58.9717.33549.300.00-
PO43−(cmolc L−1)0.880.794.070.00−0.67 ***
(mg L−1)279.07250.351287.000.00-
1 Water potential, 2 Practical salinity, 3 Pearson correlation coefficient between the water potential and other parameters (** p < 0.01, *** p < 0.001).
Table 6. Variation in the concentration of TVFAs and water potential (Ψ) during the water potential inhibition assay.
Table 6. Variation in the concentration of TVFAs and water potential (Ψ) during the water potential inhibition assay.
ParametersFermentation Time (Days)ControlTreatments 1
K1K2K3K4K5
TVFAs 2
(mg/L as acetate)
0105 (0.26) 4104
(0.22)
105
(0.25)
105
(0.09)
105
(0.21)
105
(0.08)
121461 (0.09)1474 (3.26)1448 (7.71)1261 (2.44)1789 (10.08)932
(2.87)
45109 (2.97)112
(0.81)
107
(3.72)
111
(1.04)
1393 (5.99)1357 (5.95)
78116 (0.33)123 (12.55)116
(0.56)
115
(0.42)
266
(1.33)
399 (10.35)
Ψ 3
(MPa)
0−0.37 (0.01)−0.55 (0.05)−0.86 (0.01)−1.46 (0.01)−2.84 (0.00)−3.94 (0.01)
12−0.39 (0.01)−0.55 (0.02)−0.88 (0.00)−1.46 (0.00)−2.84 (0.03)−3.91 (0.04)
45−0.42 (0.02)−0.62 (0.01)−0.92 (0.02)−1.56 (0.01)−2.86 (0.01)−3.92 (0.03)
78−0.43 (0.040−0.58 (0.03)−0.93 (0.02)−1.49 (0.01)−2.86 (0.00)−3.90 (0.04)
1 Water potential was adjusted using KCl solution in the water potential inhibition assay, 2 Total volatile fatty acids, 3 Water potential, 4 Standard deviations.
Table 7. Ultimate methane yield and the parameters calculated with the modified Gompertz model for the water potential (Ψ) inhibition assay.
Table 7. Ultimate methane yield and the parameters calculated with the modified Gompertz model for the water potential (Ψ) inhibition assay.
ParametersControlTreatments 1
K1K2K3K4K5
Ψ 2 (MPa)−0.40−0.57−0.89−1.49−2.85−3.91
VSr 3 (%)70.64 c73.63 b74.55 a74.54 a54.38 d46.21 e
Model
parameters
Bu 4 (Nm3 kg−1-VSadded)0.264 c0.275 b0.278 a0.278 a0.203 d0.172 e
p 5 (mL)127 c132 b134 a134 a98 d83 e
Rm 6 (mL day−1)8.11 a7.92 b5.84 c4.67 e5.48 d1.70 f
λ 7 (days)10.34 c10.16 c7.80 e8.65 d35.96 a25.34 b
R20.9900.9900.9840.9830.9820.985
1 Water potential was adjusted using KCl solution in the water potential inhibition assay, 2 Water potential, 3 Ratio of anaerobic degradation (Bu/Bth × 100; Bth of D-glucose = 0.373 Nm3/kg-VSadded), 4 Methane potential, 5 Maximum methane production, 6 Maximum methane production rate, 7 Lag growth phase time. Mean with different letter differs significantly between treatment (DMRT; Duncan’s multiple range test, p < 0.05).
Table 8. Stimulatory and inhibitory concentrations of K+ on anaerobic digestion.
Table 8. Stimulatory and inhibitory concentrations of K+ on anaerobic digestion.
Reactor TypeTemperatureSubstrateConcentrations (mg L−1)Reference
StimulatoryInhibitory
UASB 135 °CIW 3, SS 4200~40012,000[14]
CSTR 250 °CPS 5-3000[18]
Batch55 °CCM 6-28,000 (50% inhibition)[16]
Batch35 °CMA 7-13,000 (52~53% inhibition)[15]
Batch54 °CSS-8000 (50% inhibition)[11]
1 Upflow anaerobic sludge blanket, 2 Continuous stirred tank reactor, 3 Industrial waste, 4 Sewage sludge, 5 Pig slurry, 6 Cow manure, 7 Macroalgae.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yeo, J.; Kim, C.-G.; Lee, J.-H.; Song, E.; Yoon, Y.-M. Effects of Water Potential on Anaerobic Methane Production and a Microbial Consortium. Fermentation 2023, 9, 244. https://doi.org/10.3390/fermentation9030244

AMA Style

Yeo J, Kim C-G, Lee J-H, Song E, Yoon Y-M. Effects of Water Potential on Anaerobic Methane Production and a Microbial Consortium. Fermentation. 2023; 9(3):244. https://doi.org/10.3390/fermentation9030244

Chicago/Turabian Style

Yeo, Jin, Chang-Gyu Kim, Jun-Hyeong Lee, Eunhye Song, and Young-Man Yoon. 2023. "Effects of Water Potential on Anaerobic Methane Production and a Microbial Consortium" Fermentation 9, no. 3: 244. https://doi.org/10.3390/fermentation9030244

APA Style

Yeo, J., Kim, C. -G., Lee, J. -H., Song, E., & Yoon, Y. -M. (2023). Effects of Water Potential on Anaerobic Methane Production and a Microbial Consortium. Fermentation, 9(3), 244. https://doi.org/10.3390/fermentation9030244

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop