Driving force of soil microbial community structure in a burned area of Daxing’anling, China

Fires are an important factor impacting forest ecosystems of Daxing’anling and have a significant effect on soil microbial community structure. In this study, high-throughput sequencing for 16S rDNA and ITS rDNA were applied to analyze the changing characteristics and driving factors of bacterial and fungal community structures in burned areas with different fire severity. PICRUSt2 software was used to predict the functional characteristics of burned areas with different fire severity. The purpose was to unveil the responsive relationships among the structure and function of bacterial and fungal communities, fire severity, and post-disturbance restoration times. After high severity fires, the destruction of surface vegetation and loss of soil nutrients reduced the diversity and abundance of soil bacteria and fungi. The soil bacteria community structure, which was dominated by Acidobacteria, Proteobacteria, and Actinobacteria, changed to be dominated by Proteobacteria and Chloroflexi. As well, soil fungal community changed from domination by Helotiales, Eurotiales and Russulales to domination by Archaeorhizomycetales and Helotiales. Over time, soil bacterial community was gradually restored to pre-fire levels 30 years after the fire. Soil fungal community changed and failed to restore to pre-fire levels after 30 years. After low/intermediate severity fires, environmental factors were relatively unchanged so that soil bacteria diversity and abundance increased, optimizing community composition. The diversity and abundance of soil fungi decreased and the community structure changed slightly. Over time, both bacterial and fungal communities were gradually restored to pre-fire levels 30 years after the fire. After fire disturbance, with increasing severity, soil carbon fixation, lignin degradation, mineralization of organic nitrogen and hydrolysis of organic phosphorus are enhanced. Denitrification is weakened. Therefore, forest fires have certain positive effects on carbon, nitrogen and phosphorus cycles where soil bacteria and fungi are involved.


Introduction
Soil microbes have great complexity and diversity and are important for soil composition and nutrient cycling, and drivers of energy flow (Elsas and Boersma 2011) microbes play a critical role in plant growth and development, and in maintaining the productivity, functionality and stability of ecosystems (Felske et al. 2000;Kuramae et al. 2010). Since soil microbes are sensitive to changes in environment, they may be taken as indicators of soil restoration (Luo et al. 2020).
Abstract Fires are an important factor impacting forest ecosystems of Daxing'anling and have a significant effect on soil microbial community structure. In this study, highthroughput sequencing for 16S rDNA and ITS rDNA were applied to analyze the changing characteristics and driving factors of bacterial and fungal community structures in burned areas with different fire severity. PICRUSt2 software was used to predict the functional characteristics of burned areas with different fire severity. The purpose was to unveil the responsive relationships among the structure and function of bacterial and fungal communities, fire severity, and post-disturbance restoration times. After high severity fires, the destruction of surface vegetation and loss of soil nutrients reduced the diversity and abundance of soil bacteria and fungi. The soil bacteria community structure, which was dominated by Acidobacteria, Proteobacteria, and Actinobacteria, changed to be dominated by Proteobacteria and Chloroflexi. As well, soil fungal community changed Significant changes could result from either natural or anthropgenic disturbances to ecosystems, for which different types of soil microbes display divergent coping mechanism due to their contrasting survival characteristics and environments, thereby further changing the microbial community structure (Luo et al. 2020). Coal mining activities on the Loess Plateau have caused soil sedimentation which reduced soil bacterial diversity and resulted in significant changes in community structure (Luo et al. 2020). The removal of undergrowth vegetation caused chemical property changes in forest soils, including water contents, total nitrogen and phosphorus, hydrolyzed nitrogen, and available phosphorus, which in turn caused the reduction of relative abundance of soil bacteria of Proteobacteria, Planctomycetes, and Firmicutes, as well as a reduction in bacterial community diversity (Wu et al. 2020). Appropriate applications of nitrogen fertilizer increased the fertility of rice field soils but did not significantly change bacterial community structure (Liu et al. 2019). Fertilization in a short period increased available nutrients and soil water of a Eucalyptus forest, and the soil fungal community structure, which was dominated by Eurotiales, Archaeorhizomycetales, and Tremellales, changed to be dominated by Boletus and Eurotiales . When forest soils were degraded, the surface vegetation changed from trees to herbs, and the soil became more permeable but with lower nutrients such as total nitrogen and total carbon. In this study, the relative abundance of Eurotiales and Acidobacteria soil microbes increased dramatically but Basidiomycota levels were reduced significantly, indicating an obvious enhancement of microbe community diversity .
Among various disturbances, fires have quite different effects on soil microbial community structure which are a continuous process (Villadas et al. 2019;Luo et al. 2020). Microbial communities are first directly affected by high temperatures and indirectly by changes in physical and chemical properties (Cai et al. 2014;Wang et al. 2020). Change in organic matter quantity and quality, fluctuations in carbon and nitrogen contents (Shen et al. 2013;Kaiser et al. 2016), increases in pH (Holden et al. 2016;Wang et al. 2016;Schmidta et al. 2019), and variations in water contents are key factors changing the microbial communities (Shigeto et al. 2008;Man et al. 2010;Rasche et al. 2011;Vries 2012). Forest fires have two forms of impacts on microbial communities (Brais et al. 2000;Wallenius et al. 2004): one, the structure and functions of microbial communities can be destroyed under high severity fires, while, low severity fires optimize the forming of microbial communities, contributing to the stability of ecosystem structure (Akburak et al. 2018; Van et al. 2019).
In Daxing'anling forest, boreal and the mixed wood species are widely distributed (Tian et al. 2018). This is a region of high fire risk. In Hulunbuir City where the study area is located, fires occurred over 600 times from 1990 to 2018. Therefore, it is necessary to study the microbial restoration characteristics of the forest soils of the Daxing'anling area after fire. High-throughput sequencing for 16S rDNA and ITS rDNA were applied to analyze the changing characteristics and relative driving factors of the soil bacterial and fungal communities in burned areas under different fire severities. PICRUSt2 software was used to predict the functional characteristics of the soil element cycle in burned areas under different fire severity. The purpose of this study was to determine the relationship between structure and function of soil bacterial and fungal communities, the fire severity, and post-disturbance restoration time, to provide information for understand the restoration process of soil structure and function, and to propose targeted management measures.

