Controls on Soil Organic Matter Degradation and Subsequent Greenhouse Gas Emissions Across a Permafrost Thaw Gradient in Northern Sweden

Warming-induced permafrost thaw could enhance microbial decomposition of previously stored soil organic matter (SOM) to carbon dioxide (CO2) and methane (CH4), one of the most significant potential feedbacks from terrestrial ecosystems to the atmosphere in a changing climate. The environmental parameters regulating microbe-organic matter interactions and greenhouse gas (GHG) emissions in northern permafrost peatlands are however still largely unknown. The objective of this work is to understand controls on SOM degradation and its impact on porewater GHG concentrations across the Stordalen Mire, a thawing peat plateau in Northern Sweden. Here, we applied high-resolution mass spectrometry to characterize SOM molecular composition in peat soil samples from the active layers of a Sphagnum-dominated bog and rich fen sites in the Mire. Microbe-organic matter interactions and porewater GHG concentrations across the thaw gradient were controlled by aboveground vegetation and soil pH. An increasingly high abundance of reduced organic compounds experiencing greater humification rates due to enhanced microbial activity were observed with increasing thaw, in parallel with higher CH4 and CO2 porewater concentrations. Bog SOM however contained more Sphagnum-derived phenolics, simple carbohydrates, and organic- acids. The low degradation of bog SOM by microbial communities, the enhanced SOM transformation by potentially abiotic mechanisms, and the accumulation of simple carbohydrates in the bog sites could be attributed in part to the low pH conditions of the system associated with Sphagnum mosses. We show that Gibbs free energy of C half reactions based on C oxidation state for OM can be used as a quantifiable measure for OM decomposability and quality to enhance current biogeochemical models to predict C decomposition rates. We found a direct association between OM chemical diversity and δ13C-CH4 in peat porewater; where higher substrate diversity was positively correlated with enriched δ13C-CH4 in fen sites. Oxidized sulfur-containing compounds, produced by Sphagnum, were further hypothesized to control GHG emissions by acting as electron acceptors for a sulfate-reducing electron transport chain, inhibiting methanogenesis in peat bogs. These results suggest that warming-induced permafrost thaw might increase organic matter lability, in subset of sites that become wetlands, and shift biogeochemical processes toward faster decomposition with an increasing proportion of carbon released as CH4.

