Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 22 September 2023
Sec. Terrestrial Microbiology
This article is part of the Research Topic Anthropogenic Effects on the Microbial Communities of Terrestrial Ecosystems View all 30 articles

Elevational distribution and seasonal dynamics of alpine soil prokaryotic communities

Junpeng Rui,Junpeng Rui1,2Yuwei Zhao,Yuwei Zhao1,2Nan CongNan Cong3Fuxin WangFuxin Wang1Chao LiChao Li1Xiang LiuXiang Liu1Jingjing HuJingjing Hu1Ning LingNing Ling2Xin Jing
Xin Jing4*
  • 1State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, College of Ecology, Lanzhou University, Lanzhou, China
  • 2Center for Grassland Microbiome, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, China
  • 3Lhasa Plateau Ecosystem Research Station, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China
  • 4State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, China

The alpine grassland ecosystem is a biodiversity hotspot of plants on the Qinghai-Tibetan Plateau, where rapid climate change is altering the patterns of plant biodiversity along elevational and seasonal gradients of environments. However, how belowground microbial biodiversity changes along elevational gradient during the growing season is not well understood yet. Here, we investigated the elevational distribution of soil prokaryotic communities by using 16S rRNA amplicon sequencing along an elevational gradient between 3,200 and 4,200 m, and a seasonal gradient between June and September in the Qinghai-Tibetan alpine grasslands. First, we found soil prokaryotic diversity and community composition significantly shifted along the elevational gradient, mainly driven by soil temperature and moisture. Species richness did not show consistent elevational trends, while those of evenness declined with elevation. Copiotrophs and symbiotic diazotrophs declined with elevation, while oligotrophs and AOB increased, affected by temperature. Anaerobic or facultatively anaerobic bacteria and AOA were hump-shaped, mainly influenced by moisture. Second, seasonal patterns of community composition were mainly driven by aboveground biomass, precipitation, and soil temperature. The seasonal dynamics of community composition indicated that soil prokaryotic community, particularly Actinobacteria, was sensitive to short-term climate change, such as the monthly precipitation variation. At last, dispersal limitation consistently dominated the assembly process of soil prokaryotic communities along both elevational and seasonal gradients, especially for those of rare species, while the deterministic process of abundant species was relatively higher at drier sites and in drier July. The balance between deterministic and stochastic processes in abundant subcommunities might be strongly influenced by water conditions (precipitation/moisture). Our findings suggest that both elevation and season can alter the patterns of soil prokaryotic biodiversity in alpine grassland ecosystem of Qinghai-Tibetan Plateau, which is a biodiversity hotspot and is experiencing rapid climate change. This work provides new insights into the response of soil prokaryotic communities to changes in elevation and season, and helps us understand the temporal and spatial variations in such climate change-sensitive regions.

Introduction

As hotspots of biodiversity and climate change, mountains create a broad range of gradients in both abiotic and biotic settings (Martin and Bellingham, 2016), and provide ecologists with a natural setting for testing their hypotheses (Lomolino, 2001). Spatial distributions of flora and animals along elevational gradients have been extensively studied (Rahbek, 1995). The most prevalent elevational pattern is a hump-shaped species-richness pattern for plants and animals, followed by a decreasing pattern of species richness (Rahbek, 2005). However, the elevational patterns for soil microorganisms are not well understood that may not be the same as those of animals and plants (Fierer et al., 2011). Generally, elevational patterns of soil microbes are diverse and inconclusive, including five types (type I-V): increase, decrease, hump-shaped, U-shaped, and no discernable tread (Hendershot et al., 2017; Looby and Martin, 2020).

Many environmental factors change along elevation since it is a complex and indirect gradient (Fierer et al., 2011). For example, temperature declines with increased elevation (Körner, 2007), coupled with decreased aboveground biomass (Hendershot et al., 2017), and while precipitation is often hump-shaped because of the cloud zone at mid-elevation (Rahbek, 1995). Consequently, various environmental settings, including temperature, precipitation, vegetation, and soil characteristics, may have an impact on the diversity and community composition of soil microbes (Sundqvist et al., 2013; Singh et al., 2014; Nottingham et al., 2018). However, key factors driving the elevational distribution are still under debate. The key factor may be different owing to different vegetation types or sampling scales on the same mountain. For instance, soil pH drives bacterial community distribution along six elevations between 530 m and 2,200 m a.s.l. on Changbai Mountain (Shen et al., 2013), while the primary determinants of bacterial elevational distribution are soil carbon and nitrogen concentrations at a rather small scale from 2000 m to 2,500 m a.s.l. on the same mountain (Shen et al., 2015).

Environmental factors, including temperature, precipitation, vegetation, and soil characteristics, may change with seasons, resulting in seasonal dynamics of the soil microbes (Lazzaro et al., 2015). For example, since different soil microbial taxa may differ in their response to changes in soil moisture and temperature, the seasonal dynamics of soil temperature and moisture strongly influence the abundance of soil bacteria in a desert grassland (Bell et al., 2009). A study of the alpine ecosystem demonstrates that Acidobacteria was more abundant in spring, Actinobacteria and Bacteroides in winter, while Verrucomicrobia and Betaproteobacteria in summer (Lipson and Schmidt, 2004). Alpine ecosystems go through a remarkable variety of seasonal variations (Ernakovich et al., 2014). However, few studies have focused on seasonal dynamics of soil microbial communities along elevational gradients. According to a study along the tree line, elevational effects on soil microbial diversity outweigh seasonal influences (Shen et al., 2021). At every elevation, the bacterial and fungal community structures exhibit an annual cycle in the catchment of the Tiefen Glacier (Lazzaro et al., 2015).

The community assembly mechanism is a critical issue in microbial ecology (Stegen et al., 2012). Stochastic processes include homogenous dispersal, dispersal limitation, and drift, while deterministic processes consist of homogenous selection and heterogeneous selection (Stegen et al., 2015; Shi et al., 2022). Some studies reported that deterministic processes play critical roles in bacterial community assembly along elevation gradients in some sites of the Qinghai-Tibetan Plateau, where the range of soil pH was broad (Li et al., 2018; Shen et al., 2019). However, a transect study across the Qinghai-Tibetan Plateau showed that stochastic processes, notably dispersal limitation and drift, played a primary role in bacterial community assembly in alpine meadow soils (Kang et al., 2022).

Alpine grasslands, the dominant vegetation on the Qinghai-Tibetan Plateau, are fragile and sensitive to climate change. Soil microorganisms are mainly concentrated in surface soils of this area, and play important roles in C and N cycles (Zhao and Zhou, 1999). The mechanism of soil microbial distribution along an elevation gradient in alpine grasslands keeps unclear. Precipitation is a primary driving force on the bacterial community between 4,400 and 5,200 m a.s.l. over the Qinghai-Tibetan alpine grasslands at Damxung station (Yuan et al., 2014), while temperature may be more significant at Haibei station, Qinghai (Rui et al., 2015). Besides, the seasonal dynamics of soil microbial community along elevation may occur since climate, soil properties, and microbial biomass are seasonal various (Chen et al., 2021).

In this study, we aim to elucidate the elevational distribution and seasonal dynamics of soil prokaryotic communities along an elevational gradient (3200–4,200 m) in the Qinghai-Tibetan alpine grassland during the growing season. The main goal of this study is to answer the following questions: (i) What is the critical factor driving the elevational distribution of soil prokaryotic community? (ii) Whether soil prokaryotic community significantly changes during the growing season? and (iii) Whether the ecological processes of prokaryotic community assembly change along the elevational or seasonal gradients?

