Mycorrhizal fungal community structure in tropical humid soils under fallow and cropping conditions

Little is known to what extent soil biota, in particular, the mycorrhizae are altered through different fallow durations/types in tropical soils. We found that soil-N, -C, -Al, -K and -Ca contents significantly differed due to the fallow durations/types. Subsequently, the effects of fallow types and soil depths on the diversity, species richness and community structure of arbuscular mycorrhizal (AM) fungi were examined. A higher AM species richness was identified in the cropping than in forest fallow fields, suggesting a positive cropping feedback on the AM community composition. Distribution of the AM species was positively related to soil properties, specifically soil-pH, and soil-Pi, -Ca and -Mg contents. The soil properties conjointly accounted for 78.5% of explained variation in the AM community composition, signifying that the main factors altering the community structure under different fallow and cropping systems were the soil properties. Among the soil chemical characteristics, the soil-pH disclosed a significant explained variation in the AM community composition in the topsoil layer under the short fallow. Structural modeling equation to understand multiple predictive pathways that connect soil properties, fallow practices and AM community structures indicated that soil-C, -N and -Ca contents were highlighted as important factors influencing the AM community compositions.

further exacerbated in the actual cultivating systems of the small-scale farms, which can rarely afford organic/ inorganic inputs to support crop growth 5,14 . Therefore, soil fertility regeneration through fallow process is no longer feasible; hence, threatening the sustainability of small-scale agriculture 6,12,15 . From the soil microbiota perspective, microbes play significant roles in the productivity of lands, and are responsible for diverse ecosystem processes and their functioning 16,17 . For example, the great part of nitrogen (N) is recycled via the association between roots of legume plants and rhizobia through the biological atmospheric nitrogen (N 2 ) fixation process [18][19][20] . The soil-turning processes by soil invertebrates contribute to nutrient cycling and accumulation, affecting plant growth 1,21 . The beneficial associations between plant roots and arbuscular mycorrhizal (AM) fungi support uptake and translocation of immobile nutrients, largely N, inorganic phosphate (Pi), carbon (C) and water to host plants in exchange for fixed C source to AM fungi 22,23 . Because of their novel functions in mineral uptake and alleviation of water deprivation-related stresses, the contribution of AM fungi to smallholder agriculture in tropical areas has been recently considered the most important compared with that of other groups of the soil biota 24 . The extent of AM benefits to plants largely depends on their community structure [25][26][27] , crop species and/or genotype 28 , soil characteristics 29,30 , climatic and geographic context, and the interactions among the above-listed factors 25,27,31 . AM diversity and species richness are in turn largely influenced by the composition of associated plant communities, soil geography, and various agricultural management practices, mainly fertilizer applications, crop rotation, monocropping and tillage practices with selection of soil depths being an important part of tillage practices 32 . However, how the AM composition is reshaped under different soil horizons; this information remains markedly limited, especially under the farming agroecosystem situations in the tropical forest regions, where the cropping activities are almost restricted to the plow layer of soils with direct/ indirect consequences on soil chemical and biota changes 33,34 . Furthermore, there has been only scant information about how diverse groups of soil biota are altered by soil depths and their feedback on soil physical-chemical properties 32,35 .
As a result, investigations of AM communities and their importance in tropical ecosystems have been encouraged, especially in the farming systems where the smallholder farmers rarely use fertilizer inputs for their crops 36 . Consequences of the different fallow and cropping conditions imposed by shifting cultivation practices on microbial community and their feedback on plant growth promotion remain largely underinvestigated 5,13 . Increasing evidence depicts ample depletion of nutrients (e.g. soil-Pi and -N), following cropping with inadequate replenishment of the sources 37 . Therefore, the objectives of the present study were (i) to assess the impacts of the fallows, mainly long-term forest and short-term fallows, as well as cropping without fallow on the AM community composition, (ii) to examine how AM community composition was evolved and/or reshaped at different soil depths, and (iii) to assess the effects of soil characteristics, namely soil-pH, and soil-N, -C, -Pi, -Ca, -Mg, -Al and -K contents, and their relevant interactions on the AM community structure, diversity and species richness. Through the knowledge generated, we could gain a better insight into the factors driving the AM community composition and structure in the context of farming systems under the humid tropical soil conditions.
In each of the sampling sites, cropping fields and two kinds of fallow fields were selected, with four replicates each being sampled in two depths (0-10 and 10-20 cm) for assessment of the AM community composition and structure. The two fallow fields were the long-term forest fallow field of 20-year-old and the short-term fallow field dominated by Chromoleana odorata, an invasive weed species of 4-to 6-year-old. The cropping fields were dominated with planted food crops, mainly groundnut (Arachis hypogaea) and cassava (Manihot esculenta) grown in association with cocoyam (Colocasia esculenta), banana (Musa sp.), maize (Zea mays), yam (Discorea rotundata) and a number of leafy vegetables (Amaranthus cruentus, Corchorus olitorius, Solanum scabrum and Vernonia calvoana). The farmers traditionally managed their fields from land preparation, planting/sowing, weeding to harvesting. The soil sampling activities were conducted during vegetative growth stage of food crops, which corresponded to approximately 48-56 days after planting/sowing.
To carry out the sampling of soils from fields, sixteen core soil samples were collected from 0-10 and 10-20 cm depths, respectively, in four delimited 10 m × 10 m plots at each site. All core samples were bulked and transferred to the laboratory. In total, 72 composite soil samples were collected from all three sites and stored at 4 °C for subsequent analyses. As the first step, the soils were sieved through a 4-mm mesh, and then divided into two pools. The first pool was used to establish the bioassay culture, and the second pool was sieved to pass a 1-mm mesh and was then used for chemical analyses (e.g. measurements of soil-pH, as well as contents of total soil-C and -N, and available soil-Pi and other elements).

