Next Article in Journal
Differential Impacts of Acacia Invasion on Nutrient Fluxes in Two Distinct Bornean Lowland Tropical Rain Forests
Next Article in Special Issue
Stand Characteristics Rather than Soil Properties Contribute More to the Expansion of Moso Bamboo (Phyllostachys edulis) into Its Neighboring Forests in Subtropical Region
Previous Article in Journal
Contribution of Arbuscular Mycorrhizal Fungal Communities to Soil Carbon Accumulation during the Development of Cunninghamia lanceolata Plantations
Previous Article in Special Issue
Species Composition and Diversity of Middle-Aged Trees among Different Urban Green Space Types and Tree Age Classes in Changchun, Northeast China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Aspect on Arbuscular Mycorrhizal Fungal Diversity and Community Composition in a Natural Toona ciliata var. pubescens Forest in Subtropical China

1
Jiangxi Provincial Key Laboratory of Silviculture, Jiangxi Agricultural University, Nanchang 330045, China
2
College of Forestry, Jiangxi Agricultural University, Nanchang 330045, China
3
Nature Protected Area Construction Center of Jiangxi, Nanchang 330036, China
4
School of Agriculture, Food and Environment, Royal Agricultural University, Cirencester GL7 6JS, UK
*
Authors to whom correspondence should be addressed.
Forests 2022, 13(12), 2100; https://doi.org/10.3390/f13122100
Submission received: 31 October 2022 / Revised: 2 December 2022 / Accepted: 7 December 2022 / Published: 9 December 2022
(This article belongs to the Special Issue Forest Biodiversity and Ecosystem Stability)

Abstract

:
The aspect can affect plant communities and soil physical and chemical properties through different microclimates. However, little is known about the effect of aspect on arbuscular mycorrhizal (AM) fungal diversity and community composition, although AM fungi are beneficial for plant nutrient absorption and natural restoration. In this study, AM fungal community and diversity distribution patterns in the rhizosphere soil and roots of seven widespread plants in a natural Toona ciliata var. pubescens (Tc) forest on the north-facing (NF) aspect and south-facing (SF) aspect were investigated using Illumina PE250 high-throughput sequencing in the Guanshan National Nature Reserve, Jiangxi Province, China. Our results exhibited that aspect did not affect AM fungal diversity but significantly affected AM fungal community structure and composition. Glomeraceae was the most common and abundant family in the Tc natural forest. The Glomeromycota sequence proportion of root AM fungal community was significantly larger on NF than on SF (p < 0.05). The relative abundance of Acaulosporaceae of root AM fungal community differed significantly with aspect, being greater on NF than on SF (p < 0.05). In addition, the number of Glomeromycota sequences was significantly larger on SF than on NF, while the number of OTUs and the relative abundance of unclassified fungi in rhizosphere soil in Tc showed the opposite trend (p < 0.05). The soil properties (organic matter, nitrogen, potassium, phosphorus, and pH) were significantly correlated with these changes. These findings indicate that the habitat of NF with low insolation, high soil moisture, and high nutrient content might promote the functional realization of AM fungi; the habitat of SF with high insolation, low soil moisture, and low soil nutrient content might be beneficial for the proliferation and preservation of AM fungal groups. This study provides important information on the ecological processes of AM fungal community construction and microbiological mechanisms in natural Tc forests.

1. Introduction

Biodiversity is an important feature of life, and species composition and community construction have long been core issues in ecological research [1,2]. Indeed, the biodiversity and structural composition and distribution of communities are closely related to ecosystem multifunctionality [2]. Topography is one of the most major natural factors that affects biodiversity. Topography belongs to biogeography and affects numerous ecological functions and processes, such as soil erosion and water storage [3], weathering and nutrient leaching [4], and migration and storage of soil organic matter [5].
Aspect (shady and sunny aspects) is a feature of natural geographical mountain landscapes [6]. Shady and sunny aspects are different habitats with climate and microenvironment differences that form a natural environmental gradient [7]. The difference in climate conditions between shady and sunny aspects leads to profound changes in vegetation, soil, landform, hydrology, and other factors [8,9]. The shady aspects are darker, wetter, and colder, with less light radiation, and the climate characteristics on sunny aspects are the opposite [10]. Changes in biodiversity distribution and community construction on shady and sunny aspect gradients, as well as the relationships between environmental factors and biological communities, can reflect the roles and contributions of various biological, physical, and chemical factors.
Most ecological studies on different aspects have predominantly focused on soil stoichiometry characteristics, vegetation distribution and succession, plant community diversity and structure, aboveground plant biomass changes, and plant functional traits [11,12,13]. Zhao et al. [14] focused on vegetation–microorganism–aspects complexes and concluded that soil microbial activity in mountain ecosystems was profoundly affected by these complexes. However, previous studies have not fully solved the problem of how aspects affect biodiversity directly or indirectly, and little is known about the coexistence of species and mechanisms of community construction. Therefore, it is necessary to study different aspects in typical mountain ecosystems and reveal the relationships among environmental factors (e.g., light, temperature, and soil’s physical and chemical properties), biological factors (e.g., plant community composition, stoichiometric characteristics, and functional traits), and community composition and distribution rules of specific biotas.
Toona ciliata Roem. var. pubescens (Franch.) Hand.-Mazz. (Tc) is regarded as a national secondary protected plant in China. Moreover, Tc is an endangered species and a valuable woody plant of the Toona genus in the Meliaceae family [15]. Owing to a long period of exhaustive utilization and exploitation, habitats of Tc forest are affected by these effects and have disadvantageous factors for natural regeneration, meaning that the distribution area of natural Tc forests is reduced and fragmented [16]. Wild Tc is predominantly distributed on both sides of valleys and streams in virgin forests in subtropical China, breeding predominantly by seeds and clonal tillers [15]. Their special habitat preferences and breeding strategies may be the dominant cause of this endangerment. In recent years, scientists have addressed the problems of the natural regeneration of Tc forests [17], and found that transmission barriers [18] and microhabitat restrictions [19] might be major causes of failure of Tc regeneration. Sensitive habitat selection and lower forest regeneration capacity have caused difficulties in the protection and restoration of wild Tc forests [20]. Microhabitats involve many environmental factors, such as water, microorganisms, light, and soil. The growth of little trees is affected by these environmental factors, and there are large variations in the key factors of various plants. Soil microorganisms play a major role in seedling growth, improving soil quality, and promoting natural regeneration in Tc forests [20]. Huang et al. [21] found that soil pathogen infection was the dominant reason for seed rot and the natural regeneration failure of Tc. As an important soil microorganism, arbuscular mycorrhizal (AM) fungi has a symbiotic association with over 80% of terrestrial plants [22], and can effectively promote plant growth and natural regeneration of Tc [23,24]. The symbiotes formed by AM fungi and plants give rise to vital roles in plant nutrient absorption and resistance improvement to external adverse environments at the ecosystem or individual scale [22,25]. Wild Tc forests are largely distributed in habitats which have poor light and soil conditions, and Tc seedlings are often hindered by other plants during the initial growth stage [21]. Considering this, AM fungi may be important in wild Tc forests, as AM fungi can promote plant development and growth by increasing the chlorophyll content in plant leaves, improving photosynthesis, promoting C assimilation and metabolism of plants [26], and improving the competitiveness of plants in natural habitats.
Currently, the study of different aspect-related soil microbial community compositions and their distribution characteristics can only provide some sporadic information [27,28] and only a few studies have investigated the impact of aspect on AM fungal community structure. Chu et al. [28] reported the community composition of AM fungi on sunny and shady aspects in boreal forests of the Greater Khingan Mountains. Similarly, one study reported the distribution of AM fungi on sunny and shady aspects in the arid ecosystem of Inner Mongolia [29]. However, the above studies covered only a few ecosystem types, and studies on particular species and habitats are lacking, limiting our insight into the community-building mechanisms associated with AM fungi. Shady and sunny aspects are comprehensive natural environmental gradients, and the AM diversity and relationship among AM fungi, environmental factors, and plant communities can reveal the comprehensive effects in natural Tc forests.
Considering the particularity of wild Tc forest distribution regions and the important role of AM fungi in natural regeneration and growth of wild Tc forests, this study aimed to reveal the community characteristics and biogeographic distribution rules (sunny and shady aspects) of AM fungi in a wild Tc forests. A gradient of shady and sunny aspects in the Guanshan National Nature Reserve in Jiangxi Province was selected as the research site. Field investigations of AM fungal community and diversity distribution patterns in the roots and rhizosphere soils of seven widespread plants in natural Tc forest were analyzed using the Illumina PE250 high-throughput sequencing method. The rich, continuous environmental (biotic and abiotic) gradients and distinct vegetation types in this region provide an excellent opportunity to discover the dynamics of AM fungi on different aspects in natural Tc forests. Our results will provide important information on the ecological processes of AM fungal community construction in natural Tc forests and a reference for discussing the microbiological mechanism of adaptation of Tc to microhabitats. The null hypothesis was that Tc in shady and sunny aspects would be associated with similar AM fungi diversity and composition, but that Tc recruited significantly different AM fungi compared to other plants. Changes in soil chemical properties and tree species may be responsible for these variations. The findings of this study will improve our ability to understand ecosystem dynamics and function in wild Tc forests and facilitate scientific protection and cultivation strategy decisions on subtropical mountain ecosystems in China.

