Legume Nitrogen Fixation and Symbioses in Low-Inputs Rainfed Rice Rotations

: Cropping systems with legumes play key roles in farming systems in sub-Saharan Africa. However, how commonly legume associations perform in low input-systems is not well-known. Here, we studied four legume species used in three systems in rotation with upland rice, i.e., groundnut monocropping, sorghum–cowpea intercropping, and velvet bean–crotalaria intercropping, in two fertilization managements on the previous rice, i.e., manure alone or complemented with mineral fertilization. Legume suitability was assessed using rhizobial and mycorrhizal colonization rates, plant biomass production, shoot N and P content, and biological N 2 ﬁxation based on their δ 15 N natural abundance. Shoot and root biomasses varied signiﬁcantly between legume species and were positively correlated with nodule number (r = 0.49 and 0.74, p -value < 0.05 and <0.001, respectively) and the amount of ﬁxed N (r = 0.73 and 0.50, p -value < 0.001 and <0.05, respectively). The proportion of plant N derived from N 2 ﬁxation also varied signiﬁcantly between species, with a higher percentage for velvet bean (66%), compared to the other three species (50 to 60%). Legume roots were weakly colonized by AM fungi, with similar levels between species. Overall, fertilization management did not signiﬁcantly impact legume biomass, symbioses, or N 2 ﬁxation, yet the organo-mineral fertilization signiﬁcantly increased legume shoot P content. The lack of effect of mineral fertilization on N 2 ﬁxation and biomass could be due to other nutrient deﬁciencies (Ca, Mg, micronutrients), which can hamper symbioses with rhizobia and mycorrhizae.


Introduction
The conventional agriculture model, based on the use of synthetic inputs and nonrenewable natural resources, has made possible an increase in productivity but at the cost of damage to the environment [1]. Ecological intensification has been proposed as a nature-based alternative that sustains agricultural production while minimizing adverse effects on the environment [2]. Among the key practices to promote ecological intensification, crop diversification has shown positive impacts on both plant production and the environment [2][3][4]. It can contribute significantly to livelihoods, improved health and nutrition, household food security, climate resilience, and ecological sustainability [5].
The cropping season corresponds to the rainy season that lasts from October to April. Climatic data were recorded by an automatic meteorology station, ENERCO 404 Series (CIMEL Electronique, Paris, France). Rainfall was 1187 mm and mean temperature was 24.9 • C during the cropping season (see [26] for details on monthly rainfall and temperature during the experiment).
A field experiment was established in November 2015 (see also [26]) and hereby we focus on three two-year rotations including solely legumes or a crop mixture, namely (i) rice after groundnut (RG), (ii) rice after sorghum-cowpea (Vigna unguiculata) intercropping (RSC), and (iii) rice after velvet bean (Mucuna pruriens)-crotalaria (Crotalaria spectabilis) intercropping (RVC) (Figure 1). These were compared to rainfed rice monocropping (RR) in a factorial randomized block design with four blocks on plots of 45.9 m 2 each. The legumes and their cropping system integration were selected based on expert knowledge for different purposes: (i) RG for providing a cash crop and green manure with groundnut residues; (ii) RSC for providing forage for livestock with sorghum (grain and vegetative biomass) and grain opportunities for food and green manure with cowpea; and (iii) RVC for its potential to produce a high amount of green manure (combination of an erected plant and a climbing one) and to control plant-feeder nematodes [26]. Tillage, sowing, and harvest were performed manually. Groundnut, sorghum, cowpea, velvet bean, and crotalaria were sown at the end of November at a density of 17, 3, 8, 7, and 7 holes m −2 and of 1, 5, 5, 2, and 5 seeds per hole, respectively, in line with farmer practices. They were grown without any fertilizer or manure. However, the previous rice in the rotation was fertilized, with two modalities: (i) only organic fertilization (F1) with 5 t DM ha −1 year −1 manure was applied at sowing, directly into the rice holes, and (ii) organo-mineral fertilization (F2), with the same 5 t DM ha −1 year −1 manure addition than F1 and additional input of 400 kg ha −1 year −1 of NPK (11-22-6) and 200 kg ha −1 year −1 of urea. The manure is a low-quality farmyard manure, with an average composition for the four-year experiment of 0.68% dry, 18.6% organic matter, 0.65% N, 0.18% P, and 1.21% K [26]. NPK was applied during sowing, and urea was applied at 45 days and about 70 days after sowing, respectively, in line with recommendations by local extensionists. The F1 level represents the local practice, and F2 is a reference to achieving rice genetic potential. Some additional details on the management practices carried out during the four years of the trial are also available in a companion study [26].