Chemical analyses.
After the soil samples were shaken in water for 16 h, soil-pH was determined in the soil suspension (1:2.5, v:v) with the help of a Corning 125 pH meter (Corning, Life Sciences, Amsterdam, Netherlands). Four mg of each sample collected from the respective sites were pulverized in a ball mill (MM200, Retsch, Germany). The total soil-N and -C contents in prepared, dried samples were measured by using an elemental analyser (Flash EA 1112, Thermo Electron, Ecublens, Switzerland). For determination of the soil-Ca, -Mg, -Al and -K contents, 500 mg of soil were incinerated in a muffle furnace at 550 °C for 8 h. Subsequently, the ashes were dissolved in 2 mL of concentrated HNO 3 at 180 °C and were then made up to 50 mL with ultrapure water. The soil-Ca, -Mg, -Al and -K contents were measured in the 20-fold diluted extracts using an inductively coupled plasma mass spectrometry (ICP-MS, Agilent 7500c, Agilent, Basel, Switzerland). To determine the inorganic Pi content in soil, the Bray-I method was used. Briefly, 30 mL of Bray-I extractant were added to 3 g of air-dried soil sample. The mixture (1:10, soil:solution ratio) was shaken for 5 min and filtered, and the soil-Pi content was measured by colorimetric change 40 . Bioassay. The experiment was designed in a factorial arrangement with the following factors tested: (i) factor 1, fallow types (F) with three levels: long-term forest fallow, short fallow and cropping fields; (ii) factor 2, soil origins (S) with three levels: Tko, Tku and Rku soils; and (iii) factor 3, sampling depths (D) with two levels: 0-10 and 10-20 cm. To establish the bioassay, 200 g fresh soil of each sample were mixed together with 800 g of autoclaved sand-quartz substrate (grain size between 0.7 and 1.2 mm), and 1 L capacity of the substrate mixture was used to fill up each pot. Subsequently, one-week-old pre-germinated seedlings of leek (Allium porrum) and red clover (Trifolium pratense) were manually transplanted to the pots. The leek and red clover plants were grown for 6 months in growth chamber PGC20 (Conviron, Winnipeg, Canada), and then they were replaced by pre-germinated alfalfa (Medicago sativa) and ryegrass (Lolium multiflorum) seedlings that were subsequently grown for 6 months in a growth chamber. The following growth chamber conditions were maintained before the soils were sampled for AM microscopic observations and molecular analyses: 16 Microscopic observations. Spore extraction of the AM species was carried out using the modified wet sieving method followed by sucrose-gradient centrifugation as previously described 42 . Briefly, 25 g of wet substrate mixture were washed thoroughly with tap water through a 500-μm mesh, and then the suspension was collected using a 40-μm mesh sieve and transferred into a 50-mL centrifuge tube. Ten mL of 2.5 M sucrose solution were added followed by a centrifugation at 2000 × g for 2 min. After centrifugation, the supernatant was poured into a 40-μm mesh sieve and quickly rinsed with tap water. Turgid (viable) individuals were counted and grouped according to the species morphology. Permanent slides were prepared for each species using either polyvinyl alcohol alone or mixed with Melzer's solution (1:1), and were then observed under the Olympus AX70 microscope (Olympus Optical, Tokyo, Japan). Different AM species recorded were compared with the recent nomenclature of AM fungi as previously described 43 . Soil moisture content was determined after a subsample of fresh soil from each pot was oven-dried at 105 °C for 48 h, and the dry weight was recorded. Spore density, expressed as the number of AM spores per 25 g of fresh substrate, was adjusted to the dry matter of soil using the soil moisture content values. The Shannon-Weaver diversity index (H) was calculated for different F, S and D using the following equation Eq. 1: where pi was the spore abundance of the i th species among all N species identified in a sampling depth of a soil origin under a respective fallow or cropping field. The species abundance distribution (SAD) of each identified AM species 44,45 was determined based on the equation Eq. 2: = × SAD (Total spore number of each AM species 100) /(Total number of all identified AM species) The AM species richness, which indicates the number of the AM species found in soil, was determined as previously described 46 .

