# Impact of Wheat/Faba Bean Mixed Cropping or Rotation Systems on Soil Microbial Functionalities

Sanâa Wahbi, Yves Prin, Jean Thioulouse, Hervé Sanguin, Ezékiel Baudoin, Tasnime Maghraoui, Khalid Oufdou, Christine Le Roux, Antoine Galiana, Mohamed Hafidi, et al.

► To cite this version:

Sanâa Wahbi, Yves Prin, Jean Thioulouse, Hervé Sanguin, Ezékiel Baudoin, et al. Impact of Wheat/Faba Bean Mixed Cropping or Rotation Systems on Soil Microbial Functionalities. Frontiers in Plant Science, 2016, 7, pp.1364. 10.3389/fpls.2016.01364. hal-01911655

HAL Id: hal-01911655

https://hal.science/hal-01911655

Submitted on 18 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Impact of Wheat/Faba Bean Mixed Cropping or Rotation Systems on Soil Microbial Functionalities

Sanâa Wahbi¹,², Yves Prin³, Jean Thioulouse⁴, Hervé Sanguin³, Ezékiel Baudoin¹, Tasnime Maghraoui¹,⁴,⁵, Khalid Oufdou⁵, Christine Le Roux³, Antoine Galiana³, Mohamed Hafidi² and Robin Duponnois¹*

¹IRD, UMR LSTM, Montpellier, France, ²Laboratoire Ecologie et Environnement (Unité Associée au CNRST, URAC 32), Faculté des Sciences Semlalia, Université Cadi Ayyad, Marrakech, Morocco, ³CIRAD, UMR LSTM, Montpellier, France, ⁴Université Lyon 1, CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, Villeurbanne, France, ⁵Laboratoire de Biologie et de Biotechnologie des Microorganisms, Faculté des Sciences Semlalia, Université Cadi Ayyad, Marrakech, Morocco

Cropping systems based on carefully designed species mixtures reveal many potential advantages in terms of enhancing crop productivity, reducing pest and diseases, and enhancing ecological services. Associating cereals and legume production either through intercropping or rotations might be a relevant strategy of producing both type of culture, while benefiting from combined nitrogen fixed by the legume through its symbiotic association with nitrogen-fixing bacteria, and from a better use of P and water through mycorrhizal associations. These practices also participate to the diversification of agricultural productions, enabling to secure the regularity of income returns across the seasonal and climatic uncertainties. In this context, we designed a field experiment aiming to estimate the 2 years impact of these practices on wheat yield and on soil microbial activities as estimated through Substrate Induced Respiration method and mycorrhizal soil infectivity (MSI) measurement. It is expected that understanding soil microbial functionalities in response to these agricultural practices might allows to target the best type of combination, in regard to crop productivity. We found that the tested cropping systems largely impacted soil microbial functionalities and MSI. Intercropping gave better results in terms of crop productivity than the rotation practice after two cropping seasons. Benefits resulting from intercrop should be highly linked with changes recorded on soil microbial functionalities.

**Keywords:** arbuscular mycorrhizal fungi, cropping systems, nutrient uptake efficiency, microbial soil functions, Mediterranean region

## INTRODUCTION

The ecological key processes that warrant the productivity and stability of terrestrial ecosystems have been designated as efficient models for sustainable agricultural management. Since it is well known that the functioning and stability of terrestrial ecosystems are dependent to plant biodiversity and species composition (Naeem et al., 1994; Tilman and Downing, 1994; Tilman et al., 1996; Hooper and Vitousek, 1997), mixing plant species in cropping systems reveals many potential advantages under various conditions to contribute to modern and sustainable agriculture

## OPEN ACCESS

**Edited by:**
Khaled Masmoudi,
Center of Biotechnology of Sfax,
Tunisia

**Reviewed by:**
Oswaldo Valdes-Lopez,
National Autonomous University
of Mexico, Mexico
Zonghua Wang,
Fujian Agriculture and Forestry
University, China

*Correspondence:
Robin Duponnois
robin.duponnois@ird.fr

**Specialty section:**
This article was submitted to Plant Biotic Interactions, a section of the journal Frontiers in Plant Science

**Received:** 02 May 2016
**Accepted:** 29 August 2016
**Published:** 15 September 2016

**Citation:**
Wahbi S, Prin Y, Thioulouse J, Sanguin H, Baudoin E, Maghraoui T, Oufdou K, Le Roux C, Galiana A, Hafidi M and Duponnois R (2016) Impact of Wheat/Faba Bean Mixed Cropping or Rotation Systems on Soil Microbial Functionalities. Front. Plant Sci. 7:1364. doi: 10.3389/fpls.2016.01364
(Vandermeer, 1989). In agroecosystems, multispecies cropping systems may (i) sustain biomass production and decrease the risk of crop failure in unpredictable environments, (ii) rehabilitate disturbed ecosystem services (i.e., water and nutrient cycling), and (iii) decrease risks of invasion, pests, and diseases through enhanced biological control or direct control of pests (Gurr et al., 2003). These cultural practices have often been considered as a practical application of ecological principles related on biodiversity, plant interactions, and other natural regulation mechanisms (Malézieux et al., 2009). Numerous different multispecies cropping systems can be designed by including various criteria such the frequency of land-use rotation, the intensity of intercropping, etc. (Garcia-Barrios, 2003). For instance, the intercropping, defined as the growth of more than one crop species or cultivar simultaneously in the same field during the same growing season (Ofori and Stern, 1987; Hauggaard-Nielsen et al., 2007), enhanced the use efficiency of environmental sources for plant growth resulting in stable yields (Hauggaard-Nielsen et al., 2001; Corre-Hellou and Crozat, 2005; Jensen et al., 2015). It has been also reported that the association of cereals and legumes at the same space and time led to higher yields and improved N (via biological N₂ fixation for the legume) and P nutrition (Li et al., 2005; Betencourt et al., 2012; Latati et al., 2013, 2014).

