Next Article in Journal
The Future of European Beech in Northern Germany—Climate Change Vulnerability and Adaptation Potential
Previous Article in Journal
Effects of Long-Term Subculture on Maturation Ability and Plant Conversion in Pinus radiata: Using FT-IR Spectroscopy to Determine Biomarkers of Embryogenic Tissue Aging
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 14
Issue 7
10.3390/f14071447
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extracellular Enzyme Activity and Nutrient Characteristics of Pinus massoniana Lamb. Families with Different Growth Levels: Insights into the Ectomycorrhizal Fungal Community and Rhizosphere Soil

by
Jun Zhang
1,2,3,
Yuanxiang Zhao
1,2,3,
Hongyang He
1,2,3,
Haoyun Wang
1,2,3 and
Feng Wu
1,2,3,*
1
Institute for Forest Resources and Environment of Guizhou, Guizhou University, Guiyang 550025, China
2
Key Laboratory of Forest Cultivation in Plateau Mountain of Guizhou Province, Guizhou University, Guiyang 550025, China
3
College of Forestry, Guizhou University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Forests 2023, 14(7), 1447; https://doi.org/10.3390/f14071447
Submission received: 21 May 2023 / Revised: 10 July 2023 / Accepted: 11 July 2023 / Published: 14 July 2023
(This article belongs to the Section Forest Ecology and Management)
Browse Figures
Versions Notes

Abstract

:
The symbiosis between ectomycorrhizal (ECM) fungal species and the root system of Pinus massoniana is important for its growth. However, the characteristics of ECM fungal communities and their relationships with extracellular enzyme activities and nutrients in the rhizosphere soil of different P. massoniana genotypes have not been well characterized. In this study, P. massoniana families (groups of offspring from different parents) with different levels of growth were selected for investigating ECM fungal communities, soil nutrients, extracellular enzyme activity, and leaf nutrient concentrations to explore the relationships between P. massoniana and the composition of the ECM fungal community. The high-growth (HG) family of P. massoniana had more different ECM fungal communities than the medium-growth (MG) and low-growth (LG) families; each family had a unique and dominant genera (HG: Amphinema and Pseudoclathrosphaerina; MG: Russula and Auricularia; and LG: Russula and Amanita). Amphinema was the main contributor to the differences among the three families (contribution: HG-MG 0.225 and HG-LG 0.17) and had rich extramatrical mycelium, which favored the growth of the HG family and positively affected the accumulation of soil organic carbon. Structural equation modelling showed that the dominant genera in the HG family had significant positive effects on the activity of three extracellular enzymes (BG, NAG, and AP) (weak to moderate positive effects of Amphinema on BG, NAG, and AP and moderate positive effects of Pseudoclathrosphaerina on BG, NAG, and AP), which might have contributed to the differences in extracellular enzyme activities among the families with different growth levels. Redundancy analysis indicated that P. massoniana growth traits (tree height, diameter at breast height, and timber volume), soil total nitrogen, and the N/P ratio significantly influenced ECM fungal communities. The study revealed the characteristics of ECM fungal communities, soil extracellular enzyme activity, and nutrient features of P. massoniana with different growth levels, which help improve our understanding of the relationship between P. massoniana genotype and ECM fungal communities.

1. Introduction

Soil microbes play a critical role in driving terrestrial ecosystem processes by connecting aboveground biomass with belowground ecosystems [1]. As an important source of soil microbiota, extraradical mycelium (ECM) fungal communities are essential for plant productivity and soil nutrient cycling in forest ecosystems. ECM fungi can enhance the host plant’s nutrient acquisition, pathogen resistance, and drought tolerance and in return obtain carbon from the host [2] while producing various extracellular oxidases and hydrolytic enzymes that mobilize and release smaller organic molecules from soil organic matter (SOM) [3]. As previously mentioned, ECM fungi can enhance nutrient cycling in rhizospheric environments and their species richness and composition are also associated with soil physicochemical properties and nutrient availability [4]. Meanwhile, soil physicochemical properties can drive the community structure of fungi, thereby improving the nutrient environment of the soil [5]. For example, increased nitrogen availability in soils decreases ECM fungal colonization and alters their species richness [6]. Host plants also influence soil fungi, with recent studies showing that plants can directly affect their own root colonization through intra- and inter-specific variation [7,8,9]. While it is true that some mycorrhizal fungi have broad host ranges [10], it is important to note that within a single tree species, the composition of mycorrhizal fungal communities can be highly dynamic during stand development. Additionally, the diversity and composition of mycorrhizal fungi in the subterranean parts of the same species can vary across different developmental stages and ages. These variations are often associated with changes in soil enzyme activity [11]. Furthermore, stands of the same tree species with different growth rates exhibit differences in the range of ECM symbionts in their roots, as fast-growing stands allocate more energy and photosynthetic carbon to fine roots than slow-growing stands [12,13,14].
Soil extracellular enzymes, as crucial functional products of fungi, can be modulated by soil microorganisms to regulate the balance of necessary nutrients, thereby enabling adaptation to climate change [15] and affecting soil carbon sequestration and emission [16]. As a key factor in nutrient mobilization and acquisition [17], extracellular enzyme activity is a useful indicator for evaluating microbial nutrient demand [18,19]. As a proximate factor in the decomposition of soil organic matter and litter, soil extracellular enzymes exert effective control over the cycling of carbon, nitrogen, phosphorus, and other nutrients in terrestrial soil ecosystems [20,21,22]. Soil extracellular enzyme activity is influenced by several factors, including stand age [23], fungal community evolution time [24], and changes in both altitudinal and altitudinal gradients [25]. These variations are caused by the presence of ECM fungi, which are the key producers of extracellular enzymes [26,27], within different soil microenvironments [28,29]. Therefore, it is still unclear as to how soil extracellular enzyme activity is influenced by ECM fungal communities differ under different growth-rate genotypes of trees.
Pinus massoniana Lamb. [30], widely planted in southern China due to its adaptability, rapid growth, and high utilization rate, is used for timber and resin production and accounts for 20% of the national afforestation area [31]. Different P. massoniana families (groups of offspring from different parents) show different mechanisms of response to environmental changes [32,33,34]. This study examines the ectomycorrhizal fungal community and soil extracellular enzyme activity characteristics in the rhizosphere of P. massoniana trees with high, medium, and low growth levels selected according to progeny test data for 9-year-old stands. The aim of this study was to determine whether there are differences in the ECM fungal community and soil extracellular enzyme activity among families with the three different growth levels and to reveal the relationship between the ECM fungal community composition and the growth performance and soil nutrients of P. massoniana. The results can provide a theoretical basis for studying the interaction between P. massoniana aboveground growth and the belowground environment.

2. Materials and Methods

2.1. Study Region

The study area is located in Mochong Town, Duyun City, Guizhou Province (107°26′46″ E, 26°4′34″ N). The minimum temperature of the seed orchard is −6.9 °C, with an average annual temperature of 15.8 °C and an average annual precipitation of over 1400 mm, with the precipitation concentrated from May to October accounting for more than 80% of the annual precipitation. The average relative humidity is approximately 80%: the area is a subtropical humid area. According to the Chinese Soil Classification, the soil type in the study area is yellow soil, mainly composed of loamy sand, and with a pH of 4.77.

2.2. Mycorrhiza, Rhizosphere Soil, and Leaf Sampling

The experimental materials were obtained from the progeny test stand of P. massoniana established in the spring of 2013. The progeny test stand, which was planted using a completely randomized block design, included 46 families with 5 individuals each planted in a single row within each plot. Based on years of growth data, three families with different growth rates were selected: a high growth (HG), medium growth (MG), and low growth (LG), as shown in Table 1. Three blocks were randomly selected and then one individual from each of the three families was chosen randomly for tree samples. In August 2021, the leaves located in the southern part of the crown of each sample tree were collected and stored at −20 °C for subsequent determination of leaf carbon (leaf C), leaf nitrogen (leaf N), and leaf phosphorus (leaf P). Fine roots were excavated in the main root extension direction at the base of the tree trunk and rhizosphere soil was collected at the site of root sampling. The mixed roots and soil samples from each tree (with at least two samples per tree) were packaged and stored in an ice box to prevent loss of activity; please refer to Dong et al. [11] for specific sampling methods. After transportation to the laboratory, the soil samples were divided into two groups: one group was kept at −20 °C for soil extracellular enzyme activity analysis and identification of rhizosphere microbial community composition within one week [33], while the other group was air-dried at room temperature for determination of the soil organic carbon (SOC), total nitrogen (TN), and total phosphorus (TP).

