Skip to main content

Soil microbial functional diversity responses to different vegetation types in the Heilongjiang Zhongyangzhan Black-billed Capercaillie Nature Reserve



The soil microbial community is an important bioactive component of terrestrial ecosystems. Its structural and functional diversity directly affects carbon and nitrogen processes. This study aimed to investigate the variations in the functional diversity of soil microbial communities in forests with different types of vegetation.


We selected three typical vegetation types, larch (LG), black birch (BD), and larch and black birch mixed (LGBD) forests, located in the Heilongjiang Zhongyangzhan Black-billed Capercaillie Nature Reserve. The Biolog-Eco microplate technology was selected to perform these analyses.


Our results showed clear differences between microorganisms in the three typical forests. The average well colour development (AWCD) change rate gradually increased with incubation time. The BD type had the highest AWCD value, followed by LGBD; the LG forest type had the lowest value. The difference in the soil microbial alpha diversity index between BD and LG was significant. A principal component analysis showed that PC1 and PC2 respectively explained 62.77% and 13.3% of the variance observed. The differences in the soil microbial carbon-source utilisation patterns under different vegetation types were mainly caused by esters and carbohydrates. Redundancy analysis showed that soil microbial functional diversity was strongly affected by soil physicochemistrical properties (e.g. organic carbon, total nitrogen and pH).


These results provide a reference for further exploring the relationship between forest communities and soil microbes during the process of forest succession.


As important constituents of ecosystems and essential active components in the process of biological evolution, microorganisms are widely involved in soil biological and biochemical processes (Gao et al. 2015a, 2015b). Soil microorganisms are an important driving factor for material conversion and nutrient cycling. They play key roles in ecosystem functional maintenance and biogeochemical cycle processes, directly or indirectly promote plant growth, influence vegetation succession processes, and can significantly impact ecosystem development (Fierer and Jackson 2006; Zhang et al. 2012; Bellingham et al. 2001). Geological conditions, vegetation composition and other factors also have a feedback effect on soil microbial activity and diversity (Han et al. 2015; Kennedy et al. 2005). Changes in aboveground vegetation composition influence soil microbial biomass, activity, community structure and function (Hua et al. 2004; Sarah et al. 2010; Yoshitake et al. 2013; Glanville et al. 2012). Changes in soil microbial community and function also impact aboveground plant composition and diversity. Therefore, soil microorganisms are of great significance to the forest ecosystem and have critical roles in vegetation renewal and succession processes. The study of how vegetation impacts soil microbes is an increasingly popular topic for ecologists, especially research regarding the functional mechanisms of ecosystems.

Current research mainly focuses on studying the metabolic diversity of soil microbial communities in different soil types, climatic conditions, human disturbances and different vegetation rhizosphere environments. The activity and the function of soil microbial communities can be estimated by determining the number of bacteria and type of metabolic diversity present. The number of bacteria is a key factor to study microbial community size. Biolog-Eco 96-well micro-plate includes 31 kinds of single carbon source and blank control. The 31 single carbon sources can be roughly divided into 6 categories, which are carbohydrate, amino acids, esters, alcohols, amine and acids. In this wells of micro-plate, there is also a tetrazolium violet dye. Therefore, the free electrons generated by microorganisms using the carbon source can react with tetrazolium violet, and the change in colour reflects the microorganism’s influence on the carbon source. According to the degree of utilisation, different microbial communities have different types of utilisation of 31 single carbon sources, which reflects the diversity of microbial metabolic functions (Garland 1997; Zheng et al. 2004). It can be used to obtain microbial information that reflects the population’s metabolic characteristics and functions by testing the single energy carbon utilisation abilities of the microorganisms. The phospholipid fatty acid (PLFA) method and the various molecular biological advanced methods (e.g. high-throughput sequencing) available do not reflect the soil microbial function (Bossio and Scow 1995; Buyer and Drinkwater 1997; Yang et al. 2002; Xi et al. 2005). Conversely, the Biolog-Eco method can reflect information about the overall activity and metabolic function of microbial communities. The Biolog-Eco method has been utilised to study the impact of land use and human disturbance on soil microorganisms under environmental pollution stress. Although many studies have focused on soil microbes, few addressed the functional diversity of soil microorganisms in different forest types.

The Heilongjiang Zhongyangzhan Black-billed Capercaillie Nature Reserve is located in the transitional zone of the large Xing’an Mountains and small Xing’an Mountains. The forest types are rich, and there are obvious staggered patches. The main forest types are larch (Larix gmelinii), white birch (Betula platyphylla), black birch (Betula dahurica) and Mongolian oak (Quercus mongolica). Nevertheless, how the various types of vegetation affect soil microbial functional diversity remains unclear. In this study, Biolog-Eco microplate technology was used to investigate the variation in soil microbial functional diversity in three typical vegetation types (larch, LG; black birch, BD; and larch and black birch mixed, LGBD) in the Heilongjiang Zhongyangzhan Black-billed Capercaillie Nature Reserve. The aim was to explore the impact of forest vegetation changes and diversity on soil microbial community function and provide a basis to study the relationship between soil microbes and environmental factors in the Heilongjiang Zhongyangzhan Black-billed Grouse National Nature Reserve.