Warming-induced permafrost thaw could enhance microbial decomposition of previously stored soil organic matter (SOM) to carbon dioxide (CO 2 ) and methane (CH 4 ), one of the most significant potential feedbacks from terrestrial ecosystems to the atmosphere in a changing climate. The environmental parameters regulating microbe-organic matter interactions and greenhouse gas (GHG) emissions in northern permafrost peatlands are however still largely unknown. The objective of this work is to understand controls on SOM degradation and its impact on porewater GHG concentrations across the Stordalen Mire, a thawing peat plateau in Northern Sweden. Here, we applied high-resolution mass spectrometry to characterize SOM molecular composition in peat soil samples from the active layers of a Sphagnum-dominated bog and rich fen sites in the Mire. Microbe-organic matter interactions and porewater GHG concentrations across the thaw gradient were controlled by aboveground vegetation and soil pH. An increasingly high abundance of reduced organic compounds experiencing greater humification rates due to enhanced microbial activity were observed with increasing thaw, in parallel with higher CH 4 and CO 2 porewater concentrations. Bog SOM however contained more Sphagnum-derived phenolics, simple carbohydrates, and organic-acids. The low degradation of bog SOM by microbial communities, the enhanced SOM transformation by potentially abiotic mechanisms, and the accumulation of simple carbohydrates in the bog sites could be attributed in part to the low pH conditions of the system associated with Sphagnum mosses. We show that Gibbs free energy of C half reactions based on C oxidation state for OM can be used as a quantifiable measure for OM decomposability and quality to enhance current biogeochemical models to predict C decomposition rates. We found a direct association between OM chemical diversity and δ 13 C-CH 4 in peat porewater; where higher substrate diversity was positively correlated with enriched δ 13 C-CH 4 in fen sites. Oxidized sulfur-containing compounds, produced by Sphagnum, were INTRODUCTION Permafrost peatlands, characterized by their cold and wet conditions that promote the accumulation of soil organic matter (SOM), contain an estimated 277 Pg of carbon (within the 3 m depth) (Schuur et al., 2008), which represents over one-third of the carbon stock in the atmosphere (ca. 800 Pg) (Schuur et al., 2008;McGuire et al., 2009). In the past 100 years, the mean annual Arctic temperatures have increased at almost twice the rate of the global average (IPCC, 2007;Solomon, 2007), leaving the fate of the previously frozen landscape uncertain. Due to these increasing temperatures coupled with seasonal thawing and freezing, permafrost thaw is changing the thickness of the soil active layer (Fenner et al., 2007), surface vegetation composition (Limpens et al., 2008), carbon accumulation rates in peat profiles (McGuire et al., 2009;Frolking et al., 2011), and microbial communities that govern degradation and transformation of organic C, with direct implications on greenhouse gas (GHG) emissions (de Graaff et al., 2010;Deng et al., 2014).
Absence of oxygen, created by water-logged conditions in thawing permafrost can promote anaerobic oxidation of organic matter (OM) (McLatchey and Reddy, 1998) by microbial communities via multiple anaerobic pathways (e.g., hydrolysis, fermentation, and methanogenesis) (Keller and Bridgham, 2007;Sutton-Grier et al., 2011;Wilson et al., 2017). However accurate prediction of GHG emissions from these systems due to changing environmental conditions is limited by our understanding of substrate quality, heterogeneity, and composition (Bölscher et al., 2017), microbial-organic matter interactions and biotic and abiotic OM degradation pathways (Mikutta et al., 2019). Recently, Zhou et al. (2020) found a positive correlation between the richness and diversity of the microbial communities and the concentration and bioavailability of the dissolved organic matter in pit, headstream, and river samples collected from the northern Qinghai-Tibetan Plateau, as well as physiochemical parameters such as pH and DOC and soluble ion concentrations. While some organisms are affected mainly by the concentration of specific substrates, it was found that the overall diversity of organisms may be affected by substrate quality and heterogeneity (Allison, 2005;Hättenschwiler and Jørgensen, 2010), and these effects in turn depend on the prevalence of resource generalists or resource specialists that shapes the microbial communities (Monard et al., 2016;Muscarella et al., 2019). Recently,  investigated microbial community composition and function in the active layer of three sites across a thaw gradient at the Stordalen Mire (an intact palsa, a partially thawed bog and a fully thawed fen) and were able to link the Stordalen genome and functionality to the changing habitat. Metagenomics revealed distinct differences in diversity and overall composition of the microbial communities, with fen having the widest range of genes belonging to different phyla and thus higher microbial diversity. Metatranscriptomic data further confirmed the genomic trends, suggesting that in fen, wider range of microorganisms encode more diverse pathways in which organic substrates such as cellulose and xylan are rapidly degraded, providing substrates needed for methanogenesis and a more heterogenous organic matter composition. On the contrary, microbial communities in bog displayed less diverse functionalities in general with the lowest percentage of cellulaseand xylanase-encoding microorganisms. These results were consistent with previous studies at the same site where Hodgkins et al. (2014) and Hodgkins et al. (2016) showed that dissolved organic matter (DOM) becomes more labile with thaw in the Mire as the pH increases from bog to fen, with bog DOM being more recalcitrant, less diverse and contains more oxygen-rich compounds such as phenols and tannin-like compounds with higher oxidation state and overall lower decomposition rate due to the presence of sphagnan and other acidic inhibitors (Stalheim et al., 2009;Hájek et al., 2011). Conversely, fen DOM was more diverse and contained more low molecular weight compounds, more abundant nitrogen-containing compounds (interpreted as evidence of increased microbial activity), and less organic acids (Hodgkins et al., , 2016. Moreover, differences in microbial community composition and function and OM quality were found to alter pathways by which CO 2 and CH 4 are produced. Studies by Svensson and Rosswall (1984), Galand et al. (2005), Lupascu et al. (2012, and Singleton et al. (2018) suggested a shift from hydrogenotrophic to acetoclastic methanogenesis along the permafrost thaw gradient from bog to fen in the Stordalen Mire.
Although the ecological roles and properties of resource generalists and specialists in many discontinuous permafrost areas have been examined via extensive metagenomic and metatranscriptomics analysis (Liao et al., 2016;Sriswasdi et al., 2017;Singleton et al., 2018;Woodcroft et al., 2018), highly resolved molecular characterization of the OM is an essential foundation upon which we can understand ecosystem C cycling by both plants and microbes in peatlands for predicting future changes in CO 2 and CH 4 production potential (Qualls and Richardson, 2003). Moreover, current models for estimating permafrost soil carbon loss rates are lacking representations of key microbial processes that control OM transformations, leading to enormous uncertainty in model simulations (Burke et al., 2012;Koven et al., 2015;Chang et al., 2019). At most, such models only include a quantitative measure of OM concentration (Lucchese et al., 2010) but neglect the fact that the rate of OM degradation is dependent on the abundance and diversity of active microbes and OM quality and diversity. As such, factors accounting for OM transformation (biotically and abiotically) are poorly constrained in most models and are an important area for model improvement. To address this knowledge gap, we used gas chromatography (GC) and Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry (MS) to characterize SOM chemical composition in bog and fen soil samples collected from the Stordalen Mire. First, we hypothesized that SOM composition and heterogeneity can be used to capture shifts in the functional diversity of microbial communities with thaw gradient. Second, we hypothesized that NOSC and Gibbs free energy of C oxidation can be used to capture C decomposability and changes to the redox conditions of the surroundings (LaRowe and van Cappellen, 2011). Third, we hypothesized that sulfate reduction pathways by sulfate-reducing microorganisms (SRMs) are more active in the bog and are in competition with methanogens, rapidly changing the SOM composition, stimulating CO 2 and suppressing CH 4 emission potentials. Understanding OM energetics and metabolic transformation pathways along the thaw gradient can provide insights on resource heterogeneity and microbial diversity which can then be used to indirectly incorporate microbial processes in different biogeochemistry models and to predict how this ecosystemlevel shift will alter GHG emission under a continually changing climate.

Site Information and Soil Sampling and Preparation
Stordalen Mire (68.35 • N,19.05 • E) is a peat plateau in subarctic Sweden underlain by discontinuous permafrost, which is thawing as the region warms. The Mire has three dominant vegetation community types (Johansson et al., 2006): (1) Palsa (intact permafrost) with woody herbaceous vegetation; (2) Ombrotrophic peatland or bog (intermediate thaw) features with Sphagnum spp., sedges, and shrubs as the dominant vegetations; and (3) Minerotrophic peatland or fen (fully thawed) dominated by sedges, mainly Eriophorum spp. (Bäckstrand et al., 2010). Generally, bogs are fed only by precipitation, have acidic (pH ∼ 4) surface waters and Sphagnum-dominated vegetation. These Sphagnum mosses, are known to inhibit alternate plant growth by acidifying water and creating an environment that is harsh to other plant species (van Breemen, 1995). Fens, on the other hand, are characterized by more mineral-rich water because they are fed by surface water and groundwater, in addition to rainfall. They are usually dominated by grasses and sedges and maintain slightly acidic to alkaline pH. Due to continued warming in the Arctic, the areal extent of intact palsa (Malmer et al., 2005) and Sphagnum community (Norby et al., 2019) across the Stordalen Mire has declined significantly since 1970 (Malmer et al., 2005) while at the same time, wetter fen habitats have expanded (Kokfelt et al., 2009;Bäckstrand et al., 2010). The ground's thaw state and consequent relationship to the water table determine the plant community, which in turn determines organic matter composition. In areas where the permafrost is degraded, land subsidence leads to varying degrees of water inundation. In this study, we focused on these inundated sites, and specifically a fen and a bog site with an extensive set of historical gas flux, microbial, and geochemical measurements (previously described by Hodgkins et al., 2014Hodgkins et al., , 2016McCalley et al., 2014;Woodcroft et al., 2018). In June 2012, peat was collected from both sites with an 11-cm-diameter homemade circular push corer. Both sites were cored in triplicate to capture spatial heterogeneity. Cores were divided into three sections, and the sections were placed into plastic bags and stored frozen until analysis. Eighteen soil samples were collected in total, including nine samples (3 cores × 3 depths per core) collected from a Sphagnum-derived ombrotrophic bog and nine samples collected from a sedgedominated fen ( Table 1). Porewater for dissolved organic carbon (DOC) analysis was gathered from the same depths by syringe suction through a 1-m long, 0.5-cm-diameter stainless steel tube with holes drilled along the bottom 3 cm, then filtered through 0.7-µm Whatman GF/F glass microfiber filters into 120-mL brown borosilicate bottles. The pH was measured on-site with an Oakton Waterproof pHTestr 10. The bottles were frozen within 8 hours of collection and stored frozen until analysis. DOC concentrations were measured by high-temperature catalytic oxidation on a Shimadzu Total Organic Carbon analyzer with a non-dispersive infrared detector. Each sample was analyzed in triplicate measurements, which constantly had a coefficient of variance of <2%. The concentration was reported as millimolar (mM). Sulfate and nitrate concentrations were measured by ion chromatography on a Dionex ICS-1100 with a 4-mm IonPac AS22 column, with a 4.5 mM carbonate/1.4 mM bicarbonate eluent at a flow rate of 1.2 mL/min. The concentration was reported as micromolar (µM). In this study, we will focus on the overall SOM molecular composition differences observed in the bog vs. the fen and not the depth effect.

Porewater Gas Measurements
Porewater samples for dissolved gas analysis were collected similarly to the DOC samples, with porewater drawn up by syringe suction through perforated stainless-steel tubes that were inserted into the peat at the desired depth. Porewater was immediately filtered in the field using 2.7-µm Whatman GF/D glass-fiber filters, then stored in pre-evacuated glass vials sealed with butyl stoppers. Samples were then processed according to the method described by (Hodgkins et al., 2015). Phosphoric acid (1 mL of 20%) was added to each sample, to convert all dissolved inorganic carbon to CO 2 (thus removing differences in HCO − 3 /CO 2 ratios due to differing pH) and to preserve for shipment to Florida State University. Helium was then added to the sample headspaces to bring them to a constant pressure of 1 atm. Samples were analyzed for CH 4 and CO 2 concentrations and stable isotopic composition (δ 13 C) on a Thermo Finnigan Delta-V Isotope Ratio Mass Spectrometer. As described in (Corbett et al., 2013(Corbett et al., , 2015 and , headspace gas concentrations were then converted to dissolved concentrations (mM) in the original unacidified porewater based on previously-determined gas extraction efficiencies (defined as the proportion of formerly-dissolved gas in the headspace after acidification and equilibration between the water and headspace) (Supplementary Table 1). Each sample was analyzed twice and the average results for each sample were recorded. Analytical precision of replicate δ 13 C measurements was 0.2‰.

Organic Matter Extraction
Water was used to extract the labile fraction of SOM from the collected peat (n = 18, 9 fen and 9 bog samples) according to the protocol described by Tfaily et al. (2015) and Tfaily et al. (2017). Briefly, extracts were prepared by adding 1 mL of H 2 O to air dried and ground 100 mg bulk soil and shaking in 2 mL capped glass vials for 2 h on an Eppendorf Thermomixer. Samples were then removed from the shaker and left to stand before spinning down and pulling off the supernatant to stop the extraction. All samples were stored at −80 • C until analysis. Water extractable organic matter (WEOM; Corvasce et al., 2006) is described as a bioavailable fraction of organic matter that contains mostly polar compounds. Changes in the molecular weight and elemental composition of the WEOM fraction can provide important information regarding ecologically relevant processes in soil. Around 0.5 ml of water extracted organic matter from nine bog samples and only three fen samples (due to sample limitations) were first freeze dried using a Labconco TM FreeZone Freeze Dryer and then derivatized as described previously (Kim et al., 2015). Briefly, 20 µl of 30 mg/ml methoxyamine hydrochloride in pyridine (Sigma-Aldrich, Saint Louis, MO) was added to each freeze-dried sample for the protection of carbonyl groups. Samples were vortexed for 20s, sonicated for 60s, and were incubated at 37 • C for 90 min with generous shaking (1,000 rpm). After the incubation, 80 µL of N-methyl-N-(trimethylsilyl) trifluoroacetamide (MSTFA) (Sigma-Aldrich) with 1% trimethylchlorosilane (TMCS) (Sigma-Aldrich) was added to each sample. The samples were vortexed for 20s, sonicated for 60s, and incubated for 30 min at 37 • C with shaking (1000 rpm). After derivatization, the samples were allowed to cool to room temperature.

GC-MS Analysis
Samples were transferred to autosampler vials and analyzed by GC-MS in random order. Blanks and fatty acid methyl ester (FAME) samples (C8-28) were also included in the analysis for the background referencing and retention time calibration, respectively. Metabolites were separated using a HP-5MS column (30 m × 0.25 mm × 0.25 µm; Agilent Technologies, Santa Clara, CA) and analyzed in an Agilent GC 7890A coupled with a single quadrupole MSD 5975C (Agilent Technologies, Santa Clara, CA). For each analysis, 1 µl of sample was injected in splitless mode. The injection port temperature was kept constant at 250 • C for the duration of the analysis. The GC analysis started with an initial oven temperature fixed at 60 • C for 1 min after injection, and then at a steady rate of 10 • C/min, the temperature was increased to 325 • C with the analyses finishing with a 5 min hold at 325 • C. The GC-MS raw data files were processed with MetaboliteDetector as previously reported by Hiller et al. (2009). Briefly, retention indices (RI) of the detected polar metabolites were calibrated based on the mixture of FAME internal standards, followed by deconvolution, and subsequent chromatographic alignment of the detected features across all analyses. Identification of the metabolites was carried out by first matching GC-MS spectra and retention indices against a Pacific Northwest National laboratory (PNNL) augmented version (>900 metabolites) of FiehnLib (PMID: 19928838) library containing validated retention indices and spectral information (Kind et al., 2009). Subsequently, unidentified metabolites were screened against the NIST14 GC-MS library by spectral matching alone. To eliminate false identification and reduce the deconvolution errors created during the automated data-processing, all metabolite identifications were individually validated.
Row-wise normalization (by sum) was then performed on the metabolite data using MetaboAnalyst 4.0 to allow for general-purpose adjustment for differences among samples. Data transformation [Generalized log transformation (log 2 )] and data scaling [Auto scaling (mean-centered and divided by standard deviation of each variable)] were further performed to make features more comparable using the same platform.

FOURIER-TRANSFORM ION CYCLOTRON RESONANCE MASS SPECTROMETRY (FT-ICR-MS) Direct Inject FT-ICR-MS Analysis
A 12 Tesla Bruker FT-ICR mass spectrometer located at PNNL was used to collect high resolution mass spectra of the WEOM. A Suwannee River Fulvic Acid standard (SRFA), obtained from the International Humic Substance Society (IHCC), was used as a control to monitor potential carry over from one sample to another and ensure instrument stability in day to day operation. The instrument was flushed between samples using a mixture of Milli-Q water and HPLC grade methanol. The ion accumulation time (IAT) was varied between 0.03 and 0.05 s to account for variations in C concentrations in different samples. For each sample, 144 individual scans were averaged, and an organic matter homologous series separated by 14 Da was used as an internal calibration. The mass measurement accuracy across 100-1,200 m/z range was <1 ppm (accounted for singly charged ions) and the mass resolution at 341 m/z was ∼240 K and the transient was 0.8 s. BrukerDaltonik version 4.2 data analysis software was used to convert raw spectra (obtained from each sample) to a list of m/z values applying FT-MS peak picker module with a signal-to-noise ratio (S/N) threshold of 7 and absolute intensity threshold of 100. Formularity  software was used to assign chemical formulae based on the criteria described in Tfaily et al. (2018).

Chemical Compounds Classification and Data Processing
The molecular formulae in each sample were evaluated on van Krevelen diagrams (Van Krevelen, 1950), based on their H:C ratios (y-axis) and O:C ratios (x-axis) (Kim et al., 2003), assigning them to the major biochemical classes (e.g., lipid, protein, lignin, carbohydrate etc.). The H:C and O:C ranges for biochemical classification are provided in Table 2. Oxidation state is related to the bioenergetic potential of a specific compound based on its atomic makeup. Oxidation state indicates how a compound may be transformed in biochemical processes without necessarily knowing its structure (LaRowe and van Cappellen, 2011). For each individual compound, the average nominal oxidation state of carbon (NOSC), which represents the number of electrons transferred in C oxidation half reactions, was calculated based on Equation 1 (where C, H, N, O, P, and S represent the numbers of each element, and number coefficients represent typical oxidation states). More negative values for NOSC indicates the presence of more reduced carbon atoms.
(1) G • C-ox (kJ/mol C) indicates the thermodynamic favorability of the compound, with lower values of G • C-ox for a compound indicating greater likelihood of degradation (LaRowe and van Cappellen, 2011). G • C-ox was estimated from NOSC following Equation (2): Positive values for G • C-ox indicate that the oxidation of the C must be coupled to the reduction of a terminal electron acceptor, which in the TEA-depleted mire can include the soil organic matter itself (Keller and Takagi, 2013;Wilson et al., 2017). Organic matter degradation may become thermodynamically limited depending on the availability of terminal electron acceptors (Boye et al., 2017). Under oxic conditions the availability of high energy yielding O 2 as a terminal electron acceptor couples the decomposition of organic matter of all naturally occurring NOSC values (Keiluweit et al., 2016). However, under the inundated conditions of subsurface peat deposits, oxygen and other terminal electron acceptors become depleted thereby slowing or completely inhibiting the decomposition of organic matter with low NOSC values (Boye et al., 2017). In peat, oxygen is rapidly depleted to below detection limits (20 ppb) within 2-3 cm below the water surface (Lin et al., 2014). At Stordalen Mire, the average annual water table is typically near the peat surface in the fen  and was 3.5-7 cm above the peat surface during our sampling. Thus, even the shallow fen peat is consistently under anoxic conditions and subject to depleted TEA availability.
In the bog, the mean annual water table depth is around 17 cm , but was only 2 cm below the peat surface during sampling for this project. Thus, the bog peat is subjected to fluctuating water (and subsequently oxygen) levels throughout the growing season. This could have the effect of enhancing decomposition rates during dry periods and re-oxidizing previously reduced TEAs to support decomposition following subsequent inundation. Relative energy yields associated with the terminal electron acceptors in anoxic conditions also provide insight on the sequences in which organic matter gets decomposed (Grandy and Neff, 2008). Lower G • Cox is an indication of compounds more thermodynamically favored to undergo metabolic and chemical reactions ) (i.e., available for microbial communities), whereas higher G • C-ox reflects microbially degraded and transformed (i.e., recalcitrant) OM. For each elemental formula, double-bond equivalent (DBE) minus oxygen (DBE-O) was calculated using Equation (3): While NOSC and G • C-ox provide information on overall quality of the OM, DBE-O provides more structural information, mostly of carbon-to-carbon double bonds (i.e., lack of hydrogen saturation) without considering the oxygen atoms that mostly exist as carboxyl groups in OM (Bae et al., 2011). Using this formulation, more oxygenated molecules have lower (negative) values for DBE-O while higher values for DBE-O indicate highly unsaturated compounds such as condensed aromatic structures (Herzsprung et al., 2014).

Identification of Metabolic Transformations in Bog and Fen Samples Using Network Analysis
The rich database generated by high resolution mass spectrometry provides accurate mass information, enabling us to infer metabolite composition in complex matrices. We used network analysis to explore metabolic transformations in the bog and fen samples to further understand how two compounds are related to one another and how they might be involved in the same metabolic transformation (via a mass difference between two m/z values demonstrated as an edge in the network, and the initial and final compounds as nodes). Here, we used molecular transformation analysis to identify potential microbial transformation or decomposition pathways using FT-ICR-MS data. For a reaction to be considered and hence counted by this method, both the reactant and the product should be present, hence this approach can provide some information about SOM decomposition, uptake and accumulation (Burgess et al., 2017). We used Cytoscape (Shannon et al., 2003;Kohl et al., 2011) with the MetaNetter plug-in to generate molecular interaction networks. For each sample, the list of m/z values generated by the FT-ICR-MS analysis along with a list of 86 pre-defined m/z values and their subsequent molecular formula and IUPAC names [adopted from (Breitling et al., 2006) (Supplementary Table 2)] were imported into Cytoscape. The molecular interaction networks are consisting of nodes (m/z values) and edges (mass difference between two m/z values). Each m/z value inferred from the FT-ICR-MS analysis was defined as a node (described previously in Longnecker and Kujawinski, 2016) connected to other node (s) through the pre-defined m/z values defined as edges, representing a metabolic transformation. For example, gain or loss of CH 2 ( m = 14.01565007) between two compounds in the molecular interaction networks is depicted by two nodes and an edge. Nodes represent the two compounds with mass difference of 14.01565007 Da and edge is the mass difference (named CH 2 ). In total, 86 metabolic transformations were grouped into four different categories: amino acid-associated transformation (loss/gain of an amino acid and/or involvement of nucleoside or nucleotide), sugar -associated transformation (loss/gain of monomeric or polymeric carbohydrates), carbon transformation (reactions that can occur both biotically or abiotically depending on the environmental conditions), and other transformation that includes heteroatoms (such as loss/gain of a phosphate or an amine group). For this study, we have only focused on the first three categories since similar trends for the bog and fen samples (regardless of depth) were observed in the 'other transformation' category. The output of the Cytoscape is a list with the number of times (count) each transformation occurred between any two peaks in a particular sample as shown in Supplementary Figure S5. Normalized values (to the total number of transformations observed in each sample) and subsequently, normalized percentages are then calculated and used for comparison between samples and across sites.
Furthermore, a comprehensive set of topological parameters calculated from these undirected network plots (Assenov et al., 2008) based on transformation data generated by the Cytoscape with the MetaNetter plug-in was used to compare SOM metabolic transformation pathways between the bog and fen samples. For this study, two different parameters were chosen: clustering coefficient and average number of neighbors, all of which can be used to represent the level of population of the transformation network (Barabási and Albert, 1999;Ravasz et al., 2002). The clustering coefficient, C n , in undirected networks is a ratio of the number of edges (N) between the neighbors of a node n, and the maximum number of edges (M) between the neighbors of n ( Barabási and Oltvai, 2004). C n is calculated as C n = 2e n /(k n (k n -1)), where k n is the number (size) of neighbors of n and e n is the number of connected pairs between all neighbors of n (Watts and Strogatz, 1998). The clustering coefficient of a node is always between 0 and 1( Barabási and Oltvai, 2004). A high clustering coefficient indicates a high number of neighbors for each node and hence represents low selectivity and high diversity since a compound could be produced by multiple reactions. The average number of neighbors therefore is an average of k n of all nodes in the undirected network (Dong and Horvath, 2007). High number of neighbors indicates low degradation and high selectivity. Values for clustering coefficient and average number of neighbors for each sample are provided in Supplementary Table 3. Moreover, the network analysis parameters (clustering coefficient and average number of neighbors) were used to correlate the overall number and/or density of the metabolic reactions happening in each sample to measured biogeochemical properties in pore water such as pH, δ 13 C and DOC.

GC-MS Data
Univariate (Pearson correlation coefficient, t-test significance), multivariate (principal component analysis, PCA) and clustering analysis (heatmap, t-test significance) were performed on normalized values to visualize differences between treatments and among samples using MetaboAnalyst 4.0 and plotted using JMP 14.2.0 and R (R Core Team, 2019) package ggplot2 version 3.3.1 (Wickham, 2016).

Changes in SOM Molecular Composition Associated With Permafrost Thaw
Principal component analysis of the relative abundance of different SOM compound classes (as inferred by FT-ICR-MS), biogeochemical properties of peat porewater (TN concentration, DOC, sulfate, nitrate, and ammonia concentration and pH) and peat porewater CO 2 and CH 4 concentrations revealed clear separation between the bog and fen sites along PC1 (which explained 41% of the variance), with higher PC1 scores in rich fens (Figure 1). SOM in the Sphagnum-dominated bog sites had a higher relative abundance of tannins, lignins, unsaturated, and condensed hydrocarbons (UnSatHC and CondHC, respectively), and lipids compared to the fen sites. In contrast, fen samples were separated from bog samples along PC1 which correlated positively with the abundance of peptide-, carbohydrates (cellulose-)-, and aminosugar-like (AS) compounds. PC2 explained a small proportion of the variance (16%). Higher concentrations (measured in mM) of TN and DOC were observed in the bog porewaters compared to the fen porewater. In contrast, sulfate, nitrate and ammonia concentrations were higher in fen porewater and correlated positively with peat porewater pH. Porewater CO 2 and CH 4 concentrations were higher in the fen sites. Although some depth trends were observed, these trends were generally weak, as they were not consistently observed over the entire depth range and none of the principal components showed significant correlations with depth (Supplementary Figure 1).
To better visualize the differences in the bog and fen SOM molecular composition, unique peaks present in each group (bog FIGURE 1 | Principal components analysis on elemental composition classes from FT-ICR-MS data (black arrows)and biogeochemical properties including the concentrations of total nitrogen(TN), dissolved organic carbon (DOC), sulfate, nitrate, ammonia, and pH (measured in porewater, indicated by blue arrows) and CO 2 and CH 4 (values for each sample is described in Supplementary Table 1). Fen samples are clearly separated from bog along PC2 (which positively correlates with peptides, CO 2 , CH 4 , sulfate, pH, and nitrate) and bog samples show positive and high correlations with carbohydrates, tannins and condensed hydrocarbon and DOC.
vs. fen, regardless of the depth) were reported separately on van Krevelen diagrams (Figures 2A,B). SOM in bog samples (Figure 2A) appeared to be more rich in lignin-like (0.29 < O/C < 0.65 and 0.7 < H/C < 1.5), tannin-like (0.65 < O/C < 1 and 0.8 < H/C < 1.5), and lipid-like (0 < O/C < 0.3 and 1.5 < H/C < 2.5) compounds and fen SOM (Figure 2B) was more rich in peptide-like (0.3 < O/C < 0.6 and 1.5 < H/C < 2.3) compounds and amino sugars (0.5 < O/C < 0.7 and 0.8 < H/C < 2.5). Out of 6,906 peaks detected in the bog samples in total, 24% were unique to bog, while only 8% of the peaks (out of 5,973 in total) detected in the fen were unique to fen, suggesting a higher number of unique compounds in bog (potentially not yet degraded) compared to fen samples (Venn diagram in panel B). Figure 3A shows G • C-ox calculated for 6906 assigned molecular formula in nine bog samples and 5973 assigned molecular formulas in nine fen samples. Lower G • C-ox and high NOSC represent thermodynamically more labile organic compounds and vise-versa. For each sample, G • C-ox values ( Figure 3A; Table 3) and NOSC values ( Table 3) were averaged across all formulas and unpaired t-test was performed on the averaged values obtained from nine bog samples and nine fen samples. Average G • C-ox values in bog samples were lower than fen samples ( Figure 3A; Table 3) and subsequently, average NOSC values were higher in bog samples compared to fen Oxygen to carbon (O/C) and hydrogen to carbon (H/C) ratios were calculated for unique peaks in bog (green) and fen (blue) sites. Ven diagram (B) indicates that 24% of the peaks (out of 6,906 peaks in total) detected in bog samples are unique to bog and 8% of the peaks (out of 5,973 peaks in total) detected in fen are unique to fen.

FIGURE 3 | (A)
Average G • C-ox (Gibbs free energy of C half reactions based on C oxidation; kJ/mol C) values based on FT-ICR-MS data calculated for the bog (green) and fen (blue) samples (regardless of depth). Low G • C-ox represents thermodynamically favorable organic compounds for microbial degradation. Unpaired t-test was performed on the calculated values to determine the significant differences and P-values were reported. Additionally, (B) G • C-ox values were graphed based on compound classifications for bog and fen samples. Statistical significances were determined by unpaired t-test. The "ns" symbol indicates P > 0.05. All P < 0.05, 0.001, and 0.0001 are summarized with one, three and four asterisks, respectively, indicating G • C-ox values for those classes were significantly different between the bog and fen samples. samples (Table 3). Interestingly, different G • C-ox values were obtained for bog and fen samples within the same class of chemical compounds (Figure 3B). While amino sugars were more abundant in the fen sites compared to bog sites, they were characterized with higher average G • C-ox compared to those in the bog sites. On the other hand, while bog sites had an abundance of lignin-like and tannin-like compounds compared to fen, however, they were of lower average G • C-ox compared to those in the fen sites (more oxidized). The lower overall NOSC values of the bog SOM relative to the fen SOM could reflect the influence of oxygen at the bog site. As noted in section Chemical Compounds Classification and Data Processing, the

S/C and values for NOSC and Gibbs free energy ( G • C-ox in kJ/mol C) for unique peaks [formulas that were detected in either bog or fen samples but not both (regardless of depth and replicate number)]. For each value, average of nine samples for bog and nine samples for fen as well as standard deviations were calculated.
fen is more consistently inundated (hence anoxic), while the bog habitat experiences alternating periods of oxic and anoxic conditions. These periods of exposure to oxygenation and the potential re-oxidation of organic and inorganic TEAs may facilitate decomposition at the bog site resulting in SOM with an overall lower NOSC value relative to the fen. However, the overall lower NOSC in the bog site may also reflect different starting substrates.

Sulfur and Nitrogen Dynamics With Permafrost Thaw
To further investigate the sulfur and nitrogen dynamics in the Mire, composition of S-containing SOM (regardless of depth) was compared in the bog (A) and fen (B) samples and unique peaks were reported on a van Krevelen diagram (Figures 4A,B). Sulfur containing compounds in bog fell in the lipid-like (0 < O/C < 0.3 and 1.5 < H/C < 2.5), peptide-like (0.3 < O/C < 0.6 and 1.5 < H/C < 2.3), and amino sugar-like (0.5 < O/C < 0.7 and 0.8 < H/C < 2.5) regions of the van Krevelen diagram. In contrast, Scontaining compounds in fen fell in the condensed hydrocarbonlike (0 < O/C < 0.4 and 0.2 < H/C < 0.8) and ligninlike (0.29 < O/C < 0.65 and 0.7 < H/C < 1.5) regions of the van Krevelen diagram. Overall, a higher abundance of unique S-containing organic compounds was observed in the bog compared to the fen system (∼ 2.6 to ∼1.3% unique S containing compounds, respectively, depicted on Venn diagram in panel A) coincident with lower sulfate concentrations in bog porewater (Table 1). Furthermore, molecular weights and nitrogen to carbon ratios (N/C) for unique compounds were calculated, and both were overlaid on the van Krevelen diagram (Figures 4C,D). Higher abundance of N-containing compounds (∼5% −345 unique compounds out of 6,906 peaks in total) with molecular weight >400 m/z (>400 Da) was observed in bog compared to fen samples. Conversely, fen SOM appeared to have a lower abundance of N-containing compounds (1% −70 unique compounds out of 5,973 peaks in total) and the Ncontaining compounds in the fen tended to be of lower mass (Supplementary Figure 2). DBE-O was calculated for all unique fen and bog S-and N-containing compounds to better characterize the OM in the bog and fen samples (Supplementary Figure 3)

Shifts in Primary Metabolite Metabolism
We used multiple statistical tests to look at changes in SOM primary (central) metabolite composition between bog and fen samples (Figure 5; Table 4) where we observed statistically significant differences in bog and fen primary metabolite abundance. Principal component analysis ( Figure 5A) revealed distinct differences in metabolite abundance (explained by PC1, 47.1%) between bog and fen samples. Effect of depth was also explained by PC2 (19.9%). Correlation coefficients of 25 primary metabolite abundances between bog and fen samples (regardless of depth) using Pearson r distance measure were determined at p < 0.05 ( Figure 5B). Sugar acids such as D-galacturonic acid and simple sugars such as D-glucose, melibiose, sucrose, fructose, and D-mannose were significantly upregulated (more abundant) in bog samples (regardless of depth) whereas upregulated metabolites in fen samples were mainly classified as acids such as Erythronic acid, glycolic acid, pyruvic acid and amino acids such as glutamine and isoleucine. Another way to view this data is through the use of heatmaps. A heatmap of the same 25 identified metabolites (t-test significance) showed a similar separation (unpregulated sugar acids and simple sugars in bog and organic acid and amino acids in fen) between metabolites identified in different bog and fen samples ( Figure 5D). Figure 5C features the top four upregulated metabolites in bog and fen samples

Shifts in C Transformations With Thaw
Consistent with changes in SOM and primary metabolite composition, transformation analysis indicated that the biochemical processes associated with the SOM degradation were significantly different between the bog and fen sites ( Figure 6A). We observed higher abundances of sugar-associated transformations for the bog sites compared to the fen sites, which generally correlated with carbon transformations (Figure 6A). These data are consistent with Figures 2, 4B where bog sites had  a lower abundance of carbohydrate molecules with mass >400 Da but a higher abundance of mono and disaccharides, reflective of decomposition of complex carbohydrates into mono, and disaccharides. However, based on their higher relative abundance in the porewater (Figure 5), these sugars are likely being slowly utilized by microbial communities and/or undergoing slow decomposition. In contrast, fen samples had a higher abundance of complex carbohydrate molecules, but lower abundance of simple, mono, and disaccharides. Despite the higher abundance of complex carbohydrate compounds in the fen site and the higher abundance of sugar -associated transformations in the bog compared to the fen, principal component analysis on specific sugar -associated metabolic transformations on all elemental formulas detected in the bog and fen samples ( Figure 6B) indicated similar sugar transformations, suggesting that microbes at both sites have the capability of degrading simple and complex sugar molecules. These data suggest that simple sugars are being utilized at a higher rate in the fen sites compared to the bog sites. The slow decomposition at the bog sites and/or slow uptake of glucose by microbial communities for fermentation or anaerobic cellular respiration suggest additional controls (e.g., thermodynamic, abiotic processes) on SOM decomposition in the bog sites. The high amino acid-associated transformations in the fen samples ( Figure 6A) is consistent with a higher abundance of amino acid compounds at the fen sites, reflective of protein degradation and /or release of amino acids by fen vegetation through roots.  average number of neighbors for SOM molecules to porewater pH, DOC (mM) and δ 13 C -CH 4 values. SOM clustering coefficient increased with increasing pH (Figure 7A), decreasing DOC (Figure 7B), and enriched δ 13 C -CH 4 ( Figure 7C) values as in the case of fen and decreased with decreasing pH ( Figure 7A) and depleted δ 13 C -CH 4 values ( Figure 7C) as in the case of bog. On the other hand, SOM average number of neighbors had an opposite trend to the clustering coefficient (Supplementary Figure 4). Significant correlations (p < 0.05) were also observed between different biogeochemical properties including concentrations of DOC, TN, DIC (CO 2 ), CH 4, ammonia, nitrate, sulfate, δ 13 Cpw-DIC, δ 13 Cpw-CH 4 (Supplementary Table 4), revealing clear coupling between different biogeochemical properties and GHG emissions.

DISCUSSION
Based on our advanced analytical characterization of peat SOM along a permafrost thaw gradient, organic matter compositional differences were clear between the bog and fen sites. Fen SOM appeared to be more degraded (regardless of the depth), transformed, and reduced as evident by the more negative NOSC and higher G • C-ox values. This data is consistent with previous studies comparing bog and rich fen samples from different locations across the Stordalen Mire from different years (Hodgkins et al., , 2016. Phenol oxidase is an important enzyme in peatlands, as it is one of the few that can decompose recalcitrant phenolic compounds generally produced by Sphagnum spp. (McLatchey and Reddy, 1998;Freeman et al., 2004). Degradation of these phenolic compounds by microbes without the aid of phenol oxidase is challenging due to the high resonance energy that stabilizes the C-C bonds in aromatic rings (Harwood and Parales, 1996). Our data indicated that the bog SOM was more oxidized, dominated by tannin-like (polyphenolic) and organic acids, all with large number of carboxyl, aromatic, and methoxylated phenol groups (higher MWs). Previous reports have shown inhibition of microorganisms and the extracellular enzymes by these large polymers via crossbridged hydrogen bonds with proteins (enzymes) as well as penetration to the sensitive sites within microorganisms, hence altering their physiology (Wang et al., 1991;Field and Lettinga, 1992). These findings suggest that the high abundance of phenolic compounds in peat bogs could have a toxic effect on microorganisms in addition to lower decomposition rates by extracellular enzymes (Freeman et al., 2004). The activity of phenol oxidase is also dependent on various factors such as bimolecular oxygen availability, extent of waterlogging, and plant communities, which in the peat bog are factors controlled by Sphagnum spp. (Freeman et al., 2004). Sphagnum spp. creates waterlogged conditions, acidifies its environment using high cation exchange rates (Rudolph and Samland, 1985;Rasmussen et al., 1995;Thormann, 2006), produces lignin-like polymeric phenols that inhibit enzymatic activity (such as 5-keto-D-mannuronic acid, phydroxy-b[carboxymethyl]-cinnamic acid) (Stalheim et al., 2009;Fudyma et al., 2019), and has highly recalcitrant cell walls (Stalheim et al., 2009;Hájek et al., 2011). Sphagnum moss is also known to decay more slowly than vascular plants Tfaily et al., 2014). The preservation of labile compounds in the peat bog observed in our data is likely due to the anoxic, acidic, and nutrient sequestration traits intrinsic to Sphagnum spp. and likely impedes phenol oxidase activity and subsequent organic matter decay by microbial communities (Williams et al., 2000;Freeman et al., 2004). Diversity in microbial communities is important to the function and stability of many microbial ecosystems, such as animal-and plant-associated microbiomes (Worm et al., 2002;Muscarella et al., 2019). Uncovering mechanisms to maintain microbial diversity in these environments is important in understanding belowground metabolic transformations and subsequent GHG emissions (Muscarella et al., 2019). Substrate generalists vs. specialists can interact with different available substrates and depending upon the substrates that are used, we can partly understand microbial preferences and predict gas flux changes (Chen et al., 2020;Zhou et al., 2020). Based on our data, bog SOM exhibited higher DOC porewater concentrations, lower microbial diversity (inferred from Shannon diversity; Singleton et al., 2018;Woodcroft et al., 2018), and activity (inferred from the high average number of neighbors and lower clustering coefficient in the bog metabolic transformations analysis- Supplementary Table 3), suggesting that bog microbial communities are comprised of substrate specialists that are only capable of performing a limited set of specific functions (Monard et al., 2016;Muscarella et al., 2019) under anoxic and acidic conditions mainly created by Sphagnum spp. (Williams et al., 2000;Freeman et al., 2004). On the other hand, fen SOM was characterized by increased substrate heterogeneity (higher clustering coefficient), suggesting fen microbial communities are substrate generalists and are able to use and transform same substrates through different metabolic pathways. This observation is again consistent with Woodcroft et al. (2018) that showed microbial diversity (Shannon index) in the bog to be significantly lower than that in the fen. GC-MS data showed that the majority of upregulated metabolites in the bog were mono and disaccharides and sugar acids, in particular galacturonic acid, while the fen had higher abundances of amino acids and polyphenols (Figures 5, 6B), suggesting that in the fen, mono and disaccharides are consumed by microbial communities as soon as they are degraded by various biotic and abiotic metabolic pathways and the high abundances of amino acids and polyphenols reflect microbial byproducts and thus an active microbial community (Whiting and Chanton, 1992;Chanton et al., 1995Chanton et al., , 2008Hodgkins et al., 2014).
Recently, Fudyma et al. (2019) showed that Sphagum fallax metabolite profile is composed of flavonoid glycoside compounds (a compound formed from a simple sugar and a phenol), that could be acting as potent antimicrobial compounds. Thus, the quality of bog SOM (i.e., high abundance of lignin-like and polyphenolics and simple sugars, recalcitrant dead Sphagnum litter, and the low phenol oxidase enzymes activity) all suggest limited OM degradation by biotic factors. Here we hypothesize that the accumulation of simple sugars is either due to low rates of biotic decomposition and consumption, leading to its accumulation or that the release of simple sugars is occurring abiotically through acid catalyzed hydrolysis of lignocellulose and cellulose (BeMiller, 1967;Timell, 1967;Hopkinson, 1969;Girisuta et al., 2007;Ahmad and Pant, 2018) (Figures 2A, 5B-D). Therefore, we further postulate that low pH could be aiding the acid catalyzed hydrolysis of carbohydrates and sugar acids (glycosides) in bog systems (Ahmad and Pant, 2018) leading to the release of mono and disaccharides with lower average G • C-ox. However, accumulation of these compounds in the peat with time probably due to low pH and high abundance of galacturonic acids produced by Sphagnum, further shed light on the limited microbial activities in respiring these compounds (Figures 1, 3, 5) (Rasmussen et al., 1995;Thormann, 2006). Metabolic transformation analysis and FT-ICR-MS data confirm this trend, with higher abundance of sugar transformations present in bog samples. Because both reactant and product of the metabolic transformation are present in the bog samples, we believe in the bog biotic degradation of organic matter happens more slowly than at fen sites, and that possible byproducts of Sphagnum OM degradation in the bog (i.e., upregulated galacturonic acid-which is a by-product of pectin degradation known to have antimicrobial effects- Hájek et al., 2011) could be limiting the biotic complex carbohydrate degradation pathways, resulting in the accumulation of both complex and labile sugars. Observed correlations between the pH and the SOM clustering coefficient (positive correlation; p < 0.0001) and the average number of neighbors (negative correlation; p = 0.0002) potentially indicates the impact of pH on substrate quality and concentration. Lower pH values promote carbon accumulation and low substrate diversity whereas higher pH values are correlated with higher substrate quality and diversity. The bog porewater DOC concentration correlated positively (p = 0.0008) with the average number of neighbors, reflecting accumulation of OM in the bog sites consistent with lower decomposition rates (Supplementary Figure 3) and acidic conditions (low pH). Moreover, in fens, vascular plant roots may transport oxygen into the soil, where it activates phenol oxidase and thereby decreases the concentration of decomposition-inhibiting phenolics and organic acids (Joabsson and Christensen, 2001;Estop-Aragonés et al., 2012). However, further experiments are required to validate these assumptions.
Depleted δ 13 C -CH 4 in the bog porewater samples suggest that hydrogenotrophic rather than acetcolastic methanogenesis is likely the preferred pathway for methane production. This is consistent with low degradation potential of the bog SOM and higher DOC concentrations McCalley et al., 2014;Tfaily et al., 2014). More enriched δ 13 C -CH 4 in fen porewater, on the other hand, suggests that acetoclastic methanogenesis is the dominant pathway for CH 4 production and that methanogenisis substrates are produced by OM fermentation leading to low DOC concentration (Chanton et al., 1995;Chasar et al., 2000). Observed correlation between the clustering coefficients and the δ 13 C -CH 4 for the bog and fen samples potentially suggests that the network analysis parameters indicating the substrate diversity and heterogeneity can be used as an indication of pathways by which microbial communities are utilizing the OM and ultimately producing GHG and they indirectly provide information on the microbial diversity. However, further studies should aim to replicate these findings in a larger experimental setup.
Pathways underlying CO 2 and CH 4 productions in peatlands are also dependent on S and N cycling and the microorganisms that interact with these compounds (Granberg et al., 2001;Limpens et al., 2008;Pester et al., 2012). Sulfate (SO 4 2− ) and nitrate (NO − 3 ) can act as terminal electron acceptors in anoxic conditions, allowing microorganisms to couple the oxidization of substrates (e.g., H 2 , organosulfur and organonitrogen compounds) and reduction of TEAs to produce CO 2 and CH 4 (Muyzer and Stams, 2008). Hence sulfate reducing microorganisms (SRM) and methanogens can coexist or compete for substrates within these oxygen-limited environments (Raskin et al., 1996;Schink, 1997). Recently, Woodcroft et al. (2018) identified the presence and activation of high diversity deep-branching lineages of the dsrAB genes (dissimilatory sulfate reduction genes) belonging to sulfatereducing microorganisms (SRMs) in Sphagnum-dominated bog in the Mire. SRMs utilize sulfur containing compounds and produce CO 2. In peatlands, H 2 is generated from fermentation of carbohydrates (especially the ones with aromatic rings), while oxidized S-containing compounds come from plant matter (Rawat et al., 2012). Peat bogs generally have larger amounts of oxidized organosulfur compounds due to high S uptake potential by Sphagnum and therefore dead Sphagnum is rich in organic S compounds, which can drive sulfate reduction (Urban et al., 1989;Novák and Wieder, 1992;Fritz et al., 2014). Our data suggests that hydrogenotrophic reduction by SRM and nitrate reduction is occurring in the peat bog due to the following factors: (1) increased organosulfur compounds, potentially acting as TEAs for SRMs and suggesting that sulfate reduction is not limited by S availability (Baena et al., 1998); (2) low porewater mM concentrations of SO 4 2− and NO − 3 , suggesting rapid turnover of TEAs; (3) low pH, which is likely inhibiting methanogen activities (Tang et al., 2018); and (4) high abundance of mono and disaccharides, suggesting degradation of large carbohydrates by fermentative bacteria in addition to acid-catalyzed dissociations. Fermentation byproducts such as CO 2 , and H 2 , are then further consumed in hydrogenotrophic reduction (Conrad et al., 1987;Goodwin et al., 1991). CO 2 and CH 4 porewater concentrations measured in the bog samples were lower compared to the fen samples (Table 4), however, CO 2 :CH 4 production ratios (Table 4) were higher in the bog compared to the fen, suggesting likely consumption of fermentation and acid-catalyzed dissociation byproducts by SRMs, leading to suppressed methanogenesis potentials by methanogens in the bog (Ward et al., 2009;Rawat et al., 2012;Hausmann et al., 2018). Similarly, SRMs can oxidize CH 4 coupled with sulfate reduction, which is also likely driving the lower CH 4 production rates in bogs (Cord-Ruwisch and Ollivier, 1986;Ward et al., 2009;Rawat et al., 2012). Higher CO 2 and CH 4 porewater concentrations in the fen are likely due to a decrease in Sphagnum-derived inhibitory organic acids (Figure 5) driven by changing plant inputs and increasing pH ( Table 1 and Supplementary Table 1) from bog to fen systems. Lower C/N ratios in the fen, driven by lower C/N of plant inputs from sedges compared to Sphagnum, may also lead to greater organic matter decomposition in the fen by alleviating N limitation for decomposers Khadka et al., 2016).

CONCLUSION
Due to the contribution from aboveground vegetation and belowground metabolic transformations under various redox conditions, microbial processes both affect and are being affected by the chemical composition and quality of SOM. In order to gain better insight on the overall metabolic transformations in bog and fen and to link substrate diversity and quality with microbial activity and GHG production, we calculated the Gibbs free energy of C half reactions based on C oxidation state for OM, a quantifiable measure for OM decomposability and quality, which can be incorporated into current biogeochemical models to represent OM lability and microbial diversity. We further correlated clustering coefficients and average number of neighbors, two network analysis parameters, with different biogeochemical properties of bog and fen samples. Our data revealed that pH is a main driver of SOM decomposition. We hypothesized that the low pH of the bog, driven by Sphagnum-derived organic acids, along with phenolics, could potentially impact substrate quality and concentration and microbial diversity thereby promoting carbon accumulation and low substrate diversity. We further hypothesized that abiotic degradations such as acid catalyzed hydrolysis reactions (i.e., hydrolysis of Lignocellulose, complex carbohydrates, and sugar acids) are happening at a greater rate in bog, freeing bioavailable compounds for microbial utilization and further biotic metabolic pathways. In the fen, the presence of lower molecular weight compounds with higher Gibbs free energy, low DOC, more populated and diverse metabolic transformations, and higher net CH 4 and CO 2 fluxes and higher pH is indicative of more diverse substrate composition and higher microbial diversity than the bog, suggesting a more active system where biotic pathways are largely upregulated. The correlation observed between the clustering coefficient and the average number of neighbors with δ 13 C -CH 4 values further reflects that organic matter quality and composition has an impact on microbial activity (as revealed by FT-ICR-MS and GC-MS data) measured indirectly through CH 4 fluxes.
While many microbial pathways can be inhibited by acidic and anoxic conditions, as observed in the bog, we found evidence of activity of sulfur reducing bacteria. Lower concentrations of sulfate in the bog porewater suggest a rapid turnover of organosulfur compounds acting as TEAs. This is also confirmed by lower CH 4 production rates and lower pH in the bog, with low pH likely inhibiting methanogenesis, as well as SRM either consuming H 2 and producing CO 2 or consuming CH 4 and producing CO 2 (CH 4 + SO 2− 4 → HCO − 3 + HS − + H 2 O). Due to these results, fen ecosystems that are composed of more diverse and active microbial communities should be more stable with changing environmental conditions (e.g., similar GHG fluxes) whereas bog microbial populations may be more easily affected by a changing climate. It is expected that thaw progression in northern peatlands will rapidly change the water table level, plant communities and nutrient input sources, all contributing to a decline of Sphagnum moss (Norby et al., 2019). This could result in shifts in decomposition, as well as altered contributions from biotic and abiotic pathways of OC toward higher CH 4 and CO 2 emission which ultimately creates positive climate feedbacks. Therefore, insights into how metabolite diversity and resource heterogeneity relate to the taxonomic and functional diversity of microbial communities allows for aspects of ecosystem stability and functional redundancy to be studied and provides a unique opportunity to understand how microbial communities will respond to climate perturbations.

DATA AVAILABILITY STATEMENT
The original contributions generated for the study are included in the article/Supplementary Materials. The data sets generated from this study are currently available online at https://doi.org/ 10.17605/OSF.IO/GS2UT.

AUTHOR CONTRIBUTIONS
RW, SH, and MT designed and performed the experiment and data acquisition. RA and MT did the analysis and interpretation of data for the work. RA, MT, and RW drafted and revised the manuscript with contribution from SH, JF, JC, VR, and SS. All authors contributed to the article and approved the submitted version.