Study area
The research was located in Hulunbuir City, northwest of Daxing'anling in the northernmost part of China. It is a humid temperate zone with average annual temperatures of − 5.4 °C, and average annual precipitation of 500 mm, mainly concentrated in June to September, and average altitudes of 976.5 m a.s.l. It has an average 75% forest cover consisting of Larix gmelinii associated with Betula platyphylla, Populus davidiana and Pinus sylvestris (Tian et al. 2018). The study area has a high incidence of fires, over 677 fires from 1990 to 2018. Fires generally occur in May and June . For this study, based on consistent site conditions, five burned blocks from the southern parts of Genhe City were selected in August 2018, where fires happened all by lightning (Table 1, Fig. S1) in May 1987, June 2003, May 2010, May 2015, and June 2018. These burned areas characterize the recovery status of the forest ecosystem after 30 years, 15 years, 8 years, 3 years and 8 weeks, respectively, following the fire.

Setting of plots
Combining the forest resource features of Daxing'anling, (evolution history, forest age, tree species and forest type), the Composite Burn Index Investigation Contents and Evaluation Criteria (Table 2), which is applicable to the study area, was prepared according to the Composite Burn Index (CBI) proposed by Key and Benson (Key and Benson 2006). The field survey was conducted in August 2018. Sixty sampling sites (over 50 m intervals) were randomly established in every burned area based on the CBI of all sites and classified based on fire severity index (FSI) (Table 3), which was considered the FSI classification for burned areas of Daxing'anling. By calculating the CBI and grading the FSI of the burned areas created in different times, one low severity plot, one intermediate severity plot and one high severity plot were selected in each block. An unburned plot near the burned area was also selected which had the same site condition. There were 20 plots in total.

