Yong Zhang
Qingcheng Wang*- Key Laboratory of Sustainable Forest Ecosystem Management-Ministry of Education, School of Forestry, Northeast Forestry University, Harbin, China
In comparison with ectomycorrhizal (EM) tree species, arbuscular mycorrhizal (AM) trees have different litter quality and nitrogen cycle modes, which may affect mycorrhizal colonization and the community composition and diversity. However, available studies addressing the mycorrhizal fungal colonization rate, diversity and community composition in mixed forest stands composed of AM and EM trees are rare. In the present study, we assessed litter quality, soil physicochemical properties and correlated them with mycorrhizal community characteristics in rhizosphere soils of monoculture and mixture plantation stands of AM tree species (Fraxinus mandschurica Rupr.) and EM tree species (Larix gmelinii Rupr., Picea koraiensis Nakai) in Northeast China. We hypothesized that (1) the effect of mixture pattern on mycorrhizal colonization rate and diversity would change with tree species, (2) the effect of mixture pattern on mycorrhizal community composition would be less pronounced in comparison with that of tree species. We found that mixture did not change AMF colonization rate regardless of mixture identity, whereas mixture and tree species exerted significant effects on EMF colonization rate. For AMF community, both M-AS (Fraxinus mandschurica Rupr. and Picea koraiensis Nakai) and M-AL (Fraxinus mandschurica Rupr. and Larix gmelinii Rupr.) mixtures significantly increased Pielou index and Simpson index, whereas only M-AS significantly increased Sobs. For EMF community, mixture significantly affected examined diversity indices except for Chao1. Mixture significantly shifted AMF and EMF community, and the magnitude was tree species dependent. The dominant genera in AMF and EMF communities in plantation stands were Glomus and Tomentella, respectively. The EnvFit analysis showed that the determinant factors of EMF community are soil moisture, pH, nitrate nitrogen content, dissolved organic nitrogen content, soil organic matter content, soil organic carbon/total nitrogen and litter carbon/total nitrogen. In conclusion, mixed conifer-broadleaf trees significantly changed soil physicochemical properties, litter quality as well as mycorrhizal fungi community diversity and composition.
Introduction
A mycorrhizal association is a symbiotic relationship that occurs between plant roots and soil fungi that increases plant performance by enabling plants to extract nitrogen (N) and phosphorus (P) more efficiently from the soil, improving plant productivity and resistance to pathogens, which is extensively distributed in forest ecosystems (Van der Heijden et al., 2015). Arbuscular mycorrhizal fungi (AMF) or ectomycorrhizal fungi (EMF) form a symbiotic relationship with almost every tree species (Brundrett and Tedersoo, 2018; Averill et al., 2022). As a vital component of N cycling in ecosystems, AMF and EMF are responsible for 20 and 80% of plant N uptake, respectively, (Leigh et al., 2009; Van der Heijden et al., 2015). However, AMF and EMF have completely distinct strategies for utilizing N. For example, the AMF community, primarily relies on inorganic N as a source of N in soil (Phillips et al., 2013), whereas the EMF can utilize organic N directly from decomposing organic matter (Read and Perez-Moreno, 2003; Van der Heijden et al., 2015). In addition, AMF and EMF communities respond differently to soil N availability (Phillips et al., 2013; Aldrich-Wolfe et al., 2020).
According to the mycorrhizal type, the trees can be classified into arbuscular mycorrhizal (AM) and ectomycorrhizal (EM) trees. Different mycorrhizal tree species produce different litter qualities that affect soil N content (Phillips et al., 2013). AM tree species generally have a quality litter of a lower carbon to nitrogen ratio (C/N) than EM ones (Phillips et al., 2013; Albornoz et al., 2016; Taylor et al., 2016; Lin et al., 2017; Seyfried et al., 2021; Prada et al., 2022). As proposed by the mycorrhizal-associated nutrient economy theory, forests of AM trees dominate by an inorganic N cycle (Phillips et al., 2013), characterizing by higher soil inorganic N concentration and higher N cycling rate (Lin et al., 2017). Numerous studies found that content of N not only affects colonization rate (Johnson et al., 2003; ÖPik et al., 2010; Sun et al., 2010; Witte et al., 2017; Frater et al., 2018; Lu et al., 2020; Yang et al., 2021a), but also richness, diversity and community composition of AMF (Roy et al., 2017; Witte et al., 2017; Aldrich-Wolfe et al., 2020; Luo et al., 2020; Yang et al., 2021a; Boeraeve et al., 2022). In general, enhanced N availability induce a decline in mycorrhizal association, because of a reduction in plant carbon allocation to root and soil microorganisms (Högberg et al., 2021; Ma et al., 2021).
The plantations acreage in China is the largest in the world. Because of monoculture and high density planting, most plantations are encountering ecological problems, such as soil degradation, productivity reductions, and biodiversity loss (Wang et al., 2013; Felton et al., 2016; Pretzsch et al., 2017; Yang et al., 2017; Wu et al., 2019). In contrast, mixed plantations increase soil microbial diversity and enzyme activity, facilitate litter decomposition, soil fertility improvement and forest productivity (Seidel et al., 2013; Forrester, 2014; Lang et al., 2014; Wu et al., 2019). Earlier studies found that the spore density, community diversity, and abundance of AM fungi in mixed forests were significantly higher than those in monoculture forests (ÖPik et al., 2010; Li et al., 2021a). In contrast, mixed Juglans mandshurica and Larix gmelinii did not change AM fungal richness and diversity, but significantly affected AM fungal community structure and composition (Ji et al., 2021). To date, most of the available studies focus on assessing the effects of tree monoculture on mycorrhizal fungal colonization and communities associated with a single mycorrhiza type (Aldrich-Wolfe et al., 2020; Yang et al., 2021a; Wang and Han, 2022). In contrast to the growing interest in comparing AM and EM in terms of their functions and responses to environmental factors (Kilpeläinen et al., 2019), there is still a great knowledge gap about the potential effects of tree mixture comprising of AM and EM tree species on mycorrhizal fungal community.
The objective of this study was to determine whether and how soil nutrient differences and litter quality due to mixing AM with EM tree species would affect mycorrhizal fungi communities in 33-year-old mixture and monoculture plantation forest stands of Fraxinus mandschurica Rupr., Larix gmelinii Rupr., and Picea koraiensis Nakai (Jiang et al., 2020; Zarif et al., 2020). We hypothesized that (1) the effect of mixture pattern on mycorrhizal colonization rate and diversity would change with tree species, (2) the effect of mixture pattern on mycorrhizal community composition would be less pronounced in comparison with that of tree species. We further assumed that colonization rates, richness and diversity of AMF communities in the mixture would increase whereas the colonization rate, richness and diversity of EMF in the mixture would decline because of litter quality and soil N availability. Meanwhile, we predicted that AM and EM fungal community composition in the mixture would shift due to changes in soil physicochemical properties and litter quality.
Materials and methods
Site description
A study was conducted at the Maoershan Experimental Station of Northeast Forestry University, Heilongjiang province, China (45°25′24′′N, 127°38′56′′E, 400–500 m elevation; Supplementary Figure S1), which is located at the branch range of the Changbai Mountains, with a landform feature of low mountains and hills. This region has a temperate continental monsoon climate with annual average temperatures of 2.8°C and monthly average temperatures of −19.6°C and 20.9°C. The annual average precipitation is approximately 600–800 mm, and there are 120–140 days without frost each year (Zarif et al., 2020). Soils at the site are Hap-Boric Luvisol, which includes loamy soil with a depth of 1–10 cm and sandy loam with a depth of 10–20 cm (Yang et al., 2015). This forest is considered to be a secondary forest. The dominant species are Betula platyphylla, Fraxinus mandshurica, Picea koraiensis, Populus davidiana, Pinus sylvestris var. mongolica, Larix gmelini, Juglans mandshurica and Quercus mongolica (Yang et al., 2010).
Field sampling
Sampling was conducted in five plantation stands, including Ash (Fraxinus mandschurica Rupr.) monoculture (P-A), Larch (Larix gmelinii Rupr.) monoculture (P-L), Spruce (Picea koraiensis Nakai) monoculture (P-S), Ash and Larch mixture (M-AL), Ash and spruce mixture (M-AS). These stands were first established in 1987 by planting on a strip cutover site of a secondary forest with similar site conditions. The initial spacing between two adjacent trees was 1.5 m × 2.0 m. In the mixtures, broad-leafed trees planted in 3 rows and coniferous ones planted in 5 rows alternatively. There are three replications for each stand. The current features of the stands are given in Table 1. In July 2020, the plots were set up according to the method as described by Luo et al. (2020), with appropriate adjustments. Each replication was considered a 30 m × 50 m plot, 50 m apart from one another. In each plot, nine medium-sized trees, each 10 m apart, were randomly selected from a monoculture, and nine trees adjacent to another tree species were randomly selected for each species in the mixture. Three litter collecting nets were set up 5 m away from selected trees in each plot, and litter was collected in October 2020.
Table 1. The features of AM, EM tree species monoculture and mixture plantation stands in Maoershan Experiment Station, Northeast China.
Root and soil samples were collected from the midpoint between selected trees and the neighbors in monocultures and between selected trees and the neighboring trees of another species in mixtures. First-order root samplings and rhizosphere soil samplings were taken using the methods described previously (Phillips and Fahey, 2006; Chen et al., 2018). Before sampling, litter on the soil surface around the trees was removed and then the first-order roots were carefully traced along the main root of the tree by a hand shovel. At least 100 first-order roots were taken from each sample tree. After gently shaking the fine roots, the soil loosely and tightly attached to the root system was taken by handshaking and brushing, respectively, and was defined as rhizosphere soil. Rhizosphere soil samples were taken from each sample tree within a 0–10 cm soil depth. Bulk soil, i.e., soil not adhering to root, was also sampled within the same soil depth. Nine root samples, nine rhizosphere soil samples and nine bulk soil samples of each tree species in each plot (3 replicated plots) were pooled to produce plot a level sample, respectively (Han et al., 2020; Luo et al., 2020). The root and soil samples were cold stored in a cooler and taken to the lab in no time.
Root samples were rinsed by deionized water, the living roots were selected and fixed in formalin-aceto-alcohol (FAA, 90 ml of 70% alcohol, 5 ml of 38% formaldehyde and 5 ml of glacial acetic acid) solution for subsequent determination of mycorrhizal colonization rate (Yang et al., 2021b). Soil samples were sieved with a 2-mm mesh to remove roots and litter. Rhizosphere soil samples were put in a −80°C refrigerator for later DNA detection, and fresh bulk soil of each sample was divided into two parts, one was stored at −20°C in a refrigerator for nitrate N (NO3−-N), ammonium N (NH4+-N), and dissolved organic N (DON) assaying, and the other was air-dried and ground for the other variables’ determination.
Soil chemical analyses
Soil moisture content (SM) was measured using the drying and weighing method. In brief, SM was measured by weighting before and after drying approximate 25 g fresh soils at 105°C overnight. Soil pH was measured in a soil-water suspension at a soil:water ratio of 1:2.5 (Zarif et al., 2020). The contents of nitrate N (NO3−-N) and ammonium N (NH4+-N) were extracted by 1 M KCl extraction and measured using an AA3 continuous flow analyzer (SEAL AA3, Norderstedt, Germany; Zarif et al., 2020). Litter C/N ratio and soil total N (TN) content were determined by an Elemental Analyzer (Elementar Analysensysteme GmbH, Hanau, Germany). Available P (AP) was extracted by 0.5 M NaHCO3 and determined using the molybdenum blue method (Watanabe and Olsen, 1965). The contents of dissolved organic N (DON) and organic material (OM) were determined using the methods described by Lu (1999).
Mycorrhizal colonization rate
First-order root tips of AM and EM tree species were randomly selected from the FAA solution in each plot and stained using the previous method (Phillips, 1970). In brief, root tips were softened using 5% KOH and then heated in water at a constant temperature of 65°C for 5 min. The root tips were cleaned with 1% HCl for 10 min and then stained with 0.05% trypan blue. The AMF colonization rate was estimated according to the gridline intercept method as described by Giovannetti and Mosse (1980). The roots of EMF tree species were observed under the microscope to calculate the EMF mycorrhizal colonization rate, which was calculated by the number of colonized root tips/total number of root tips (Cline et al., 2005; Suz et al., 2014; Nagati et al., 2019).
DNA extraction and PCR amplification
All genomic DNA samples were distilled using OMEGA Soil DNA Kit (M5635-02; Omega Bio-Tek, Norcross, GA, United States) following the manufacturer’s instructions, and placed at –20°C for the next step of assay. The number and quality of DNA obtained in the previous step were measured using a NanoDrop NC2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States) and agarose gel electrophoresis, respectively.
Nested PCR was conducted to amplify the AMF gene fragments with high specific amplification. In the first round of amplification, the AML1F (5′-ATCAACTTTCGATGGTAGGATAGA-3′) and AML2R (5′-GAACCCAAACACTTTGGTTTCCTTGGTTTCC-3′) were used as the primers, whereas in the second amplification step, the AMV4.5NF (5′-AAGCTCGTAGTTGAATTTCG-3′) and AMDGR (5′-CCCAACTATCCCTATTAATCAT-3′) (Luo et al., 2020) were used as the primers. For the amplification of EMF, the highly variable ITS1 region of the fungal ITS gene was selected as the target fragment, and fungal ITS5F (5′-GGAAGTAAAAGTCAACAAGG-3′) and ITS1R (5′-GCTGCGTTCTTCATCGATGC-3′) were selected as specific primers (Ma et al., 2022).
PCR amplification: firstly, denature the mold at 98°C for 5 min, followed by denaturation at 98°C for 30 s, then 53°C for 30 s, 72°C for 45 s, allowing the primers to extend on the template and synthesize DNA. The cycle was maintained for 25 times, so that a large amount of amplified DNA fragments accumulated. Final extension was at 72°C for 5 min to make the primer extension complete (Zong et al., 2022). Vazyme V AHTSTM DNA Clean Beads (Vazyme, Nanjing, China) and the Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, United States) were used to purify and quantify PCR amplicons. When individual quantification step finished, amplicons were pooled in equal amounts, and pair-end 2 × 250 bp sequencing was performed using the Illumina MiSeq platform with MiSeq Reagent Kit v3 at Shanghai Personal Biotechnology Co., Ltd. (Shanghai, China; Zong et al., 2022).
Sequence analysis
Following the official course, with minor modifications, microbiome bioinformatics was carried out by QIIME2 2019.41 (Bolyen et al., 2019). In short, the original sequence data were dissected repeatedly using the demux plugin, and then the primers were cut using the cut adapt plugin (Martin, 2011). Then, Sequences were merged, quality filtered and de-replicated with fastq_mergepairs, fastq_filter and derep_fullength in Vsearch plugen. All the unique sequences were clustered at 98% (via cluster_size), and then chimera was removed (via uchime_denovo). Finally, the non_chimera sequences were clustered at 97% to generate OTU representative sequences and OUT Table. Non-singleton amplicon sequence variants (OTUs) were aligned using mafft (Katoh, 2002). After the completion of sequence clustering, RDP classifier algorithm was used to classify the taxonomy of each AMF and EMF OTU from the domain level to the OTU level by comparing with MaarjAM database (Maarjam 081; ÖPik et al., 2010) and Unite v8 database (Nilsson et al., 2019), respectively, with the default confidence threshold of 0.7, and OTUs of EMF were distinguished against the FUNGuild database (Nguyen et al., 2016). Sequence data analyses were mainly performed using QIIME2 and R (v3.2.0) packages. Dilution curves were generated to compare the effects of sequencing depth on the diversity of observation samples (QIIME2). The OTU-level diversity, including number of OTUs observed (Sobs), Chao1, Shannon diversity index, Pielou index and Simpson index, were calculated using the OUT table in QIIME2.
Statistical analyses
The differences in AMF community characteristics, soil chemical properties and litter quality associated with AMF trees across treatments were examined with one-way analysis of variance (one-way ANOVA) or Kruskal-Wallis test. The normality and variance homogeneity were examined by using the Shapiro–Wilk normality test and Levene’s test, respectively. When the assumption met, one-way ANOVA followed by Turkey’s HSD for pair-wise comparisons was used. Otherwise, the Kruskal-Wallis test was used. The differences in EMF community characteristics soil chemical properties and litter quality associated with EMF trees across treatments were examined with Two-way analysis of variance (two-way ANOVA) or Scheirer-Ray-Hare test. Spearman correlations between colonization rates and environmental factors as well as those between alpha (α) diversity indices and environmental factors were investigated using the “pheatmap” package. Venn diagrams were generated to compare the shared and unique OTUs between different stands with the R package, based on the occurrence of OTUs in a sample or group regardless of their relative abundance using VennDiagram. Beta (β) diversity were measured by calculated the Bray-Curtis dissimilarities. In addition, the Bray–Curtis dissimilarity metrics were fit on the nonmetric multidimensional scaling (NMDS) graph using the “vegan” package to show the differences across treatments. Statistical significance was determined by permutational multivariate analysis of variance (PERMANOVA) on Bray–Curtis dissimilarities with 999 permutations to study the effect of different factors on the structure of fungal communities using the “vegan” package. Finally, the major environmental factors shaping soil mycorrhizal community were selected with ‘EnvFit’ function in the ‘vegan’ package of R. All statistical analysis was performed in R (version 4.2.0).
Results
Soil physicochemical properties and litter quality
In AM tree stands, soil nitrate nitrogen content, DON, OM, soil C/N and litter C/N were significant changed by mixture identity (Table 2). Specifically, the M-AS significantly reduced soil nitrate nitrogen content, while significantly increased OM content, soil C/N and litter C/N in comparison with the P-A monoculture. DON content was significantly reduced in mixture stands, and its reduction in the M-AS stands was more pronounced than that in the M-AL stands. AP in the M-AS stands was significantly higher than that in the M-AL stands. In EM tree stands, soil pH and OM content were significantly affected by planting pattern, tree species and their interaction, soil nitrate nitrogen content, DON content, soil C/N and litter C/N were significantly affected by planting pattern and tree species, AP was significantly affected by the interaction of planting pattern and tree species (Supplementary Table S1).
Table 2. Soil physicochemical properties and litter quality of ash monoculture, Ash-Larch mixture and Ash-Spruce mixture stands at Maoershan Experiment Station, Northeast China.
Colonization rate of mycorrhizal fungi
AMF colonization rate in roots of forest plantation stands ranging from 55.86 to 64.81%, and no significant difference in the AMF colonization rate across examined stands (Figure 1A, Kruskal-Wallis test, χ2 = 2.489, df = 2, p = 0.34). The EMF colonization rate of forest plantation stands ranging from 40 to 57%. Both mixture (F = 50.417, df = 1, p = 0.0001) and tree species (F = 6.017, df = 1, p = 0.0398) exerted significant effects of colonization rate of EMF community. Specifically, the colonization rate of EMF community in mixture was significantly lower than that of monoculture for the same tree species, and the colonization rate of EMF community in spruce stand was significantly higher than that of larch stand under identical plantation pattern. Further analysis showed that EMF colonization rate was significantly positively correlated to OM content and litter C/N, whereas negatively correlated to soil pH, nitrate nitrogen content and DON content (Table 3; Supplementary Table S1).
Figure 1. Colonization rate of arbuscular mycorrhizal fungi (A) and ectomycorrhizal fungi (B) in stands of monocultures and mixtures at the Maoershan Experiment Station, Northeast China. M-AL, Ash-Larch mixture; M-AS, Ash-Spruce mixture; P-L, Larch monoculture; P-S, Spruce monoculture. Mon, monoculture; Mix, mixture. All values represent mean ± standard error (SE), with three replicates per treatment. Different lowercase letters indicate significant differences between different stands according to Tukey’s honestly significant difference test (α = 0.05).
Diversity of mycorrhizal fungi
For AMF community, both M-AS and M-AL mixture significantly increased Pielou index and Simpson index, whereas only M-AS significantly increased Sobs (Table 4). All diversity indices significantly positively correlated to litter C/N, while significantly negatively correlated to DON and soil nitrate nitrogen content (Supplementary Figure S2A). For EMF community, mixture significantly affected examined diversity indices except for Chao1, regardless of tree species (Figure 2; Supplementary Table S2). The Sobs significantly positively correlated to litter C/N, while significantly negatively correlated to soil pH. The Pielou index significantly positively correlated to litter C/N. The Shannon index significantly positively correlated to litter C/N and soil OM content, whereas significantly negatively correlated to soil pH and soil nitrate nitrogen content. The Simpson index significantly positively correlated to litter C/N and soil C/N, whereas significantly negatively correlated to soil nitrate nitrogen content, DON content and pH (Supplementary Figure S2B).
Figure 2. The relative abundance of genera of arbuscular mycorrhizal fungi (A) and ectomycorrhizal fungi (B) communities, venn diagram displaying the shared OTUs in arbuscular mycorrhizal fungi (C) and ectomycorrhizal fungi (D) communities as well as NMDS analysis of arbuscular mycorrhizal fungi (E), ectomycorrhizal fungi (F) communities in rhizosphere soils of plantation stands at the Maoershan Experiment Station, Northeast China. P-A, Ash monoculture; M-AL, Ash-Larch mixture; M-AS, Ash-Spruce mixture; P-L, Larch monoculture; P-S, Spruce monoculture; pat, pattern; spe, species.
Community structure and composition
By clustering OTU at a 97% similarity threshold, we identified 770 AMF OTUs and 146 EMF OTUs. The representative AMF and EMF OTUs were associated with 1 phylum, 3 classes, 4 orders, 8 families, 8 genera and 2 phyla, 4 classes, 11 orders, 12 families, 12 genera, respectively. As indicated by the gentle dilution curve, the sequencing depth was the sequencing depth was adequate for evaluating the structure and diversity of mycorrhizal fungi across all samples (Supplementary Figure S3). For AMF community in the rhizosphere soils, the dominant genus of all plantation stands was the Glomus (Figure 2A). Its relative abundances were 90.89 ± 8.58 (%), 60.89 ± 19.2 (%) and 59.84 ± 0.86 (%), respectively. In addition, the relative abundance of the Paraglomus was relatively high in some rhizosphere soil samples of the M-AL and M-AS mixture stands. For EMF community, the dominant genera and their relative abundances varied with plantation pattern and tree species. In rhizosphere soils of all examined plantation stands, the Tomentella was the dominant genus (Figure 2B). Its relative abundances in rhizosphere soils of the M-AL and M-AS stands were 83.41 ± 8.78 (%) and 81.65 ± 4.63 (%), respectively (Figure 2B). In contrast, the relative abundances of the Tomentella genus in rhizosphere soils of P-L and P-S were 57.19 ± 4.62 (%) and 51.72 ± 13.29 (%), respectively. As shown in the Veen plots, 19 AMF OTUs were shared by the monoculture and mixture stands, and the number of unique OTU in mono stands was also much lower than that in the mixed stands (Figure 2C). A similar result was observed for EMF. Specifically, the number of EMF OTUs shared by M-AL and P-L, M-AS and P-S stands were 23 and 9, respectively, while there were only 2 shared EMF OTUs in four stands (Figure 2D). The NMDS displayed that AMF community composition in the rhizosphere soil under AM monoculture and AM-EM mixtures differed significantly (Figure 2E), the same is true for EMF community composition (Figure 2F). Further analysis showed that the determinant factors of EMF community are soil moisture, pH, nitrate nitrogen content, DON content, OM content, soil C/N and litter C/N (Supplementary Table S3).
Discussion
Since AM and EM tree species usually coexist (Kubisch et al., 2016), this study investigated the effects of mixed conifer-broadleaf trees on mycorrhizal community in rhizosphere soils of AM and EM tree species, which differ in litter quality and N utilizing strategy in Northeast, China by the Illumina MiSeq sequencing. Our results showed that mixed conifer-broadleaf differentially affected the soil physicochemical properties and the diversity, community structure and composition of AMF and EMF.
Litter quality and soil physicochemical properties
As hypothesized, mixture of conifer-broadleaf trees changed soil physicochemical properties and litter quality (Tables 2, 3). Previous studies suggested that difference in plant species composition and litter diversity contributes to changes in patterns and rates of litter decomposition by microbial community in forests (ÖPik et al., 2010). In general, litter C/N ratio of EM trees is significantly higher than that of AM trees (Phillips et al., 2013), and thus difficult for microorganisms to utilize and decompose (Zeng et al., 2014), resulting in a low contribution to soil N availability content (Phillips et al., 2013; Lin et al., 2017; Witte et al., 2017). Moreover, Cremer and Prietzel (2017) found that mixed conifer-broadleaf tree litter had a significant effect on soil pH, which depends on organic acid production resulting from litter decomposition. In agreement with these findings, the litter C/N ratio of the mixtures was higher than that of EM monocultures and lower than that of AM monoculture, whereas the opposite held true for soil N availability in the present study (Tables 2, 3). In combination with Li et al. (2021b), we observed mixture significantly reduced pH but only in spruce monoculture stands (Table 3), indicating the effect of mixture on soil pH depends on mixture identity and tree species.
Table 3. Soil physicochemical properties and litter quality of Larch monoculture, Spruce monoculture, Ash-Larch mixture and Ash-Spruce mixture stands at Maoershan Experiment Station, Northeast China.
Colonization rate of mycorrhizal fungi
Previous studies suggested that differences in soil N content will lead to significant changes in mycorrhizal colonization (Blanke et al., 2005), and increased soil N content are expected to reduce AMF colonization rate (e.g., Johnson et al., 2003; Johnson, 2010; Frater et al., 2018) because host plants will allocate less carbon to AMF when belowground resources no longer constrain plants (Johnson, 2010; Chen et al., 2017). Furthermore, AMF abundance was negatively correlated with soil nitrate nitrogen content (Cui et al., 2016). Similarly, enhanced N availability is also detrimental for EMF colonization (ÖPik et al., 2010; Witte et al., 2017; Yang et al., 2021a). In the present study, we found no significant effect of mixture on AMF colonization rate (Figure 1A). Additionally, we found mixed conifer-broadleaf trees significantly changed soil nitrate nitrogen content and DON content (Table 2), whereas no significant association between the AMF colonization rate and environmental factors including litter C/N and soil physicochemical properties were found in our study (Supplementary Figure S2A). This implies there maybe another important factors influencing AMF colonization, further study should be conducted. However, we found EMF colonization was significantly lower in mixtures than in monoculture for the same tree species and negatively correlated with soil nitrate nitrogen content and DON content in this study (Figure 1B; Supplementary Figure S2B). With an increase in soil N availability, host plants obtain enough nitrogen with less carbon input, resulting in a decrease in EMF colonization (Sun et al., 2010).
Diversity of mycorrhizal fungi
In this study, we found mixture did not affect the Chao1 for AMF and EMF communities, whereas the diversity change magnitude depends on the target index and mixture identity (Table 4; Supplementary Table S2). Our findings are contrast to the reports that addressing the mixed conifer-broadleaf trees on soil fungal diversity and richness (Wu et al., 2019; Ji et al., 2021; Li et al., 2021b). Wu et al. (2019) found that diversity indices were higher in soil samples from mixed plantations than in those from corresponding pure forests. Ji et al. (2021) found that the diversity and richness in soil samples from mixed plantations were comparable to those from corresponding pure forests. Li et al. (2021b) reported that the Chao1 index of EMF was responsive to mixture. The discrepancy in the effects of mixture on mycorrhizal fungal community richness and diversity across studies maybe correlated to tree species, mycorrhizal type, soil nutrient, soil fauna and plant–soil organism interactions. In agreement with previous studies that EMF richness and diversity are negatively correlated to soil nutrient content (Corrales et al., 2016; Erlandson et al., 2016), particularly the N availability (Lilleskov et al., 2011; Hasselquist and Högberg, 2014; Yang et al., 2021a; Boeraeve et al., 2022).
Table 4. AMF alpha diversity indices of Ash monoculture, Ash-Larch mixture and Ash-Spruce mixture stands at Maoershan Experiment Station, Northeast China.
Community structure and composition of mycorrhizal fungi
In agreement with our prediction that AM and EM fungal community composition in the mixture would shift due to changes in soil physicochemical properties and litter quality. Previous studies demonstrate that the identities of mixed species are determinant of AMF community in rhizosphere soil (Mummey and Rillig, 2006; Turrini et al., 2016; Ji et al., 2021) and changes in soil properties result in shifts in mycorrhizal fungal communities (Avis et al., 2003; Roy et al., 2017; Chen et al., 2018; Lu et al., 2020; Gu et al., 2022). We found strikingly differences in the AMF communities between M-AL and P-A stands, as well as in the EMF communities between M-AS and P-S stands (Figure 2).
In line with previous studies (Wang et al., 2018; Zeng et al., 2021; Wang and Han, 2022), the Glomus and the Tomentella were the most abundant genera in AMF and EMF communities, respectively. The Glomus was primarily distributed in the upper soil layer, particularly where the available N content was high (Wang et al., 2018; Chen et al., 2021). We found the relative abundance of Glomus was higher in AM monoculture, and the relative abundance was negatively correlated to nitrate nitrogen. What accounts for this discrepancy is unknown but warrants further studying. The Tomentella has a strong adaptability to the environment and strong colonization ability (Nara, 2006; ÖPik et al., 2010). Meanwhile, more abundance and species of the Tomentella was found in soils with high N than that with low N (Kranabetter et al., 2015; Wang and Han, 2022). Besides, EM root tips colonized by Tomentella are capable of producing extracellular enzymes used to break down proteins, polysaccharides, and organic forms of P (Courty et al., 2005; Tedersoo et al., 2012) and mobilizing organic N (Pena et al., 2013). We found that mixture enhanced the relative abundance of Tomentella (Figure 2B), indicating it would increase the nutrient availability.
Finally, we found that the determinant factors of EMF community are soil moisture, pH, nitrate nitrogen content, dissolved organic nitrogen content, soil organic matter content, soil organic carbon/total nitrogen and litter carbon/total nitrogen. In earlier studies, soil moisture (Erlandson et al., 2016), soil organic carbon (Yang et al., 2019; Li et al., 2021a), nitrate nitrogen (He et al., 2016), soil C/N (Zhang et al., 2013; Sterkenburg et al., 2015), pH (Bahnmann et al., 2018; Craig et al., 2018) and litter C/N (Hewitt et al., 2017) are found to be the critical factors shaping fungal community composition. In addition, soil pH, soil organic matter, total nitrogen, C/N ratio, and total phosphorus influenced the EMF community structure (Li et al., 2021b). In general, our findings again highlight the linkage between soil physicochemical properties and EMF community composition.
The major limitation of the present studies is that we did not consider the seasonal variation. In addition, we only compared the AMF and EMF with a few tree species with limited samples. To draw a general conclusion about responses of mycorrhizal community to conifer-broadleaf mixture, future studies should be conducted under different conditions, such as host plants, climates, soil texture and altitudes. Finally, we did not consider the effects of changes in biological factors, including species and diversity of aboveground plants, root exudates and other soil fauna belowground.
Conclusion
In the present study, changes in soil physicochemical properties and litter quality due to mixed planting influenced the mycorrhizal colonization rate, composition, and diversity of the rhizosphere mycorrhizal community, especially the EMF community. The determinant factors of EMF community were soil moisture, pH, nitrate nitrogen content, dissolved organic nitrogen content, soil organic matter content, soil organic carbon/total nitrogen and litter carbon/total nitrogen. These results will assist in understanding and predicting variations in mycorrhizae in forest ecosystems composed of different tree species associated with different mycorrhizae. These results provide a better understanding of the ecological functions of the soil–fungal community between mixed and pure plantations, and cast a new light on underlying mechanisms between soil nutrients and soil–fungal functional community, especially for those tree species used in forest plantation.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/genbank/, PRJNA846491 https://www.ncbi.nlm.nih.gov/genbank/, PRJNA846556.
Author contributions
All authors contributed to the study conception and design. The field and laboratory work were carried out by YZ, LX, SM, DC, and KZ. The first draft of the manuscript was written by YZ and revised by QW and WF. All authors commented on previous versions of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Fundamental Research Funds for the Central Universities of China (2572017PZ03, 2572020DR04, and 2572019CP16).
Acknowledgments
We acknowledge Limei Yan, Jiajing Duan, Hongli Li, and Hongwei Zhang for their help in sampling and laboratory analysis. We appreciate Zarif Nowsherwan for his help in the manuscript revision process. We are also grateful to the staff of the Maoershan Experimental Station and the Shanghai Personal Biotechnology Co., Ltd. for their support.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.986515/full#supplementary-material
Footnotes
1. ^https://docs.qiime2.org/2019.4/tutorials/ (accessed September 13, 2022).
References
Albornoz, F. E., Lambers, H., Turner, B. L., Teste, F. P., and Laliberté, E. (2016). Shifts in symbiotic associations in plants capable of forming multiple root symbioses across a long-term soil chronosequence. Ecol. Evol. 6, 2368–2377. doi: 10.1002/ece3.2000
Aldrich-Wolfe, L., Black, K. L., Hartmann, E. D. L., Shivega, W. G., Schmaltz, L. C., McGlynn, R. D., et al. (2020). Taxonomic shifts in arbuscular mycorrhizal fungal communities with shade and soil nitrogen across conventionally managed and organic coffee agroecosystems. Mycorrhiza 30, 513–527. doi: 10.1007/s00572-020-00967-7
Averill, C., Fortunel, C., Maynard, D. S., Van den Hoogen, J., Dietze, M. C., Bhatnagar, J. M., et al. (2022). Alternative stable states of the forest mycobiome are maintained through positive feedbacks. Nat. Ecol. Evol. 6, 375–382. doi: 10.1038/s41559-022-01663-9
Avis, P. G., Mclaughlin, D. J., Dentinger, B. C., and Reich, P. B. (2003). Long-term increase in nitrogen supply alters above-and below-ground ectomycorrhizal communities and increases the dominance of Russula spp. in a temperate oak savanna. New Phytol. 160, 239–253. doi: 10.1046/J.1469-8137.2003.00865.X
Bahnmann, B., Mašínová, T., Halvorsen, R., Davey, M. L., Sedlak, P., Tomšovský, M., et al. (2018). Effects of oak, beech and spruce on the distribution and community structure of fungi in litter and soils across a temperate forest. Soil Biol. Biochem. 119, 162–173. doi: 10.1016/j.soilbio.2018.01.021
Blanke, V., Renker, C., Wagner, M., Füllner, K., Held, M., Kuhn, A. J., et al. (2005). Nitrogen supply affects arbuscular mycorrhizal colonization of Artemisia vulgaris in a phosphate-polluted field site. New Phytol. 166, 981–992. doi: 10.1111/j.1469-8137.2005.01374.x
Boeraeve, M., Kohout, P., Ceulemans, T., Cajthaml, T., Tedersoo, L., and Jacquemyn, H. (2022). Changes in the root microbiome of four plant species with different mycorrhizal types across a nitrogen deposition gradient in ombrotrophic bogs. Soil Biol. Biochem. 169:108673. doi: 10.1016/j.soilbio.2022.108673
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9
Brundrett, M. C., and Tedersoo, L. (2018). Evolutionary history of mycorrhizal symbioses and global host plant diversity. New Phytol. 220, 1108–1115. doi: 10.1111/nph.14976
Chen, N., Gregory, M., Louise, E. W., Andrew, W., Yan, W. D., and Xiang, W. H. (2018). Diversity and enzyme activity of ectomycorrhizal fungal communities following nitrogen fertilization in an urban-adjacent pine plantation. Forests 9:99. doi: 10.3390/f9030099
Chen, L., Li, W., and Xiao, Y. (2021). Biochar and nitrogen fertilizer increase glomus synergism and abundance and promote trifolium pratense growth while inhibiting pollutant accumulation. Ecol. Indic. 133:108377. doi: 10.1016/j.ecolind.2021.108377
Chen, Y. L., Xu, Z. W., Xu, T. L., Veresoglou, S. D., Yang, G. W., and Chen, B. D. (2017). Nitrogen deposition and precipitation induced phylogenetic clustering of arbuscular mycorrhizal fungal communities. Soil Biol. Biochem. 115, 233–242. doi: 10.1016/j.soilbio.2017.08.024
Cline, E. T., Ammirati, J. F., and Edmonds, R. L. (2005). Does proximity to mature trees influence ectomycorrhizal fungus communities of Douglas-fir seedlings? New Phytol. 166, 993–1009. doi: 10.1111/j.1469-8137.2005.01387.x
Corrales, A., Arnold, A. E., Ferrer, A., Turner, B. L., and Dalling, J. W. (2016). Variation in ectomycorrhizal fungal communities associated with Oreomunnea mexicana (Juglandaceae) in a Neotropical montane forest. Mycorrhiza 26, 1–17. doi: 10.1007/s00572-015-0641-8
Courty, P. E., Pritsch, K., Schloter, M., Hartmann, A., and Garbaye, J. (2005). Activity profiling of ectomycorrhiza communities in two forest soils using multiple enzymatic tests. New Phytol. 167, 309–319. doi: 10.1111/j.1469-8137.2005.01401.x
Craig, M. E., Turner, B. L., Liang, C., Clay, K., Johnson, D. J., and Phillips, R. P. (2018). Tree mycorrhizal type predicts within-site variability in the storage and distribution of soil organic matter. Glob. Change Biol. 24, 3317–3330. doi: 10.1111/gcb.14132
Cremer, M., and Prietzel, J. (2017). Soil acidity and exchangeable base cation stocks under pure and mixed stands of European beech, Douglas fir and Norway spruce. Plant Soil 415, 393–405. doi: 10.1007/s11104-017-3177-1
Cui, X., Hu, J., Wang, J., Yang, J., and Lin, X. (2016). Reclamation negatively influences arbuscular mycorrhizal fungal community structure and diversity in coastal saline-alkaline land in Eastern China as revealed by Illumina sequencing. Appl. Soil Ecol. 98, 140–149. doi: 10.1016/j.apsoil.2015.10.008
Erlandson, S. R., Savage, J. A., Cavender-Bares, J. M., and Peay, K. G. (2016). Soil moisture and chemistry influence diversity of ectomycorrhizal fungal communities associating with willow along an hydrologic gradient. FEMS Microbiol. Ecol. 92:fiv148. doi: 10.1093/femsec/fiv148
Felton, A., Nilsson, U., Sonesson, J., Felton, A. M., Roberge, J.-M., Ranius, T., et al. (2016). Replacing monocultures with mixed-species stands: ecosystem service implications of two production forest alternatives in Sweden. Ambio 45, 124–139. doi: 10.1007/s13280-015-0749-2
Forrester, D. I. (2014). The spatial and temporal dynamics of species interactions in mixed-species forests: from pattern to process. Forest Ecol. Manag. 312, 282–292. doi: 10.1016/j.foreco.2013.10.003
Frater, P. N., Borer, E. T., Fay, P. A., Jin, V., Knaeble, B., Seabloom, E., et al. (2018). Nutrients and environment influence arbuscular mycorrhizal colonization both independently and interactively in Schizachyrium scoparium. Plant Soil 425, 493–506. doi: 10.1007/s11104-018-3597-6
Gessner, M. O., Swan, C. M., Dang, C. K., Mckie, B. G., Bardgett, R. D., Wall, D. H., et al. (2010). Diversity meets decomposition. Trends Ecol. Evol. 25, 372–380. doi: 10.1016/j.tree.2010.01.010
Giovannetti, M., and Mosse, B. (1980). An evaluation of techniques for measuring vesicular arbuscular mycorrhizal infection in roots. New Phytol. 84, 489–500. doi: 10.1111/j.1469-8137.1980.tb04556.x
Gosling, P., Mead, A., Proctor, M., Hammond, J. P., and Bending, G. D. (2013). Contrasting arbuscular mycorrhizal communities colonizing different host plants show a similar response to a soil phosphorus concentration gradient. New Phytol. 198, 546–556. doi: 10.1111/nph.12169
Gu, T., Mao, Y., Chen, C., Wang, Y., Lu, Q., Wang, H., et al. (2022). Diversity of arbuscular mycorrhiza fungi in rhizosphere soil and roots in Vetiveria zizanioides plantation chronosequence in coal gangue heaps. Symbiosis 86, 111–122. doi: 10.1007/s13199-022-00829-0
Han, M., Sun, L., Gan, D., Fu, L., and Zhu, B. (2020). Root functional traits are key determinants of the rhizosphere effect on soil organic matter decomposition across 14 temperate hardwood species. Soil Biol. Biochem. 151:108019. doi: 10.1016/j.soilbio.2020.108019
Hasselquist, N. J., and Högberg, P. (2014). Dosage and duration effects of nitrogen additions on ectomycorrhizal sporocarp production and functioning: an example from two N-limited boreal forests. Ecol. Evol. 4, 3015–3026. doi: 10.1002/ece3.1145
He, F., Yang, B., Wang, H., Yan, Q., Cao, Y., and He, X. (2016). Changes in composition and diversity of fungal communities along Quercus mongolica forests developments in Northeast China. Appl. Soil Ecol. 100, 162–171. doi: 10.1016/j.apsoil.2015.12.014
Herrmann, L., Lesueur, D., Bräu, L., Davison, J., Jairus, T., Robain, H., et al. (2016). Diversity of root-associated arbuscular mycorrhizal fungal communities in a rubber tree plantation chronosequence in Northeast Thailand. Mycorrhiza 26, 863–877. doi: 10.1007/s00572-016-0720-5
Hewitt, R. E., Chapin, F. S. III, Hollingsworth, T. N., and Taylor, D. L. (2017). The potential for mycobiont sharing between shrubs and seedlings to facilitate tree establishment after wildfire at Alaska arctic treeline. Mol. Ecol. 26, 3826–3838. doi: 10.1111/mec.14143
Högberg, M. N., Högberg, P., Wallander, H., and Nilsson, L. O. (2021). Carbon-nitrogen relations of ectomycorrhizal mycelium across a natural nitrogen supply gradient in boreal forest. New Phytol. 232, 1839–1848. doi: 10.1111/nph.17701
Ji, L., Yang, Y. C., Yang, N., Khan, A., and Yang, L. X. (2021). Seasonal variation of diversity and co-occurrence patterns of arbuscular mycorrhizal fungal communities in mixed broadleaf-conifer forests. Appl. Soil Ecol. 158:103782. doi: 10.1016/j.apsoil.2020.103782
Jiang, F., Zhu, K., Cadotte, M. W., and Jin, G. (2020). Tree mycorrhizal type mediates the strength of negative density dependence in temperate forests. J. Ecol. 108, 2601–2610. doi: 10.1111/1365-2745.13413
Johnson, N. C. (2010). Resource stoichiometry elucidates the structure and function of arbuscular mycorrhizas across scales. New Phytol. 185, 631–647. doi: 10.1111/j.1469-8137.2009.03110.x
Johnson, N. C., Rowland, D. L., Corkidi, L., Egerton-Warburton, L. M., and Allen, E. B. (2003). Nitrogen enrichment alters mycorrhizal allocation at five Mesic to semiarid grasslands. Ecology 84, 1895–1908. doi: 10.1890/0012-9658(2003)084[1895:NEAMAA]2.0.CO;2
Katoh, K. (2002). Mafft: a novel method for rapid multiple sequence alignment based on fast fourier transform. Nucleic Acids Res. 30, 3059–3066. doi: 10.1093/nar/gkf436
Kilpeläinen, J., Barbero-López, A., Adamczyk, B., Aphalo, P. J., and Lehto, T. (2019). Morphological and ecophysiological root and leaf traits in ectomycorrhizal, arbuscular-mycorrhizal and non-mycorrhizal alnus incana seedlings. Plant Soil 436, 283–297. doi: 10.1007/s11104-018-03922-w
Kranabetter, J. M., Hawkins, B. J., Jones, M. D., Robbins, S., Dyer, T., and Li, T. (2015). Species turnover (β-diversity) in ectomycorrhizal fungi linked to uptake capacity. Mol. Ecol. 24, 5992–6005. doi: 10.1111/mec.13435
Kubisch, P., Hertel, D., and Leuschner, C. (2016). Fine root productivity and turnover of ectomycorrhizal and arbuscular mycorrhizal tree species in a temperate broad-leaved mixed forest. Front. Plant Sci. 7:1233. doi: 10.3389/fpls.2016.01233
Lang, A. C., von Oheimb, G., Scherer-Lorenzen, M., Yang, B., Trogisch, S., Bruelheide, H., et al. (2014). Mixed afforestation of young subtropical trees promotes nitrogen acquisition and retention. J. Appl. Ecol. 51, 224–233. doi: 10.1111/1365-2664.12157
Leigh, J., Hodge, A., and Fitter, A. H. (2009). Arbuscular mycorrhizal fungi can transfer substantial amounts of nitrogen to their host plant from organic material. New Phytol. 181, 199–207. doi: 10.1111/j.1469-8137.2008.02630.x
Li, W. Q., Huang, Y. X., Chen, F. S., Liu, Y. Q., Lin, X. F., Zong, Y. Y., et al. (2021a). Mixing with broad-leaved trees shapes the rhizosphere soil fungal communities of coniferous tree species in subtropical forests. Forest Ecol. Manag. 480:118664. doi: 10.1016/j.foreco.2020.118664
Li, X., Kang, W., Liu, S., Yin, H., Lyu, Q., Su, Y., et al. (2021b). Diversity of ectomycorrhizal fungal communities in four types of stands in Pinus massoniana plantation in the west of China. Forests 12:719. doi: 10.3390/f12060719
Lilleskov, E. A., and Bruns, T. D. (2003). Root colonization dynamics of two ectomycorrhizal fungi of contrasting life history strategies are mediated by addition of organic nutrient patches. New Phytol. 159, 141–151. doi: 10.1046/j.1469-8137.2003.00794.x
Lilleskov, E. A., Hobbie, E. A., and Horton, T. R. (2011). Conservation of ectomycorrhizal fungi: exploring the linkages between functional and taxonomic responses to anthropogenic N deposition. Fungal Ecol. 4, 174–183. doi: 10.1016/j.funeco.2010.09.008
Lin, G. G., McCormack, M. L., Ma, C. G., and Guo, D. L. (2017). Similar below-ground carbon cycling dynamics but contrasting modes of nitrogen cycling between arbuscular mycorrhizal and ectomycorrhizal forests. New Phytol. 213, 1440–1451. doi: 10.1111/nph.14206
Lu, R. K. (1999). Analytical Methods of Agricultural Chemistry in Soil. Beijing: China Agricultural Scinetech Press.
Lu, Y., Liu, X., Chen, F., and Zhou, S. R. (2020). Shifts in plant community composition weaken the negative effect of nitrogen addition on community-level arbuscular mycorrhizal fungi colonization. P. Roy. Soc. B-Biol. Sci. 287:20200483. doi: 10.1098/rspb.2020.0483
Luo, Y., Wang, Z., He, Y., Li, G., and Zhuang, L. (2020). High-throughput sequencing analysis of the rhizosphere arbuscular mycorrhizal fungi (AMF) community composition associated with ferula sinkiangensis. BMC Microbiol. 20:335. doi: 10.1186/s12866-020-02024-x
Ma, L. W., Liu, L., Lu, Y. S., Lin, C., Zhang, Z. C., Zhang, H. W., et al. (2022). When microclimates meet soil microbes: temperature controls soil microbial diversity along an elevational gradient in subtropical forests. Soil Biol. Biochem. 166:108566. doi: 10.1016/j.soilbio.2022.108566
Ma, X., Zhu, B., Nie, Y., Liu, Y., and Kuzyakov, Y. (2021). Root and mycorrhizal strategies for nutrient acquisition in forests under nitrogen deposition: a meta-analysis. Soil Biol. Biochem. 163:108418. doi: 10.1016/j.soilbio.2021.108418
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. Embnet J. 17, 10–12. doi: 10.14806/ej.17.1.200
Mummey, D. L., and Rillig, M. C. (2006). The invasive plant species Centaurea maculosa alters arbuscular mycorrhizal fungal communities in the field. Plant Soil 288, 81–90. doi: 10.1007/s11104-006-9091-6
Nagati, M., Roy, M., Desrochers, A., Manzi, S., and Gardes, M. (2019). Facilitation of balsam fir by trembling aspen in the boreal forest: do ectomycorrhizal communities matter? Front. Plant Sci. 10:932. doi: 10.3389/fpls.2019.00932
Nara, K. (2006). Ectomycorrhizal networks and seedling establishment during early primary succession. New Phytol. 169, 169–178. doi: 10.1111/j.1469-8137.2005.01545.x
Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2016). FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. doi: 10.1016/j.funeco.2015.06.006
Nilsson, R. H., Larsson, K. H., Taylor, A. F. S., Bengtsson-Palme, J., Jeppesen, T. S., Schigel, D., et al. (2019). The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications. Nucleic Acids Res. 47, D259–D264. doi: 10.1093/nar/gky1022
Oehl, F., Laczko, E., Bogenrieder, A., Stahr, K., Bisch, R., van der Heijden, M., et al. (2010). Soil type and land use intensity determine the composition of arbuscular mycorrhizal fungal communities. Soil Biol. Biochem. 42, 724–738. doi: 10.1016/j.soilbio.2010.01.006
ÖPik, M., Vanatoa, A., Vanatoa, E., Moora, M., Davison, J., Kalwij, J. M., et al. (2010). The online database Maarj AM reveals global and ecosystemic distribution patterns in arbuscular mycorrhizal fungi (Glomeromycota). New Phytol. 188, 223–241. doi: 10.1111/j.1469-8137.2010.03334.x
Pena, R., Tejedor, J., Zeller, B., Dannenmann, M., and Polle, A. (2013). Interspecific temporal and spatial differences in the acquisition of litter-derived nitrogen by ectomycorrhizal fungal assemblages. New Phytol. 199, 520–528. doi: 10.1111/nph.12272
Phillips, J. M. (1970). Improved procedures for clearing roots and staining parasitic and vesicular-arbuscular mycorrhizal fungi for rapid assessment of infection. Trans. Br. Mycol. Soc. 55:158-IN18. doi: 10.1016/S0007-1536(70)80110-3
Phillips, R. P., Brzostek, E., and Midgley, M. G. (2013). The mycorrhizal-associated nutrient economy: a new framework for predicting carbon–nutrient couplings in temperate forests. New Phytol. 199, 41–51. doi: 10.1111/nph.12221
Phillips, R. P., and Fahey, T. J. (2006). Tree species and mycorrhizal associations influence the magnitude of rhizosphere effects. Ecology 87, 1302–1313. doi: 10.1890/0012-9658(2006)87[1302:TSAMAI]2.0.CO;2
Prada, C. M., Turner, B. L., and Dalling, J. W. (2022). Growth responses of ectomycorrhizal and arbuscular mycorrhizal seedlings to low soil nitrogen availability in a tropical montane forest. Funct. Ecol. 36, 107–119. doi: 10.1111/1365-2435.13901
Pretzsch, H., Forrester, D.I., and Bauhus, J. (2017). Mixed-Species Forests. Berlin, Heidelberg: Springer
Read, D. J., and Perez-Moreno, J. (2003). Mycorrhizas and nutrient cycling in ecosystems-a journey towards relevance? New Phytol. 157, 475–492. doi: 10.1046/j.1469-8137.2003.00704.x
Roy, J., Reichel, R., Brüggemann, N., and Rillig, M. (2017). Succession of arbuscular mycorrhizal fungi along a 52-year agricultural recultivation chronosequence. FEMS Microbiol. Ecol. 93:9. doi: 10.1093/femsec/fix102
Seidel, D., Leuschner, C., Scherber, C., Beyer, F., Wommelsdorf, T., Cashman, M. J., et al. (2013). The relationship between tree species richness, canopy space exploration and productivity in a temperate broad-leaf mixed forest. Forest Ecol. Manag. 310, 366–374. doi: 10.1016/j.foreco.2013.08.058
Seyfried, G. S., Dalling, J. W., and Yang, W. H. (2021). Mycorrhizal type effects on leaf litter decomposition depend on litter quality and environmental context. Biogeochemistry 155, 21–38. doi: 10.1007/s10533-021-00810-x
Sterkenburg, E., Bahr, A., Brandström Durling, M., Clemmensen, K. E., and Lindahl, B. D. (2015). Changes in fungal communities along a boreal forest soil fertility gradient. New Phytol. 207, 1145–1158. doi: 10.1111/nph.13426
Sun, Y., Gu, J. C., Zhuang, H. F., and Wang, Z. Q. (2010). Effects of ectomycorrhizal colonization and nitrogen fertilization on morphology of root tips in a Larix gmelinii plantation in northeastern China. Ecol. Res. 25, 295–302. doi: 10.1007/s11284-009-0654-x
Suz, L. M., Barsoum, N., Benham, S., Dietrich, H. P., Fetzer, K. D., Fischer, R., et al. (2014). Environmental drivers of ectomycorrhizal communities in Europe's temperate oak forests. Mol. Ecol. 23, 5628–5644. doi: 10.1111/mec.12947
Tatsumi, C., Hyodo, F., Taniguchi, T., Shi, W., Koba, K., Fukushima, K., et al. (2021). Arbuscular mycorrhizal community in roots and nitrogen uptake patterns of understory trees beneath ectomycorrhizal and non-ectomycorrhizal overstory trees. Front. Plant Sci. 11:583585. doi: 10.3389/fpls.2020.583585
Taylor, M. K., Lankau, R. A., and Wurzburger, N. (2016). Mycorrhizal associations of trees have different indirect effects on organic matter decomposition. J. Ecol. 104, 1576–1584. doi: 10.1111/1365-2745.12629
Tedersoo, L., Bahram, M., Toots, M., Diedhiou, A. G., Henkel, T. W., Kjøller, R., et al. (2012). Towards global patterns in the diversity and community structure of ectomycorrhizal fungi. Mol. Ecol. 21, 4160–4170. doi: 10.1111/j.1365-294X.2012.05602.x
Tedersoo, L., and Nara, K. (2010). General latitudinal gradient of biodiversity is reversed in ectomycorrhizal fungi. New Phytol. 185, 351–354. doi: 10.1111/j.1469-8137.2009.03134.x
Turrini, A., Sbrana, C., Avio, L., Njeru, E. M., Bocci, G., Bàrberi, P., et al. (2016). Changes in the composition of native root arbuscular mycorrhizal fungal communities during a short-term cover crop-maize succession. Biol. Fertil. Soils 52, 643–653. doi: 10.1007/s00374-016-1106-8
Van der Heijden, M. G. A., Martin, F. M., Selosse, M., and Sanders, I. R. (2015). Mycorrhizal ecology and evolution: the past, the present, and the future. New Phytol. 205, 1406–1423. doi: 10.1111/nph.13288
Wang, X., and Han, Q. (2022). A closer examination of the ‘abundant-center’ for ectomycorrhizal fungal community associated with Picea crassifolia in China. Front. Plant Sci. 13:759801. doi: 10.3389/fpls.2022.759801
Wang, H., Liu, S., Wang, J., Shi, Z., Lu, L., Zeng, J., et al. (2013). Effects of tree species mixture on soil organic carbon stocks and greenhouse gas fluxes in subtropical plantations in China. Forest Ecol. Manag. 300, 4–13. doi: 10.1016/j.foreco.2012.04.005
Wang, X., Sui, X., Liu, Y., Xiang, L., Zhang, T., Fu, J., et al. (2018). N-P fertilization did not reduce AMF abundance or diversity but alter AMF composition in an alpine grassland infested by a root hemiparasitic plant. Plant Divers. 40, 117–126. doi: 10.1016/j.pld.2018.05.001
Wardle, D. A., Nilsson, M. C., Zackrisson, O., and Gallet, C. (2003). Determinants of litter mixing effects in a Swedish boreal forest. Soil Biol. Biochem. 35, 827–835. doi: 10.1016/S0038-0717(03)00118-4
Watanabe, F. S., and Olsen, S. R. (1965). Test of an ascorbic acid method for determining phosphorus in water and NaHCO3 extracts from soil. Soil Sci. Soc. Am. J. 29, 677–678. doi: 10.2136/sssaj1965.03615995002900060025x
Witte, L. D., Rosenstock, N. P., Linde, S., and Braun, S. (2017). Nitrogen deposition changes ectomycorrhizal communities in Swiss beech forests. Sci. Total Environ. 605-606, 1083–1096. doi: 10.1016/j.scitotenv.2017.06.142
Wu, D., Zhang, M., Peng, M., Sui, X., Li, W., and Sun, G. (2019). Variations in soil functional fungal community structure associated with pure and mixed plantations in typical temperate forests of China. Front. Microbiol. 10:1636. doi: 10.3389/fmicb.2019.01636
Yang, W., Jeelani, N., Xia, L., Zhu, Z., Luo, Y., Cheng, X., et al. (2019). Soil fungal communities vary with invasion by the exotic Spartina alternifolia Loisel. In coastal salt marshes of eastern China. Plant Soil 442, 215–232. doi: 10.1007/s11104-019-04184-w
Yang, B., Pang, X. Y., Hu, B., Bao, W. K., and Tian, G. L. (2017). Does thinning-induced gap size result in altered soil microbial community in pine plantation in eastern Tibetan plateau? Ecol. Evol. 7, 2986–2993. doi: 10.1002/ece3.2714
Yang, L., Wang, P., and Kong, C. (2010). Effect of larch (Larix gmelini Rupr.) root exudates on Manchurian walnut (Juglans mandshurica maxim.) growth and soil juglone in a mixed-species plantation. Plant Soil 329, 249–258. doi: 10.1007/s11104-009-0149-0
Yang, N., Wang, B., Liu, D., Wang, X., Li, X., Zhang, Y., et al. (2021a). Long-term nitrogen deposition alters ectomycorrhizal community composition and function in a poplar plantation. J Fungi 7:791. doi: 10.3390/jof7100791
Yang, Y., Zhang, X. H., Hartley, I. P., Dungait, J., Wen, X. F., Li, D. D., et al. (2021b). Contrasting rhizosphere soil nutrient economy of plants associated with arbuscular mycorrhizal and ectomycorrhizal fungi in karst forests. Plant Soil 470, 81–93. doi: 10.1007/s11104-021-04950-9
Yang, K., Zhu, J. J., Gu, J. C., Yu, L. Z., and Wang, Z. Q. (2015). Changes in soil phosphorus fractions after 9 years of continuous nitrogen addition in a Larix gmelinii plantation. Ann. Forest Sci. 72, 435–442. doi: 10.1007/s13595-014-0444-7
Zarif, N., Khan, A., and Wang, Q. C. (2020). Linking soil acidity to P fractions and exchangeable base cations under increased N and P fertilization of mono and mixed plantations in Northeast China. Forests 11:1274. doi: 10.3390/f11121274
Zeng, Y. L., Xiang, W. H., Deng, X. W., Fang, X., Liu, C., and Peng, C. H. (2014). Soil N forms and gross transformation rates in Chinese subtropical forests dominated by different tree species. Plant Soil 384, 231–242. doi: 10.1007/s11104-014-2206-6
Zeng, H., Yu, L., Liu, P., Wang, Z., Chen, Y., and Wang, J. (2021). Nitrogen fertilization has a stronger influence than cropping pattern on AMF community in maize/soybean strip intercropping systems. Appl. Soil Ecol. 167:104034. doi: 10.1016/j.apsoil.2021.104034
Zhang, J., Taniguchi, T., Tateno, R., Xu, M., Du, S., Liu, G. B., et al. (2013). Ectomycorrhizal fungal communities of Quercus liaotungensis along local slopes in the temperate oak forests on the loess plateau. China. Ecol. Res. 28, 297–305. doi: 10.1007/s11284-012-1017-6
Keywords: arbuscular mycorrhizal fungi, ectomycorrhizal fungi, mixed plantation, colonization rate, diversity, community dissimilarity
Citation: Zhang Y, Wang Q, Xu L, Ma S, Cui D, Zhu K and Feng W (2022) Mixed conifer-broadleaf trees on arbuscular mycorrhizal and ectomycorrhizal communities in rhizosphere soil of different plantation stands in the temperate zone, Northeast China. Front. Microbiol. 13:986515. doi: 10.3389/fmicb.2022.986515
Edited by:
Rentao Liu, Ningxia University, ChinaReviewed by:
Muhammad Arif, Southwest University, ChinaWei Qiang, Chinese Academy of Sciences, China
Copyright © 2022 Zhang, Wang, Xu, Ma, Cui, Zhu and Feng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qingcheng Wang, wqcnefu@163.com