Molecular analysis.
To complement the morphological identification, sequencing of the 28S ribosomal DNA (rDNA) of the respective AM isolates observed under the microscope was conducted according to a previously published method 42 . Briefly, DNA was recovered from fresh and juvenile spores of each AM species, and amplified by the polymerase chain reaction (PCR) using the LR1 (5′-GCATATCAATAAGCGGAGGA-3′) and LR2 (5′-GTCGTTTAAAGCCATTACGTC-3′) primers 47 . The PCR reactions were performed in an automated thermal cycler (Biometra T-Gradient, Leusden, Netherlands) with an initial denaturation step at 94 °C for 3 min, amplification with 36 cycles (denaturation for 15 s at 94 °C, annealing for 90 s at 50 °C, elongation for 90 s at 72 °C), and final elongation for 10 min at 72 °C. The PCR products were visualized on 0.8% agarose gel with ethidium bromide under UV light to confirm their integrity. The PCR products were then purified using a gel extraction kit (Qiagen AG, Hombrechtikon, Switzerland), following the manufacturer's procedure, and cloned into the pGEM-T plasmid (Promega, Dübendorf, Switzerland). Sequencing of the cloned PCR products was carried out using a Sanger DNA sequencer (Microsynth, Balgach, Switzerland), and each sequence was manually edited and blasted against GenBank for homology search using the BLASTN algorithm. Multiple sequence alignment was conducted using BioEdit (v.7.2.1) 48 . Phylogenetic analysis was performed by the Neighbour-joining method with a bootstrap value of 1000 using the MEGA5 software 49 . We chose the AM species from the GenBank to represent as broadly as possible the phylogenetic diversity within the phylum Glomeromycota and using Morchella esculenta and Fusarium sambucinum as outgroup species. Data analysis. A three-way analysis of variance (ANOVA) was performed using the General Linear Procedure 'Proc GLM' in the Statistical Analysis System (SAS) software (v. 9.2) 50 to test the effects of F, S, D, and their interactions on the soil parameters and the AM community structure and diversity. On the basis of the significant Fisher's (f) values, a Least-Means test (Least-Means/PDIFF option) was performed to separate means among the F, S and D. This was carried out exclusively when the Student's t-test from the ANOVA resulted in a significant p-value (p < 0.05). Levels of significance are given by 'ns' (not significant, p > 0.05), *p < 0.05, **p < 0.01 and ***p < 0.001. A principal component analysis (PCA) was conducted to explore the interrelated effects of soil characteristic data (e.g. soil-pH, and soil-N, -C, -Pi, -Al, -K, -Ca and -Mg contents) obtained in the two depths of the forest, short fallow and cropping fields. A redundancy analysis (RDA)-based variation-partitioning method was used to extract the effects of explanatory variables (i.e. soil-pH, and soil-N, -C, -Pi, -Al, -K, -Ca and -Mg contents), and factorial predictors (i.e: F, S and D) on the AM community composition using the 'rda' function in the 'vegan' library of the R program. Prior to RDA, a detrended correspondence analysis (DCA) was used to extract the axes of maximum variation in the AM species composition, and to identify the gradient lengths from original data. Biplots of the DCA and RDA were displayed using the CANOCO software (v.5.2) 51 . The response curves of individual AM species against the DCA axis 1 (DCA1) were performed, and the model fitted for the respective AM species was presented. The significant levels for both canonical axes were determined, and their f-values were reported following the Monte Carlo permutation test with 499 permutations 52 . Furthermore, the correlations between individual soil characteristics (SC) and different AM species were analyzed by generalized additive model using the Gaussian distribution and log (x + 1)-transformed SAD values. When the effects of explanatory variables were significantly correlated with the AM species composition, a variation-partitioning analysis was performed adopting the approach previously described 53 . A prior stepwise selection of explanatory variables in each of the two groups of predictors (quantitative and factor estimators) was carried out with candidate predictor inclusion based on the significant values 54 . Venn diagram of the variation-partitioning results was then created. The adjusted p-values for multiple comparison tests reporting statistically significant influences among the analyzed soil properties (i.e. soil-pH, and soil-N, -C, -Pi, -Al, -K, -Ca and -Mg contents) on the AM species composition were determined based on the Bonferroni correction test.
The structural equation modeling (SEM) was used to estimate the multiple cause-effect relationships between the fallow practices, SC and AM composition 32 . The SEM model allows the specification of multiple predictive pathways among model variables to account for their influence on each other based on the two sub-model approach: the structural model and measurement model. The structural model was built from the analyzed soil parameters, results of the RDA-based variation-partitioning and ANOVA analyses, taking into account the published literature on the AM studies. An explanatory stepwise procedure was conducted through the use of modification indices to increase the fit of the models by relaxing imposed constraints from the measured vairables. The final models were reached after their validation using the maximum criteria required for the chi-square (χ 2 ), the standardized root mean square residual (SRMR), and the comparative fit index (CFI) based on the maximum-likelihood analysis. The path coefficients indicated the degree (strong/weak) of relationship, while the thick and thin arrows indicated direct or indirect effects from the measured variables that were well fitted into the models.