Soil physical-chemical properties

The physical and chemical properties of the three different soil forest types tested are summarised in Table 1. Significant differences in soil pH, organic carbon, total nitrogen and carbon-to-nitrogen ratio (P < 0.05) were observed amongst the three forest types. BD forest soil had the highest pH value (5.86) and LG forest had the lowest (4.61). BD also had the highest value of soil organic carbon (SOC) (60.3−1), whilst LG had the lowest (49.9−1). The same trend was observed for total nitrogen (TN), with BD having the highest (4.12−1) and LG the lowest (2.76−1) values. Finally, LG and BD had the highest (18.1) and lowest (14.7) soil carbon to nitrogen ratio (C/N), respectively.

Table 1 Soil chemical and physical properties in three different forest types

Difference in average colour change rate of soil microorganisms

The AWCD is used to characterise microbial communities. It is an important indicator of the soil microbial community’s ability to use a single carbon source, and it reflects soil microbial activity and diversity of physiological functions (Kenarova et al. 2014; Magdalena et al. 2012). A high AWCD can lead to high activity and density of soil microorganisms; conversely, low AWCD decreases microorganism density and activity. The metabolic activity for soils of the three forest types was analysed from 24 to 240 h. As seen in Fig. 1, the AWCD value for the three forest types increased with increasing the incubation time. LD and LGBD had a significant increasing trend within the first 24 h, but the value strongly increased for BD. This rapid increase in AWCD indicated that the soil microbial community could metabolise organic substrates on the Biolog-Eco microplates. The increasing trend was maintained from 24 to 144 h (Fig. 1). The highest AWCD value occurred at 96 h, which indicated the soil microorganism had the highest ability to utilise carbon sources at that time point. After 96 h, the AWCD increase decelerated and tended to stabilise after 216 h of incubation, indicating all culturable microorganisms can use the carbon source stably during the stable period (Miyake et al. 2016).

Fig. 1

Average well colour development (AWCD) of the soil community for three forest types with incubation time

BD reached the highest AWCD value amongst the three forest-type soil samples. The AWCD increased from 0 to 1.43, at 216 h, for BD, and approximately to 1.14 for the mixed LGBD, indicating that the broad-leaved BD forest had the highest soil microbial activity. However, the AWCD of LG tended to stabilise around 1.12, indicating that the coniferous forest had the lowest soil microbial activity. In summary, different vegetation types caused differences in the microbial AWCD utilisation rate, with BD having the highest value, followed by LGBD, and with LG having the lowest AWCD.

Changes in soil microbial functional diversity

In this study, the Biolog-Eco board had different types of carbon sources that were used to treat the soil microbial communities. Therefore, the alpha diversity could be used as an indicator of the soil microbial communities’ diversity in various perspectives. The absorbance values of the soil microorganisms at 96 h were used to determine Shannon-Wiener (H), Simpson (D) and Mclntosh (U) values. As shown in Table 2, the diversity index reflects the impact of different plant compositions on microbial community functions. Overall, there were significant differences in the Shannon-Wiener, Simpson and McIntosh indexes between the three forest types (Table 2, P < 0.05). The Shannon-Wiener index was the highest for BD, at 3.1. The McIntosh index was the highest for LG, at 9.3, and the lowest for BD, at 7.1.

Table 2 Soil microbial functional diversity of carbon utilisation

Carbon source metabolisation response rates of microbial communities from various vegetation types

Based on physiological and biochemical characteristics, the 31 carbon substrates of the Biolog-Eco microplate were categorised into 6 types: carbohydrates, amino acids, esters, alcohols, amines and acids (Cao et al. 2018) (Table S1). The microbial metabolic fingerprint represents the ability of microorganisms to utilise a wide variety of carbon sources. The utilisation capacity of the 6 carbon sources varied with the microbial communities. As seen in Fig. 2, for BD, there was 1 type of carbon source with an AWCD higher than 2.0 in the metabolic fingerprint, l-asparagine (2.2242). Six kinds of sources were observed in LG: glycogen, d-cellobiose, Tween 40, d-mannitol, phenylethylamine, N-acetyl-d-glucosamine and d-cellobiose, which had a highest AWCD of 2.5686. Finally, there were 5 types in LGBD: l-arginine, l-asparagine, Tween 40, d-mannitol, γ-hydroxybutyric acid and l-asparagine, which had a highest AWCD of 2.7193.

Fig. 2

Metabolic fingerprint of carbon physiological profiles of the microbial community. Note: LG, Larix gmelinii; BD, Betula dahurica; LGBD, mixed Larix gmelinii and Betula dahurica. B2, H1, A2, G2, E1, F1 and G1: carbohydrates; A4, B4, C4, D4, E4 and F4: amino acids; B1, C1, D1 and A3: esters; C2, D2 and H2: alcohols; G4, H4 and E2: amines; B3, F2, C3, D3, E3, F3, G3 and H3: acids

