Absence of oxygen effect on microbial structure and methane production during drying and rewetting events

Scientific Reports volume 12, Article number: 16570 (2022) Cite this article

964 Accesses

5 Altmetric

Metrics details

Natural environments with frequent drainage experience drying and rewetting events that impose fluctuations in water availability and oxygen exposure. These relatively dramatic cycles profoundly impact microbial activity in the environment and subsequent emissions of methane and carbon dioxide. In this study, we mimicked drying and rewetting events by submitting methanogenic communities from strictly anaerobic environments (anaerobic digestors) with different phylogenetic structures to consecutive desiccation events under aerobic (air) and anaerobic (nitrogen) conditions followed by rewetting. We showed that methane production quickly recovered after each rewetting, and surprisingly, no significant difference was observed between the effects of the aerobic or anaerobic desiccation events. There was a slight change in the microbial community structure and a decrease in methane production rates after consecutive drying and rewetting, which can be attributed to a depletion of the pool of available organic matter or the inhibition of the methanogenic communities. These observations indicate that in comparison to the drying and rewetting events or oxygen exposure, the initial phylogenetic structure and the organic matter quantity and quality exhibited a stronger influence on the methanogenic communities and overall microbial community responses. These results change the current paradigm of the sensitivity of strict anaerobic microorganisms to oxygen exposure.

Natural environments with frequent drainage are subjected to drying and rewetting events that impose a fluctuation in oxygen (O2) exposure caused by drastic changes in the water table. Strict anaerobes are generally believed to have a low tolerance to O2, but they can apply different strategies to cope with such O2 exposure, such as by forming resting stages or producing enzymes deactivating reactive oxygen species (ROS)1. Additionally, methanogens, although considered O2-sensitive microorganisms, have shown to have a recovery capacity to O2 exposure promoted by desiccation. Methanogens have often been found in oxidative environments such as rice-paddy anoxic soils that are frequently air drained, and antioxidant genes such as the catalase (KAT) gene have also been detected in their genomes, especially in those of class-II methanogens, e.g., Methanosarcina barkeri, Methanobrevibacter arboriphilus, and Methanocellales, also demonstrating aerotolerance2,3. Furthermore, the higher resistance to oxidative stress in Methanocellales, Methanomicrobiales and Methanosarcinales, compared to class-I methanogens (e.g., Methanopyrales, Methanobacteriales and Methanococcales) has been reported2 This resistance might be attributed to their ability to counteract the damage caused by O2 and reactive oxygen species (ROS) via enhanced defence/repair mechanisms, such as a more diverse and higher expression level of antioxidant enzymes and the development of a new O2/ROS elimination pathway2.Thus, it has been recently proposed that the theory concerning O2 intolerance in microbiology should be modified, i.e., that some taxa of strict anaerobic methanogens can survive oxic conditions4. In addition, their function has been shown to rapidly recover, as evidenced by enhanced organic matter degradation and methane production after the drying of freshwater sediments followed by rewetting5. It is unclear what environmental factor might contribute to this recovery, but changes in sedimentary microbial communities were reported after drying and rewetting6, which were distinct for sediments of different origins. This observation led to postulations on the potential effect of different organic sources on the microbial response to desiccation4, as air drying and rewetting could affect the characteristics of organic matter7. Similarly, it has been observed that the mineralization of organic matter was enhanced upon drying and rewetting of soil due to an increased solubilization of organic matter and disruption of soil aggregates8. In contrast, Hernandez et al. (2019) observed that in comparison to the drying and rewetting effect, the soil organic matter effect on the microbial community was minor9. Conrad (2020) proposed that methanogenic insensitivity to desiccation and O2 exposure could be a plausible explanation for these contrasting results and hypothesized that the oxidant resistance of an anaerobic microbial community could be gained or enhanced by their seasonally/periodically drying and rewetting native eco-environments4.

However, understanding how the oxidant resistance of an anaerobic microbial community is developed can be limited by studying the native eco-environments since these possible oxidant resistance mechanisms should have already been established in such environments. It is thus reasonable to assume that natural environments subjected to frequent drainage should contain more O2-tolerant and/or desiccation-tolerant microorganisms4, while anaerobic communities growing in environments that rarely experience aerobic situations should be negatively affected by air-drying and rewetting. Anaerobic digesters (ADs) are engineered systems designed to maintain an anaerobic environment for a relatively long period of time and provide continuous availability of organic matter to maximize organic matter degradation and methanogenesis10. Microbial communities in ADs are not exposed to periodic O2 (or water loss) and therefore are not expected to be adapted to cope with oxidative stress or desiccation. Thus, the drying and high O2 exposure microbial communities experience from ADs should, at least according to the current established knowledge, promote dramatic changes in microbial community structure and functional damage, as shown by their difficulties in resuming methane production. Our results challenged the abovementioned understanding regarding the severe toxicity of O2/oxidants to strict anaerobes, revealing some O2-tolerant methanogens as well as the yet-to-be-explored mechanism of anaerobe oxidant resistance.

Anaerobic digester sludge from the (a) Henriksdal wastewater treatment plant (WWTP) (Stockholm, Sweden), (b) Åby biogas plant with food waste (Linköping, Sweden), and (c) Gasum Jordberga with agricultural waste (Linköping, Sweden) were sampled in October 2017 and transported to the laboratory and kept at room temperature (approximately 25 °C) for a few hours. The sludge was then used as both the medium and inoculum source to investigate the effect of drying (under air and N2) and rewetting on methane production recovery. Control groups without the drying and rewetting treatment for the three sludges were included. Batch bottles for each type of sludge and treatment were set in triplicate.