The positive effects of species diversity in intercropping systems result from two main processes: complementarity and facilitation (Fridley, 2001; Hinsinger et al., 2011) whereas in rotation systems (i.e., legumes/cereals rotation), they occur through indirect feedback interactions (Schnitzer et al., 2011). It has been reported that these biological processes were mainly driven by soil microbe activities (Klironomos, 2002; Tang et al., 2014). For instance, the enhancement of P acquisition in the context of cereal/legume intercrops could occur as a consequence of microbial mediated processes involving soil fungi and bacteria (Wang et al., 2007). Positive feedback often involves changes in abundance and diversity of symbiotic mutualists such as nitrogen-fixing rhizobacteria and mycorrhizal fungi (van der Putten et al., 2013). These biological processes have been particularly studied for their impacts on primary production, nutrient retention, and resilience after stress in natural prairie ecosystems (Tilman et al., 1996, 1997) or in natural forest ecosystems (Altieri, 1999). In contrast, few studies have been conducted in agricultural systems to determine the impacts of multispecies cropping systems on microbial soil functionalities and mycorrhizal soil infectivity in relation with crop productivity.

It has been suggested that Faba bean could enable diversification of the agrosystems (Garofalo et al., 2009; Köpke and Nemecek, 2010) and it has been highlighted the importance of the N and P contributed by faba bean in intercropping and rotation systems (Li et al., 2003). However, the Faba bean impacts on soil microbial functionalities in field conditions have been less investigated (Köpke and Nemecek, 2010). In the present study, two cereal/legume systems (intercropping and rotation) were investigated in a field experiment during two growing seasons. The specific aims of this study were (i) to evaluate the agronomic performance of wheat/Faba bean in intercrop and in rotation and (ii) to monitor the impacts of these cultural

practices on the mycorrhizal soil infectivity and on the microbial soil functionalities.

# MATERIALS AND METHODS

## Field Conditions and Experimental Design

The field experiment was performed in the 2011–2013 growing seasons in the Haouz plain at about 30 km at the East of Marrakech (31°4′60″ N and 7°3′0″ W, Morocco). Soil chemical properties (0–0.10 m layer) were as follows: pH (H₂O) 7.2; carbon (%) 1.54; nitrogen (%) 0.08, C/N 19.2; Total P (mg·kg⁻¹) 502.9 and Olsen P (mg·kg⁻¹) 22.1. This soil exhibited a high content of available P presumably non-limiting for plant growth and resulting from large-P fertilizer applications during the last decades, until this part of the field was dedicated to organic farming 1 year before the beginning of this experiment. The regional climate of the experimental site is typical Mediterranean with surface soils regularly undergoing drying-rewetting cycles from the irregular distribution of rainfall. The annual average rainfall was 282 mm, mostly in Autumn/Winter (59%) and in Spring (22%). The dry season is from May to September. The mean air temperature is 17.9°C in autumn, 12.8°C in winter, 18.5°C in spring, and 24.7°C in summer. Soil, cropped in the previous growing season with durum wheat (Triticum durum Desf.), was plowed to a depth of 0.30 m in Summer and then shallowly harrowed to control weeds. No herbicide nor chemical fertilizers were applied. The experimental design had a randomized block design with two factors and four replication blocks. The factor was the cropping system (wheat/faba bean intercropping, rotation, or wheat monoculture). Plots were 3.0 m × 3.0 m; each main plot was spaced 1.0 m out from the next. The crops were sown in December 2011 and 2012 at a rate of 400 viable seeds m⁻² in rows 0.18 m apart for wheat and at a rate of 200 kg·ha⁻¹ for faba bean. When intercropped, the two species were sown in the same row in order to maximize root proximity and plant-plant interactions. The experimental plot consisted of four rows 3 m long. Weeds were controlled by hand during the experiment. Hence, 3 treatments were examined, namely durum wheat as sole crop (W) and two cropping systems: durum wheat/Faba bean intercrop (WF) and wheat/Faba bean rotation (F+W).

## Plant Analyses

Wheat was considered as the main crop and faba bean as an intercrop component. Hence, the expected benefits of the cultural practices on yield productivity were only examined on wheat plants. After one and two growing seasons, at wheat tillering and at the same time, the total number of wheat plant and spike per plot were counted. Then 10 randomly chosen plants in the middle of the plot were harvested. The seeds from each plant were collected, counted, and weighed to determine the dry weight of 1000 seeds (grain yield). These measurements recorded on a plot basis were converted to hectare for statistical analysis. The aerial parts of each plant were then oven dried at 70°C during 2 weeks and weighed. After drying, shoot tissues were ground,
ashed (500°C), digested in 2 ml HCl 6N and 10 ml HNO₃ N for nitrogen and then analyzed by colorimetry for phosphorus (John, 1970). For nitrogen determination (Kjeldahl method), they were digested in 15 ml H₂SO₄ (36 N) containing 50 g l⁻¹ of salicylic acid. Roots from five other randomly chosen plants in the middle of each plot were sampled and root subsamples of about 3 g each were taken. Each subsample was stained with 0.05% trypan blue in lactic acid according to Phillips and Hayman (1970); root colonization by AMF was then measured with the grid intersect method according to Giovannetti and Mosse (1980).

## Soil Microbial Analysis