Results
Variations in soil properties across fallow types, soil origins and sampling depths. Results from the one-way ANOVA showed significant different effects due to the fallow types (F) on six of the eight examined soil chemical properties, with the exception of the soil-Mg and -Pi contents ( Table 1). The soil-pH values ranged from 4.3 to 7.2, and were acidic regarding the Tko and Tku, and were neutral with reference to the Rku soils (Table 2). Furthermore, the soil-pH exhibited higher values in the 0-10 cm depth of the Tko soil under the cropping and short fallow as compared with that of the forest fallow. Both the cropping and short fallow soils displayed greater soil-pH values relative to that of the forest fallow in both top and lower layers (0-10 and 10-20 cm depths) of the Tku soil. We also found that the cropping displayed the highest soil-pH value in comparison with the forest and short fallows in the 0-10 cm sampling layer of the Rku soil. In addition, the cropping fields recorded greater soil-C, -N and -Al contents than the forest and short fallows in both the top and subsoil layers of the Tko soil. We also noticed that the cropping soils displayed higher soil-C, -N, -Pi, -Mg, -Al and -Ca contents than the forest and short fallows in both two sampling depths of the Rku soil. Additionally, the cropping and short fallow exhibited greater soil-C contents than the forrest fallow in both sampling depths of the Rku soil. With respect to the Tku, the cropping fields showed an increasing tendency in soil-Mg and -Al contents as compared with the other two fallows in the 0-10 cm soil layer. We also found that the cropping and short fallow fields displayed the lowest soil-N contents as compared with that of the forest fallow fields in the 10-20 cm depth of Tku soil ( Table 2).
With regard to the sampling depths (D), the ANOVA exhibited significant f-values for the majority of analyzed soil characteristics (SC), except the soil-Al and -K contents (Table 1); with the 0-10 cm sampling depth displaying the highest values of the soil-pH, and soil-C, -N and -Ca contents, irrespective of the F and soil origins (S) ( Table 2). Furthermore, the Rku soil generally displayed greater values of soil-pH than the Tko and Tku soils in both soil depths (Table 2). Among the interactions tested, the F × S interaction resulted in significantly different effects on the soil-pH, and soil-Ca and -Pi contents (Table 1). Similarly, the S × D interaction displayed statistically significant f-values with respect to the soil-pH, and soil-Mg, -Al, -K, -Ca and -Pi contents, while the F × D interaction showed significantly different influences on the soil-N, -C and -Pi levels (Table 1). Finally, the F × S × D interaction exhibited significantly different effects on the soil-pH, and soil-K, -Ca and -Pi contents following a 3-way ANOVA ( Table 1).
The PCA investigating how the SC were interrelated under the different F and D allowed us to identify two major components accounting for the total variance explained (TVE) with cumulative values ranging from 66 Table 1 Table 1). Finally, we noticed positive and high loading scores for PC1 from the soil-Mg variable in the 0-10 cm soil depth under the short fallow and cropping conditions (Fig. 1c,e; Supplementary Table 1). As for the PC2, the soil-N and -C variables exhibited high and positive loadings in the 0-10 cm soil depth of the forrest fallow, while high and negative loadings in the 0-10 cm soil layer under the short fallow and cropping conditions (Fig. 1a,c,e; Supplementary Table 1). The loadings from the soil-Al variable into the direction of PC2 indicated high and positive scores in the 0-10 cm depth of the forest fallow and in the 10-20 cm depth of the forest and short fallows ( Fig. 1a,b,d; Supplementary Table 1). In addition, the soil-Mg content displayed positive and high loading values in the 10-20 cm soil depth of all examined forrest fallow, short fallow and cropping practices in the direction to the PC2 (Fig. 1b,

Identification and classification of the AM species in soils after the fallows and cropping. The
28S rDNA fragments were amplified from the spores of different AM species, which were detected under a microscope, and were sequenced. Subsequently, the sequence data (accession numbers provided in Supplementary  Fig. 1) were blasted with available sequences in GenBank for homology search using the BLASTN. A phylogenic tree was then constructed using the sequence data of the 28S rDNA fragments to determine the phylogenic relationship of the individually identified AM species ( Supplementary Fig. 1). Up to 12 different AM species were detected from various fallow types, soil origins and sampling depths. Using the nomenclature of AM fungi previously reported 43 , the identified AM fungi could be classified into the Glomeraceae, Diversporaceae, Acaulosporaceae, Archeosporaceae, Gigasporaceae, Ambisporaceae and Paraglomaraceae families ( Supplementary Fig. 1). Impact of fallow types, soil origins and sampling depths on the AM species abundance distribution, richness, diversity and number of their spores. Results of the 3-way ANOVA testing the effects of fallow types (F), soil origins (S) and sampling depths (D) on the species abundance distribution (SAD) of the AM fungi are presented in Table 3. Significantly different effects of the F were detected on the SAD of four out of the 12 recorded AM species, namely Acaulospora scrobiculata, A. mellea, Rhizophagus intraradices and Ambispora sp., while the effects of S on the SAD of AM species displayed significant f-values for A. mellea, Cetraspora pellucida, Racocetra tropicana and Glomus microaggregatum and Ambispora sp. (Table 3). Significantly different effects of the D were also observed on the SAD of five out of 12 detected AM species, including A. scrobiculata, A. mellea, Paraglomus occultum, Gigaspora decipiens and R. intraradices. Concerning their interaction effects, the F × S interaction revealed significant f-values for A. mellea and P. occultum, whereas F × D and S × D interactions showed significantly different effects on A. scrobiculata, C. pellucida and R. tropicana, and A. scrobiculata, A. mellea and A. paulineae, respectively. Additionally, the triple F × S × D interaction exhibited critical influence on SAD of six AM communities, namely A. scrobiculata, A. mellea, C. pellucida, R. tropicana, A. gerdemannii and Ambispora sp. (Table 3).
The 3-way ANOVA also demonstrated significantly different effects of the F on the AM species richness ( Table 3). The cropping fields generally recorded greater AM species richness (α-diversity) as compared with the forest fallow in both two sampling depths of the Tko and Rku soils, according to a one-way ANOVA  ( Fig. 2a-c). Furthermore, the S also displayed significantly different influences on the AM species richness (Table 3). Moreover, the Fisher's test indicated that the effects of the D on the AM species richness were also highly significant (Table 3), with the 10-20 cm depth of the Rku soil exhibiting greater values of AM species richness than the 0-10 cm depth under all fallow and cropping conditions (Fig. 2c). Overall, the AM species richness ranged from 1.5 to 7.5; with the 0-10 cm sampling depth of the Rku under the forest fallow displaying the lowest, and the 0-10 cm sampling depth of the cropping field of Tku soil disclosing the highest value of the AM species richness (Fig. 2a-c). Among the tested interactions, F × S, S × D and F × S × D interactions showed significantly different influences on the AM species richness (Table 3).
With respect to the AM species diversity, the Shannon-Weaver diversity index ranged from 0.2 to 12.3 ( Fig. 2d-f). The 10-20 cm depth of Tko soil displayed the lowest, and the same sampling depth of the Rku soil disclosed the highest Shannon-Weaver diversity value under the short fallow (Fig. 2d,f). The 3-way ANOVA testing the effect of the F on the Shannon-Weaver diversity index noted a significant f-value (Table 3). For instance; within the Tko soil, the forest fallow and cropping fields showed the highest diversity index in 0-10 and 10-20 cm sampling depth, respectively (Fig. 2d). The 3-way ANOVA also revealed significantly different effects of the S on the Shannon-Weaver diversity index (Table 3). In addition, the effects of D on the Shannon-Weaver diversity index were greatly and significantly different according to the 3-way ANOVA (Table 3), displaying remarkable differences between the 0-10 and the 10-20 cm depths (Fig. 2d-f). For instance, the 10-20 cm depth of the Rku soil recorded higher Shannon-Weaver diversity values than the 0-10 cm depth of the same soil type in all forest fallow, short fallow and cropping fields (Fig. 2f). Two out of the four tested interactions, namely the S × D and F × S × D interactions were significantly different in terms of the Shannon-Weaver diversity index ( Table 3). The effects of the F were significantly different on the AM spore numbers as shown by the 3-way ANOVA ( Table 3). The short fallow recorded remarkably higher AM spore numbers than the forest fallow and cropping fields with respect to the 10-20 cm sampling depth of the Tko and Rku soils (Fig. 2g,i). The effects of D on the AM spore numbers were also significantly different (Table 3), with the 10-20 cm depth displaying higher AM spore numbers than the 0-10 cm depth with respect to the Rku soil (Fig. 2i). Additionally, greater AM spore numbers were observed in the 10-20 cm than in the 0-10 cm sampling depth of the Tko soil under the short fallow conditions (Fig. 2g). When comparing the AM spore numbers across the F, S and D, the values ranged from 49 to 1406 (Fig. 2g-i); with the 0-10 cm sampling depth of the Tku soil under the forest fallow displaying the lowest value, while the 10-20 cm depth of the Rku soil under the short fallow showing the highest value ( Fig. 2h-i). None of the four tested interactions among the F, S and D was significant in terms of the AM spore numbers (Table 3).  Fig. 2). The DCA interpreting original variations extracted from the soil parameter predictors indicated that the longest DCA1 axes ranged from 39.8 to 49.3% of the AM species variabilities for the two sampling depths under the three forest fallow, short fallow and cropping conditions ( Supplementary Fig. 3). Next, we performed a RDA to explain the relationships between the AM community structure and various soil characteristics (SC) under various combinations of D and F (Fig. 3). In the 0-10 cm soil depth of the forest fallow, the soil-pH and soil-Ca content exhibited significant influences on the AM community composition (Supplementary Table 2), whereas three (e.g. soil-Mg, -Al and -Ca contents) out of eight analyzed SC displayed significant roles in explaining the AM community composition in the 10-20 cm subsoil layer of the forest fallow (Supplementary Table 2). We further noticed remarkable contributions from soil-pH, -Mg and -K variables to the AM variability in the 0-10 cm depth under the short fallow (Supplementary Table 2), while none of the SC used as the predictors showed significant effect on the AM community composition in the 10-20 cm soil depth under the short fallow (Supplementary Table 2). The soil-Mg and -Ca contents displayed significant contributions to the variation in the AM community structure in the 0-10 cm soil depth of the cropping fields, while the soil-C and soil-Pi variables revealed their significant proportion in the AM composition variation in the 10-20 cm soil depth of the cropping fields (Supplementary Table 2). It should be noticed that among the p-adjusted values reporting the differential effects among the different soil parameters (soil-pH, -N, -C, -Mg, -Al, -K, -Ca or -Pi) on the AM species variabilities, only that of the soil-pH variable with respect to the 0-10 cm soil depth under the short fallow conditions was significant (Supplementary Table 2). Table 4 depicts the correlation patterns between the individual soil property and the species abundance distribution (SAD) of each of the AM species under different soil depths, and fallow and cropping conditions. In the 0-10 cm soil depth of the forrest fallow, the relationships between the soil-pH and the SAD of Ambispora sp.; between the soil-Mg content and the SADs of Ambispora sp. and C. pellucida; and between the soil-Al content and the SAD of A. scrobicultata were significant. Furthermore, the soil-Ca content exhibited positive relationships with the SADs of A. scrobiculata, A. paulineae and R. tropicana in the 0-10 cm soil depth of the forest fallow. In the 10-20 cm soil layer of the forest fallow, both soil-pH and soil-Ca content showed significant correlations with the SADs of C. pellucida and R. tropicana (Table 4).