Heatmap analysis visually showed the difference in the metabolic rates of different carbon sources for the microbial communities in all soil samples (Fig. 3), highlighting differences for the various vegetation types. The heatmap results were separated into three categories (Fig. 3). In group I, the AWCD values of D, l-α-glycerol phosphate, glycogen, α-cyclodextrin, l-phenylalanine, d-xylose and d-cellobiose were much higher for the LG forest than for LGBD and BD. In group II, the AWCD values of putrescine, l-serine, Tween 80, γ-hydroxybutyric acid, l-arginine, l-asparagine, d-galacturonic acid, 4-hydroxy benzoic acid and glycyl-l-glutamic acid were also higher for LG than for LGBD and BD. In group III, the AWCD values of l-phenylalanine, l-threonine, α-ketobutyric acid and pyruvic acid methyl ester, and of i-erythritol and itaconic acid were not significantly different between forest types. This indicated that LG had higher density and enhanced soil microorganism activity than the BD forest and LGBD forests.

Fig. 3

Heatmap and hierarchical cluster analysis based on the AWCD of 72 h for the soil microbial communities with different revegetation types. Note: The samples are grouped based on the similarity to each other, and the clustering results are arranged horizontally according to the clustering results. In the figure, red represents the higher AWCD in the corresponding sample, and green represents the lower AWCD. LG, Larix gmelinii; BD, Betula dahurica; LGBD, mixed Larix gmelinii and Betula dahurica

PCA of carbon source metabolisation

For the Biolog-Eco, the load value of the 31 carbon sources in the main component increased with the impact of the carbon source on the main component. As seen in Table S1, 24 carbon sources constituted the first main subdivision (PC1). These included 5 carbohydrates, 6 amino acids, 3 alcohols, 2 amines and 8 acids, of which 2-hydroxy benzoic acid, γ-hydroxybutyric acid and itaconic acid were the most relevant carbon sources to PC1 (1.000). Therefore, acids had a major effect on PC1. In addition, 14 kinds of carbon sources constituted the secondary major subdivision (PC2). These included 3 carbohydrate types, 2 amino acid types, 4 ester types, 2 amines types and 3 acid types, amongst which d-galacturonic acid was the most important carbon source. Therefore, carbohydrates, esters and acids had important effects on PC2.

A PCA analysis was performed to find the variance of carbon utilisation of soil microorganisms in different samples (Felipe-Sotelo et al. 2008; Illian et al. 2009; Sui et al. 2016a, 2016b). The variance contribution rate of principal component 1 was 62.77%, and the contribution rate of principal component 2 was 13.38% (Fig. 4). After a dimensionality reduction, the location images of the points in the main carrier space of the carbon source directly reflected the differences in soil carbon utilisation ability of microorganisms in soils from different forest types (Kim et al. 2017). In this study, PCA of carbon source metabolism of microbial communities indicated that the soil microbial community function in LG coniferous forests tended to be far from that of LGBD and BD broad-leaved forests, especially along the PCA1 axis (Fig. 4), whilst the soil microbial community function of LGBD and BD forests tended to be similar (Fig. 4). This means that soil microbial carbon source metabolic capacity varied significantly between three forest types.

Fig. 4

Principal component analysis of soil carbon source metabolism of the different microbial communities. Note: BD, Betula dahurica; LG, Larix gmelinii; LGBD, mixed Larix gmelinii and Betula dahurica

Differences and factors influencing soil microbial carbon source utilisation patterns in different vegetation types

The environmental factors that affected soil microbial community activity and soil physicochemical properties were analysed through canonical correlation analyses. The two CCA axes explained 56.43% and 12.89% of the variance (Fig. 5). The first CCA axis was close to the SOC (r = 0.98) and TN (r = 0.96). The second CCA axis was negatively correlated with soil pH (r = −0.4820) and C/N (r = −0.6165). Soil microbial community activities significantly correlated with SOC and TN, whilst the BD and LG forests significantly correlated with C/N (Fig. 5).

Fig. 5

Canonical correspondence analysis (CCA) ordination plots of the functional microbial community structure under different revegetation types. Note: SOC, soil organic carbon; TN, total nitrogen; C/N, soil carbon to nitrogen ratio. BD, Betula dahurica; LG, Larix gmelinii; LGBD, mixed Larix gmelinii and Betula dahurica

We analysed the relationship of soil physical and chemical properties and soil microbial alpha diversity (Table 3). As shown in Table 3, for the soil physical and chemical properties, the Shannon index was positively correlated with soil organic carbon (r = 0.678, P < 0.05), and total nitrogen (r = 0.828, P < 0.01) was extremely and positively related with C/N (r = −0.758, P < 0.01).

Table 3 Correlation coefficients between soil physicochemical properties and soil microbial functional diversities


Different vegetation types correlate with different soil physicochemical characteristics

Climate, vegetation and soil physical and chemical properties are important environmental factors related to soil microbial characteristics. In this study, the three forest types investigated (Larix gmelinii coniferous forest, LG; Betula dahurica broad-leaved forest, BD; and Betula dahurica and Larix gmelinii mixed forest, LGBD) had the same soil texture and the same climate. Therefore, differences in soil physicochemical properties and microbial alpha diversity can significantly affect the quality and quantity of organic carbon in the different vegetation types. In this study, the soil pH ranged from 4.61 to 5.86. The soil pH was the lowest for LG and the highest for BD. The differences may come from the chemical composition of different litters (Waid 1999). Litter quality and total nitrogen content were lower in LG than in other forest types, whereas the C/N, lignin content and lignin/N were higher in LG. Therefore, the decomposition rate of the LG forest litter and the release rate of plant nutrients should be reduced (Gao et al. 2016). This may be the cause for the low soil pH in the LG forest.

In addition, we found that the total carbon and total nitrogen content of the BD ecosystem were higher than in LG and LGBD forests (Table 1). Our results are consistent with previously published findings that broad-leaved forests have higher soil organic matter and nitrogen content than coniferous forests (Ding et al. 2017). Other studies have also found that the values for these soil physicochemistrical properties were lower in coniferous forests than in broad-leaved forests (Gao et al. 2015a, 2015b; Geng et al. 2013). The high content of soil organic carbon and total nitrogen in broad-leaved forests may be due to the continuous input of litter, root effects and turnover of soil nutrients (Ge et al. 2012). Litter in the larch forests contains high levels of difficult-to-decompose lignin and organic acids (Ge et al. 2012). Research also demonstrated that the soil carbon-to-nitrogen ratio in various forest types can drastically change the structure and function of soil microbes (Jiang et al. 2014). In this study, the soil C/N ratio was higher in LG than in the BD forest, which contradicts previously published results (Yoshimura et al. 2008; Stéphane et al. 2008). This may be because the LG forest in this study has a single community type and a single type of litter. The chemical composition of the main litter has a significant impact on soil chemical properties (Fu et al. 2019). A large number of previous studies have shown that the composition ratio of aboveground tree species significantly impacts the soil, which is consistent with the results of this study (Chodak and Niklinska 2010; Huang et al. 2012; Wu et al. 2013).

Different vegetation types contribute to different soil microbial functional diversities

In this study, the soil alpha diversity was significantly different in the three forest types (Table 2). These results indicate that the type of forest vegetation correlates with soil microbial functional diversity. The metabolic functional diversity of the soil microbial community is reflected by the functional diversity index (Sui et al. 2016a, 2016b). Our research demonstrated that the soil microbial functional diversity index was higher in BD than in the LG and LGBD forests. Previous research has indicated that the Shannon-Wiener index, Simpson index and McIntosh index of soil microorganisms are higher in broad-leave than in coniferous forests (Hu et al. 2010). Other studies have shown this using different method (Wu et al. 2013). In this study, the soil microorganisms’ functional diversity was the lowest in coniferous LG forests and the highest for broad-leaf BD forests. The soil microbial Biolog-Eco functional diversity Shannon index was lower in LG than in LGBD and BD. These analyses found that the soil microbial functional diversity is correlated to the forest’s dominant tree species. Indeed, various forest types impact microbial functional diversity by influencing soil physical and chemical properties (Anderson et al. 2009).

Our research found that microbial community functional activities were significantly different in the three forest types (Figs. 2 and 3), which may be closely related to the soil physical and chemical characteristics (Fig. 5, Table 3). As evident in the heatmap and PCA results, the beta diversity of soil microbial communities was significantly different in the three forests. It is known that regions with different vegetation types also differ in the soil microorganism composition (Ushio et al. 2010a, 2010b; Eilers et al. 2010). The organic carbon and total nitrogen contents were positively correlated with the Shannon index in our data, and soil C/N was negatively correlated with the same index. These results are inconsistent with what was reported by Deng et al. (2019). Different vegetation types produce different amounts and qualities of plant litter input (Ushio et al. 2010a, 2010b; Zhang et al. 2006), promoting different types of microorganisms in the soil (Sariyildiz and Anderson 2003; Tripathi et al. 2013). Therefore, these results indicate that the difference in litter composition of each vegetation composition leads to changes in nutrient content and pH in the soil, affecting the structure and function of microorganisms present.


In this paper, we analysed the changes in metabolic activity and functional diversity of soil microbial communities in three different forest types. The results showed that different vegetation types led to significant differences in the functional diversity of soil microbial communities. The metabolic activities were significantly higher in the coniferous and broad-leaved mixed and broad-leaved forests than in the coniferous forest type. Additionally, the broad-leaved forests had the highest soil microbial functional diversity. Soil organic carbon and total nitrogen are the main parameters that affect the functional diversity of soil microbial communities. Therefore, forest composition can change the functional diversity of microorganisms by altering the physical and chemical properties of the soil.

Due to the limitation of carbon source types available, the Biolog-Eco method is not sufficient to completely explore the function and diversity of microorganisms. In the future, it will be necessary to combine enzyme kinetics and molecular biology approaches, especially high-throughput sequencing, in this type of research.


Study area