After one and two growing seasons, soil cores (1 kg) were collected at 0- to 20-cm depth in each plot. About 10 soil samples were taken from each plot and pooled together. Soil samples were crushed and passed through a 2-mm sieve. Then hyphal length that is considered as a main component of the mycorrhizal soil infectivity (Kisa et al., 2007) was measured by the filtration-grid-line method on membrane filters according to Jakobsen and Rosendahl (1990). Patterns of *in situ* catabolic potential (ISCP) were designed to assess the functions of soil microbial communities and the microbial functional diversity in soil treatments after one and two growing seasons. The diversity of the catabolic potentials of the total soil bacterial community was evaluated according to Campbell et al. (2003) by a microrespirometry method performed in 96-well microtiter plates. In order to ensure the resumption of microbial activity, sterile distilled water was added to reach 30% of the water-holding capacity and plates were incubated 3 days in the dark at 28°C. Then soil wells received 28 organic substrate solutions (three wells per substrate). Stock solutions for thirteen carbohydrates (D-mannose, D-mannitol, D-trehalose, L-arabinose, D-xylose, D-sucrose, D-galactose, meso-inositol, D-sorbitol, L-rhamnose, L-arabitol, mesoerythriol, D-Glucose), eight carboxylic acids (citric acid, maleic acid, D,L-malic acid, oxalic acid, Na-gluconate, α-ketoglutaric acid, L-ascorbic acid) and seven amino acids (L-asparagine, D,L-valine, L-methionine, L-glutamine, D,L-alanine, D,L-serine, N-acetyl-D-Glucosamine) were prepared with distilled water and their concentrations were calculated to lead, respectively, 0.03, 0.04, and 0.004 mmol g⁻¹ soil. Basal respiratory activity was calculated in triplicate with distilled water. The colorimetric detection plates were assembled and used according to MicroResp™ (Aberdeen, UK) recommendations. Absorbance was measured at 572 nm with a Tecan infinite M200 reader before substrate spiking (t0) and after 6 h of incubation at 28°C (t6). For each well, absolute respiratory activity was calculated by subtracting the absorbance value at t0 from the value at t6. The average basal respiration value was then subtracted from all the individual substrate respiration values. For each carbon source, this substrate-specific respiratory activity was averaged and the value was finally divided by the sum of all the mean substrate-specific respiratory activities (pi value). The catabolic evenness (E) was calculated to determine the catabolic diversity of soil treatments. It represents the variability of catabolized substrates amongst the range of the targeted substrates and is calculated using the Simpson-Yule index $E = 1/\sum p_i^2$ with $p_i = (\text{respiration}$

response to individual substrates)/(Total respiration activity induced by all substrates for a soil treatment; Magurran, 1988). The catabolic evenness could be used to evaluate the ability of microbial communities to resist against environmental stress or disturbance (Magurran, 1988). Data were calculated for the individual responses to substrates but also for the average responses with carbohydrates, carboxylic acids, and amino acids.

## Statistical Analysis

All the data were subjected to a two-way analysis of variance and comparisons among means were made using the Newman-Keuls test (*P* < 0.05). The percentages of the mycorrhizal colonization were transformed by arcsin(sqrt) before the statistical analysis. The relationships between ISCP profiles table (2012 data) and yield variables table (2013 data) were analyzed using Between-Group CO-Inertia Analysis (BGCOIA). BGCOIA is a Co-Inertia Analysis on the two tables of group means obtained after a Between-Group Analysis (BGA, Thioulouse et al., 2012). As a first step, a BGA is therefore computed on the two data sets, considering each treatment as a group. The technical details of BGCOIA are given in Franquet et al. (1995). Examples of use and a comparison with other methods are presented in Thioulouse (2011) in the framework of k-tables data analysis methods. Let g be the number of groups (treatments here). The table of group means for SIR profiles is obtained by computing the means of each substrate within each treatment. This gives a new table, with g rows and p columns (p substrates). The same computations are done for the yield data table, leading to a second new table with g rows and q columns (q yield variables). A Co-Inertia Analysis is then performed on these two new tables. The rows of the initial tables can be projected into this analysis to help interpret the results (Lebart et al., 1984). Computations and graphical displays can be produced with the ade4 package for the R software (Thioulouse and Dray, 2007).

## RESULTS

### Grain Yield and Biomass Production