2.3. Soil and Leaf Nutrient Determination

After being air-dried at room temperature in a shaded area for one week and sieved, the pine leaves were dried and ground with a ball mill and 0.1 g samples (accurate to 0.0001 g) were taken for analysis. The total nitrogen of rhizosphere soil and pine leaves were determined by the sulfuric acid digestion and sodium salicylate method; the total phosphorus was determined by the sulfuric acid digestion and molybdate-antimony anti-method. Total nitrogen and total phosphorus analyses were performed using a fully automated chemical sequential analyzer (De Chem-Tech. GmbH, Hamburg, Germany). The organic carbon of rhizosphere soil and total carbon in leaves were determined by potassium dichromate oxidation; the soil pH was measured by potentiometry. The specific methods were consistent with those described by Bao [35].

2.4. Determination of Soil Extracellular Enzyme Activity

β-1,4-Glucosidase (BG, EC 3.2.1.21), β-1,4-N-acetylglucosaminidase (NAG, EC 3.2.1.14), and acid phosphatase (AP, EC 3.1.3.1) activities were determined in this study. These three enzymes are resistant to denaturation and have relatively long half-lives [36]. After soil sample processing, the enzyme activity was measured based on the fluorescence method using a standard set reported by German [37] with a mesh size of 30–50.

2.5. Molecular Identification of the Ectomycorrhizal Community

2.5.1. DNA Extraction, PCR Amplification, and Illumina Sequencing

DNA extraction, PCR amplification, and Illumina sequencing of mycorrhizal fungi were completed by Novogene Co., Ltd. (Beijing, China). The genomic DNA of the samples was extracted by the CTAB method [38] and then the extracted DNA purity and concentration were detected by agarose gel electrophoresis. Appropriate amounts of DNA were then placed into a centrifuge tube and diluted with sterile water to 1 ng/μL. The ITS1-1F region of fungi was amplified by PCR using the diluted genomic DNA as a template. The primers were ITS5-1737F (F: CTTGGTCATTTAGAGGAAGTAA and R: GCTGCGTTCTTCATCGATGC) [39]. The PCR products were examined by agarose gel electrophoresis at a 2% concentration, mixed thoroughly in equal amounts according to the concentration of the PCR products, and then examined again using agarose gel electrophoresis at a 2% concentration after which the target bands were recovered. Libraries were constructed using the NEBNext® Ultra™ IIDNA Library Prep Kit and the libraries were quantified by a Qubit instrument and Q-PCR; after the libraries were qualified, they were sequenced using NovaSeq6000.

2.5.2. Sequencing Data Processing