The experimental site is located in the Heilongjiang Zhongyangzhan Black-billed Capercaillie National Nature Reserve (Fig. 6). The reserve is located in south-western Heihe City, Heilongjiang Province, and this reserve is the transition zone between the south-western mountains of Xiaoxing’anling and the Songnen Plain. The geographical location is 126° 00′–126° 45′ E, 48° 30′–48° 50′ N, and the total area is 988.6 km2. This area is located at a temperate continental monsoon climate zone, with severe cold in winter and cool and short summers. The annual average temperature is −0.5 °C, the frost-free period is 121 days, the annual average rainfall is 476.33 mm, and the average relative humidity for many years has been 69.2%. The different vegetation types are randomly distributed and are affected by the surrounding environment. The plant communities are in different succession stages. The research site (nature reserve) is a typical alpine forest ecosystem with forests, shrubs, wetlands and meadows. Forests include coniferous and broad-leaved mixed forests, coniferous forests and broad-leaved forests. Wetlands include forest swamps, shrub swamps and herb marshes and swampy meadows and other wetland types. Most of the nature reserve is covered by forests, with a forest coverage rate of 82.4% (Yin et al. 2019).

Fig. 6

Location of the study area

The main plants present are Larix gmelinii, Tilia amurensis, Chosenia arbutifolia, Betula platyphylla, Betula dahurica, Quercus mongolica, Populus davidiana, Prunus padus, Salix raddeana and Alnus mandshurica.

Plot setting and sample collection

In July 2019, plot selection and sample collection were conducted in a day with sunshine and a temperature of 20 °C. Three typical vegetation types, Larix gmelinii (coniferous larch forest, LG), Betula dahurica (broad-leaved black birch forest, BD) and Betula dahurica and Larix gmelinii mixed forest (LGBD) were selected in larger Xing’an. Three plots were used for each vegetation type. The area of each plot was approximately 400 m2. In each forest type, three plot (10 m × 10 m) were established and the distance between any two plots was > 50 m. Soil samples were collected using a sterile soil drill from 5 randomly selected locations inside each plot and were pooled to obtain a mixed soil specimen (approximately 1 kg of fresh soil) for each plot. After removing plant debris, the samples were passed through a 2-mm sieve and placed in sterile bags. The samples used for the determination of soil microbial diversity and other indicators were placed in an incubator with an ice bag and kept in the laboratory at 4 °C. Another part of the sample was dried naturally, and soil nutrients and other indicators were measured after sieving.

Determination of soil physical and chemical factors

Soil pH was analysed by the water extraction potential method (water-soil ratio 2.5:1). Soil organic carbon (soil organic carbon, SOC) and total nitrogen (total nitrogen, TN) were determined with an elemental analyser (Arvato EA3000).

Determination of the functional diversity of soil microbial community

The Biolog-Eco micro-plate was used to test the soil microbial carbon metabolic function for three different forest types. The first well works as a control without a carbon source, and the others contain a carbon source and redox dye tetrazolium blue. Microorganisms use carbon sources to respire changing the redox potential and reducing tetrazoles (TV) from colourless to purple. Ten grammes of fresh soil was filtered through a 2-mm sieve into an autoclaved Erlenmeyer flask, and 100 ml of 0.85% sterile sodium chloride solution was added. The sample was sealed, shaken at 180 r/min for 30 min at 25 °C and statically incubated for 2 min. The suspensions were diluted (10−3) in NaCl solution and inoculated onto Biolog-Eco plates (150 μL per well) using a multichannel pipette. All treatments had three triplicates. The absorbance was measured at 590 nm (0.13 ± 0.02). The plate was incubated at 25 °C, and the absorbance of the ecological plate was measured once every 12 h, continuously for 10 d (Classen et al. 2003; Liao et al. 2013; Jin et al. 2014).

Statistical analysis

AWCD was used to describe the soil microbial metabolic activity. The calculation formula is as follows:

$$ AWCD=\sum \left({C}_i-R\right)/31 $$

In the formula, Ci is the absorbance value of the well with a carbon source, R is the absorbance value of the control well (the control is empty), and n is the number of carbon sources (31 × 3 repetitions) (Velasco et al. 2009). If Ci − R ≤ 0, then the whole value is recorded as 0 (Kong et al. 2008).

The 96-h cultivation data were used to calculate the functional diversity of the soil microbial communities. The calculation formula is as follows:

Shannon-Wiener diversity index (Kong et al. 2008):

$$ H={P}_i\ \ln\ {P}_i $$

Simpson diversity index (Kong et al. 2008):

$$ D=1-\sum {\left({P}_i\right)}^2 $$

McIntosh diversity index (Xie et al. 2011):

$$ U=\sqrt{\sum {n_i}^2} $$

In the formula, Pi is the ratio of the relative absorbance value of the ith to the sum of the relative absorbance values of all.

Data analysis was performed with Excel, SPSS 20.0, and R. SPSS was used for one-way ANOVA analyses with the detection level set at 0.05. Heatmap, principal component analysis (PCA) and canonical correspondence analysis (CCA) were performed with R (Vegan package).



Betula dahurica


Larix gmelinii


Larix gmelinii and Betula dahurica mixed


Average well colour development


Canonical correspondence analysis


Soil organic carbon


Total nitrogen