Plant Sampling and Analyses
All legumes were sampled during the fourth cropping year, from 5 to 7 March 2019 (first sampling date), i.e., 14 weeks after sowing (WAS). Plants were collected between the full development and ripening stage of legumes: groundnut at the ripening stage (Stage 8 on the BBCH growth stage scale [27]); cowpea at the end of Stage 8; velvet bean at the full development stage (between Stages 4 and 5), and crotalaria at the flowering Stage 6.
Plants from four seed holes from the middle of each plot were taken. The plants were collected together with their roots and carefully washed. The plants were separated between shoot, root, and nodules, and the nodules were counted. The shoot, root, and nodule samples were oven dried at 60 • C for 48 h, weighed, and stored for analysis.
Shoot biomass was ground into a fine powder for chemical analysis. Total N content was determined by dry combustion using a Vario Pyro Cube CHN Analyzer (Elementar, Langenselbold, Germany). Total P content was determined after calcination of 0.5 g of ground material at 550 • C [28]. After cooling, the ash was dissolved in warm 2% HCl and analyzed by colorimetry using the blue method of Murphy and Riley [29]. The δ 15 N isotopic abundance of each legume plant, the rice grown in the trials, and of Digitaria sp.
(an abundant weed present in all plots of the trial) were analyzed on ground plant residues using a mass spectrometer (Isoprime Precision, Elementar, Langenselbold, Germany) coupled with a CN elemental analyzer (Vario Pyro Cube, Elementar). The air N 2 was used as reference to determine δ 15 N. The air N 2 scale for nitrogen isotopes is formally defined by atmospheric N 2 , which has an assigned δ 15 N value of 0‰.
These data were then used to determine legume %Ndfa, as follows [30]: where δ 15 Nleg correspond to the legume δ 15 N, and δ 15 Nref is the non-fixing reference plants δ 15 N, which corresponded to the average between rice and Digitaria sp. δ 15 N, as described by Samago et al. [31]. The "B" value is the δ 15 N of legume grown in a glasshouse on N-free substrate. The B values of groundnut (−2.27‰), crotalaria (−2.25‰), and velvet bean (−2.62‰) were taken from Okito et al. [32], and that of the cowpea (−2.20‰) was taken from Nyemba and Dakora [33]. The amount of N fixed by the crop was deduced by multiplying %Ndfa by shoot N content and shoot biomass. Plants from four other holes were collected as previously for mycorrhization rate determination. Roots were carefully washed with tap water, cut into 2 cm fragments, and stored at 4 • C in 70% alcohol. The root samples were stained with Schaeffer black ink as described in Vierheilig et al. [34]. Thirty root segments selected randomly were mounted on microscope slides for observation ( Figure 2). The frequency of mycorrhiza in the root system (F. Myc) and the root mycorrhization rate (M. Myc) were calculated according to Trouvelot et al. [35].
The transition from the individual plant scale (g m −2 ) to the intercropping plot scale (t ha −1 ) was performed in two steps. First, biomass of the association at the time of the first sampling was obtained by summing species individual biomass relatively to the surface areas they occupied in each plot. The amount of N fixed per ha was then deduced by the sum of each species N 2 fixation within the association (i.e., by multiplying the amount of N fixed by the shoot biomass of each species). Second, another subset of individual legume plants was harvested at maturity (April-May), 20 WAS for groundnut and 26 WAS for cowpea, crotalaria, and velvet bean (second sampling date). Fresh biomass was collected in two 1 m 2 quadrats and weighted. Sub-samples of approximately 200 g of fresh biomass were oven dried at 65 • C for 72 h to determine the dry weight (DW) produced at the harvest stage.