2. Materials and Methods

2.1. Study Site

The Guanshan National Nature Reserve (28°30′–28°40′ N, 114°29′–114°45′ E) is located in the western part of the Jiuling Mountains in the northwest Jiangxi Province, China. It was one of the first national nature reserves in Jiangxi Province and belongs to the subtropical evergreen broad-leaved forest ecosystem. The mean annual temperature is 17.2 °C. The mean annual precipitation is 1680.2 mm. The total area is 11,500.5 ha, comprising 3621.1 ha in the core area, 1466.4 ha in the buffer zone, and 6413.0 ha in the experimental area. The soil type of this study area was typical red soil.

2.2. Experimental Design and Sampling

The natural Tc forest community is predominantly distributed in a subtropical evergreen broad-leaved forest (Figure 1). The Tc is dominantly distributed sporadically on both sides of valleys and streams and on both sides of the aspects with different orientations and illumination and hydrothermal conditions. The photophilous tree species is mainly distributed in the arborous layer. A field survey was started at 2017, and eight experimental sample plots of 20 × 20 m were set up in a natural Tc forest (north-facing [NF] aspect—4, and south-facing [SF] aspect—4). In the sample plot, we selected seven dominant plants that were distributed in NF and SF, as follows: Toona ciliates (Tc), Padus buergeriana (Pb), Maesa japonica (Mj), Melliodendron xylocarpum (Mx), Ilex chinensis (Ic), Alniphyllum fortunei (Af), and Mallotus japonicas (Mj2). At least one individual sample from each plant was selected in each plot. Five individual root samples and five individual rhizosphere soil samples were randomly selected from each plant, respectively. Fine fibrous roots in the 0–20 cm soil depth were collected from the root samples, and 1–2 kg of soil around the root was collected with a clean small steel shovel to extract AM fungal spores. Finally, 70 root samples and 70 rhizosphere soil samples per plot (2 aspects × 7 plants × 5 individuals) were collected. The collected fine roots were rinsed with water and stored at −80 °C for DNA extraction. The collected soil samples were passed through a 2 mm sieve to eliminate stones, roots, and litter for DNA extraction (−80 °C) and soil experiments (air-dried soil, 0.149 mm mesh).

2.3. Soil Chemical Property Determination

The soil pH was measured using a glass electrode at 1:2.5 (w/v). The soil organic matter (SOM, g·kg−1) content was measured using the K2Cr2O7 heating method. Soil available phosphorus (AP) content was determined after extraction with 0.5 mol·L−1 NaHCO3. Total potassium (TK) content was extracted using a 1 M NH4OAc solution (pH 7.0) and measured by flame photometry. Total phosphorus (TP) content was determined using the H2SO4-HClO4 heating digestion method, and the total nitrogen (TN) content was measured using the Kjeldahl method. The detailed steps can be found in SD Bao [30].

2.4. DNA Extraction Method

Soil DNA was extracted using a Fast DNA SPIN Kit (MP Biomedicals LLC, Santa Ana, CA, USA) according to the instructions and using 0.5 g fresh soil samples. Root DNA was extracted using a Fast Plant Kit (Beijing Tiangen) according to the instructions and using 0.05 g frozen root samples. The total concentration of the extracted soil/root DNA in each soil or root sample was determined using a NanoDrop ND-8000 spectrophotometer (NanoDrop, Wilmington, DE, USA). The DNA quality was checked by 1% agarose gel electrophoresis, diluted to 10–20 ng·μL−1 with ultrapure water, and preserved at −20 °C in a refrigerator for molecular biological analysis.

2.5. PCR Reaction Method

To guarantee the reliability and accuracy of the further data analysis, the PCR reaction needs to meet the following two requirements: (1) use a small number of amplification cycles and (2) ensure that the number of cycles amplified by each sample is consistent. The SSU rDNA of AM fungi in soil and plant roots was amplified using AM FUNGI-specific primers AMV4.5NF-F (5′-AAGCTCGTAGTTGAATTTCG-3′) and AMDGR-R (5′-CCCAACTATCCCTATTAATCAT-3′). Then, PCR was performed using TransStart Fastpfu DNA Polymerase, with a 20 μL reaction system. Each sample was repeated three times. The reaction system is shown in Table 1.
The PCR reaction system was mixed mildly and put into the PCR instrument (ABI GeneAmp® 9700 type), and the reaction conditions were as follows: (1) 5 min of incipient denaturation at 95 °C; (2) 30 s of denaturation at 95 °C; (3) 30 s of annealing at 55 °C; (4) 45 s of elongation at 72 °C; (5) 27 cycles from step 4 to step 2; and (6) incubation at 72 °C for 10 min. Each sample was subjected to three replicates, and the PCR products of the same sample were blended and detected using 2% agarose gel electrophoresis.

2.6. Illumina High-Throughput Sequencing and Bioinformatics