Materials and methods

Study site description and experimental design

The study site is located at the Haibei National Alpine Grassland Ecosystem Research Station (101°12′E, 37°37’N, ~ 3,200 m above sea level) on the northeastern Qinghai-Tibetan Plateau in Qinghai, China. The annual mean air temperature is −1.1°C, and the annual mean precipitation is 485 mm (Ma et al., 2017). The soil is classified as Mat-Gryic Cambisol, and is dominated by perennial plants (Ma et al., 2017). The growing season lasts from May to September, when more than 80% of the annual rainfall falls (Yang et al., 2014).

Study sites were located at six elevations along an elevational gradient (3,200, 3,400, 3,600, 3,800, 4,000, and 4,200 m). Six 1 m × 1 m quadrats were selected as replicates at each elevation. The geographic coordinates and vegetation information for each elevation were reported previously (Rui et al., 2022). We collected soil samples at the beginning of June, July, August, and September in 2021, respectively. Five soil cores (5 cm in diameter) with 0–10 cm depth were taken randomly within each quadrat, mixed as one sample, and sieved with a 2 mm mesh. The aboveground plant biomass and richness were surveyed at the same time (Supplementary Figure S1). Soil temperature was monitored using EM-50 devices (Decagon devices, United States) at 5 cm, and the mean monthly temperature was calculated. Methods of measuring soil moisture, pH, total nitrogen (TN), total organic carbon (TOC), nitrate, nitrite, and ammonium contents were reported previously (Rui et al., 2022). Monthly precipitations at 3200 m between June and September in 2021 were 90, 64, 125, and 104 mm, respectively. The data of soil and plant properties are in the Supplementary materials.

DNA extraction, qPCR, and sequencing

PowerSoil DNA Isolation kit (MO BIO Laboratories, United States) was used to extract genomic DNA in soils following the manufacturer’s instructions (Rui et al., 2022). The universal primers 515F and 806R were used for the amplification of the 16S rRNA gene for qPCR and sequencing (Borton et al., 2017). The amplifications for sequencing were conducted following the “16S Illumina Amplicon Protocol” on the website of the Earth Microbiome Project.1 Real-time PCR was performed on the Applied Biosystems 7,300 Sequence Detection System (ABI, United States). The qPCR reaction system was composed of 10 μL of 2X Taq Plus Master Mix (Vazyme Biotech Co., Ltd., China), 0.8 μL of forward primer (5 μmol L−1), 0.8 μL of reverse primer, 1 μL of DNA, and 7.4 μL of ddH2O. The PCR program included an initial denaturation at 95°C for 5 min, followed by 35 cycles of 95°C for 30 s, 58°C for 30 s, and 72°C for 1 min. Sequencing libraries were generated with TruSeq® DNA Kit (Illumina, United States), and sequenced with the Reagent Kit v2 (2 × 250 bp) on the Illumina NovaSeq 6,000 platform (Novogene, Beijing, China) (Rui et al., 2022).

Sequencing data analysis