Soil carbon to nitrogen ratio


  1. Anderson JAH, Hooper MJ, Zak JC, Cox SB (2009) Characterization of the structural and functional diversity of indigenous soil microbial communities in smelter-impacted and nonimpacted soils. Environ Toxicol Chem 28(3):534–540.

    CAS  Article  PubMed  Google Scholar 

  2. Bellingham PJ, Walker LR, Wardle DA (2001) Differential facilitation by a nitrogen-fixing shrub during primary succession influences relative performance of canopy tree species. J Ecol 89(5):861–875.

    Article  Google Scholar 

  3. Bossio DA, Scow KM (1995) Impact of carbon and flooding on the metabolic diversity of microbial communities in soils. Appl Environ Microbiol 61(11):4043–4050.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. Buyer JS, Drinkwater LE (1997) Comparison of substrate utilization assay and fatty acid analysis of soil microbial communities. J Microbiol Methods 30(1):1–11.

    Article  Google Scholar 

  5. Cao HJ, Wang LM, Xu MY, Huang QY, Xie LH, Luo CY, Ni HW (2018) Functional diversity of soil microbial communities of different vegetation types in the new volcanic lava platform of Wudalianchi. Acta Eco Sini 39(21):7927–7937.

    Article  Google Scholar 

  6. Chodak M, Niklinska M (2010) Effect of texture and tree species on microbial properties of mine soils. Appl Soil Ecol 46(2):268–275.

    Article  Google Scholar 

  7. Classen AT, Boyle SI, Haskins KE, Overby ST, Hart SC (2003) Community-level physiological profiles of bacteria and fungi: plate type and incubation temperature influences on contrasting soils. FEMS Microbiol Ecol 44(3):319–328.

    CAS  Article  PubMed  Google Scholar 

  8. Deng JJ, Zho YB, Bai X, Luo JY, Zhu WX (2019) Soil microbial functional diversity responses to different revegetation types in Baishilazi Nature Reserve. Pol J Environ Stud 28(5):3675–3686.

    Article  Google Scholar 

  9. Ding GC, Wang XH, Yang QF, Lin QX, Huang ZQ (2017) Effects of tree species transition on soil microbial community composition and functions in subtropical China. J Appl Ecol 28(11):3751–3758.

    Article  Google Scholar 

  10. Eilers KG, Lauber C, Fierer N (2010) Shifts in bacterial community structure associated with inputs of low molecular weight carbon compounds to soil. Soil Biol Biochem 42(6):896-903.

  11. Felipe-Sotelo M, Tauler R, Vives I, Grimalt JO (2008) Assessment of the environmental and physiological processes determining the accumulation of organochlorine compounds in European mountain lake fish through multivariate analysis. Sci Total Environ 404(1):148–161.

    CAS  Article  PubMed  Google Scholar 

  12. Fierer N, Jackson RB (2006) The diversity and biogeography of soil bacterial communities. Proc Natl Acad Sci U S A 103(3):626–631.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. Fu Q, Xing YJ, Yan GY, Dong XD, Zhang JH, Wang QG (2019) Response of litter dynamics of boreal forest to long-term nitrogen deposition. Ecol Environ Sci 28(7):1341–1350

    Google Scholar 

  14. Gao C, Zhang Y, Shi NN, Zheng Y, Chen L, Wubet T (2015b) Community assembly of ectomycorrhizal fungi along a subtropical secondary forest succession. New Phytol 205(2):771–785.

    Article  PubMed  Google Scholar 

  15. Gao F, Jiang H, Cui XY (2015a) Soil organic carbon pools and their turn over under two different types of forest in Xiaoxing’an Mountains, Northeast China. J Appl Ecol 26(7):1913–1920

    CAS  Google Scholar 

  16. Gao J, Kang F, Han H (2016) Effect of litter quality on leaf-litter decomposition in the context of home-field advantage and non-additive effects in temperate forests in China. Pol J Environ Stud 25(5):1911–1920.

    Article  Google Scholar 

  17. Garland JL (1997) Analysis and interpretation of community-level physiological profiles in microbial ecology. FEMS Microb Ecol 24(4):289–300.

    CAS  Article  Google Scholar 

  18. Ge XG, Huang ZL, Cheng RM, Zeng LX, Tan BW (2012) Effects of litter fall and root input on soil physical and chemical properties in Pinus Massoniana plantations in Three Gorges Reservoir Area, China. J Appl Ecol 23(12):3301–3308

    CAS  Google Scholar 

  19. Geng SC, Chen ZJ, Zhang JH, Lou X, Yu DD (2013) Soil methane fluxes of three forest types in Changbai Mountain of northeast China. Chin J Ecol 32(5):1091–1096

    Google Scholar 

  20. Glanville HC, Hill PW, Maccarone LD, Golyshin P, Murphy DV, Jones DL (2012) Temperature and water controls on vegetation emergence, microbial dynamics, and soil carbon and nitrogen fluxes in a high arctic tundra ecosystem. Funct Ecol 26(6):1366–1380.

    Article  Google Scholar 

  21. Han DX, Wang N, Wang NN, Sun X, Feng FJ (2015) Soil microbial functional diversity of different altitude Pinus Koraiensis forests. Chin J Appl Ecol 26(12):3649–3656

    CAS  Google Scholar 

  22. Hu C, Fu B, Liu G, Jin T, Guo L (2010) Vegetation patterns influence on soil microbial biomass and functional diversity in a hilly area of the loess plateau, China. J Soils Sediments 10(6):1082–1091.

    CAS  Article  Google Scholar 

  23. Hua Z, Zhiyun O, Xiaoke W, Zhiguo F, Tongqian Z, Hong M (2004) Effects of forest restoration patterns on soil microbial communities. Chin J Appl Ecol 15(11):2019–2024.

    Article  Google Scholar 

  24. Huang YY, Qu LY, Qu XC, Du SQ (2012) Soil microbial community characteristics under different vegetation types at the Holocene-basalt Platform, Jingpo Lake area, northeast China. Acta Ecol Sin 32(9):2827–2836.

    Article  Google Scholar 

  25. Illian JB, Prosser JI, Baker KL, Rangel-Castro JI (2009) Functional principal component data analysis: a new method for analysing microbial community fingerprints. J Microbiol Methods 79(1):89–95.

    CAS  Article  PubMed  Google Scholar 

  26. Jiang YS, Xiao J, Luo CL, Li J, Zhang G (2014) The microbial community of Three Mountains Pingding, Changbai and Saihanwula based on the analysis of phospholipid fatty acids (PLFAs). Ecolo Environ Sci 23(8):1272–1279

    Google Scholar 

  27. Jin Z, Ji FY, Xu X, Xu XY, Chen QK, Li Q (2014) Microbial and metabolic characterization of a denitrifying phosphorus-uptake/side stream phosphorus removal system for treating domestic sewage. Biodegradation 25(6):777–786.

    CAS  Article  PubMed  Google Scholar 

  28. Kenarova A, Radeva G, Traykov I, Boteva S (2014) Community level physiological profiles of bacterial communities inhabiting uranium mining impacted sites. Ecotox Environl Safe 100:226–232.

    CAS  Article  Google Scholar 

  29. Kennedy N, Brodie E, Connolly J, Clipson N (2005) Impact of lime, nitrogen and plant species on fungal community structure in grassland microcosms. Environ Microbiol 7(6):780–788.

    CAS  Article  PubMed  Google Scholar 

  30. Kim JW, Rehmann L, Ray MB (2017) Development of microalgal bioassay based on the community level physiological profiling (clpp). Algal Res 25:47–53.

    Article  Google Scholar 

  31. Kong WD, Zhu YG, Fu BJ, Han XZ, Zhang L, He JZ (2008) Effect of long-term application of chemical fertilizers on microbial biomass and functional diversity of a black soil. Pedosphere 18(6):801–808.

    CAS  Article  Google Scholar 

  32. Liao M, Xie XM, Peng Y, Chai JJ, Na C (2013) Characteristics of soil microbial community functional and structure diversity with coverage of Solidago Canadensis L. J Cent South Univ 20(3):749–756.

    Article  Google Scholar 

  33. Magdalena F, Oszust K, Lipiec J (2012) Community level physiological profiles (CLPP), characterization and microbial activity of soil amended with dairy sewage sludge. Sensors 12(3):3253–3268.

    Article  Google Scholar 

  34. Miyake H, Maeda Y, Ishikawa T, Tanaka A (2016) Calorimetric studies of the growth of anaerobic microbes. J Biosci Bioeng 122(3):364–369.

    CAS  Article  PubMed  Google Scholar 

  35. Sarah R, Sattin CC, Cleveland (2010) Functional shifts in unvegetated, perhumid, recently-deglaciated soils do not correlate with shifts in soil bacterial community composition. J Microbiol 47(6):673–681.

    Article  Google Scholar 

  36. Sariyildiz T, Anderson JM (2003) Interactions between litter quality, decomposition and soil fertility: a laboratory study. Soil Biol Biochem 35(3):391–399.

    CAS  Article  Google Scholar 

  37. Stéphane B, Frédéric MG, Balandier P (2008) Influence of tree species on understory vegetation diversity and mechanisms involved—a critical review for temperate and boreal forests. Forest Ecol Manag 254(1):1–15.

    Article  Google Scholar 

  38. Sui X, Zhang RT, Yang LB, Xu N, Chai CR, Wang JF (2016a) Effect of simulation nitrogen depositions on bacterial diversity of Deyeuxia angustifolia in wetland of Sanjiang Plain. Pratacult Sci 33(4):589–598

    CAS  Google Scholar 

  39. Sui X, Zhang RT, Yang LB, Xu N, Zhong HX, Wang JF, Ni HW (2016b) Functional diversity of soil bacterial communities in different types of Calamagrostis angustifolia wetlands in Sanjiang Plain. Res Environ Sci 29(10):1479–1486

    CAS  Google Scholar 

  40. Tripathi G, Deora R, Singh G (2013) The influence of litter quality and micro-habitat on litter decomposition and soil properties in a silvopasture system. Acta Oecol 50:40–50.

    Article  Google Scholar 

  41. Ushio M, Kitayama K, Balser TC (2010a) Tree species effects on soil enzyme activities through effects on soil physicochemical and microbial properties in a tropical montane forest on Mt. Kinabalu, Borneo. Pedobiologia 53(4):230–233.

    CAS  Article  Google Scholar 

  42. Ushio M, Kitayama K, Balser TC (2010b) Tree species-mediated spatial patchiness of the composition of microbial community and physicochemical properties in the topsoils of a tropical montane forest. Soil Biol Biochem 42(9):1588–1595.

    CAS  Article  Google Scholar 

  43. Velasco GV, Probanza A, Ero FJG, Trevi OAC, Moreno JM, Garcia JAL (2009) Effect of fire and retardant on soil microbial activity and functional diversity in a Mediterranean pasture. Geoderma 153(1-2):185–193.

    CAS  Article  Google Scholar 

  44. Waid JS (1999) Does soil biodiversity depend upon metabiotic activity and influences? Appl Soil Ecol 13(2):151–158.

    Article  Google Scholar 

  45. Wu ZY, Lin WX, Chen ZF, Fang CX, Zhang ZX, Wu LK, Zhou MM, Shen LH (2013) Characteristics of soil microbial communities of different vegetation types in Wuyishan National Nature Reserve. Chin J Appl Ecol 24(8):2301–2309

    CAS  Google Scholar 

  46. Xi JY, Hu HY, Jiang J, Qian Y (2005) Metabolic characteristics of the microbial community in the biological filter tower. Environ Sci 26(4):167–172

    Google Scholar 

  47. Xie XM, Liao M, Ma AL, Zhang HJ (2011) Effects of contamination of single and combined cadmium and mercury on the soil microbial community structural diversity and functional diversity. Chin J Geochem 30(3):366–374.

    CAS  Article  Google Scholar 

  48. Yang YG, Paterson E, Campbell C (2002) Application of Biolog method to distinguish microbial characteristics between urban soil and rural soil. Acta Pedol Sin 39(4):582–589

    CAS  Google Scholar 

  49. Yin WP, Li GF, Zhao ZZ, Chen FY, Qiu GL, Liu WS (2019) Species diversity and distribution of orchids in Heilongjiang Zhongyangzhan Black-billed Capercaillie National Nature Reserve. J Qiqihar Uni 35(5):52–65

    Google Scholar 

  50. Yoshimura C, Gessner MO, Furumai TH (2008) Chemical properties, microbial respiration, and decomposition of coarse and fine particulate organic matter. J N Am Benthol Soc 27(3):664–673.

    Article  Google Scholar 

  51. Yoshitake S, Fujiyoshi M, Watanabe K, Takehiro M (2013) Successional changes in the soil microbial community along a vegetation development sequence in a subalpine volcanic desert on Mount Fuji, Japan. Plant Soil 364(1-2):261–272.

    CAS  Article  Google Scholar 

  52. Zhang QR, Zhou QX, Ren LP, Zhu YG, Sun SL (2006) Ecological effects of crude oil residues on the functional diversity of soil microorganisms in three weed rhizospheres. J Environ Sci 18(6):1101–1106.

    Article  Google Scholar 

  53. Zhang X, Xuan Z, Zhao M (2012) Functional diversity changes of microbial communities along a soil aquifer for reclaimed water recharge. FEMS Microb Ecol 80(1):9–18.

    CAS  Article  Google Scholar 

  54. Zheng H, Ouyang ZY, Fang ZG, Zhao TQ (2004) Application of Biolog in functional diversity of soil microbial community. Acta Pedol Sin 1(3):456–461

    Google Scholar 