The PCR products were subjected to 1% agarose gel electrophoresis in 1 × TAE buffer (added with Biotake producing a green as blue nucleic acid dye). After electrophoresis at 125 V for 30 min, the PCR results of the gel were observed under UV light. Then, the target band for recovering the PCR products was cut off using a AxyPrepDNA gel recovery kit (AXYGEN), which was eluted with Tris-HCl and detected by 2% agarose electrophoresis. The PCR products were observed and identified using 2% agarose gel, while the electrophoresis condition was 120 V for 40 min. The target bands were separated by electrophoresis, cut, and purified (Aidlab Biotechnologies Co., Ltd., Tianjin, China). The purified product was quantified, homogenized using a TBS-380 fluorescence spectrometer, and sequenced using the Illumina PE250 high-throughput sequencing platform (Shanghai BIOZERON Co., Ltd., Shanghai, China) referring to standard methods. A 97% identity threshold was used to group sequences into operational taxonomic units (OTUs) [31]. The largest ample sequence from each OTU was regarded as the typical sequence for the OTU. The Usearch (version 7.1 http://drive5.com/uparse/, accessed on 25 December 2017) software was used to cluster sequences. Taxonomy was allotted to the fungal OTUs against a subset of the Silva 104 database (http://www.arb-silva.de/download/archive/qiime/, accessed on 28 December 2017). The GenBank (http://www.ncbi.nlm.nih.gov/, accessed on 31 December 2017) for the OTUs was used to confirm the typical sequences that could not be recognized at the family or class level in the fungal database. The original sequences described here can be accessed through GenBank with SRP277481 and accession numbers SAMN15815414 to SAMN15815425 and SAMN15815492 to SAMN15815503.

2.7. Calculation Methods of Diversity Correlation Indexes Data Analysis

We used the Chao index to calculate the total number of species in an ecosystem [32]. Here, Chao 1 was used to calculate the number of OTUs, as follows:
S c h a o 1 = S o b s + n 1 ( n 1 1 ) 2 ( n 2 + 1 )
where Schao1 is the number of OTUs, Sobs is the number of observed OTUs, n1 is the number of OTUs containing only one sequence (singletons), and n2 is the number of OTUs containing only two sequences (doubletons).
The Simpson index was used to quantitatively calculate the microbial diversity [33]. The larger the Simpson index, the lower the community diversity, as follows:
D s i m p s o n = i = 1 S o b s n i ( n i 1 ) N ( N 1 )
where Sobs is the number of observed OTUs, ni is the sequence number of OTU-i, and N is the total number of sequences.
The Shannon index (i.e., the Shannon–Weiner index) was also applied to quantitatively calculate the microbial diversity. The larger the Shannon index, the greater the community diversity, as follows:
H s h a n n o n = i = 1 S o b s n i N ln n i N
where Sobs is the number of observed OTUs, ni is the sequence number of OTU-i, and N is the total number of sequences.
Subsequently, the coverage index was calculated as follows:
C = 1 n 1 N
where C is the coverage of OTUs, n1 is the number of OTUs containing only one sequence (singletons), and N is the total sequence.

2.8. Data Analysis

Firstly, the normal distribution test is carried out on all the data used in the analysis. If the data does not conform to the normal distribution, the natural logarithm or trigonometric function is used for data conversion. Secondly, the Duncan’s multiple comparative analysis (version 26.0; IBM, Armonk, NY, USA) was conducted to analyze the variations in soil chemical properties (pH, SOM, AP, TK, TP, and TN), soil and root AM fungi community diversity (the Chao, Shannon, Simpson, and coverage indices), and sequence parameters (the total sequence number, the number of Glomeromycota sequences, the proportion of Glomeromycota sequences, the number of OTUs, and the relative abundance of the belowground AM fungi community sequence) in Tc natural forest in different aspects. Pairwise comparative analysis was used to compare the difference of the corresponding value between the north and south aspects. Thirdly, there were, respectively, 18 response factors for AM fungal community characteristics in soil and root. Through factor analysis, 6 types of comprehensive factors (which mainly reflected the information of 11 key factors) were used to reflect the information of 18 response factors, so as to reduce the dimension, and the key response factors were determined for further analysis. Therefore, the redundancy analysis (RDA) was applied to analyze the relationships between soil environmental variables and the key characteristics of the AM fungal communities by using the standardized data (Canoco 5.0, Biometrics Wageningen, The Netherlands).

3. Results

3.1. Soil Chemical Properties

In general, soil TN, TP, and AP differed significantly with aspect, being greater on NF than on SF (p < 0.05, Table 2); moreover, the significant effects mainly showed in soils of Tc, Pb, Ic, and Mj (Figure 2). However, there was no significant effect of species on soil TN and TP (p > 0.05, Table 2). Although soil TK showed no obvious difference with aspect (p > 0.05), significant species effect and interaction were found in Table 2 and Figure 2 (p < 0.05). For example, the soil TK of MJ2 was significantly higher than other trees on SF. Additionally, regarding soil pH and SOM, there was no significant difference between SF and NF (p > 0.05, Table 2 and Figure 2).

3.2. Sequence Characteristics of AM Fungi

The basic sequence information of soil AM fungal communities showed no significant difference with aspect (p > 0.05) (Table 2 and Table 3). Due to the Glomeromycota sequences making up the largest proportion of the total sequences, we conducted the related analysis (Table 3). Although the number of Glomeromycota sequences showed no obvious difference with aspect (p > 0.05), a significant interaction between aspect and species was found in Table 2 and Table 3 (p < 0.05). The number of Glomeromycota sequences of Tc and Mj2 was significantly higher on SF than on NF (p < 0.05, Table 3).
The sequence information of the root AM fungal communities in Tc natural forest at different aspects is shown in Table 4. On the whole, the number of sequence differed significantly with aspect (p < 0.05, Table 2 and Table 4), being greater on SF than on NF. However, there was no significant effect of species on the number of sequence (p > 0.05, Table 4). The Glomeromycota sequence proportion showed significant differences between SF and NF (p < 0.05, Table 2 and Table 4). Moreover, the significant effect of species on the Glomeromycota sequence proportion was found in Mx (p < 0.05, Table 4).

3.3. OTUs Analysis of AM Fungi

Although the relative OTU abundance of AM fungal communities of soil and root showed no difference between SF and NF (p > 0.05), the significant effect of species was found in soil (p < 0.05) (Table 2 and Table 5). Moreover, we found that the relative OTU abundance in rhizosphere soil of Tc was significantly higher for NF than for SF (p < 0.05, Table 5).

3.4. Analysis of AM Fungal Community Composition

At the family level, the relative abundance of each composition changed obviously in the soil AM fungal community in the Tc natural forest at different aspects (Figure 3). Glomeraceae was the most common and plentiful family in the Tc natural forest soil, but there was no significant difference between SF and NF (p > 0.05, Table A1 and Table A2). The relative abundance of unclassified fungi showed a significant interaction between aspect and species (p < 0.05, Table A2).
Similarly, Glomeraceae was the most common and abundant AM fungi group in the Tc natural forest root (Figure 3). The relative abundance of Glomeraceae in the Tc root was 92.68% in SF, and it was 90.10% in NF. Generally, the relative abundance of Acaulosporaceae differed significantly with aspect, being greater on NF than on SF (p < 0.05, Table A1 and Table A2, Figure 3). However, there was no significant effect of species on the relative abundance of Acaulosporaceae (p > 0.05, Table A2). In addition, the relative abundance of Archaeosporaceae showed significant differences among different species (p < 0.05, Table A1 and Table A2, Figure 3).

3.5. Diversity of AM Fungi Community

The diversity of the soil and root AM fungal communities in the Tc natural forest at different aspects is shown in Table 6 and Table 7. The coverage index of the root AM fungal communities differed significantly with aspect (p < 0.05, Table 2 and Table 7), being greater on SF than on NF. Moreover, there was significant effect of species on the coverage index of the soil and root AM fungal communities (p < 0.05, Table 2, Table 6 and Table 7). Although the Chao index of the root AM fungal communities showed no obvious difference with aspect (p > 0.05, Table 2), a significant species effect was found in Table 2 and Table 7 (p < 0.05). In addition, the Shannon and Simpson indices of the soil and root AM fungal community found no significant difference between SF and NF (p > 0.05).

3.6. Redundancy Analysis in Different Aspects

Through the factor analysis, we extracted the key response factors for RDA analysis. The environmental variables (aspects and soil physicochemical properties) and response variables (characteristics of AM fungal communities of dominant tree species in Tc natural forest) were analyzed using RDA (Figure 4). The results showed that soil AM fungi communities had a higher total sequence numbers, Glomeromycota sequence numbers, Chao index, Shannon index, Archaeosporaceae abundance, and Acaulosporaceae abundance in SF, while the opposite trend was shown in NF (Figure 4A). Furthermore, SOM, TK, TN, TP, and AP were negatively correlated with total sequence numbers, Glomeromycota sequence numbers, Chao index, Shannon index and Archaeosporaceae abundance, and negatively correlated with the relative abundance of unclassified fungi (Figure 4A). The SOM and TK were positively correlated with OTUs (Figure 4A).
The RDA results (Figure 4B) showed that the proportion of Glomeromycota sequences, Glomeromycota sequence numbers, Diversisporaceae abundance, and Shannon index of the root AM fungi community were relatively high in NF, while total sequence numbers, coverage index, Gigasporaceae abundance, and unclassified fungi abundance of the root AM fungi community were relatively high on SF. The SOM was positively correlated with Glomeromycota sequence numbers and Glomeraceae abundance (Figure 4B). The AP was positively correlated with coverage index, total sequence numbers, and Gigasporaceae abundance, while it was negatively correlated with unclassified fungi abundance, the proportion of Glomeromycota sequences, Diversisporaceae abundance, and Shannon index (Figure 4B). Moreover, TN, TK, and pH showed the opposite trend with AP.

4. Discussion

4.1. Soil AM Fungal Community Diversity and Composition in Tc Forest at Different Aspects

An Illumina PE250 high-throughput sequencing method was applied to explore the features of the AM fungi community in a natural Tc forest at different aspects and to analyze the AM fungal community characteristics and diversity changes. The AM fungal community structure refers to the composition and quantitative characteristics of the fungal community in each habitat, and is mainly characterized by abundance, density, frequency, dominance, importance value, and community coefficient [34]. Stability is an important factor in achieving ecological balance. Therefore, studying the structural changes in AM fungal communities in specific ecosystems and their influencing factors can provide a basis for optimizing the natural Tc forest ecosystem structure and regulating ecological functions. Indeed, AM fungal species diversity is an important component of soil microbial diversity and has a vital ecological and physiological role in providing plant nutrient uptake, enhancing soil structure stability, and changing the global C and N cycles [35]. Many factors influence the diversity of AM fungal communities, including nutrient inputs [36], anthropogenic disturbances [37], transmission limitations [38], plant communities [39], and other environmental factors [40].
Aspect, one of the most common environmental impact factors of natural mountain forests, influences the diversity of AM fungal communities in plant soils [41]. Soil microorganisms are seen as the driving force of productivity diversity in terrestrial ecosystems [42]. In this study, species richness and diversity indices were used to describe the diversity characteristics of AM fungal communities in different aspects of Tc natural forests, including Chao, Shannon, and Simpson. In terms of species richness, the number of OTUs of AM fungal communities in the inter-rhizosphere soil of Tc was significantly higher on NF (16.4%) than on SF (13.7%). In terms of diversity indices, the Chao, Shannon, and Simpson indices of AM fungal communities in the rhizosphere soil and root of Tc were not significantly different between NF and SF. This phenomenon may be because the composition of the AM fungi community in the roots and rhizosphere soil primarily consisted of Glomeraceae, reaching 45.07%–86.39% and 74.82%–93.08%, respectively. Some scholars have studied the AM fungal diversity and community structure of farmland, wetland, forest, and grassland, and found that the AM fungal community composition is mainly determined by Glomeraceae [43,44], which is the same as the result of this study. The reason is that the Glomeraceae has a unique reproductive mode and can tolerate different environmental changes [45]. However, the rhizosphere soil contained more rare AM fungal taxa with low relative abundance, which lowered the community diversity index (Figure 3; Table 5). We found that there were a large number of rare AM fungal taxa with very few sequences in the Tc natural forest, and these rare taxa may be the main reason for the functional redundancy in the Tc natural forest, which should be investigated in a follow-up study. The results suggest that the effects of various environmental factors on AM fungal communities in natural ecosystems, especially in complex forest ecosystems, should not be evaluated simply in terms of diversity or species richness but should be considered comprehensively to obtain more accurate and objective analysis results.

4.2. Effects of Aspects and Other Factors on the Characteristics of AM Fungal Community

The AM fungi are mainly distributed in the soil around the roots; therefore, they are closely related to the soil environment. Soil type, depth, water status, fertility, pH, and other factors have important effects on the AM fungal community [46]. In a past study, ectomycorrhizal fungal diversity was higher on sunny aspects than on shady aspects, probably due to differences in the distribution of plants in different aspect orientations caused by differences in insolation, water, and heat conditions [47]. This study found that soil nutrient characteristics were relatively better on NF than on SF. The differences may be due to the variations of space and climate which were observed in the selected Tc community distribution area.
This is similar to the results of other studies, which found that the taxonomy and spectral composition of AM fungal communities significantly differed between sunny and shady aspects through NMDS sorting and PERMANOVA analysis, with sunny aspects having higher light intensity, soil temperature, and soil pH, and shady aspects having higher soil moisture and soil fertility index [48]. Community genealogy analysis also showed that AM fungal communities were spectrally aggregated on sunny aspects and spectrally random on shady aspects, indicating that the primary ecological process of AM fungal community construction is a shift from ecological niche-dominated filtering to a synthesis of competitive exclusion and ecological niche filtering [49]. The diversity of mycorrhizal fungal communities (evenness index, Simpson index, Shannon–Wiener index, and richness) in different areas differed according to aspect, altitude, soil moisture, and human disturbance [50].

4.3. Effects of Unique Tree Species Characteristics of Tc on AM Fungal Community Characteristics

The characteristics of Tc and its symbiotic relationship with AM fungi jointly affected the soil AM fungal community. Symbiotics play a vital role in plant nutrient uptake, improving ecosystem or individual resistance to adverse environments [22,25], which can help break the regeneration dilemma caused by transmission barriers [21] and microhabitat restrictions in Tc forests [51], and are conducive to the protection and restoration of wild Tc forests [52]. The Tc is the main tree species in the Tc natural forest, but its seeds are susceptible to soil pathogens, which affects its natural forest regeneration ability. Therefore, it does not have an absolute advantage in community competition [53]. In different aspects in the natural Tc forest, the diversity and relative abundance of AM fungal communities in the roots and rhizosphere soil of Tc were only at the medium level compared to those of the other six species. This is because there is a certain selectivity between different mycorrhizal fungi and host plants, which affects the binding ability of plants to mycorrhizal fungi [54], and different plants have different adaptabilities to different AM fugal community [55] and aspect habitats. In addition, the unclassified fungi were found less in roots (0%–1%) than in rhizosphere soil (6%–17%). Published research recognized more new species of Glomeromycota in the rhizospheres of endemic plants [56]. Therefore, Tc still has a certain development scope for AM fungal infections. It can enhance its competition for environmental resources by increasing AM fungal infection, which may be of great value for the natural regeneration of Tc and protection of the Tc natural community.

4.4. AM Fungal Library in SF and ‘Functional’ AM Fungi in NF

The AM fungi in soil are often considered a library, and their significance is to preserve and maintain the diversity of AM fungal resources. The AM fungi in roots are often considered functional AM fungi, and their presence helps realize some plant functions. In this study, the number of OTUs (16.4%) of AM fungal communities in rhizosphere soil on NF of Tc was significantly larger than that on SF (13.7%), whereas the number of OTUs (11.3%) of AM fungal communities in roots on NF of Tc was lower than on SF (15.1%). This difference is precisely why taxonomic and functional diversities are not synchronized [57] with aspect. Moreover, abiotic factors, such as soil physicochemical properties, had an important influence on AM fungal community composition [58]. In this study, the low insolation, high moisture, and high nutrient content (TP, TK, and TN) conditions on NF might be beneficial for the functional realization of AM fungi. In contrast, the high insolation, low soil moisture, and low soil nutrient content conditions (TP, TK, and TN) on SF might be conducive to the proliferation and preservation of AM fungal groups. For example, one study showed that when SOM content was high, the associated microorganisms of AM fungi promote mass-production due to more nutrients, thus, inhibiting the infiltration of roots of AM fungi [59]. The changes in vegetation and environment factors did not lead to the extinction of a large number of soil AM fungi, but only caused the change in their relative proportion. Therefore, the results of this study found that the soil AM fungal pool might be relatively stable in the Tc natural forest.

4.5. Implications

The rhizosphere is the main medium for the interaction between plant and soil. A complex feedback process occurs at the interface between plant and soil. Soil fungi form a dense underground mycelium bridge network in soil through coupling with plant roots, thus, promoting water and nutrient absorption by the plant [60]. Most terrestrial plants can form mycorrhizal symbionts with AM fungi, and the fungi cannot grow and complete their life cycle without host plants [22]. In this study, Illumina high-throughput sequencing technology was used to study the composition and functional characteristics of the AM fungal communities in the natural forest of Tc. The results of this study can provide a theoretical basis for better understanding of the natural regeneration law of the community of Tc. We provide a reference for ecological restoration, protection, and utilization of Tc from the perspective of underground microbial diversity, as well as for the protection of rare and endangered plants, sustainable use of forests, maintenance of ecosystem balance and protection of biodiversity.

5. Conclusions

We analyzed the diversity and composition of AM fungal communities in a natural Tc forest using Illumina high-throughput sequencing of rhizosphere soil and roots from north-facing and south-facing aspects. The number of Glomeromycota sequences was significantly larger on SF than on NF, while the number of OTUs on NF was significantly higher than on SF in Tc soil. The composition of the AM fungi community in roots and rhizosphere soil was primarily consisted of Glomeraceae, reaching 45.07%–86.39% and 74.82%–93.08%, respectively. Moreover, the relative abundance of unclassified fungi in rhizosphere soil of Tc on NF was significantly higher than that on SF. However, the Chao, Shannon, and Simpson indices of the soil and root AM fungal community found no significant difference between SF and NF. Furthermore, our study revealed that the low insolation, high soil moisture, and high nutrient content conditions on NF might be beneficial for the functional realization of AM fungi. In contrast, the high insolation, low soil moisture, and low soil nutrient content conditions on SF might be beneficial for the proliferation and preservation of AM fungal groups. In the future, how topography and various environmental factors to regulate AM fungal function should be given attention due to the fact that AM symbiosis plays non-negligible roles in ecosystems.

Author Contributions

Q.W., J.P., W.L. and L.Z. conceptualized the study, and J.P. identified, processed, and collected samples with help from S.Y. and T.L.; Q.W. and Y.K. primarily interpreted the data with help from P.J.M. and L.Z.; Q.W., J.P., L.Z. and W.L. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (31860202) and the Forestry Science and Technology Innovation Special Project of Jiangxi Provincial Department of Forestry (Innovation Special Project [2019] No. 25).

Data Availability Statement

Not applicable.

Acknowledgments

We would especially like to thank the Administration of Guanshan National Natural Reserve for assisting in the investigation.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Analysis of relative abundance of belowground AM fungal communities’ sequence in dominant tree species in T. ciliate natural forest in different aspects.
Table A1. Analysis of relative abundance of belowground AM fungal communities’ sequence in dominant tree species in T. ciliate natural forest in different aspects.
SpeciesAspectsUnclassifiedParaglomeraceaePacisporaceaeGlomeraceaeGigasporaceaeGeosiphonaceaeDiversisporaceaeArchaeosporaceaeAmbisporaceaeAcaulosporaceae
Soil
Toona ciliatesSouthbaaabaaaabaa
NorthA *BAAAAAAAA
Padus buergerianaSouthaaacaaaabaa
NorthAABAAAAAAAA
Maesa japonicaSouthbaaabcaaaabaa
NorthAABAAAAAAAA
Melliodendron xylocarpumSouthbaaabcaaaabaa
NorthAABAAAAAAAA
Ilex chinensisSouthbaabcaaaabaa
NorthABAAAAAAAA
Alniphyllum fortuneiSouthbaabcaaaaaa
NorthAABAAAAAAAA
Mallotus japonicusSouthbaaaaaabaa
NorthA *A *AA *AAAAAA
Root
Toona ciliatesSouthaaaaaaabaa
NorthAAAABBBAABA
Padus buergerianaSouthaaaaaaaaaa
NorthAAAABBBABA
Maesa japonicaSouthaaaaaaaabaa
NorthAAAABABABA
Melliodendron xylocarpumSouthaaaaaaabaa
NorthAAAABBAAABA
Ilex chinensisSouthAAaaaaaabaa
NorthA *AAABBBAAA
Alniphyllum fortuneiSouthaaaaaaaabaa
NorthAAAAAABBA *BA *
Mallotus japonicusSouthaaaaaaaabaa
NorthAAAABBBABA
Note: Data with different capital letters means significant difference in NS at p < 0.05. Data with different lowercase letters means significant difference in SS at p < 0.05; * means that there is significant difference between SS and NS in Tc.
Table A2. Results of variance analysis for the effects of aspect and species on the relative abundance of belowground AM fungal communities’ sequence.
Table A2. Results of variance analysis for the effects of aspect and species on the relative abundance of belowground AM fungal communities’ sequence.
FactorsAspectSpeciesAspect × Species
FpFpFp
S-Acaulosporaceae2.939NS1.299NS0.606NS
S-Ambisporaceae2.171NS1.326NS0.443NS
S-Archaeosporaceae1.503NS0.938NS1.292NS
S-Diversisporaceae0.468NS0.755NS0.865NS
S-Geosiphonaceae0.274NS0.381NS1.170NS
S-Gigasporaceae1.338NS0.757NS0.425NS
S-Glomeraceae0.419NS1.478NS1.789NS
S-Pacisporaceae1.284NS0.754NS0.754NS
S-Paraglomeraceae1.582NS0.416NS1.975NS
S-Unclassified fungi0.465NS1.092NS2.3730.049
R-Acaulosporaceae7.5980.0090.903NS0.705NS
R-Ambisporaceae1.934NS1.779NS2.057NS
R-Archaeosporaceae2.516NS3.9960.0360.640NS
R-Diversisporaceae0.743NS0.811NS1.175NS
R-Geosiphonaceae0.276NS1.331NS1.866NS
R-Gigasporaceae0.535NS1.299NS0.652NS
R-Glomeraceae0.112NS1.319NS1.042NS
R-Pacisporaceae1.274NS1.318NS1.336NS
R-Paraglomeraceae1.337NS1.028NS1.427NS
R-Unclassified fungi3.986NS0.676NS0.800NS
Abbreviations are as follows: S, soil; R, root; NS means that there is no significant difference (p > 0.05).

References

  1. Cardinale, B.J.; Duffy, J.E.; Gonzalez, A.; Hooper, D.U.; Perrings, C.; Venail, P.; Narwani, A.; Mace, G.M.; Tilman, D.; Wardle, D.A.; et al. Biodiversity loss and its impact on humanity. Nature 2012, 486, 59–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Jing, X.; Sanders, N.J.; Shi, Y.; Chu, H.; Classen, A.T.; Zhao, K.; Chen, L.; Shi, Y.; Jiang, Y.; He, J.S. The links between ecosystem multifunctionality and above- and belowground biodiversity are mediated by climate. Nat. Commun. 2015, 6, 8159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Mvondo Owono, F.; Ntamak-Nida, M.-J.; Dauteuil, O.; Guillocheau, F.; Njom, B. Morphology and long-term landscape evolution of the South African plateau in South Namibia. Catena 2016, 142, 47–65. [Google Scholar] [CrossRef]
  4. Bennie, J.; Huntley, B.; Wiltshire, A.; Hill, M.O.; Baxter, R. Slope, aspect and climate: Spatially explicit and implicit models of topographic microclimate in chalk grassland. Ecol. Model. 2008, 216, 47–59. [Google Scholar] [CrossRef]
  5. Liu, S.L.; An, N.N.; Yang, J.J.; Dong, S.K.; Wang, C.; Yin, Y.J. Prediction of soil organic matter variability associated with different land use types in mountainous landscape in southwestern Yunnan province, China. Catena 2015, 133, 137–144. [Google Scholar] [CrossRef]
  6. Zhang, Z.W.; Yin, H.Y.; Chang, J.J.; Xue, J. Spatial variability of surface soil water content and its influencing factors on shady and sunny slopes of an alpine meadow on the Qinghai–Tibetan Plateau. Glob. Ecol. Conserv. 2022, 34, e02035. [Google Scholar] [CrossRef]
  7. Wei, T.F.; Shangguan, D.H.; Tang, X.L.; Qin, Y. The role of glacial gravel in community development of vascular plants on the glacier forelands of the Third Pole. J. Arid Land 2022, 14, 1022–1037. [Google Scholar] [CrossRef]
  8. Conforti, M.; Lucà, F.; Scarciglia, F.; Matteucci, G.; Buttafuoco, G. Soil carbon stock in relation to soil properties and landscape position in a forest ecosystem of southern Italy (Calabria region). Catena 2016, 144, 23–33. [Google Scholar] [CrossRef]
  9. Huang, Y.M.; Liu, D.; An, S.S. Effects of slope aspect on soil nitrogen and microbial properties in the Chinese Loess region. Catena 2015, 125, 135–145. [Google Scholar] [CrossRef]
  10. Small, C.J.; McCarthy, B.C. Spatial and temporal variability of herbaceous vegetation in an eastern deciduous forest. Plant Ecol. 2003, 164, 37–48. [Google Scholar] [CrossRef]
  11. Måren, I.E.; Karki, S.; Prajapati, C.; Yadav, R.K.; Shrestha, B.B. Facing north or south: Does slope aspect impact forest stand characteristics and soil properties in a semiarid trans-Himalayan valley? J. Arid Environ. 2015, 121, 112–123. [Google Scholar] [CrossRef] [Green Version]
  12. Pereira, P.; Cerda, A.; López, A.J.; Zavala, L.M.; Mataix-Solera, J.; Arcenegui, V.; Misiūnė, I.; Keesstra, S.D.; Novara, A.J.L.D. Short-Term Vegetation Recovery after a Grassland Fire in Lithuania: The Effects of Fire Severity, Slope Position and Aspect. Land Degrad. Dev. 2016, 27, 1523–1534. [Google Scholar] [CrossRef]
  13. Wang, Y.Y.; Yu, J.H.; Xiao, L.; Zhong, Z.L.; Wang, Q.; Wang, W.J. Dominant species abundance, vertical structure and plant diversity response to nature forest protection in Northeastern China: Conservation effects and implications. Forests 2020, 11, 295. [Google Scholar] [CrossRef] [Green Version]
  14. Zhao, N.N.; Guggenberger, G.; Shibistova, O.; Thao, D.T.; Shi, W.J.; Li, X.G. Aspect-vegetation complex effects on biochemical characteristics and decomposability of soil organic carbon on the eastern Qinghai-Tibetan Plateau. Plant Soil 2014, 384, 289–301. [Google Scholar] [CrossRef]
  15. Liu, J.; Jiang, J.; Chen, Y. Genetic diversity of central and peripheral populations of Toona ciliata var. pubescens, an endangered tree species endemic to China. Genet. Mol. Res. GMR 2014, 13, 4579–4590. [Google Scholar] [CrossRef]
  16. Zhan, X.; Li, P.; Hui, W.K.; Deng, Y.W.; Gan, S.M.; Sun, Y.; Zhao, X.H.; Chen, X.Y.; Deng, X.M. Genetic diversity and population structure of Toona ciliata revealed by simple sequence repeat markers. Biotechnol. Biotechnol. Equip. 2019, 33, 214–222. [Google Scholar] [CrossRef] [Green Version]
  17. Huang, H.L.; Zhang, L.; Jia, L.M.; Liang, Y.L.; Cai, J.H. Fitness of sexual reproduction of Toona ciliata var. pubescens natural populations and their sexual reproduction and regeneration. J. Appl. Ecol. 2018, 29, 1098–1106. [Google Scholar] [CrossRef]
  18. Münzbergová, Z.; Herben, T. Seed, dispersal, microsite, habitat and recruitment limitation: Identification of terms and concepts in studies of limitations. Oecologia 2005, 145, 1–8. [Google Scholar] [CrossRef]
  19. Marques, M.C.M.; Burslem, D.F.R.P. Multiple stage recruitment limitation and density dependence effects in two tropical forests. Plant Ecol. 2015, 216, 1243–1255. [Google Scholar] [CrossRef]
  20. Guo, X.Y.; Fu, L.; Zhang, L.; Su, H.; Liang, Y.L. Effects of forest soil and soil-borne fungi on seed germination and seedling survival of Toona ciliata var. pubescens. For. Res. 2017, 30, 871–877. [Google Scholar]
  21. Huang, H.L.; Zhang, L.; Liao, C.K. Seed rain, soil seed bank, and natural regeneration of natural Toona ciliata var. pubescens forest. Chin. J. Appl. Ecol. 2012, 23, 972–978. [Google Scholar]
  22. Smith, S.E.; Read, D.J. Mycorrhizal Symbiosis; Academic Press: Cambridge, MA, USA, 2008. [Google Scholar]
  23. Mridha, M.A.U.; Dhar, P.P. Biodiversity of arbuscular mycorrhizal colonization and spore population in different agroforestry trees and crop species growing in Dinajpur, Bangladesh. J. Res. 2007, 18, 91–96. [Google Scholar] [CrossRef]
  24. Pan, J.; Wang, Q.; Guo, X.; Jiang, X.; Cheng, Q.; Fu, L.; Liu, W.; Zhang, L. Local patterns of arbuscular mycorrhizal fungal diversity and community structure in a natural Toona ciliata var. pubescens forest in South Central China. PeerJ 2021, 9, e11331. [Google Scholar] [CrossRef]
  25. Hoeksema, J.D.; Chaudhary, V.B.; Gehring, C.A.; Johnson, N.C.; Karst, J.; Koide, R.T.; Pringle, A.; Zabinski, C.; Bever, J.D.; Moore, J.C.; et al. A meta-analysis of context-dependency in plant response to inoculation with mycorrhizal fungi. Ecol. Lett. 2010, 13, 394–407. [Google Scholar] [CrossRef] [PubMed]
  26. Liu, R.J.; Li, X.L. Arbuscular Mycorrhiza and Its Application; Science Press: Beijing, China, 2000. [Google Scholar]
  27. Ehrlich, R.; Schulz, S.; Schloter, M.; Steinberger, Y. Effect of slope orientation on microbial community composition in different particle size fractions from soils obtained from desert ecosystems. Biol. Fertil. Soils 2015, 51, 507–510. [Google Scholar] [CrossRef]
  28. Chu, H.Y.; Xiang, X.J.; Yang, J.; Adams, J.M.; Zhang, K.P.; Li, Y.T.; Shi, Y. Effects of slope aspects on soil bacterial and arbuscular fungal communities in a Boreal forest in China. Pedosphere 2016, 26, 226–234. [Google Scholar] [CrossRef]
  29. Liu, M.; Zheng, R.; Bai, S.; Bai, Y.; Wang, J. Slope aspect influences arbuscular mycorrhizal fungus communities in arid ecosystems of the Daqingshan Mountains, Inner Mongolia, North China. Mycorrhiza 2017, 27, 189–200. [Google Scholar] [CrossRef]
  30. Bao, S.D. Soil and Agricultural Chemistry Analysis; China Agriculture Press: Beijing, China, 2000. [Google Scholar]
  31. Mahmoudi, N.; Cruz, C.; Mahdhi, M.; Mars, M.; Caeiro, M.F. Arbuscular mycorrhizal fungi in soil, roots and rhizosphere of Medicago truncatula: Diversity and heterogeneity under semi-arid conditions. PeerJ 2019, 7, e6401. [Google Scholar] [CrossRef] [Green Version]
  32. Chao, A. Nonparametric-estimation of the number of classes in a population. Scand. J. Stat. 1984, 11, 265–270. [Google Scholar]
  33. Simpson, E.H. Measurement of Diversity. Nature 1949, 163, 688. [Google Scholar] [CrossRef]
  34. Liu, M.; Zheng, R.; Bai, S.L.; Wang, J.G.; Li, L.; Duan, G.Z. Advances of species diversity of arbuscular mycorrhizal fungi. Microbiol. China 2016, 43, 1836–1843. [Google Scholar]
  35. Yang, H.S.; Xiong, Y.Q.; Wang, Q.; Guo, Y.; Dai, Y.J.; Xu, M.M. Arbuscular mycorrhizal fungal species diversity: Ecological functioning, determinants and assembling mechanisms. Acta Ecol. Sin. 2016, 36, 2826–2832. [Google Scholar]
  36. Veresoglou, S.D.; Chen, B.; Rillig, M.C. Arbuscular mycorrhiza and soil nitrogen cycling. Soil Biol. Biochem. 2012, 46, 53–62. [Google Scholar] [CrossRef]
  37. Kabir, Z. Tillage or no-tillage: Impact on mycorrhizae. Can. J. Plant Sci. 2005, 85, 23–29. [Google Scholar] [CrossRef] [Green Version]
  38. Dumbrell, A.J.; Nelson, M.; Helgason, T.; Dytham, C.; Fitter, A.H. Idiosyncrasy and overdominance in the structure of natural communities of arbuscular mycorrhizal fungi: Is there a role for stochastic processes. J. Ecol. 2010, 98, 419–428. [Google Scholar] [CrossRef]
  39. Řezáčová, V.; Řezáč, M.; Gryndler, M.; Hršelová, H.; Gryndlerová, H.; Michalová, T. Plant invasion alters community structure and decreases diversity of arbuscular mycorrhizal fungal communities. Appl. Soil Ecol. 2021, 167, 104039. [Google Scholar] [CrossRef]
  40. Hawkes, C.V.; Kivlin, S.N.; Rocca, J.D.; Huguet, V.; Thomsen, M.; Suttle, K.B. Fungal community responses to precipitation. Glob. Chang. Biol. 2011, 17, 1637–1645. [Google Scholar] [CrossRef]
  41. Guo, Y.N.; Zhang, H.D.; Bao, Y.Y.; Tan, H.Z.; Liu, X.H.; ur Rahman, Z. Distribution characteristics of soil AM fungi community in soft sandstone area. J. Environ. Manag. 2022, 316, 115193. [Google Scholar] [CrossRef]
  42. 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]
  43. Jeffries, P.; Gianinazzi, S.; Perotto, S.; Turnau, K.; Barea, J.-M. The contribution of arbuscular mycorrhizal fungi in sustainable maintenance of plant health and soil fertility. Biol. Fertil. Soils 2003, 37, 1–16. [Google Scholar] [CrossRef]
  44. Scheublin, T.R.; Ridgway, K.P.; Young, J.P.; van der Heijden, M.G. Nonlegumes, legumes, and root nodules harbor different arbuscular mycorrhizal fungal communities. Appl. Environ. Microbiol. 2004, 70, 6240–6246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Brito, I.; Goss, M.J.; de Carvalho, M.; Chatagnier, O.; van Tuinen, D. Impact of tillage system on arbuscular mycorrhiza fungal communities in the soil under Mediterranean conditions. Soil Till. Res. 2012, 121, 63–67. [Google Scholar] [CrossRef]
  46. Canini, F.; Zucconi, L.; Pacelli, C.; Selbmann, L.; Onofri, S.; Geml, J. Vegetation, pH and water content as main factors for shaping fungal richness, community composition and functional guilds distribution in soils of Western Greenland. Front. Microbiol. 2019, 10, 2348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Voortman, B.R.; Fujita, Y.; Bartholomeus, R.P.; Aggenbach, C.; Witte, J.-P.M. How the evaporation of dry dune grasslands evolves during the concerted succession of soil and vegetation. Ecohydrology 2017, 10, e1848. [Google Scholar] [CrossRef]
  48. Chai, Y.X. Diversity and Community Structure of Arbuscular Mycorrhizal Fungi in Yinyang Slope of Qilian Mountains. Ph.D. Thesis, Lanzhou University, Lanzhou, China, 2018. [Google Scholar]
  49. Li, R.; Liu, M.X.; Zhang, C.; Zhao, R.D.; Shao, P. Distribution characteristics of soil microbial communities along different slope direction of Gannan sub-alpine meadows. Ecol. Environ. Sci. 2017, 26, 1884–1891. [Google Scholar]
  50. Sun, X.F.; Su, Y.Y.; Zhang, Y.; Wu, M.Y.; Zhang, Z.; Pei, K.Q.; Sun, L.F.; Wan, S.Q.; Liang, Y. Diversity of arbuscular mycorrhizal fungal spore communities and its relations to plants under increased temperature and precipitation in a natural grassland. ChSBu 2013, 58, 4109–4119. [Google Scholar] [CrossRef]
  51. Guo, Y.; He, Y.J.; Wu, P.; Wu, B.L.; Lin, Y.; He, M.H.; Han, X.; Xia, T.T.; Shen, K.P.; Kang, L.L.; et al. The interspecific competition presents greater nutrient facilitation compared with intraspecific competition through AM fungi interacting with litter for two host plants in karst soil. J. Plant Ecol. 2021, 15, 399–412. [Google Scholar] [CrossRef]
  52. Liang, Y.M.; He, X.Y.; Chen, C.Y.; Feng, S.Z.; Liu, L.; Chen, X.B.; Zhao, Z.W.; Su, Y.R. Influence of plant communities and soil properties during natural vegetation restoration on arbuscular mycorrhizal fungal communities in a karst region. Ecol. Eng. 2015, 82, 57–65. [Google Scholar] [CrossRef]
  53. Huang, H.L.; Zhang, L.; Guo, X.Y.; Liang, Y.L.; Liu, Z.G. Fruiting characteristics and sexual fecundity of Toona ciliata var. pubescens populations in Jiulianshan national nature reserve. Sci. Silvae Sin. 2013, 49, 170–174. [Google Scholar]
  54. Helgason, T.; Merryweather, J.W.; Denison, J.; Wilson, P.; Young, J.P.W.; Fitter, A.H. Selectivity and functional diversity in arbuscular mycorrhizas of co-occuring fungi and plants from a temperate deciduous woodland. J. Ecol. 2002, 90, 371–384. [Google Scholar] [CrossRef] [Green Version]
  55. Lovelock, C.E.; Andersen, K.; Morton, J.B. Arbuscular mycorrhizal communities in tropical forests are affected by host tree species and environment. Oecologia 2003, 135, 268–279. [Google Scholar] [CrossRef] [PubMed]
  56. López-García, Á.; Varela-Cervero, S.; Vasar, M.; Öpik, M.; Barea, J.M.; Azcón-Aguilar, C. Plant traits determine the phylogenetic structure of arbuscular mycorrhizal fungal communities. Mol. Ecol. 2017, 26, 6948–6959. [Google Scholar] [CrossRef] [PubMed]
  57. Van Diepen, L.; Lilleskov, E.; Pregitzer, K. Simulated nitrogen deposition affects community structure of arbuscular mycorrhizal fungi in northern hardwood forests. Mol. Ecol. 2011, 20, 799–811. [Google Scholar] [CrossRef] [PubMed]
  58. Jamiołkowska, A.; Księżniak, A.; Gałązka, A.; Hetman, B.; Kopacki, M.; Skwaryło-Bednarz, B. Impact of abiotic factors on development of the community of arbuscular mycorrhizal fungi in the soil: A Review. Int. Agrophys. 2018, 32, 133–140. [Google Scholar] [CrossRef]
  59. Wang, Q.; Wang, J.; Li, Y.; Chen, D.; Ao, J.; Zhou, W.; Shen, D.; Li, Q.; Huang, Z.; Jiang, Y. Influence of nitrogen and phosphorus additions on N(2)-fixation activity, abundance, and composition of diazotrophic communities in a Chinese fir plantation. Sci. Total Environ. 2018, 619–620, 1530–1537. [Google Scholar] [CrossRef] [PubMed]
  60. Corrales, A.; Ovrebo, C.L. Fungi of the Fortuna Forest Reserve: Taxonomy and ecology with emphasis on ectomycorrhizal communities. Ecol. Environ. Conserv. 2020, 4, 1–29. [Google Scholar]