FLASh v1.2.11 was used to merge the paired-end reads (Magoč and Salzberg, 2011). Primers were removed from the merged reads using the Perl script trim_primer_in fq.pl (All Perl scripts mentioned in this article are available on GitHub at https://github.com/PeterRui/perl4amplicon). Low-quality reads were filtered using Trimomatic version 0.39 (Bolger et al., 2014). Chimeras were checked using Usearch 11.0.667.2 Amplicon sequence variants (ASVs) were obtained using the Unoise3 algorithm (Edgar, 2016a). The representative reads of ASVs were aligned with PyNAST, and the phylogenetic tree was constructed with Qiime 1.9.1 (Caporaso et al., 2010). The taxonomy of ASV-representative sequences was predicted with the RDP training set v16 database using the SINTAX algorithm (Edgar, 2016b). Moreover, chloroplasts and mitochondria from eukaryotic cells were checked using Greengenes v13.5 and SILVA v123 database using SINTAX. Sequences of each sample in the ASV-table were randomly normalized to the same sequence depth using the Perl script subsample_in_table.pl. After subsample, a total of 3,947,600 high-quality 16S rRNA sequences were obtained from 142 samples, affiliated to 41,750 ASVs. Alpha diversity indices, such as observed species, Chao1 richness, Shannon diversity, reciprocal Simpson diversity (1/D), and Pielou evenness index were calculated with the script alpha_diversity.pl. The Supplementary material contains information on using Perl scripts.

Statistical analyses

We performed statistical analyses based on the ASV-table and/or environmental factors using R (Version 4.1.0). Details of the statistical analysis can be found in a previous publication (Rui et al., 2022). For example, the function varclus in package Hmisc and vif.cca in package vegan were used to estimate collinearities among environmental factors. Using the function HSD.test in package agricolae, ANOVA with Tukey’s test was calculated to check the differences of prokaryotic groups or environmental factors among elevations/months. Bray-Curtis distance-based PCoA was conducted using cmdscale in vegan. Environmental factors were fitted onto the PCoA plots using envfit in vegan. Adonis (i.e., PerMANOVA) was estimated using adonis to evaluate community differences between sample groups. To calculate the multivariate homogeneity of community dispersions, PERMDISP was estimated using betadisper. Partial Mantel test and multiple regression on distance matrices (MRM) were conducted with the package ecodist to evaluate the effects of environmental factors on the prokaryotic community. The package rdacca.hp was used to estimate the independent contribution of each factor to prokaryotic community variation (Lai et al., 2022).

Community assembly processes of subcommunities

Abundant and rare subcommunities may have different ecological roles and community assembly mechanisms (Logares et al., 2014; Jiao et al., 2017). The definition of abundant, intermediate, and rare species were modified based on Jiao’s method (Jiao and Lu, 2020). Specifically, across all samples, ASVs with average relative abundances <0.001% and > 0.01% were defined as rare and abundant species, respectively. Those between 0.001% and 0.01% were defined as intermediate species. The community assembly processes of abundant, intermediate, and rare species were evaluated with iCAMP, using the confidence index for the null model significance test. ASVs were divided into different “bins” based on the phylogenetic tree, and their ecological processes were identified using null model analysis (Ning et al., 2020).

Results

Elevational shift of prokaryotic diversity and community composition

Alpha diversity indices changed significantly along the elevational gradient (two-way ANOVA p < 0.001). In general, elevational changes of evenness were more obvious than that of species richness (Observed species and Chao1 index). Elevational patterns of species richness were almost no trend in most cases (Figures 1A,B), while those indices including Shannon diversity, reciprocal Simpson diversity, and Pielou evenness declined with elevation, although not always monotonic (Figures 1CE). Soil temperature was the most important factor in elevational shift, positively correlated to alpha diversity, especially evenness (p < 0.01).

FIGURE 1
www.frontiersin.org

Figure 1. Elevational changes of prokaryotic alpha diversity indices during the growing season. Error bars represent standard deviations (n = 6).

The community structures at low elevations (3200–3,400 m) were substantially different from those at high elevations (3600–4,200 m), which could be distinguished along axis 1 of PCoA (Figure 2A). In the same month, most elevations we sampled had substantially distinct prokaryotic community structures (Supplementary Figure S2). Soil temperature was the critical factor driving the prokaryotic elevational distribution in each month (Tables 13). Temperature was strongly correlated to prokaryotic community structure, best fitted in PCoA patterns in each month (R2 ≥ 0.75, Supplementary Figure S3). Moisture was another important factor to elevational distribution, particularly in July (Tables 13). The effects of temperature and moisture decreased in August, in which the community dispersions were higher (Supplementary Figure S4). In addition, soil pH was important to elevational distribution in June, while plant richness in July (Tables 1, 2). It seems that aboveground biomass was less important because of its collinearity with temperature. The impact of precipitation along the elevation gradient was unclear due to lack of data.

FIGURE 2
www.frontiersin.org

Figure 2. PCoA plots of prokaryotic community structure based on Bray-Curtis distance. Plots include all samples with (A) axis 1 vs. axis 2, and (B) with axis 1 vs. axis 3; samples at 3200 m (C) with axis 1 vs. 2, and (D) with axis 2 vs. 3, respectively. Important environmental factors were fitted onto the patterns. Values on axes indicated the percentages of total variation explained by each axis.

TABLE 1
www.frontiersin.org

Table 1. Spearman’s correlation coefficients of prokaryotic community structure with environmental factors based on partial Mantel test.

TABLE 2
www.frontiersin.org

Table 2. MRM regression coefficients of prokaryotic community structure with environmental factors.

TABLE 3
www.frontiersin.org

Table 3. The independent contribution of each factor to prokaryotic community variation based on hierarchical partitioning in canonical analysis.

Across all samples, 58 phyla and 943 genera were detected. Proteobacteria, Actinobacteria, and Acidobacteria were dominant phyla, accounting for >60% of relative abundance in most samples (Figure 3). Besides, Verrucomicrobia, Bacteroidetes, Gemmatimonadetes, Planctomycetes, Nitrospirae, Chloroflexi, Firmicutes, and Thaumarchaeota were also important phyla with high relative abundance.

FIGURE 3
www.frontiersin.org

Figure 3. Elevational shifts of important phyla (A–K) and proteobacterial classes (L–O) in each month. Error bars represent standard deviations (n = 6).

Elevational patterns of community composition were diverse at the genus/phylum level, including increase (type I, genus level 2.3%/phylum level 4.9%), decrease (type II, 7.8%/3.3%), hump-shaped (type III, 3.2%/13.1%), U-shaped (type IV, 1.4%/1.6%) and no trend (type V, 85%/77%) (Supplementary Figures S5–S9). For example, Verrucomicrobia and Betaproteobacteria (mainly Nitrosomonadaceae) increased with elevation, while Alphaproteobacteria decreased (Figure 3). Patterns of Nitrospirae, Thaumarchaeota, and Deltaproteobacteria were hump-shaped, while Gemmatimonadetes was opposite. Besides, elevational patterns of many phyla were no trend, such as Acidobacteria, Actinobacteria, Firmicutes, and Chloroflexi. Besides, many rare taxa did not show any trends because of low abundance and low frequency.

At the genus level, type I genera were mainly composed of Nitrospira, Rhodoplanes, Rhodomicrobium, Chthoniobacter, Bauldia, and Terrimicrobium (Supplementary Figure S5). Several actinobacterial genera were type II, such as Oryzihumus, Blastococcus, Mycobacterium, and Arthrobacter (Supplementary Figure S6). Moreover, some rhizobia also belonged to type II, such as Mesorhizobium, Rhizobium, and Devosia. Type III genera were hump-shaped, such as Nitrososphaera and Candidatus Solibacter (Supplementary Figure S7). Type IV genera were lower at middle elevations, such as acidobacterial genera Terrimonas and Aridibacter (Supplementary Figure S8). Type I and II genera were significantly correlated to soil temperature, aboveground biomass, and altitude (Supplementary Table S1). Most type III and IV genera were significantly correlated to soil moisture.

Seasonal dynamics of prokaryotic diversity and community composition

Species richness exhibited more pronounced seasonal changes than evenness (Supplementary Figure S10). Both richness and evenness were lowest in September. Richness indices declined from July to September (Supplementary Figures S10A,B). However, richness indices in June were not significantly different from those in July. Precipitation was important to seasonal dynamics of alpha diversity at 3200 m, negatively correlated to these indices, especially to richness (p < 0.01, Supplementary Table S2).

The community structures in June were quite different from those in other months, which was obvious along axis 2 of PCoA (Figure 2A). Moreover, samples in July were separated from those in August and September by axis 3 (Figure 2B). Community dissimilarities between sample groups (i.e., between elevations and months) were higher than those within each group (Supplementary Figure S2). Therefore, prokaryotic community succession during the growing season could be divided into 3 stages, that is early stage (June), middle stage (July), and late stage (August and September). Community structures (i.e., beta-diversity) were significantly different among months, even between the same month in different years (Figure 2D). According to Adonis’ analysis, altitude, sampling season and their interaction could explain 19%, 22%, and 12% of community variations in all samples, respectively, while 22%, 14%, and 12% if samples in June were excluded, indicating that the difference between early stage and other stages contributed more to community variation.

Seasonal dynamics of the prokaryotic community were mainly affected by aboveground biomass, precipitation, and temperature (Tables 1, 2, Figures 2C,D). Aboveground biomass was more important to seasonal dynamics at low elevations, since it increased from June to August (Supplementary Figure S1E), while temperature played more important roles at higher elevations (Tables 13). However, the seasonal changes of aboveground biomass were less significant at high elevations (Supplementary Figure S1E), resulting in less importance of aboveground biomass. The community dispersions were high at 3800 and 4,200 m (Supplementary Figure S4), weakening the impacts of temperature. Soil nitrate was significantly correlated to prokaryotic community structures at 3800 m (Tables 1, 2), indicating the importance of nitrogen cycling-related microorganisms there. Results at 3200 m revealed that precipitation was more important than soil moisture to seasonal dynamics (Tables 13; Figures 2C,D).

Seasonal patterns of community composition were also diverse during the growing season in 2021, including increase (type 1, 44%/36%), decrease (type 2, 9.1%/28%), hump-shaped (type 3, 8.6%/6.9%), V-shaped (type 4, 1.3%/3.4%), and no trend (type 5, 37%/26%) at the genus/phylum level. For example, Thaumarchaeota and Firmicutes increased and were most abundant at the late stage (i.e., August and September, type 1, Supplementary Figure S11), while Acidobacteria, Nitrospirae, and Deltaproteobacteria (mainly Syntrophobacteraceae) were most abundant in June and decreased during the growing season (type 2). Actinobacteria were most abundant in July (type 3), while Proteobacteria (such as beta and gamma classes), Gemmatimonadetes, and Verrucomicrobia indicated lowest in July (type 4). Those phyla with no trend were mainly rare taxa.

Community assembly processes of abundant, intermediate, and rare species

Across all samples, abundant, intermediate, and rare species consisted of 116, 10,568, and 31,065 ASVs, respectively, and their total relative abundances were 24.8%, 65.6%, and 9.6%, respectively (Figure 4). The PCoA plots of abundant and intermediate species were similar to those of all samples, that is, the subcommunity structure in June was quite different from other months (Supplementary Figures S12A,B). Key factors affecting these subcommunities were also similar to that of all samples, especially the intermediate species (Tables 13, Supplementary Table S3). However, the subcommunity structure of rare species in September was quite different from other months (Supplementary Figures S12C,D). The key factors influencing rare subcommunities were nitrate, C: N ratio, plant richness, and aboveground biomass (Supplementary Table S3).

FIGURE 4
www.frontiersin.org

Figure 4. Relative abundances of phyla for abundant, intermediate and rare subcommunities (A) along the elevational gradient and (B) during the growing season. Error bars represent standard deviations (n = 6).

The primary assembly mechanism, particularly in rare subcommunities, was dispersal limiting (Figure 5). The proportions of contribution from different assembling processes were relatively steady in rare and intermediate subcommunities. The relative importance of homogeneous selection was higher in intermediate subcommunities than that in rare subcommunities. Besides, ecological drift was almost only observed in intermediate subcommunities.

FIGURE 5
www.frontiersin.org

Figure 5. Relative importance of different ecological processes for (A) abundant, (B) intermediate, and (C) rare subcommunities along the elevational gradient during the growing season.

In abundant subcommunities, the relative importance of different processes shifted variously, especially homogeneous selection and heterogeneous selection. The elevation pattern of the deterministic process was opposite to that of moisture, and this process was higher in drier July. Thus, water conditions (precipitation/moisture) might be the key factor affecting the balance between deterministic and stochastic processes in abundant subcommunities. Generally, dispersal limitation was lower at low elevations but higher at middle elevations, and the lowest in July (Figure 5A). Its relative importance was contributed by almost all abundant species (Supplementary Table S4). The homogeneous selection was higher in drier July, and it was also higher at 3200 and 4,200 m, indicating drought might increase its relative importance. Its relative importance was mainly contributed by Betaproteobacteria in June, Actinobacteria (“bin2”) in July, while Acidobacteria (classes Vicinamibacteria and Holophagae) and Verrucomicorbia in September (Supplementary Table S4). The heterogeneous selection was higher at low elevations, while it was lower at middle elevations, especially in July. Its relative importance was mainly contributed by Acidobacteria (classes Blastocatellia and Gp16) in August and September, while by Actinobacteria (“bin2”) in June (Supplementary Table S4). Interestingly, the roles of “bin2” changed, and its relative abundance increased dramatically from wetter June to drier July (Supplementary Figure S13). Therefore, community assemblies of different taxa might be different even though in the same place.

Discussion

The main finding of this study was that soil prokaryotic diversity and community composition significantly shifted along the elevational gradient, mainly driven by soil temperature, followed by soil moisture. Seasonal patterns of community composition were also diverse and driven mainly by aboveground biomass, precipitation, and temperature. Water conditions might be the key drivers for the balance between deterministic and stochastic processes in abundant subcommunities along the elevational and seasonal gradients.

Elevational shift of prokaryotic diversity and community composition

Our results clearly showed that temperature was the critical factor driving prokaryotic elevational distribution. Temperature is one of the most important factors affecting soil microbial community. First, high temperatures may directly alter the microbial community by increasing the metabolism and growth rates of microorganisms (Brown et al., 2004). Second, a higher temperature may increase aboveground biomass and alter plant community composition, then indirectly affecting the soil microbial community (Wang et al., 2009). Third, a higher temperature may increase the degradation rate of organic matter such as litter, and then increase the substrates for microbes (Zhou et al., 2016). In our study, aboveground biomass and richness were higher at low elevations with higher temperatures (Supplementary Figure S1). Our previous study showed that temperature indirectly affected prokaryotic communities through vegetation (Rui et al., 2015). Some important taxa were affected by temperature. For instance, Verrucomicrobia and Betaproteobacteria increased with elevation, while Alphaproteobacteria decreased (Figure 3). Verrucomicrobia is a common oligotroph in soils, adaptable to the high elevation where available nutrients are limited (Lipson and Schmidt, 2004; Yuan et al., 2021). Nitrosomonadaceae, the ammonia-oxidizing bacterium (AOB) (Garrity et al., 2015), was the most abundant family in Betaproteobacteria in this study, which increased with elevation. Many members of Alphaproteobacteria are copiotrophs, resulting in the decrease pattern of this class, in line with our previous study (Rui et al., 2015). Moreover, alphaproteobacterial genera Mesorhizobium, Rhizobium, and Devosia declined with elevation, which are symbiotic rhizobia associated with legumes (Franche et al., 2009). Our previous study showed that legume biomass is higher at low elevation, resulting in higher symbiotic nitrogen fixation and a higher abundance of symbiotic rhizobia (Rui et al., 2022).

Soil moisture was another important factor driving prokaryotic elevation distribution, although its effect was smaller than temperature (Tables 13). The physiological state of microorganisms could be directly influenced by soil moisture (Prado and Airoldi, 1999), and indirectly affect them by altering vegetation and soil properties (McHugh et al., 2014; Chen et al., 2015). Excessive soil moisture may lead to an anaerobic environment unfavorable to aerobic microbes (Unger et al., 2009). Due to the influence of topography on rainfall in this study, the soil moisture got maximum at 3600 m during the growth season at the study location (Xu et al., 2010). Many soil properties were significantly correlated to moisture, such as TN, TOC, nitrite contents, and pH (p < 0.001). Our previous study showed that anaerobic free-living diazotrophs are most abundant at 3600 m, indicating the anaerobic microenvironment there (Rui et al., 2022). In this study, some type III taxa (hump-shaped) were anaerobic or facultatively anaerobic, such as Defluviicoccus, Hyphomicrobium, and Phaselicystis (Supplementary Figure S7). Nitrososphaera, known as common ammonia-oxidizing archaea (AOA) in soils (Tourna et al., 2011), was higher at mid-elevation. The positive correlation between Nitrososphaera and soil moisture/precipitation is identical to other studies (Gu et al., 2017; Yang et al., 2018). In contrast, type IV phylum Gemmatimonadetes was U-shaped. Gemmatimonadetes prefers low moisture (DeBruyn et al., 2011), and is abundant in dry soils (Singh et al., 2014). Moisture also indirectly affected microbial taxa through soil pH. For example, Candidatus Solibacter, a member of Acidobacteria subdivision 3, was hump-shaped and displayed highest at 3600 m because low pH was preferred (Halamka et al., 2022). The U-shaped Aridibacter and Blastocatella are members of Acidobacteria subdivision 4, which prefers high pH, common in arid soils and crusts (Foesel et al., 2013).

Soil pH is often recognized as the critical driver of microbial elevational distribution in the literature (Bryant et al., 2008; Looby and Martin, 2020). However, it is not the case in this work. The importance of soil pH depends on the range of soil pH. A very narrow range of pH perhaps lead to less importance of soil pH (Singh et al., 2014). Interestingly, substantial changes in soil pH are often observed in elevational gradients where soil pH is the key factor driving bacterial elevational distribution (Shen et al., 2013; Xu et al., 2014; Li et al., 2018). In our study, the range of soil pH was narrow (Supplementary Figure S1), resulting in less importance of soil pH.

In summary, soil temperature and moisture drove the elevational distribution of soil prokaryotic community. Copiotrophs and symbiotic diazotrophs declined with elevation, while oligotrophs and AOB increased, affected by temperature. Anaerobic or facultatively anaerobic bacteria and AOA were hump-shaped, mainly influenced by moisture.

Seasonal dynamics of prokaryotic diversity and community composition

In this study, we found obvious seasonal dynamics of soil prokaryotic community at each elevation, which were mainly affected by aboveground biomass, precipitation, and soil temperature. Plants primarily influenced the soil microbial community by adding litter and root exudates (Knelman et al., 2012). Root exudates may change during the plant growth cycle, and alter soil microbial community structures, especially in rhizospheres (Dunfield and Germida, 2003; Li et al., 2014). At the early stage of the growing season, aboveground biomass was much lower than in the mid and late stages, indicating available carbon input from plants might be also lower. In this study, relative abundances of Acidobacteria, Nitrospirae, and Deltaproteobacteria were relatively higher in June (Supplementary Figure S11). Acidobacteria is a keystone oligotroph in soils, abundant in bulk soils with low resource availability (Fierer et al., 2007). It plays an important role in the degradation of soil organic matter (Banerjee et al., 2016). Nitrospira, the dominant genus of Nitrospirae in our study, is able to oxidize ammonia to nitrite (Daims and Wagner, 2018), perhaps related to high ammonium and low nitrite contents in June (Supplementary Figure S1). Deltaproteobacteria were dominated by Syntrophobacteraceae, related to anaerobic degradation of organic substrates and sulfate reduction (Kuever et al., 2015). Interestingly, relative abundances of rhizobia, such as Bradyrhizobium, Mesorhizobium, and Rhizobium, were the lowest in June, owing to the low biomass of legumes. Thus, symbiotic nitrogen fixation might be less important in June, and the high ammonium content could be from other pathways, such as free-living nitrogen fixation (Rui et al., 2022). It is necessary to investigate the seasonal dynamics of the diazotrophic community and measure the nitrogen fixation rate here in the future.

Actinobacteria plays important roles in carbon cycling, especially in the solubilization of fungal and plant cell walls (Lacombe-Harvey et al., 2018). The relative abundance of Actinobacteria was the lowest in June due to the lowest aboveground biomass. When aboveground biomass increased with lower precipitation in July, Actinobacteria increased dramatically and reached the highest level (Supplementary Figure S11). The relative abundance of Actinobacteria was negatively correlated to precipitation in soils during July and September (p < 0.001), indicating precipitation became the main factor for Actinobacteria, in line with another study (Yao et al., 2017). Seasonal dynamics of Actinobacteria in grassland soils were also observed in another study, caused by seasonal variation of rainfall (Jenkins et al., 2009). In our study, several actinobacterial genera such as Arthrobacter were positively correlated to aboveground biomass, which prefer soils receiving organic inputs (Jenkins et al., 2009). Compared to other phyla, the relative abundance of Actinobacteria changed dramatically during the growing season. Thus, Actinobacteria was the keystone indicator of seasonal dynamics of the soil prokaryotic community in the alpine meadow.

Our study showed that several important phyla changed almost linearly during the growing season. Community structures were more similar between the same months in different years than those between different stages (Figure 2). The seasonal dynamics indicated that the soil prokaryotic community was sensitive to short-term climate change, such as monthly variation of precipitation. Thus, seasonal variation should be paid attention to the time of soil sampling for bacterial community sequencing.

Community assembly processes of abundant, intermediate, and rare species

We found that dispersal limitation was the main assembly process for the prokaryotic community in alpine grassland soils, in accordance with the transect study across the Qinghai-Tibetan Plateau (Kang et al., 2022). Two reasons could explain it. First, soil microbes produce dormancy to deal with stressors such as low temperature, and microbial dormancy serves as a seed bank that may weaken the selection effects of community assembly, but increase the relative importance of stochastic processes (Nemergut et al., 2013). Second, soil freezing is a physical barrier that can limit the movement of microbes (Steven et al., 2008), increasing dispersal limitation (Kang et al., 2022). Permafrost and seasonally frozen zones cover 96% of the areas on the Qinghai-Tibetan Plateau (Zou et al., 2017). Therefore, it is not surprising that dispersal limitation is the primary process for the prokaryotic community in alpine grassland soils (Kang et al., 2022).

Different ecological processes were dominant in different subcommunities. First, the relative importance of dispersal limitation was the highest in rare subcommunities. We suspect that the reason is that species with higher abundance may have a greater probability of dispersal than rare species (Nemergut et al., 2011; Zhou and Ning, 2017). In this study, abundant species were found in most samples, while rare species were in only 10 samples on average. Second, the relative importance of drift was the highest in intermediate subcommunities. Drift refers to random changes due to birth, death, and reproduction (Zhou and Ning, 2017). A large fraction of a rare biosphere consists of nongrowing cells, such as dormancy cells and spores (Pedrós-Alió, 2012). Therefore, drift was seldom observed in rare subcommunities in our study. Third, heterogeneous selection was more important to abundant subcommunities, while the homogeneous selection was more important to intermediate subcommunities. Heterogeneous selection (i.e., variable selection) refers to the heterogeneous changes of communities in selective pressure from environmental heterogeneous, while homogeneous selection refers to the homogeneous changes in a consistent selective pressure if environmental conditions are homogeneous (Dini-Andreote et al., 2015; Stegen et al., 2015). The environmental conditions were relatively homogeneous at the same elevation and at the same time, causing strong homogeneous selection. Abundant species usually grow fast and contribute more to nutrient cycling and carbon flow (Pedrós-Alió, 2012). Thus, the response of abundant species to selective pressure might be more active than immediate and rare species.

A transect research reported that precipitation can balance deterministic and stochastic processes in grassland soils (Yang et al., 2022). We found that water conditions might be the key factor affecting the balance between deterministic and stochastic processes in abundant subcommunities. Actinobacteria played an important role here (Supplementary Table S4). From wetter June to drier July, the relative abundance of Actinobacteria increased dramatically. As a result, communities at the same elevation became similar (Supplementary Figure S4), in line with the definition of homogeneous selection. Therefore, water conditions might be the critical factor influencing community assembly.

Conclusion

In conclusion, soil prokaryotic diversity and community composition significantly shifted along the elevational gradient, mainly driven by soil temperature, followed by soil moisture. Elevational patterns of community composition were diverse, including increase, decrease, hump-shaped, U-shaped, and no trend. Moreover, seasonal patterns of prokaryotic community composition were also diverse and driven mainly by aboveground biomass, precipitation, and temperature. Actinobacteria was the keystone indicator of seasonal dynamics, more abundant in drier July. Dispersal limitation consistently dominated the assembly process of soil prokaryotic communities along both elevational and seasonal gradients. The deterministic process of abundant species was relatively higher at drier sites and in drier July. Our findings suggest that rapid changes in growing season climate can alter the patterns of soil prokaryotic biodiversity along both elevational and seasonal environmental gradients, especially carbon/nitrogen cycling microorganisms. This work highlights that temperature, water conditions, and plant productivity are critical factors driving soil prokaryotic community along elevations during the growing season, and provides new insights into the response of soil prokaryotic community to changes in elevation and season, and helps us understand the temporal and spatial fluctuations of the climate-sensitive Qinghai-Tibetan Plateau.

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 at: https://ngdc.cncb.ac.cn/gsa/browse/CRA009192, CRA009192.

Author contributions

JR: Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Software, Validation, Visualization, Writing – original draft, review & editing. YZ: Investigation, Resources, Writing – original draft. NC: Methodology, Writing – review & editing. FW: Investigation, Writing – original draft. CL: Investigation, Writing – original draft. XL: Resources, Writing – review & editing. JH: Investigation, Writing – original draft. NL: Funding acquisition, Writing – review & editing. XJ: Project administration, Supervision, Writing – original draft, review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by grants from the Natural Science Foundation for Distinguished Young Scholars of Gansu Province (grant number 20JR5RA308), the National Natural Science Foundation of China (grant numbers 41771293 and 42277100), and the Fundamental Research Funds for the Central Universities (grant numbers lzujbky-2019-ct01, lzujbky-2022-ct04 and lzujbky-2023-ct04).

Acknowledgments

We thank Juan Du (Lanzhou University) for helping with soil property measurement. We also thank Professor Jin-Sheng He and the Haibei National Alpine Grassland Ecosystem Research Station for helping with environmental data collection.

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.1280011/full#supplementary-material

Footnotes

References

Banerjee, S., Kirkby, C. A., Schmutter, D., Bissett, A., Kirkegaard, J. A., and Richardson, A. E. (2016). Network analysis reveals functional redundancy and keystone taxa amongst bacterial and fungal communities during organic matter decomposition in an arable soil. Soil Biol. Biochem. 97, 188–198. doi: 10.1016/j.soilbio.2016.03.017

CrossRef Full Text | Google Scholar

Bell, C. W., Acosta-Martinez, V., McIntyre, N. E., Cox, S., Tissue, D. T., and Zak, J. C. (2009). Linking microbial community structure and function to seasonal differences in soil moisture and temperature in a Chihuahuan desert grassland. Microb. Ecol. 58, 827–842. doi: 10.1007/s00248-009-9529-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Borton, M. A., Sabag-Daigle, A., Wu, J., Solden, L. M., O’Banion, B. S., Daly, R. A., et al. (2017). Chemical and pathogen-induced inflammation disrupt the murine intestinal microbiome. Microbiome 5:47. doi: 10.1186/s40168-017-0264-8

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Bryant, J. A., Lamanna, C., Morlon, H., Kerkhoff, A. J., Enquist, B. J., and Green, J. L. (2008). Microbes on mountainsides: contrasting elevational patterns of bacterial and plant diversity. Proc. Natl. Acad. Sci. USA 105, 11505–11511. doi: 10.1073/pnas.0801920105

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, D., Mi, J., Chu, P., Cheng, J., Zhang, L., Pan, Q., et al. (2015). Patterns and drivers of soil microbial communities along a precipitation gradient on the Mongolian plateau. Landsc. Ecol. 30, 1669–1682. doi: 10.1007/s10980-014-9996-z

CrossRef Full Text | Google Scholar

Chen, Y., Han, M., Yuan, X., Cao, G., and Zhu, B. (2021). Seasonal changes in soil properties, microbial biomass and enzyme activities across the soil profile in two alpine ecosystems. Soil Ecol. Lett. 3, 383–394. doi: 10.1007/s42832-021-0101-7

CrossRef Full Text | Google Scholar

Daims, H., and Wagner, M. (2018). Nitrospira. Trends Microbiol. 26, 462–463. doi: 10.1016/j.tim.2018.02.001

CrossRef Full Text | Google Scholar

DeBruyn, J. M., Nixon, L. T., Fawaz, M. N., Johnson, A. M., and Radosevich, M. (2011). Global biogeography and quantitative seasonal dynamics of Gemmatimonadetes in soil. Appl. Environ. Microbiol. 77, 6295–6300. doi: 10.1128/AEM.05005-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Dini-Andreote, F., Stegen, J. C., Van Elsas, J. D., and Salles, J. F. (2015). Disentangling mechanisms that mediate the balance between stochastic and deterministic processes in microbial succession. Proc. Natl. Acad. Sci. U. S. A. 112, E1326–E1332. doi: 10.1073/pnas.1414261112

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunfield, K. E., and Germida, J. J. (2003). Seasonal changes in the rhizosphere microbial communities associated with field-grown genetically modified canola (Brassica napus). Appl. Environ. Microbiol. 69, 7310–7318. doi: 10.1128/AEM.69.12.7310-7318.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2016a). UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. bio Rxiv :081257. doi: 10.1101/081257

