Abstract
Glacier foreland soils have long been considered as methane (CH4) sinks. However, they are flooded by glacial meltwater annually during the glacier melting season, altering their redox potential. The impacts of this annual flooding on CH4 emission dynamics and methane-cycling microorganisms are not well understood. Herein, we measured in situ methane flux in glacier foreland soils during the pre-melting and melting seasons on the Tibetan Plateau. In addition, high-throughput sequencing and qPCR were used to investigate the diversity, taxonomic composition, and the abundance of methanogenic archaea and methanotrophic bacteria. Our results showed that the methane flux ranged from −10.11 to 4.81 μg·m−2·h−1 in the pre-melting season, and increased to 7.48–22.57 μg·m−2·h−1 in the melting season. This indicates that glacier foreland soils change from a methane sink to a methane source under the impact of glacial meltwater. The extent of methane flux depends on methane production and oxidation conducted by methanogens and methanotrophs. Among all the environmental factors, pH (but not moisture) is dominant for methanogens, while both pH and moisture are not that strong for methanotrophs. The dominant methanotrophs were Methylobacter and Methylocystis, whereas the methanogens were dominated by methylotrophic Methanomassiliicoccales and hydrogenotrophic Methanomicrobiales. Their distributions were also affected by microtopography and environmental factor differences. This study reveals an alternative role of glacier foreland meadow soils as both methane sink and source, which is regulated by the annual glacial melt. This suggests enhanced glacial retreat may positively feedback global warming by increasing methane emission in glacier foreland soils in the context of climate change.
Introduction
Methane (CH4) is the second most important greenhouse gas in the atmosphere, with 28 times the global warming potential of CO2 on the centennial-scale (IPCC, 2013). Methanogens and methanotrophs are responsible for methane-cycling in the natural environment, and their relative activities determine the global methane dynamics (Le Mer and Roger, 2001), while anaerobic methanotrophs, coupling methane oxidation to the reduction of nitrate, nitrite, iron, manganese, and sulfate, also play a role in the methane cycle (Knittel and Boetius, 2009; Ettwig et al., 2010, 2016; Haroon et al., 2013; Hu et al., 2014; Leu et al., 2020; Guerrero-Cruz et al., 2021). Methane flux in glacier foreland soils has attracted increased attention due to the rapid retreat of mountain glaciers under global warming (Bárcena et al., 2010). Methane dynamics in the natural environment depend on nutrient availability (Le Mer and Roger, 2001), vegetation coverage (Adachi et al., 2006), terrain topography (Wei et al., 2015b), soil moisture content (Coles and Yavitt, 2002; Wei et al., 2015a), and temperatures (Avery et al., 2003; Metje and Frenzel, 2005; Bárcena et al., 2010). Global warming has greatly enhanced glacier retreat, which accelerates the expansion of glacier foreland. Thus, the methane dynamics at glacier foreland soils can enhance or mitigate the impact of climate change.
Glacier foreland soils are typically regarded as a sink for methane, and most previous studies on methane-cycling microorganisms have focused on the community succession of methanotrophs (Nauer et al., 2012; Chiri et al., 2015, 2017; Mateos-Rivera et al., 2018). For example, a previous study measured the abundance of methanotrophs and revealed a stable CH4 uptake (−0.082 to −2.2 mg CH4 m−2 d−1) during the snow-free season in Damma and Griessfirm glacier forefields in Switzerland (Chiri et al., 2015). In contrast, a net CH4 production has been reported for some glacier foreland (e.g., the Swiss Alps glacier foreland; Bárcena et al., 2010; Nauer et al., 2012). These results suggest that glacier foreland soil may also serve as a potential role of methane sources. However, the driving factors that regulate the methane flux in glacier foreland are largely unknown.
Accelerated glacier retreat has occurred on the Tibetan Plateau since the early 20th Century (Yao et al., 2007), and the exposed barren soil is expected to develop into meadow within 80 years (Eichel, 2019). The alpine meadow soils on the Tibetan Plateau are considered an important sink for atmospheric CH4 (Wei et al., 2015b). In comparison, wetland and swamp meadows are major methane source regions, which could emit 0.56 ~ 1 Tg CH4 a−1 (Jin et al., 1999; Ding and Cai, 2007). Alpine meadow and swamp meadow are inter-convertible by the alternation of hydrological conditions. The enhanced glacial retreat can increase glacial meltwater discharge, which affects the hydrological processes in the glacier foreland regions (Yao et al., 2007; Bradley et al., 2014). During the melting seasons, a large amount of glacial meltwater flows into glacier foreland soils, transforming alpine meadow into swamp meadow, and reducing the oxygen availability, which favors methanogens over methanotrophs (Chen et al., 2008). This may transform the glacier foreland soils from a methane sink to a methane source, but the microbial process underlying this transformation is poorly understood.
The overall goals of this study were to: assess the impact of glacier melting on the dynamics of methane flux and identify the environmental and biological drivers of the methane dynamics. Based on the impact of glacier meltwater on the hydrological conditions in glacier foreland soils, we hypothesized that the methanogens would be more abundant than methanotrophs, and turn glacier foreland meadow into a methane source.
Materials and Methods
Foreland of the Longxiazailongba Glacier
The Longxiazailongba (LXZ) Glacier is located in the middle of the Tanggula Mountain on the Tibetan Plateau, China (Figure 1). The glacier is 7 km in length, 19.3 km2 in area, 118 ± 10 m in thickness, and the highest altitude is 6,000 m above sea level (a.s.l.; Yao, 2014). The altitude of the glacier terminus is about 5,240 m a.s.l. The area near the terminus (about 0.5 km in length and 1.1 km in width) was moraine without any vegetation. Vegetation was developed at approximately 600 m from the glacier terminus, and meadows (dominated by Kobresia tibetica) were developed at approximately 1.5 km from the glacier terminus (Figure 1A).
Figure 1
The glacial meltwater enters the foreland from Mid-June, transforming the alpine meadow into a swamp meadow until October. Therefore, June is defined as pre-melting season, and July to October is defined as melting season in the present research. During the glacial melting season, the swamp meadow exhibits both hummock and hollow microtopography (Figure 1B), which are commonly observed in wetlands. Hummocks are typically located higher above the water table, and hollows are closer to the water table and may occasionally be flooded during glacial melt season (Belyea and Clymo, 1998). At the LXZ glacier foreland, hummock soils were covered by K. tibetica, while hollows were inundated with glacier meltwater and no vegetation was observed (Figure 1B).
Soil Sampling
In August 2017, soil samples were collected at three sites at the glacier foreland during the melting season (Figure 1A). These meadows feature both hummock and hollow microtopography. At each site, soil samples were collected from three hollow or hummock microtopography as three independent replicates. In each microtopography, five surface soil samples (0–10 cm) were collected randomly using a sterile shovel and were then mixed thoroughly to form a composite hollow or hummock sample. The fresh soil samples were stored in Whirl-Pak bags at approximately 4°C in a portable refrigerator and delivered to the laboratory within 48 h. In the laboratory, soil samples were frozen at −80°C for further analysis.
Soil Characteristics Measurement
Moisture was measured by using the gravimetric method. Approximately 10 g 2 mm sieved fresh soil was weighted, and oven-dried at 105°C until no further mass loss was observed and then reweighted. The moisture content is expressed as the mass of water per mass of dry soil. The pH values were measured after mixing wet soil with distilled water at the soil-to-water ratio of 1:5 (g/g). Soil organic matter (OM) was determined by the external heating-potassium dichromate volumetric method with air-dried soils (Ji, 2005). The total nitrogen content (TN) was determined using the Kjeldahl method (Bremner, 1960). Soil ammonium nitrogen (NH+ 4) and nitrate-nitrogen (NO-3) were extracted from wet soils with 2 M KCl (soil/solution, 1:5) using Smartchem200 Discrete Auto Analyzer (Alliance, France). The available phosphorus contents (AP) were determined using the acid digestion method (Kuo, 1996).
DNA Extraction and Quantitative PCR
Genomic DNA was extracted from 0.5 g soil using Fast DNA®SPIN Kit for Soil (MP Biomedicals, Santa, CA, United States). The Quality and quantity of the extracted DNA were measured using a NanoDrop 2000 Spectrophotometer (Thermo-Scientific). In order to minimize the potential inhibitory effects of co-extracted substrates from soil (e.g., humic acid), soil DNA was diluted 10 times for qPCR. The mcrA gene encodes the α-subunit of methanogenic methyl coenzyme M reductase, and the pmoA gene encodes the α-subunit of methane monooxygenase, which is used for the initial conversion of CH4 to methanol. The copy numbers of mcrA and pmoA genes were quantified using qPCR with primer mals-mod-F (5′-GGYGGTGTMGGDTTCACMCARTA-3′)/mcrA-rev-R (5′-CGTTCATBGCGTAGTTVGGRTAGT-3′; Angel et al., 2012) and A189f (5′-GGNGACTGGGACTTCTGG-3′)/mb661r (5′-CCGGMGCAACGTCYTTACC-3′; Kolb et al., 2003), respectively. Quantitative PCR amplification was performed on a CFX-96 Optical Real-Time PCR System (Bio-Rad Inc. Hercules, CA, United States). The 20 μl reaction mixture contained 2 μl of the DNA template, 10 μl of TB GreenTM Premix EX TaqTM II (TaKaRa), 0.4 μl of each forward and reverse primer (20 μM), and 7.2 μl ddH2O. qPCR conditions for mcrA is as follows: initial denaturation (94°C, 6 min), followed by 50 cycles of denaturation (94°C, 25 s), annealing (65.5°C, 20 s), and elongation (72°C, 45 s). While qPCR condition for pmoA is as follows: initial denaturation at 95°C for 30 s, followed by 40 cycles of 5 s at 95°C, 30s at 54°C, 40s at 72°C, and 30s at 80°C. A melting curve analysis was conducted to confirm the specificity of the PCR products.
In situ CH4 Flux Measurement
Methane flux was measured in situ in June (pre-melting season) and August (melting season) 2020. Methane flux was measured using polyvinylchloride cylindrical (PVC) tubes at site 1 to site 3 (Figure 1). Tubes consist of two parts, a PVC base (25 cm in diameter and 10 cm in height) installed permanently into the soil of about 5 cm, and a cylindrical box (without bottom, 25 cm in diameter and 30 cm in height; Figure 1B). In the field, six chambers were placed simultaneously, with three chambers placed on hollow soils and the other three on hummock soils. Chambers were closed for 20 min before the gas was collected. The gas was collected every 15 min using a plastic syringe (100 ml) and stored in a sterile gas package (Haide, Dalian).
The gas samples were analyzed within 48 h using a gas chromatograph (Agilent GC-7890B, United States) at the Naqu Ecological and Environmental Observation and Research Station, China (31°17’N, 92°06′E; 4,501 m a.s.l.). The gas chromatograph was equipped with a flame ionization detector (FID) and electron capture detector (ECD), using N2 as the carrier gas to remove O2 and water vapor (Chen et al., 2017).
The CH4 flux was calculated as the following (Wang and Wang, 2003):
Where, F is the methane flux (mg·m−2·h−1), M is the molecular mass of methane (16.12 g/mol), V0 is the gas molecular mass (22.41 L/mol), P is the atmospheric pressure at the sampling site, P0 is standard atmospheric pressure (1013.25 mbar), is the slope of the linear regression for the methane concentration gradient through time (m3·m−3·h−1), and H is the chamber height above the soil (m).
High-Throughput Sequencing and Data Analysis
The compositions of methanogenic and methanotrophic communities were evaluated using high-throughput sequencing targeting the mcrA and pmoA genes, respectively. The primer pairs for mcrA and pmoA genes were mlas-mod-F/mcrA-rev-R and A189f/mb661r, respectively (Costello and Lidstrom, 1999; Angel et al., 2012). PCR amplification was performed with primers with sample-specific barcodes. Around 35 cycles were used to amplify both the mcrA gene and pmoA gene PCR products. The PCR products of each sample were purified and quantified using a Qubit instrument (Life Technologies). The PCR products were pooled in an equimolar concentration, and sequenced on an Illumina MiSeq system using 2 × 300 cycle combination mode at Shanghai Meiji Biotechnology Co., Ltd.
For the mcrA and pmoA gene sequence analyses, the paired-end sequences were merged and quality checked. Candidate OTUs sequences of mcrA and pmoA genes were obtained using the “unoise3” command in Usearch v11.0.667 (Edgar, 2013). Chimeras were removed during the clustering process. These candidate OTUs sequences were uploaded to FrameBot (Wang et al., 2013) to check the frameshift and delete frameshift sequences. In addition, the candidate OTUs sequences were imported into ARB (Ludwig et al., 2004) to calculate the distance matrix based on amino acid sequences. Using cluster command in Mothur (Schloss et al., 2009), the final mcrA and pmoA OTUs at approximate species-level were assigned with amino acid dissimilarity levels of 0.11 (Yang et al., 2014) and 0.07 (Lüke and Frenzel, 2011), respectively. Meantime, we double-checked the number of new mcrA OTUs based on the distance level of 0.11 amino acids sequences, which was equal to the number obtained based on 0.84 nucleotide distance (Yang et al., 2014). Finally, based on the representative sequences obtained above, the OTU table was generated using the “otutab” command in Usearch.
To establish the phylogenetic trees for methanogens and methanotrophs, representative OTUs sequences and the reference sequences were first aligned using mafft v7.464 (Rozewicki et al., 2019). AliView (Larsson, 2014) was used to check and realign the aligned sequences based on the length of sequences that automatically trimmed by trimAI (Capella-Gutiérrez et al., 2009). FastTree v2.1 (Price et al., 2010) was used to generate the approximate maximum likelihood tree. Finally, these tree files were visualized in iTOL (Letunic and Bork, 2016).1 Meanwhile, the taxonomic information of the representative OTUs was assigned according to the relationship between the representative OTUs and the reference sequences.
Statistical Analyses
The Kruskal–Wallis test was used to identify the significant differences in soil physicochemical properties, methane-cycling microorganisms, and methane fluxes among different microtopography using the software of PAST v3.0 (Hammer et al., 2001). Alpha diversity of Chao1 richness and Shannon indexes were processed with vegan package v2.5-7 (Dixon, 2003). Spearman correlations were used to determine the relationships of methane flux and methane cycling communities with soil properties. The heatmap was performed with the package of “heatmap” (Kolde, 2013) in R 4.0.4.
Results
Soil Physiochemical Properties
The physicochemical properties of glacier foreland soils are shown in Supplementary Table S1. All soils were neutral to slight alkaline with pH values ranging from 6.94 to 8.14. The air temperature was 10.43 ± 5.06°C and 10.96 ± 4.71°C when methane flux was measured in June and August, and no significant difference was detected (Kruskal–Wallis test, p = 0.282). The pH values of hollow soils were significantly higher than those of hummock soils (Kruskal–Wallis test, p < 0.001). In contrast, the available phosphorous (AP) was significantly higher in hummock soils than in hollow soils (p = 0.013). Other soil properties, including soil water content, values of total nitrogen (TN), ammonium (NH+ 4), nitrate (NO-3), and organic matter (OM) had no significant differences between hollow and hummock samples.
The Abundances of Methane-Cycling Microorganisms and Methane Flux
The abundance of methanogens was inferred by the mcrA gene copy number, which ranged from 1.19 × 106 to 5.17 × 108 copies g−1 dry soil. The abundance of methanotrophs was inferred using the pmoA gene copy number, which ranged from 4.76 × 105 to 1.09 × 108 copies g−1 dry soil (Supplementary Table S2). Except for the hummock soils in site 2, the mcrA-to-pmoA gene copy number ratios were consistently greater than 1 (z-test, all p < 0.05; Figure 2A). Correlation analysis indicated that both mcrA and pmoA gene abundance were significantly and positively correlated to soil pH (Figure 3). In addition, pmoA gene abundance was also significantly and positively correlated to soil water content, while no significant correlation was observed between soil water content and mcrA gene abundance (Figure 3).
Figure 2
Figure 3
To evaluate the impact of glacier meltwater on methane dynamics in glacier foreland soils, methane fluxes during the pre-melting season (June) and glacier melting season (August) were measured. In the pre-melting season, the methane flux was −3.76 ± 9.84 μg·m−2·h−1, and significantly increased to 13.95 ± 12.44 μg·m−2·h−1 in the glacier melting season (Kruskal–Wallis test, p < 0.001; Figure 2B). Methane flux variations were also observed by the microtopographic effects, the methane flux was −8.29 ± 9.44 μg·m−2·h−1 in hollow soils in the pre-melting season, which was significantly higher than that in hummock soils (0.78 ± 7.97 μg·m−2·h−1, p = 0.027). However, such variation diminished during the glacier melting season, i.e., the methane flux in hollow soils was 14.27 ± 12.15 μg·m−2·h−1, which was no longer significantly different from that in hummock soils (13.64 ± 12.71 μg·m−2·h−1, p = 0.964). There is a negative and statistically significant correlation between the methane flux and soil pH (r = −0.48, p = 0.003; Figure 2C). However, the correlation between methane flux and pH was insignificant when the methane flux in the pre-melting season and melting season were investigated separately (p = 0.08 and 0.75, respectively). Furthermore, the methane flux significantly correlated with soil moisture in the pre-melting season (p = 0.02), but not in the melting season (p = 0.17; Figure 2D).
The Diversity of Methane-Cycling Microorganisms in Glacier Foreland Soils
For methanogens, the Chao1 richness and Shannon indexes ranged from 5.0 to 13.0, and 1.3 to 3.1, respectively. Both indexes were significantly higher in hollow soil samples than those in hummock soil samples (Kruskal–Wallis test, all p < 0.001; Supplementary Figures S1A,C). The Chao1 richness and Shannon indexes of the methanotrophs ranged from 4 to 11, and 0.58 to 2.27, respectively, and no significant differences in both alpha diversity indexes were observed between the hollow and hummock soils (p = 0.19 and p = 0.44, respectively; Supplementary Figures S1B,D).
Spearman correlation analyses demonstrated that both Chao1 richness and Shannon indexes of the methanogens were significantly and positively correlated with soil pH (r = 0.82, p < 0.001; Figure 3). The Chao1 richness index of the methanotrophs was positively correlated with pH (r = 0.58, p = 0.02), and the Shannon index was positively correlated with OM (r = 0.54, p = 0.02) and TN (r = 0.51, p = 0.03).
Methanogenic and Methanotrophic Community Compositions
In the present study, a total of 15 mcrA genes OTUs were obtained and classified. Within these OTUs, three were affiliated with the Methanomassiliicoccales (Figure 4A), which together accounted for about 50% of the sequences retained (Figure 4C). Four additional OTUs were affiliated with Methanomicrobiales and accounted for another 25%. The other OTUs were affiliated with the Methanobacteriales (three OTUs), Methanosarcinales (two OTUs), Methanotrichales (two OTUs), and Methanocellales (one OTUs), which accounted for the remaining 25% of the community. No significant differences were identified between hollow and hummock soils in Methanobacteriales, Methanocellales, Methanomicrobiales, and Methanotrichales (Kruskal–Wallis test, p = 0.413, 0.269, 0.075, and 0.962, respectively; Supplementary Figure S2A). The relative abundances of Methanomassiliicoccales and Methanosarcinales were found to be more abundant in hummock soil samples than in hollow soil samples (p = 0.005 and p < 0.001, respectively). Furthermore, the relative abundance of three OTUs (OTUs 1, 3, and 9) was significantly higher in the hummock soils, and their relative abundance all exhibited negative correlations with pH (Figure 5A). In comparison, the relative abundance of eight OTUs was significantly higher in the hollow soils (OTUs 4, 6, 7, 8, 11, 12, 14, and 15), and their relative abundance all exhibited positive correlations with pH value.
Figure 4
Figure 5
Taxonomic analysis of pmoA genes revealed the dominance of type I methanotrophs (Figure 4B). Type Ia methanotroph Methylobacter was the most abundant, accounting for approximately 42% of the sequences retained (Figure 4D). Other type Ia (LP20, Methylomonas and Aquifer-cluster) and type Ib group (FWs) were also identified but with much lower abundance. In addition, type II methanotrophs (Methylocystis) were detected with an average relative abundance of 21.6%. Hollow soil samples exhibited a significantly higher relative abundance of Methylobacter (Kruskal–Wallis test, p = 0.010) and Methylomonas (p < 0.001) than that in hummock soils, while a higher relative abundance of Methylocystis (p = 0.013) and FWs (p = 0.017) were detected in hummock soils (Supplementary Figure S2B). The relative abundance of Methylobacter significantly correlated with soil water content and organic matter negatively (p = 0.012 and p = 0.040, respectively), while the relative abundance of Methylomonas significantly correlated with pH and available phosphorous positively and negatively, respectively (p < 0.001 and p = 0.005, respectively; Figure 5B). The relative abundance of Methylocystis positively correlated with soil moisture (p = 0.008), while FWs negatively correlated with pH (p = 0.036) and positively correlated with available phosphorous (p = 0.005).
Discussion
Glacier Foreland Soil as an Atmospheric CH4 Source During the Melting Season
The copy number of the mcrA gene was typically higher than that of the pmoA gene in the glacier melting season (Figure 2A). Although the qPCR result does not reflect the actual number and activity of microorganisms, the mcrA-to-pmoA gene copy number ratio has commonly been used to estimate the relative abundance between methanogens and methanotrophs (Morris et al., 2016; Yuan et al., 2018; Kong et al., 2019). Thus, our results suggest that the methanogens were dominant over methanotrophs in the glacier melting season, and the LXZ glacier foreland soils could be a methane source. This was confirmed by the positive in situ methane flux, which showed that glacial meltwater turned glacier foreland meadow soils from a methane sink to a methane source (Figure 2B).
The LXZ glacier foreland meadow was a weak methane sink during the pre-melting season (−3.76 ± 9.84 μg·m−2·h−1; Figure 2B). However, its methane oxidation capacity was much smaller compared with other alpine meadows, which typically ranged from 28 to 71.5 μg·m−2·h−1 (Lin et al., 2009; Wei et al., 2015b; Wu et al., 2020). Glacier meltwater turned glacier foreland soils into a methane source. This also differed from other glacier foreland soils, which are predominately identified as atmospheric methane sinks (or as weak methane sources) during the glacier melting season, such as observed in Greenland (−0.76 ~ −0.14 μg·m−2·h−1; Bárcena et al., 2010), Switzerland (Damma and Griessfirn; −2,200 ~ −82 μg·m−2·h−1; Chiri et al., 2015, 2017), and Svalbard (−110 ~ 300 μg·m−2·h−1; Adachi et al., 2006). These contrastive differences could be attributed to the soil water content difference, which is a major driver of methane flux (Le Mer and Roger, 2001; Xu et al., 2010). The soil water content of the glacier foreland soils in the Arctic and the Alps (Chiri et al., 2015; Hofmann et al., 2016) is typically less than 30% even during the melting season. This is much lower compared with the soil water content in the present study (>72%), which could be due to the rapid glacier melting on the Tibetan Plateau (Yao et al., 2004; Qiu, 2008). Thus, the enhanced glacial meltwater discharge increased soil water content, which subsequently changed the redox potential and favored methanogens over methanotrophs (Hofmann et al., 2013).
Our results revealed a hidden methane source associated with annual glacial melting on the Tibetan Plateau. This seems to be unique compared with the glacier foreland soils in the Arctic and the Alps (Adachi et al., 2006; Bárcena et al., 2010; Chiri et al., 2017). This uniqueness could be associated with the faster glacier mass loss and greater glacial meltwater discharge (Adachi et al., 2006; Chiri et al., 2015, 2017), as the glacier mass change rate of the TP is faster than those in the Alps and Svalbard regions (Wouters et al., 2019). The Tibetan Plateau has the largest number of mid-latitude glaciers, thus the impact of glacier melting on the glacier foreland methane emission could be overlooked, which needs to be carefully evaluated in future studies.
Abundance and Diversity of Methanogenic and Methanotrophic Communities
The copy numbers of marker genes (mcrA) for methanogens in LXZ glacier foreland soils ranged from 1.19 × 106 to 5.17 × 108 copies g−1 dry soil, which is higher than those in the foreland soils of the Alps (4.6 × 104–2.5 × 106 copies g−1 dry soil; Hofmann et al., 2016) and the meadow soils of the Tibetan Plateau (6.0 × 105–6.7 × 106 copies g−1 dry soils; Yang et al., 2017). The mcrA gene copy number increased with elevated soil water content (Figure 3), which is consistent with the results in rice paddy fields (Ma et al., 2012) and Arctic soils (Høj et al., 2006). This finding indicates that the increase in soil water content can enhance the abundance of methanogens, and subsequently increase methane production.
Our results showed that the methanogens in LXZ glacier foreland soils were dominated by Methanomassiliicoccales (49%) and Methanomicrobiales (24%), which are consistent with those recovered from the NamCo wetland (Deng et al., 2019), Zoige wetland (Zhang et al., 2008), and Qilian Mountain alpine permafrost soils (Wang et al., 2020). Methanomassiliicoccales has been identified from a wide range of environments, such as tropical peat swamp forest soils (Too et al., 2018), lake sediment (Emerson et al., 2020), and alpine cave soils (Jurado et al., 2020). Methanomassiliicoccales is phylogenetically distant from other methanogen orders, and belongs to a large evolutionary branch composed of many non-methanogenic archaea (Borrel et al., 2014). The wide distribution of these archaea suggests that they could be adapted to psychrophilic/mesophilic conditions, and their methanogenic activity in LXZ glacier foreland soils may be further enhanced under global warming. In comparison, the other dominant hydrogenotrophic methanogen Methanomicrobiales were also dominant in the foreland soils of the Alps, and have also been identified from Antarctica and Alaska marine sediment (Dong and Chen, 2012), indicating their adaptation to the psychrophilic environment (Hofmann et al., 2016).
The pmoA gene abundance varied from 4.76 × 105 to 1.09 × 108 copies g−1 dry soil, which is higher than the siliceous and calcareous foreland soils in Switzerland (8.2 × 102 to 5.5 × 105 copies g−1 soil; Chiri et al., 2015, 2017), foreland soils in Norway (~103 g−1 soil; Mateos-Rivera et al., 2018), and foreland soils in the Central Alps (~103–105 g−1 dry soil; Hofmann et al., 2016). Type I methanotrophs Methylobacter and type II methanotrophs Methylocystis dominated the LXZ glacier foreland soils. These lineages have been identified in the wetland soils of the Tibetan Plateau (Deng et al., 2014; Zhang et al., 2019), foreland soils of Griessfirn Glacier (Chiri et al., 2017), and Swiss foreland soils (Nauer et al., 2012). In contrary to our results, USCγ cluster is the dominant type I methanotrophs in upland soil ecosystem, such as those reported in glacier foreland upland soils (Nauer et al., 2012; Chiri et al., 2017), Tibetan upland grassland soils (Deng et al., 2019), and Canadian upland tundra soils (Martineau et al., 2014). Members of the USCγ cluster exhibit high methane uptake affinity, and are capable of atmospheric methane scavenging (Knief et al., 2003). The observed difference in the USCγ cluster could be ascribed to the different environmental conditions, such as soil water content. In general, USCγ is more abundant in regions with low precipitation, whereas low-affinity methanotrophs (such as Methylobacter) dominate soils with high water content (Deng et al., 2019). The high-water content favors anaerobic methanogens, subsequently leading to enhanced methane production (Huttunen et al., 2003), which supports the growth of low-affinity methanotrophs (such as Methylobacter and Methylocystis).
Methane Flux and Methane-Cycling Microorganisms Within Microtopography
Our results illustrated that the hollow soils were a sink for atmospheric methane in the pre-melting season, while hummock soils were a weak source (Figure 2B). This could be due to the presence of vegetation in hummocks (Laine et al., 2007), as root respiration can enhance methanogenic activity and methane emissions by creating an anaerobic environment and providing nutrients via root extrudes (Sutton-Grier and Megonigal, 2011; Minick et al., 2019). Glacial meltwater discharge during glacier melting season made both hollow and hummock soils into methane sources (Figure 2B). This is evidenced by the high soil water content (Figure 2D), which can alter the redox potential by reducing oxygen availability (Wadham et al., 2007; Hofmann et al., 2013).
Our results further demonstrated the distinct compositions of methane-cycling microorganisms in glacier foreland hollow and hummock soils. Hummock soils enriched methanogens that negatively correlated with soil pH, whereas hollow soils enriched methanogens that positively correlated with pH. Our results also illustrated pH was the driving factor that regulates methanogens in glacier foreland hollow and hummock soils (Figure 5A). Hummock soils exhibited lower pH values than the hollow soils, thus suggesting hummock soil-dwelling methanogens could have a lower pH preference than those in hollow soils. This is consistent with a previous study that pH is the primary factor influencing methanogen distribution (Deng et al., 2014). Further to this, certain OTUs within the same order (i.e., Methanobacteriales and Methanomassiliicoccales) exhibited distinct pH preferences (Figure 5A). This finding may indicate the distinct adaptation strategies of the methanogens in the same order. In contrast, Methanomicrobiales, Methanosarcinales, and Methanotrichales were only enriched in hollow soils. Both Methanomicrobiales and Methanosarcinales are versatile regarding carbon sources compared with those more specialized species (such as the Methanobacteriales and Methanomassiliicoccales; Angelidaki et al., 2011). Methanomicrobiales and Methanosarcinales (Conrad, 2020) are also known to be oxygen tolerant, attributed to the presence of superoxide dismutase, which can protect them from oxygen toxicity (Guerrero-Cruz et al., 2018). This can be particularly important for methanogens in hollow soils, where the soil water content can be much lower and the soils can become aerated during the pre- and post-melting season.
The relative abundance of methanotrophs Methylobacter (Ia) and Methylomonas (Ia) was higher in hollow soils, while Methylocystis (IIa) and FWs (Ib) dominant in hummock soils (Figure 5B). This finding is supported by previous observations derived from arctic peat soils (Tveit et al., 2014), lake sediment (Rahalkar et al., 2007), and peat bog soils (Danilova et al., 2016). Methylobacter and Methylomonas have been found to present in habitats with very low oxygen concentration. During melting season, glacier foreland soils are flooded by the meltwater, providing a very low oxygen condition that allows the more appropriately adapted Methylobacter and Methylomonas taxa to grow, dominating in hollow soils. In comparison, Type II methanotrophs such as Methylocystis have frequently been identified in plant rhizosphere (Bender and Conrad, 1995), partially due to their ability to fix nitrogen (Mäkipää et al., 2018). Furthermore, Methylocystis could grow on multi-carbon compounds such as acetate (Hakobyan and Liesack, 2020), which is available through root exudates (Girkin et al., 2018).
Conclusion
This study shows that meadow soils in the glacier forefield of a receding glacier are an overlooked methane sink in the non-melting season and can be transferred into a methane source in the melting season, and this transformation is regulated by glacier melting and also influenced by local vegetation coverage. Furthermore, the methane production and oxidation balance can be further changed under the projected global warming and accelerate glacier retreat.
Funding
This research was funded by National Key R&D Program of China (Grant No. 2019YFC1509103), the National Natural Science Foundation of China (Grant Nos. 41971077 and 91851207), the Second Tibetan Plateau Scientific Expedition and Research (STEP) program (Grant No. 2019QZKK0503), and the Strategic Priority Research Program (A) of the Chinese Academy of Sciences (Grant No. XDA20050101).
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented for this study can be found in the NCBI Sequence Read Archive under study accession number PRJNA679195. The soil properties and methane fluxes in present study have been archived at the PANGAEA dataset (https://issues.pangaea.de/browse/PDI-31033 and https://issues.pangaea.de/browse/PDI-31032).
Author contributions
TX, PL, YD, and MJ performed the experiments, analyzed the data, and prepared figures and tables. TX and YD drafted the work. YL, MJ, YD, KL, and PL conceived and designed the experiments and approved the final draft. WW and YL contributed to field work. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interests.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.862242/full#supplementary-material
Footnotes
References
1
AdachiM.OhtsukaT.NakatsuboT.KoizumiH. (2006). The methane flux along topographical gradients on a glacier foreland in the high Arctic, Ny-Alesund, Svalbard. Polar Biosci.20, 131–139. doi: 10.15094/00006266
2
AngelR.ClausP.ConradR. (2012). Methanogenic archaea are globally ubiquitous in aerated soils and become active under wet anoxic conditions. ISME J.6, 847–862. doi: 10.1038/ismej.2011.141
3
AngelidakiI.KarakashevD.BatstoneD. J.PluggeC. M.StamsA. J. (2011). Biomethanation and its potential. Methods Enzymol.494, 327–351. doi: 10.1016/B978-0-12-385112-3.00016-0
4
AveryG. B.ShannonR. D.WhiteJ. R.MartensC. S.AlperinM. J. (2003). Controls on methane production in a tidal freshwater estuary and a peatland: methane production via acetate fermentation and CO2 reduction. Biogeochemistry62, 19–37. doi: 10.1023/A:1021128400602
5
BárcenaT. G.YdeJ. C.FinsterK. W. (2010). Methane flux and high-affinity methanotrophic diversity along the chronosequence of a receding glacier in Greenland. Ann. Glaciol.51, 23–31. doi: 10.3189/172756411795932001
6
BelyeaL. R.ClymoR. S. (1998). “Do hollows control the rate of peat bog growth,” in Patterned Mires and Mire Pools. eds. StandenV.TallisJ. H.MeadeR. (London: British Ecological Society), 55–65.
7
BenderM.ConradR. (1995). Effect of CH4 concentrations and soil conditions on the induction of CH4 oxidation activity. Soil Biol. Biochem.27, 1517–1527. doi: 10.1016/0038-0717(95)00104-M
8
BorrelG.ParisotN.HarrisH. M.PeyretailladeE.GaciN.TotteyW.et al. (2014). Comparative genomics highlights the unique biology of Methanomassiliicoccales, a Thermoplasmatales-related seventh order of methanogenic archaea that encodes pyrrolysine. BMC Genomics15, 1–24. doi: 10.1186/1471-2164-15-679
9
BradleyJ. A.SingarayerJ. S.AnesioA. M. (2014). Microbial community dynamics in the forefield of glaciers. Proc. R. Soc. B Biol. Sci. U. S. A.281:20140882. doi: 10.1098/rspb.2014.0882
10
BremnerJ. (1960). Determination of nitrogen in soil by the Kjeldahl method. J. Agric. Sci.55, 11–33. doi: 10.1017/S0021859600021572
11
Capella-GutiérrezS.Silla-MartínezJ. M.GabaldónT. (2009). trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics25, 1972–1973. doi: 10.1093/bioinformatics/btp348
12
ChenX.WangG.ZhangT.MaoT.WeiD.HuZ.et al. (2017). Effects of warming and nitrogen fertilization on GHG flux in the permafrost region of an alpine meadow. Atmos. Environ.157, 111–124. doi: 10.1016/j.atmosenv.2017.03.024
13
ChenH.YaoS.WuN.WangY.LuoP.TianJ.et al. (2008). Determinants influencing seasonal variations of methane emissions from alpine wetlands in Zoige plateau and their implications. J. Geophys. Res. Atmos.113:D12. doi: 10.1029/2006JD008072
14
ChiriE.NauerP. A.HennebergerR.ZeyerJ.SchrothM. H. (2015). Soil–methane sink increases with soil age in forefields of alpine glaciers. Soil Biol. Biochem.84, 83–95. doi: 10.1016/j.soilbio.2015.02.003
15
ChiriE.NauerP. A.RainerE. M.ZeyerJ.SchrothM. H. (2017). High temporal and spatial variability of atmospheric-methane oxidation in alpine glacier forefield soils. Appl. Environ. Microbiol.83, e01139–e02117. doi: 10.1128/AEM.01139-17
16
ColesJ. R.YavittJ. B. (2002). Control of methane metabolism in a forested northern wetland, New York state, by aeration, substrates, and peat size fractions. Geomicrobiol J.19, 293–315. doi: 10.1080/01490450290098397
17
ConradR. (2020). Methane production in soil environments-anaerobic biogeochemistry and microbial life between flooding and desiccation. Microorganisms8:881. doi: 10.3390/microorganisms8060881
18
CostelloA. M.LidstromM. E. (1999). Molecular characterization of functional and phylogenetic genes from natural populations of methanotrophs in lake sediments. Appl. Environ. Microbiol.65, 5066–5074. doi: 10.1128/AEM.65.11.5066-5074.1999
19
DanilovaO. V.SuzinaN. E.Van De KampJ.SvenningM. M.BodrossyL.DedyshS. N. (2016). A new cell morphotype among methane oxidizers: a spiralshaped obligately microaerophilic methanotroph from northern low-oxygen environments. ISME J.10, 2734–2743. doi: 10.1038/ismej.2016.48
20
DengY.CuiX.HernándezM.DumontM. G. (2014). Microbial diversity in hummock and hollow soils of three wetlands on the Qinghai-Tibetan plateau revealed by 16S rRNA pyrosequencing. PLoS One9:e103115. doi: 10.1371/journal.pone.0103115
21
DengY.LiuP.ConradR. (2019). Effect of temperature on the microbial community responsible for methane production in alkaline NamCo wetland soil. Soil Biol. Biochem.132, 69–79. doi: 10.1016/j.soilbio.2019.01.024
22
DingW. X.CaiZ. C. (2007). Methane emission from natural wetlands in China: summary of years 1995-2004 studies. Pedosphere17, 475–486. doi: 10.1016/S1002-0160(07)60057-5
23
DixonP. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci.14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x
24
DongX.ChenZ. (2012). Psychrotolerant methanogenic archaea: diversity and cold adaptation mechanisms. Sci. China Life Sci.55, 415–421. doi: 10.1007/s11427-012-4320-0
25
EdgarR. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods10, 996–998. doi: 10.1038/nmeth.2604
26
EichelJ. (2019). “Vegetation succession and biogeomorphic interactions in glacier forelands,” in Geomorphology of Proglacial Systems. Geography of the Physical Environment. eds. HeckmannT.MorcheD. (Cham: Springer), 327–349.
27
EmersonJ. B.VarnerR. K.WikM.ParksD. H.NeumannR. B.JohnsonJ. E.et al. (2020). Diverse Arctic lake sediment microbiota shape methane emission temperature sensitivity. BioRxiv. [Preprint]. doi: 10.1101/2020.02.08.934661
28
EttwigK. F.ButlerM. K.Le PaslierD.PelletierE.MangenotS.KuypersM. M. (2010). Nitrite-driven anaerobic methane oxidation by oxygenic bacteria. Nature464, 543–548. doi: 10.1038/nature08883
29
EttwigK. F.ZhuB.SpethD.KeltjensJ. T.JettenM. S.KartalB. (2016). Archaea catalyze iron-dependent anaerobic oxidation of methane. Proc. Natl. Acad. Sci. U. S. A.113, 12792–12796. doi: 10.1073/pnas.1609534113
30
GirkinN. T.TurnerB. L.OstleN.CraigonJ.SjögerstenS. (2018). Root exudate analogues accelerate CO2 and CH4 production in tropical peat. Soil Biol. Biochem.117, 48–55. doi: 10.1016/j.soilbio.2017.11.008
31
Guerrero-CruzS.CremersG.van AlenT. A.Op den CampH. J.JettenM. S.RasigrafO.et al. (2018). Response of the anaerobic methanotroph “Candidatus Methanoperedens nitroreducens” to oxygen stress. Appl. Environ. Microbiol.84, e01832–e01918. doi: 10.1128/AEM.01832-18
32
Guerrero-CruzS.VaksmaaA.HornM. A.NiemannH.PijuanM.HoA. (2021). Methanotrophs: discoveries, environmental relevance, and a perspective on current and future applications. Front. Microbiol.12:678057. doi: 10.3389/fmicb.2021.678057
33
HakobyanA.LiesackW. (2020). Unexpected metabolic versatility among type II methanotrophs in the Alphaproteobacteria. Biol. Chem.401, 1469–1477. doi: 10.1515/hsz-2020-0200
34
HammerØ.HarperD.RyanP. (2001). PAST: paleontological statistics software package for education and data analysis. Palaeontol. Electron.4, 1–9.
35
HaroonM. F.HuS.ShiY.ImelfortM.KellerJ.HugenholtzP.et al. (2013). Anaerobic oxidation of methane coupled to nitrate reduction in a novel archaeal lineage. Nature500, 567–570. doi: 10.1038/nature12375
36
HofmannK.PauliH.PraegN.WagnerA. O.IllmerP. (2016). Methane-cycling microorganisms in soils of a high-alpine altitudinal gradient. FEMS Microbiol. Ecol.92:fiw99. doi: 10.1093/femsec/fiw009
37
HofmannK.ReitschulerC.IllmerP. (2013). Aerobic and anaerobic microbial activities in the foreland of a receding glacier. Soil Biol. Biochem.57, 418–426. doi: 10.1016/j.soilbio.2012.08.019
38
HøjL.RustenM.HaugenL. E.OlsenR. A.TorsvikV. L. (2006). Effects of water regime on archaeal community composition in Arctic soils. Environ. Microbiol.8, 984–996. doi: 10.1111/j.1462-2920.2006.00982.x
39
HuB. L.ShenL. D.LianX.ZhuQ.LiuS.HuangQ.et al. (2014). Evidence for nitrite-dependent anaerobic methane oxidation as a previously overlooked microbial methane sink in wetlands. Proc. Natl. Acad. Sci. U. S. A.111, 4495–4500. doi: 10.1073/pnas.1318393111
40
HuttunenJ. T.NykänenH.TurunenJ.MartikainenP. J. (2003). Methane emissions from natural peatlands in the northern boreal zone in Finland, Fennoscandia. Atmos. Environ.37, 147–151. doi: 10.1016/S1352-2310(02)00771-9
41
IPCC (2013). Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, UK: Cambridge University Press.
42
JiT. W. (2005). Comparison on determining the organic matter contents in the soils by different heating methods in the potassium dichromate-volumetric method. Acta Agricult. Zhejiangensis17, 311–313.
43
JinH.WuJ.ChengG.NakanoT.SunG. (1999). Methane emissions from wetlands on the Qinghai-Tibet plateau. Chin. Sci. Bull.44, 2282–2286. doi: 10.1007/BF02885940
44
JuradoV.Gonzalez-PimentelJ. L.MillerA. Z.HermosinB.D’AngeliI. M.TogniniP.et al. (2020). Microbial communities in vermiculation deposits from an alpine cave. Front. Earth Sci.8:586248. doi: 10.3389/feart.2020.586248
45
KniefC.LipskiA.DunfieldP. F. (2003). Diversity and activity of methanotrophic bacteria in different upland soils. Appl. Environ. Microbiol.69, 6703–6714. doi: 10.1128/AEM.69.11.6703-6714.2003
46
KnittelK.BoetiusA. (2009). Anaerobic oxidation of methane: progress with an unknown process. Annu. Rev. Microbiol.63, 311–334. doi: 10.1146/annurev.micro.61.080706.093130
47
KolbS.KniefC.StubnerS.ConradR. (2003). Quantitative detection of methanotrophs in soil by novel pmoA-targeted real-time PCR assays. Appl. Environ. Microbiol.69, 2423–2429. doi: 10.1128/AEM.69.5.2423-2429.2003
48
KoldeR. (2013). Pheatmap: Pretty Heatmaps. R Package Version 0.7.7. Available at: https://cran.r-project.org/web/packages/pheatmap/index.html
49
KongD.LiS.JinY.WuS.ChenJ.HuT.et al. (2019). Linking methane emissions to methanogenic and methanotrophic communities under different fertilization strategies in rice paddies. Geoderma347, 233–243. doi: 10.1016/j.geoderma.2019.04.008
50
KuoS. (1996). “Phosphorus” in Methods of Soil Analysis. Part 3. SSSA Book Series. 5. ed. SparksD. L. (Madison, WI: SSSA), 869–919.
51
LaineA.WilsonD.KielyG.ByrneK. A. (2007). Methane flux dynamics in an Irish lowland blanket bog. Plant Soil299, 181–193. doi: 10.1007/s11104-007-9374-6
52
LarssonA. (2014). AliView: a fast and lightweight alignment viewer and editor for large datasets. Bioinformatics30, 3276–3278. doi: 10.1093/bioinformatics/btu531
53
Le MerJ.RogerP. (2001). Production, oxidation, emission and consumption of methane by soils: a review. Eur. J. Soil Biol.37, 25–50. doi: 10.1016/S1164-5563(01)01067-6
54
LetunicI.BorkP. (2016). Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res.44, W242–W245. doi: 10.1093/nar/gkw290
55
LeuA. O.CaiC.McIlroyS. J.SouthamG.OrphanV. J.YuanZ.et al. (2020). Anaerobic methane oxidation coupled to manganese reduction by members of the Methanoperedenaceae. ISME J.14, 1030–1041. doi: 10.1038/s41396-020-0590-x
56
LinX.WangS.MaX.XuG.LuoC.LiY.et al. (2009). Fluxes of CO2, CH4, and N2O in an alpine meadow affected by yak excreta on the Qinghai-Tibetan plateau during summer grazing periods. Soil Biol. Biochem.41, 718–725. doi: 10.1016/j.soilbio.2009.01.007
57
LudwigW.StrunkO.WestramR.RichterL.MeierH.Yadhukumaret al. (2004). ARB: a software environment for sequence data. Nucleic Acids Res.32, 1363–1371. doi: 10.1093/nar/gkh293
58
LükeC.FrenzelP. (2011). Potential of pmoA amplicon pyrosequencing for methanotroph diversity studies. Appl. Environ. Microbiol.77, 6305–6309. doi: 10.1128/aem.05355-11
59
MaK.ConradR.LuY. (2012). Responses of methanogen mcrA genes and their transcripts to an alternate dry/wet cycle of paddy field soil. Appl. Environ. Microbiol.78, 445–454. doi: 10.1128/AEM.06934-11
60
MäkipääR.LeppänenS. M.MunozS. S.SmolanderA.TiirolaM.TuomivirtaT.et al. (2018). Methanotrophs are core members of the diazotroph community in decaying Norway spruce logs. Soil Biol. Biochem.120, 230–232. doi: 10.1016/j.soilbio.2018.02.012
61
MartineauC.PanY.BodrossyL.YergeauE.WhyteL. G.GreerC. W. (2014). Atmospheric methane oxidizers are present and active in Canadian high Arctic soils. FEMS Microbiol. Ecol.89, 257–269. doi: 10.1111/1574-6941.12287
62
Mateos-RiveraA.ØvreåsL.WilsonB.YdeJ. C.FinsterK. W. (2018). Activity and diversity of methane-oxidizing bacteria along a Norwegian sub-Arctic glacier forefield. FEMS Microbiol. Ecol.94:fiy059. doi: 10.1093/femsec/fiy059
63
MetjeM.FrenzelP. (2005). Effect of temperature on anaerobic ethanol oxidation and methanogenesis in acidic peat from a northern wetland. Appl. Environ. Microbiol.71, 8191–8200. doi: 10.1128/AEM.71.12.8191-8200.2005
64
MinickK. J.KelleyA. M.MiaoG.LiX.NoormetsA.MitraB.et al. (2019). Microtopography alters hydrology, phenol oxidase activity and nutrient availability in organic soils of a coastal freshwater forested wetland. Wetlands39, 263–273. doi: 10.1007/s13157-018-1107-5
65
MorrisR. L.TaleV. P.MathaiP. P.ZitomerD. H.MakiJ. S. (2016). mcrA gene abundance correlates with hydrogenotrophic methane production rates in full-scale anaerobic waste treatment systems. Lett. Appl. Microbiol.62, 111–118. doi: 10.1111/lam.12515
66
NauerP. A.DamB.LiesackW.ZeyerJ.SchrothM. H. (2012). Activity and diversity of methane-oxidizing bacteria in glacier forefields on siliceous and calcareous bedrock. Biogeosciences9, 2259–2274. doi: 10.5194/bg-9-2259-2012
67
PriceM. N.DehalP. S.ArkinA. P. (2010). FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One5:e9490. doi: 10.1371/journal.pone.0009490
68
QiuJ. (2008). China: the third pole. Nat. News454, 393–396. doi: 10.1038/454393a
69
RahalkarM.BussmannI.SchinkB. (2007). Methylosoma difficile gen. Nov., sp. nov., a novel methanotroph enriched by gradient cultivation from littoral sediment of lake constance. Int. J. Syst. Evol. Microbiol.57, 1073–1080. doi: 10.1099/ijs.0.64574-0
70
RozewickiJ.LiS.AmadaK. M.StandleyD. M.KatohK. (2019). MAFFT-DASH: integrated protein sequence and structural alignment. Nucleic Acids Res.47, W5–W10. doi: 10.1093/nar/gkz342
71
SchlossP. D.WestcottS. L.RyabinT.HallJ. R.HartmannM.HollisterE. B.et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol.75, 7537–7541. doi: 10.1128/AEM.01541-09
72
Sutton-GrierA. E.MegonigalJ. P. (2011). Plant species traits regulate methane production in freshwater wetland soils. Soil Biol. Biochem.43, 413–420. doi: 10.1016/j.soilbio.2010.11.009
73
TooC. C.KellerA.SickelW.LeeS. M.YuleC. M. (2018). Microbial community structure in a Malaysian tropical peat swamp forest: the influence of tree species and depth. Front. Microbiol.9:2859. doi: 10.3389/fmicb.2018.02859
74
TveitA. T.UrichT.SvenningM. M. (2014). Metatranscriptomic analysis of arctic peat soil microbiota. Appl. Environ. Microbiol.80, 5761–5772. doi: 10.1128/AEM.01030-14
75
WadhamJ. L.CooperR. J.TranterM.BottrellS. (2007). Evidence for widespread anoxia in the proglacial zone of an Arctic glacier. Chem. Geol.243, 1–15. doi: 10.1016/j.chemgeo.2007.04.010
76
WangY.CuiH.SuX.WeiS.ZhuY.LuZ.et al. (2020). Diversity and distribution of methanogenic community between two typical alpine ecosystems on the Qinghai–Tibetan plateau. Curr. Microbiol.77, 1061–1069. doi: 10.1007/s00284-020-01891-x
77
WangQ.QuensenJ. F.FishJ. A.LeeT. K.SunY.TiedjeJ. M.et al. (2013). Ecological patterns of nifH genes in four terrestrial climatic zones explored with targeted metagenomics using FrameBot, a new informatics tool. mBio4, e00592–e00513. doi: 10.1128/mBio.00592-13
78
WangY.WangY. (2003). Quick measurement of CH4, CO2 and N2O emissions from a short-plant ecosystem. Adv. Atmos. Sci.20, 842–844. doi: 10.1007/BF02915410
79
WeiD.Xu-RiTarchenT.DaiD.WangY.WangY. (2015a). Revisiting the role of CH4 emissions from alpine wetlands on the Tibetan plateau: evidence from two in situ measurements at 4758 and 4320m above sea level. J. Geophys. Res. Biogeosci.120, 1741–1750. doi: 10.1002/2015JG002974
80
WeiD.Xu-RiTarchenT.WangY.WangY. (2015b). Considerable methane uptake by alpine grasslands despite the cold climate: in situ measurements on the central Tibetan plateau, 2008-2013. Glob. Chang. Biol.21, 777–788. doi: 10.1111/gcb.12690
81
WoutersB.GardnerA. S.MoholdtG. (2019). Global glacier mass loss during the GRACE satellite mission (2002–2016). Front. Earth Sci.7:96. doi: 10.3389/feart.2019.00096
82
WuH.WangX.GanjurjavH.HuG.QinX.GaoQ. (2020). Effects of increased precipitation combined with nitrogen addition and increased temperature on methane fluxes in alpine meadows of the Tibetan plateau. Sci. Total Environ.705:135818. doi: 10.1016/j.scitotenv.2019.135818
83
XuX. F.TianH. Q.ZhangC.LiuM. L.RenW.ChenG. S.et al. (2010). Attribution of spatial and temporal variations in terrestrial methane flux over North America. Biogeosciences7, 3637–3655. doi: 10.5194/bg-7-3637-2010
84
YangS.LiebnerS.AlawiM.EbenhöhO.WagnerD. (2014). Taxonomic database and cut-off value for processing mcrA gene 454 pyrosequencing data by MOTHUR. J. Microbiol. Methods103, 3–5. doi: 10.1016/j.mimet.2014.05.006
85
YangS.LiebnerS.WinkelM.AlawiM.HornF.DörferC.et al. (2017). In-depth analysis of core methanogenic communities from high elevation permafrost-affected wetlands. Soil Biol. Biochem.111, 66–77. doi: 10.1016/j.soilbio.2017.03.007
86
YaoP. (2014). Climate change during the past 2000 years recoeded by Tanggula ice core in central Tibetan plateau. Dissertation/Doctoral thesis. Beijing: University of Chinese Academy of Sciences.
87
YaoT.PuJ.LuA.WangY.YuW. (2007). Recent glacial retreat and its impact on hydrological processes on the Tibetan plateau, China, and surrounding regions. Arct. Antarct. Alp. Res.39, 642–650. doi: 10.1657/1523-0430(07-510)[YAO]2.0.CO;2
88
YaoT.WangY.LiuS.PuJ.ShenY.LuA. (2004). Recent glacial retreat in high Asia in China and its impact on water resource in Northwest China. Sci. China Ser. D Earth Sci.47, 1065–1075. doi: 10.1360/03yd0256
89
YuanJ.YuanY.ZhuY.CaoL. (2018). Effects of different fertilizers on methane emissions and methanogenic community structures in paddy rhizosphere soil. Sci. Total Environ.627, 770–781. doi: 10.1016/j.scitotenv.2018.01.233
90
ZhangY.CuiM.DuanJ.ZhuangX.ZhuangG.MaA. (2019). Abundance, rather than composition, of methane-cycling microbes mainly affects methane emissions from different vegetation soils in the Zoige alpine wetland. MicrobiologyOpen8:e00699. doi: 10.1002/mbo3.699
91
ZhangG.TianJ.JiangN.GuoX.WangY.DongX. (2008). Methanogen community in Zoige wetland of Tibetan plateau and phenotypic characterization of a dominant uncultured methanogen cluster ZC-I. Environ. Microbiol.10, 1850–1860. doi: 10.1111/j.1462-2920.2008.01606.x
Summary
Keywords
methane flux, methanogens, methanotrophs, glacier foreland, Tibetan Plateau
Citation
Xing T, Liu P, Ji M, Deng Y, Liu K, Wang W and Liu Y (2022) Sink or Source: Alternative Roles of Glacier Foreland Meadow Soils in Methane Emission Is Regulated by Glacier Melting on the Tibetan Plateau. Front. Microbiol. 13:862242. doi: 10.3389/fmicb.2022.862242
Received
25 January 2022
Accepted
01 March 2022
Published
21 March 2022
Volume
13 - 2022
Edited by
Peng Xing, Nanjing Institute of Geography and Limnology (CAS), China
Reviewed by
Ang Hu, Hunan Agricultural University, China; Michael Pester, German Collection of Microorganisms and Cell Cultures GmbH (DSMZ), Germany
Updates
Copyright
© 2022 Xing, Liu, Ji, Deng, Liu, Wang and Liu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pengfei Liu, liupf@lzu.edu.cn; liupfskygre@gmail.comYongcui Deng, dengyongcui@njnu.edu.cn
This article was submitted to Terrestrial Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.