Figure 1. Location of Guanshan National Nature and distribution of study sites.
Figure 1. Location of Guanshan National Nature and distribution of study sites.
Forests 13 02100 g001
Figure 2. Soil chemical properties (AP, TK, SOM, pH, TP, and TN) in Tc natural forest in different aspects. The error line is the standard error. Data with different capital letters means significant difference in NF at p < 0.05. Data with different lowercase letters means significant difference in SF at p < 0.05; * means that there is significant difference between SF and NF. Abbreviations are as follows: Tc, Toona ciliates; Pb, Padus buergeriana; Mj, Maesa japonica; Mx, Melliodendron xylocarpum; Ic, Ilex chinensis; Af, Alniphyllum fortunei; Mj2, Mallotus japonicas.
Figure 2. Soil chemical properties (AP, TK, SOM, pH, TP, and TN) in Tc natural forest in different aspects. The error line is the standard error. Data with different capital letters means significant difference in NF at p < 0.05. Data with different lowercase letters means significant difference in SF at p < 0.05; * means that there is significant difference between SF and NF. Abbreviations are as follows: Tc, Toona ciliates; Pb, Padus buergeriana; Mj, Maesa japonica; Mx, Melliodendron xylocarpum; Ic, Ilex chinensis; Af, Alniphyllum fortunei; Mj2, Mallotus japonicas.
Forests 13 02100 g002
Figure 3. The AM fungal community composition in Tc natural forest (at the family level) in different aspects. Tc represents Toona ciliates; Pb represents Padus buergeriana; Mx represents Melliodendron xylocarpum; Mj represents Maesa japonica; Ic represents Ilex chinensis; Af represents Alniphyllum fortunei; Mj2 represents Mallotus japonicas.
Figure 3. The AM fungal community composition in Tc natural forest (at the family level) in different aspects. Tc represents Toona ciliates; Pb represents Padus buergeriana; Mx represents Melliodendron xylocarpum; Mj represents Maesa japonica; Ic represents Ilex chinensis; Af represents Alniphyllum fortunei; Mj2 represents Mallotus japonicas.
Forests 13 02100 g003
Figure 4. Redundancy ordination (RDA) results (A, soil; B, root) of the characteristics of the AM fungi community and aspects and soil physicochemical properties in Tc natural forest. Abbreviatons are as follows: GSP, the proportion of Glomeromycota sequences; NTS, the number of total sequences; NGS, the number of Glomeromycota sequences; Div, Diversisporaceae; Gig, Gigasporaceae; Arc, Archaeosporaceae; Unc, unclassified fungi; Glo, Glomeraceae; Pac, Pacisporaceae; Aca, Acaulosporaceae.
Figure 4. Redundancy ordination (RDA) results (A, soil; B, root) of the characteristics of the AM fungi community and aspects and soil physicochemical properties in Tc natural forest. Abbreviatons are as follows: GSP, the proportion of Glomeromycota sequences; NTS, the number of total sequences; NGS, the number of Glomeromycota sequences; Div, Diversisporaceae; Gig, Gigasporaceae; Arc, Archaeosporaceae; Unc, unclassified fungi; Glo, Glomeraceae; Pac, Pacisporaceae; Aca, Acaulosporaceae.
Forests 13 02100 g004aForests 13 02100 g004b
Table 1. PCR amplification system used in Illumina PE250 high-throughput sequencing.
Table 1. PCR amplification system used in Illumina PE250 high-throughput sequencing.
PCR Amplification System (20 μL)
5 × FastPfu Buffer4 μL
2.5 mM dNTPs2 μL
Forward Primer (5 μM)0.8 μL
Reverse Primer (5 μM)0.8 μL
FastPfu Polymerase0.4 μL
Template DNA10 μL
Add ddH2O2 μL
Table 2. Results of variance analysis for the effects of aspect and species on the variables.
Table 2. Results of variance analysis for the effects of aspect and species on the variables.
VariablesAspectSpeciesAspect × Species
FpFpFp
TN4.5190.0450.363NS1.074NS
TK0.117NS2.8760.0232.5550.039
pH0.694NS2.291NS1.147NS
SOM0.208NS1.107NS0.099NS
TP3.7010.0441.063NS0.449NS
AP9.1500.0054.1020.0042.6890.031
S-Sequence3.332NS0.634NS0.846NS
S-Glomeromycota sequence0.465NS1.092NS2.3730.049
S-Glomeromycota sequence proportion0.577NS0.872NS1.314NS
S-OTUs0.203NS5.6830.0001.033NS
S-Chao2.716NS1.758NS0.184NS
S-Shannon2.576NS1.964NS0.122NS
S-Simpson0.672NS2.097NS0.226NS
S-Coverage0.065NS2.7640.0260.551NS
R-Sequence6.4720.0151.969NS0.792NS
R-Glomeromycota sequence3.986NS0.676NS0.800NS
R-Glomeromycota sequence proportion10.1160.0031.664NS1.109NS
R-OTUs3.725NS1.536NS0.171NS
R-Chao0.462NS5.1000.0011.134NS
R-Shannon0.426NS1.961NS0.989NS
R-Simpson0.179NS0.839NS0.744NS
R-Coverage3.3230.0315.0030.0090.720NS
Abbreviations are as follows: S, soil; R, root; NS means that there is no significant difference (p > 0.05).
Table 3. Basic sequence information of soil AM fungal communities in Tc natural forest in different aspects.
Table 3. Basic sequence information of soil AM fungal communities in Tc natural forest in different aspects.
SpeciesAspectsThe Number of SequenceThe Number of Glomeromycota SequenceGlomeromycota Sequence Proportion (%)
Toona ciliatesSouth50,890 (±2465) a15,143 (±241) a29.86 (±1.08) a
North40,187 (±3655) A12,062 (±2219) A *30.63 (±6.19) A
Padus buergerianaSouth40,673 (±2257) a12,076 (±1251) b29.46 (±1.78) a
North37,203 (±2745) A14,939 (±387) A40.93 (±3.59) A *
Maesa japonicaSouth44,598 (±4338) a14,748 (±345) a33.66 (±3.08) a
North39,575 (±7654) A15,104 (±200) A39.75 (±8.19) A
Melliodendron xylocarpumSouth49,182 (±5272) a14,960 (±224) a31.13 (±3.39) a
North35,266 (±7225) A15,320 (±96) A51.71 (±13.96) A
Ilex chinensisSouth49,826 (±3474) a14,698 (±302) a30.00 (±2.45) a
North40,715 (±5318) A15,044 (±418) A37.91 (±3.75) A
Alniphyllum fortuneiSouth42,443 (±5429) a14,584 (±402) a36.43 (±5.48) a
North48,785 (±5089) A13,194 (±1122) A28.59 (±5.07) A
Mallotus japonicusSouth39,661 (±12,721) a15,486 (±65) a53.09 (±22.60) a
North38,330 (±3760) A13,652 (±405) A *36.42 (±2.97) A
Note: Data are exhibited as the mean ± standard error. Data with different capital letters means significant difference in NF at p < 0.05. Data with different lowercase letters means significant differences in SF at p < 0.05; * means that there is significant difference between SF and NF.
Table 4. Basic sequence information of root AM fungal communities in Tc natural forest in different aspects.
Table 4. Basic sequence information of root AM fungal communities in Tc natural forest in different aspects.
SpeciesAspectsThe Number of SequenceThe Number of Glomeromycota SequenceGlomeromycota Sequence Proportion (%)
Toona ciliatesSouth49,755 (±4673) a30,186 (±70) a62.28 (±5.72) a
North50,058 (±3202) AB30,068 (±122) A60.88 (±4.28) AB
Padus buergerianaSouth50,257 (±2626) a29,920 (±157) a60.06 (±3.35) a
North52,982 (±2217) A30,126 (±89) A57.16 (±2.40) B
Maesa japonicaSouth53,405 (±655) a29,985 (±136) a56.18 (±0.91) a
North45,872 (±5774) AB30,119 (±97) A68.95 (±8.83) AB
Melliodendron xylocarpumSouth48,572 (±4766) a28,840 (±1099) a61.66 (±7.98) a
North42,850 (±5793) AB30,185 (±30) A74.13 (±9.20) AB *
Ilex chinensisSouth53,724 (±1835) a29,754 (±69) a55.59 (±2.00) a
North45,667 (±3392) AB30,200 (±83) A *67.20 (±4.82) AB
Alniphyllum fortuneiSouth48,351 (±1445) a28,635 (±1425) a59.30 (±3.02) a
North39,375 (±3677) AB30,124 (±118) A78.25 (±6.26) AB
Mallotus japonicusSouth45,732 (±1246) a29,803 (±302) a65.31 (±1.88) a
North38,492 (±4235) B30,015 (±166) A80.64 (±8.23) A
Note: Data are exhibited as the mean ± standard error. Data with different capital letters means significant difference in NF at p < 0.05. Data with different lowercase letters means significant differences in SF at p < 0.05; * means that there is significant difference between SF and NF.
Table 5. Relative OTU abundance of AM fungal communities on different aspects in Tc natural forest.
Table 5. Relative OTU abundance of AM fungal communities on different aspects in Tc natural forest.
SpeciesSouthNorth
RootSoilRootSoil
Toona ciliata15.1% a13.7% bc11.3% B16.4% AB *
Maesa japonica21.7% a17.9% a23.0% A18.6% A
Padus buergeriana12.3% a16.0% ab14.3% AB16.2% ABC
Melliodendron xylocarpum14.0% a10.6% c16.6% AB12.5% BC
Alniphyllum fortunei12.0% a13.9% bc13.6% B11.0% C
Ilex chinensis9.8% a13.1% bc8.3% B12.5% BC
Mallotus japonicus15.1% a14.8% ab12.9% B12.8% BC
Note: Data with different capital letters means significant difference in NF at p < 0.05. Data with different lowercase letters means significant difference in SF at p < 0.05; * means that there is significant difference between SF and NF.
Table 6. The diversity of soil AM fungal communities in Tc natural forest in different aspects.
Table 6. The diversity of soil AM fungal communities in Tc natural forest in different aspects.
SpeciesAspectsChaoShannonSimpsonCoverage
Toona ciliatesSouth421 (±156) a3.67 (±0.67) ab0.07 (±0.05) a0.9986 (±0.0005) ab
North236 (±73) B3.03 (±0.36) A0.10 (±0.03) A0.9990 (±0.0004) AB
Padus buergerianaSouth356 (±139) a3.81 (±0.36) ab0.05 (±0.02) a0.9985 (±0.0006) ab
North312 (±65) AB3.28 (±0.46) A0.11 (±0.05) A0.9985 (±0.0004) AB
Maesa japonicaSouth630 (±86) a4.35 (±0.10) a0.03 (±0.00) a0.9972 (±0.0000) b
North505 (±98) A3.91 (±0.07) A0.04 (±0.01) A0.9971 (±0.0014) B
Melliodendron xylocarpumSouth367 (±86) a3.85 (±0.10) ab0.04 (±0.01) a0.9989 (±0.0003) a
North353 (±93) AB3.56 (±0.31) A0.06 (±0.01) A0.9977 (±0.0010) AB
Ilex chinensisSouth272 (±170) a2.80 (±0.75) b0.18 (±0.10) a0.9992 (±0.0005) a
North161 (±26) B2.72 (±0.48) A0.15 (±0.06) A0.9995 (±0.0001) A
Alniphyllum fortuneiSouth337 (±72) a3.70 (±0.18) ab0.04 (±0.01) a0.9987 (±0.0003) a
North301 (±37) AB3.41 (±0.14) A0.08 (±0.02) A0.9989 (±0.0001) AB
Mallotus japonicusSouth382 (±83) a3.85 (±0.10) ab0.04 (±0.01) a0.9985 (±0.0002) ab
North277 (±45) B3.68 (±0.25) A0.05 (±0.01) A0.9985 (±0.0003) AB
Note: Data are exhibited as mean ± standard error. Data with different capital letters means significant difference in NF at p < 0.05. Data with different lowercase letters means significant difference in SF at p < 0.05.
Table 7. The diversity of root AM fungal communities in Tc natural forest in different aspects.
Table 7. The diversity of root AM fungal communities in Tc natural forest in different aspects.
SpeciesAspectsChaoShannonSimpsonCoverage
Toona ciliatesSouth409 (±52) bc3.52 (±0.29) a0.08 (±0.02) a0.9988 (±0.0002) ab
North466 (±63) AB3.93 (±0.05) A0.05 (±0.01) A0.9983 (±0.0002) A
Padus buergerianaSouth471 (±18) ab3.88 (±0.08) a0.05 (±0.01) a0.9985 (±0.0002) ab
North463 (±40) AB3.80 (±0.18) A0.05 (±0.01) A0.9984 (±0.0002) A
Maesa japonicaSouth534 (±34) a4.11 (±0.17) a0.04 (±0.01) a0.9983 (±0.0001) b
North517 (±59) A3.91 (±0.15) A0.05 (±0.01) A0.9980 (±0.0002) A
Melliodendron xylocarpumSouth308 (±31) c3.70 (±0.20) a0.05 (±0.01) a0.9990 (±0.0002) a
North358 (±30) B3.79 (±0.28) A0.07 (±0.04) A0.9986 (±0.0002) A
Ilex chinensisSouth387 (±31) bc3.63 (±0.13) a0.07 (±0.02) a0.9988 (±0.0002) ab
North374 (±25) B3.47 (±0.15) A0.07 (±0.01) A0.9985 (±0.0001) A
Alniphyllum fortuneiSouth448 (±32) ab3.57 (±0.11) a0.06 (±0.01) a0.9983 (±0.0002) b
North335 (±35) B3.44 (±0.17) A0.05 (±0.01) A0.9983 (±0.0002) A
Mallotus japonicusSouth455 (±37) ab3.83 (±0.18) a0.05 (±0.01) a0.9983 (±0.0002) b
North399 (±35) AB3.47 (±0.09) A0.07 (±0.01) A0.9979 (±0.0001) A *
Note: Data are exhibited as mean ± standard error. Data with different capital letters means significant difference in NF at p < 0.05. Data with different lowercase letters means significant difference in SF at p < 0.05. * means that there is significant difference between SF and NF.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Q.; Pan, J.; Ke, Y.; Yu, S.; Murray, P.J.; Luo, T.; Zhang, L.; Liu, W. Impact of Aspect on Arbuscular Mycorrhizal Fungal Diversity and Community Composition in a Natural Toona ciliata var. pubescens Forest in Subtropical China. Forests 2022, 13, 2100. https://doi.org/10.3390/f13122100

AMA Style

Wang Q, Pan J, Ke Y, Yu S, Murray PJ, Luo T, Zhang L, Liu W. Impact of Aspect on Arbuscular Mycorrhizal Fungal Diversity and Community Composition in a Natural Toona ciliata var. pubescens Forest in Subtropical China. Forests. 2022; 13(12):2100. https://doi.org/10.3390/f13122100

Chicago/Turabian Style

Wang, Qiong, Jianfeng Pan, Yunfang Ke, Suqin Yu, Philip J. Murray, Tianyu Luo, Lu Zhang, and Wei Liu. 2022. "Impact of Aspect on Arbuscular Mycorrhizal Fungal Diversity and Community Composition in a Natural Toona ciliata var. pubescens Forest in Subtropical China" Forests 13, no. 12: 2100. https://doi.org/10.3390/f13122100

APA Style

Wang, Q., Pan, J., Ke, Y., Yu, S., Murray, P. J., Luo, T., Zhang, L., & Liu, W. (2022). Impact of Aspect on Arbuscular Mycorrhizal Fungal Diversity and Community Composition in a Natural Toona ciliata var. pubescens Forest in Subtropical China. Forests, 13(12), 2100. https://doi.org/10.3390/f13122100

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