In this study, methane production was an indirect measure of methanogenic function and was therefore also interpreted as a microbial function in the studied systems. However, the term methane production was chosen for clarity and to avoid misunderstandings. Methane production performance after one drying and rewetting event (cycle 1) was assessed using an Automatic Methane Potential Test System II (AMPTS II, Bioprocess Control, Sweden). 400 mL of sludge was added to each batch bottle while being flushed with N2, leaving 250 mL headspace. All batch bottles were then placed in a thermostat water bath (37 °C) and connected to the AMPTS II system. Mixers were set as stirring in a cycle of 20 min of mixing and 5 min of resting. The produced methane volume was automatically converted to standard temperature and pressure by the software. For air and N2 drying, two trays were used for each type of sludge and each treatment. Each tray was supplied with 1.2 L sludge, closed with a lid, connected to an input hose with N2 or air and an output hose for the ventilation system in an acclimatized room (37 °C). The hoses were placed inside samples, so air/N2 could penetrate deep into the digestates during drying. The drying process lasted around 14 and 20 days for air- and N2-dry treatment, respectively, until the sludge samples were completely dried (evaluated by the decrease in TS weight due to water loss). After drying, all batch bottles were filled with certain amount of dried sludge and rewetted with O2-free Milli-Q water to reach a 400 mL volume of sludge, and they were then flushed with N2 to ensure anaerobic conditions and connected to the AMPTS II as previously described. In parallel, one tray of batch bottles (triplicates) had the same setup except the N2 or air-drying treatment was set as a control group. To facilitate subsequent repeated drying-rewetting and incubating operations, cycles 2 and 3 were performed with a similar drying and rewetting process as described in cycle 1 in a traditional batch system. The batch bottle setup and methane content calculation were described previously by Sun et al.11, except that 160 mL sludge was added to each 320-mL bottle. For each drying treatment, the bottles were opened and flushed with air or N2 in a 37 °C incubation room. For the cycle 2 and 3, sludge samples were completely dry in 5 to 7 days gas flush. There was no stirring during the drying process in all three cycles. Same as cycle 1, the sludge samples were then rewetted with O2-free Milli-Q water and reincubated at 37 °C. One tray of batch bottles (triplicate) without the N2 or air-drying treatment was set up in parallel as a control group. The produced gas in the control group was released to reset the bottle pressure between each cycle. Approximately 2 mL of sludge was sampled from each batch bottle after each drying-rewetting and incubation event as well as from the original inoculum, and these samples were maintained at −20 °C for later molecular analysis.

For cycle 1, total biogas production was determined by using gas meters working based on water displacement. For cycles 2 and 3, the gas pressure in each bottle was regularly measured by a pressure meter (Testo 312, Germany) to quantify the volume of biogas production. Gas samples were collected for headspace methane content analysis by gas chromatography (5880A series, Hewlett Packard, USA). The methane production was normalized to the sludge content of each bottle (i.e., norm. ml g−1 digestate) and reported at standard atmospheric pressure and 0 °C. Before and after each treatment, several parameters of the sludge samples were determined. Dissolved organic carbon (DOC) was analysed by a TOC analyser (Shimadzu TOC-VCPH, Japan) after the samples were filtered through a 0.45 μm polyethersulfone membrane (USA). VFA concentrations, including acetate, propionate, butyrate, isobutyrate, valerate, and isovalerate, were analysed using a gas chromatograph (6890 Series, Hewlett Packard, USA) according to Jonsson and Borén12. The pH was determined using a pH electrode (InfoLab pH 7310, Germany). Student t-test and Welch's ANOVA was performed in R software (version 4.0.2) to analyse the difference of methane production amount and rate across the samples.

Aliquots of 200 mg in triplicate of each saved sample were used to extract total genomic DNA with the FastDNA spin kit for soil (MP Biomedicals, Santa Ana, CA, USA), as described previously11. DNA concentrations were quantified by a Qubit 3.0 Fluorometer (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA). The 16S rRNA gene sequence library was prepared by polymerase chain reaction (PCR) using the primer pair 515′F(GTGBCAGCMGCC GCG GTAA)/805R(GAC TAC HVGGG TAT CTA ATC C) for the bacterial community and 516F(TGY CAG CCG CCG CGG TAA HACCVGC)/915R(GTG CTC CCC CGC CAA TTC CT) for the archaeal community13,14. The detailed PCR procedure and the sequence library preparation steps were described previously15. Next-generation amplicon sequencing was performed using Illumina MiSeq technology at the SNP&SEQ Technology Platform of SciLifeLab in Uppsala, Sweden. The raw data were analysed through the open-source bioinformatics pipeline DADA2 (version 1.16) in R software (version 4.0.2)16. Forward and reverse sequences were cut to lengths of 259 and 150 bp, respectively, for bacteria with the quality threshold of maxEE = 2 and truncQ = 2 and to 280 and 240 bp for archaea with maxEE = 2 and truncQ = 2, according to in silico calculation by FIGARO17. Taxonomic profiles were assigned based on amplicon sequence variants (ASVs) by using the rRNA database SILVA, release 13218.