Relationships between the soil properties and
As for the short fallow conditions, the SADs of six out 12 AM species, namely A. scrobicultata, A. mellea, Ambispora sp., C. pellucida, R. tropicana and G. microagregatum were significantly related to the soil-pH in the    (Table 4). Similarly, the SADs of A. scrobiculata, A. mellea and Ambispora sp. and the soil-Mg content were significantly interrelated; while that of G. decipiens and the soil-K content, and that of R. intraradices and the soil-Ca content showed a positive correlation in the 0-10 cm soil depth of the short fallow. In comparing the relationships between the soil properties and the SAD of the AM fungi in the 10-20 cm soil depth under the short fallow conditions, we found that both soil-pH and soil-Mg content exhibited significant correlations with the SADs of A. scrobiculata and G. microagregatum, while both soil-N and -C contents showed significant relationships to the SADs of A. paulineae and R. intraradices. Furthermore, significant relationships were also observed between soil-Pi content and the SAD of A. scrobiculata; soil-Al content and the SAD of G. margarita and; and soil-Ca content and the SAD of A. mellea in the 10-20 cm soil depth under the short fallow conditions (Table 4). With regard to the cropping fields, the soil-C content displayed significant relationships with the SADs of A. mellea, C. pellucida, R. tropicana and G. microagregatum in the 0-10 cm depth, while the soil-Mg content exhibited significant correlations with the SADs of A. scrobiculata and A. paulineae ( Table 4). The soil-Ca content was significantly related to the SADs of A. paulineae, P. ocultum and A. gerdermannii, whereas the soil-K content showed significant relationships to the SADs of A. gerdermannii and G. margarita in the 0-10 cm sampling depth of the cropping fields. Significant correlations were also observed between soil-Pi content and G. decipiens, and between soil-Al content and R. intraradices under the same conditions. After the cropping, in the 10-20 cm soil depth, the soil-Pi content showed significant correlation with the SAD of R. intraradices, while the soil-N content exhibited significant relationships with the SADs of A. paulineae, Ambispora sp. and G. margarita. Furthermore, we recorded significant relationships between the soil-C content and the SADs of Ambispora sp. and G. microaggregatum, soil-Al content and the SAD of P. occultum, and soil-K content and the SAD of R. intraradices in the 10-20 cm soil depth of the cropping fields (Table 4).   Identification of factors altering the AM community composition after the fallows and cropping.
Results of the RDA-based variation-partitioning approach dissecting the contribution of each of the predictors on the AM community composition revealed that the soil characteristics (SC), soil origins (S) and fallow types (F) were responsible for 78.5, 31.0 and 10.0% of the variations in the AM community composition, respectively (Fig. 4a). Among the SC, soil-pH, and soil-Ca, -Mg, -C and -N contents exhibited significant effects on the variations of the AM composition structure across the sampled locations (Supplementary Table 3). The p-adjusted values explained the significant variations for the soil-pH, and soil-Ca, -Mg, -C and -N contents among the SC (Supplementary Table 3). Among the various combinations of F, S and SC, the combination S × SC had the lowest effect on the variation of the AM composition structure as factorial predictor, followed by the F × S (Fig. 4a). In analyzing the influences of the factorial F, S and D predictors and their combinations, we found that the D, S and F were responsible for 61.9, 20.6 and 20.4% of the variations in the AM community structure, respectively, while their combinations F × S, F × D, S × D, F × S × D exhibited neglible effects on the AM community composition structure variability (Fig. 4b). These data indicated that the effects of the individual predictors were stronger than that of their combinations in controlling the AM species composition and distribution under the analyzed environmental conditions (Fig. 4a,b). Results of the SEM using the data of analyzed soil properties under the different F, S and D conditions as a single group supported our previous observations obtained from the RDA, and further indicated the different pathways and effects in which each predictor exercised its influence on the AM community composition (Fig. 4c). The final built models fitted well to our original data with a non-significant χ 2 of 52.3, a SRMR value < 0.08 and a CFI value > 0.9, indicating that the hypothesized measurement models properly fitted to the data. Soil properties included as latent variables in the models indicated their direct effects on the AM composition, with the measured soil-C, -N and -Ca variables significantly connecting the soil properties and AM composition (Fig. 4c). The fallow practices, when included in the model, also showed direct effects on soil characteristics, with the measured soil-C, -N and -Ca variables displaying strong relationships to the fallow practices on the basis of the observed path coefficients (Fig. 4c). On the other hand, we discovered that the fallow practices exhibited indirect influences on the AM composition (Fig. 4c). When the models were constructed separately using the datasets obtained from the two soil depths, we found that the fallow practices and soil characteristics disclosed direct effects on each other, and the soil-Ca variable exhibited relatively strong relationships to fallow practices under both sampling depths (Fig. 4c). The soil origins showed strong direct effects on the AM community compositions in the upper 0-10 cm soil depth only (Fig. 4c). With respect to the fallow types, they exhibited indirect weak effects on the AM composition under both sampling depths as observed from their path coefficients (Fig. 4c).