Sample collection
Three sampling points were randomly selected in each plot. With the sampling point as center, and within 3-m diameter, an auger with a 3-cm inner diameter was used to collect 10 samples at 1-10 cm depths. The samples were mixed and screened through a 2 mm-sieve, packed in self-sealing bags and stored at − 20 °C. These samples were used for measuring the chemical properties, bacterial diversity, and fungal diversity of the soil.

Measurement of soil chemical properties
An aqueous extract of potassium chloride solution was used to determine pH (measurement standard: NY/T 1377-2007); an potassium dichromate volumetric method to measure total organic carbon (TOC) (measurement standard: NY/T 1121.6-2006); the Kjeldahl distillation was used to measure total nitrogen (TN) (measurement standard: NY/T 53-1987); the sodium hydroxide fusion-molybdenum stibium resist colorimetric method was used to measure total phosphorus (TP) (measurement standard: NY/T 88-1988) and total potassium (TK) (measurement standard: NY/T 87-1988); the alkaline hydrolysis diffusion method was used to measure available nitrogen (AN) (measurement standard: DB13T 843-2007); and, the sodium bicarbonate extraction-molybdenum stibium resist colorimetric method was used to measure available phosphorus (AP) (measurement standard: NY/T 1121.7-2014) ( Table 4).

DNA extraction and PCR amplification
Total DNA of the microbial communities were extracted according to the instructions of E.Z.N.A. ® soil DNA kit (Omega Bio-tek, Norcross, GA, USA.). A 1% agarose gel electrophoresis was applied to detect the DNA extraction quality while the NanoDrop2000 determined the DNA  concentration and purity. For bacterial samples, the 338F (5′-ACT CCT ACG GGA GGC AGC AG-3′) and 806R (5′-GGA CTA CHVGGG TWT CTAAT-3′) were used to carry out PCR amplification for the V3-V4 variable region of 16S rRNA gene; for fungal samples, ITS1 (5′-TCC GTA GGT GAA CCT GCG G-3′) and ITS4 (5′-TCC TCC GCT TAT TGA TAT GC-3′) were applied to perform PCR amplification for the ITS region of rRNA gene. Each operation was repeated three times on each sample.

High-throughput sequencing
The PCR products of the same sample were mixed and recycled using 2% agarose gel, and then purified the recycled products by AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Ltd, Union City, CA, USA). A 2% agarose gel electrophoresis was carried out and the recycled products quantified by Quantus™ Fluorometer (Promega, Ltd, Madison City, WI, USA). NEXTFLEX Rapid DNA-Seq Kit was next applied to set up a library and sequencing was performed by Miseq PE300 platform (Illumina, Ltd, San Diego City, CA, USA) of the company. The original data were uploaded to the NCBI SRA (Sequence Read Archive) database.