After one cropping season, no significant effect on the wheat development has been recorded between the wheat monoculture and the wheat/faba bean intercropping (Table 1). After the second cropping season, crop data evidenced that total biomass yield, spike number, spike dry weight were higher in the intercropping treatment (WF) compared to the wheat monoculture (W) and Faba bean/wheat rotation (W+F) treatments with some enhancements resulting from the intercropping vs. monoculture of +84.2, +24.8, and 122.7%, respectively (Table 1). The thousand-seed weight was significantly higher in the WF and W+F treatments compared to the monoculture (W; Table 1). After the second growing season, the shoot P content was significantly increased in the intercropping and rotation treatments whereas nitrogen enhancement was only recorded in the WF treatment with (Table 1).
<table><caption>TABLE 1 | Effects of cropping systems (W: wheat monoculture; WF: wheat/Faba bean intercropping; W+F: Wheat/Faba bean rotation) on total biomass yield (kg.ha⁻<fnref n="1" />), spike number per ha, spike dry weight (kg.ha⁻<fnref n="1" />), thousand-seed weight (TSW), mineral nutrition, mycorrhizal colonization of durum wheat and on soil catabolic evenness, and standardized average substrate-induced respiration (SIR) responses with each substrate group (carboxylic acids, amino-acids, and carbohydrates) in 2012 and 2013 in the field experiment located in the Haouz plain at about 30 km at the East of Marrakech (Morocco).</caption><thead><tr><th rowspan="3"></th><th colspan="6">Treatments</th></tr><tr><th colspan="3">2012</th><th colspan="3">2013</th></tr><tr><th>W</th><th>WF</th><th>W+F</th><th>W</th><th>WF</th><th>W+F</th></tr></thead><tbody><tr><td>Total biomass yield (kg.ha<sup>-1</sup>)</td><td>4662 (78.5)<sup>1</sup> a<sup>2</sup></td><td>3897 (49.3) a</td><td>-</td><td>4434 (66.1) a</td><td>8167 (97.3) b</td><td>5038 (94.6) a</td></tr><tr><td>Spike number per ha (× 10<sup>4</sup>)</td><td>219.3 (12.2) a</td><td>213.7 (24.4) a</td><td>-</td><td>192.3 (37.5) a</td><td>240.1 (10.2) b</td><td>179.1 (39.9) a</td></tr><tr><td>Spike dry weight (kg.ha<sup>-1</sup>)</td><td>2290 (12.6) a</td><td>2283 (33.4) a</td><td>-</td><td>2224 (17.1) a</td><td>4953 (58.3) b</td><td>2661 (41.5) a</td></tr><tr><td>Thousand-seed weight (g)</td><td>42.7 (2.9) a</td><td>41.1 (1.9) a</td><td>-</td><td>42.7 (2.9) a</td><td>53.9 (2.7) b</td><td>52.8 (4.5) b</td></tr><tr><td>Shoot N content (%)</td><td>nd<sup>3</sup></td><td>nd<sup>3</sup></td><td>-</td><td>5.08 (0.31) a</td><td>5.55 (0.17) b</td><td>4.71 (0.41) a</td></tr><tr><td>Shoot P content (g.kg<sup>-1</sup>)</td><td>nd<sup>3</sup></td><td>nd<sup>3</sup></td><td>-</td><td>6.01 (0.59) a</td><td>7.47 (0.38) b</td><td>7.56 (0.34) b</td></tr><tr><td>Mycorrhizal colonization (%)</td><td>46.7 (5.8) a</td><td>62.4 (6.4) b</td><td>-</td><td>54.7 (4.7) a</td><td>74.7 (4.8) b</td><td>69.9 (3.2) b</td></tr><tr><td>Hyphal length (m g<sup>-1</sup> dry soil)</td><td>1.66 (0.08) a</td><td>2.85 (0.07) c</td><td>2.14 (0.06) b</td><td>1.72 (0.07) a</td><td>2.98 (0.05) b</td><td>2.85 (0.04) b</td></tr><tr><td>Catabolic evenness</td><td>15.2 (0.83) b</td><td>12.3 (0.23) b</td><td>17.4 (0.48) c</td><td>15.1 (1.13) b</td><td>12.7 (1.01) a</td><td>13.5 (1.23) ab</td></tr><tr><td>Carbohydrates</td><td>0.025 (0.01) a</td><td>0.029 (0.02) b</td><td>0.028 (0.02) b</td><td>0.024 (0.02) a</td><td>0.021 (0.01) a</td><td>0.023 (0.02) a</td></tr><tr><td>Amino-acids</td><td>0.008 (0.01) a</td><td>0.006 (0.01) a</td><td>0.008 (0.01) a</td><td>0.007 (0.01) a</td><td>0.005 (0.01) a</td><td>0.006 (0.01) a</td></tr><tr><td>Carboxylic acids</td><td>0.067 (0.04) a</td><td>0.102 (0.07) b</td><td>0.060 (0.03) a</td><td>0.064 (0.04) a</td><td>0.075 (0.05) b</td><td>0.078 (0.05) b</td></tr></tbody></table>

## Mycorrhizal Soil Infectivity and Soil Microbial Functionalities

The mycorrhizal colonization of wheat root systems was significantly higher in the WF treatment compared to the wheat monoculture after one cropping season (Table 1). After two cropping seasons, wheat root systems were highly colonized by the native mycorrhizal fungi in the WF and W+F treatments compared to the wheat monoculture (Table 1). The development of the mycorrhizal hyphal network was also stimulated in the intercropping and rotation treatments with highest values recorded in the WF treatment after the first growing season (Table 1).

In 2012, the catabolic evenness ranged as follow W+F > W = WF whereas it ranged as W > W+F > WF in 2013 (Table 1). No significant differences were recorded for the average amino-acid induced respiration within the treatments in 2012 and 2013 and for the average carbohydrate induced respiration in 2013 (Table 1). In contrast, the average SIR was significantly higher with carbohydrates in the WF and W+F treatments after the first growing season (Table 1). For the carboxylic acids, highest values were measured in 2012 with the WF treatment and in 2013 with the WF and W+F treatments (Table 1).

After one cropping season, the highest SIRs have been obtained with citric acid, maleic acid, malic acid, oxalic acid, ketoglutaric acid, and ascorbic acid in the WF treatment (Figure 1) whereas and after two cropping seasons, a significant highest SIR response was only recorded with maleic acid in the WF treatment (Figure 2). The BGA permutation test of ISCP profiles in the three treatments (W, W+F, WF) is highly significant (p ~1.10 10⁻⁴) in 2012, but not in 2013. The BGA

permutation test for the comparison of SIR profiles in 2012 vs. 2013 is also significant (p ~0.03). There is no relationship between SIR profiles and yield variables in 2012 or in 2013 (the permutation tests of Co-Inertia Analyses are not significant). But there is a strong relationship between 2012 SIR profiles and 2013 yield variables (p ~0.004; Figure 3). Samples with higher yields are on the left, they correspond to intercropping (WF), and to high consumption of the following organic substrates: trehalose, glucosamine, glucose, sucrose, and organic acids (Figure 3). Lower yield levels are on the right, they correspond to wheat monoculture (W), and to the following organic substrates: galactose, glutamate, asparagine (Figure 3). On the vertical axis, rotation system (W+F) correspond to intermediate yield levels, and to particular organic substrates: gluconic acid, glutamine, xylose, mannose, mannitol, sorbitol (Figure 3).

## DISCUSSION

This study conducted during two growing seasons clearly shows that intercropping could lead to highest performances in terms of crop yield than the rotation and the monoculture after two cropping seasons. Benefits resulting from intercrop seem to be linked with changes recorded on the mycorrhizal soil infectivity and on soil microbial functionalities.

These results are in accordance with previous studies showing that productivity advantages of intercropping may arise from complement use of growth resources such as N and water in either space or time (Akter et al., 2004; Chu et al., 2004). In the same way, intercropping legume and cereals may lead in highest nitrogen content in the cereal grain, improving that quality criterion (Bulson et al., 1997; Gooding et al., 2007).