Statistical Analysis
The influence of legume species and fertilization levels on individual plant traits was first assessed with two-way ANOVAs with interactions. The nodulation variables (abundance and dry weight) that did not meet the criteria of homoscedasticity were analyzed with a linear model coupled with ANOVA with White correction [36] with the ANOVA function of the {car} R package. Since there were no legume species × fertilization interactive effects for all variables (excepted shoot N, which was analyzed separately), differences between legume species and fertilization managements were tested separately a second time with a post hoc Tukey test. Relationships between legume properties were then assessed with a principal component analysis (PCA) and a Pearson correlation matrix with Benjamini Hochberg correction for multiple comparisons. All statistical analyses were performed using R software (R-4.0.4) and the following packages: agricolae, car, factoextra, FactoMineR, multcompView, psych, and stats.

Legume Growth and Biological N2 Fixation
Shoot and root biomasses varied significantly between legume species, whereas fertilization had no effect on these variables ( Table 1). The highest values of shoot biomass were found for groundnut and crotalaria (respectively, 185 ± 65 and 181 ± 52 g m −2 ), while the lowest ones were observed for velvet bean and cowpea (respectively, 87 ± 36 and 109 ± 19 g m −2 ). On the other hand, all species significantly differed in their root biomass, and were, from the highest to the lowest values: groundnut, crotalaria, cowpea and velvet bean. These differences in shoot and root biomasses resulted in variable shoot: root ratio, with the highest and lowest values for velvet bean (14.4 ± 4.0) and groundnut (3.3 ± 0.6), respectively, while both crotalaria and cowpea presented intermediate values ( Table 1). The root biomass represented 18-29% of the shoot biomass for the various species, except for velvet bean, where it represented only 7% of the aboveground biomass.

Statistical Analysis
The influence of legume species and fertilization levels on individual plant traits was first assessed with two-way ANOVAs with interactions. The nodulation variables (abundance and dry weight) that did not meet the criteria of homoscedasticity were analyzed with a linear model coupled with ANOVA with White correction [36] with the ANOVA function of the {car} R package. Since there were no legume species × fertilization interactive effects for all variables (excepted shoot N, which was analyzed separately), differences between legume species and fertilization managements were tested separately a second time with a post hoc Tukey test. Relationships between legume properties were then assessed with a principal component analysis (PCA) and a Pearson correlation matrix with Benjamini Hochberg correction for multiple comparisons. All statistical analyses were performed using R software (R-4.0.4) and the following packages: agricolae, car, factoextra, FactoMineR, multcompView, psych, and stats.

Legume Growth and Biological N 2 Fixation
Shoot and root biomasses varied significantly between legume species, whereas fertilization had no effect on these variables ( Table 1). The highest values of shoot biomass were found for groundnut and crotalaria (respectively, 185 ± 65 and 181 ± 52 g m −2 ), while the lowest ones were observed for velvet bean and cowpea (respectively, 87 ± 36 and 109 ± 19 g m −2 ). On the other hand, all species significantly differed in their root biomass, and were, from the highest to the lowest values: groundnut, crotalaria, cowpea and velvet bean. These differences in shoot and root biomasses resulted in variable shoot: root ratio, with the highest and lowest values for velvet bean (14.4 ± 4.0) and groundnut (3.3 ± 0.6), respectively, while both crotalaria and cowpea presented intermediate values ( Table 1). The root biomass represented 18-29% of the shoot biomass for the various species, except for velvet bean, where it represented only 7% of the aboveground biomass. Shoot N and P contents exhibited contrasting responses to legume species and fertilization levels. Shoot N content differed significantly between species only in F2, where crotalaria presented the highest N content (3.36 ± 0.22% dry matter), groundnut the lowest shoot N content (2.59 ± 0.10% dry matter), and velvet bean and cowpea showed intermediate values (Table 1). In contrast, shoot P content did not differ between species but was significantly impacted by fertilization levels, with higher content in F2 than F1 (1.83 ± 0.30‰ against 1.55 ± 0.32‰ dry matter, respectively), leading to higher shoot N: shoot P ratio in F1 than F2 modality.
The highest shoot δ 15 N was observed for crotalaria (3.14 ± 1.29‰) and the lowest for the three other species, with values between 0.8 and 2‰. The percentage of plant N derived from N 2 fixation varied significantly between species, with a higher %Ndfa for velvet bean (66 ± 8%) compared to the other three species (50% to 60%). However, we observed no effect on the percentage of plant N using the different fertilization methods. This resulted in an average amount of fixed N in shoot biomass of 2.25 ± 0.82 g m −2 at the first sampling date, all modalities taken into account, without significant differences between species or fertilization management.