The samples were split from the downstream data according to the PCR amplification primer sequences and the reads were spliced using FLASH (V1.2.11, http://ccb.jhu.edu/software/FLASH/, accessed on 11 January 2023) software to obtain raw tags [40]. Finally, the clean tags were compared with the database using Vsearch software (version 2.8.1) to detect chimaera and removed to obtain the final effective tags [41].

2.5.3. Amplicon Sequence Variation (ASV) Noise Reduction and Species Annotation

For the effective tags obtained above, the DADA2 module in QIIME2 software (version 2) was used to reduce noise and filter out sequences with abundances less than five to obtain the final amplicon sequence variants (ASVs) and the feature list [42]. The obtained ASVs were subjected to species annotation by aligning them with a pretrained naive Bayes classifier using the classify-sklearn module in the QIIME2 software [43], thereby obtaining species-level information for each ASV based on the UNITEv8.2 database.

2.6. Data Analyses

Differences in rhizospheric soil nutrients, leaf nutrients, and extracellular enzyme activity among the three families of P. massoniana were analysed using one-way analysis of variance (ANOVA) in SPSS software (version 22.0, SPSS Inc. 2018) followed by Duncan’s multiple range test at a significance level of 0.05. Pearson correlation analysis (p < 0.05, p < 0.01, and p < 0.001) was performed using the corrplot package in R software (version 4.0) to examine the correlations among the rhizospheric soil nutrients, leaf nutrients, extracellular enzyme activity, and growth performance of P. massoniana. Simper analysis (similarity percentage) [44] was conducted using the vegan package in R software (version 4.0) to quantify the contribution of ECM fungal species to the differences among the three families of P. massoniana. A structural equation model was constructed using AMOS software (version 4.0) to examine the relationship between the dominant ECM genera of P. massoniana and extracellular enzyme activity. Bar charts were drawn using Origin 8.9.0 (version 2021, Origin Inc. 2021) software. The results of denoised ASVs were visualized using the Draw Venn Diagram platform (https://bioinformatics.psb.ugent.be/webtools/Venn/, accessed on 11 January 2023) to generate a Venn diagram illustrating the shared and unique ASVs of ECM fungal communities in the three families of P. massoniana. Redundancy analysis (RDA) was conducted using Canoco5.0 [45] software (Microcomputer Power, Inc., Ithaca, NY, USA) to study the effects of P. massoniana growth performance and rhizosphere soil nutrient content on the composition of ectomycorrhizal fungal communities.

3. Results

3.1. Rhizosphere Soil and Leaf Nutrients, Extracellular Enzyme Activities, and Ratio Characteristics of P. massoniana Families

The results revealed that the nitrogen and phosphorus content in the leaves of the HG family was significantly higher than that of the MG and LG families (Figure 1B,C). Specifically, the leaf nitrogen content (leaf N) in all three families ranged from 9.53 to 12.27 g/kg. Moreover, the patterns of soil organic carbon (SOC) and total nitrogen (TN) in the rhizosphere soil exhibited a similar trend, with the HG family displaying significantly higher content compared to the MG and LG families (Figure 1A,B). However, no significant differences were observed in terms of the content of leaf carbon (Leaf C) and total phosphorus (TN) in rhizosphere soil (Figure 1A,C).
In terms of nutrient ratios, only the carbon-to-phosphorus ratio (C/P) exhibited significant differences among the three families, as depicted in Figure 1D–F. Specifically, the rhizosphere soil carbon-to-phosphorus ratio (Soil C/P) was found to be higher in the HG family compared to the other two families. Conversely, the carbon-to-phosphorus ratio in the leaves (Leaf C/P) was considerably lower in the HG family in comparison to the MG family, with no significant difference observed in relation to the LG family. It was noteworthy that the nitrogen-to-phosphorus ratio (N/P) in the leaves of all three families were found to be below 10 (Figure 1F). This observation implied the existence of potential variations in nutrient allocation strategies among the families of P. massoniana.
Furthermore, the acid phosphatase (AP) activity was found to be the highest among the extracellular enzymes analyzed in this study, with values averaging around 2000 u−1·g−1·h−1 (Figure 2). Notably, the AP activity in the rhizosphere soil of the HG family exhibited a significantly higher level compared to the MG and LG families, with an increase of approximately 40%. A similar pattern was observed in the activities of β-1,4-N-acetylglucosaminidase (NAG), with a much greater increase observed in the HG family ranging from 77% to 129% (Figure 2). It is worth mentioning that the HG family also exhibited higher activity of β-1,4-glucosidase (BG) compared to the LG family but showed no significant difference compared to the MG family (Figure 2).

3.2. Community Composition of ECM Fungi

Based on the ASVs’ results, a Venn diagram was constructed to illustrate the distribution of ASVs among different families of P. massoniana. The diagram revealed that, in addition to 81 shared ASVs, the HG, MG, and LG families possessed 148, 132, and 132 unique ASVs, respectively (Figure 3). These findings highlighted the presence of unique microbial populations within each family. Additionally, the most dominant ectomycorrhizal (ECM) fungal genera in the HG family were Amphinema P. Karst. (30.5%) and Pseudoclathrosphaerina Voglmayr. (13.4%) (Figure 4). In the MG family, the ECM fungal genera Auricularia Bull. and Russula Pers. were found to be the most abundant, representing 14.0% and 10.2% of the total community, respectively. Russula was also dominant in ECM fungal community in the LG family, followed by Amanita Pers. (15.6%) (Figure 4).
Using simper analysis, the Bray–Curtis dissimilarity index of ECM fungal community composition among the three families of P. massoniana was decomposed. The results indicated that the main contributor to this dissimilarity index was Amphinema (contribution: HG-MG 0.225, HG-LG 0.17). Additionally, other genera with an abundance of less than 0.1% and unidentified genera also made substantial contributions (contribution: HG-MG 0.25 and HG-LG 0.225) (Figure 5 and Figure 6). These findings underscored the distinctive composition and contribution patterns exhibited by ECM fungal genera within the HG family.

3.3. Relationships between the Soil Extracellular Enzyme Activities, Nutrients, and Growth Performance of P. massoniana

Pearson correlation analysis showed significant correlations between the three soil enzyme activities and the growth performance, rhizosphere soil nutrients, and leaf nutrient contents (Figure 7). Extracellular enzyme activities in the rhizosphere soil displayed a moderate to strong positive correlation with the growth performance of P. massoniana, such as Height, DBH and V. There was no significant correlation observed between BG and NAG or between BG and AP. However, NAG exhibited a significant correlation with AP. Extracellular enzyme activities in the rhizosphere soil exhibited a non-significant correlation with leaf and rhizosphere soil nutrients for BG. However, AP and NAG demonstrated a significant positive correlation with Leaf P, SOC, and Soil N while displaying a negative correlation with Leaf C/P and Leaf N/P. Similar trends were also observed in the growth performance parameters. These findings implied that NAG and AP may have a more prominent role in nutrient dynamics.

3.4. Relationships between ECM Fungal Communities and Extracellular Enzyme Activity, Soil Nutrients, and Growth Performance

Redundancy analysis (RDA) was conducted with the ECM fungal abundance of P. massoniana as the response variable and soil chemical properties and tree growth performance as explanatory variables. The results showed that the first two axes of the RDA explained 44.16% of the total microbial variation (Figure 8). The Monte Carlo test indicated significant correlations (p < 0.05) between tree growth performance (Height, DBH, and V), soil total nitrogen (TN), soil N/P, and ECM fungal community composition (Table 2).
Specifically, the growth performance of P. massoniana exhibited a positive correlation with the following genera: Amphinema, Archaeorhizomyces Rosling and T.Y. James, Clavulina J. Schroöt., Pseudoclathrosphaerina, and Russula (Figure 8). Furthermore, soil TN and soil N/P demonstrated a positive correlation with Amanita, Tylospora Donk., Amphinema, Archaeorhizomyces, and Clavulina but displayed a negative correlation or no correlation with other genera.
Structural equation modeling was used to analyze the relationship between the dominant ECM genera and three extracellular enzyme activities. The results showed that the dominant genera Amphinema and Pseudoclathrosphaerina had a positive effect on the activities of the three enzymes (Figure 9).

4. Discussion

4.1. Effects of Different Growth Families of P. massoniana on Rhizosphere Soil and Leaf Nutrients

In this study, the range of leaf total nitrogen (leaf N) contents in P. massoniana was 9.53–12.27 g/kg (Figure 1B), lower than the global average of 20.06 g/kg [46] and the average of Chinese plants measured at 20.20 g/kg [47]. The N/P ratio has been widely used to predict nutrient limitation of plant growth [48]. All N/P ratios of P. massoniana leaves in this study (Figure 1F) were less than 10 (ranging from 2.56 to 2.95), indicating a N limitation [49] consistent with the findings of previous studies [50]. Although N/P ratios have been widely used for determining nutrient limitations in plants, an accurate prediction of the optimal N limitation for P. massoniana requires further investigation due to variations in plant N nutrient levels over time and differences in optimal N/P ratios among different tree species [51]. Among the three families of P. massoniana, the leaf total nitrogen and phosphorus contents of the HG family were significantly higher than those of the MG and LG families (Figure 1B,C) and positively correlated with tree growth performance to varying degrees. It has been previously reported that plants with a higher growth rate have lower leaf C/P ratios than slower-growing individuals [52]. An increase in the leaf C/P ratio indicates a decrease in the plant growth rate [53]. The significantly lower leaf C/P ratio in the HG family (Figure 1E) suggests that the higher growth rate of the HG family results in a large amount of nitrogen and phosphorus being stored in the leaves to support rapid growth, supporting the growth rate hypothesis [54].
In this study, the rhizosphere soil SOC content and three extracellular enzyme activities were significantly higher in the HG family than in the MG and LG families (Figure 1A and Figure 2) and the growth performance (height, DBH, and volume) of P. massoniana was significantly positively correlated with rhizosphere soil nutrients (SOC and TN) and the activities of three soil extracellular enzymes (BG, NAG and AP) (Figure 8). Previous studies have shown that there are significant differences in root biomass allocation patterns among different genotypes (provenances, clones, or cultivars) of the same species [55,56] and, as an important component of soil organic carbon (SOC) [57,58], root biomass plays an important role in the input and storage of soil carbon [46]. The three different growth families of P. massoniana selected in this study showed differences in height, DBH, and volume (Table 1) and, according to the “quality ratio effect”, taller and more extensive plant individuals have higher productivity and more complex underground functions than shorter and less extensive individuals, producing greater biomass [59,60,61]. Therefore, the high-growth family of P. massoniana could input more carbon into the rhizosphere soil, providing a large amount of substrate for extracellular enzymes, which in turn can affect their activities. The increase in extracellular enzyme activity can promote the release of nutrients in the rhizosphere soil and benefit the growth of P. massoniana. This study also showed that the rhizosphere soil C/P ratio of the HG family was significantly higher than that of the other two families and that the soil total nitrogen content of the HG family was significantly higher than that of the LG family. Previous studies have shown that a higher nitrogen content and C/P ratio can promote the accumulation of soil organic matter (SOM) [62]. SOM in forest soil mainly comes from plant inputs on the surface (plant litter) or underground (plant roots); soil organic carbon (SOC) is one of its main sources [63]. Therefore, it is speculated that the differences in rhizosphere soil organic carbon among the families with different growth levels are likely to be caused by their growth differences. The difference in rhizosphere soil C/P ratio between the HG family and the other two families also indicates that the mineralization ability of phosphorus varies between high-growth and low-growth families of P. massoniana [64]. The significant positive correlations between the rhizosphere soil C/P ratio and the height and DBH of P. massoniana confirm this inference.

4.2. Effects of Different Growth Families of P.massoniana on ECM Community Composition and Extracellular Enzyme Activity

Figure 4 revealed differences in the dominant ECM fungal genera and their abundances at the genus level among the three families, with different dominant genera for each. The dominant genus in the ECM fungal community of the high-growth family (HG) was Amphinema (30.5%) followed by Pseudoclathrosphaerina (13.4%); Auricularia (14.0%) was the dominant genus in the medium-growth family (MG) and Russula (17.2%) in the low-growth family (LG) of P. massoniana. Although the ECM fungal communities had similarities in their function as mycorrhizae, there were differences in the colonization habits and decomposition and utilization of soil nutrients among different ECM fungal species [65]. Moreover, differences in extracellular enzyme activities can also reflect differences in the composition of ECM fungal communities [66]. In this study, the NAG and AP extracellular enzyme activities in the rhizosphere soil of the HG family were significantly higher than those in the MG and LG families, while the BG activity was significantly higher in the LG family than in the MG family but not significantly different from that in the HG family (Figure 2), indicating a distinct ECM fungal community composition and stronger nutrient decomposition and release ability in the rhizosphere of the HG family which were beneficial for its growth. Furthermore, beyond the top 10 abundant ECM fungal genera shared among the three families of P. massoniana, the MG and LG families harbored other fungi at abundances lower than 0.1%, accounting for 63.53% and 46.80%, respectively, while the representation of other fungi with abundances less than 0.1% in the HG family was 29.30% (Figure 4), possibly related to the carbon source supply of the host plant to its own ECM mycorrhizae [2]. The superior growth performance of the HG family was attributed to its unique ECM fungal community, which was different from that of the MG and LG families, enabling it to achieve greater aboveground and belowground biomass. Therefore, compared with the MG and LG families, the HG family provided more colonization space and carbon sources (for fungi such as Amphinema and Pseudoclathrosphaerina) for the ECM fungal community, resulting in differences in the shared and dominant genera of ECM fungi among the three families with different growth levels. The RDA also showed that the growth performance (volume, height, and DBH) of P. massoniana was an important factor affecting the abundance of ECM fungi (Figure 8).
Differences in ECM fungal communities and extracellular enzyme activity in the rhizosphere soil of P. massoniana may reflect different nutrient acquisition modes. Prior research indicates that fungi with widespread hyphal growth provide less nitrogen to host plants [67], which may explain the higher NAG activity in the rhizosphere soil of the HG family. As a dominant fungal genus in the HG family, Amphinema possesses abundant external mycelia [68] and produces a large number of nitrogen-acquiring enzymes by oxidizing SOM to meet its own and the host plant’s nitrogen requirements. The positive effect of Amphinema on NAG in the structural equation model also confirms this mechanism (Figure 9). Compared to those of the MG and LG families, the ECM fungal community and extracellular enzyme activity in the HG family appear to better match the nutrient acquisition model proposed by Hupperts et al. [69].
In addition, the rhizosphere soil SOC of the HG family was significantly higher than that of the other two families (Figure 1), which may also be related to the impact of ECM fungal changes on soil carbon accumulation. In the study by Weight et al. [70], long-distance exploring ECM fungal groups, namely those capable of producing abundant external mycelia, can provide more biomass than short-distance exploring ECM fungi and their hyphae typically have hydrophobic surfaces and relatively low rhizosphere soil contact areas [71], resulting in their hyphae being more persistent than those of short-distance groups [72]. Therefore, in this study, the differences in the ECM fungal community structure of P. massoniana families with different growth levels may have had an impact on the input of fungal necromass into the soil carbon pool, where the ECM fungal community of the HG family, dominated by Amphinema, may be more conducive to SOC accumulation in the soil compared to the ECM fungal communities of the MG and LG families.
This study also revealed evidence of carbon sources for P. massoniana ECM fungal communities other than host plants. In the correlation analysis (Figure 7), rhizosphere soil organic carbon had varying degrees of positive correlation with three extracellular enzymes, as enzyme activity has a strong dependence on soil carbon and is sensitive to changes in soil SOC [23,73]. Although there was a certain positive correlation between organic carbon and the BG enzyme, it was not significantly associated with NAG and AP enzymes which may be related to the carbon source of soil microorganisms. The soil carbon pool is one of the carbon sources for soil microorganisms [74] where the extracellular enzymes BG and NAG can help microorganisms obtain carbon by decomposing different substrates. When cellulose decomposed by BG is the main carbon resource in the soil, microorganisms use more cellulose than chitin and peptidoglycan degraded by NAG. However, under conditions limited by environmental factors, such as chitin and peptidoglycan, cellulose would become the main substrate for microorganisms to utilize [75]. In this study, we explored the impact of nitrogen limitation on P. massoniana. We found that the NAG enzyme plays a critical role in catalyzing chitin decomposition into carbohydrates and inorganic nitrogen [76]. As a result, P. massoniana may release more NAG enzymes through root exudates, influencing ECM fungal communities and obtaining more nitrogen to alleviate its own nitrogen stress. Additionally, this alleviation mechanism provides a source of carbon for the ECM fungal community by decomposing a significant amount of chitin and peptidoglycan that NAG can break down.

5. Conclusions

This study explored the ECM fungal community composition, rhizosphere soil extracellular enzyme activity and nutrient characteristics, and growth performance of different growth families of P. massoniana. The high-growth family of P. massoniana displayed a higher leaf nutrient content (leaf N and leaf P), allowing it to cope with the high nutrient demand from its rapid growth. The application of appropriate nitrogen fertilizer can help alleviate P. massoniana’s nitrogen limitations in the study area. Furthermore, the high-growth family of P. massoniana displayed a higher abundance of ECM fungal communities that favor its growth and enable soil SOC accumulation. This family allocated more carbon sources to promote the secretion of more extracellular enzymes by ECM fungal communities, resulting in higher enzyme activity in the rhizosphere and an increase in phosphorus mineralization. In this way, the high-growth family provided a source of carbon for the ECM fungal communities beyond the host plants. Overall, this study offers important insights into the ECM fungal community structure, rhizosphere soil extracellular enzyme activity, and nutrient characteristics of P. massoniana with different growth levels. The findings provide guidance for a deeper understanding of the relationship between growth performance and the belowground environment of P. massoniana.

Author Contributions

J.Z.: field and lab experiments, formal analysis, writing—original draft, and writing—review and editing. Y.Z.: field and lab experiments and formal analysis. H.H.: collected the data. H.W.: revised and improved the manuscript. F.W.: methodology, project administration, writing—original draft, writing—review and editing, and supervision. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by the National Key Research and Development Program for Young Scientists (2022YFD2201700) and the National Natural Science Foundation of China under Grant (32260406, 31660201).

Data Availability Statement

The data supporting the conclusions of this article will be made available by authors, without undue reservation.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflict of interest.

References

  1. Kaiser, C.; Koranda, M.; Kitzler, B.; Fuchslueger, L.; Schnecker, J.; Schweiger, P.; Rasche, F.; Zechmeister-Boltenstern, S.; Sessitsch, A.; Richter, A. Belowground carbon allocation by trees drives seasonal patterns of extracellular enzyme activities by altering microbial community composition in a beech forest soil. New Phytol. 2010, 187, 843–858. [Google Scholar] [CrossRef] [Green Version]
  2. Van Der Heijden, M.G.; Bardgett, R.D.; Van Straalen, N.M. The unseen majority: Soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 2008, 11, 296–310. [Google Scholar] [CrossRef] [PubMed]
  3. Courty, P.E.; Pritsch, K.; Schloter, M.; Hartmann, A.; Garbaye, J. Activity profiling of ectomycorrhiza communities in two forest soils using multiple enzymatic tests. New Phytol. 2005, 167, 309–319. [Google Scholar] [CrossRef]
  4. Kyaschenko, J.; Clemmensen, K.E.; Karltun, E.; Lindahl, B.D. Below-ground organic matter accumulation along a boreal forest fertility gradient relates to guild interaction within fungal communities. Ecol. Lett. 2017, 20, 1546–1555. [Google Scholar] [CrossRef]
  5. Morrison, E.W.; Frey, S.D.; Sadowsky, J.J.; Van Diepen, L.T.; Thomas, W.K.; Pringle, A. Chronic nitrogen additions fundamentally restructure the soil fungal community in a temperate forest. Fungal Ecol. 2016, 23, 48–57. [Google Scholar] [CrossRef]
  6. Hasselquist, N.J.; Högberg, P. Dosage and duration effects of nitrogen additions on ectomycorrhizal sporocarp production and functioning: An example from two N-limited boreal forests. Ecol. Evol. 2014, 4, 3015–3026. [Google Scholar] [CrossRef] [PubMed]
  7. Pérez-Izquierdo, L.; Zabal-Aguirre, M.; González-Martínez, S.C.; Buée, M.; Verdú, M.; Rincón, A.; Goberna, M. Plant intraspecific variation modulates nutrient cycling through its below ground rhizospheric microbiome. J. Ecol. 2019, 107, 1594–1605. [Google Scholar] [CrossRef]
  8. Arfi, Y.; Buee, M.; Marchand, C.; Levasseur, A.; Record, E. Multiple markers pyrosequencing reveals highly diverse and host-specific fungal communities on the mangrove trees Avicennia marina and Rhizophora stylosa. FEMS Microbiol. Ecol. 2012, 79, 433–444. [Google Scholar] [CrossRef] [PubMed]
  9. Hagenbo, A.; Piñuela, Y.; Castaño, C.; Martínez De Aragón, J.; De-Miguel, S.; Alday, J.G.; Bonet, J.A. Production and turnover of mycorrhizal soil mycelium relate to variation in drought conditions in Mediterranean Pinus pinaster, Pinus sylvestris and Quercus ilex forests. New Phytol. 2021, 230, 1609–1622. [Google Scholar] [CrossRef]
  10. Molina, R.J.M.F.a.I.P.-F.P. Specificity phenomena in mycorrhizal symbioses: Community-ecological consequences and practical implications. In Mycorrhizal Functioning: An Integral Plant–Fungal Process; Allen, M.F., Ed.; Chapman & Hall: New York, NY, USA, 1992; pp. 357–423. [Google Scholar]
  11. Dong, H.; Ge, J.; Sun, K.; Wang, B.; Xue, J.; Wakelin, S.A.; Wu, J.; Sheng, W.; Liang, C.; Xu, Q. Change in root-associated fungal communities affects soil enzymatic activities during Pinus massoniana forest development in subtropical China. For. Ecol. Manag. 2021, 482, 118817. [Google Scholar] [CrossRef]
  12. Toju, H.; Sato, H.; Tanabe, A.S. Diversity and spatial structure of belowground plant–fungal symbiosis in a mixed subtropical forest of ectomycorrhizal and arbuscular mycorrhizal plants. PLoS ONE 2014, 9, e86566. [Google Scholar] [CrossRef]
  13. Carriconde, F.; Gardes, M.; Bellanger, J.-M.; Letellier, K.; Gigante, S.; Gourmelon, V.; Ibanez, T.; Mccoy, S.; Goxe, J.; Read, J. Host effects in high ectomycorrhizal diversity tropical rainforests on ultramafic soils in New Caledonia. Fungal Ecol. 2019, 39, 201–212. [Google Scholar] [CrossRef]
  14. Tedersoo, L.; Bahram, M. Mycorrhizal types differ in ecophysiology and alter plant nutrition and soil processes. Biol. Rev. 2019, 94, 1857–1880. [Google Scholar] [CrossRef] [PubMed]
  15. Deforest, J.L.; Moorhead, D.L. Effects of elevated pH and phosphorus fertilizer on soil C, N and P enzyme stoichiometry in an acidic mixed mesophytic deciduous forest. Soil Biol. Biochem. 2020, 150, 107996. [Google Scholar] [CrossRef]
  16. Zhou, L.; Liu, S.; Shen, H.; Zhao, M.; Xu, L.; Xing, A.; Fang, J. Soil extracellular enzyme activity and stoichiometry in China’s forests. Funct. Ecol. 2020, 34, 1461–1471. [Google Scholar] [CrossRef]
  17. Pellitier, P.T.; Zak, D.R. Ectomycorrhizal fungi and the enzymatic liberation of nitrogen from soil organic matter: Why evolutionary history matters. New Phytol. 2018, 217, 68–73. [Google Scholar] [CrossRef] [Green Version]
  18. Sinsabaugh, R.L.; Hill, B.H.; Follstad Shah, J.J. Ecoenzymatic stoichiometry of microbial organic nutrient acquisition in soil and sediment. Nature 2009, 462, 795–798. [Google Scholar] [CrossRef]
  19. Guo, K.; Zhao, Y.; Liu, Y.; Chen, J.; Wu, Q.; Ruan, Y.; Li, S.; Shi, J.; Zhao, L.; Sun, X.J.G. Pyrolysis temperature of biochar affects ecoenzymatic stoichiometry and microbial nutrient-use efficiency in a bamboo forest soil. Geoderma 2020, 363, 114162. [Google Scholar] [CrossRef]
  20. Sinsabaugh, R.L.; Lauber, C.L.; Weintraub, M.N.; Ahmed, B.; Allison, S.D.; Crenshaw, C.; Contosta, A.R.; Cusack, D.; Frey, S.; Gallo, M.E. Stoichiometry of soil enzyme activity at global scale. Ecol. Lett. 2008, 11, 1252–1264. [Google Scholar] [CrossRef] [Green Version]
  21. Burns, R.G.; Deforest, J.L.; Marxsen, J.; Sinsabaugh, R.L.; Stromberger, M.E.; Wallenstein, M.D.; Weintraub, M.N.; Zoppini, A. Soil enzymes in a changing environment: Current knowledge and future directions. Soil Biol. Biochem. 2013, 58, 216–234. [Google Scholar] [CrossRef]
  22. Nannipieri, P.; Trasar-Cepeda, C.; Dick, R.P. Soil enzyme activity: A brief history and biochemistry as a basis for appropriate interpretations and meta-analysis. Biol. Fertil. Soils 2018, 54, 11–19. [Google Scholar] [CrossRef]
  23. Wang, C.; Xue, L.; Jiao, R. Soil organic carbon fractions, C-cycling associated hydrolytic enzymes, and microbial carbon metabolism vary with stand age in Cunninghamia lanceolate (Lamb.) Hook plantations. For. Ecol. Manag. 2021, 482, 118887. [Google Scholar] [CrossRef]
  24. Albornoz, F.E.; Teste, F.P.; Lambers, H.; Bunce, M.; Murray, D.C.; White, N.E.; Laliberté, E. Changes in ectomycorrhizal fungal community composition and declining diversity along a 2-million-year soil chronosequence. Mol. Ecol. 2016, 25, 4919–4929. [Google Scholar] [CrossRef]
  25. Ning, P.; Cheng, X.M.; Yang, X.F.; Huang, X.X. Chemical properties and extracellular enzymatic activity in the rhizosphere soil of Abies fabri at different altitudes on Mount Gongga. Chin. J. Appl. Environ. Biol. 2021, 27, 1138–1146. (In Chinese) [Google Scholar]
  26. Welc, M.; Frossard, E.; Egli, S.; Bünemann, E.K.; Jansa, J. Rhizosphere fungal assemblages and soil enzymatic activities in a 110-years alpine chronosequence. Soil Biol. Biochem. 2014, 74, 21–30. [Google Scholar] [CrossRef]
  27. Nicholson, B.A.; Jones, M.D. Early-successional ectomycorrhizal fungi effectively support extracellular enzyme activities and seedling nitrogen accumulation in mature forests. Mycorrhiza 2017, 27, 247–260. [Google Scholar] [CrossRef] [PubMed]
  28. Maillard, F.; Didion, M.; Fauchery, L.; Bach, C.; Buée, M. N-Acetylglucosaminidase activity, a functional trait of chitin degradation, is regulated differentially within two orders of ectomycorrhizal fungi: Boletales and Agaricales. Mycorrhiza 2018, 28, 391–397. [Google Scholar] [CrossRef] [PubMed]
  29. Nicolás, C.; Martin-Bertelsen, T.; Floudas, D.; Bentzer, J.; Tunlid, A. The soil organic matter decomposition mechanisms in ectomycorrhizal fungi are tuned for liberating soil organic nitrogen. ISME J. 2018, 13, 977–988. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, Z.M.; Zhang, J.C.; He, Y.S.; Zhang, Y.W.; Liu, Y.Y. The research progress on ectomycorrhizal fungi in Pinus massoniana. Tillage Cultiv. 2016, 2, 66–68+72. (In Chinese) [Google Scholar]
  31. Luo, Y.H.; Sun, D.J.; Lin, J.Y.; Guo, W.F.; Lu, L.H.; Wen, Y.G. Effect of Close-to-Nature management on the natural regeneration and species diversity in a massson pine plantation. Acta Ecol. Sin. 2013, 33, 6154–6162. (In Chinese) [Google Scholar]
  32. Huang, S.Y.; Wu, T.G.; Chu, X.L.; Wang, B.; Wang, X.H.; Zhang, D.B.; Zhou, Z.C. Effects of phosphorus addition and inoculation of mycorrhizal fungi on the growth and phosphorus utilization of masson pine container seedlings from different families. For. Res. 2021, 34, 142–151. (In Chinese) [Google Scholar]
  33. Zhao, Y.; Zhou, Z.C.; Wu, J.F.; Lan, Y.Z. Growth response of superior provenance of masson pine to phosphorus supply and persistence of phosphorus effect. Sci. Silvae Sin. 2007, 43, 64–70. (In Chinese) [Google Scholar]
  34. Jian, Z.; Ni, Y.; Zeng, L.; Lei, L.; Xu, J.; Xiao, W.F.; Li, M.-H. Latitudinal patterns of soil extracellular enzyme activities and their controlling factors in Pinus massoniana plantations in subtropical China. For. Ecol. Manag. 2021, 495, 119358. [Google Scholar] [CrossRef]
  35. Bao, S. Soil and Agricultural Chemistry Analysis; China Agriculture Press: Beijing, China, 2000. [Google Scholar]
  36. Mayor, Á.G.; Goirán, S.B.; Vallejo, V.R.; Bautista, S. Variation in soil enzyme activity as a function of vegetation amount, type, and spatial structure in fire-prone Mediterranean shrublands. Sci. Total Environ. 2016, 573, 1209–1216. [Google Scholar] [CrossRef]
  37. German, D.P.; Weintraub, M.N.; Grandy, A.S.; Lauber, C.L.; Rinkes, Z.L.; Allison, S.D. Optimization of hydrolytic and oxidative enzyme methods for ecosystem studies. Soil Biol. Biochem. 2011, 43, 1387–1397. [Google Scholar] [CrossRef]
  38. Allen, G.C.; Flores-Vergara, M.; Krasynanski, S.; Kumar, S.; Thompson, W. A modified protocol for rapid DNA isolation from plant tissues using cetyltrimethylammonium bromide. Nat. Protoc. 2006, 1, 2320–2325. [Google Scholar] [CrossRef]
  39. White, T.J.; Bruns, T.; Lee, S.; Taylor, J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. PCR Protoc. A Guide Methods Appl. 1990, 18, 315–322. [Google Scholar]
  40. Magoč, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Haas, B.J.; Gevers, D.; Earl, A.M.; Feldgarden, M.; Ward, D.V.; Giannoukos, G.; Ciulla, D.; Tabbaa, D.; Highlander, S.K.; Sodergren, E. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011, 21, 494–504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Callahan, B.J.; Mcmurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [Green Version]
  43. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Gregory Caporaso, J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 2018, 6, 90. [Google Scholar] [CrossRef] [PubMed]
  44. Warton, D.I.; Wright, S.T.; Wang, Y. Distance-based multivariate analyses confound location and dispersion effects. Methods Ecol. Evol. 2012, 3, 89–101. [Google Scholar] [CrossRef]
  45. Šmilauer, P.; Lepš, J. Multivariate Analysis of Ecological Data Using Canono 5, 2nd ed.; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  46. Elser, J.J.; Fagan, W.F.; Denno, R.F.; Dobberfuhl, D.R.; Folarin, A.; Huberty, A.; Interlandi, S.; Kilham, S.S.; Mccauley, E.; Schulz, K.L. Nutritional constraints in terrestrial and freshwater food webs. Nature 2000, 408, 578–580. [Google Scholar] [CrossRef] [PubMed]
  47. Han, W.; Fang, J.; Guo, D.; Zhang, Y. Leaf nitrogen and phosphorus stoichiometry across 753 terrestrial plant species in China. New Phytol. 2005, 168, 377–385. [Google Scholar] [CrossRef]
  48. Koerselman, W.; Meuleman, A.F. The vegetation N: P ratio: A new tool to detect the nature of nutrient limitation. J. Appl. Ecol. 1996, 33, 1441–1450. [Google Scholar] [CrossRef]
  49. Güsewell, S. N: P ratios in terrestrial plants: Variation and functional significance. New Phytol. 2004, 164, 243–266. [Google Scholar] [CrossRef] [PubMed]
  50. Kuang, Y.W.; Wen, D.Z.; Yan, J.H.; Liu, S.Z.; Chu, G.W.; Zhou, C.Y.; Wang, G.Q.; Zhang, Q.M. Characteristics of element contents in leaves of 3 dominant species in karst forest in Puding, Guizhou, China. Chin. J. Appl. Environ. Biol. 2010, 16, 158–163. (In Chinese) [Google Scholar] [CrossRef]
  51. Zhang, X.Y.; Ning, Q.R.; Li, S.Z.; Jiang, L.C.; Liu, Y.; Zhang, H.Z.; Tao, J.J. Stoichiometric characteristics of Pinus massoniana plantation in the subtropical red soil erosion region. Res. Soil Water Conserv. 2017, 24, 156–161. (In Chinese) [Google Scholar]
  52. Sterner, R.W.; Elser, J.J. Ecological Stoichiometry: The Biology of Elements from Molecules to the Biosphere; Princeton University Press: Princeton, NJ, USA, 2002. [Google Scholar]
  53. Bai, X.; Dippold, M.A.; An, S.; Wang, B.; Zhang, H.; Loeppmann, S. Extracellular enzyme activity and stoichiometry: The effect of soil microbial element limitation during leaf litter decomposition. Ecol. Indic. 2021, 121, 107200. [Google Scholar] [CrossRef]
  54. Aerts, R. Nutrient resorption from senescing leaves of perennials: Are there general patterns? J. Ecol. 1996, 84, 597–608. [Google Scholar] [CrossRef]
  55. Li, L.X.; Zheng, S.S.; Xu, J.W.; Wu, P.F. Research advance in influence mechanism of tree root biomass allocation. World For. Res. 2022, 35, 15–20. (In Chinese) [Google Scholar]
  56. Wu, P.F.; Ma, X.Q.; Chen, Y.L.; Lin, W.J.; Huang, S.Y.; Liu, L.Q. Comparison of phosphorus use efficiency among clonal test plantations of Chinese fir. J. Fujian Agric. For. Univ. (Nat. Sci. Ed.) 2012, 41, 40–45. (In Chinese) [Google Scholar]
  57. Mao, Y.L. Advance of research on fractionation methods of soil organic matter. J. Quanzhou Teach. Univ. (Nat. Sci.) 2008, 26, 70–75. (In Chinese) [Google Scholar]
  58. Christensen, B.T. Physical fractionation of soil and organic matter in primary particle size and density separates. In Advances in Soil Science; Springer: New York, NY, USA, 1992; Volume 20, pp. 1–90. [Google Scholar]
  59. Erfanzadeh, R.; Bahrami, B.; Motamedi, J.; Pétillon, J. Changes in soil organic matter driven by shifts in co-dominant plant species in a grassland. Geoderma 2014, 213, 74–78. [Google Scholar] [CrossRef] [Green Version]
  60. Gruba, P.; Socha, J. Exploring the effects of dominant forest tree species, soil texture, altitude, and pHH2O on soil carbon stocks using generalized additive models. For. Ecol. Manag. 2019, 447, 105–114. [Google Scholar] [CrossRef]
  61. Mariotte, P.; Buttler, A.; Kohler, F.; Gilgen, A.K.; Spiegelberger, T. How do subordinate and dominant species in semi-natural mountain grasslands relate to productivity and land-use change? Basic Appl. Ecol. 2013, 14, 217–224. [Google Scholar] [CrossRef]
  62. Gong, X.; Liu, Y.; Li, Q.; Wei, X.; Guo, X.; Niu, D.; Zhang, W.; Zhang, J.; Zhang, L. Sub-tropic degraded red soil restoration: Is soil organic carbon build-up limited by nutrients supply. For. Ecol. Manag. 2013, 300, 77–87. [Google Scholar] [CrossRef]
  63. Whalen, E.D.; Lounsbury, N.; Geyer, K.; Anthony, M.; Morrison, E.; Van Diepen, L.T.; Le Moine, J.; Nadelhoffer, K.; Vanden Enden, L.; Simpson, M.J. Root control of fungal communities and soil carbon stocks in a temperate forest. Soil Biol. Biochem. 2021, 161, 108390. [Google Scholar] [CrossRef]
  64. Xiang, Y.X.; Chen, S.K.; Pan, P.; Ouyang, X.Z.; Ning, J.K.; Li, Q. Stoichiometric characteristics of carbon, nitrogen and phosphorus in leave-litter-soil of Pinus massoniana. J. For. Environ. 2019, 39, 120–126. [Google Scholar]
  65. Smith, A.J.; Potvin, L.R.; Lilleskov, E.A. Fertility-dependent effects of ectomycorrhizal fungal communities on white spruce seedling nutrition. Mycorrhiza 2015, 25, 649–662. [Google Scholar] [CrossRef]
  66. Ning, C.; Xiang, W.; Mueller, G.M.; Egerton-Warburton, L.M.; Yan, W.; Liu, S. Differences in ectomycorrhizal community assembly between native and exotic pines are reflected in their enzymatic functional capacities. Plant Soil 2020, 446, 179–193. [Google Scholar] [CrossRef]
  67. Corrêa, A.; Gurevitch, J.; Martins-Loução, M.; Cruz, C. C allocation to the fungus is not a cost to the plant in ectomycorrhizae. Oikos 2012, 121, 449–463. [Google Scholar] [CrossRef]
  68. Huang, L.L.; Wang, Y.L.; Yu, F.Q. Exogenous competitive fungi and their mycorrhizal morphology during ectomycorrhizal synthesis. Mycosystema 2020, 39, 955–962. [Google Scholar]
  69. Hupperts, S.F.; Karst, J.; Pritsch, K.; Landhäusser, S.M. Host phenology and potential saprotrophism of ectomycorrhizal fungi in the boreal forest. Funct. Ecol. 2017, 31, 116–126. [Google Scholar] [CrossRef]
  70. Weigt, R.B.; Raidl, S.; Verma, R.; Agerer, R. Exploration type-specific standard values of extramatrical mycelium—A step towards quantifying ectomycorrhizal space occupation and biomass in natural soil. Mycol. Prog. 2012, 11, 287–297. [Google Scholar] [CrossRef]
  71. Koide, R.T.; Fernandez, C.; Malcolm, G. Determining place and process: Functional traits of ectomycorrhizal fungi that affect both community structure and ecosystem function. New Phytol. 2014, 201, 433–439. [Google Scholar] [CrossRef]
  72. Certano, A.K.; Fernandez, C.W.; Heckman, K.A.; Kennedy, P.G. The afterlife effects of fungal morphology: Contrasting decomposition rates between diffuse and rhizomorphic necromass. Soil Biol. Biochem. 2018, 126, 76–81. [Google Scholar] [CrossRef]
  73. Fanin, N.; Mooshammer, M.; Sauvadet, M.; Meng, C.; Alvarez, G.; Bernard, L.; Bertrand, I.; Blagodatskaya, E.; Bon, L.; Fontaine, S. Soil enzymes in response to climate warming: Mechanisms and feedbacks. Funct. Ecol. 2022, 36, 1378–1395. [Google Scholar] [CrossRef]
  74. Kramer, C.; Trumbore, S.; Fröberg, M.; Dozal, L.M.C.; Zhang, D.; Xu, X.; Santos, G.M.; Hanson, P.J. Recent (<4 year old) leaf litter is not a major source of microbial carbon in a temperate forest mineral soil. Soil Biol. Biochem. 2010, 42, 1028–1037. [Google Scholar]
  75. Mori, T.; Aoyagi, R.; Kitayama, K.; Mo, J. Does the ratio of β-1, 4-glucosidase to β-1, 4-N-acetylglucosaminidase indicate the relative resource allocation of soil microbes to C and N acquisition? Soil Biol. Biochem. 2021, 160, 108363. [Google Scholar] [CrossRef]
  76. Boerner, R.E.; Giai, C.; Huang, J.; Miesel, J.R. Initial effects of fire and mechanical thinning on soil enzyme activity and nitrogen transformations in eight North American forest ecosystems. Soil Biol. Biochem. 2008, 40, 3076–3085. [Google Scholar] [CrossRef]
Forests 14 01447 g001 550
Figure 1. Rhizosphere soil and leaf nutrients and their ratios in P. massoniana families with different growth levels. Note: mean ± standard error, different lowercase letters indicated significant differences (p < 0.05) based on Duncan’s multiple range test. Note: TN, total nitrogen of the rhizosphere soil; TP, total phosphorus of the rhizosphere soil; C, carbon; C/P, carbon-to-phosphorus ratio; N/P, nitrogen-to-phosphorus ratio; C/N, carbon-to-nitrogen ratio; HG, high growth; MG, medium growth; and LG, low growth. The subgraphs letters denote (A) SOC and C (carbon); (B) TN, total nitrogen of the rhizosphere soil; (C) TP, total phosphorus of the rhizosphere soil; (D) C/N, carbon-to-nitrogen ratio; (E) C/P, carbon-to-phosphorus ratio; (F) N/P, nitrogen-to-phosphorus ratio.
Figure 1. Rhizosphere soil and leaf nutrients and their ratios in P. massoniana families with different growth levels. Note: mean ± standard error, different lowercase letters indicated significant differences (p < 0.05) based on Duncan’s multiple range test. Note: TN, total nitrogen of the rhizosphere soil; TP, total phosphorus of the rhizosphere soil; C, carbon; C/P, carbon-to-phosphorus ratio; N/P, nitrogen-to-phosphorus ratio; C/N, carbon-to-nitrogen ratio; HG, high growth; MG, medium growth; and LG, low growth. The subgraphs letters denote (A) SOC and C (carbon); (B) TN, total nitrogen of the rhizosphere soil; (C) TP, total phosphorus of the rhizosphere soil; (D) C/N, carbon-to-nitrogen ratio; (E) C/P, carbon-to-phosphorus ratio; (F) N/P, nitrogen-to-phosphorus ratio.
Forests 14 01447 g001
Forests 14 01447 g002 550
Figure 2. Extracellular enzyme activities in the rhizosphere soil of P. massoniana families with different growth levels. Note: mean ± standard error, where different lowercase letters indicate significant differences (p < 0.05) based on Duncan’s multiple range test. Note: BG, β-1,4-glucosidase; NAG, β-1,4-N-acetylglucosaminidase; AP, acid phosphatase; HG, high growth; MG, medium growth; and LG, low growth.
Figure 2. Extracellular enzyme activities in the rhizosphere soil of P. massoniana families with different growth levels. Note: mean ± standard error, where different lowercase letters indicate significant differences (p < 0.05) based on Duncan’s multiple range test. Note: BG, β-1,4-glucosidase; NAG, β-1,4-N-acetylglucosaminidase; AP, acid phosphatase; HG, high growth; MG, medium growth; and LG, low growth.
Forests 14 01447 g002
Forests 14 01447 g003 550
Figure 3. ASV Venn diagram of P. massoniana families with different growth levels. Note: HG, high growth; MG, medium growth; and LG, low growth.
Figure 3. ASV Venn diagram of P. massoniana families with different growth levels. Note: HG, high growth; MG, medium growth; and LG, low growth.
Forests 14 01447 g003
Forests 14 01447 g004 550
Figure 4. Relative abundance distribution of ECM fungal communities (genus level) of P. massoniana families with different growth levels. Note: All were the top 10 genera with average abundance > 0.1%. Note: HG, high growth; MG, medium growth; and LG, low growth.
Figure 4. Relative abundance distribution of ECM fungal communities (genus level) of P. massoniana families with different growth levels. Note: All were the top 10 genera with average abundance > 0.1%. Note: HG, high growth; MG, medium growth; and LG, low growth.
Forests 14 01447 g004
Forests 14 01447 g005 550
Figure 5. Contribution map of ECM fungal community composition in the HG and MG families of P. massoniana. Note: HG, high growth and MG, medium growth.
Figure 5. Contribution map of ECM fungal community composition in the HG and MG families of P. massoniana. Note: HG, high growth and MG, medium growth.
Forests 14 01447 g005
Forests 14 01447 g006 550
Figure 6. Contribution map of ECM fungi community composition in the HG and LG families of P. massoniana. Note: HG, high growth and LG, low growth.
Figure 6. Contribution map of ECM fungi community composition in the HG and LG families of P. massoniana. Note: HG, high growth and LG, low growth.
Forests 14 01447 g006
Forests 14 01447 g007 550
Figure 7. Pearson correlation test. Note: V, timber volume; Height, tree height; DBH, diameter at breast height; BG, β-1,4-glucosidase; NAG, β-1,4-N-acetylglucosaminidase; AP, acid phosphatase; TN, total nitrogen of the rhizosphere soil; TP, total phosphorus of the rhizosphere soil; C, carbon; C/P, carbon-to-phosphorus ratio; N/P, nitrogen-to-phosphorus ratio; C/N, carbon-to-nitrogen ratio. * represents p < 0.05, ** represents p < 0.01, *** represents p < 0.001.
Figure 7. Pearson correlation test. Note: V, timber volume; Height, tree height; DBH, diameter at breast height; BG, β-1,4-glucosidase; NAG, β-1,4-N-acetylglucosaminidase; AP, acid phosphatase; TN, total nitrogen of the rhizosphere soil; TP, total phosphorus of the rhizosphere soil; C, carbon; C/P, carbon-to-phosphorus ratio; N/P, nitrogen-to-phosphorus ratio; C/N, carbon-to-nitrogen ratio. * represents p < 0.05, ** represents p < 0.01, *** represents p < 0.001.
Forests 14 01447 g007
Forests 14 01447 g008 550
Figure 8. Relationship between the top 10 ECM fungi (>0.01%) in terms of relative abundance, soil nutrients, and growth performance of P. massoniana. Note: Clavulina (Clavulin), Amphinema (Amphinem), Auricularia (Auriculr), Tylospora (Tylospor), Archaeorhizomyces (Archaeor), Lactarius (Lactariu), Pseudoclathrosphaerina (Pseudocl); V, timber volume; Height, tree height; DBH, diameter at breast height; TN, total nitrogen of the rhizosphere soil; TP, total phosphorus of the rhizosphere soil; C/P, carbon-to-phosphorus ratio; N/P, nitrogen-to-phosphorus ratio; C/N, carbon-to-nitrogen ratio.
Figure 8. Relationship between the top 10 ECM fungi (>0.01%) in terms of relative abundance, soil nutrients, and growth performance of P. massoniana. Note: Clavulina (Clavulin), Amphinema (Amphinem), Auricularia (Auriculr), Tylospora (Tylospor), Archaeorhizomyces (Archaeor), Lactarius (Lactariu), Pseudoclathrosphaerina (Pseudocl); V, timber volume; Height, tree height; DBH, diameter at breast height; TN, total nitrogen of the rhizosphere soil; TP, total phosphorus of the rhizosphere soil; C/P, carbon-to-phosphorus ratio; N/P, nitrogen-to-phosphorus ratio; C/N, carbon-to-nitrogen ratio.
Forests 14 01447 g008
Forests 14 01447 g009 550
Figure 9. Structural equation model of the dominant genera and extracellular enzyme activities of P. massoniana. Note: x2 = 8.752, DF = 9, p = 0.460, GFI = 0.925, RMSEA = 0.000. The red and blue arrows represent positive and negative pathways, respectively. The numbers adjacent to the arrows are standardized path coefficients. Note: BG, β-1,4-glucosidase; NAG, β-1,4-N-acetylglucosaminidase; AP, acid phosphatase.
Figure 9. Structural equation model of the dominant genera and extracellular enzyme activities of P. massoniana. Note: x2 = 8.752, DF = 9, p = 0.460, GFI = 0.925, RMSEA = 0.000. The red and blue arrows represent positive and negative pathways, respectively. The numbers adjacent to the arrows are standardized path coefficients. Note: BG, β-1,4-glucosidase; NAG, β-1,4-N-acetylglucosaminidase; AP, acid phosphatase.
Forests 14 01447 g009
Table
Table 1. Growth statistics of family sample trees with different growth rates. Note: mean ± standard error where different lowercase letters indicate significant differences (p < 0.05) based on Duncan’s multiple range test. V = 0.000094602 × (D(1.88156−0.0030651×(D+H))) × (H(0.7684+0.0046574×(D+H))).
Table 1. Growth statistics of family sample trees with different growth rates. Note: mean ± standard error where different lowercase letters indicate significant differences (p < 0.05) based on Duncan’s multiple range test. V = 0.000094602 × (D(1.88156−0.0030651×(D+H))) × (H(0.7684+0.0046574×(D+H))).
Parameter ListHigh Growth
(HG)		Medium Growth
(MG)		Low Growth
(LG)
Height/(m)9.37 ± 0.57 a8.90 ± 0.13 a4.77 ± 0.20 b
DBH/(cm)17.62 ± 0.39 a16.60 ± 0.15 a9.59 ± 0.94 b
V/(m3)0.122 ± 0.007 a0.105 ± 0.003 b0.023 ± 0.004 c
year999
Table
Table 2. Table 2. Conditional effects of soil nutrients and growth performance on ECM fungal community composition. Explanation (%): the explanation rate of a single variable relative to the total variation. Note: V, timber volume; Height, tree height; DBH, diameter at breast height; TN, total nitrogen of the rhizosphere soil; TP, total phosphorus of the rhizosphere soil; C/P, carbon-to-phosphorus ratio; N/P, nitrogen-to-phosphorus ratio; C/N, carbon-to-nitrogen ratio. * p < 0.05, ** p < 0.01.
Table 2. Table 2. Conditional effects of soil nutrients and growth performance on ECM fungal community composition. Explanation (%): the explanation rate of a single variable relative to the total variation. Note: V, timber volume; Height, tree height; DBH, diameter at breast height; TN, total nitrogen of the rhizosphere soil; TP, total phosphorus of the rhizosphere soil; C/P, carbon-to-phosphorus ratio; N/P, nitrogen-to-phosphorus ratio; C/N, carbon-to-nitrogen ratio. * p < 0.05, ** p < 0.01.
VariablesExplanation %F Valuep Value
V195.90.002 **
Height18.87.30.002 **
DBH6.83.10.026 *
TN83.40.012 *
TP3.51.90.11
SOC0.1<0.10.992
N/P7.33.80.008 *
C/P1.30.70.5
C/N0.60.30.82
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, J.; Zhao, Y.; He, H.; Wang, H.; Wu, F. Extracellular Enzyme Activity and Nutrient Characteristics of Pinus massoniana Lamb. Families with Different Growth Levels: Insights into the Ectomycorrhizal Fungal Community and Rhizosphere Soil. Forests 2023, 14, 1447. https://doi.org/10.3390/f14071447

AMA Style

Zhang J, Zhao Y, He H, Wang H, Wu F. Extracellular Enzyme Activity and Nutrient Characteristics of Pinus massoniana Lamb. Families with Different Growth Levels: Insights into the Ectomycorrhizal Fungal Community and Rhizosphere Soil. Forests. 2023; 14(7):1447. https://doi.org/10.3390/f14071447

Chicago/Turabian Style

Zhang, Jun, Yuanxiang Zhao, Hongyang He, Haoyun Wang, and Feng Wu. 2023. "Extracellular Enzyme Activity and Nutrient Characteristics of Pinus massoniana Lamb. Families with Different Growth Levels: Insights into the Ectomycorrhizal Fungal Community and Rhizosphere Soil" Forests 14, no. 7: 1447. https://doi.org/10.3390/f14071447

APA Style

Zhang, J., Zhao, Y., He, H., Wang, H., & Wu, F. (2023). Extracellular Enzyme Activity and Nutrient Characteristics of Pinus massoniana Lamb. Families with Different Growth Levels: Insights into the Ectomycorrhizal Fungal Community and Rhizosphere Soil. Forests, 14(7), 1447. https://doi.org/10.3390/f14071447

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop