Xiaorong Zhou1,2Xianke Chen1,3,4Xiangning Qi1,2Yiyuan Zeng1,2
Xiaowei Guo5
Guoqiang Zhuang1,2
Anzhou Ma1,2*- 1Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing, China
- 2College of Resources and Environment, University of Chinese Academy of Sciences, Beijing, China
- 3Sino-Danish College of University of Chinese Academy of Sciences, Beijing, China
- 4Sino-Danish Center for Education and Research, Beijing, China
- 5Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
Introduction: The functions of terrestrial ecosystems are mainly maintained by bacteria, as a key component of microorganisms, which actively participate in the nutrient cycling of ecosystems. Currently, there are few studies have been carried out on the bacteria contributing to the soil multi-nutrient cycling in responding to climate warming, which hampers our obtainment of a comprehensive understanding of the ecological function of ecosystems as a whole.
Methods: In this study, the main bacteria taxa contributing to the soil multi-nutrient cycling under the long-term warming in an alpine meadow was determined based onphysichemical properties measurement and high-throughput sequencing, and the potential reasons that warming altered the main bacteria contributing to the soil multi-nutrient cycling were further analyzed.
Results: The results confirmed that the bacterial β-diversity was crucial to the soil multi-nutrient cycling. Furthermore, Gemmatimonadetes, Actinobacteria, and Proteobacteria were the main contributors to the soil multi-nutrient cycling, and played pivotal roles as keystone nodes and biomarkers throughout the entire soil profile. This suggested that warming altered and shifted the main bacteria contributing to the soil multi-nutrient cycling toward keystone taxa.
Discussion: Meanwhile, their relative abundance was higher, which could make them have the advantage of seizing resources in the face of environmental pressures. In summary, the results demonstrated the crucial role of keystone bacteria in the multi-nutrient cycling under the climate warming in the alpine meadow. This has important implications for understanding and exploring the multi-nutrient cycling of alpine ecosystems under the global climate warming.
1. Introduction
Ecosystems maintain multiple functions simultaneously. Nutrient cycling is an important part of ecosystem function. The index of multi-nutrient cycling defined by humans, which is not a single measurable index that can quantify more than one process of the ecosystem (Byrnes et al., 2014). Microbes are the engines driving the biogeochemical cycle of the earth, especially in the terrestrial ecosystems (Falkowski et al., 2008), and play a vital role in the function of ecosystems. Numerous studies have focused on the role of microorganisms in the multifunctionality of ecosystems (Wagg et al., 2014; Jiao et al., 2019, 2021). It has been found that the biodiversity of microbes clearly influenced the ecosystem multifunctionality (Jing et al., 2015). The multifunctionality was driven by soil microbial diversity, which has been demonstrated by previous studies in terrestrial ecosystems (Jing et al., 2015; Delgado-Baquerizo et al., 2016). In addition, species themselves also affect the function of ecosystems. The functional importance of a single species does not directly translate into the functional importance of the microorganisms in the community as a whole (Pasari et al., 2013). This is because there are interactions between species. The research demonstrated that the positive and negative associations between species were tightly correlated with the ecosystem function (Jiao et al., 2021). Thus, it is necessary to pay attention to the role of microorganisms in the microbial network, especially the keystone nodes in the network which may have a more dramatic impact on the ecosystem multifunctionality. Soil microbial communities account for the majority of belowground biodiversity, particularly bacterial communities, make significant contributions to ecosystem biomass, biodiversity, element cycling, and energy flow (Nielsen et al., 2015; Jiao et al., 2019). However, little attention has been paid to the impact of bacterial communities on the multi-nutrient cycling under climate change. Due to the importance of multifunctionality, it is necessary to determine the contribution of bacteria to the soil multi-nutrient cycling in terrestrial ecosystems.
In recent years, evidence has emerged that terrestrial ecosystems have a positive feedback on climate warming (Heimann and Reichstein, 2008). Specifically, climate warming effected the activities of soil microorganisms directly and indirectly, and greenhouse gases (e.g., carbon dioxide, methane) produced by microbial activities were released into the atmosphere, which made a contribution to global warming (Petchey et al., 1999; Bardgett et al., 2008; Vincent, 2010; Dove et al., 2021). The soil multi-nutrient cycling mainly mediated by microbes, so the effect of climate warming on microbiomes will ultimately affect soil multi-nutrient cycling. The bacterial diversity has responded differently to climate warming. It was found that bacterial α-diversity in subtropical planted forest topsoil (0–10 cm; Wang et al., 2019) and alpine meadow topsoil (0–20 cm; Li et al., 2016) increased under climate warming. However, there was no significant change in bacterial α-diversity in the topsoil (0–10 cm) of the forest-alpine cross zone after warming (Zheng et al., 2020). Thus, differences in the α-diversity of bacteria responding to warming may affect the contribution of bacteria populations to soil multi-nutrient cycling in the ecosystem, owing to previous studies that have proved that the diversity of bacteria was tightly correlated with the multi-nutrient cycling. Warming also affected the bacterial population changes besides the bacterial diversity (Guo et al., 2018; Dove et al., 2021; Zhou et al., 2021). The composition of soil bacteria have been altered by warming (DeAngelis et al., 2015; Cheng et al., 2017). For instance, there has been a decline in the relative abundance of Planctomycetes and Actinobacteria (Zhou et al., 2021). Climate warming enhanced the interactions among species, owing to the increased complexity and stability of microbial networks (Yuan et al., 2021). Connectors in microbial networks were more connected to other modules (Guimerà and Nunes Amaral, 2005). Therefore, these species were more resistant to environmental stresses. However, few researches focused on the changes in bacterial community composition affecting the roles of bacteria in soil multi-nutrient cycling. Furthermore, bacterial communities are present throughout soil profile. There is little further attention has been paid to the impact of population changes on the multi-nutrient cycling under warming across the entire soil profile. It remains unclear which bacterial groups contribute to the multi-nutrient cycling, particularly under long-term warming.
More than 85% of the Qinghai-Tibetan Plateau is covered by alpine meadows, shrub-meadows, and alpine steppes, and the rate of increase in its temperature is approximately three times that of global warming (Qiu, 2008). Particularly, the alpine meadow ecosystem is extremely sensitive and vulnerable to climate change, which is very important for the dynamic changes in the element cycling (Chen et al., 2013). Thus, the Qinghai-Tibetan Plateau is a typical experimental area for exploring the response and feedback of subsurface bacterial communities to climate warming (Heimann and Reichstein, 2008). It is very necessary to explore the effects of climate warming on bacteria involved in soil multi-nutrient cycling in alpine ecosystems that are extremely sensitive to climate warming. The aim of this study was to identify the potential association between the bacterial community and soil multi-nutrient cycling, and to investigate whether warming affected the main bacteria contributors to the multi-nutrient cycling in the alpine meadow. The importance of bacterial diversity to the soil multi-nutrient cycling was first explored. Then, the main bacteria involved in the multi-nutrient cycling were identified under control and warming. This study deepened our understanding of the bacteria involved in soil multi-nutrient cycling under climate warming, especially in alpine ecosystems.
2. Materials and methods
2.1. Study site and sampling
This study was conducted at the Haibei National Alpine Grassland Ecosystem Research Station (37°36′38.53″ N, 101°18′49.31″ E), in northwestern China at an average elevation of 3,195 m. During the summer, the station area has a southeast monsoon climate, and in winter, the northwest Asian cold effects it. The warm seasons are short and cool, and the cold seasons are long and cold. The annual maximum and minimum extreme temperatures are 27.6°C and −37.1°C, respectively. The annual average temperature is −1.7°C, and the annual precipitation is 426–860 mm, of which 80% is distributed in the plant growth period from May to September. The parent material of the soil is loess with alluvial sediments underneath. The open top chambers (OTCs) were often used to explore the responses of ecosystems to climate warming (Marion et al., 1997). In this study, these OTCs were 1.5 m in diameter and 40 cm in height. The temperature of the OTCs was measured by online monitoring device. During the growing season, OTCs increased the average daily temperature by 0.6–2.0°C, and more information about OTCs referred to the reference (Klein et al., 2004). In August 2021, the depth of the soil core collected in the study area was 0–80 cm. Nine soil cores (core diameter of 38 mm) were randomly collected at one time in the warming (soil in OTCs) and control (soil outside OTCs) plots using a soil sampling device, respectively (Supplementary Figure S1). Meanwhile, the number of soil samples replicates is nine. The cores were manually segmented as follows: 0–15 cm, 15–30 cm, 30–45 cm, 45–60 cm, and 60–80 cm (corresponding to L1, L2, L3, L4, and L5, respectively). Consequently, there were 90 (2 × 9 × 5) samples in total. After the soil samples were transported to the laboratory, some were stored at 4°C to measure the physicochemical properties, while other parts were stored at −20°C for other experiments.
2.2. The assessment of multi-nutrient cycling index
The contents of ammonium nitrogen (NH4+-N) and nitrate nitrogen (NO3−-N) were measured using an AA3 continuous flow analytical system (Seal, German). The soil organic carbon (SOC) concentration was determined using the potassium dichromate oxidation method (Walkley and Black, 1934). The soil dissolved organic carbon (DOC) was extracted by adding 10 mL of water to 1 g of soil, shaking for 24 h at 30°C and filtering through a 0.45 μm filter (Millipore). The sample was analyzed using an elementary total organic carbon (TOC) analyzer (Vario TOC; Elementar, Germany). The total carbon (TC) and total nitrogen (TN) content were measured using an elementary analyzer (Vario MACRO Cube, Elementar, Germany), while the total phosphorus (TP) content was measured using the acid fusion-Mo-Sb anti spectrophotometric method. The available nitrogen (AN) concentration was determined using the alkali diffusion method. Finally, the available phosphorus (AP) was extracted using sodium bicarbonate and was measured using the molybdenum-blue method (Huang et al., 2021).
The concept of soil multi-nutrient cycling has been applied in many studies (Pasari et al., 2013; Delgado-Baquerizo et al., 2016; Gross et al., 2017; Le Bagousse-Pinguet et al., 2019; Zheng et al., 2019). First, we standardized each of the nine measured functions (NH4+-N, NO3−-N, DOC, TN, TC, TP, SOC, AN, and AP) using the Z-score transformation, and then, we averaged them to obtain an index of the multi-nutrient cycling for each site (Delgado-Baquerizo et al., 2016). We chose these functions because they provided some of the essential functions needed to support and regulate ecosystems.
2.3. DNA extraction and sequencing
The total deoxyribonucleic acid (DNA) in the soil was extracted from 0.5 g of thoroughly mixed soil using a MOBIO PowerSoil DNA isolation kit (MOBIO Laboratories, Carlsbad, CA, United States) according to the instructions of the kit manufacturer. The integrity and purity of the DNA were monitored on a 1% agarose gel. The V3–V4 hypervariable regions of 16S ribosomal ribonucleic acid (rRNA) genes were PCR amplified with primers 338F (5′-ACTCCTACGGGAGGCAGCA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). Polymerase chain reactions (PCRs) were conducted in a 50 μL reaction system, containing 25 μL 2Χ Premix Taq (TaKaRa), 1 μL each of 10 μM forward and reverse primers, and 50 ng of the DNA. The thermal cycle conditions were as follows: 5 min at 94°C for initialization; followed by 30 cycles of 30 s of denaturation at 94°C, 30 s of annealing at 52°C, and 30 s of extension at 72°C, 10 min final elongation at 72°C, and then held at 4°C. The PCR products were tested using 1% agarose gels and purified using an E.Z.N.A. Gel Extraction Kit (Omega). The concentrations were measured using NanoDrop 2000 and Qubit 3.0 (Thermo Fisher Scientific, Waltham, United States) and were stored at −80°C. The DNA samples were then sequenced using an Illumina Nova6000 platform and 250 bp paired-end reads were generated (Guangdong Magigene Biotechnology Co., Ltd. Guangzhou, China).
2.4. Bioinformatics and statistical analyses
A pipeline was used to generate an operational taxonomic unit (OTU) table1 (Feng et al., 2017). Wilcoxon test was used to compare whether α-diversity was significantly different between treatments. Principal coordinates analysis (PCoA) was conducted to assess the community similarity among the depth intervals. Significant differences in the composition of the community between depth intervals were assessed using permutation multivariate analysis of variance (PERMANOVA), with 999 permutations calculated per test.
After analysis by the online pipeline,2 phylogenetic molecular ecological network (pMEN) graphs were created using the Gephi 0.9.2 software program (Zhou et al., 2010; Deng et al., 2012). The OTUs appearing in more than half of the samples were retained, and the RMT threshold were consistent (0.830). Then the Spearman correlation matrix were used to construct the microbial network. The roles of the network nodes were identified based on their values of within-module connectivity (Zi) values and among-module connectivity (Pi) values (Guimerà and Nunes Amaral, 2005). Furthermore, the Levins’ niche breadth index was estimated according to the reference (Jiao et al., 2021).
In addition, the random forest (RF) algorithm was used to identify the taxa that represented the most under warming and control using the rfPermute package. The A3 package was applied to calculate the significance level and cross-validated R2 value of the model with 5,000 permutations of the response variable (Jiao et al., 2020).
All of the statistical analyses were carried out in the R environment (v4.0.4; https://www.r-project.org/), using vegan, ggplot2, and spaa.
3. Results
3.1. The diversity of bacteria correlated with nutrient variables
There was a total of 11,463,321 reads were obtained, and they were grouped into 34,192 OTUs with 97% similarity. The species accumulation curves were close to the asymptotes, suggesting that the sequencing recovered much of the local species diversity (Supplementary Figure S2). First, warming had no significant effect on the α-diversity of the bacterial community (Figures 1A–D). The Shannon index of bacteria in L1 increased after warming, indicating that the diversity of the community increased (Supplementary Figure S3D). Then, the analysis of PCoA was used to explore the overall variability of the bacteria composition (Figures 1E,F), which that about 54% of the variance could be explained by depth, indicating that warming significantly increased the differences among the bacterial communities in the soil profile. Meanwhile, warming significantly increased the differences between bacterial communities in L1, L2, and L3 layers (Supplementary Figures S4A–C).
Figure 1. The bacterial diversity throughout the entire soil profile under control and warming. The α-diversity of bacteria between control and warming (using the Mann–Whitney U-test) (A–D). PCoA analysis was performed on differences in community composition at different depths under control (E) and warming (F).
To explore the relationship between bacterial diversity and the individual nutrient variables, the results of their correlation coefficients were shown by a correlation heatmap (Figures 2A,B). Both α-diversity (PD, Chao 1, Richness, Shannon) and β-diversity (PC1) of bacteria were significantly correlated with all individual variables of nutrient in the whole soil profile (Figures 2A,B), indicating that both diversity and community composition of bacteria play crucial roles in soil multi-nutrient cycling. Most indexes of α-diversity were significantly positively correlated with individual nutrient variables except AP, while the β-diversity was significantly negatively correlated. Warming appeared to have no effect on the correlation of bacterial diversity with individual nutrient variables.
Figure 2. The correlation between nutrient variables and the diversity of bacteria in the alpine ecosystem. Heatmaps of the Spearman’s correlation coefficients between community diversity and individual nutrient variables under control (A) and warming (B). *p < 0.05; **p < 0.01. The importance of bacterial α-and β-diversity driving the cycling of soil multi-nutrient in the entire profile under control (C) and warming (D), was assessed by RF analysis.
Bacterial diversity was significantly correlated with multiple soil nutrient variables (Figures 2A,B), but the ecosystems maintain multiple functions at the same time, the index of multi-nutrient cycling could better reflect the functions of terrestrial ecosystems. Therefore, it is of great significance to explore the importance of bacterial diversity to the cycling of soil multi-nutrient. Based on the random forest analysis, the main bacterial predictors of the soil multi-nutrient cycling index were identified, which elucidated the major factors driving the soil multi-nutrient cycling. The results found that most essential variable for predicting multi-nutrient cycling across the entire soil profile was the β-diversity, followed by the α-diversity (Figures 2C,D). It was worth noting that the bacterial diversity index related to the change in soil multi-nutrient cycling index differed after warming, although this difference was little obvious.
3.2. Potential bacteria contribute to the soil multi-nutrient cycling in the alpine meadow
The biological contribution of the bacteria to the soil properties was further assessed at the phylum level (dominant phyla) based on the RF analysis (Figure 3; Supplementary Figure S5). It is evident that rather than all bacterial phyla made the same contribution to various soil variables. Specifically, the most important variable to predict the soil properties was Planctomycetes (TC, SOC, TN, and AN), indicating that they were crucial to soil nutrient cycling (p < 0.001). In addition, other important variables for predicting the soil properties were the Gemmatimonadetes for TP (p = 0.001), and the Verrucomicrobia for DOC (p < 0.001; Figure 3). In addition, after warming, Gemmatimonadetes was the most important to predict the soil properties (DOC, SOC, TN, and AN), followed by Actinobacteria. Proteobacteria was a more important variable for predicting TC and AP under warming (p < 0.001). Planctomycetes and Verrucomicrobia played potential roles in multi-nutrient cycling in the alpine meadow (Supplementary Figure S5). However, these roles were filled by Gemmatimonadetes, Actinobacteria, and Proteobacteria after warming (Figure 3), suggesting that warming shaped the main bacteria contributing to the soil multi-nutrient cycling throughout the soil entire profile.
Figure 3. Main potential bacteria contribute to soil multi-nutrient cycling in the alpine ecosystem under warming. The importance of these predictors was estimated using the percent increase in the MSE of the variables, with higher MSE% values indicating more important predictors. The p < 0.05 (*) and p < 0.01 (**) indicate significance at the levels of 5% and 1%, respectively.
3.3. Species coexistence of bacteria and their biomarkers in the entire soil profile
Phylogenetic molecular ecological network (pMEN) often used to study the interaction between species (Deng et al., 2012). The coexistence of species in soil was investigated by constructing microbial co-occurrence networks (Figure 4A), owing to the existence of bacteria throughout the entire profile. The number of nodes for control and warming were 120 and 192, respectively; and the number of edges were 431 and 831. In addition, the total nodes, links, and average degree increased after warming (Supplementary Table S1), indicating that the bacteria network became more complex and connected. The proportion of negative interactions also increased (Supplementary Table S1). Due to environmental stress caused by warming, competitive interactions between species has increased, which may lead to a decline in community diversity. Furthermore, the roles of bacterial taxa were also determined based on their Zi and Pi values in the bacteria networks. The proportion and number of keystone nodes (i.e., connectors, network, and module hubs) increased after warming, and were associated with Acidobacteria, Actinobacteria, Bacteroidetes, Proteobacteria (Figure 4B; Supplementary Table S2).
Figure 4. The bacteria network and biomarkers under control and warming. (A) Microbial co-occurrence networks under control and warming. The thickness of the edges is proportional to the value of the Spearman’s correlation coefficient, and the size of the nodes is proportional to the degree of the OTUs. (B) Roles of nodes based on their Zi and Pi values in the co-occurrence networks.
Generally, there are some microorganisms that are particularly sensitive to changes in the environment, which could be regarded as biomarkers. The biomarkers were identified throughout the whole soil profile based on the analysis of RF. There were 29 OTUs were biomarkers, and most of them belonged to Acidobacteria, Actinobacteria, Bacteroidetes, and Proteobacteria (Figure 5). In addition, Verrucomicrobia, Planctomycetes, Firmicutes, Gemmatimonas, Roseisolibacter, Dissulfurirhabdus, Solimonas, Sulfuricaulis, and Thiobacter could be considered to be biomarkers of the entire soil profile after warming (Supplementary Table S3). Moreover, some specific biomarkers (OTUs) were also found to play important roles in the network under warming. Some bacteria in Acidobacteria (OTU_104, OTU_1760, OTU_256, and OTU_36), Gemmatimonadetes (OTU_17 and OTU_20), Actinobacteria (OTU_598), Proteobacteria (OTU_1, OTU_43 and OTU_51), and Verrucomicrobia (OTU_131) were the key connectors and module hubs in the bacterial network under warming (Supplementary Tables S2, S3). Therefore, these OTUs with dual identities could be regarded as species that are more sensitive and even more strongly affected by climate warming.
Figure 5. Biomarkers of bacteria between control and warming throughout the entire profile. The top 29 biomarker bacterial OTUs were accurately identified by RF analysis in the control (A), and 32 were identified in the warming profile (B).
4. Discussion
Microorganisms actively participate in the cycling of various nutrient elements, which is of great significance to the soil multi-nutrient cycling in terrestrial ecosystems. The diversity of microorganisms is closely related to the nutrient variables, whether they are single variables or multi-nutrient cycling index. Multiple elemental cycles occur simultaneously in the soil ecosystems. Therefore, the index of multi-nutrient cycling was used to explore the functions of terrestrial ecosystems. Soil microbial communities account for the majority of belowground biodiversity, particularly bacterial communities. Therefore, the roles of bacteria in the soil multi-nutrient cycling were first determined by correlating the bacterial diversity with the multi-nutrient cycling. The results indicated that the bacterial β-diversity had a stronger effect on the soil multi-nutrient cycling throughout the entire soil profile. Most importantly, Gemmatimonadetes, Actinobacteria, and Proteobacteria, which were identified as biomarkers and keystone nodes, were important contributors to the cycling of multi-nutrient in the soil profile after warming. Long-term warming has resulted in keystone taxa being the main bacteria involved in the soil multi-nutrient cycling. This study shed light on the contribution of bacteria to soil multi-nutrient cycling, especially provided insights into the impact of changes in bacterial populations on the soil multi-nutrient cycling in the context of a warming climate.
4.1. Main bacteria contributing to the soil multi-nutrient cycling were altered in the whole soil profile after warming
Bacteria are important components in driving biogeochemical cycles in terrestrial ecosystems (Wagg et al., 2014), such as the degradation of organic matter (Xue et al., 2016), and nitrogen cycling (Zhou et al., 2015; Li et al., 2019). The microbial diversity is tightly correlated with the soil multi-nutrient cycling in terrestrial ecosystems (Delgado-Baquerizo et al., 2016; Zheng et al., 2019). In this study, it was found that the bacterial β-diversity was the most important in predicting the soil multi-nutrient cycling index across the entire soil profile (Figures 2C,D). This was supported by the results of a previous study, which reported that the β-diversity was more important when multiple functions were considered (Pasari et al., 2013). Moreover, %IncMSE (percent increase of the mean square error) in PC1 increased from 21.84% to 23.20% after warming (Figures 2C,D), which was consistent with the fact that differences between the bacterial communities along the soil profile increased after warming (Figures 1E,F). The bacterial diversity associated with changes in the soil multi-nutrient cycling index differed after warming (Figures 2C,D), which indicated the importance of bacterial diversity on soil multi-nutrient cycling differed after warming.
Warming altered the bacteria involved in the multi-nutrient cycling, although there was seemed to have less effect on the proportion of bacteria in the whole profile (Supplementary Figure S6A). Planctomycetes and Verrucomicrobia were important contributors to the soil multi-nutrient cycling, but Gemmatimonadetes, Actinobacteria and Proteobacteria became the important contributors after warming (Figure 3; Supplementary Figure S5). Studies have demonstrated that the warming accelerated lignin decomposition caused by Proteobacteria (Tao et al., 2020), which supported the results in this study that Proteobacteria played a crucial role in the multi-nutrient cycling under warming. Actinobacteria exhibited an increasing trend under climate warming (Rui et al., 2015), which meant Actinobacteria could be better able to adapt to warming. In addition, Gemmatimonadetes, as dominant phylum, have a stronger tolerance to environmental stress.
4.2. Species coexistence of bacteria and their biomarkers in the whole soil profile
The activities and interspecific interactions of microbiomes affected the soil nutrient cycling and many ecosystem processes (Cardinale et al., 2012), so the role of bacteria was explored in multi-nutrient cycling based on the coexistence of species. Warming increased the complexity of bacterial networks, including the network size, and average degree (Supplementary Table S2), which was consistent with a study in a tallgrass prairie ecosystem (Yuan et al., 2021). Previous studies have found that climate warming has enhanced network stability over time (Yuan et al., 2021). The higher complexity (larger network, higher connectivity) can render the system more resistant (Okuyama and Holland, 2008; Landi et al., 2018) but less resilient (Boucot, 1985). Therefore, it could be speculated that the microorganisms have been adapted to changes in the environment after long-term warming. However, Hernandez et al. found that the increasing stress through anthropogenic perturbations could destabilize microbiomes (Hernandez et al., 2021), suggesting that climate warming has less impact on microbial networks than anthropogenic disturbances. Warming could select some fast-growing bacteria as a deterministic filter (Guo et al., 2018) based on the metabolic theory of ecology (MTE; Brown et al., 2004). The bacteria developed thermal adaptation (Bradford, 2013) to this environmental stress during the long-term warming in this study. This further confirmed the previous conclusion that the microbial community could adapt to the changes in the environment (Walther et al., 2002). The proportion of negative interactions increased after warming, which may lead to a decrease in community diversity or even reduced the multi-nutrient cycling of soil (Delgado-Baquerizo et al., 2016). Most of the keystone nodes were abundant taxa [OTUs with relative abundance above 0.1% of the total sequences (Liu et al., 2015)], and the members of the Betaproteobacteria, Thermoleophilia and Acidobacteria Gp4 classes were the most prominent keystone nodes under warming (Supplementary Table S3), accounting for 52.3% of all of the module hub connectors and network hubs. Some of the biomarkers were keystone nodes in the network under warming (Supplementary Tables S2, S3). It may be the hub species in the network that initially respond to the changes in the environment, sensing and then passing on this information to other species in the population.
Gemmatimonadetes, Actinobacteria, and Proteobacteria were important contributors to the cycling of soil multi-nutrient after warming. Based on this, their roles were further explored in microbial networks and biomarkers. The proportion of Gemmatimonadetes, Actinobacteria, and Proteobacteria in the keystone nodes reached 63.63% under warming (Supplementary Table S2), and their proportion also reached 56.25% in biomarkers after warming (Supplementary Table S3). These results suggested that the main contributors (Gemmatimonadetes, Actinobacteria, and Proteobacteria) serve dual role under warming. They may have been the first to perceive and even adapt to the changes in the environment. In addition, they transmitted this change to their partners through interactions as keystone nodes in the network. Finally, they affected the multi-nutrient cycling of the terrestrial ecosystems.
4.3. Bacterial phyla with high relative abundance after warming were the main contributors to multi-nutrient cycling
The potential reasons for the change in bacteria participating in the multi-nutrient cycling after warming were further analyzed from the perspective of relative abundance (Supplementary Figures S6A,B). Planctomycetes and Verrucomicrobia were important contributors to the soil multi-nutrient cycling, and their relative abundance were 0.56% and 2.14%, respectively (Supplementary Figure S6A). The main bacteria that contributed to the soil multi-nutrient cycling became Gemmatimonadetes, Actinobacteria, and Proteobacteria after warming, with the relative abundance 12.83%, 7.06%, and 32.88%, respectively (Supplementary Figure S6A). This indicated that the relative abundance of the three phyla was higher. In order to eliminate the difference in soil heterogeneity caused by depth, the relative abundance of bacteria in each layer was studied (Supplementary Figure S6B). The relative abundance of Gemmatimonadetes, Actinobacteria, and Proteobacteria was higher in each layer. Therefore, the phyla with lower relative abundance are more stressed by warming, which inhibited their functions in soil multi-nutrient cycling, while the phyla with higher relative abundance were more adaptable and better copied with the environmental stress caused by warming. In addition, the community-level niche breadths (Bcom) of bacteria in the whole profile showed an increasing trend after warming (Supplementary Figure S6C). It could be explained that the species with high relative abundance have more advantages to occupy resources than the species with low relative abundance, with the increase of niche breadth under warming.
In conclusion, this study filled a gap in the impact of climate warming on bacteria involved in the multi-nutrient cycling. In this study, the diversity and contribution of bacteria were explored to the soil multi-nutrient cycling in an alpine meadow under long-term warming. The results indicated that the bacterial β-diversity played a dominant role in the multi-nutrient cycling in the alpine meadow. In addition, climate warming altered and resulted in keystone taxa being the major bacteria contributing to the soil multi-nutrient cycling. In summary, this study deepens our understanding of bacteria involved in soil multi-nutrient cycling in alpine meadow ecosystems, and shows that keystone bacteria are more essential for the soil multi-nutrient cycling under climate warming. Future studies could further focus on whether metabolic exchange between keystone bacteria affects their roles in soil multi-nutrient cycling under environmental stress of warming.
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
XZ, GZ, and AM: conceptualization and design. XZ, XC, and XQ: methodology, data collection and formal analysis. XZ, YZ, and XG: validation, investigation, and resources. XZ: writing the original manuscript. GZ and AM: supervision, reviewing and editing the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Second Tibetan Plateau Scientific Expedition and Research Program (2019QZKK0402 and 2019QZKK0307), the National Key Research and Development Program of China (2018YFA0901200), and the National Natural Science Foundation of China (41671270 and 41673082).
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.
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.2023.1136187/full#supplementary-material
Footnotes
References
Bardgett, R. D., Freeman, C., and Ostle, N. J. (2008). Microbial contributions to climate change through carbon cycle feedbacks. ISME J. 2, 805–814. doi: 10.1038/ismej.2008.58
Boucot, A. J. (1985). The complexity and stability of ecosystems. Nature 315, 635–636. doi: 10.1038/315635c0
Bradford, M. (2013). Thermal adaptation of decomposer communities in warming soils. Front. Microbiol. 4:333. doi: 10.3389/fmicb.2013.00333
Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M., and West, G. B. (2004). Toward a metabolic theory of ecology. Ecology 85, 1771–1789. doi: 10.1890/03-9000
Byrnes, J. E. K., Gamfeldt, L., Isbell, F., Lefcheck, J. S., Griffin, J. N., Hector, A., et al. (2014). Investigating the relationship between biodiversity and ecosystem multifunctionality: challenges and solutions. Methods Ecol. Evol. 5, 111–124. doi: 10.1111/2041-210X.12143
Cardinale, B. J., Duffy, J. E., Gonzalez, A., Hooper, D. U., Perrings, C., Venail, P., et al. (2012). Biodiversity loss and its impact on humanity. Nature 486, 59–67. doi: 10.1038/nature11148
Chen, H., Zhu, Q., Peng, C., Wu, N., Wang, Y., Fang, X., et al. (2013). The impacts of climate change and human activities on biogeochemical cycles on the Qinghai-Tibetan plateau. Glob. Chang. Biol. 19, 2940–2955. doi: 10.1111/gcb.12277
Cheng, L., Zhang, N. F., Yuan, M. T., Xiao, J., Qin, Y. J., Deng, Y., et al. (2017). Warming enhances old organic carbon decomposition through altering functional microbial communities. ISME J. 11, 1825–1835. doi: 10.1038/ismej.2017.48
DeAngelis, K. M., Pold, G., Topçuoğlu, B. D., van Diepen, L. T. A., Varney, R. M., Blanchard, J. L., et al. (2015). Long-term forest soil warming alters microbial communities in temperate forest soils. Front. Microbiol. 6:104. doi: 10.3389/fmicb.2015.00104
Delgado-Baquerizo, M., Maestre, F. T., Reich, P. B., Jeffries, T. C., Gaitan, J. J., Encinar, D., et al. (2016). Microbial diversity drives multifunctionality in terrestrial ecosystems. Nat. Commun. 7:10541. doi: 10.1038/ncomms10541
Deng, Y., Jiang, Y. H., Yang, Y. F., He, Z. L., Luo, F., and Zhou, J. Z. (2012). Molecular ecological network analyses. BMC Bioinform. 13:113. doi: 10.1186/1471-2105-13-113
Dove, N. C., Torn, M. S., Hart, S. C., and Taş, N. (2021). Metabolic capabilities mute positive response to direct and indirect impacts of warming throughout the soil profile. Nat. Commun. 12:2089. doi: 10.1038/s41467-021-22408-5
Falkowski, P. G., Fenchel, T., and Delong, E. F. (2008). The microbial engines that drive Earth’s biogeochemical cycles. Science 320, 1034–1039. doi: 10.1126/science.1153213
Feng, K., Zhang, Z. J., Cai, W. W., Liu, W. Z., Xu, M. Y., Yin, H. Q., et al. (2017). Biodiversity and species competition regulate the resilience of microbial biofilm community. Mol. Ecol. 26, 6170–6182. doi: 10.1111/mec.14356
Gross, N., Bagousse-Pinguet, Y. L., Liancourt, P., Berdugo, M., Gotelli, N. J., and Maestre, F. T. (2017). Functional trait diversity maximizes ecosystem multifunctionality. Nat. Ecol. Evol. 1:0132. doi: 10.1038/s41559-017-0132
Guimerà, R., and Nunes Amaral, L. A. (2005). Functional cartography of complex metabolic networks. Nature 433, 895–900. doi: 10.1038/nature03288
Guo, X., Feng, J. J., Shi, Z., Zhou, X. S., Yuan, M. T., Tao, X. Y., et al. (2018). Climate warming leads to divergent succession of grassland microbial communities. Nat. Clim. Chang. 8, 813–818. doi: 10.1038/s41558-018-0254-2
Heimann, M., and Reichstein, M. (2008). Terrestrial ecosystem carbon dynamics and climate feedbacks. Nature 451, 289–292. doi: 10.1038/nature06591
Hernandez, D. J., David, A. S., Menges, E. S., Searcy, C. A., and Afkhami, M. E. (2021). Environmental stress destabilizes microbial networks. ISME J. 15, 1722–1734. doi: 10.1038/s41396-020-00882-x
Huang, D. H., Zhou, L. L., Fan, H. M., Jia, Y. F., and Liu, M. B. (2021). Responses of aggregates and associated soil available phosphorus, and soil organic matter in different slope aspects, to seasonal freeze-thaw cycles in Northeast China. Geoderma 402:115184. doi: 10.1016/j.geoderma.2021.115184
Jiao, S., Peng, Z. H., Qi, J. J., Gao, J. M., and Wei, G. H. (2021). Linking bacterial-fungal relationships to microbial diversity and soil nutrient cycling. mSystems. 6, e01052–e01020. doi: 10.1128/mSystems.01052-20
Jiao, S., Xu, Y. Q., Zhang, J., Hao, X., and Lu, Y. H. (2019). Core microbiota in agricultural soils and their potential associations with nutrient cycling. mSystems 4, e00313–e00318. doi: 10.1128/mSystems.00313-18
Jiao, S., Yang, Y. F., Xu, Y. Q., Zhang, J., and Lu, Y. H. (2020). Balance between community assembly processes mediates species coexistence in agricultural soil microbiomes across eastern China. ISME J. 14, 202–216. doi: 10.1038/s41396-019-0522-9
Jing, X., Sanders, N. J., Shi, Y., Chu, H. Y., Classen, A. T., Zhao, K., et al. (2015). The links between ecosystem multifunctionality and above-and belowground biodiversity are mediated by climate. Nat. Commun. 6:8159. doi: 10.1038/ncomms9159
Klein, J. A., Harte, J., and Zhao, X. Q. (2004). Experimental warming causes large and rapid species loss, dampened by simulated grazing, on the Tibetan plateau. Ecol. Lett. 7, 1170–1179. doi: 10.1111/j.1461-0248.2004.00677.x
Landi, P., Minoarivelo, H. O., Brännström, Å., Hui, C., and Dieckmann, U. (2018). Complexity and stability of ecological networks: a review of the theory. Popul. Ecol. 60, 319–345. doi: 10.1007/s10144-018-0628-3
Le Bagousse-Pinguet, Y., Soliveres, S., Gross, N., Torices, R., Berdugo, M., and Maestre, F. T. (2019). Phylogenetic, functional, and taxonomic richness have both positive and negative effects on ecosystem multifunctionality. Proc. Natl. Acad. Sci. U. S. A. 116, 8419–8424. doi: 10.1073/pnas.1815727116
Li, Y. M., Lin, Q. Y., Wang, S. P., Li, X. Z., Liu, W., Luo, C. Y., et al. (2016). Soil bacterial community responses to warming and grazing in a Tibetan alpine meadow. FEMS Microbiol. Ecol. 92:fiv152. doi: 10.1093/femsec/fiv152
Li, Z. L., Tian, D. S., Wang, B. X., Wang, J. S., Wang, S., Chen, H. Y. H., et al. (2019). Microbes drive global soil nitrogen mineralization and availability. Glob. Chang. Biol. 25, 1078–1088. doi: 10.1111/gcb.14557
Liu, L. M., Yang, J., Yu, Z., and Wilkinson, D. M. (2015). The biogeography of abundant and rare bacterioplankton in the lakes and reservoirs of China. ISME J. 9, 2068–2077. doi: 10.1038/ismej.2015.29
Marion, G. M., Henry, G. H. R., Freckman, D. W., Johnstone, J., Jones, G., Jones, M. H., et al. (1997). Open-top designs for manipulating field temperature in high-latitude ecosystems. Glob. Chang. Biol. 3, 20–32. doi: 10.1111/j.1365-2486.1997.gcb136.x
Nielsen, U. N., Wall, D. H., and Six, J. (2015). Soil biodiversity and the environment. Annu. Rev. Env. Resour. 40, 63–90. doi: 10.1146/annurev-environ-102014-021257
Okuyama, T., and Holland, J. N. (2008). Network structural properties mediate the stability of mutualistic communities. Ecol. Lett. 11, 208–216. doi: 10.1111/j.1461-0248.2007.01137.x
Pasari, J. R., Levi, T., Zavaleta, E. S., and Tilman, D. (2013). Several scales of biodiversity affect ecosystem multifunctionality. Proc. Natl. Acad. Sci. U. S. A. 110, 10219–10222. doi: 10.1073/pnas.1220333110
Petchey, O. L., McPhearson, P. T., Casey, T. M., and Morin, P. J. (1999). Environmental warming alters food-web structure and ecosystem function. Nature 402, 69–72. doi: 10.1038/47023
Rui, J. P., Li, J. B., Wang, S. P., An, J. X., Liu, W. T., Lin, Q. Y., et al. (2015). Responses of bacterial communities to simulated climate changes in alpine meadow soil of the Qinghai-Tibet plateau. Appl. Environ. Microbiol. 81, 6070–6077. doi: 10.1128/AEM.00557-15
Tao, X. Y., Feng, J. J., Yang, Y. F., Wang, G. S., Tian, R. M., Fan, F. L., et al. (2020). Winter warming in Alaska accelerates lignin decomposition contributed by Proteobacteria. Microbiome 8:84. doi: 10.1186/s40168-020-00838-5
Vincent, W. F. (2010). Microbial ecosystem responses to rapid climate change in the Arctic. ISME J. 4, 1087–1090. doi: 10.1038/ismej.2010.108
Wagg, C., Bender, S. F., Widmer, F., and van der Heijden, M. G. A. (2014). Soil biodiversity and soil community composition determine ecosystem multifunctionality. Proc. Natl. Acad. Sci. U. S. A. 111, 5266–5270. doi: 10.1073/pnas.1320054111
Walkley, A., and Black, I. A. (1934). An examination of the Degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method. Soil Sci. 37, 29–38. doi: 10.1097/00010694-193401000-00003
Walther, G. R., Post, E., Convey, P., Menzel, A., Parmesan, C., Beebee, T. J. C., et al. (2002). Ecological responses to recent climate change. Nature 416, 389–395. doi: 10.1038/416389a
Wang, H., Liu, S. R., Schindlbacher, A., Wang, J. X., Yang, Y. J., Song, Z. C., et al. (2019). Experimental warming reduced topsoil carbon content and increased soil bacterial diversity in a subtropical planted forest. Soil Biol. Biochem. 133, 155–164. doi: 10.1016/j.soilbio.2019.03.004
Xue, K., Yuan, M. M., Shi, Z, J., Qin, Y., Deng, Y., Cheng, L., et al. (2016). Tundra soil carbon is vulnerable to rapid microbial decomposition under climate warming. Nat. Clim. Chang. 6, 595–600. doi: 10.1038/nclimate2940
Yuan, M. M., Guo, X., Wu, L. W., Zhang, Y., Xiao, N. J., Ning, D. L., et al. (2021). Climate warming enhances microbial network complexity and stability. Nat. Clim. Chang. 11, 343–348. doi: 10.1038/s41558-021-00989-9
Zheng, Q., Hu, Y. T., Zhang, S. S., Noll, L., Böckle, T., Dietrich, M., et al. (2019). Soil multifunctionality is affected by the soil environment and by microbial community composition and diversity. Soil Biol. Biochem. 136:107521. doi: 10.1016/j.soilbio.2019.107521
Zheng, H. F., Liu, Y., Chen, Y. M., Zhang, J., Li, H. J., Wang, L. F., et al. (2020). Short-term warming shifts microbial nutrient limitation without changing the bacterial community structure in an alpine timberline of the eastern Tibetan plateau. Geoderma 360:113985. doi: 10.1016/j.geoderma.2019.113985
Zhou, J. Z., Deng, Y., Luo, F., He, Z. L., Tu, Q. C., and Zhi, X. Y. (2010). Functional molecular ecological networks. MBio 1, e00169–e00110. doi: 10.1128/mBio.00169-10
Zhou, X., Fornara, D., Wasson, E. A., Wang, D. M., Ren, G. D., Christie, P., et al. (2015). Effects of 44 years of chronic nitrogen fertilization on the soil nitrifying community of permanent grassland. Soil Biol. Biochem. 91, 76–83. doi: 10.1016/j.soilbio.2015.08.031
Keywords: climate warming, profile, multi-nutrient cycling, keystones, bacteria
Citation: Zhou X, Chen X, Qi X, Zeng Y, Guo X, Zhuang G and Ma A (2023) Soil bacterial communities associated with multi-nutrient cycling under long-term warming in the alpine meadow. Front. Microbiol. 14:1136187. doi: 10.3389/fmicb.2023.1136187
Edited by:
Congrui Jin, Binghamton University, United StatesReviewed by:
Kashif Rahim, Cholistan University of Veterinary and Animal Sciences, PakistanSudeep Tiwari, Ben-Gurion University of the Negev, Israel
Copyright © 2023 Zhou, Chen, Qi, Zeng, Guo, Zhuang and Ma. 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: Anzhou Ma, ✉ azma@rcees.ac.cn