Legume Symbioses
The analysis of variance showed significant differences in nodulation between the species of legumes (p-values < 0.001), while fertilization had no effect on this symbiosis (p-values > 0.05, Table 1). Groundnut presented the highest nodule number out of all other species (10 to 19 times higher than the other species, all fertilization modes taken into account). However, both groundnut and velvet bean presented high level of total nodule dry weight, corresponding to many small nodules for groundnut (average weight of 5.5 mg per nodule) and few big nodules for velvet bean (average weight of 35.8 mg per nodule). On the other hand, both crotalaria and cowpea presented low nodule number and dry weight.
Legumes roots were weakly colonized by AM fungi, without significant effect of legume species and fertilization. The frequency of mycorrhiza in the root system (F. Myc) varied between 5% to 11%, and the root mycorrhization rate (M. Myc) was only about 1% (Table 1).

Legume Response to Symbioses
The PCA biplot indicated strong positive covariations between shoot and root biomass production, and to a lesser extent with the nodulation and the amount of fixed N on the first axis of the PCA (Figure 3). In contrast, shoot: root ratio and Shoot N content negatively covariated on the same axis. Contrastingly, shoot P content and to a lesser extent mycorhization frequency covariated positively on the second axis, while shoot N:P ratio covariated negatively on this axis. Traits contrasted between species were mainly displayed on the first axis, with globally higher biomass production and nodulation levels for groundnut and crotalaria than cowpea and velvet bean (Figure 3). The relationship between the biomass produced (shoot and root) and the amount of fixed N was highlighted by Pearson correlations (Table 2), with significant correlations between root or shoot biomass and nodule number (r = 0.74 and 0.49, respectively) and significant correlation between shoot biomass with the amount of fixed N (r = 0.73).

Biomass Production and N Restitution at the Plot Scale
At the intercropping plot scale, shoot biomass of legumes nearly at the flowering stage was the highest for velvet bean and crotalaria intercropping (2.68 ± 0.69 t ha −1 ) and the lowest for cowpea intercropped with sorghum (1.26 ± 0.55 t ha −1 ), while groundnut

Biomass Production and N Restitution at the Plot Scale
At the intercropping plot scale, shoot biomass of legumes nearly at the flowering stage was the highest for velvet bean and crotalaria intercropping (2.68 ± 0.69 t ha −1 ) and the lowest for cowpea intercropped with sorghum (1.26 ± 0.55 t ha −1 ), while groundnut monocropping presented intermediate values (1.85 ± 0.65 t ha −1 ) (Figure 4a). The corresponding amounts of fixed N in shoot biomass were 25.3 ± 9.8 kg N ha −1 for groundnut monocropping, 21.0 ± 9.0 kg N ha −1 for sorghum-cowpea intercropping, and 47.2 ± 11.5 kg ha −1 Sustainability 2021, 13, 12349 9 of 14 for velvet bean-crotalaria intercropping (Figure 4b). At harvest, shoot biomass decreased for species that were at a more advanced development stage at the first sampling date (groundnut, cowpea: Stage 8) by 51% and 22%, respectively (Figure 4c). However, velvet bean and crotalaria that were between Stage 4 to 6 at the sampling date continued to grow by 123% (from 2.68 to 5.98 t ha −1 ) until harvest.  (Figure 4a). The corresponding amounts of fixed N in shoot biomass were 25.3 ± 9.8 kg N ha −1 for groundnut monocropping, 21.0 ± 9.0 kg N ha −1 for sorghum-cowpea intercropping, and 47.2 ± 11.5 kg ha −1 for velvet bean-crotalaria intercropping (Figure 4b). At harvest, shoot biomass decreased for species that were at a more advanced development stage at the first sampling date (groundnut, cowpea: Stage 8) by 51% and 22%, respectively (Figure 4c). However, velvet bean and crotalaria that were between Stage 4 to 6 at the sampling date continued to grow by 123% (from 2.68 to 5.98 t ha −1 ) until harvest. Significance of the differences between intercropping and fertilization were tested with two-way ANOVAs with interactions. Since only intercropping presented a significant effect, ANOVAs were coupled with post-hoc Tukey tests. Treatments with significant differences (p-values < 0.05) are shown with different lower case letters.