<footnote n="1">¹Standard error. ²For each year, data in the same line followed by the same letter are not significantly different according to the Newman-Keul's test (p < 0.05). ³nd, not determined.</footnote>
<figure><img src="image_3.png" /></figure>

**FIGURE 1** | *In situ* catabolic potential (ISCP) profiles after the first growing season (2011–2012). Error bars represent standard errors (n = 4). Blue bars: wheat monoculture (W); Red bars: wheat/Faba bean intercropping (WF); Green bars: Faba bean/wheat rotation (W+F). An asterisk indicates a significant difference between the WF treatment and the others according to the Newman-Keul's test (*P* < 0.05).

<figure><img src="image_4.png" /></figure>

**FIGURE 2** | *In situ* catabolic potential profiles after the second growing season (2012–2013). Error bars represent standard errors (n = 4). Blue bars: wheat monoculture (W); Red bars: wheat/Faba bean intercropping (WF); Green bars: Faba bean/wheat rotation (W+F). An asterisk indicates a significant difference between the WF treatment and the others according to the Newman-Keul's test (*P* < 0.05).
<figure><img src="image_5.png" /></figure>

**FIGURE 3** | Between-Group CO-Inertia Analysis (BGCOIA) on patterns of ISCP profiles (2012) and yield variables (2013). The three red ellipses represent the samples of the yield variables table in 2013 for the three treatments (W = wheat monoculture, W+F = wheat/Faba bean rotation, WF = wheat/Faba bean intercropping). The four yield variables (in red) are represented by the four arrows pointing to the left (W1000, weight of 1000 grains, PB, plant biomass, DW, dry weight of spikes, SN, spike number). The three blue ellipses represent the samples of the ISCP profiles table in 2012 for the same three cultivation modes. The names of the organic substrates are given in blue.

The importance of above-belowground interactions in both natural and agricultural systems has been highlighted during the last decade (Wolters et al., 2000; Soler et al., 2012; Orrell and Bennett, 2013). Mixing plant species will create new habitats for associated species and more particularly through its impact on the soil microbiota composition (Bartelt-Ryser et al., 2005). In particular changes in plant cover composition alter the composition of Arbuscular Mycorrhizal (AM) fungal communities (Lovelock et al., 2003). AM fungi facilitate plant uptake and transport of less mobile soil nutrients (Jakobsen et al., 2001), enhance drought tolerance (Kaya et al., 2003) and reduce pathogenic infections (Abdallah and Abdel-Fattah, 2000). These fungal symbionts are also involved in the biological mechanisms that influence plant community productivity and plant-plant interactions (Wagg et al., 2011). It has been also reported that an AM fungal diversity increase could relax competition in species network (Bastolla et al., 2009). Although AM fungi have been traditionally believed to be non-host specific in their ability to infect and to promote the host plant growth, the benefits expected from the mycorrhizal symbiosis to enhance the plant host development may highly depend on the particular

species involved (Bever et al., 2001; Jansa et al., 2005). For instance, it is well known that plant species, highly dependent to the mycorrhizal symbiosis for their growth (i.e., legumes) will promote the development of the mycorrhizal fungal growth, spore production, and hyphal network extend (Duponnois et al., 2001). Our results corroborate these previous studies as higher plant development; better mineral nutrition and higher mycorrhizal development were recorded in the intercropping treatment after one growing season.

Interactions of AM fungi with soil microbiota are an important driver of plant growth (Linderman, 1988; Johansson et al., 2004; Artursson et al., 2006). AM symbiosis is known to promote root exudation (Rambelli, 1973), and to influence rhizosphere microbial communities (Johansson et al., 2004). All these interactions constitute a plant root-mycorrhizal-bacteria continuum named “mycorrhizosphere.” In the current study, the soil compartment can be assimilated to the mycorrhizosphere. After one growing season, soil microbial functionalities were highly impacted by the plant cover composition (monoculture vs. intercropping) highlighting the importance of these aboveground-belowground interactions.
The main differences were observed between the treatments with SIR response to carboxylic acids recorded with the intercropping treatment. AM fungi and their associated bacteria may excrete carboxylic acids as well as faba bean (Zhou et al., 2009; Hafidi et al., 2013). These organic compounds could also exert a selective influence on soil microbial communities by enhancing the multiplication of microorganisms able to catabolise organic acids (Ouahmane et al., 2007; Hafidi et al., 2013).

In the current study, a strong relationship between 2012 SIR profiles and 2013 yield variables was recorded. This result showed that intercropping strongly impacted soil microbial functionalities resulting in a positive effect of this cultural practice recorded again in the second cropping year.

## CONCLUSION

This field experiment demonstrates that the benefits of intercropping are highly subjected to the mycorrhizal symbiosis establishment and its impact on the soil microflora functionalities. In the present study, intercropping gave better results than the rotation practice in terms of crop productivity. However, this study has been only performed during two growing years and other studies have to be undertaken for more than two growing seasons to evaluate the long-term impact of intercropping compared to that resulting from rotations. These results emphasize the need to develop crop diversity

in agroecosystems and to include the management of AM fungal communities in agro-ecological strategies in order to sustainably maintain the crop productivity. Another point of major interest would be to check, in different situations and over several crop rotations, the predictive validity of the SIR approaches as a bioindicator of soil high productivity with wheat and other cereals. In other words, in the absence of any chemical inputs, are SIR values valid indicators of soil resilience, across several successive legume/crops associations? Limiting SIR substrates to carboxylic acids would be a way to simplify and reduce the cost of such a biotest, as to maintain its reliability.

## AUTHOR CONTRIBUTIONS

SW, YP, JT, HS, EB, and RD: These authors contributed to defining the objectives of the experiment, the interpretation of results, and the redaction of the article. TM, KO, CL, AG, and MH: These authors contributed to the implementation of the experiences and proofreading of the article.

## ACKNOWLEDGMENT