CrossRef Full Text | Google Scholar

Edgar, R. C. (2016b). SINTAX: a simple non-Bayesian taxonomy classifier for 16S and ITS sequences. bio Rxiv :074161. doi: 10.1101/074161

CrossRef Full Text | Google Scholar

Ernakovich, J. G., Hopping, K. A., Berdanier, A. B., Simpson, R. T., Kachergis, E. J., Steltzer, H., et al. (2014). Predicted responses of arctic and alpine ecosystems to altered seasonality under climate change. Glob. Chang. Biol. 20, 3256–3269. doi: 10.1111/gcb.12568

PubMed Abstract | CrossRef Full Text | Google Scholar

Fierer, N., Bradford, M. A., and Jackson, R. B. (2007). Toward an ecological classification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839

CrossRef Full Text | Google Scholar

Fierer, N., McCain, C. M., Meir, P., Zimmermann, M., Rapp, J. M., Silman, M. R., et al. (2011). Microbes do not follow the elevational diversity patterns of plants and animals. Ecology 92, 797–804. doi: 10.1890/10-1170.1

CrossRef Full Text | Google Scholar

Foesel, B. U., Rohde, M., and Overmann, J. (2013). Blastocatella fastidiosa gen. Nov., sp. nov., isolated from semiarid savanna soil–the first described species of Acidobacteria subdivision 4. Syst. Appl. Microbiol. 36, 82–89. doi: 10.1016/j.syapm.2012.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Franche, C., Lindström, K., and Elmerich, C. (2009). Nitrogen-fixing bacteria associated with leguminous and non-leguminous plants. Plant Soil 321, 35–59. doi: 10.1007/s11104-008-9833-8