Download references


We are grateful to Mr. Zhenzhu Zhao, leader of Heilongjiang Zhongyangzhan Black-billed Capercaillie Nature Reserve, for allowing us access to the Nature Reserve; Mr. Zhenzhu Zhao, Fuyuan Chen and Xianda Li for help with sampling, and Mr. Xianda Li for help with sample processing.


The work was funded by National Natural Sciences Foundation of China (32011530425, 32001126), Natural Sciences Foundation of Heilongjiang Province (LH2020C088), Outstanding Youth Foundation of Heilongjiang University (JCL202006) and supported by the State Forestry and Grassland Administration Programme (AS2017) and Heilongjiang Finance Fundamental Financial Special Project (ZNJZ2020R01) .

Author information




Sui Xin and Weng Xiaohong performed this experiment, analysed the data and wrote this MS, and Li Jiayi performed this experiment and help to write this MS; Yin Weiping and Li Guofu helped to collect the soil samples; Ma Weichao helped to analysed the data; Yang Libin and Mu Liqiang designed this experiment and revised this MS. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Xin Sui, Libin Yang or Liqiang Mu.

Ethics declarations

Ethics approval and consent to participate

The study did not violate ethics, and all participants agreed to publish the paper.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplement Table 1.

Loading factors of the principal components of the 31 individual carbon sources.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Weng, X., Li, J., Sui, X. et al. Soil microbial functional diversity responses to different vegetation types in the Heilongjiang Zhongyangzhan Black-billed Capercaillie Nature Reserve. Ann Microbiol 71, 26 (2021).

Download citation


  • Soil microorganisms
  • Vegetation types
  • Carbon source utilisation patterns
  • Biolog microplate technology