This study was funded by the project Fabatropimped supported by Agropolis Foundation under the reference ID 1001-009.

## REFERENCES

* Abdallah, M. E., and Abdel-Fattah, G. M. (2000). Influence of the endomycorrhizal fungus *Glomus mosseae* on the development of peanut pod rot disease in Egypt. *Mycorrhiza* 10, 29–35. doi: 10.1007/s005720050284

* Akter, N., Alim, A. M., Islam, M. M., Naher, Z., Rahman, M., and Hossain, L. M. (2004). Evaluation of mixed and intercropping lentil and wheat. *J. Agron.* 3, 48–51. doi: 10.3923/ja.2004.48.51

* Altieri, M. A. (1999). The ecological role of biodiversity in agroecosystems. *Agric. Ecosyst. Environ.* 74, 19–31. doi: 10.1016/S0167-8809(99)00028-6

* Artursson, V., Finlay, R. D., and Janson, J. K. (2006). Interactions between arbuscular mycorrhizal fungi and bacteria and their potential for stimulating plant growth. *Environ. Microbiol.* 8, 1–10. doi: 10.1111/j.1462-2920.2005.00942.x

* Bartelt-Ryser, J., Joshi, J., Schmid, B., Brandl, H., and Balser, T. (2005). Soil feedbacks of plant diversity on soil microbial communities and subsequent plant growth. *Perspect. Plant Ecol. Evol. Syst.* 7, 27–49. doi: 10.1016/j.ppees.2004.11.002

* Bastolla, U., Fortuna, M. A., Pascal-Garcia, A., Ferrera, A., Luque, B., and Bascompte, J. (2009). The architecture of mutualistic networks minimizes competition and increases biodiversity. *Nature* 458, 1018–1021. doi: 10.1038/nature07950

* Betencourt, E., Duputel, M., Colomb, B., Desclaux, D., and Hinsinger, P. (2012). Intercropping promotes the ability of durum wheat and chickpea to increase rhizosphere phosphorus availability in low P soil. *Soil Biol. Biochem.* 46, 21–33. doi: 10.1016/j.soilbio.2011.11.015

* Bever, J. D., Schultz, P. A., Pringle, A., and Morton, J. B. (2001). Arbuscular mycorrhizal fungi: more diverse than meets the eye, and the ecological tale of why. *Bioscience* 51, 923–931. doi: 10.1641/0006-3568(2001)051[0923:AMFMDT]2.0.CO;2

* Bulson, H. A. J., Snydon, R. W., and Stopes, C. E. (1997). Effects of plant density on intercropped wheat and field beans in an organic farming system. *J. Agric. Sci.* 128, 59–71. doi: 10.1017/S0021859696003759

* Campbell, C. D., Chapman, S. J., Cameron, C. M., Davidson, M. S., and Potts, J. M. (2003). A rapid microtiter plate method to measure carbon dioxide evolved from carbon substrate amendments so as to determine the physiological profiles of soil microbial communities by using whole soil. *Appl. Environ. Microbiol.* 69, 3593–3599. doi: 10.1128/AEM.69.6.3593-3599.2003

* Chu, G. X., Shen, Q. R., and Cao, J. L. (2004). Nitrogen fixation and N transfer from peanut to rice cultivated in aerobic soil in intercropping system and its effect on soil N-fertility. *Plant Soil* 263, 17–27. doi: 10.1023/B:PLSO.0000047722.49160.9e

* Corre-Hellou, G., and Crozat, Y. (2005). Assessment of root system dynamics of species grown in mixtures under field conditions using herbicide injection and 15N natural abundance methods: a case study with pea, barley and mustard. *Plant Soil* 276, 177–192. doi: 10.1007/s11104-005-4275-z

* Duponnois, R., Plenchette, C., and Bâ, A. M. (2001). Growth stimulation of seventeen fallow leguminous plants inoculated with *Glomus aggregatum* in Senegal. *Eur. J. Soil Biol.* 37, 181–186. doi: 10.1016/S1164-5563(01) 01077-9

* Franquet, E., Doledec, S., and Chessel, D. (1995). Using multivariate analyses for separating spatial and temporal effects within species-environment relationships. *Hydrobiologia* 300, 425–431. doi: 10.1007/BF000 24484

* Fridley, J. D. (2001). The influence of species diversity on ecosystem productivity: how, where and why? *Oikos* 93, 514–526. doi: 10.1034/j.1600-0706.2001.930318.x

* Garcia-Barrios, L. (2003). "Plant-plant interactions in tropical agriculture," in *Tropical Agroecosystems*, ed. J. Vandermeer (Boca Raton, FL: CRC Press), 11–58.
* Garofalo, P., Di Paolo, E., and Rinaldi, M. (2009). Durum wheat (*Triticum durum* Desf.) in rotation with faba bean (*Vicia faba* var. minor L.): long-term simulation case study. *Crop Pasture Sci.* 60, 240–250. doi: 10.1071/CP08208

* Giovannetti, M., and Mosse, B. (1980). An evaluation of techniques for measuring vesicular arbuscular mycorrhizal infection in roots. *New Phytol.* 84, 489–500. doi: 10.1111/j.1469-8137.1980.tb04556.x
Gooding, M. J., Kasyanova, E., Ruske, R., Hauggaard-Nielsen, H., Jensen, E. S., Dahlmann, C., et al. (2007). Intercropping with pulses to concentrate nitrogen and sulphur in wheat. *J. Agric. Sci.* 145, 469–479.
doi: 10.1017/S0021859607007241

Gurr, G. M., Wratten, S. D., and Luna, J. M. (2003). Multi-function agricultural biodiversity: pest management and other benefits. *Basic Appl. Ecol.* 4, 107–116.
doi: 10.1078/1439-1791-00122