Biological information extraction
Trimmomatic software (version 0.39, http://www.usade llab. org/cms/?page=trimm omati c) was applied to control the quality of original sequencing, and the FLASH software used for splicing. UPARSE software (version 7.1, http:// drive 5.com/upars e/) was used to carry out OTU clustering and to eliminate chimera on sequence according to 97% similarity. The RDP classifier (http://rdp.cme.msu.edu/) was then used to carry out species classification and annotation for each sequence, and the bacterial and fungal samples were compared with Silva database (Release132, http://www.arbsilva .de) and Unite database (Release 7.2, http://unite .ut.ee/ index .php), respectively, for which the comparison threshold was set as 70%.

Data processing
SPSS (IBM, USA) software was used to calculate the mean and standard deviation of each chemical factor and to perform significance analysis. The mothur software package was used to calculate the diversity index and the R programming language (AT&T Bell Laboratories, Ltd, Hurray Hill City, NJ, USA) vegan package to carry out redundancy analysis (RDA) and canonical correspondence analysis (CCA)

Chemical properties of soils
The mixed forest soil of Daxing'anling has pH 4.5-5.2 with total organic carbon (TOC) content of 23.2-80.6 g kg −1 , total nitrogen (TN) of 2.2-6.9 g kg −1 , and total phosphorus (TP) of 0.2-0.5 g kg −1 (Table 4). Eight weeks after fire in soils subjected to intermediate/high severity burning, the TOC, TN, and AP contents dropped significantly, but soil pH rose dramatically. TP, TK, and AN levels increased slightly. In soil under a low severity fire, pH rose slightly and TOC, TN, TP, TK, AN, and AP levels increased slightly. Fifteen years after the fire, TP, TK, and AN contents were restored to pre-fire levels. Thirty years after fire, TOC and TN contents as well as pH were at pre-fire levels. Over time, the AP content gap between the burned area and the control plot increased significantly (Table 4).

Structural characteristics
The 16SrRNA high-throughput sequencing was carried out on all 60 soil samples of 20 groups. There were 3,302,646 high quality sequences with 432 bp mean reading length obtained after quality control and chimera removed, and were clustered into 4771 OTUs (97% similarity), belonging to 33 phylum, 234 orders, and 604 genera.

Bacterial diversity characteristics
Chao and Shannon indexes, respectively, characterize the abundance ratio and diversity of bacterial species (Fig. 2). Eight weeks after the fire, both the richness index and diversity indexes were reduced dramatically; 3 years after the fire, the two indexes increased beyond pre-fire levels; 15 years later, the two indexes dropped to pre-fire levels. The β diversity characteristics of soil bacteria in all plots is shown by Fig. 1 Soil bacterial community composition at the phylum level RDA analysis. All plots are located in five agglomeration areas which take the sampling site as the center (Fig. 3). Soil bacteria are capable of resisting disturbance and able to maintain a relatively stable community structure and restore rapidly after fire.

Influence of environmental factors on bacterial communities
From the redundancy analysis (RDA) (Fig. 3), it is clear that the explained variable in the first and second axes are 17.7% and 7.7%, respectively, and the cumulative contribution rate is 25.45% (p < 0.01). Soil pH, TOC, and AP are explained by 19.31%, 18.21%, and 17.69% of the variance, respectively.
The cover vegetation has less effect on soil bacterial community (Fig. 3); the explained variable in the first and second axes are 17.75% and 7.7%, respectively, and the cumulative contribution rate is 25.45% (p < 0.01).

Fungal diversity characteristics
The soil fungal diversity index (value range: 1.5-3.5) was significantly lower than the soil bacterial diversity index (value range: 5.0-6.5) (Fig. 5). After fire, fungal richness and diversity indexes were dramatically reduced. With increased recovery time, abundance ratios and diversity index gradually recovered; the abundance ratio was restored to pre-fire levels after 30 years while the diversity index did not (Fig. 5).
From the β diversity characteristics of soil fungi, it is apparent that: the high severity burned and the low/intermediate severity burned areas have fungi in two clusters (Fig. 6). Fungi are more sensitive to environmental changes than bacteria after fire. This is because soil fungi have a homogeneous structure and are mainly distributed over a few predominant phyla. They are unstable and have poor adaptability to disturbance, and easily show significant changes in community structure after fire, requiring a length period for restoration.

Influences of environmental factors on fungal communities
The results of the canonical correspondence analysis (CCA) (Fig. 6) show that soil factors have little influence on the fungal community. The relative explained variable in the first and second axes are 7.1% and 6.4%, respectively. The cumulative contribution rate is 13.6% (p < 0.01). Soil pH, TOC and AP contents explain 9.0%, 8.9% and 6.7% variance, respectively, which indicates that these are the main factors affecting soil fungal communities. The vegetation factors such as tree biomass, herbaceous biomass, and herbaceous diversity can greatly affect the fungal community. The explained variable in the first and second axes are 14.4% and 11.9%, respectively. The cumulative contribution rate is 26.4% (p < 0.01). Tree biomass, herbaceous biomass, and herbaceous diversity respectively explain 18.7%, 16.0%, and 15.6% variance, respectively, showing that they are the main influencing factors on soil fungal communities (Fig. 6).

Prediction on soil functions
PICRUSt2 was applied to identify the functional characteristics of soil microbes related to element cycling. Among the functions that the soil bacteria are involved in, the dominating genes were those relating to the anaerobic acetyl-CoA pathway and the 3-hydroxypropionic acid pathway in carbon fixation, and those involved in cellulose degradation, and organic phosphorus hydrolysis. Relative  The analysis on carbon cycling functions showed that: after fire, with increasing severity, Calvin fixation was enhanced to an ever-increasing extent. Among them, the relative abundance ratios (relative abundance ratio refers to the ratio of the abundance of the soil enzyme to the sum of the abundance of all enzymes in the soil) of 1,5-ribulose bisphosphate carboxylase [EC: 4.  (Fig. 7).
Among the processes that soil fungi are involved in, the dominating genes are those relating to the oxaloacetate pathway and organic phosphorus hydrolysis in carbon fixation, while other functions were comparatively low abundances. The analysis of carbon cycling showed that with increasing fire severity, the relative abundance ratios of pyruvate carboxylase [EC:6.4.  (Fig. 8).

Characteristics and driving factors of soil bacterial community structure after fire
The soils of the Daxing'anling mixed forest contained similar bacterial communities, dominated by Acidobacteria, Proteobacteria, and Actinobacteria. After high severity fire disturbance, the two bacterial indexes of abundance ratio and diversity both declined dramatically. The structure also changed considerably to a community dominated by Proteobacteria and Chloroflexi. Over time, the indexes of abundance and diversity increased significantly to pre-fire levels 30 years after the fire. After low/intermediate-severity fire, Acidobacteria decreased significantly, while other bacteria increased slightly. Over time, the bacterial community which was originally dominated as noted above, was dominated by Proteobacteria, Actinobacteria and Chloroflexi. Mild fire disturbance only resulted in small changes to the soil bacteria.
When the environment changes, different bacteria show different coping strategies (k-strategy, r-strategy) due to their different survival characteristics. Acidobacteria, Actinobacteria, and Firmicutes are gram-positive bacteria (Fierer et al. 2007), which follow a "k-strategy (Kapazitätsgrenzestrategistis)" for survival, and were characterized by good growth in a stable environment but with a weak capability in adapting to a changed environment (Mickan et al. 2017); Proteobacteria, Gemmatimonadetes, and most Chloroflexi gram-negative bacteria (Chaudhry et al. 2012) with "r-strategy (rate-strategistis)" for survival, which features rapid growth and capable of adapting to new environment (Mickan et al. 2017). Fontaine et al. (2003) considered that fast-growing "r-strategy" species dominated in the early stages of environmental change, and later, these species were replaced by "k-strategy" species when easy-to-degrade components of nutrition were depleted. In the study area, the soil environment was relatively stable. The "k-strategy" species like Acidobacteria and Actinobacteria dominated the soil bacteria communities. After fire disturbance, nutrition from organic matter were greatly reduced (Table 4), and the bacteria species changed. In this event, "r-strategy" species such as Proteobacteria and Chloroflexi, multiplied and became dominant. Then, when the soil environment tended to be stable again, "k-strategy" species gradually resumed a dominating status.
Among the environmental factors, vegetation has little effect on soil bacteria; pH, TOC (total organic carbon) and AP (available phosphorous) have obvious effects. Soil pH was negatively correlated with the relative abundance of Acidobacteriales of the Acidobacteria phylum, Actinobacteria and Corynebacteriales of the Actinobacteria phylum, and Gemmatales and Isosphaerales of the Planctomycetes phylum It was positively correlated with Gaiellales and IMCC26256 of the Actinobacteria phylum, Betaproteobacteriales, Elsterales, Myxococcales and Rhizobiales of the Proteobacteria phylum, AD3 and KD4-96 of the Chloroflexi phylum, Chloroplast of the Cyanophyta phylum, Nitrospira of the Nitrospirae phylum, Chitinophagales of the Bacteroidetes phylum, and Chthoniobacterales of the Verrucobacteria phylum (Fig. 9). Austin and Smith (1990) found that pH had a direct physiological gradient on organisms, i.e., the growth of organisms is subject to the restraining level of pH. The ability of different organisms to adapt to low pH is an acid tolerance response (ATR). The optimal pH for Acidobacteria growth is 4-5, having high ATR, but most bacteria grow better at a pH 7-8 environment, having comparatively low acid tolerance. Soil pH in the study area was 4.2-5.2 so that, compared to other bacteria, Acidobacteriales and other bacteria of the same kind, which possessed higher ATR, were more competitive. But due to the alkaline hydrolysis effect of ash contents after fire, pH rose slightly to 4.6-5.6. At this pH, the competitiveness of Acidobacteria was weakened while the relative abundance of other bacteria increased accordingly.
Soil TOC (total organic carbon) was positively correlated with the relative abundance of Subgroup_11 of the Acidobacteria phylum, Gaiellales and Corynebacteriales of the Actinobacteria phylum, and Bacillales of the Firmicutes phylum (Fig. 9). Soil TOC is the main carbon source for heterotrophic microorganisms such as Acidobacteria, Actinobacteria, and Firmicutes. After fire, the loss of soil TOC showed significant depressing effects on these bacteria.
Available soil phosphorous was positively correlated with the relative abundance ratios of Acidobacteriales of the Acidobacteria phylum, Solirubrobacterales and Actinobacteria of the Actinobacteria phylum, Bacillales of the Firmicutes phylum, and Gemmatales of the Planctomycetes phylum, but negatively correlated to KD4-96 of the Chloroflexi phylum, Chloroplast of the Cyanobacteria phylum, and Elsterales and Myxococcales of the Proteobacteria phylum ( Fig. 9). Phosphorus is an important element of biological cells. Available soil phosphorous is directly absorbed and utilized by soil microorganisms, so that its level in the soil has an important influence on the growth of microorganisms. Total phosphorus (0.21-0.49 g kg −1 ) is lower than the total phosphorous (0.41-1.5 g kg −1 ) of the upper 10-cm layer in generally all Chinese forest soils (Gholz et al. 2000;Mcgroddy et al. 2004;Han et al. 2005;Parton et al. 2007;Li et al. 2018) which confirms that soils of the study area lack phosphorus. In the meantime, after fire disturbance and due to appropriate environmental conditions, bacteria such as Chloroflexi, Cyanobacteria, and Proteobacteria multiplied and consumed any available phosphorous and allowed for the production of weak bacteria like Acidobacteria and Actinobacteria (Parton et al. 2007).
According to previous studies, total soil nitrogen and available nitrogen are strongly correlated to Planctomycetes, Firmicuates, and Proteobacteria bacteria (Mao and Qi 2015). However, the results of this study show that the levels of total and available nitrogen had little influence on bacteria, which is different from previous research findings. Total soil nitrogen (2.24-6.93 g kg −1 ) in the study area is relatively higher than the total nitrogen levels (1.2-6.35 g kg −1 ) of the upper 10 cm layer in Chinese forests (Gholz et al. 2000;Mcgroddy et al. 2004;Han et al. 2005;Parton et al. 2007;Li et al. 2018), which indicates that the soils of the study area are not lacking in nitrogen; after the fire, however, total soil Fig. 9 Correlation coefficients of soil bacterial community and environmental factors (**represents a significant difference between the factors at the 0.01 level; *a significant difference at the 0.05 level) nitrogen decreased, the available nitrogen which is directly absorbed by organisms increased. Therefore, the level of nitrogen resources was not restricted by soil bacteria.

Characteristics and driving factors of fungal community structure after fire
The dominant soil fungi orders are Helotiales, Agaricales, and Russulales. Eight weeks after a high intensity fire, the abundance ratios and diversity decreased significantly and the community structure showed significant change from the original fungi to fungi dominated by Archaeorhizomycetales, Helotiales, and Eurotiales. Over time, the community structure enlarged. Thirty years after fire, the fungal community changed to a community dominated by Archaeorhizomycetales, Russulales, and Helotiales. After low/intermediate fire, the abundance ratios and diversity decreased but the community structure was only slightly changed from being dominated by Helotiales, Agaricales, and Russulales to being dominated by Helotiales. With greater recovery time, the bacterial community gradually recovered and was restored to pre-fire levels.
Different species of soil fungi have various survival characteristics and specific relationships with vegetation. Therefore, changes in vegetative cover could be a significant influence on soil fungi. Archaeorhizomycetales was negatively correlated with tree biomass, and positively with herbaceous diversity. Eurotiales and Helotiales are negatively correlated with tree biomass while positive with herbaceous diversity. Agaricales and Russulales show positive correlation with biomass and negative correlation with herbaceous diversity (Fig. 10). Agaricomycetes of the Agaricales phylum and Russula of the Russulales phylum are the main ectomycorrhizal fungi of pines and also able to use saprophytic fungi using litter (Zhang et al. 2016). Archaeorhizomycetales is commonly distributed in the rhizosphere of Pinus, Rhododendron, and other species which can either coexist with mycorrhizal fungi or provide organic matter for saprophytic fungi (Rosling et al. 2011). The mixed forest of Daxing'anling is largely homogeneous plant species dominated by Chinese larch associated with birch and aspen. Their roots are densely distributed in the soil which is rich in organic matter. Agaricales and Russulales fungi multiply rapidly through symbiosis and saprophytic relationship with roots. After fire, most vegetation was killed and roots decreased accordingly. Soil organic matter was largely consumed and further reduced the quantity Agaricales and Russulales due to the loss of nutrients. For low/intermediateseverity fire, the surface vegetation structure had not been greatly changed and recovered quickly, so that the number of Agaricales and Russulales could gradually restore to prefire levels. For high-severity fire, the restoration process of surface vegetation should be: open herb community → shrub and grass community → mixed forest. Agaricales and Russulales are specific to ecosystems and host habitats and unable to adapt to frequent changes to surface vegetation and failed to restore to pre-fire level. However, Archaeorhizomycetales, due to its powerful adaption capability dominated among soil fungi (Fig. 10).
The influence of the soil environment on soil fungi cannot be ignored. Archaeorhizomycetales is positively correlated with pH, while Helotiales, Agaricales, and Russulales are negatively correlated (Fig. 10). Soil fungi also show a prominent physiological gradient to pH. The optimal pH for Agaricales and Russulales is 3-5, while the optimal for Archaeorhizomycetales is around 7 (Xiong et al. 2019). Archaeorhizomycetales was negatively correlated with available phosphorous and total organic carbon; Eurotiales was negatively correlated with available phosphorous; Russulales and Agaricales were positively correlated with total organic carbon; Helotiales was positively correlated with available phosphorous (Fig. 10). After fire disturbance, the pH rose which promoted the growth of Archaeorhizomycetales which grows best in a neutral soil. Organic matter in the soil was consumed by the fire so that the growth of Agaricales and Russulales were significantly reduced due to the lack of saprophytic nutrition. Soil available phosphorous was consumed and reduced by Archaeorhizomycetales and Eurotiales, which dramatically contained the growth of Helotiales.

Prediction of soil functions after fire disturbance
According to the prediction of PICRUSt2, soil microbe functions related to carbon fixation, cellulose degradation and organic phosphorus hydrolysis were in relatively high abundance. In contrast, hemicellulose and lignin degradation, nitrogen fixation, nitrification and denitrification were low. This indicates that, in the soil of the study area, microbes are capable in decomposing cellulose but weak in decomposing hemicellulose and lignin, and weak in functions related to nitrogen cycling.
The carbon cycling process in the soil occurs with both bacteria and fungi simultaneously. Carbon fixation by the soil converts inorganic carbon to organic carbon, thus enhancing total organic carbon. In the soil of the study area, the Calvin cycle pathway, involving participation of photoautotrophic microorganisms like Cyanobacteria, the 3-hydroxypropionate pathway involving input by Chloroflexi, the anaerobic acetyl-CoA pathway with anaerobic bacteria, and the oxaloacetate pathway involving of symbiotrophism fungi like ectomycorrhizal fungi and arbuscular mycorrhizal fungi, were the main ways to increase total organic carbon (Yuan et al. 2011). The degrading cellulose, hemicellulose and lignin represented microbe capability in decomposing and re-utilizing organic matter in animal and plant residues (Yuan et al. 2011), among which, Acidobacteria helped degrade plant carbon like cellulose and hemicellulose efficiently (Fierer et al. 2007). Actinobacteria played an important role in lignin degradation (Pan et al. 2011). After fire disturbance, the rise of soil pH led to an increase in the relative abundance of Cyanobacteria, Chloroflexi, and Acidobacteria, which further enhanced microbe carbon fixation and lignin degradation; the relative abundance of Acidobacteria decreased, which weakened microbe capability in degrading cellulose and hemicellulose. In fact, the enhanced capability in carbon fixation and lignin degradation were favorable for total organic carbon accumulation. The forest fire played a positive role in soil carbon cycling.
Soil nitrogen cycling was accomplished by soil bacteria. Soil nitrogen cycle mainly includes nitrogen fixation, nitrification, and denitrification, etc. Nitrogen fixation converts nitrogen in the air into ammonia, which is beneficial for the accumulation of nitrogen in the soil (Sierra 1992;Chen et al. 2006). Nitrogen fixation in the study area soil was achieved through symbiotic nitrogen fixation bacteria like Frankia. Nitrification is an important step for organic nitrogen mineralization, for soil microbes in decomposing and reutilizing nitrogen in animal and plant residues (Weier et al. 1993;Galloway et al. 2004). In the study area, nitrification is mainly achieved by nitrifying bacteria of Nitrospira and Nitrosomonadaceae. Denitrification converts nitrates and nitrites into nitrogen to return it to the atmosphere (Weier et al. 1993;Galloway et al. 2004). In the study area, denitrification is mainly achieved by denitrifying bacteria. Research indicates that most denitrifying bacteria are heterotrophic which take organic matter as carbon and energy sources (Mao and Qi2015). After fire, Frankia decreased due reduction in available phosphorous and weakened nitrogen fixation. The increase in pH resulted in an increase in nitrifying bacteria such as Nitrospira and Nitrosomonadaceae which enhanced nitrification. The decrease total organic carbon levels resulted in an inadequate supply of denitrifying bacteria, weakening the denitrification effect. The strengthening of soil nitrification and nitrogen fixation were important for the accumulation of nitrogen; the weakening of denitrification reduced nitrogen consumption and fire played a positive role in soil nitrogen cycling.
Organic phosphorus hydrolysis hydrolyzes macromolecular organic phosphorus into available phosphorus that can be directly absorbed by organisms (Podgorskii 1988;Liu et al. 2015). This occurs in both soil bacteria and fungi simultaneously. After fire, soil pH (4.6-5.6) rose slightly but was still weakly acidic. This had little effect on acid phosphatase but and weakened the inhibitory effect of alkaline phosphatase. Overall levels of phosphatase increased, the hydrolysis of organic phosphorus was enhanced, and the availability of phosphorus in the soil was also increased. Fire had positive influence on phosphorus cycling.

Conclusions
Fire of different intensities has various effects on soil bacterial communities. Severe fires negatively impact soil bacterial communities. Fire changes the nutrients levels, damages original bacterial communities, and replaces them by bacteria with more powerful adaption capabilities. Such kinds of changes require at least 30 years to restore to the original communities. Low/intermediate fires enhances diversity and abundances of bacterial communities and optimizes composition by improving the soil environment.
Changes brought by high intensity fires on soil fungal communities are irreversible in the short-term because surface vegetation is destroyed and there is a loss of soil nutrients. This results in the decrease or loss of dominant fungi like Agaricales and Russulales while Archaeorhizomycetales becomes dominant fungi by its strong adaptability. Over time, these changes become worse and are not restored in 30 years. However, low/intermediate-severity fires only affects the environment to a small extent. The relative amounts of Agaricales and Russulales are slightly reduced but are restored essentially after 30 years.
After fire disturbance, soil carbon fixation, lignin degradation, mineralization of organic nitrogen and hydrolysis of organic phosphorus are enhanced, while denitrification is weakened. This positively promotes soil carbon, nitrogen, and phosphorus cycles to some extent.