CrossRef Full Text | Google Scholar

Garrity, G. M., Bell, J. A., and Lilburn, T. (2015). “Nitrosomonadaceae fam. nov” in Bergey's Manual of Systematics of Archaea and Bacteria, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer), 901–908.

Google Scholar

Gu, Y., Wang, Y., Lu, S., Xiang, Q., Yu, X., Zhao, K., et al. (2017). Long-term fertilization structures bacterial and archaeal communities along soil depth gradient in a paddy soil. Front. Microbiol. 8:1516. doi: 10.3389/fmicb.2017.01516

PubMed Abstract | CrossRef Full Text | Google Scholar

Halamka, T. A., Raberg, J. H., McFarlin, J. M., Younkin, A. D., Mulligan, C., Liu, X. L., et al. (2022). Production of diverse br GDGTs by Acidobacterium Solibacter usitatus in response to temperature, pH, and O2 provides a culturing perspective on br GDGT proxies and biosynthesis. Geobiology 21, 102–118. doi: 10.1111/gbi.12525

PubMed Abstract | CrossRef Full Text | Google Scholar

Hendershot, J. N., Read, Q. D., Henning, J. A., Sanders, N. J., and Classen, A. T. (2017). Consistently inconsistent drivers of microbial diversity and abundance at macroecological scales. Ecology 98, 1757–1763. doi: 10.1002/ecy.1829

