Gone and forgotten: facilitative effects of intercropping combinations did not carry over to affect barley performance in a follow-up crop rotation

Intercropping often leads to improved productivity of individual species compared to monocultures. We have practically little knowledge of facilitation effects in different intercropping systems and their importance in creating soil legacies that can indirectly affect the succeeding crop in a crop rotation through plant-soil feedback (PSF) effects. To test this, we used a two-phased field experiment where we combined intercropping and crop rotation. During intercropping, we grew maize, faba bean, and lupine in monocultures or two-species crop combinations. The following season, we grew winter barley on the soil previously used for intercropping to test PSF effects under field conditions. We found evidence for facilitative effects on aboveground biomass production that were species-specific with faba bean and maize biomass benefitting when intercropped compared to their expected biomasses in monocultures. Lupine, in contrast, performed best in monocultures. After the intercropping phase, total soil mineral nitrogen was higher in legume monocultures creating soil legacies but this did not affect soil microbial parameters and barley biomass production in the follow-up rotation phase. We found support for species-specific positive and negative interactions in intercropping. Our results also demonstrated that soil legacies play no significant role under moderately high nutrient environments.


Introduction
To meet the projected food demand by 2050, agricultural production must increase by 60-110% and this increase should be environment-friendly through reduced usage of synthetic pesticides and fertilizers and increased ecological intensification (Tilman et al. 2011;Wezel et al. 2014). In this regard, Gurr et al. (2016) showed evidence that ecological intensification can be promoted by crop diversification. Through crop diversification, increasing the positive biodiversity effects, that is, higher productivity in mixed cultures than the corresponding monocultures, may help us enhance the ecological intensification. We have demonstrated knowledge where increasing plant species richness has shown to increase multiple ecosystem functions in forest (Huang et al. 2018) and grassland (Isbell et al. 2017) ecosystems. This knowledge may be applied in cropping systems to boost agricultural productivity. However, differences in experimental designs and management practices in both forest & grassland ecosystems and cropping systems make it difficult to apply the knowledge gained from the former to cropping systems. For instance, in cropping systems, different intercropping types (relay-, strip-and mixed intercropping) and planting densities have been used, whereas, in most of the biodiversity-ecosystem functioning experiments, planting densities remain constant whereas the proportion of plant species vary. Further, the cropping systems remain intensively managed as compared to forest and less intensive grasslands. Therefore, we require more evidence from cropping systems on relative biomass production with different crop combinations as regards the strength of facilitative interactions leading to enhanced productivity. In situations where facilitation is particularly strong, one could envisage that its effects may even carry over into a subsequent crop.
In cropping systems, crop diversification can be achieved either spatially by growing more than one cultivar or crop simultaneously in close proximity (intercropping) or temporally by growing different consecutive crops (crop rotation). The positive effects of intercropping are mediated through trait complementarity and plasticity as well as the facilitative effects of interacting plant species (summarized in Li et al. 2014) whereas, for crop rotation, such effects are mediated through indirect feedback interactions (Schnitzer et al. 2011;Mushonga et al. 2020). There is more evidence on the positive spatio-temporal crop diversification on plant productivity (Li et al. 2014;Gaudin et al. 2015;Zhang et al. 2017;Dong et al. 2018), but negative effects have also been observed (Polley et al. 2003;Bukowski et al. 2018). This suggests that such effects appear to be species-specific and to a larger extent depend on soil biotic and abiotic properties as well as environmental conditions (Van der Putten et al. 2013;Craven et al. 2016;Png et al. 2019). For instance, by growing 4 barley cultivars (Prague, Spire, Waggon, and Krystal) and 3 legumes (Trifolium subterranaeum, Ornithopus sativus, and Medicago trunculata) in monocultures and possible intercropping combinations, Darch et al. (2018) showed that, compared to monocultures, barley-legume intercropping resulted in an up to 40% increase in overall biomass production (combined of both crops in intercropping). This increase was dependent on soil P availability, with the highest gain occurring at or below the sub-critical P demand for barley. They further showed that intercropping of different cultivars of barley did not change their productivity compared to when growing in monocultures. Su et al. (2014) showed that even though the total chlorophyll content (chlorophyll a + chlorophyll b) of two soyabean cultivars increased in a relay intercropping with maize, the photosynthetic activity decreased as compared to their monocultures. This decreased photosynthetic activity was attributed to shading effects of maize. In an another wheat/maize relay intercropping system, the SPAD values (measure of leaf greenness) of maize decreased when intercropped with wheat . This suggests that it is not always the bigger plant in the mixed cultures that suppress the growth of the 'subordinate' plant. The underlying mechanisms still need to be identified.
The cornerstone of crop rotation practices lies on the assumption of plant-soil feedback (PSF) effects, that is, a preceding plant alters the soil abiotic and biotic components that may ultimately affect succeeding plant performance (Bever 1994;Ehrenfeld et al. 2005). It has been shown how PSF effects contribute to overyielding in intercropping and the succeeding crops by altering soil microbial communities (Wang et al. 2017(Wang et al. , 2020. As microbiome assemblages in the soil appear to be generally plant species-dependent (Panke-Buisse et al. 2015;Uroz et al. 2019), it is believed that having phylogenetically distinct preceding and succeeding plant species may disrupt the species-specific pathogen accumulation in soil, thereby resulting in negative PSF effects (better plant performance in soil previously grown with different species) (Bever 2003;Miller et al. 2019;Heinen et al. 2020). This ideology is not strongly supported by either empirical (Fitzpatrick et al. 2017;Ingerslew and Kaplan 2018;Kaplan et al. 2020) or synthesis (Mehrabi and Tuck 2015) evidence. For example, Ingerslew and Kaplan (2018) demonstrated using the PSF approach that the succeeding plant biomass (tomato) strongly depended on the identity of 36 plant species that previously trained the soil. However, this effect was independent of phylogenetic relatedness of tomato with the preceding plant species.
These findings urge us to identify optimal crop species combinations in intercropping as well as in the rotation with an overall positive interaction effect on both above-and belowground yields and processes, irrespective of their phylogenetic relatedness. As most of the biodiversity-ecosystem functioning and PSF effects knowledge is derived from noncropping systems, we have limited knowledge if these ecological interactions and underlying mechanisms can also be utilized in cropping systems to enhance productivity through ecological intensification. More specifically, there have been limited attempts to combine spatio-temporal diversity (intercropping together with crop rotation) in cropping systems (Karpenstein-Machan and Stuelpnagel 2000;Scalise et al. 2015;Wang et al. 2017Wang et al. , 2020Kaplan et al. 2020). It is also not clear if different crop species that are performing better when intercropped would also create a positive soil legacy by improving soil nutrient contents, dilution of soil borne pathogens, and increased abundance of mutualists (e.g. arbuscular mycorrhizal fungi) that later would benefit the next crop in the rotation. It has recently been shown that identity of previous crop may lead to changes in AMF communities in soil that may persist over time to affect the follow-up crop (Roy et al. 2021). To fill this knowledge gap, we performed a field experiment comprising of two phases: an intercropping phase followed by a crop rotation phase. The intercropping phase consisted of monocultures and intercrops (a combination of two crops) of maize, faba bean, and lupine. The rotation phase had barley monocultures grown on soils from intercropping phase. The overall aim was two-fold: (1) identify the crop combinations in intercropping with overall enhanced biomass production relative to their expected biomasses in monocultures, and (2) to test the PSF effects of intercropping on soil biochemical parameters and barley biomass production in the rotation. Therefore, we hypothesized that: 1. Compared to monocultures, intercropped species will have a greater aboveground biomass production. 2. Intercropping would alter soil properties and create soil legacies which, in turn, affect the micro-bial parameters & the performance of the next crop in the rotation through PSF effects.

Experimental design and management
The field experiment started in May 2019 in an agricultural field in Lüneburg (53° 12′ N and 10° 22′ E), Germany. The climate is typical of temperate regions with mild summers and cold winters. The daily mean temperature and precipitation are shown in Supplementary Fig. 1. The agricultural field was under conventional practices with 800 kg ha −1 chalk lime, 470 HAS (Harnstoff-Ammonsulfat) solution containing 20% N and 6% S, and 300 kg ha −1 Caralonkali containing 12% P 30% K, 6% Mg, and 4% S applied for summer barley in 2018. Soil was slightly acidic (pH H20 6) and classified as Cambisol and contained around 2.1% total C and 0.2% total N. The experiment comprised of block design in which five blocks were placed parallel to each other and six plots of 2 × 2 m were randomly placed inside each block, yielding a total of five replicates per monoculture and intercrop combination. In each block, plots were 1 m apart from each other to avoid edge effects. The experiment consisted of two phases: an intercropping phase and a rotation phase.

SPAD measurements
Approximately on the 80th day after sowing, leaf greenness was measured as a proxy of chlorophyll content using a SPAD 502 Plus Chlorophyll Meter (SPAD-502-, Minolta Camera, Tokyo, Japan). We chose this time period as all the crop species were fully developed and were in their reproductive phase. For this, we randomly selected five plants from each crop species from both monocultures and intercrops. SPAD values were taken from two youngest yet fully developed healthy leaves at 10 points along the leaf length by avoiding edges and mid ribs. For faba bean and lupine, the measurements were distributed over the leaflets per leaf. Doing this, we had 4500 measurement points.
Harvest and soil sampling On 20th and 30th August 2019, we harvested faba bean and lupine, respectively, whereas maize was harvested on 26th September 2019. This differential harvest date was chosen to allow complete maturity of each crop at harvest. We had initially planned to separate the grain yield from the total aboveground biomass at harvest, but due to a rust pathogen infection on the faba bean, we had to harvest before grain maturity. Since we did not have this separation in the faba bean, we decided in order to be consistent to treat all three species in the same way by measuring aboveground biomass. We randomly harvested 5, 10, and 10 plants of maize, faba bean, and lupine, respectively by cutting stem from soil surface towards the center of each plot to avoid edge effects at plot level. The harvested biomass was dried at 60 °C for 5 days to measure dry biomass and extrapolated to kg m −2 . The total aboveground biomass in intercrops for each species was calculated as difference between observed and expected values in intercrops compared to their respective monocultures to identify either positive or negative effects of intercropping on biomass production. For the expected aboveground biomass estimation for each crop species in intercrops, we halved their respective biomasses in monocultures to correct for planting density using paired monocultures and intercrops per block to account for block effects. After maize harvest on 26th September 2019, all the plots experienced a fallow period of 12-13 days before rotation phase started ( Fig. 1 lower panel). Prior to sowing winter barley in the rotation phase, we collected soil samples for loss-on-ignition and soil mineral N measurements to assess soil legacies created by intercropping phase.
Loss-on-ignition (LOI) was used as a proxy for soil organic matter (SOM). Pre-weighed fresh soil samples were first oven dried at 105 °C for overnight to remove the moisture content. Pre-and post-ignition (500 °C for 24 h) soil weight was recorded. Percent LOI was calculated as below: For ammonium (NH 4+ ) and nitrate (NO 3− ), 5 g of fresh soil were extracted in 20 ml of a 0.01 M CaCl 2 solution. After horizontal shaking for 30 min and subsequent centrifugation (for 5 min at 4500 rpm) and filtration (through a Whatman 595 filter paper), ammonium and nitrate concentrations were immediately determined using ion-selective electrodes (Nico 2000 Ltd, UK).

Phase 2: rotation phase
We grew winter barley (Hordeum vulgare vr. Meridian) in the same plots which were used for the previous intercropping phase to investigate if barley performance is affected by soil legacy created by intercropping through PSF effects. For this, barley seeds were hand sown in rows on 10th and 11th October 2019 at a planting density of 300 seeds m −2 . As it was impractical to hand-sow 30 plots of 2 × 2 m 2 , we reduced the sowing area to 1 × 2 m (2 m 2 ) per plot in the rotation phase. To facilitate sowing, we superficially ploughed all plots (~ 10 cm deep) and barley seeds were placed at 4-5 cm soil depth.

Harvest and soil sampling
Barley was harvested on 27th May 2020 from an area of 0.5 × 0.25 m 2 . After harvesting barley, we randomly took 4 soil cores (4 cm inner diameter and 10 cm depth) from the harvested area by placing the soil cores on the cut stem. This allowed us to collect . 1 Upper panel shows the experimental design and layout with all the crop combinations and planting densities in monocultures and in the intercropping phase and that of barley in the subsequent rotation phase. The lower panel shows sowing and harvesting times for both the phases ◂ soil and barley roots. Four cores were then pooled together to make one composite sample per plot and stored at 4 °C overnight before sieving (2 mm sieve) the next day. After the sieving process, the roots were transferred to 250 ml plastic bottles containing distilled water and shaken overnight to remove soil adhering to roots. Afterwards, roots were carefully washed and stored in 90% glycerol for later counting for root length colonization by arbuscular mycorrhizal fungi.

Microbial biomass and potential enzyme activities
Sieved soil samples were used to measure microbial biomass C and N by chloroform-fumigationextraction with modifications (Vance et al. 1987;Witt et al. 2000). Two sets of subsamples (5 g) were taken from fresh samples. One set was horizontally shaken in 25 ml of 0.5 M K 2 SO 4 for 1 h and thereafter centrifuged for 5 min at 4500 rpm. Subsequently, 3 ml of the supernatant were transferred to another plastic vessel and stored frozen until they were analysed for dissolved organic C (DOC) and total dissolved N (TDN) with a TOC analyser (multi N/C 2100S, Analytik Jena, Germany). The other set of samples was fumigated with 50 ml of ethanol-free chloroform for 24 h. After fumigation, soil extractions and C and N measurements were performed as described above.
Soil microbial biomass C and N were determined as the difference of fumigated and non-fumigated DOC and TDN, respectively. Microbial biomass C and N were corrected by extraction efficiency factors of 0.45 (Vance et al. 1987) and 0.54 (Brookes et al. 1985), respectively.
Potential activities of leucine aminopeptidase (LAP), N-acetyl-ß-D-glucosaminidase (NAG), ß-glucosidase (GLU), and phosphomonoesterase (PHO) were measured fluorometrically according to the method described in Marx et al. (2001) and German et al. (2011). Briefly, 0.5 g of soil was suspended in 50 ml sterile deionized water, homogenized for 1 min in a sonication bath, and aliquots of 200 μl were subsequently pipetted under constant stirring into black 96-well microplates (Puregrade, Germany). Optimal substrate concentrations and incubation times for substrates were evaluated ahead. 50 μl of substrate solution were added to each well, followed by a 120 min incubation in the dark at 20 °C. Fluorescence was measured using a Perkin Elmer EnSpire multiplate reader with an excitation of 365 nm and an emission of 450 nm. Potential enzyme activity was expressed in units of nmol MUB/AMC cleaved g −1 dry soil h −1 .
Barley root length colonization by AMF AMF abundance was determined as root length colonization in percent. Fresh roots stored in 90% glycerol were cut into 1-1.5 cm fragments and cleared in 10% KOH for 20 min in a water bath at 80 °C. Afterward, roots were washed 4 times with distilled water and acidified for 10 min with 1% HCl and placed in a 2% blue ink in 1% HCl for 30 min at 80 °C before clearing them overnight in lactoglycerol (1:1:1) (Phillips and Hayman 1970;Vierheilig et al. 1998). Cleared root fragments were mounted on glass slides and the percent of root length colonization was quantified with the intersection method (McGonigle et al. 1990).

Statistics
All the statistical analyses were performed within R environment (Team 2020) and graphs were prepared with the 'ggplot2' (Wickham 2016) and 'ggpurb' (Kassambara 2020) libraries. The measured variables are presented as means with confidence intervals (CIs) of 95% that were computed by using non-parametric bootstrap resampling with 10,000 iterations. To avoid common statistical errors, we followed the step-wise protocol for data exploration (Zuur et al. 2010 (Mangiafico 2020). We refer to significant differences at the p < 0.05 level but based on recent discussion on the significance and null hypothesis testing using α = 0.05, we refrain from using the word 'significant' and mostly mention the mean differences between the treatments and effect sizes wherever possible (Ho et al. 2019;Rillig et al. 2019).

Aboveground biomass production and SPAD values
Total aboveground biomass was affected in intercrops relative to their monocultures and this effect was crop-specific. Aboveground biomass of faba bean and maize increased when intercropped whereas that of lupine decreased in intercropping irrespective of crop combinations. This increase was 65% and 47% for faba bean when intercropped with maize and lupine, respectively (Fig. 2, Supplementary table 2). Similarly, maize aboveground biomass increased by 135% and 131% in intercropping with faba bean and lupine, respectively. On the contrary, the aboveground biomass of lupine decreased by 28% and 36% when intercropped with maize and faba bean, respectively. At the species level, we found greatest SPAD values for lupine and the values remained similar in both monoculture and when intercropped. On the other hand, faba bean and maize had greater SPAD values in their monocultures compared to their intercropping independent of species combination (Supplementary Figure 2).

Legacy effects on soil properties
The total SOM content was unaffected after intercropping phase (Fig. 3a). The mineral N in the soil, however, showed the legacy effects created from the intercropping phase (Fig. 3b). For instance, compared to maize monoculture, legume monocultures had higher mineral N content in the soil (38% higher in Fb-Fb monoculture and 46% higher in L-L monoculture), whereas that of intercrops was in between.

Microbial biomass and potential enzyme activities
In the rotation phase, microbial biomass C and N remained similar and did not vary depending on the intercropping phase. The same was true for microbial biomass C:N ratios (data not shown). Similarly, the potential activities of four measured enzymes were not dependent on the soil legacies from the intercropping phase (Fig. 4).

Root length colonization by arbuscular mycorrhizal fungi
The percent root length of barley colonized by AMF (as determined by staining technique) was affected by the intercropping phase (Fig. 5). Barley grown in soil previously trained by maize-faba bean (M-Fb) crop combination had the highest root length colonization (61%) followed by faba bean-lupine (Fb-L)

Barley aboveground biomass production
Barley aboveground biomass varied from 0.25 to 0.37 kg m −2 but the soil feedback effects from the intercropping phase had no effect on barley biomass production (Fig. 6).

Positive effects of intercropping on aboveground biomass production are species-specific
In support of the first hypothesis, we showed positive effects of intercropping on aboveground biomass production and such effects were crop-specific and dependent on the exact combinations of species grown together (Fig. 2). Maize benefitted the most from facilitative interactions with the legumes. The mineral N accumulation that occurred only in legume monoculture plots underlines the importance of legume-grass interactions as strong candidates for creating facilitative interactions. Such crop-specific effects of intercropping on aboveground biomass production have previously been shown to be likely mediated by interspecific interactions and soil type (Dissanayaka et al. 2015;Gou et al. 2016;Chen et al. 2019), with the underlying mechanisms varying with crop species identity. In intercropping systems of two crops growing simultaneously in close proximity, multiple scenarios may arise in terms of biomass production. For instance, in intercropping, (1) both crop species may benefit from each other thereby increasing their biomasses, (2) one crop species may benefit without affecting the performance of other species, (3) one species may benefit on the expense of other species, and (4) no benefit of intercropping on biomass production of both crop species. In our study, we found that faba bean and maize had greater aboveground biomass production in in intercropping than their expected biomasses in monocultures (Fig. 2). Such stimulated productivity of maize and faba bean biomass has been attributed to inter-specific rhizosphere interactions, in which, root exudates from maize act as signaling molecules to induce faba bean root nodulation and consequently higher rates of biological N fixation (Li et al. 2016). Maize, on the other hand, gets access to soil nutrients such as N that is spared by faba bean but also to increased P availability through faba bean mediated by rhizosphere acidification (Li et al. 2007;Zhang et al. 2016). It should also be noted that even though maize is found to be benefitting when intercropped with legumes (Sileshi et al. 2008;Chai et al. 2014;Latati et al. 2014), the aboveground biomass of maize was exceptionally high in intercropping in the present study. Along with facilitation, this increase may additionally be attributed to the lowered competition for resources as a result of early harvest of faba bean and lupine than maize as well as reduced maize planting density when intercropped. Such temporal differentiation due to different harvesting period has been shown to significantly contribute to  (Yu et al. 2015;Dong et al. 2018).
Next, we showed that when intercropped with lupine, the aboveground biomass production for faba bean and maize was greater (Fb-L and M-L) than their corresponding monocultures, whereas, the aboveground biomass of lupine decreased as compared to its monoculture. Lupines in nature tend to grow in large stands that dominate the surrounding vegetation, and hence evolutionarilyspeaking our results seem to underline this habit, in that it did not benefit from intercropping. The lower performance of lupine when intercropped also hints toward antagonistic inter-specific interactions negatively affecting lupine growth and is very likely that the competition for resources severely constrained lupine growth when intercropped. This notion is supported by the smaller SPAD values (a proxy for chlorophyll content) for faba bean and maize leaves in the intercropping than their corresponding monocultures whereas SPAD values were the highest and remained similar for lupine in both monoculture and intercropping with maize and faba bean. This finding indicated that energy and resource investments for photosynthetic activity was generally greater for lupine than both maize and faba bean and did not change depending on monocultures and intercropping. On the contrary, smaller SPAD values for maize and faba bean when intercropped may suggest reduction in inter-specific competition for resources and efficient resource investments in biomass production for both maize and faba bean. We are aware that SPAD values are only the measure for leaf greenness and the actual rate of photosynthesis in both monocultures and intercrops may vary. It has been shown that intercropping maize with lupine resulted in higher maize biomass production but there was a tendency of lower biomass for lupine (although not significant) compared to their monocultures (Dissanayaka et al. 2015). In support of our results, Hauggaard-Nielsen et al. (2008) found that narrow-leafed lupine (L.angustifolius L.) performance was lowered in intercropping with barley with a reduction in atmospheric N-fixation from 15 to 5-6 g N m −2 . Further investigations are required to quantify C costs for resource acquisition and biomass production for lupine before adopting lupine as a viable companion crop in intercropping and mixed cultures.  Soil legacies from intercropping phase did not affect soil microbial parameters and barley aboveground biomass production in the rotation phase We showed that variation in soil mineral N was dependent on the intercropping phase, thereby, created soil N legacies (Fig. 3). Mineral N was greater in legume monocultures (both Fb-Fb and L-L) compared to maize monoculture, with in-between effect for intercropped combinations. This is very likely a result of residual N in soil from decomposition of high-N plant residues that was reported previously (Freschet et al. 2012). Lower C:N ratios of legume residues make them faster to decompose by microbes. Legumes are also known to increase soil N availability through rhizodeposition (Fustec et al. 2009) and biological N fixation (Jensen and Hauggaard-Nielsen 2003) that leads to facilitative effects on neighbors (Temperton et al. 2007). Temperton et al. (2007) found that legume presence across a gradient of grassland plant diversity in the Jena Experiment facilitated a grass and a forb species, but the exact effect was largest for the grass, with the forb only increasing leaf N but not growing larger with legumes. Contrary to our expectation (i.e. higher mineral N in intercropped combinations than maize monoculture) we found comparable amounts of mineral N in maize monoculture (M-M) and all the intercrop combinations. This contrasts with previous findings showing that plants in mixed cultures extracted more nutrients from soil than those in monocultures due to complementarity in resource acquisition (Levine and HilleRisLambers 2009;Yang et al. 2013;Hacker et al. 2015;Wang et al. 2015). However, in the present study, this could be an artifact as we only measured the mineral N (NO 3 − and NH 4 + ) after the intercropping phase which does not represent all N pools in soil. Depending on cropping systems, discrepancies in total and (in)organic N pools have been reported suggesting alteration in soil N pools after intercropping (Cong et al. 2015;Wang et al. 2015). Future studies would need to also measure organic N pools and mineralization rates to better understand soil N dynamics in relation to relative importance of organic and mineral N in intercropping settings. Similar to mineral N, we expected that faster decomposition of legume residues would increase the fraction of their residues that becomes a part of the SOM thereby increasing the total SOM content in soil. However, we found that the total SOM remained similar after intercropping phase. This finding suggests that the legume-derived SOM fraction decomposed quickly without affecting the total pool of SOM. This is plausible as one would expect an increased pool of particulate organic matter during early decomposition stages which is characterized by faster decomposition than the mineral associated organic matter pool. We suggest that future studies directly quantify the litter decomposition and its contribution in the formation of stable SOM from different cropping systems.
Soil legacies from intercropping phase neither affected the microbial parameters measured nor the aboveground biomass of barley in the rotation phase, thereby rejecting our second hypothesis. We found that, in the rotation phase, microbial biomass C and N, and their potential enzyme activities remained unchanged, suggesting an absence of strong PSF effects. In agreement with our findings, Wang et al. (2015) showed in a decade long mixed cropping experiment that even though the soil chemical parameters such as soil pH, exchangeable potassium, and cation exchange capacity varied depending upon cropping systems (monocultures versus continuous and rotational mixed cultures), the soil biological parameters such as activities of urease, phosphomonoesterase, and nitrate reductase remained largely unaffected. In an another experiment under rainfed conditions, Scalise et al. (2015) showed that legume-cereal intercropping had rather low impact compared to soil type and environmental factors on succeeding durum wheat productivity. Our findings are in contrast with results from Barel et al. (2019), where the identity of preceding crop affected the microbial biomass in the succeeding cropping phase. These discrepancies may arise from different plant species and nutrient availability in different soil types under investigation. For example, total N and total P content was higher in the present study than that found in Barel et al. (2019), and soil nutrient availability has strong regulation on microbial community composition and their activities (Olander and Vitousek 2000;Bell et al. 2015;Kumar et al. 2018) as well as the PSF effects (in't Zandt et al. 2019;Klinerová and Dostál 2019).
No change in barley aboveground biomass production in the rotation phase suggests that PSF effects are context dependent and edaphic factors may, in part, play a significant role. It further suggests that the generally observed positive plant-microbial interaction and plant performance in nutrient-limited soils may fade with higher nutrient availability. Our results are supported by a recent study (in't Zandt et al. 2019), where PSF effects on shoot biomass production of four grassland species were neutralized under increased nutrient availability. Interestingly, we found variation in the barley root length colonized by AMF but this did not lead to a measurable benefit (higher biomass production) for barley. This is in line with the long-held view that under higher nutrient availability, plants are less dependent on AMF for nutrient acquisition (Treseder 2004;Camenzind et al. 2016). Altered AMF colonization of barley roots may be due to changes in their community composition from intercropping phase. AMF communities have been shown to co-vary with their host plant community composition and diversity (Schmid et al. 2020;Smilauer et al. 2020). Therefore, it is very likely that soil harbored differential AMF communities from intercropping phase, which may have varied in the degree of root colonization potential.

Conclusions
We found evidence for good intercropping species combinations (maize and faba bean) as well as not so effective intercropping combinations (with lupine), with species-specific increases in biomass production in intercropping despite relatively high nutrient content in the agricultural soil. This suggests that inter-specific interactions overwhelmed the soil nutrient availability. Density-dependent relaxation in competition with maize in the intercropped combinations may have further resulted in increased biomass production but this was the case only for faba bean and not for lupine. Although belonging to the same plant functional group (i.e. legume), our study underlines that faba bean and lupine have a different potential in intercropping systems for biomass production. Further, we showed that the feedback effects of intercropping did not lead to improved barley biomass production even if there were changes in residual mineral N after the intercropping phase. These effects were also similar for other biological parameters (microbial biomass and their potential enzyme activity) in the rotation phase. Even though we showed that intercropping did not lead to significant PSF effects in our study, such effects may become important in management practices promoting reduced external mineral inputs or in soils with low fertility. of Education and Research (BMBF) within the 'INPLAMINT -Increasing agricultural nutrient-use efficiency by optimizing plant-soil-microorganisms interactions' project with the grant numbers: 031B0508H and 031B0508F. Authors would like to sincerely thank the editor and three anonymous reviewers for their insightful comments and suggestions.
Author contribution AK and VMT designed the experiment. AK collected field data and led lab measurements. CR and HC measured mineral N, microbial biomass, and potential enzyme activities. SP took SPAD measurements for leaf greenness. AK analyzed all the data and prepared the first draft. All coauthors read and contributed to finalize the manuscript, and gave their permission for submission.
Funding Open Access funding enabled and organized by Projekt DEAL.
Data availability Should the manuscript be accepted for publication, the presented data will be submitted to BonaRes data repository (https:// daten zentr um. bonar es. de/ data-portal. php).

Conflict of interest Authors declare no conflict of interest for this submission.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.