Importances of Rhizobial and AMF Symbionts on Plant N2 Fixation in Acid Ferralsols
Nodulation of all studied legume species was observed with a quite large number of nodules per plant (Table 1), similar or higher than that reported in other studies in SSA (e.g., [37,38]). Efficient native strains seem to be present in the soil and able to nodulate the legume species. For some other species (e.g., common bean), efficient native strains could be lacking in Malagasy soils and inoculation with rhizobium could have a positive effect [24]. Our results also showed an overall significant relationship between shoot and root DW, nodule number and DW, and fixed-N ( Figure 3 and Table 2). This results in a mean of 16.5 kg fixed N t −1 DM (fixed N = 16.49 × Shoot biomass − 0.64; R 2 = 0.85), in line with the data of Unkovich et al. [39]. Regression analysis between shoot dry matter and amounts of shoot fixed N are generally strong for most legumes [39,40]. The meta-analysis performed by Unkovich et al. [39] for Australian studies shows that the relationship between N 2 fixation and shoot dry matter had slopes between 20.2-24.3 kg N 2 fixed t −1 DW for pasture legumes and between 10.7 to 23.0 kg N 2 fixed t −1 DW for crop legumes, depending upon species.
Although rice in rotation with legume was fertilized, a legume P limitation could be suspected, as suggested by their high shoot N: shoot P ratio of shoot biomass, that was 17-19 (Table 1). Indeed, a N:P ratio above 16 generally suggests P-limited growth, while an N: P ratio below 13.5 suggests N-limited growth [41]. It is well known that N 2 fixation is restricted by low P availability [20]. The lack of adequate levels of available P in tropical soils, as is found in Madagascar, is one of the major constraints for crop production [42] due to the high P sorption on Fe-and Al-oxides and clay minerals [43]. The organo-mineral fertilization (F2) provided to the previous rice only resulted in a significant increase in shoot P content, without any effect on shoot biomass or nodulation (Table 1). According to Kouas et al. [44], nodule numbers are significantly diminished by low P availability, whereas the intrinsic characteristics of the nodules (individual biomass and size, P concentration, and efficiency of N 2 fixation) do not depend on P availability. The increase in P content with F2 fertilization showed that the bioavailability of P was improved. Therefore, the lack of response of nodulation to increased P bioavailability could be linked to other macro-or micro-nutrient deficiencies, likely to occur in these multi-deficient soils [42]. In addition, the nutrient deficiencies that affect the soil (N, P, Ca, and Si) [42] lead to a large decrease in plant growth. Indeed, as shown by Santachiara et al. [45] on soybeans, legumes can be affected by various nutrient deficiencies, such as potassium, sulfur, molybdenum, iron, etc., in addition to P, most of them having a critical control on BNF. In such degraded soils, the poor soil fertility often results in a large gap between the potential and observed yields (e.g., [7]).
It is also well-known that N 2 fixation is restricted by high N availability [46]. The %Ndfa of legumes declines with high concentrations of soil nitrate, but concentrations <20-25 kg N ha −1 are generally considered too low to have an impact on %Ndfa [46]. Results of a companion study on the same trial [26] show that total mineral N content was only 13 kg N ha −1 at the harvest of the previous rice. However, nitrate levels were only measured in rice and not in legume plots. The soil mineral N had probably no or a low inhibitory effect on nodulation under these conditions. The low nitrogen content of the soil after rice cultivation may also explain the lack of effect of mineral fertilization on the nitrogen content of legume stems.
Unlike nodulation, the mycorrhization rates were very low (Table 1). Our findings for F. Myc (5-9%) were generally lower than the observations from studies published in SSA on the plants we were studying (mycorrhization rate of 2-35% for velvet bean [47,48]; 20-60% for cowpea [49]). This low mycorrhization rate is not well-understood, and different hypotheses could be proposed, such as a very low amount of soil spores or the inhibitory effects of some practices (e.g., fertilizer inputs) on mycorrhization. According to Lambers et al. [50], mycorrhizal strategies are relatively inefficient in highly weathered soils, presumably because AMF are not able to mobilize P from poorly available pools for plant roots. Moreover, the mycorrhization rate is usually restricted when P availability is low, whereas AM colonization increases after P additions in P-limited soils [48,51]. Indeed, we showed a positive relationship on the PCA between shoot P and F. Myc (Figure 3). A recent study performed on a Ferralsol typical of the hills of the Malagasy Highlands has also shown that the inoculation of a legume (common bean) with efficient native strains of AM fungi collected on the same soil promotes a significant increase in shoot biomass, plant P content, and soil phosphatase activity [24].