PubMed Abstract | CrossRef Full Text | Google Scholar

Jenkins, S. N., Waite, I. S., Blackburn, A., Husband, R., Rushton, S. P., Manning, D. C., et al. (2009). Actinobacterial community dynamics in long term managed grasslands. Antonie Van Leeuwenhoek 95, 319–334. doi: 10.1007/s10482-009-9317-8

CrossRef Full Text | Google Scholar

Jiao, S., Chen, W., and Wei, G. (2017). Biogeography and ecological diversity patterns of rare and abundant bacteria in oil-contaminated soils. Mol. Ecol. 26, 5305–5317. doi: 10.1111/mec.14218

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiao, S., and Lu, Y. (2020). Abundant fungi adapt to broader environmental gradients than rare fungi in agricultural fields. Glob. Chang. Biol. 26, 4506–4520. doi: 10.1111/gcb.15130

CrossRef Full Text | Google Scholar

Kang, L., Chen, L., Zhang, D., Peng, Y., Song, Y., Kou, D., et al. (2022). Stochastic processes regulate belowground community assembly in alpine grasslands on the Tibetan plateau. Environ. Microbiol. 24, 179–194. doi: 10.1111/1462-2920.15827

PubMed Abstract | CrossRef Full Text | Google Scholar

Knelman, J. E., Legg, T. M., O’Neill, S. P., Washenberger, C. L., González, A., Cleveland, C. C., et al. (2012). Bacterial community structure and function change in association with colonizer plants during early primary succession in a glacier forefield. Soil Biol. Biochem. 46, 172–180. doi: 10.1016/j.soilbio.2011.12.001