Hafidi, M., Ouahmane, L., Thioulouse, J., Sanguin, H., Boumezzough, A., Prin, Y., et al. (2013). Managing mediterranean nurse plants-mediated effects on soil microbial functions to improve rock phosphate solubilization processes and early growth of *Cupressus atlantica* G. Ecol. Eng.* 57, 57–64. doi: 10.1016/j.ecoleng.2013.04.006

Hauggaard-Nielsen, H., Ambus, P., and Jensen, E. S. (2001). Temporal and spatial distribution of roots and competition for nitrogen in pea-barley intercrops - a field study employing 32P technique. *Plant Soil* 236, 63–74. doi: 10.1023/A:1011909414400

Hauggaard-Nielsen, H., Jornsgaard, B., Kinane, J., and Jensen, E. S. (2007). Grain legume-cereal intercropping: the practical application of diversity, competition and facilitation in arable and organic cropping systems. *Renew. Agric. Food Syst.* 23, 3–12. doi: 10.1017/S1742170507002025

Hinsinger, P., Bettencourt, E., Bernard, L., Brauman, A., Plassard, C., Shen, J., et al. (2011). P for two sharing a scarce resource: soil phosphorus acquisition in the rhizosphere of intercropped species. *Plant Physiol.* 156, 1075–1086. doi: 10.1104/pp.111.175331

Hooper, D. U., and Vitousek, P. M. (1997). The effects of plant composition and diversity on ecosystem processes. *Science* 277, 1302–1305. doi: 10.1126/science.277.5330.1302

Jakobsen, I., Gazey, C., and Abbott, L. K. (2001). Phosphate transport by communities of arbuscular mycorrhizal fungi in intact soil cores. *New Phytol.* 149, 95–103. doi: 10.1046/j.1469-8137.2001.00006.x

Jakobsen, I., and Rosendahl, L. (1990). Carbon flow into soil and external hyphae from roots of mycorrhizal cucumber plants. *New Phytol.* 115, 77–83. doi: 10.1111/j.1469-8137.1990.tb00924.x

Jansa, J., Mozafar, A., and Frossard, E. (2005). Phosphorus acquisition strategies within arbuscular mycorrhizal fungal community of a single field site. *Plant Soil* 276, 163–176. doi: 10.1007/s11104-005-4274-0

Jensen, E. S., Bedoussac, L., Carlsson, G., Journet, E. P., Justes, E., and Hauggaard-Nielsen, H. (2015). Enhancing yields in organic crop production by eco-functional intensification. *Sustain. Agric. Res.* 4, 42–50. doi: 10.5539/sar.v4n3p42

Johansson, J. F., Paul, L. R., and Finlay, R. D. (2004). Microbial interactions in the mycorrhizosphere and their significance for sustainable agriculture. *FEMS Microbiol. Ecol.* 48, 1–13. doi: 10.1016/j.femsec.2003.11.012

John, M. K. (1970). Colorimetric determination in soil and plant material with ascorbic acid. *Soil Sci.* 68, 171–177.

Kaya, C., Higgs, D., Kirnak, H., and Tas, I. (2003). Mycorrhizal colonization improves fruit yield and water use efficiency in water melon (*Citrullus lanatus* Thumb) grown under well-watered and water-stressed conditions. *Plant Soil* 253, 287–292. doi: 10.1023/A:1024843419670

Kisa, M., Sanon, A., Thioulouse, J., Assigbetse, K., Sylla, S., Spichiger, R., et al. (2007). Arbuscular mycorrhizal symbiosis can counterbalance the negative influence of the exotic tree species *Eucalyptus camaldulensis* on the structure and functioning of soil microbial communities in a sahelian soil. *FEMS Microbiol. Ecol.* 62, 32–44. doi: 10.1111/j.1574-6941.2007.00363.x

Klironomos, J. N. (2002). Feedback with soil biota contributes to plant rarity and invasiveness in communities. *Nature* 417, 67–70. doi: 10.1038/417067a

Köpke, U., and Nemecek, T. (2010). Ecological services of faba bean. *Field Crops Res.* 115, 217–233. doi: 10.1016/j.fcr.2009.10.012

Latati, M., Blavet, D., Alkama, N., Laoufi, H., Drevon, J. J., Gérard, F., et al. (2014). The intercropping cowpea–maize improves soil phosphorus availability and maize yields in an alkaline soil. *Plant Soil* 85, 181–191. doi: 10.1007/s11104-014-2214-6

Latati, M., Pansu, M., Drevon, J. J., and Ounane, S. M. (2013). Advantage of intercropping maize (*Zea mays* L.) and common bean (*Phaseolus vulgaris* L.) on yield and nitrogen uptake in Northeast Algeria. *Int. J. Res. Appl. Sci.* 01, 1–7.

Lebart, L., Morineau, A., and Warwick, K. M. (1984). *Multivariate Descriptive Statistical Analysis*. New York, NY: Wiley.

Li, L., Tang, C., Rengel, Z., and Zhang, F. (2005). Chickpea facilitates phosphorus uptake by intercropped wheat from an organic phosphorus source. *Plant Soil* 248, 297–303. doi: 10.1023/A:1022389707051

Li, W., Li, L., Sun, J., Zhang, F. S., and Christie, P. (2003). Effects of nitrogen and phosphorus fertilizers and intercropping on uptake of nitrogen and phosphorus by wheat, maize, and faba bean. *J. Plant Nutr.* 26, 629–642. doi: 10.1081/PLN-120017670

Linderman, R. G. (1988). Mycorrhizal interactions with the rhizosphere microflora: the mycorrhizosphere effect. *Phytopathology* 78, 366–371.

Lovelock, C. E., Andersen, K., and Morton, J. B. (2003). Arbuscular mycorrhizal communities in tropical forests are affected by host tree species and environment. *Oecologia* 135, 268–297. doi: 10.1007/s00442-002-1166-3

Magurran, A. E. (1988). *Ecological Diversity and its Measurement*. London: Croom Helm.