Software STAMP (version 2.1.3)19 was used to perform ANOVA and Post-hoc tests for analysing the difference in relative abundance on genus level across samples. Nonmetric multidimensional scaling (NMDS) was applied to assess the dissimilarity of the microbial community among the samples (vegan package in R, permutations = 999). Bacterial community diversity was evaluated using the Hill diversity index 0D and 1D20 in the R package ‘hillR’ (https://CRAN.R-project.org/package=hillR). A combined bacterial and archaeal ASV table with a relative abundance ≥ 0.1% of the total sequence threshold was used to conduct co-occurrence network analysis (RCy3 and igraph package in R). Spearman's rank correlation was calculated between each pair of ASVs, and p values were corrected by the Benjamini–Hochberg correction to control the false discovery rate upon multiple comparisons21. Correlations with significance p ≤ 0.001 and coefficients ≥ 0.5 were present in the co-occurrence network. Degree, betweenness, and closeness centrality indices were calculated to evaluate the co-occurrence network structure and potential keystone microorganisms22. The co-occurrence network figure was plotted using Cytoscape (version 3.7.2).

To differentiate methanogenic production and microbial community alterations caused by drying from the effects of O2 exposure, we dried anaerobic sludge under air and N2 atmospheres and rewetted it with O2-free water three consecutive times. Sludge samples were retrieved from three different categories of ADs, i.e., sewage sludge, agricultural wastes, and food waste digesters, which are known to potentially have methanogenic communities with distinctively different phylogenetic structures23,24. After the drying-rewetting events, the methanogenic production and community structure were evaluated. As expected, methane production rates decreased over time and after each drying-rewetting event, with higher rates from agricultural waste followed by those from food waste and sewage sludge (Table 1). However, the overall trend was constant among all experiments, with a quick recovery of methanogenic production after rewetting for all three drying-rewetting events and with an overall lack of a clear difference between drying with air and drying with N2 (Table 1).

While the first drying-rewetting event had a negative effect on the methane production of the three types of sludge, the next two drying-rewetting events generally had significant positive effects on the methane production amount and rate compared to those of the control groups (Table 1). For the food waste and sewage sludge samples after the last treatment, methane production performance became similar for both treatments (air and N2-dry), possibly due to the depletion of organic matter. The ‘exhausted’ sludge (after one drying event) produced more methane than the undried control, suggesting that drying could enhance the availability of organic matter and favour methanogenic production (Table 1). This result is in line with those of several previous studies on the availability and composition of organic matter, including the evaluation of a range of fractionation procedures, such as sample air-drying followed by rewetting and wet-sieving, recognising that air-drying and rewetting of soil samples affect the characteristics of organic matter (OM)7, and increasing OM solubilization and disrupting soil aggregates8, which could increase OM availability to microorganisms. Similar results for methane production have also been reported, i.e., an increase in methane production after drying-rewetting in Amazonian oxbow lake sediments5 and rain-fed paddy soil from Thailand25.

The recovery time of methane production was different between the three different materials and between the first, second, and third drying-rewetting events (Welch's ANOVA, p < 0.01), indicating that the observed different lag phases of methanogenesis were related to the origin of the microbial community and the availability of organic matter. Similar differences have been reported for natural sediments where methanogens have different lag times, depending on their metabolic stage and the presence of electron acceptors other than O25. However, in the present study, a surprisingly similar response between the air and N2 treatments (except for sewage sludge, which produces much less methane than the other two types of sludge and showed an inconsistent trend) suggests that O2 exposure had a limited, if any, influence on the lag phase and recovery of methane production (Table 1).

To assess the effect of drying and rewetting on the solubilization of organic matter, we evaluated the changes in the dissolved organic carbon (DOC) concentrations. The DOC levels varied among the three anaerobic sludge samples, with the highest DOC level occurring in samples from agricultural waste digesters (2526 ± 99 mg/L), followed by those in food waste (1862 ± 18 mg/L) and sewage sludge (293 ± 5 mg/L) digesters. Nevertheless, the DOC levels in the three different sludge types occurred in the same descending order (i.e., agricultural wastes > food waste > sewage sludge) during all drying-rewetting events (Table S1). It is well known that drying can substantially alter the quality and physicochemical characteristics of DOC, e.g., in soil and aquatic sediment samples26,27, but information regarding the drying effect on anaerobic sludge is limited to a few studies28,29,30, with no study evaluating the drying effect on DOC in terms of microbes. Nevertheless, Knoop et al. (2018) observed a decrease in DOC during drying in anaerobic digesters operating with municipal organic waste. Their results also showed a significant positive correlation between the DOC level and dissolved Mg, Ca, Ni, and Zn29.

Mg, Zn, and Ni concentrations affect microbial enzyme activity31, while Ca is involved in the regulation of microbial membrane permeability32; thus, a decreased accessibility of these elements might reduce microbial tolerance to drying. As expected, we observed a DOC decrease after each drying and rewetting event across all sludge types independent of N2 and air drying. A plausible hypothesis for this decrease is that drying stabilizes micropores due to the cohesion of slightly soluble components and strengthens the aggregates, thus causing a reduction in DOC, as was observed for the soil samples33,34. Furthermore, drying can also change the hydrophobicity of particles, and an increase in hydrophobicity in soils can slow the increase in water potential in the following rewetting, thus decreasing the ability of organic matter to redissolve27,35,36.

The archaeal community structures of the three sludge types prior to the experiments (1_AW_I, 1_FW_I, and 1_SS_I) were different, as expected due to their different substrate sources (Fig. 1). Archaeal communities from agricultural waste and sewage sludge digesters were dominated by the obligate acetolactic genus Methanosaeta (43 and 48% of the archaeal communities, respectively), while those from the food sludge digester had the highest relative abundance of Methanoculleus (66.5% of the archaeal communities), followed by Methanosarcina (24.0% of the archaeal communities) (Fig. 1). Methanosaeta can only use acetate, and Methanoculleus can use H2 and formate for methane production, while Methanosarcina is very versatile in its substrate range but cannot use formate for methane production37,38. The first drying-rewetting event promoted an enrichment of Methanobacterium in relation to Methanosarcina in the sludge from the agricultural AD. However, the relative abundance of Methanosarcina significantly recovered and replaced Methanobacterium and Methanosaeta (Fig. S1) after the second and third drying-rewetting events, where an unclassified Bathyarchaeia archaea also became dominant. Among the archaea genera, Methanosarcina had comparatively high relative abundances across all samples from the agricultural and food waste digesters even after the drying and rewetting events, while the relative abundance of Methanobacterium largely decreased. In the archaeal community from the sewage sludge digester, Methanobacterium prevailed, followed by Methanosarcina, after the drying-rewetting events (Fig. 1). Overall, although archaeal communities with different origins developed differently, the drying-rewetting events favoured the establishment of Methanosarcina in relation to other methanogens. Interestingly, no clear difference was seen regarding the compositions of the methanogen communities between the air and N2-dried treatments.

Relative abundance of the archaeal community at the genus level of three types of sludge (AW agricultural waste, FW food waste, and SS sewage sludge) separated after the drying-rewetting (sample name ends with "I") and incubation events (sample name ends with "F"). From left to right of each figure, blocks corresponding to the first, second, and third drying-rewetting and incubation events (sample names start with 1, 2, and 3, respectively). Where the genus name could not be assigned to the sequences, the closest classified taxonomic level is depicted: class (C), order (O), and family (F).

The most unexpected finding of our study was that O2 exposure during the drying and rewetting treatments, scarcely affected the methane production in the later anaerobic incubation step by the three types of sludge, as demonstrated by the similar amounts and rates of methane production after drying with air or N2. A slightly longer lag phase for recovery of methane production by the air-dried archaeal community from the food waste in the first rewetting event and sewage sludge in the third rewetting event were exceptions for the whole experiment, whereas similar lag phases were observed for the other treatment conditions (Table 1). These combined results challenge our current understanding of the O2 sensitivity of methanogens. In addition, the absence of a clear distinction in the archaeal community structure that would be expected between air- and N2-dried anaerobic sludge suggests that the O2 resistance capacity was native even for the methanogenic communities living in strict anaerobic environments. A few groups of methanogens, in this study represented by members of Methanosarcina, showed a higher stress tolerance than that of other methanogens of the archaeal community (Fig. 1). In ADs, compared to other methanogens, Methanosarcina has often been identified as the dominant methanogen due to its high adaptability to variable digester operating conditions, such as low pH and high ammonia and salt content39. A dominance of Methanosarcina also occurs in natural ecosystems, for instance, in Amazonian oxbow lake sediments that undergo periodic rainy and dry seasons5 and biological soil crusts of arid regions experiencing dry and oxic conditions40. A recent study grouped Methanosarcina in class-II among 40 other methanogens regarding oxidant tolerance according to their genome sequences. Furthermore, the study indicated that this high-oxidant tolerance might be attributed to a more diverse and higher expression level of antioxidant enzymes and the development of a new O2/ROS elimination pathway compared to those of class-I methanogens2. Methanosarcina is also flexible in its substrate utilization and grows faster than Methanosaeta41, which contribute to its high tolerance under stressed conditions. Furthermore, the high stress tolerance of Methanosarcina may be further enhanced when harboured in a community instead of as a pure culture42,43. The drying event, however, did not explain the decrease in Methanoculleus (class-II) in the food waste sludge, considering that Methanoculleus has also been shown to be predominant in dry and oxic soils40. A change in the organic matter characteristics after the drying and rewetting events leading to an increase in the availability of acetate in relation to H2 or formate is a possible explanation. This scenario would have favoured the revival of acetolactic Methanosarcina instead of formate- and hydrogenotrophic Methanoculleus in the agricultural and food waste sludge. At the comparably low level of organic matter in the sewage sludge, this effect might have been limited.

The bacterial community structure of the three sludge types had relatively similar compositions but with varying relative abundances (Fig. 2a and b). The phylum Firmicutes accounted for 52 and 48% of the total bacteria from the agricultural and food waste digesters, respectively, followed by Bacteroidetes (25 and 14% of bacteria, respectively). The phyla Cloacimonetes, Chloroflexi, Atribacteria, and Thermotogae were also dominant in the agricultural and food waste digesters, with a few differences, including a higher relative abundance of Thermotogae (16%) in the food waste digester than in the agricultural waste digester and the presence of Synergistets in the agricultural waste digester. The bacterial community from the sewage sludge digester was dominated by the phyla Proteobacteria (25%), Bacteroidetes (22%), Chloroflexi (21%) and Firmicutes (19%). Actinobacteria (6%), Acidobacteria (3%), Aegiribacteria (2%) and Spirochaetes (2%) were also present at low abundances, while they were absent in the samples from agricultural and food waste digesters (Fig. 2b).

Bacterial community structure of the three types of sludge (AW: agricultural waste, FW: food waste, and SS: sewage sludge). (a) Nonmetric multidimensional scaling (NMDS) analysis based on ASV read counts displaying bacterial community distribution coloured by sludge types. (b) Relative abundance of the bacterial community at the phylum level of three types of sludge separated after the drying-rewetting (sample name end with "I") and incubation events (sample name end with "F"). From top to bottom, the bar plots correspond to the agricultural waste, food waste, and sewage sludges. From left to right, the bar plots correspond to the first, second, and third drying-rewetting and incubation events (sample names start with 1, 2, and 3, respectively). Bacteria with relative abundances ≥ 1.0% in at least one sample are depicted.

The first drying-rewetting event altered the structure of the three bacterial communities, with a one order of magnitude substantial increase in the relative abundance of Proteobacteria from < 1 to 53% (air-dried) and 68% (N2-dried) in agricultural waste sludge; from 4% to up to 34% (air-dried) and 13% (N2-dried) for food waste sludge; and from 25% to up to 46% (air-dried) and 37% (N2-dried) for sewage sludge (Fig. 2b). However, the relative abundance of Proteobacteria dropped (< 10%) after the first incubation event, where Firmicutes, particularly in the agricultural and food waste digestate samples, dominated the bacterial community. The same pattern, i.e., alternating dominance between Proteobacteria and Firmicutes during the drying-rewetting and incubating events, was repeated in the second and third rounds of treatments, and at the end, they outcompeted most of the other bacteria [also shown as decreased population diversity indices (Table S2)]. At the genus level, this process occurred as a unified increase in the relative abundance of Hydrogenispora and Alkaliphilus within the phylum Firmicutes (Fig. S2), while for Proteobacteria, the increase in relative abundance differed somewhat according to the origin of the microbial communities, represented by the genus Pusillimonas in the agricultural and food waste samples and an unclassified member of the family Burkholderiaceae in the sewage sludge samples (Fig. S3). In addition, the family Halobacteroidaceae, also belonging to the phylum Proteobacteria, increased at the end of the experiment in the food waste and sewage sludge samples. The genus Hydrogenispora can ferment various sugars to produce hydrogen and acetate, while Alkaliphilus can ferment sugars to acetate and formate44,45. Members of the family Burkholderiaceae and Halobacteroidaceae have also been shown to ferment various sugars46,47, while the genus Pusillimonas can use acetate as a carbon source and potentially compete with methanogens48. The increase in the relative abundance of Pusillimonas after the drying-rewetting events, especially later for the air-dried samples, suggests that acetolactic bacteria might have competed temporally with acetoclastic methanogens (in our case Methanosarcina) for acetate, especially when the latter was inhibited.

Microbial community dynamics under dry and oxic stress in ADs have garnered limited interest, and most studies have focused on the responses in natural environments such as soil and lake sediments where an overall trend of increasing gram-positive phyla (e.g., Firmicutes) in relation to gram-negative phyla (e.g., Proteobacteria and Bacteroidetes) has been reported5,49. The tolerance of bacteria to drying/drought stress has also been relatively well studied. The physiological characteristics influencing tolerance mainly include the thickness of cell walls and the ability of sporulation49. Thus, it is not surprising that Hydrogenispora and Alkaliphilus (belonging to Firmicutes), both spore-forming and gram-positive bacteria, exhibited an overall higher stress tolerance than that of Bacteroidetes (represented by order Bacteroidales, mostly gram-negative) in our study44,50,51,52,53,54,55. However, these arguments do not explain the increase in the abundance of the phylum Proteobacteria, represented by the genera Pusillimonas, Pseudomonas, and Nitrincola and the families Burkholderiaceae and Halobacteroidaceae in the present study, which are all gram-negative and nonspore-forming bacteria, with the exception of a few species in Halobacteroidaceae that can produce endospores46,47,56,57,58. Thus, there must be other bacterial strategies that counter drying effects. One possible explanation could be the ability of these bacteria to produce and accumulate osmolytes, contributing to the maintenance of cellular turgor and protection of macromolecular structure in microorganisms59. Another possibility could be an established network between different microbial groups, i.e., a beneficial cooperative interaction in microbial communities might potentially enhance their environmental stress resistance60.

To assess whether drying-rewetting stress might affect cooperative interactions in microbial communities, we conducted a co-occurrence network analysis on microbial communities after the drying-rewetting and incubation events. Figure 3 displays signification correlations (p ≤ 0.001, R ≥ 0.5) in a microbial network constituted by archaea and bacteria containing at least three nodes (≥ 0.1% in relative abundance) at the genus level based on their degree, betweenness, and closeness centrality. The network pattern clearly differed between the drying-rewetting and recovery/incubation events (Fig. 3a and b). A slightly less clustered and overall lower R value occurred for the drying-rewetting group (Fig. 3a) than for the recovery/incubation group (Fig. 3b). An unclassified genus belonging to the phylum Aegiribacteria and genus Methanolinea was positively linked to each other and showed a relatively high degree of centrality in both networks (Fig. 3, Tables S3 and S4). These taxa, however, occurred with low relative abundances in the samples (Fig. 2). For the drying- and rewetting-treated samples (Fig. 3a), Aegiribacteria was also positively correlated with the genus Candidatus Methanofastidiosum via Bacteriodes vadinHA17 and Anaerolineaceae, which was further correlated with a subcluster of methanogens, including Candidatus Methanomethylicus and Methanosaeta; however, it was negatively correlated with Methanosarcina. Interestingly, the competition relationship between Methanosarcina and the other two methanogens was absent and replaced by a positive interaction between two bacterial genera Hydrogenispora and Alkaliphilus, particularly increasing after the incubation stage (Fig. 3b).

Co-occurrence network analysis based on the correlation of the relative abundance of amplicon sequence variant (ASV) reads for the microbial profiles at the genus level or the closest classified taxonomic level. (a) after the drying-rewetting and (b) incubation events. Bacterial and archaeal groups are shown by nodes in green and orange, respectively. Each edge represents significant correlations between pairs of nodes (p ≤ 0.001), where positive and negative correlations are coloured in green and red, respectively. The thickness of the edge is proportional to the R value of the correlation (R ≥ 0.5).

In general, the network analysis suggested an overall weakened network of interactions between microorganisms after the drying-rewetting events, while a new positive correlation was identified after incubation events among the increased bacteria Hydrogenispora and Alkaliphilus, together with a subcluster of methanogens. Hydrogenispora and Alkaliphilus have also both been found in rice paddies, another environment exposed to periodical drying and rewetting events61,62. Since Hydrogenispora and Alkaliphilus are both identified as fermentative bacteria, producing acetate as well as hydrogen (Hydrogenispora) and formate (Alkaliphilus)44,45, our observations implied that their increase in relative abundance could potentially serve as new substrate suppliers for both formate-dependent, hydrogenotrophic, and acetotrophic methanogens. The correlations between Aegiribacteria and Methanolinea have yet to be explored, except for a recent study reporting their cooccurrence in anaerobic digesters with high sulfide content63. Moreover, Aegiribacteria and Methanolinea showed a relatively high degree of centrality, closeness and connectance in both the drying-rewetting and recovering groups, while their betweenness drastically decreased after the incubation stage (Table S4), suggesting that they might act as keystone species in methanogenic function recovery according to Berry and Widder22.

In contrast to the current established knowledge regarding the severe toxicity of O2/oxidants to methanogens, our observations indicated unexpectedly similar methanogen community structure and methane production recovery after drying/wetting with and without the presence of O2. It is also not clear whether drying promotes the anaerobic microbial community in establishing their strategy to resist O2/oxidants. Therefore, future research should carefully examine genes that are regulated in methanogens under these conditions to further reveal their self-protecting mechanisms under O2/oxidants or drying stress. Our study showed that the original characteristics of the AD sludge, including DOC content and the composition of microbial communities, seemed to be key factors determining drying and rewetting tolerance, as revealed by methane production. Repeatedly recovered methane production after drying and rewetting events implies that microorganisms may form new metabolic interaction networks leading to more stress-tolerant communities.

The raw DNA sequencing data obtained in the current study are available in the National Center for Biotechnology Information database (NCBI) under accession number PRJNA736312, BioSample: SAMN19641899 (Bacteria), SAMN19643011 (Archaea).

Madigan, M. T., Bender, K. S., Buckley, D. H., Sattley, W. M. & Stahl, D. A. (Pearson Benjamin Cummings, San Francisco, CA94111, 2021).

Lyu, Z. & Lu, Y. Metabolic shift at the class level sheds light on adaptation of methanogens to oxidative environments. ISME J. 12, 411–423. https://doi.org/10.1038/ismej.2017.173 (2018).

Article CAS PubMed Google Scholar

Brioukhanov, A. L., Netrusov, A. I. & Eggen, R. I. L. The catalase and superoxide dismutase genes are transcriptionally up-regulated upon oxidative stress in the strictly anaerobic archaeon Methanosarcina barkeri. Microbiology 152, 1671–1677. https://doi.org/10.1099/mic.0.28542-0 (2006).

Article CAS PubMed Google Scholar

Conrad, R. Methane production in soil environments-anaerobic biogeochemistry and microbial life between flooding and desiccation. Microorganisms 8, 881. https://doi.org/10.3390/microorganisms8060881 (2020).

Article CAS PubMed Central Google Scholar

Conrad, R. et al. Response of the methanogenic microbial communities in Amazonian oxbow lake sediments to desiccation stress. Environ. Microbiol. 16, 1682–1694. https://doi.org/10.1111/1462-2920.12267 (2014).

Article CAS PubMed Google Scholar

Ji, Y. et al. Structure and function of methanogenic microbial communities in sediments of Amazonian lakes with different water types. Environ. Microbiol. 18, 5082–5100. https://doi.org/10.1111/1462-2920.13491 (2016).

Article CAS PubMed Google Scholar

Kaiser, M., Kleber, M. & Berhe, A. A. How air-drying and rewetting modify soil organic matter characteristics: An assessment to improve data interpretation and inference. Soil Biol. Biochem. 80, 324–340. https://doi.org/10.1016/j.soilbio.2014.10.018 (2015).

Article CAS Google Scholar

Wu, J. & Brookes, P. C. The proportional mineralisation of microbial biomass and organic matter caused by air-drying and rewetting of a grassland soil. Soil Biol. Biochem. 37, 507–515. https://doi.org/10.1016/j.soilbio.2004.07.043 (2005).

Article CAS Google Scholar

Hernández, M. et al. Structure, function and resilience to desiccation of methanogenic microbial communities in temporarily inundated soils of the Amazon rainforest (Cunia Reserve, Rondonia). Environ. Microbiol. 21, 1702–1717. https://doi.org/10.1111/1462-2920.14535 (2019).

Article CAS PubMed Google Scholar

Schnürer, A., Bohn, I. & Moestedt, J. in Hydrocarbon and Lipid Microbiology Protocols: Bioproducts, Biofuels, Biocatalysts and Facilitating Tools (eds Terry J. McGenity, Kenneth N. Timmis, & Balbina Nogales) 171–200 (Springer Berlin Heidelberg, 2017).

Sun, L., Liu, T., Müller, B. & Schnürer, A. The microbial community structure in industrial biogas plants influences the degradation rate of straw and cellulose in batch tests. Biotechnol. Biofuels 9, 1–20. https://doi.org/10.1186/s13068-016-0543-9 (2016).

Article CAS Google Scholar

Jonsson, S. & Borén, H. Analysis of mono- and diesters of o-phthalic acid by solid-phase extractions with polystyrene–divinylbenzene-based polymers. J. Chromatogr. A 963, 393–400. https://doi.org/10.1016/S0021-9673(02)00647-7 (2002).

Article CAS Google Scholar

Hugerth, L. W. et al. DegePrime, a program for degenerate primer design for broad-taxonomic-range PCR in microbial ecology studies. Appl. Environ. Microbiol. 80, 5116–5123. https://doi.org/10.1128/AEM.01403-14 (2014).

Article ADS CAS PubMed PubMed Central Google Scholar

Takai, K. & Horikoshi, K. Rapid detection and quantification of members of the archaeal community by quantitative PCR using fluorogenic probes. Appl. Environ. Microbiol. 66, 5066–5072 (2000).

Article ADS CAS Google Scholar

Liu, T., Schnürer, A., Björkmalm, J., Willquist, K. & Kreuger, E. Diversity and abundance of microbial communities in UASB reactors during methane production from hydrolyzed wheat straw and lucerne. Microorganisms 8, 1394 (2020).

Article CAS Google Scholar

Callahan, B. J. et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581 (2016).

Article CAS Google Scholar

Weinstein, M. M., Prem, A., Jin, M., Tang, S. & Bhasin, J. M. FIGARO: An efficient and objective tool for optimizing microbiome rRNA gene trimming parameters. bioRxiv https://doi.org/10.1101/610394 (2019).

Article Google Scholar

Quast, C. et al. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucl. Acids Res. 41, D590–D596. https://doi.org/10.1093/nar/gks1219 (2013).

Article CAS PubMed Google Scholar

Parks, D. H., Tyson, G. W., Hugenholtz, P. & Beiko, R. G. STAMP: Statistical analysis of taxonomic and functional profiles. Bioinformatics 30, 3123–3124. https://doi.org/10.1093/bioinformatics/btu494 (2014).

Article CAS PubMed PubMed Central Google Scholar

Lucas, R. et al. A critical evaluation of ecological indices for the comparative analysis of microbial communities based on molecular datasets. FEMS Microbiol. Ecol. https://doi.org/10.1093/femsec/fiw209 (2017).

Article PubMed Google Scholar

Haynes, W. in Encyclopedia of Systems Biology (eds Werner Dubitzky, Olaf Wolkenhauer, Kwang-Hyun Cho, & Hiroki Yokota) 78–78 (Springer New York, 2013).

Berry, D. & Widder, S. Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Front. Microbiol. https://doi.org/10.3389/fmicb.2014.00219 (2014).

Article PubMed PubMed Central Google Scholar

Liu, T., Sun, L. & Schnürer, B. M. A. Importance of inoculum source and initial community structure for biogas production from agricultural substrates. Bioresour. Technol. https://doi.org/10.1016/j.biortech.2017.08.213 (2017).

Article PubMed Google Scholar

Sundberg, C. et al. 454 pyrosequencing analyses of bacterial and archaeal richness in 21 full-scale biogas digesters. FEMS Microbiol. Ecol. 85, 612–626 (2013).

Article CAS Google Scholar

Reim, A. et al. Response of methanogenic microbial communities to desiccation stress in flooded and rain-fed paddy soil from Thailand. Front. Microbiol. 8, 785 (2017).

Article Google Scholar

del Campo, R., Gómez, R. & Singer, G. Dry phase conditions prime wet-phase dissolved organic matter dynamics in intermittent rivers. Limnol. Oceanogr. 64, 1966–1979. https://doi.org/10.1002/lno.11163 (2019).

Article ADS CAS Google Scholar

Borken, W. & Matzner, E. Reappraisal of drying and wetting effects on C and N mineralization and fluxes in soils. Glob. Change Biol. 15, 808–824. https://doi.org/10.1111/j.1365-2486.2008.01681.x (2009).

Article ADS Google Scholar

Awiszus, S., Meissner, K., Reyer, S. & Müller, J. Ammonia and methane emissions during drying of dewatered biogas digestate in a two-belt conveyor dryer. Bioresour. Technol. 247, 419–425. https://doi.org/10.1016/j.biortech.2017.09.099 (2018).

Article CAS PubMed Google Scholar

Knoop, C., Dornack, C. & Raab, T. Effect of drying, composting and subsequent impurity removal by sieving on the properties of digestates from municipal organic waste. Waste Manage. 72, 168–177. https://doi.org/10.1016/j.wasman.2017.11.022 (2018).

Article CAS Google Scholar

Salamat, R., Scaar, H., Weigler, F., Berg, W. & Mellmann, J. Drying of biogas digestate: A review with a focus on available drying techniques, drying kinetics, and gaseous emission behavior. Dry. Technol. https://doi.org/10.1080/07373937.2020.1781879 (2020).

Article Google Scholar

Šafarič, L. et al. Effect of cobalt, nickel, and selenium/tungsten deficiency on mesophilic anaerobic digestion of chemically defined soluble organic compounds. Microorganisms 8, 598 (2020).

Article Google Scholar

Matheri, A. N., Belaid, M., Seodigeng, T. & Ngila, J. C. In Proceedings of the World Congress on Engineering, London, UK.

Kemper, W. & Rosenau, R. Soil cohesion as affected by time and water content. Soil Sci. Soc. Am. J. 48, 1001–1006 (1984).

Article ADS Google Scholar

Peng, X. & Horn, R. Modeling Soil shrinkage curve across a wide range of soil types. Soil Sci. Soc. Am. J. 69, 584–592. https://doi.org/10.2136/sssaj2004.0146 (2005).

Article ADS CAS Google Scholar

Halverson, L. J., Jones, T. M. & Firestone, M. K. Release of intracellular solutes by four soil bacteria exposed to dilution stress. Soil Sci. Soc. Am. J. 64, 1630–1637. https://doi.org/10.2136/sssaj2000.6451630x (2000).

Article ADS CAS Google Scholar

Schimel, J., Balser, T. C. & Wallenstein, M. Microbial stress-response physiology and its implications for ecosystem function. Ecology 88, 1386–1394 (2007).

Article Google Scholar

Smith, K. S. & Ingram-Smith, C. Methanosaeta, the forgotten methanogen?. Trends Microbiol. 15, 150–155. https://doi.org/10.1016/j.tim.2007.02.002 (2007).

Article CAS PubMed Google Scholar

Thauer, R. K., Kaster, A.-K., Seedorf, H., Buckel, W. & Hedderich, R. Methanogenic archaea: Ecologically relevant differences in energy conservation. Nat. Rev. Microbiol. 6, 579–591. https://doi.org/10.1038/nrmicro1931 (2012).

Article CAS Google Scholar

De Vrieze, J., Hennebel, T., Boon, N. & Verstraete, W. Methanosarcina: The rediscovered methanogen for heavy duty biomethanation. Bioresour. Technol. 112, 1–9. https://doi.org/10.1016/j.biortech.2012.02.079 (2012).

Article CAS PubMed Google Scholar

Angel, R., Claus, P. & Conrad, R. Methanogenic archaea are globally ubiquitous in aerated soils and become active under wet anoxic conditions. ISME J. 6, 847–862. https://doi.org/10.1038/ismej.2011.141 (2012).

Article CAS PubMed Google Scholar

Schnürer, A. Biogas production: Microbiology and technology. Adv Biochem Eng Biotechnol. 156, 195–234. https://doi.org/10.1007/10_2016_5 (2016).

Article CAS Google Scholar

Kato, M. T., Field, J. A. & Lettinga, G. High tolerance of methanogens in granular sludge to oxygen. Biotechnol. Bioeng. 42, 1360–1366. https://doi.org/10.1002/bit.260421113 (1993).

Article CAS PubMed Google Scholar

Kiener, A. & Leisinger, T. Oxygen sensitivity of methanogenic bacteria. Syst. Appl. Microbiol. 4, 305–312. https://doi.org/10.1016/S0723-2020(83)80017-4 (1983).

Article CAS PubMed Google Scholar

Liu, Y., Qiao, J.-T., Yuan, X.-Z., Guo, R.-B. & Qiu, Y.-L. Hydrogenispora ethanolica gen. nov., sp. Nov., an anaerobic carbohydrate-fermenting bacterium from anaerobic sludge. Int. J. Syst. Evolut. Microbiol. 64, 1756–1762. https://doi.org/10.1099/ijs.0.060186-0 (2014).

Article CAS Google Scholar

Wu, X.-Y. et al. Alkaliphilus halophilus sp. Nov., a strictly anaerobic and halophilic bacterium isolated from a saline lake, and emended description of the genus Alkaliphilus. Int. J. Syst. Evolut. Microbiol. 60, 2898–2902. https://doi.org/10.1099/ijs.0.014084-0 (2010).

Article CAS Google Scholar

Coenye, T. in The Prokaryotes: Alphaproteobacteria and Betaproteobacteria (eds Eugene Rosenberg et al.) 759–776 (Springer Berlin Heidelberg, 2014).

Oren, A. in Bergey's Manual of Systematics of Archaea and Bacteria 1–4.

Stolz, A., Bürger, S., Kuhm, A., Kämpfer, P. & Busse, H.-J. Pusillimonas noertemannii gen. nov., sp. nov., a new member of the family Alcaligenaceae that degrades substituted salicylates. Int. J. Syst. Evolut. Microbiol. 55, 1077–1081. https://doi.org/10.1099/ijs.0.63466-0 (2005).

Article CAS Google Scholar

Naylor, D. & Coleman-Derr, D. Drought stress and root-associated bacterial communities. Front. Plant. Sci. https://doi.org/10.3389/fpls.2017.02223 (2018).

Article PubMed PubMed Central Google Scholar

Acosta-Martínez, V. et al. Predominant bacterial and fungal assemblages in agricultural soils during a record drought/heat wave and linkages to enzyme activities of biogeochemical cycling. Appl. Soil. Ecol. 84, 69–82. https://doi.org/10.1016/j.apsoil.2014.06.005 (2014).

Article Google Scholar

Bouskill, N. J. et al. Pre-exposure to drought increases the resistance of tropical forest soil bacterial communities to extended drought. ISME J. 7, 384–394. https://doi.org/10.1038/ismej.2012.113 (2013).

Article CAS PubMed Google Scholar

Chodak, M., Gołębiewski, M., Morawska-Płoskonka, J., Kuduk, K. & Niklińska, M. Soil chemical properties affect the reaction of forest soil bacteria to drought and rewetting stress. Ann. Microbiol. 65, 1627–1637. https://doi.org/10.1007/s13213-014-1002-0 (2015).

Article CAS PubMed Google Scholar

Curiel Yuste, J. et al. Strong functional stability of soil microbial communities under semiarid Mediterranean conditions and subjected to long-term shifts in baseline precipitation. Soil Biol. Biochem. 69, 223–233. https://doi.org/10.1016/j.soilbio.2013.10.045 (2014).

Article CAS Google Scholar

Hartmann, M. et al. A decade of irrigation transforms the soil microbiome of a semi-arid pine forest. Mol. Ecol. 26, 1190–1206. https://doi.org/10.1111/mec.13995 (2017).

Article PubMed Google Scholar

Zakharyuk, A. et al. Alkaliphilus namsaraevii sp. nov., an alkaliphilic iron- and sulfur-reducing bacterium isolated from a steppe soda lake. Int. J. Syst. Evolut. Microbiol. 67, 1990–1995. https://doi.org/10.1099/ijsem.0.001904 (2017).

Article CAS Google Scholar

Jin, L. et al. Pusillimonas caeni sp. nov., isolated from a sludge sample of a biofilm reactor. Antonie Van Leeuwenhoek 110, 125–132 (2017).

Article CAS Google Scholar

Joshi, A. et al. Nitrincola tapanii sp. nov., a novel alkaliphilic bacterium from An Indian Soda Lake. Int. J. Syst. Evol. Microbiol. 70, 1106–1111. https://doi.org/10.1099/ijsem.0.003883 (2020).

Article CAS PubMed Google Scholar

Palleroni, N. J. Pseudomonas. Bergey's Manual of Systematics of Archaea and Bacteria. https://doi.org/10.1002/9781118960608.gbm01210 (2015).

Welsh, D. T. Ecological significance of compatible solute accumulation by micro-organisms: From single cells to global climate. FEMS Microbiol. Rev. 24, 263–290. https://doi.org/10.1111/j.1574-6976.2000.tb00542.x (2000).

Article CAS PubMed Google Scholar

Morris, B. E. L., Henneberger, R., Huber, H. & Moissl-Eichinger, C. Microbial syntrophy: Interaction for the common good. FEMS Microbiol. Rev. 37, 384–406. https://doi.org/10.1111/1574-6976.12019 (2013).

Article CAS PubMed Google Scholar

Li, H., Peng, J., Weber, K. A. & Zhu, Y. Phylogenetic diversity of Fe(III)-reducing microorganisms in rice paddy soil: Enrichment cultures with different short-chain fatty acids as electron donors. J. Soils Sedim. 11, 1234. https://doi.org/10.1007/s11368-011-0371-2 (2011).

Article CAS Google Scholar

Yi, X.-Y. et al. Coupling metabolisms of arsenic and iron with humic substances through microorganisms in paddy soil. J. Hazard. Mater. 373, 591–599. https://doi.org/10.1016/j.jhazmat.2019.03.113 (2019).

Article CAS PubMed Google Scholar

Shakeri Yekta, S. et al. Effluent solids recirculation to municipal sludge digesters enhances long-chain fatty acids degradation capacity. Biotechnol. Biofuels 14, 56. https://doi.org/10.1186/s13068-021-01913-1 (2021).

Article CAS PubMed PubMed Central Google Scholar

Download references

The authors would like to thank Professor Bo H. Svensson and Professor Ralf Conrad for providing valuable comments and reviews for this manuscript. This study was funded by the Swedish Energy Agency [Grant number: 35624-2] and conducted within the Biogas Research Centre hosted by Linköping University, Sweden.

Open access funding provided by Linköping University.

These authors contributed equally: Tong Liu and Xiaoxiao Li.

Department of Molecular Science, Uppsala BioCenter, Swedish University of Agricultural Science, 750 07, Uppsala, Sweden

Tong Liu & Anna Schnürer

Laboratory Medicine Department, Putuo Hospital, Shanghai University of Traditional Chinese Medicine, Shanghai, 200062, People's Republic of China

Xiaoxiao Li

Department of Thematic Studies - Environmental Change, Linköping University, 581 83, Linkoping, Sweden

Xiaoxiao Li, Sepehr Shakeri Yekta, Annika Björn & Alex Enrich-Prast

State Key Laboratory of Bioreactor Engineering and Institute of Applied Chemistry, East China University of Science and Technology, Shanghai, 200237, People's Republic of China

Xiaoxiao Li & Bo-Zhong Mu

Multiuser Unit for Environmental Analysis, Universidade Federal do Rio de Janeiro – UFRJ, Av. Carlos Chagas Filho, 373, Bloco A, Ilha do Fundão, Rio de Janeiro, RJ, CEP 21941-971, Brazil

Laura Shizue Moriga Masuda & Alex Enrich-Prast

Biogas Research Center, Linköping University, 581 83, Linköping, Sweden

Tong Liu, Xiaoxiao Li, Sepehr Shakeri Yekta, Annika Björn, Anna Schnürer & Alex Enrich-Prast

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Conceptualization, A.P. and B.M.; methodology, T.L. and A.S. (microbiology), X.L., L.M., and A.P. (anaerobic culturing and drying-rewetting treatment); validation, T.L. and A.P.; formal analysis, T.L. and X.L.; investigation, T.L. and X.L.; resources, A.S. and A.P.; data curation, T.L. and X.L.; writing—original draft preparation, A.P. (abstract and introduction), T.L. (microbiology, methods, results and discussion), X.L. (methods, anaerobic culturing and drying-rewetting treatment); writing—review and editing, T.L., X.L., S.Y., A.B., A.S. and A.P.; visualization, T.L. and X.L.; supervision, A.S. (microbiology) and A.P. (anaerobic culturing and drying-rewetting treatment); project administration, A.P., A.B. and A.S.; funding acquisition, A.P., A.B. and A.S.. All authors have read and agreed to the published version of the manuscript.

Correspondence to Tong Liu or Alex Enrich-Prast.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Liu, T., Li, X., Yekta, S.S. et al. Absence of oxygen effect on microbial structure and methane production during drying and rewetting events. Sci Rep 12, 16570 (2022). https://doi.org/10.1038/s41598-022-20448-5

Download citation

Received: 11 July 2022

Accepted: 13 September 2022

Published: 04 October 2022

DOI: https://doi.org/10.1038/s41598-022-20448-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.