CrossRef Full Text | Google Scholar

Körner, C. (2007). The use of ‘altitude’ in ecological research. Trends Ecol. Evol. 22, 569–574. doi: 10.1016/j.tree.2007.09.006

CrossRef Full Text | Google Scholar

Kuever, J., Rainey, F. A., and Widdel, F. (2015). “Syntrophobacteraceae fam. Nov” in Bergey's Manual of Systematics of Archaea and Bacteria, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer), 289–299.

Google Scholar

Lacombe-Harvey, M.-È., Brzezinski, R., and Beaulieu, C. (2018). Chitinolytic functions in actinobacteria: ecology, enzymes, and evolution. Appl. Microbiol. Biotechnol. 102, 7219–7230. doi: 10.1007/s00253-018-9149-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Lai, J., Zou, Y., Zhang, J., and Peres-Neto, P. R. (2022). Generalizing hierarchical and variation partitioning in multiple regression and canonical analyses using the rdacca.hp R package. Methods Ecol. Evol. 13, 782–788. doi: 10.1111/2041-210X.13800

CrossRef Full Text | Google Scholar

Lazzaro, A., Hilfiker, D., and Zeyer, J. (2015). Structures of microbial communities in alpine soils: seasonal and elevational effects. Front. Microbiol. 6:1330. doi: 10.3389/fmicb.2015.01330

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Rui, J., Mao, Y., Yannarell, A., and Mackie, R. (2014). Dynamics of the bacterial community structure in the rhizosphere of a maize cultivar. Soil Biol. Biochem. 68, 392–401. doi: 10.1016/j.soilbio.2013.10.017

CrossRef Full Text | Google Scholar

Li, J., Shen, Z., Li, C., Kou, Y., Wang, Y., Tu, B., et al. (2018). Stair-step pattern of soil bacterial diversity mainly driven by pH and vegetation types along the elevational gradients of Gongga Mountain, China. Front. Microbiol. 9:569. doi: 10.3389/fmicb.2018.00569

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipson, D. A., and Schmidt, S. K. (2004). Seasonal changes in an alpine soil bacterial community in the Colorado Rocky Mountains. Appl. Environ. Microbiol. 70, 2867–2879. doi: 10.1128/AEM.70.5.2867-2879.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Logares, R., Audic, S., Bass, D., Bittner, L., Boutte, C., Christen, R., et al. (2014). Patterns of rare and abundant marine microbial eukaryotes. Curr. Biol. 24, 813–821. doi: 10.1016/j.cub.2014.02.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Lomolino, M. V. (2001). Elevation gradients of species-density: historical and prospective views. Glob. Ecol. Biogeogr. 10, 3–13. doi: 10.1046/j.1466-822x.2001.00229.x

CrossRef Full Text | Google Scholar

Looby, C. I., and Martin, P. H. (2020). Diversity and function of soil microbes on montane gradients: the state of knowledge in a changing world. FEMS Microbiol. Ecol. 96:fiaa 122. doi: 10.1093/femsec/fiaa122

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Z., Liu, H., Mi, Z., Zhang, Z., Wang, Y., Xu, W., et al. (2017). Climate warming reduces the temporal stability of plant community biomass production. Nat. Commun. 8:15378. doi: 10.1038/ncomms15378

PubMed Abstract | CrossRef Full Text | Google Scholar

Magoč, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, P. H., and Bellingham, P. J. (2016). Towards integrated ecological research in tropical montane cloud forests. J. Trop. Ecol. 32, 345–354. doi: 10.1017/S0266467416000432

CrossRef Full Text | Google Scholar

McHugh, T. A., Koch, G. W., and Schwartz, E. (2014). Minor changes in soil bacterial and fungal community composition occur in response to monsoon precipitation in a semiarid grassland. Microb. Ecol. 68, 370–378. doi: 10.1007/s00248-014-0416-3

CrossRef Full Text | Google Scholar

Nemergut, D. R., Costello, E. K., Hamady, M., Lozupone, C., Jiang, L., Schmidt, S. K., et al. (2011). Global patterns in the biogeography of bacterial taxa. Environ. Microbiol. 13, 135–144. doi: 10.1111/j.1462-2920.2010.02315.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Nemergut, D. R., Schmidt, S. K., Fukami, T., O'Neill, S. P., Bilinski, T. M., Stanish, L. F., et al. (2013). Patterns and processes of microbial community assembly. Microbiol. Mol. Biol. Rev. 77, 342–356. doi: 10.1128/MMBR.00051-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Ning, D., Yuan, M., Wu, L., Zhang, Y., Guo, X., Zhou, X., et al. (2020). A quantitative framework reveals ecological drivers of grassland microbial community assembly in response to warming. Nat. Commun. 11:4717. doi: 10.1038/s41467-020-18560-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Nottingham, A. T., Fierer, N., Turner, B. L., Whitaker, J., Ostle, N. J., McNamara, N. P., et al. (2018). Microbes follow Humboldt: temperature drives plant and soil microbial diversity patterns from the Amazon to the Andes. Ecology 99, 2455–2466. doi: 10.1002/ecy.2482

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedrós-Alió, C. (2012). The rare bacterial biosphere. Annu. Rev. Mar. Sci. 4, 449–466. doi: 10.1146/annurev-marine-120710-100948

CrossRef Full Text | Google Scholar

Prado, A. G. S., and Airoldi, C. (1999). The influence of moisture on microbial activity of soils. Thermochim. Acta 332, 71–74. doi: 10.1016/S0040-6031(99)00062-3

CrossRef Full Text | Google Scholar

Rahbek, C. (1995). The elevational gradient of species richness: a uniform pattern? Ecography 18, 200–205. doi: 10.1111/j.1600-0587.1995.tb00341.x

CrossRef Full Text | Google Scholar

Rahbek, C. (2005). The role of spatial scale and the perception of large-scale species-richness patterns. Ecol. Lett. 8, 224–239. doi: 10.1111/j.1461-0248.2004.00701.x

CrossRef Full Text | Google Scholar

Rui, J., Hu, J., Wang, F., Zhao, Y., and Li, C. (2022). Altitudinal niches of symbiotic, associative and free-living diazotrophs driven by soil moisture and temperature in the alpine meadow on the Tibetan plateau. Environ. Res. 211:113033. doi: 10.1016/j.envres.2022.113033

PubMed Abstract | CrossRef Full Text | Google Scholar

Rui, J., Li, J., Wang, S., An, J., Liu, W. T., Lin, Q., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, C., He, J.-Z., and Ge, Y. (2021). Seasonal dynamics of soil microbial diversity and functions along elevations across the treeline. Sci. Total Environ. 794:148644. doi: 10.1016/j.scitotenv.2021.148644

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, C., Ni, Y., Liang, W., Wang, J., and Chu, H. (2015). Distinct soil bacterial communities along a small-scale elevational gradient in alpine tundra. Front. Microbiol. 6:582. doi: 10.3389/fmicb.2015.00582

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, C., Shi, Y., Fan, K., He, J.-S., Adams, J. M., Ge, Y., et al. (2019). Soil pH dominates elevational diversity pattern for bacteria in high elevation alkaline soils on the Tibetan plateau. FEMS Microbiol. Ecol. 95:fiz 003. doi: 10.1093/femsec/fiz003

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, C., Xiong, J., Zhang, H., Feng, Y., Lin, X., Li, X., et al. (2013). Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biol. Biochem. 57, 204–211. doi: 10.1016/j.soilbio.2012.07.013