Malézieux, E., Crozat, Y., Dupraz, C., Laurans, M., Makowski, D., Ozier-Lafontaine, H., et al. (2009). Mixing plant species in cropping systems: concepts, tools and models. A review. *Agron. Sustain. Dev.* 29, 43–62. doi: 10.1051/agro:2007057

Naeem, S., Thompson, L. J., Lawler, S. P., Lawton, J. H., and Woodfin, R. M. (1994). Declining biodiversity can alter the performance of ecosystems. *Nature* 368, 734–737. doi: 10.1038/368734a0

Ofori, F., and Stern, W. R. (1987). Cereal-legume intercropping systems. *Adv. Agron.* 41, 41–90. doi: 10.1016/S0065-2113(08)60802-0

Orrell, P., and Bennett, E. (2013). How can we exploit above-belowground interactions to assist in addressing the challenges of food security. *Front. Plant Sci.* 4:432. doi: 10.3389/fpls.2013.00432

Ouahmane, L., Thioulouse, J., Hafidi, M., Prin, Y., Ducouso, M., Galiana, A., et al. (2007). Soil functional diversity and P solubilization from rock phosphate after inoculation with native or exotic arbuscular mycorrhizal fungi. *Forest Ecol. Manage.* 241, 200–208. doi: 10.1016/j.foreco.2007.01.015

Phillips, J. M., and Hayman, D. S. (1970). Improved procedures for clearing roots and staining parasitic and vesicular-arbuscular mycorrhizal fungi for rapid assessment of infection. *Trans. British Mycol. Soc.* 55, 158–161. doi: 10.1016/S0007-1536(70)80110-3

Rambelli, A. (1973). "The rhizosphere of mycorrhizae," in *Ectomycorrhizae: Their Ecology and Physiology*, eds G. C. Marks and T. T. Kozlowski (New York, NY: Academic Press), 299–343.

Schnitzer, S. A., Klironomos, J. N., Hille Ris Lambers, J., Kinkel, L. L., Reich, P. B., Xiao, K., et al. (2011). Soil microbes drive the classic plant diversity-productivity pattern. *Ecology* 92, 296–303. doi: 10.1890/10-0773.1

Soler, R., van der Putten, W. H., Harvey, J. A., Vet, L. E. M., Dickie, M., and Bezemer, T. M. (2012). Root herbivore effects on aboveground multitrophic interactions: patterns, processes and mechanisms. *J. Chem. Ecol.* 38, 755–767. doi: 10.1007/s10886-012-0104-z

Tang, X., Bernard, L., Brauman, A., Daufresne, T., Deleporte, P., Desclaux, D., et al. (2014). Increase in microbial biomass and phosphorus availability in the rhizosphere of intercropped cereal and legumes under field conditions. *Soil Biol. Biochem.* 75, 86–93. doi: 10.1016/j.soilbio.2014.04.001

Thioulouse, J. (2011). Simultaneous analysis of a sequence of paired ecological tables: a comparison of several methods. *Annals Appl. Stat.* 5, 2300–2325. doi: 10.1214/10-AOAS372

Thioulouse, J., and Dray, S. (2007). Interactive multivariate data analysis in R with the ade4 and ade4TkGUI packages. *J. Stat. Soft.* 5, 1–14.

Thioulouse, J., Prin, Y., and Duponnois, R. (2012). Multivariate analyses in soil microbial ecology: a new paradigm. *Environ. Ecol. Stat.* 19, 499–520. doi: 10.1007/s10651-012-0198-z

Tilman, D., and Downing, J. A. (1994). Biodiversity and stability in grasslands. *Nature* 367, 363–365. doi: 10.1038/367363a0

Tilman, D., Knops, J., Wedin, D., Reich, P., Ritchie, M., and Siemann, E. (1997). The influence of functional diversity and composition on ecosystem processes. *Science* 277, 1300–1302. doi: 10.1126/science.277.5330.1300

Tilman, D., Wedin, D., and Knops, J. (1996). Productivity and sustainability influenced by biodiversity in grassland ecosystems. *Nature* 379, 718–720. doi: 10.1038/379718a0

van der Putten, W. H., Bardgett, R. D., Bever, J. D., Bezemer, T. M., Casper, B. B., Fukami, T., et al. (2013). Plant-Soil Feedbacks: the past, the present and future challenges. *J. Ecol.* 101, 265–276. doi: 10.1111/1365-2745.12054
Vandermeer, J. H. (1989). *The Ecology of Intercropping*. Cambridge: Cambridge University Press.

Wagg, C., Jansa, J., Stadler, M., Schmid, B., and van der Heijden, M. G. A. (2011). Mycorrhizal fungal identity and diversity relaxes plant-plant competition. *Ecology* 92, 1303–1313. doi: 10.1890/10-1915.1

Wang, D., Marschner, P., Solaiman, Z., and Rengel, Z. (2007). Growth, P uptake and rhizosphere properties in intercropped wheat and chickpea in soil amended with iron phosphate or phytate. *Soil Biol. Biochem.* 39, 249–256. doi: 10.1016/j.soilbio.2006.07.013

Wolters, V., Silver, W. L., Bignell, D. E., Coleman, D. C., Lavelle, P., van der Putten, W. H., et al. (2000). Effects of global changes on above- and belowground biodiversity in terrestrial ecosystems: implications for ecosystem functioning. *Bioscience* 50, 1089–1098. doi: 10.1641/0006-3568(2000)050[1089:EOGCOA]2.0.CO;2

Zhou, L. L., Cao, J., Zhang, F. S., and Li, L. (2009). Rhizosphere acidification of faba bean, soybean and maize. *Sci. Tot. Environ.* 407, 4356–4362. doi: 10.1016/j.scitotenv.2009.02.006

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

Copyright © 2016 Wahbi, Prin, Thioulouse, Sanguin, Baudoin, Maghraoui, Oufdou, Le Roux, Galiana, Hafidi and Duponnois. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.