Discussion
Changes in soil characteristics by fallow types across soil origins and sampling depths. In this study, we first dissected the contributions of the fallow types (F), including forest fallow, short-fallow and cropping to the changes in chemical characteristics of three soil types (S; Tko, Tku and Rku) in the two sampling depths (D; 0-10 and 10-20 cm). Our results demonstrated a clear tendency that the examined F exhibited significantly different effects on the chemical properties of all three types of soils collected from the two sampling depths (Tables 1 and 2). Additionally, the S and D significantly impacted the analyzed soil parameters (Tables 1  and 2). Remarkable differences in soil-pH of different soil types have been reported earlier in the soils sampled in these areas of southern Cameroon 39 . Similarly, significant differences among the soil types regarding the available soil-Pi content were also previously noted, when the earthworm cast production was examined in different forest and short fallow fields at the same locations 39 . The deterioration in available soil-Pi content in the cropping fields as compared with the forest fallow fields was observed in the Tko soils for both sampling depths (Table 2), corroborating the hypothesis of nutrient deterioration due to cropping which is partly restored through long-term fallowing process under shifting cultivation system 5 . The PCA approach allowed us to summarize the different soil data into two orthogonal components from the various F, S and D conditions ( Fig. 1; Supplementary Table 1). The first PC component showed highly positive loading scores for soil-pH and -Ca in the 0-10 cm soil depth of the forest fallow and cropping conditions, and in the 10-20 cm soil depth under the short fallow and cropping conditions, irrespective of the S (Supplementary Table 1). It is worth mentioning that the Rku soil showed a closer relationship with soil-pH than Tko and Tku soils in both two soil depths under all three forest fallow, short fallow and cropping conditions (Fig. 1), which might be associated with the higher soil-Ca content observed in the Rku than in Tko and Tku soils ( Fig. 1; Table 2).    25,59,60 . Possible reasons of the relatively low AM diversity observed in the present study could be attributed to the use of the trap culture approach that could potentially inhibit the sporulation of rare AM species, and thus might underestimate the α-diversity 61 . We opted to use the trap culture method to recover fresh AM spores that were further used as isolates for AM inoculant production for the planned assays in the frame of this project. Another possible reason of the low AM diversity could be due to the incompatibilities between the primers used and the target regions, which might result in low matching of the original sequence and lack of possibility to detect some important rare AM species 62 . In addition, the limited availability of reference sequence data for various AM species could contribute to the reduction in recovery of more AM species. We further noticed that the species abundance distributions (SADs) of several AM species, mainly A. scrobiculata, A. mellea, Ambispora sp. and R. intraradices, were greatly influenced by the F (Table 3), with the cropping revealing the highest AM species richness (Fig. 2a-f). The effects of the cropping practices, such as tillage, The various pathways drawn from the models used the dataset as single group. Path coefficient values shown above to the pathways are for the soil depths and listed in their order (0-10 cm/10-20 cm). Only the significant path coefficient values (p < 0.05) from the models are presented. The thick and thin arrows indicate the direct and indirect effects, respectively. Dotted frame indicates the soil properties, namely soil-C, -N and -Ca contents that strongly influenced the model groups. Values in parentheses indicate path coefficient values for the main models using the dataset as single group. ns, not significant. crop rotation and other agricultural management approaches, on the AM community compositions have been extensively studied 5,25,32,63 . Cropping practices comprising of tillage, crop rotational systems, and/with organic manuring incorporation and inorganic fertilizer applications negatively affect the AM community composition, and lead to shifts in composition and reduction in the AM diversity 16,25,64 . In addition, the effects of the S were found to be significant on the SADs of several AM species, particularly A. mellea, Ambispora sp., C. pellucida, R. tropicana and G. microaggregatum (Table 3). This result is supported by a previous report that showed that S is among the determinant factors for the AM diversity in the Swiss agroecosystem 58 .
To explore the extent to which top and subsoil layers influence the AM species richness and diversity in a community under different fallow and cropping conditions, we conducted different hierarchical analyses. We found that the D remarkably impacted the AM species richness and Shannon-Weaver diversity index (Table 3). Interestingly, the subsoil layer (10-20 cm) displayed higher AM species richness and Shannon-Weaver diversity index than the topsoil layer (0-10 cm) under all fallow and cropping conditions of the Rku soils (Fig. 2c,f). Recent studies conducted under the temperate climate conditions have also shown different impacts of the topsoil and subsoil layers on the species richness and Shannon-Weaver diversity index of the AM fungi 32,65 . Factors driving the AM community compositions across the different fallow and cropping conditions, as well as the sampling soil layers remain to be determined, particularly in the tropical soils. Plausible hypotheses are linked to the differences in the composition of native plant species used, which play an important role in shaping the AM community composition in different soil depths as suggested by previous reports 32, 65 . The reason that the cropping practice revealed higher AM species richness than both the forest and short fallows in almost all combinations of D and S still remains unknown as well (Fig. 2a-c).
Factors influencing the AM community composition. Using the RDA-based variation-partitioning method and the SEM approach, we identified predictors of importance in structuring the AM community composition (Figs 3 and 4). Under different fallow type (F) and sampling depth (D) conditions, several members of the examined soil characteristics (SC) were shown to affect the AM community composition. The effect of soil-pH on the AM composition was particularly observed in the topsoil layer (0-10 cm) under the forest and short fallow conditions (Fig. 3a,c; Supplementary Table 2). Previous studies investigating the effects of soil-pH on the AM and bacterial community compositions suggested that low soil-pH conditions alter the microbe community structure 66,67 . Furthermore, the variation-partitioning method allowed us to define the contributions from different SC, F, soil origins (S) and D, as well as their combinations, to the AM composition variations (Fig. 4a,b). Among the SC, F and S, the SC and the F had the highest (78.5%) and lowest (10.0%) contribution, respectively, to the AM variation (Fig. 4a). A significant contribution of SC to the variation of AM community composition has also been reported in a previous investigation of soils on the North American Great Plains 32 . Furthermore, among the soil properties examined, the soil-pH, and soil-Ca, -Mg, -C, and -N contents showed significant effects on the AM composition structure (Supplementary Table 3). For instance, we observed that the soil-C content was highly related to the SADs of four out of 12 AM species [e.g. A. mellea, C. pellucida, R. tropicana and G. microaggregatum (Table 4)] in the topsoil layer under cropping conditions, suggesting a promoting role of the organic-C for growth and development of both plants and soil microorganisms 68,69 .
Previous investigations have demonstrated that the labile C fractions in soils provide energy source for soil microorganisms, and the proportion of organic matter in these fractions is a valuable indicator as for how the soils are biologically fertile 70 . Our findings are in agreement with the results of previous reports that showed that the community compositions of soil microbes, including bacteria and AM fungi, are under strong influences of the soil organic constituents 23,71,72 . Recent evidence has depicted the role of mycelia from AM fungi in decomposing the organic matter via the release of compounds like ammonia oxidizers, and secondary metabolites like phenolic compounds from the fungi to produce organic N for their uptake and utilization 23,72 . In addition, among the factorial predictors F, S and D, the D accounted for the highest variation of the AM community composition (61.9%) (Fig. 4b). These results are in agreement with recent observations that also indicated the influences of D on both the AMF species richness and community composition 32,65 .
Results from the SEM displayed various pathways in which each predictor exercised its effect (direct/indirect) on the AM community composition (Fig. 4c). The direct effect of the soil characteristics (SC) on structuring the AM community composition is consistent with previous observation 58 . The effects of fallow practices on the AM community compositions were indirect, while they were direct on the SC, with the observed strong path coefficients from the soil-C, -N and -Ca variables (Fig. 4c), suggesting that the fallow practices might affect the AM diversity by modifying the soil chemical properties, particularly the soil-C, -N and -Ca contents, to impose their influence on the surrounding microbial communities 32 . Interestingly, although soil-pH did not show a strong direct path from the model built by the SEM, it showed significant relationship with the AM community compositions in 0-10 cm soil depth under the forest and short fallows according to the results of RDA ( Fig. 3; Supplementary Table 3).
In conclusion, the present study showed the differential effects of the F, S and D, and their interactions on the soil properties. Furthermore, the F, S and D, and their interactions also significantly and differently affected the AM species richness and diversity, with the cropping practices revealing the highest effect, supporting the hypothesis of a positive feedback of the cropping practices on the AM structures under the conditions of the shifting agriculture in tropical Africa. Our findings also suggest that the fallow practices indirectly affect the AM diversity by altering the soil chemical properties via the production of high amount of C and N, which will subsequently influence the development of different AM species in soils.