Abstract
Wetlands are the largest natural sources of methane (CH4) emissions worldwide, with methanogenic archaea serving as the primary drivers of CH4 production. Nevertheless, the influences of biotic factors (e.g., methanogen abundance, community diversity and composition) and abiotic factors (e.g., soil properties) on potential CH4 production rates remain insufficiently understood in Qinghai-Tibet Plateau. In this study, we examined soil properties, potential methane production rates (PMPRs), methanogenic archaeal abundance, diversity, community structure, and co-occurrence networks across four wetlands (two desert wetlands and two peatlands) with contrasting soil conditions on the northeastern edge of the Qinghai-Tibet Plateau. We found no significant differences in methanogen abundance and PMPRs among the four wetlands, but the two were significantly positively correlated. The structure of methanogenic communities varied markedly among wetlands and was mainly shaped by soil pH. The complexity of co-occurrence networks was positively correlated with both methanogen diversity and PMPRs. Further analysis using partial least squares path modeling (PLS-PM) revealed that PMPRs were closely associated with soil nutrition (soil total organic carbon and total nitrogen; standardized path coefficient = 0.307), methanogenic abundance (0.570) and network complexity (0.238). It indicated that biotic factors may exert a greater influence than abiotic factors on soil PMPRs in wetland ecosystems. Additionally, complex microbial interaction networks may play a more crucial role in regulating PMPRs than methanogenic diversity and community structure. Our study highlights a strong link between methanogenic network complexity and methane-producing potential, offering a novel perspective on the relationship between community interactions and ecosystem function.
1 Introduction
The concentration of methane (CH4) in the troposphere has increased by 150% since the onset of industrialization (Etminan et al., ). Compared to carbon dioxide, methane is 72 times more potent in trapping heat radiation over a 20-year period (Stavert et al., 2022). As a result, it is expected to play a significant role in driving future global warming caused by the greenhouse effect. Wetlands are the largest natural sources of CH4 (Saunois et al., 2020). Ground-penetrating radar and pore water-sampling efforts have shown that wetland belowground topography can shape hydrology in ways which enhance downward transport of labile organic carbon into deep peat, encouraging methanogenesis meters below the surface (Abdalla et al., ). Among various wetland ecosystems, freshwater wetlands (138–165 Tg CH4 yr−1) and lakes (23–142 Tg CH4 yr−1) constitute the dominant methane emitters. These are followed by rice paddies (25–32 Tg CH4 yr−1), reservoirs (9–28 Tg CH4 yr−1), coastal oceans (<200 m depth; 5–28 Tg CH4 yr−1), and rivers/streams (2–21 Tg CH4 yr−1) (Rosentreter et al., 2021). Aquatic CH4 is predominantly generated by methanogenic archaea in soils during the terminal stage of anaerobic organic matter degradation (Hao et al., ; Nazaries et al., 2013; Gruca-Rokosz et al., ; Rosentreter et al., 2021; Stavert et al., 2022).
Archaea-mediated methanogenesis accounts for approximately 69% of global atmospheric methane emissions (Conrad, ). Methanogens are key microorganisms influencing methane production, as their diversity, abundance, and metabolic activity directly determine the rate and total amount of methane generated. Some studies have revealed that soil properties and methanogenic microbial communities can influence the function of methanogens (Brooker et al., ; Deng et al., ; Yang et al., 2022). For example, pH plays an important role in habitat filtering, which shapes the methanogenic biogeographic pattern in paddy soils, lakes, and dry lands (Hu et al., ). In coastal wetland systems such as mangroves, the abundance of labile organic matter is considered a major driver of methane emissions (Xiang et al., 2015). Meanwhile, biotic factors also drive methane emissions in soils. Archaeal methanogens genes and bacterial methylphosphonate degradation genes were both positively correlated with methane flux in salt ponds (Zhou et al., 2022). In an ombrotrophic peat bog in Maine, USA, researchers observed variation in methanogen abundance explaining 70% of the flux heterogeneity in a subset of plots (Arnold et al., ). These suggested that the effects of biotic and abiotic factors on the function of methanogenic communities remain controversial.
Obtaining detailed knowledge of complex species interactions in natural environments through empirical studies remains a significant challenge, particularly for the most abundant and diverse microbial taxa (Layeghifard et al., 2017). In recent years, co-occurrence network analysis has become a widely used approach in ecology for inferring potential microbial interactions (Barberán et al., ). Methanogens engage in complex associations that encompass both interspecies and intraspecies syntrophic relationships, as well as competitive interactions (Li et al., 2021). Although coexistence cannot be strictly equated with co-occurrence, co-occurrence patterns offer valuable insights into potential coexistence relationships, spanning from taxon pairs to complex multi-taxon communities across diverse ecosystems (Eiler et al., ; Tavella and Cagnolo, 2019). While previous studies have made important contributions in describing methanogenic community composition and diversity, as well as in elucidating the effects of biotic factors on CH4 fluxes in wetlands (Feng et al., ; Dong et al., ; Qin et al., 2024; Yang et al., 2024), few have examined the interactions among soil methanogens themselves, which are likely to play a more critical role in the functioning of complex ecosystems.
The Qinghai-Tibet Plateau plays a critical role in both regional and global climate systems and ecological processes. It is also a sensitive indicator of climate change (Chen et al., ), and a critical component of the global carbon cycle. It has been called the “third pole” of the Earth, and is a major CH4 emission hotspot (Ding and Cai, ). We collected 128 soil samples from four different wetland ecosystems, which vary in habitat characteristics and soil properties on the northeastern edge of the Qinghai-Tibet Plateau (Supplementary Figure 1). We analyzed soil properties, measured soil potential CH4 production rates (PMPRs) in lab, assessed the abundance and community composition of methanogens. Studying methanogenic communities in diverse wetland habitats provides a broader understanding of their variation and enhances the generalizability of findings, and offers new perspectives on the mechanisms driving methanogen function and adaptation. Specifically, we address the following questions: (i) Whether abiotic factors or biotic factors have a greater influence on potential CH4 production rates. (ii) More complex networks may enhance ecosystem functioning (Morriën et al., 2017; Zhang et al., 2024), though whether this applies to methanogen networks in promoting soil methanogenesis remains unclear.
2 Materials and methods
2.1 Study site description and sampling
Our samples were collected from four wetland ecosystems on the northeastern edge of the Qinghai-Tibet Plateau. The four wetlands are located in Maqu (MQ), Luqu (LQ), Zhangye (ZY) and Sugan lake (SGH), respectively (Supplementary Table 1 and Figure 1). MQ and LQ are peatlands located in the eastern Qinghai-Tibet Plateau, while ZY and SGH are desert wetlands situated in the arid northern region, characterized by low precipitation and a dry climate (Zhou et al., 2019).
Figure 1
On July and August 2020, eight 2 × 2-m2 plots (>50 m interval) were randomly arranged from each region, respectively. The soil samples were collected from five points (the four corners and center) in each plot and subsequently blended to obtain a composite sample. At each plot, soil samples were collected at four depths, 0–5, 5–10, 10–20, 20–40 cm, respectively. In total, 128 soil samples were collected in quadruplicate from eight plots and four depths. The samples were then transported to the laboratory and divided into two parts as follows: one part was air dried and sieved through a 2 mm mesh for the soil properties analysis; the other part was stored at −20°C before DNA extraction and the incubation experiment.
2.2 Soil properties analysis
Weighed 10 g air-dried soil samples, added 50 ml of 1 mol/L KCl solution, then measured the soil pH by a pH meter (Sartorius PB-10; Sartorius AG, Göttingen, Germany) (Xi et al., 2023). Soil electrical conductivity (EC) was measured in water (1:5 w/v) by conductivity meter (CT-3031) (Fang et al., ). The soil moisture content (SWC) was determined by drying the fresh soil to a constant weight at 105°C, 5°C. The soil total organic carbon (TOC) and soil total nitrogen (TN) were determined by Elementar analysis system CHNS (Elementar Analysensysteme GmbH, Langenselbold, Germany) (Yang et al., 2020). The soil total phosphorus (TP) was analyzed by digestion with sulfuric acid at 375°C, 5°C. Soil available ( and ) were extracted using 2 M KCl (1:5 w/v) and analyzed by FIAstar 5000 Analyzer (FOSS, Hillerød, Denmark) (Peng et al., 2024). Due to logistical constraints, soil available nitrogen was measured using air-dried samples. While fresh soil is generally preferred for such measurements, all samples were processed under the same conditions, ensuring relative comparability across samples. The results of soil properties are presented in Supplementary Table 2.
2.3 Potential methane production rates (PMPRs)
Anaerobic incubation experiments were conducted to measure sediment CH4 production potential (Yang et al., 2022; Mu et al., 2025). Approximately 5 g of each soil sample was put into a 120 mL serum vials capped with two-way valves and in deionized water was added in 1:1 v/v. Before the start of incubation, the bottles were flushed with pure nitrogen gas (N2) for 5–8 min to create an anoxic condition (Vizza et al., 2017). From each sample, we prepared three replicates for quality control and one control filled only with N2. CH4 concentrations were measured on the first and fifteenth day use a gas chromatograph (Agilent 7890, Santa Clara, CA, USA) equipped with a flame ion detector (FID) using 10 ml gas samples from the bottle. The furnace temperature, FID, and ECD detector temperature were 55, 200, and 300°C, respectively. 99.999% high purity nitrogen was selected as the carrier gas, and the flow rate was 2 mL/min. High purity hydrogen and air were used as the gas with flow rates of 40 and 400 mL/min, respectively. Soil potential methane production rates (PMPRs) were calculated according to soil incubation time and gas concentration as the following equation:
PP is the flux of CH4 in soil samples [ng CH4 g−1 (dry weight) day−1]; dc/dt is the rate of change in headspace CH4 in the incubation bottle over time (mmol mol−1 d−1); MM is the molar mass (g mol−1) of CH4 (g); VH is the volume of serum bottles headspace (L); PA is the atmospheric pressure (kPa); R is the gas constant (m3 Pa °K−1 mol−1); WS is the dry weight of soil sample (g); TST and T are the standard temperature (°K) and the incubation temperature (°K), respectively (Li et al., 2021; Yang et al., 2022).
2.4 DNA extraction and gene qPCR
DNA extraction was performed from 0.25 g soil samples with DNeasy Power Soil kit (Qiagen, Germantown, MD, United States). The methyl coenzyme M reductase (mcrA) gene, which catalyzes the final step in all methanogenic pathways by reducing the methyl group attached to coenzyme M, is widely utilized as a functional gene marker in the characterization of methanogenic communities (Luton et al., 2002). The copy numbers of mcrA gene were determined by real-time PCR with primer set mlas (55″- GGTGGTGTMGGDTTCACMCARTA) – rev (55″ - CGTTCATBGCGTAGTTVGGRTAGT) (Ma et al., 2010; Angel et al., ) and a SYBR Green System (Takara Bio Inc., Shiga Japan) as described previously (Costello and Lidstrom, ; Kolb et al., 2003), Technical replicates were performed in triplicate. The assays were performed using a Q5 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) and the associated software. The 20 μl reaction mixtures contained: 2 μl template DNA, 10 μl SYBR Green, 0.4 μl ROX mixture (2×, Takara Bio Inc., Shiga, Japan), 0.4 μl forward primer (10 μmol), 0.4 μl reverse primer (10 μmol l−1), and 6.8 μl nuclease-free water. Standard curves were constructed using plasmids harboring the gene fragment. PCR runs started with an initial denaturation and enzyme activation step at 98°C8°C for 2 min, followed by 40 cycles of 10 s at 98°C, 8°C, 30 s at 60°C, 0°C and 40 s at 72°C, 2°C, and 10 min at 72°C, 2°C. We recorded the fluorescence signal at 80°C, 0°C to attenuate influences of primer dimers. The specificities of PCR products were tested by melting curve analysis. The R2 and amplification efficiency of the standard curve were 0.99 and 96.2%, respectively.
2.5 Illumina sequencing and bioinformatic analysis
We used high-throughput sequencing to characterize the community compositions of methanogensis. The soil DNA were extracted and then sent to the Majorbio Biotechnology Company (Shanghai, China) for Illumina MiSeq sequencing. For each DNA sample, amplified the mcrA gene in triplicate utilizing specific primers (mlas-rev) equipped with a unique barcode. Subsequently, the triplicate amplicons were combined and subjected to a purification process. The purified PCR products from each sample were quantified and pooled in an equimolar fashion, followed by the construction of Illumina libraries using the MiSeq Reagent Kit v3 (Illumina, USA). Sequencing was conducted as paired-end or single-direction format by Majorbio Company (Shanghai, China) on an Illumina MiSeq PE300 platform. Kanzhege, due to unsuccessful library preparation of certain DNA samples, the final available sample count is 121 (MQ:31 samples, LQ:32 samples, ZY: 29 samples, SGH: 29 samples).
A total of 20,71,622 sequences were obtained from 121 soil samples. QIIME 2 (2,022.8) was used to analyze the sequencing data. Specifically, paired-end reads were trimmed to a minimum Q-score of 20. These sequences mcrA gene were processed to generate amplicon sequence variants (ASVs) by DADA2. Insertions and deletions caused the frameshifts of the mcrA gene sequences were corrected using FrameBot (Wang et al., 2013) from FunGene database (http://fungene.cme.msu.edu). The resulting ASVs were assigned to the database from (http://doi.org/10.5880/GFZ.4.5.2014.001) (Yang et al., 2014), to further proofread the taxonomical assignments, the representative sequences of dominant ASVs were compared with the GenBank database using BLAST (https://blast.ncbi.Nlm.nih.gov). The ASVs that contained less than 10 reads were removed from the datasets, and all samples were rarified to 1,016 sequences based on the lowest sequencing depth. After processing, the total number of representative ASVs across all samples was 3,147. The sequencing depth of amplicon sequencing was estimated using rarefaction analyses (Supplementary Figure 2). All raw sequencing data of mcrA were submitted to the Sequence Read Archive of NCBI under the accession numbers PRJNA1227586 (https://www.ncbi.nlm.nih.gov/sra/PRJNA1227586).
2.6 Microbial network construction and keystones
The co-occurrence patterns observed in the networks are indicative of the interactions among microorganisms in the ecosystems (Ma et al., 2016). Network structures were calculated in R4.2.1 using the “picante” package (Kembel et al., ) and visualized using the interactive platform Gephi 0.10.1 using directed network and the Fruchterman-Reingold layout. We considered a valid co-occurrence event to have a Spearman's correlation coefficient (R > 0.6, P < 0.01). A set of metrics including the number of nodes [NN], number of edges[NE], connectance, average degree [AD], global clustering coefficient [GCC], average clustering coefficient [ACC], average neighborhood [AN], average clustering centralization [ADC], degree centralization [DC], number of positive nodes and connectance was calculated to describe the network.
The Zi-Pi thresholds were based on topology heterogeneous properties of the network structure in “igraph” package (Csardi and Nepusz, ), and we sorted all species into four groups: peripherals (zi ≤ 2.5; pi ≤ 0.62), connector hubs (zi ≤ 2.5; pi > 0.62), module hubs (zi > 2.5; pi ≤ 0.62), and network hubs (zi > 2.5; pi > 0.62) (Olesen et al., 2007).
2.7 Network complexity index
A comprehensive index was established to reflect the complexity of microbial network, and was calculated by averaging the standardized scores of the topological properties including the NN, NE, connectance, AD, GCC, ACC, AN, ADC, DC and number of positive edges, as follows:
where Xraw, Xmin, and Xmax represent the raw topological properties, the minimum and maximum values across all samples, respectively (Zhang et al., 2024).
2.8 Statistical analysis
Since there were no significant differences between soil samples from different depths at the same site (Supplementary Table 3), they were treated as independent samples for subsequent analysis. All statistical analyses were performed using R4.2.1 (https://www.R-project.org/). One-way analysis of variance (ANOVA) with Tukey's test was uesd to evaluate differences in mcrA gene copies and PMORs among regions. The data derived from qPCR (i.e., gene copy numbers) were log10 transformed and used in the following analyses. Based on the rarified ASV tables of methangentic group, we calculated the ASV richness and the relative abundance of each ASV in each sample. The relative abundance of taxonomical taxa at the levels of genus, family or class were also calculated on the basis of the ASV tables with taxonomical annotations. Regional differences in microbial communities were assessed using constrained analysis of principal coordinates (CAP) based on Bray-Curtis distance with the “vegan” package (Jiao et al., ). CAP is particularly suitable for analyzing non-normally distributed ecological data. Compared to traditional linear methods like RDA, CAP offers improved resolution in detecting and visualizing group-level differences, making it ideal for examining microbial community variation across distinct wetland habitats. Therefore, to evaluate environmental variables significantly contributing to ASV composition variation among samples, we performed CAP with stepwise forward selection and 999 Monte Carlo permutations (using the “ordistep” function from the “vegan” package), and confirmed feature significance via permutational MANOVA (PERMANOVA). Variation partitioning was performed with adjusted R2 to determine the proportion of variation in bacterial communities explained purely by environmental factors and shared variation of environmental. The “pheatmap” package in R (https://github.com/raivokolde/pheatmap) was used to generate heatmaps of soil properties, biotic factors and topo properties. Random forest analysis (RFA) was performed using the ‘randomForest' package in R (Jin et al., ). Linear correlation between variables was tested by calculating Pearson's correlation coefficient. Variation partitioning analysis (VPA) to discern the contributions of these variables to the overall variations of the PMPRs (using “vegan” package).
Partial least squares path modeling (PLS-PM) is an effective statistical method for studying cause and effect relationships among biotic and abiotic variables (Barber et al., ; Wagg et al., 2014; Fan et al., ). Based on the previous CAP screening results for soil properties, the nutritional and non-nutritional factors in the PLS-PM model were determined. We conducted principal component analysis (PCA) to create a new index to represent the methanogenic community of soils before PLS-PM (Qin et al., 2021). The latent variable “diversity” includes methanogensis richness and shannon, “community” includes PC1 of methanogenic community composition. Then, we used PLS-PM analysis to evaluate the potential causal relationships between the variables to PMPRs. Path coefficients were assessed for difference from 0 by bootstrapping using 1,000 resamples; this allowed calculation of the precision of each path and the direction and strength of the linear relationships between variables (direct effects). Indirect effects are the multiplied path coefficients between a predictor and a response variable, adding the product of all possible paths excluding the direct effect. For PLS-PM, the function “inner plot” was used and calculated by the package “plspm” in R (Monecke and Leisch, 2012). The model reliability was evaluated using the Goodness of Fit (GoF) statistic (Monecke and Leisch, 2012).
3 Results
3.1 Correlations of methanogenic abundance and PMPRs
In our study, the log-transformed mcrA gene abundance ranged from 1.648 ± 0.037 to 1.970 ± 0.125. The potential methane production rates (PMPRs) ranged from 49.668 ± 5.557 ng g−1 dry soil d−1 to 83.361 ± 13.377 ng g−1 dry soil d−1, with no significantly differences among the sampled sites (Figure 1a). Linear regression analysis was conducted to examine the relationship between mcrA copies and PMPRs across four sites and all samples (Figure 1b). A significant positive correlation was observed in MQ, LQ, SGH and all samples, while no significant correlation was found in ZY.
3.2 Methanogenic communities associated with environmental factors
Taxonomic annotations of these ASVs offer limited information, 99.1% of ASVs belong to the Archaea domain, and 33.2% of total ASVs taxonomically assigned only up to the phylum level (see the Supplementary material). The top ten dominant ASVs in the four regions exhibited significant differences, with no shared ASVs among them (Figure 2a). The most dominant methanogen ASVs in the four wetlands were ASV2300, ASV2815, ASV1415, and ASV1380, respectively. The alpha diversity of soil methanogenic communities in ZY respectively was strikingly higher than other sites (Figure 2b).
Figure 2
The Constrained analysis of principal coordinates (CAP) based on Bray-Curtis distance was performed to investigate the relationship between methanogenic community composition and environmental factors. TN, TOC, pH and EC were found to significantly influence methanogenic community distribution (P < 0.05) after forward selection procedure (Figure 2c, Table 1). The first axis was positively correlated with EC but negatively correlated with pH, TOC and TN. The second axis was positively correlated with TOC and TN, but negatively correlated with pH and EC. Among the variables, pH had the strongest influence on the difference among the bacterial communities of all samples.
Table 1
| Factor | Df | Pseudo-F | P-value | Significance |
|---|---|---|---|---|
| pH | 1 | 1.94 | 0.001 | *** |
| EC | 1 | 1.75 | 0.001 | *** |
| TN | 1 | 1.70 | 0.001 | *** |
| TOC | 1 | 1.68 | 0.001 | *** |
ANOVA of the environmental factors correlated to the methangens beta-diversity in CAP analysis.
Significance: 0 = “***”; 0.001 = “**”; 0.01 = “*”; 0.05 = “.”; 0.1 =.”
Permutations = 999.
Overall model significance: pseudo-F = 3.6249, P < 0.001***.
3.3 Co-occurrence networks of methanogenic community and network complexity
We analyzed microbial networks and key ASVs that connect community nodes and maintain stability across different regions. Molecular ecological networks of methanogenic community were constructed at the ASVs level for MQ, LQ, ZY, and SGH regions (Supplementary Figure 3). Supplementary Table 5 summarizes the network parameters. Correlation analysis showed predominantly positive relationships among methanogens in all regions. ZY had the most nodes and edges, the smallest geodesic length, and the highest degree, indicating a more complex and interconnected network compared to others, while SGH had the fewest nodes and edges, suggesting a simpler structure. In addition, we conducted further analysis on the Zi-Pi relationships among ASVs and found 34 provincial hubs and 29 connectors were identified as keystones in all site (Supplementary Table 4). We also calculated the topological properties of the methanogen subnetwork in each sample and assessed the network complexity. The network complexity index for methanogens exhibited a creasing trend with increasing PMPRs (Figure 3). Correlation analysis indicated a significant positive relationship between these properties and methanogenic community richness, pH and EC, a negative correlation with soil nutrients (Supplementary Figure 4).
Figure 3
3.4 Correlations among biotic factors, abiotic factors, and PMPRs
We conducted principal component analysis (PCA) to create a new index to represent the methanogenic community of soils (Supplementary Figure 5). Random forest analysis indicates that soil properties (TOC, SWC, TN, TP, EC, pH, , and ), diversity (richness), microbial community structure (PC1 of PCA analysis) and network complexity to PMPRs (Figure 4). The abundance of mcrA is the most significant factor contributing to PMPRs, with the variable importance exceeding 40%. Following closely is network complexity. Variation partitioning analysis (VPA) revealed that biotic factors independently explained 28% of the variation in PMPRs, while abiotic factors accounted for only 4%. An additional 18% was jointly explained by the interaction between biotic and abiotic factors. In total, these variables explained 50% of the observed variation in PMPRs (Supplementary Figure 6).
Figure 4
We further explored the direct and indirect effects of soil properties, mcrA abundance, the diversity of community and the network complexity on PMPRs by PLS-PM. The goodness of fit of this module is 0.57. The solid line in the figure shows a significant relationship. Generally, PMPRs are mainly directly influenced by mcrA abundance (0.57), soil nutrition (0.307) and the network complexity (0.238) of methanogens. methanogenic diversity indirectly affected PMPRs (Figure 5). Soil nutrition had significant effects on mcrA abundance (0.228), methanogenic diversity (−0.31) and community structure (0.498). Soil non-nutrition, as well as the methanogenic diversity (−0.403) and community structure (−0.303). methanogenic diversity and community structure do not have a significant direct effect on PMPRs.
Figure 5
4 Discussion
4.1 Biotic and abiotic drivers of PMPRs
In this study, we observed significant differences in microbial community composition among the sampling sites, with no shared top 10 dominant methanogenic ASVs across locations (Figure 2a). The differences in community composition were primarily influenced by pH, EC, TOC, and TN, with pH having the most significant effect (Figure 2c). It should be noted that pH can be regarded as a cumulative and integrative parameter that can potentially mask and/or reflect unmeasured environmental drivers and their interactions, such as vegetation, hydrography or chemical concentrations. Soil pH can lead to differences in the community composition of methanogens. Methanoregulaceae are found across a wide pH range (Juottonen, ). Among them, Methanoregula predominated as a bioindicator of acidic niches (pH < 4.9) (Seppey et al., 2023). The methylotrophic methanogens Methanomassiliicoccales served as important bioindicators for only one niche—the niche with a pH ≥ 4.9 (Seppey et al., 2023). In studies with a pH range of 4.0 to 10.0, extreme pH conditions reduced the relative abundance of acetoclastic methanogens, while the abundance of obligate hydrogenotrophic and facultative acetoclastic/hdrogenotrophic methanogens increased (Qiu et al., 2023). It has thus been suggested that the low pH of bogs causes a fundamental disconnect between acetogenesis and acetoclastic methanogenesis (Yavitt and Seidman-Zager, 2006). In addition, since the study sites are located in different habitats, endemism or dispersal effects may lead to the formation of distinct communities within different habitats (Von Eggers et al., 2024), geographic distance may still contribute to variations in microbial community structure across sites. Therefore, the observed differences are likely driven by the combined effects of soil physicochemical properties and geographic isolation.
Our results indicate that soil nutrient promotes PMPRs. Soil TN and TOC can exert a positive influence on PMPRs, either directly or indirectly (Figure 5). The anoxic environment inhibits the aerobic decomposition of organic matter, resulting in SOC accumulation, which increases the diversity of CH4-cycling microbial communities and mcrA gene abundance (Yang et al., 2024). The relationships between availability of organic carbon provide important substrates for microbes and play an important role in CH4 production (Wik et al., 2016). Abundant organic matters in the sediments can reduce competition between sulfate reducers and methanogens by providing more competitive substrate and/or providing more noncompetitive substrate for methanogenesis (Zhuang et al., 2016). Furthermore, Ammonia nitrogen can also serve as a nitrogen resource for microbes to deliver methanogenic substrates, thereby increasing the availability of organic carbon for methanogens (Banger et al., ). Recent studies on CH4 emission prediction from thermokarst lakes have also shown that the DOC and ratios of SOC and TN in sediments had direct positive effects on CH4 release (Mu et al., 2023). This suggests that soil organic matter quality directly promates PMPRs in wetlands.
Although soil nutrients contribute to enhancing potential CH4 production rates, these results suggest that biotic factors may exert a stronger influence on PMPRs than abiotic factors. We found significant positive correlation was observed between mcrA gene abundance and PMPRs (Figure 1b), suggesting that mcrA gene abundance can regulate soil methanogenic potential. The random forest and PLS-PM analyses yielded consistent results, indicating that mcrA gene abundance had the most significant impact on PMPRs among all research factors (Figures 4, 5). The influence of mcrA gene abundance and network complexity on PMPRs is significantly greater than that of soil nutrients. The results of the VPA analysis indicate that biotic factors alone explained 28% of the variation in PMPRs (Supplementary Figure 6), which is significantly higher than the explanatory power of abiotic factors. It is precisely for this reason that, despite the significant differences in soil properties among the four wetlands, there was no significant difference in PMPRs.
4.2 Network complexity promotes PMPRs
Building on previous research into community composition and structure, we examined the interactions among soil methanogens—an aspect often overlooked in earlier studies. Interestingly, we found that network complexity showed a stronger positive effect on PMPRs than methanogen diversity and composition, a pattern consistently observed across different wetland habitats (Figures 4, 5). In this study on co-occurrence networks, weaker correlations were filtered out. Correlation analysis revealed predominantly positive relationships among methanogensis across all habitats (see Supplementary Table 5), signifying cooperative rather than competitive interactions within these communities. We found that the richness of methanogens was positively correlated with specific topological features such as AD, number of nodes and edges, complexity and AN (Supplementary Figure 4). Conversely, methanogenic network might be more susceptible to external interference, potentially leading to a less stable structure (Faust et al., ), sensitive to changes or disruptions in SGH. Similar results were observed in rice paddies (Li et al., 2021). This result suggested that the microbially derived ecological processes are not necessarily captured by the sum of its coexisting individuals. Rather, these are a consequence of integrated metabolic pathways that are conducted by a myriad of interactions among taxa. More tightly connected microbial members supported a higher level of ecosystem functions, which could be associated with a higher efficiency of resource use and metabolic regulation of ecological processes (Morriën et al., 2017). The combined effects of mechanisms that alter methanogenesis, methanotrophy, nitrogen cycling, and ammonium release, along with those that enhance decomposition and promote the growth of syntrophic bacterial populations, could collectively contribute to increased net CH4 flux in wetlands (Hartman et al., ). These network structures provide us with deeper insights into the interconnections among microbes and the ecological assembly rules, far beyond the mere understanding of diversity and community composition (Ziegler et al., 2018).
Additionally, our findings reveal that network complexity and specific topological properties (such as GCC, ACC, connectance, AD, and AN) are negatively correlated with soil nutrients (Supplementary Figure 4). As soil nutrients decline, resource input becomes limited, prompting microbes to intensify their interactions both within and between species (e.g., cooperation, mutualism) and enhance their metabolic capacity to acquire these limited resources (Brown et al., ; Jansson and Hofmockel, ), particularly in relation to carbon decomposition (Ren et al., 2021).
5 Conclusion
In summary, this study advances the existing framework for assessing the relationships between biotic and abiotic factors and potential CH4 production rates by incorporating an evaluation of methanogenic microbial co-occurrence networks—an essential yet previously underexplored dimension of methanogenic microbial diversity. Our findings clearly indicate that the complexity of methanogenic microbial networks is a key role of soil methanogenic microbial diversity in methane production processes. Integrating this perspective enables a more mechanistic understanding of how interactions among methanogenic microorganisms shape potential methane production and broader ecosystem functioning in natural environments. These results highlight the need to shift focus from the total number of methanogenic species to their association networks, which may help reduce uncertainty in biodiversity–ecosystem function relationships. The strong correlations between methanogenic microbial networks and PMPRs underscore the importance of expanding both empirical and theoretical investigations into the network structures of methanogenic communities. Future predictions of potential CH4 production rates from wetlands should take into account not only the methanogenic microbial communities in soil environments, but also the surrounding catchment vegetation-derived inputs and hydrochemical conditions, as well as conduct more in-depth investigations into the mechanisms influencing methane production.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
KH: Software, Data curation, Visualization, Formal analysis, Writing – review & editing, Supervision, Writing – original draft, Methodology. JZ: Writing – original draft, Writing – review & editing, Resources, Methodology. JP: Writing – review & editing. QZ: Conceptualization, Writing – review & editing, Formal analysis. HF: Writing – review & editing, Conceptualization, Supervision, Funding acquisition.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the National Natural Science Foundation of China (U21A20186 and 32371592) and Natural Science Foundation of Gansu Province (23JRRA1034).
Acknowledgments
The authors thank the Core Facility of School of Life Sciences, State Key Laboratory of Grassland Agro-ecosystems of Lanzhou University and Supercomputing Center of Lanzhou University for samples and data analysis. We also give our sincere gratitude to Yongjun Liu, Chenyu Li, Yue Yang, Yong Zhang, Haohao Zhang and Weizhong Zhang for their help during the fieldwork and useful suggestions on the manuscript.
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 interest.
Generative AI statement
The author(s) declare that no Gen AI was used in the creation of this manuscript.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1616051/full#supplementary-material
References
1
AbdallaM.HastingsA.TruuJ.EspenbergM.ManderÜ.SmithP.et al. (2016). Emissions of methane from northern peatlands: a review of management impacts and implications for future management options. Ecol. Evol.6, 7080–7102. 10.1002/ece3.2469
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. 10.1038/ismej.2011.141
3
ArnoldW.TaylorM.BradfordM.RaymondP.PecciaJ. (2023). Microbial activity contributes to spatial heterogeneity of wetland methane fluxes. Microb. Spectr.11, e02714–23. 10.1128/spectrum.02714-23
4
BangerK.TianH.LuC. (2012). Do nitrogen fertilizers stimulate or inhibit methane emissions from rice fields?Glob. Change Biol.18, 3259–3267. 10.1111/j.1365-2486.2012.02762.x
5
BarberA.LeffJ. W.BradfordM. A.WallD. H.FiererN. (2014). Why are some microbes more ubiquitous than others?Predicting the habitat breadth of soil bacteria. Ecol. Lett.17, 794–802. 10.1111/ele.12282
6
BarberánA.BatesS. T.CasamayorE. O.FiererN. (2012). Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J.6, 343–351. 10.1038/ismej.2011.119
7
BrookerM. R.BohrerG.MouserP. J. (2014). Variations in potential CH4 flux and CO2 respiration from freshwater wetland sediments that differ by microsite location, depth and temperature. Ecol. Eng.72, 84–94. 10.1016/j.ecoleng.2014.05.028
8
BrownJ. H.GilloolyJ. F.AllenA. P.SavageV. M.WestG. B. (2004). Toward a metabolic theory of ecology. Ecology85, 1771–1789. 10.1890/03-9000
9
ChenH.ZhuQ.PengC.WuN.WangY.FangX.et al. (2013). The impacts of climate change and human activities on biogeochemical cycles on the Qinghai-Tibetan Plateau. Glob. Change Biol.19, 2940–2955. 10.1111/gcb.12277
10
ConradR. (2009). The global methane cycle: recent advances in understanding the microbial processes involved. Environ. Microbiol. Rep.1, 285–292. 10.1111/j.1758-2229.2009.00038.x
11
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. 10.1128/AEM.65.11.5066-5074.1999
12
CsardiG.NepuszT. (2005). The igraph software package for complex network research. Available online at: https://www.researchgate.net/publication/221995787 (Accessed May 21, 2025).
13
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. 10.1016/j.soilbio.2019.01.024
14
DingW-. X.CaiZ-. C. (2007). Methane emission from natural wetlands in China: summary of years 1995–2004 studies. Pedosphere17, 475–486. 10.1016/S1002-0160(07)60057-5
15
DongY.YuanJ.LiJ.LiuD.QiuY.ZhangX.et al. (2023). Conversion of natural coastal wetlands to mariculture ponds dramatically decreased methane production by reducing substrate availability. Agric. Ecosyst. Environ. 356:108646. 10.1016/j.agee.2023.108646
16
EilerA.HeinrichF.BertilssonS. (2012). Coherent dynamics and association networks among lake bacterioplankton taxa. ISME J.6, 330–342. 10.1038/ismej.2011.113
17
EtminanM.MyhreG.HighwoodE. J.ShineK. P. (2016). Radiative forcing of carbon dioxide, methane, and nitrous oxide: a significant revision of the methane radiative forcing. Geophys. Res. Lett.43, 12,614–12,623. 10.1002/2016GL071930
18
FanX.YinC.ChenH.YeM.ZhaoY.LiT.et al. (2019). The efficacy of 3,4-dimethylpyrazole phosphate on N2O emissions is linked to niche differentiation of ammonia oxidizing archaea and bacteria across four arable soils. Soil Biol. Biochem.130, 82–93. 10.1016/j.soilbio.2018.11.027
19
FangZ.PengZ.ZhangY.HanJ.PanJ.ZhangQ.et al. (2024). Drivers of denitrification and nitrification in a dryland agroecosystem: the role of abiotic and biotic factors. Appl. Soil Ecol.204:105691. 10.1016/j.apsoil.2024.105691
20
FaustK.SathirapongsasutiJ. F.IzardJ.SegataN.GeversD.RaesJ.et al. (2012). Microbial co-occurrence relationships in the human microbiome. PLoS Comput. Biol.8:e1002606. 10.1371/journal.pcbi.1002606
21
FengH.GuoJ.MalghaniS.HanM.CaoP.SunJ.et al. (2021). Effects of soil moisture and temperature on microbial regulation of methane fluxes in a poplar plantation. Forests12, 1–16. 10.3390/f12040407
22
Gruca-RokoszR.SzalD.BartoszekL.PekalaA. (2020). Isotopic evidence for vertical diversification of methane production pathways in freshwater sediments of Nielisz reservoir (Poland). Catena195:104803. 10.1016/j.catena.2020.104803
23
HaoW.ScharffeD.CrutzenP.SanhuezaE. (1988). Production of N2O, CH4, and CO2 from soils in the tropical savanna during the dry season. J. Atmos. Chem. 7, 93–105. 10.1007/BF00048256
24
HartmanW. H.Bueno De MesquitaC. P.TherouxS. M.Morgan-LangC.BaldocchiD. D.TringeS. G.et al. (2024). Multiple microbial guilds mediate soil methane cycling along a wetland salinity gradient. mSystems9, e00936–23. 10.1128/msystems.00936-23
25
HuH-. W.ZhangL-. M.YuanC-. L.HeJ-. Z. (2013). Contrasting Euryarchaeota communities between upland and paddy soils exhibited similar pH-impacted biogeographic patterns. Soil Biol. Biochem.64, 18–27. 10.1016/j.soilbio.2013.04.003
26
JanssonJ. K.HofmockelK. S. (2020). Soil microbiomes and climate change. Nat. Rev. Microbiol.18, 35–46. 10.1038/s41579-019-0265-7
27
JiaoS.LiuZ.LinY.YangJ.ChenW.WeiG. (2016). Bacterial communities in oil contaminated soils: biogeography and co-occurrence patterns. Soil Biol. Biochem.98, 64–73. 10.1016/j.soilbio.2016.04.005
28
JinZ.ShangJ.ZhuQ.LingC.XieW.QiangB.et al. (2020). “RFRSF: employee turnover prediction based on random forests and survival analysis,” in Web Information Systems Engineering – WISE 2020. WISE 2020. Lecture Notes in Computer Science, vol. 12343, eds. Z. Huang, W. Beek, H. Wang, R. Zhou, and Y. Zhang (Cham: Springer). 503–515. 10.1007/978-3-030-62008-0_35
29
JuottonenH. (2020). Disentangling the effects of methanogen community and environment on peatland greenhouse gas production by a reciprocal transplant experiment. Funct. Ecol.34, 1268–1279. 10.1111/1365-2435.13536
30
KembelS. W.CowanP. D.HelmusM. R.CornwellW. K.MorlonH.AckerlyD. D.et al. (2010). Picante: R tools for integrating phylogenies and ecology. Bioinformatics26, 1463–1464. 10.1093/bioinformatics/btq166
31
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. 10.1128/AEM.69.5.2423-2429.2003
32
LayeghifardM.HwangD. M.GuttmanD. S. (2017). Disentangling Interactions in the Microbiome: a network perspective. Trends Microbiol.25, 217–228. 10.1016/j.tim.2016.11.008
33
LiD.NiH.JiaoS.LuY.ZhouJ.SunB.et al. (2021). Coexistence patterns of soil methanogens are closely tied to methane generation and community assembly in rice paddies. Microbiome9, 1–13. 10.1186/s40168-020-00978-8
34
LutonP. E.WayneJ. M.SharpR. J.RileyP. W. (2002). The mcrA gene as an alternative to 16S rRNA in the phylogenetic analysis of methanogen populations in landfill. Microbiology148, 3521–3530. 10.1099/00221287-148-11-3521
35
MaB.WangH.DsouzaM.LouJ.HeY.DaiZ.et al. (2016). Geographic patterns of co-occurrence network topological features for soil microbiota at continental scale in eastern China. ISME J.10, 1891–1901. 10.1038/ismej.2015.261
36
MaK. E.QiuQ.LuY. (2010). Microbial mechanism for rice variety control on methane emission from rice field soil. Global Change Biol.16, 3085–3095. 10.1111/j.1365-2486.2009.02145.x
37
MoneckeA.LeischF. (2012). SemPLS: structural equation modeling using partial least squares. J. Statis. Softw.48, 1–32. 10.18637/jss.v048.i03
38
MorriënE.HannulaS. E.SnoekL. B.HelmsingN. R.ZweersH.de HollanderM.et al. (2017). Soil networks become more connected and take up more carbon as nature restoration progresses. Nat. Commun.8:14349. 10.1038/ncomms14349
39
MuC.MuM.WuX.JiaL.FanC.PengX.et al. (2023). High carbon emissions from thermokarst lakes and their determinants in the Tibet Plateau. Global Change Biol. 2732–2745. 10.1111/gcb.16658
40
MuM.MuC.LiuH.LeiP.GeY.ZhouZ.et al. (2025). Thermokarst lake drainage halves the temperature sensitivity of CH4 release on the Qinghai-Tibet Plateau. Nat. Commun.16:1992. 10.1038/s41467-025-57356-x
41
NazariesL.MurrellJ. C.MillardP.BaggsL.SinghB. K. (2013). Methane, microbes and models: fundamental understanding of the soil methane cycle for future predictions. Environ. Microbiol.15, 2395–2417. 10.1111/1462-2920.12149
42
OlesenJ. M.BascompteJ.DupontY. L.JordanoP. (2007). The modularity of pollination networks. Proc. Natl. Acad. Sci. U.S.A. 104, 19891–19896. 10.1073/pnas.0706375104
43
PengZ.JohnsonN. C.JansaJ.HanJ.FangZ.ZhangY.et al. (2024). Mycorrhizal effects on crop yield and soil ecosystem functions in a long-term tillage and fertilization experiment. New Phytol.242, 1798–1813. 10.1111/nph.19493
44
QinL.NiB.ZouY.FreemanC.PengX.YangL.et al. (2024). Deciphering soil environmental regulation on reassembly of the soil bacterial community during wetland restoration. Sci. Total Environ.954:176586. 10.1016/j.scitotenv.2024.176586
45
QinS.KouD.MaoC.ChenY.ChenL.YangY.et al. (2021). Temperature sensitivity of permafrost carbon release mediated by mineral and microbial properties. Sci. Adv.7:eabe3596. 10.1126/sciadv.abe3596
46
QiuS.ZhangX.XiaW.LiZ.WangL.ChenZ.et al. (2023). Effect of extreme pH conditions on methanogenesis: methanogen metabolism and community structure. Sci. Total Environ.877:162702. 10.1016/j.scitotenv.2023.162702
47
RenC.ZhangX.ZhangS.WangJ.XuM.GuoY.et al. (2021). Altered microbial CAZyme families indicated dead biomass decomposition following afforestation. Soil Biol. Biochem.160:108362. 10.1016/j.soilbio.2021.108362
48
RosentreterJ. A.BorgesA. V.DeemerB. R.HolgersonM. A.LiuS.SongC.et al. (2021). Half of global methane emissions come from highly variable aquatic ecosystem sources. Nat. Geosci. 14, 225–230. 10.1038/s41561-021-00715-2
49
SaunoisM.StavertA. R.PoulterB.BousquetP.CanadellJ. G.JacksonR. B.et al. (2020). The global methane budget 2000–2017. Earth Syst. Sci. Data12, 1561–1623. 10.5194/essd-12-1561-2020
50
SeppeyC. V. W.CabrolL.ThalassoF.GandoisL.LavergneC.Martinez-CruzK.et al. (2023). Biogeography of microbial communities in high-latitude ecosystems: contrasting drivers for methanogens, methanotrophs and global prokaryotes. Environ. Microbiol. 25, 3364–3386. 10.1111/1462-2920.16526
51
StavertA. R.SaunoisM.CanadellJ. G.PoulterB.JacksonR. B.RegnierP.et al. (2022). Regional trends and drivers of the global methane budget. Global Change Biol.28, 182–200. 10.1111/gcb.15901
52
TavellaJ.CagnoloL. (2019). Does fire disturbance affect ant community structure? Insights from spatial co-occurrence networks. Oecologia189, 475–486. 10.1007/s00442-018-4320-2
53
VizzaC.WestW. E.JonesS. E.HartJ. A.LambertiG. A. (2017). Regulators of coastal wetland methane production and responses to simulated global change. Biogeosciences14, 431–446. 10.5194/bg-14-431-2017
54
Von EggersJ. M.WisnoskiN. I.CalderJ. W.CapoE.GroffD. V.KristA. C.et al. (2024). Environmental filtering governs consistent vertical zonation in sedimentary microbial communities across disconnected mountain lakes. Environ. Microbiol.26:e16607. 10.1111/1462-2920.16607
55
WaggC.BenderS. F.WidmerF. (2014). Soil biodiversity and soil community composition determine ecosystem multifunctionality. PNAS111, 5266–5270. 10.1073/pnas.1320054111
56
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–13. 10.1128/mBio.00592-13
57
WikM.VarnerR. K.AnthonyK. W.MacIntyreS.BastvikenD. (2016). Climate-sensitive northern lakes and ponds are critical components of methane release. Nat. Geosci.9, 99–105. 10.1038/ngeo2578
58
XiH.JiaM.KuzyakovY.PengZ.ZhangY.HanJ.et al. (2023). Key decomposers of straw depending on tillage and fertilization. Agric. Ecosyst. Environ.358:108717. 10.1016/j.agee.2023.108717
59
XiangJ.LiuD.DingW.YuanJ.LinY. (2015). Invasion chronosequence of Spartina alterniflora on methane emission and organic carbon sequestration in a coastal salt marsh. Atmos. Environ.112, 72–80. 10.1016/j.atmosenv.2015.04.035
60
YangK.LuoS.HuL.ChenB.XieZ.MaB.et al. (2020). Responses of soil ammonia-oxidizing bacteria and archaea diversity to N, P and NP fertilization: relationships with soil environmental variables and plant community diversity. Soil Biol. Biochem.145:107795. 10.1016/j.soilbio.2020.107795
61
YangP.TangK. W.TongC.LaiD. Y. F.WuL.YangH.et al. (2022). Changes in sediment methanogenic archaea community structure and methane production potential following conversion of coastal marsh to aquaculture ponds. Environ. Pollut.305:119276. 10.1016/j.envpol.2022.119276
62
YangR.JiM.ZhangX.HeF.YuZ.ZengJ.et al. (2024). Methane emissions and microbial communities under differing flooding conditions and seasons in littoral wetlands of urban lake. Environ. Res. 250:118390. 10.1016/j.envres.2024.118390
63
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. 10.1016/j.mimet.2014.05.006
64
YavittJ. B.Seidman-ZagerM. (2006). Methanogenic Conditions in Northern Peat Soils. Geomicrobiol. J.23, 119–127. 10.1080/01490450500533957
65
ZhangC.LeiS.WuH.LiaoL.WangX.ZhangL.et al. (2024). Simplified microbial network reduced microbial structure stability and soil functionality in alpine grassland along a natural aridity gradient. Soil Biol. Biochem.191:109366. 10.1016/j.soilbio.2024.109366
66
ZhouJ.TherouxS. M.Bueno De MesquitaC. P.HartmanW. H.TianY.TringeS. G.et al. (2022). Microbial drivers of methane emissions from unrestored industrial salt ponds. ISME J.16, 284–295. 10.1038/s41396-021-01067-w
67
ZhouJ.XiangJ.WangL.ZhongG.ZhuG.LiuC.et al. (2019). The impacts of groundwater chemistry on wetland vegetation distribution in the northern Qinghai-Tibet Plateau. Sustainability11, 1–15. 10.3390/su11185022
68
ZhuangG-. C.EllingF. J.NigroL. M.SamarkinV.JoyeS. B.TeskeA.et al. (2016). Multiple evidence for methylotrophic methanogenesis as the dominant methanogenic pathway in hypersaline sediments from the Orca Basin, Gulf of Mexico. Geochim. Cosmochim. Acta187, 1–20. 10.1016/j.gca.2016.05.005
69
ZieglerM.EguíluzV. M.DuarteC. M.VoolstraC. R. (2018). Rare symbionts may contribute to the resilience of coral-algal assemblages. ISME J.12, 161–172. 10.1038/ismej.2017.151
Summary
Keywords
methanogens, coexistence network complexity, CH4 emission, wetlands, mcrA gene
Citation
He K, Zhao J, Pan J, Zhang Q and Feng H (2025) Coexistence networks of soil methanogens are closely tied to methane generation in wetlands on the northeastern of the Qinghai–Tibet Plateau. Front. Microbiol. 16:1616051. doi: 10.3389/fmicb.2025.1616051
Received
22 April 2025
Accepted
13 June 2025
Published
01 July 2025
Volume
16 - 2025
Edited by
Shuai Liu, Qingdao University of Science and Technology, China
Reviewed by
Ying Wang, Anhui Normal University, China
Dejie Li, Shandong Normal University, China
Updates
Copyright
© 2025 He, Zhao, Pan, Zhang and Feng.
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: Huyuan Feng fenghy@lzu.edu.cn
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.