Amount of N Fixed by Legumes in SSA
Estimates of N inputs by BNF in cropping systems in SSA are of the uppermost importance. Nevertheless, there is a limited number of studies on the extent of N 2 fixation (%Ndfa, amount of N fixed) by legumes and their contribution to N input in cropping systems in SSA, with most studies focusing on a few grain legumes (cowpea, groundnut, Bambara groundnut, and soybean) (e.g., [9]). Moreover, to our knowledge, except the pioneering study of Becker and Johnson [52], none of them focus on forage or fallow legumes (e.g., Cajanus, Crotalaria, Dolichos, Mucuna, Stylosanthes, etc.).
Our results of N 2 fixation (Table 1) were based on selected B values from the literature and measurements of above-ground biomass N at various stages of plant development during the growing season. The B values for a given plant species can vary due to several parameters, such as the growth stage of the plant [53], the legume cultivar [54], the rhizobium strain(s) that colonizes the roots [32], and differences between field and greenhouse results [55]. Concerning measurements on biomass, we did not take into account the contribution of fixed N derived from the roots nor the contribution of total biomass reached at maximum plant development. An estimate of the quantity of N fixed considering these two factors can be made. Since the %Ndfa of nodulated roots and stems were found to be almost identical [46], a 'root factor' was applied to shoot-based data to estimate total plant N (i.e., shoot N + nodulated root N) [39]. Considering the root: shoot ratio at the sampling stage, the amount of N fixed would increase by 7-29% according to legume species to take into account the N of nodulated roots. An estimate of biomass production can also be difficult given the senescence of plants and leaf decay after flowering. For example, the biomass was lower at harvest for groundnut and cowpea than that at flowering, whereas it increased by 123% for velvet bean-crotalaria intercropping (Figure 3). For our calculations, we used the maximum biomass obtained of the two sampling dates. The estimated amounts of fixed N in total biomass (shoot and root) at maximum development stage of the legume species were 33 kg ha −1 for groundnut monocropping, 24 kg ha −1 for sorghum-cowpea intercropping, and 121 kg ha −1 for velvet bean-crotalaria intercropping. The corresponding amounts of soil-derived N in legumes were 28 kg ha −1 for groundnut monocropping, 17 kg ha −1 for sorghum-cowpea intercropping, and 89 kg ha −1 for velvet bean-crotalaria intercropping.
In the trial we studied, the hierarchy observed for the amount of fixed N by the different cropping systems, i.e., a higher amount for velvet bean-crotalaria intercropping compared to sorghum-cowpea and groundnut intercropping, is mirrored by the differences in rice yield and rice N content observed in the following rice crop cycle (see [26]). This highlights the greater contribution of forage, fodder, and green manure legumes with regard to grain legumes to provide a large amount of fixed N in agricultural systems [18].

Conclusions
Symbiotic N 2 fixation by legumes can add a significant amount of nitrogen for the next crop in the rotation through the incorporation of residues into the soil. The overall inputs of N 2 were estimated at 25-30 kg ha −1 for groundnut monocropping and sorghumcowpea intercropping and 120 kg ha −1 for velvet bean-crotalaria intercropping. The choice of legume species to integrate into the cropping system will depend on the N balance but also on the profitability of the cropping system over the entire rotation cycle, as well as constraints linked to the production system (e.g., mulching of crop residues in conservation agriculture system or harvesting of forage in crop-livestock integrated systems). In addition, the lack of effect of NPK fertilization on highly weathered soils with low nutrient availability could be due to multiple nutrient deficiencies (K, Ca, Mg, or micronutrients). which can affect symbiosis with rhizobia or mycorrhizae. Then, the increase in the productivity of legumes requires managing fertilization in order to improve the efficiency of rhizobial or mycorrhizal symbioses. In the highland of Madagascar, inputs of P and Ca seem essential for a greater efficiency of the symbionts. The presence of native rhizobial strains in the soil may allow the nodulation of the various species studied, whereas the mycorrhization rate remained weak. Further studies should investigate if the mycorrhization rate could be increased following AMF inoculation in the soil with more efficient native or commercial mycorrhizal strains or with additional fertilization.