CrossRef Full Text | Google Scholar

Shi, Y., Xu, M., Zhao, Y., Cheng, L., and Chu, H. (2022). Soil pH determines the spatial distribution, assembly processes and co-existence networks of microeukaryotic community in wheat fields of the North China plain. Front. Microbiol. 13:911116. doi: 10.3389/fmicb.2022.911116

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, D., Lee-Cruz, L., Kim, W.-S., Kerfahi, D., Chun, J.-H., and Adams, J. M. (2014). Strong elevational trends in soil bacterial community composition on Mt. Halla, South Korea. Soil Biol. Biochem. 68, 140–149. doi: 10.1016/j.soilbio.2013.09.027

CrossRef Full Text | Google Scholar

Stegen, J. C., Lin, X., Fredrickson, J. K., and Konopka, A. E. (2015). Estimating and mapping ecological processes influencing microbial community assembly. Front. Microbiol. 6:370. doi: 10.3389/fmicb.2015.00370

PubMed Abstract | CrossRef Full Text | Google Scholar

Stegen, J. C., Lin, X., Konopka, A. E., and Fredrickson, J. K. (2012). Stochastic and deterministic assembly processes in subsurface microbial communities. ISME J. 6, 1653–1664. doi: 10.1038/ismej.2012.22

PubMed Abstract | CrossRef Full Text | Google Scholar

Steven, B., Pollard, W. H., Greer, C. W., and Whyte, L. G. (2008). Microbial diversity and activity through a permafrost/ground ice core profile from the Canadian high Arctic. Environ. Microbiol. 10, 3388–3403. doi: 10.1111/j.1462-2920.2008.01746.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sundqvist, M. K., Sanders, N. J., and Wardle, D. A. (2013). Community and ecosystem responses to elevational gradients: processes, mechanisms, and insights for global change. Annu. Rev. Ecol. Evol. Syst. 44, 261–280. doi: 10.1146/annurev-ecolsys-110512-135750

CrossRef Full Text | Google Scholar

Tourna, M., Stieglmeier, M., Spang, A., Könneke, M., Schintlmeister, A., Urich, T., et al. (2011). Nitrososphaera viennensis, an ammonia oxidizing archaeon from soil. Proc. Natl. Acad. Sci. U. S. A. 108, 8420–8425. doi: 10.1073/pnas.1013488108

PubMed Abstract | CrossRef Full Text | Google Scholar

Unger, I. M., Kennedy, A. C., and Muzika, R.-M. (2009). Flooding effects on soil microbial communities. Appl. Soil Ecol. 42, 1–8. doi: 10.1016/j.apsoil.2009.01.007

CrossRef Full Text | Google Scholar

Wang, Z., Brown, J. H., Tang, Z., and Fang, J. (2009). Temperature dependence, spatial scale, and tree species diversity in eastern Asia and North America. Proc. Natl. Acad. Sci. U. S. A. 106, 13388–13392. doi: 10.1073/pnas.0905030106

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, G., Hu, Y., Wang, S., Zhang, Z., Chang, X., Duan, J., et al. (2010). Effects of litter quality and climate change along an elevation gradient on litter mass loss in an alpine meadow ecosystem on the Tibetan plateau. Plant Ecol. 209, 257–268. doi: 10.1007/s11258-009-9714-0

CrossRef Full Text | Google Scholar

Xu, M., Li, X., Cai, X., Gai, J., Li, X., Christie, P., et al. (2014). Soil microbial community structure and activity along a montane elevational gradient on the Tibetan plateau. Eur. J. Soil Biol. 64, 6–14. doi: 10.1016/j.ejsobi.2014.06.002

CrossRef Full Text | Google Scholar

Yang, Y., Gao, Y., Wang, S., Xu, D., Yu, H., Wu, L., et al. (2014). The microbial gene diversity along an elevation gradient of the Tibetan grassland. ISME J. 8, 430–440. doi: 10.1038/ismej.2013.146

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Ning, D., Yang, Y., He, N., Li, X., Cornell, C. R., et al. (2022). Precipitation balances deterministic and stochastic processes of bacterial community assembly in grassland soils. Soil Biol. Biochem. 168:108635. doi: 10.1016/j.soilbio.2022.108635

CrossRef Full Text | Google Scholar

Yang, Y.-D., Ren, Y.-F., Wang, X.-Q., Hu, Y.-G., Wang, Z.-M., and Zeng, Z.-H. (2018). Ammonia-oxidizing archaea and bacteria responding differently to fertilizer type and irrigation frequency as revealed by Illumina Miseq sequencing. J. Soils Sediments 18, 1029–1040. doi: 10.1007/s11368-017-1792-3

CrossRef Full Text | Google Scholar

Yao, M., Rui, J., Niu, H., Heděnec, P., Li, J., He, Z., et al. (2017). The differentiation of soil bacterial communities along a precipitation and temperature gradient in the eastern Inner Mongolia steppe. Catena 152, 47–56. doi: 10.1016/j.catena.2017.01.007

CrossRef Full Text | Google Scholar

Yuan, C., Li, F., Yuan, Z., Li, G., and Liang, X. (2021). Response of bacterial communities to mining activity in the alpine area of the Tianshan Mountain region, China. Environ. Sci. Pollut. Res. 28, 15806–15818. doi: 10.1007/s11356-020-11744-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, Y., Si, G., Wang, J., Luo, T., and Zhang, G. (2014). Bacterial community in alpine grasslands along an altitudinal gradient on the Tibetan plateau. FEMS Microbiol. Ecol. 87, 121–132. doi: 10.1111/1574-6941.12197

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, X.-Q., and Zhou, X.-M. (1999). Ecological basis of alpine meadow ecosystem management in Tibet: Haibei alpine meadow ecosystem research station. Ambio 28, 642–647.

Google Scholar

Zhou, J., Deng, Y., Shen, L., Wen, C., Yan, Q., Ning, D., et al. (2016). Temperature mediates continental-scale diversity of microbes in forest soils. Nat. Commun. 7:12083. doi: 10.1038/ncomms12083

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J., and Ning, D. (2017). Stochastic community assembly: does it matter in microbial ecology? Microbiol. Mol. Biol. Rev. 81, e00002–e00017. doi: 10.1128/MMBR.00002-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, D., Zhao, L., Sheng, Y., Chen, J., Hu, G., Wu, T., et al. (2017). A new map of permafrost distribution on the Tibetan plateau. Cryosphere 11, 2527–2542. doi: 10.5194/tc-11-2527-2017

CrossRef Full Text | Google Scholar

Keywords: soil microbial community, altitudinal gradient, temporal dynamic, alpine grassland, Actinobacteria, community assembly

Citation: Rui J, Zhao Y, Cong N, Wang F, Li C, Liu X, Hu J, Ling N and Jing X (2023) Elevational distribution and seasonal dynamics of alpine soil prokaryotic communities. Front. Microbiol. 14:1280011. doi: 10.3389/fmicb.2023.1280011

Received: 19 August 2023; Accepted: 06 September 2023;
Published: 22 September 2023.

Edited by:

Xingjia Xiang, Anhui University, China

Reviewed by:

Shuo Jiao, Northwest A&F University, China
Yue Lu, Hunan University, Changsha, China
Congcong Shen, Chinese Academy of Sciences (CAS), China

Copyright © 2023 Rui, Zhao, Cong, Wang, Li, Liu, Hu, Ling and Jing. 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: Xin Jing, amluZ3hAbHp1LmVkdS5jbg==

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.