The potential of corn-soybean intercropping to improve the soil health status and biomass production in cool climate boreal ecosystems

Intercropping (IC) is a promising approach used to improve soil health and sustainable crop production. However, it is unknown whether IC improve the soil health status and biomass productivity of crops cultivated in podzols under cool climate in boreal ecosystems. Two silage corn and three forage soybean genotypes were cultivated either as inter or monocrop (MC) treatments in a randomized complete block design. IC resulted in 28% increase in total forage production (FP). A reduction in rhizosphere soil pH (RS-pH) was observed in the IC treatments. Conversely, the rhizosphere soil acid phosphatase (RS-APase) activity was significantly higher (26–46%) in the IC treatments and occurred concomitant with a significant increase in available phosphorus (RS-Pavailable) (26–74%) in the rhizosphere. Furthermore, IC enhanced the active microbial composition and strong positive correlations were observed between RS-Pavailable, RS-APase, microbial biomass and FP; while RS-pH was negatively correlated with FP, RS-APase and RS-Pavailable. These findings suggested silage corn intercropped with forage soybean could be a viable approach to enhance FP through improved active microbial community structure, RS-APase activity and RS-Pavailable when cultivated on podzols in cool climate boreal ecosystem.

The global population is anticipated to increase by 33% from 7.2 billion to 9.6 billion by 2050 1 . As such, the demand for sustainable food production systems to support this expanding population under different climate change and land use scenarios will be a challenge to food security 2 . This is of particular relevance in cool climatic regions where the current population are expanding, but there are limitations to increasing agriculture productions in these regions due to short growing season, limited agricultural crop varieties suitably adapted for production in these regions, low availability of arable land and poor soils (podzols) 3 . Therefore, there is a need to have innovative or improved cropping systems that will permit intensive, but sustainable crop production to provide food for livestock and the expanding population in cool climatic regions across the globe.
Intercropping (IC) is defined as the growing of two or more crops simultaneously on the same piece of land [4][5][6] , and is increasingly being adopted as a more sustainable practice in modern agricultural production systems throughout the world [7][8][9] . Conversely, monocropping [MC] (growing of one crop on the same piece of land) is less sustainable though the dominant production system in industrial agriculture. A common and popular index used to evaluate the yield advantage or agronomic performance of IC over MC is the land equivalent ratio (LER) [10][11][12][13][14] . The LER is defined as the relative land area required in MC to produce the same yield as in an IC production system. One advantage of LER is that it can be used to assess the agronomic performance, as well as determine the

Results plant performance indicators (forage production, chlorophyll contents, and plant height).
Total eleven experimental treatments for two corn and three soybean genotypes is presented in Table 1. Out of total 37 observed PLFAs (Table 2), 27 were used to as biomarkers to identify the various microbial groups living under prevailing climatic conditions (Fig. 1). The forage production (FP), chlorophyll contents and plant height were used as indicators of agronomic performance of the cropping systems evaluated in this study. The results showed that all the agronomic parameters except the soybean chlorophyll contents were significantly affected by IC corn with forage soybeans (Tables 3-5).
In general, the chlorophyll contents of corn intercropped with soybean significantly increased compared to corn MC, while a reduction in the soybean plant chlorophyll contents was observed after IC with corn (Table 3). Upright soybean (US) varieties had significant decreased chlorophyll contents when intercropped with corn compared to the vine type soybean (VS) chlorophyll contents. The highest corn chlorophyll contents were measured in S2C2 (46.67), while the lowest was measured in C2 (41.17). For soybean, the highest chlorophyll contents were observed in S2 (31.66), while the lowest was observed in S2C1 (27.34). The overall average corn plant height was significantly different between the monocropped and intercropped treatments (Table 4). However, we observed a significant reduction in the soybean plant height when both upright and vine type soybeans were intercropped with corn (Table 4). In general, the plant height in upright soybeans MC was 54 cm compared to 49 cm in the IC. Conversely, the plant height of S3 (vine soybean) cultivated as MC was 125 cm compared to 112 cm in the IC treatments ( Table 4). The trend was opposite for corn plants when compared to soybeans. IC significantly (p < 0.05) increased the corn plant height as compared to MC treatments ( Table 4). The highest and lowest corn plant heights were recorded in S1C1 (206.57 cm) and C2 (181.80 cm), respectively during the growing season.
In general, IC treatments significantly increased the total FP as compared to corn and forage soybean MC (Table 5). Overall, IC treatments produced 28% higher FP compared to corn MC treatments. It is important to note that the S3 treatment is a vine soybean that were intercropped with silage corn. Overall, IC increased the corn FP but decrease the FP of soybean during the growing seasons ( Table 5).
The highest FP was observed in the S2C1 (15.43 Mg ha −1 ) IC treatment, while the lowest yields were observed in C2 (11.09 Mg ha −1 ) and S3 (1.05 Mg ha −1 ) MC treatments. The LER values were greater than 1 for all IC treatments compared to MC treatments during both years (Supplementary Table S2), that means, more land is required in MC to produce the same crop yield as was obtained in the IC treatments. The observed LER values for the IC treatments ranged from 1.49-1.58 (Supplementary Table S2).
Rhizosphere soil acid phosphatase activity. The (RS-APase) activity was significantly higher (p < 0.05) in the IC treatments compared to MC (Fig. 2). However, no significant difference was observed in (RS-APase) activity between upright (S1 and S2) and vine (S3) soybean varieties in either IC or MC treatments (Fig. 3). The highest RS-APase was observed in S2C1 (70.13 µmol pNP g −1 soil 30 min −1 ). Conversely, the lowest RS-APase was recorded in C2 (34.90 µmol pNP g −1 soil 30 min −1 ). changes in rhizosphere available p and pH. The changes in RS-P available and RS-pH due to different IC treatments are given in Table 6. There was no difference in RS-P available or pH between the VS and US cultivated as MC. However, when both the US and VS were intercropped with corn, there was a significant increase in the content of RS-P available compared to the MC treatments.  Table 2. Phospholipid fatty acids (PLFA) biomarkers used to characterize the active microbial community structure. G+: gram positive bacteria; G−: gram negative bacteria; F: fungi; P: protozoa.
The increase in RS-P available was significantly greater in the US varieties intercropped with corn, than when IC was done with the VS. Furthermore, the increase in RS-P available was significantly (p < 0.05) different for the IC treatments compared to the corn MC treatments. The highest RS-P available was measured in S1C2 (91.44 mg kg −1 ), and the lowest in C2 (41.27 mg kg −1 ). Similarly, the highest P availability was observed in S1C2 (102.71 mg kg −1 ), but the lowest was in S2 (37.66 mg kg −1 ).
Overall corn and soybean IC decreased the RS-pH as compared to when cultivated as MC, but the decrease was not significant. Similarly, there was no significant difference observed between RS-pH of US or VS cultivated as monocrops. The highest pH was measured in S2 (5.89), and the lowest in S1C1 (5.27) of all treatments evaluated in this study (Table 6).
Soil microbial community composition. The effect of corn-soybean IC on rhizosphere active soil microbial community structure (gram-positive (G+) or gram-negative (G−) bacteria, fungi, and protozoa) was assessed using microbial membrane PLFAs (Tables 2 and 7). The total bacterial population was significantly higher in all treatments than the fungi and protozoa population. Both G+ and G− bacteria contributed (48% and 52%, respectively) to the overall observed total bacterial population. Overall, the active microbial population (G+, G− and protozoa) was significantly (p < 0.05) higher in IC compared to the corn and soybean MC; except for the fungal population which was not significantly different than that of the corn MC treatments ( Table 7). The overall trend showed that the microbial community biomass was higher when corn was intercropped with upright soybean varieties than when intercropped with vine soybean. The highest G+ and G− bacterial populations were found in S1C1 (27.83 nmol g −1 ), while the lowest population was found in S2 (20.66 nmol g −1 ) treatment. The fungi and protozoa populations also followed similar trends. The highest populations were observed in S1C1 (7.45 and 1.86 nmol g −1 ), while the lowest populations of both fungi and protozoans were observed in S3 (4.57, 1.52 nmol g −1 ) treatment respectively. The maximum G+: G− and fungi: bacteria ratios were recorded in C2, but overall these ratios were not significantly different between the MC and IC treatments ( Table 7). The total PLFAs content indicate the overall active microbial population was higher in IC treatments compared to corn and soybean MC treatments.
Relationship between soil health status and agronomic performance following corn-soybean cultivation under cool climatic conditions in boreal ecosystem. The active microbial community structure, RS-pH, RS-APase activity, and RS-P available were used as indicators to assess the soil health status; while  Table 8).
Several significant associations were observed between the active microbial community structure, soil biochemical properties and agronomic performance (Figs 5-7; Table 8). These associations were very influential in clustering the MC treatments in different quadrants of the biplot from the IC treatments following redundancy analysis (Fig. 4). For example, RS-APase, RS-P available and FP were the most influential factors clustering the IC together (Fig. 4). We also observed that vine type forage soybean intercropped with silage corn [Inter-(VS)] clustered in separate quadrants of the biplots compared to US intercropped with silage corn [Inter (US)]. The output from the RDA analysis revealed that axis 1 and axis 2 explained 52.55% and 22.97% of the total variance (Fig. 4). Significant positive correlations were observed between the soil microbial PLFAs biomarkers (G+, G−, P, F, T PLFAs), RS-P available and FP. Conversely, all these parameters were negatively correlated with RS-pH (Figs 5 and 6; Table 8). Furthermore, the G+: G− and F:B ratios were negatively correlated with RS-APase, RS-P available and FP; except for the F:B ratio which showed a positive correlation with FP (Fig. 4). Further analysis by Pearson correlation confirmed the significant associations between the active soil microbes, RS-P available , and FP (Figs 5 and 6, Table 8). For example, FP was observed to be significantly associated with not only RS-APase activity and RS-P available but also the following estimated soil microbial PLFAs biomarkers (B, P, F:B, G+: G−, total PLFAs). However, FP was inversely associated with the RS-pH. A consistent pattern also emerged that indicates the RS-APase activity was significantly associated with the RS-P available (r = 0.57) of podzols present in cool climate boreal ecosystems. The RS-P available was also observed to be significantly associated with most of the active soil microbial population present in the root rhizosphere, except for fungi ( Table 8). The agronomic parameters (plant height, chlorophyll content and forage production) were also significantly, and positively correlated (r = 0.87-0.94) with each other (Fig. 7). Specifically, the chlorophyll content was highly associated with increase forage production. This trend was similar between plant height and total forage production in the inter and mono cropping systems evaluated in this study.

Discussion
The growing seasons in cool climate boreal environments are short and characterized by low growth temperatures averaging between 14-17 °C 3 , consistent with the low growth temperatures observed in Newfoundland during the study (Fig. 1). In the current study, we attempted to evaluate the effects of intercropping silage corn with vine type or upright forage soybeans on the crop agronomic performance and soil health status following field cultivation under cool climatic conditions in boreal ecosystems. Intercropping (IC) can enhance crop growth   and productivity due to superior utilization of above and below ground resources 4,44 . Plant chlorophyll content is an important indicator of agronomic performance (enhanced crop growth and productivity) and play a key role in plant photosynthesis 45 . Chlorophyll content and plant height of corn plants increased (by 8.5% and 4% respectively) in the IC treatments compared to MC (Tables 3 and 4). Conversely, the soybean chlorophyll content and plant height decreased (−9.4% and −9.7% respectively) in the IC treatments compared to the MC treatments (Tables 3 and 4). These findings are in agreement with previous studies reported in the literature [46][47][48] , and suggest that both the chlorophyll content and plant height could be the possible drivers for the changes observed in forage production. This is consistent with the high correlations observed between chlorophyll content, plant height and total forage production ( Fig. 7). IC can result in either increased forage production (LER > 1), decreased forage production (LER < 1), or have no effects on forage production (LER = 1) 12,49 . In our experiments, IC increased the FP compared to either corn or soybean MC following cultivation on podzols under cool climatic conditions in boreal ecosystems. The increase observed in FP suggested that IC has an advantage over MC in terms of plant growth, which is also supported by the LER values which were greater than 1 in all IC treatments; further confirming that corn and soybean IC in cool climate boreal environment production systems was superior compared to MC in regards to increasing FP (Supplementary Table S2). Our LER findings are consistent with the results obtained by several other researchers who reported LER values greater than 1 in IC following crop production under warmer climatic conditions 14,23,50 .
The LER values in our study suggest that 49-58% more land would be required for the MC production system to produce a crop yield equal to that of the IC production system in boreal ecosystem. These findings indicate that the IC system appears to use land resources (nutrients) more efficiently than the MC system 17 , and this may account for the superior agronomic performance (enhanced forage production) observed when the crops were cultivated under cool climatic conditions in boreal ecosystem.
The study results demonstrated that IC increased the RS-APase activity between 46% and 26% compared to corn and soybean MC treatments, respectively (Fig. 2). These findings suggest IC may utilize organic P more efficiently than MC corn when cultivated on podzols under cool climatic conditions in boreal ecosystem. Organic P comprises 30-80% of total P in most agricultural soils, and can be converted into RS-P available forms after hydrolyzation by phosphate enzymes 51,52 , consistent with the high association observed between RS-P available and the RS-APase activity (Figs 4 and 5) in this study. The higher RS-APase activity in IC could also be attributed to compatibility or suitability of the silage corn and soybean combinations as companion plants in an IC production system under cool climatic conditions. In arid soils, the RS-P available was significantly correlated with RS-APase activity 53 because of an association between mobilization of organic P and RS-APase activity 54    www.nature.com/scientificreports www.nature.com/scientificreports/ are similar with this and other studies, which reported that cereal-legume IC increased the RS-APase activity compared to when either silage corn or forage soybean were cultivated as monocrops 7,47 . The legume species have been considered as the major contributor to the increased RS-APase activity observed in IC due to the fact that large amounts of acid phosphatase are known to be released from their roots into the root rhizosphere 26,56 . This may account in part for the increase forage production observed in IC considering the significant correlations observed between RS-APase activity, the RS-P available and forage production (Fig. 5).      www.nature.com/scientificreports www.nature.com/scientificreports/ Rhizosphere acidification have been previously reported in several IC systems following cultivation in warm climate production systems 7,12,28,57 . Similar to these findings, we observed that IC decreased the RS-pH compared to when corn and soybean were cultivated as monocrops under cool climatic conditions in boreal ecosystem ( Table 6). The rhizosphere acidification was shown in previous studies to be due to the release of large quantities of protons or organic acids in the root rhizosphere from the crop roots during IC 28,58 . Similar to the RS-pH, IC had positive effects on RS-P available when silage corn and forage soybeans are cultivated as the companion crops in cool climate boreal ecosystem. In fact, IC increased the RS-P available between 74% and 26% in the plant root rhizosphere compared to when corn and soybeans were cultivated as monocrops (Table 6). Consistent with our findings, increased P availability in the rhizosphere have also been reported in intercropping garlic-cucumber 59 ; and maize-chickpea 26 . Increase in RS-APase activity and acidification of the rhizosphere via release of protons and organic acids have been suggested to be responsible for the enhanced RS-P available observed in IC production systems 28 ; and maybe related to the significant inverse relationships observed in our study between RS-P available or RS-APase activity and RS-pH (Fig. 5). www.nature.com/scientificreports www.nature.com/scientificreports/ The active microbial community was investigated by performing PLFAs analysis present in the root rhizosphere following corn and soybean cultivation as either mono or intercrops in podzolic soils under cool climatic conditions in boreal ecosystem. PLFAs are present in the membranes of living cells, but not in dead cells because of rapid degradation during cell death. As such, they can give an accurate estimate of the living (active) microbial community present in the root rhizosphere, and how these community composition change in response to factors such as crop management systems, environmental conditions, and production inputs 60 . Different microbial www.nature.com/scientificreports www.nature.com/scientificreports/ groups present in the soil are comprised of phospholipid fatty acids in their membranes that are diagnostics of their presence and rate of change in the soil habitat 61 . Consequently, the diversity of the active microbial community is referred to as an imperative indicator of soil quality or the health status of the soil 62 . As such, a combination of PLFA markers were used consistent with convention in the literature to delineate changes in the microbial community composition or structure in response to the cropping systems evaluated in this study. Our results showed that IC in general, increased the total microbial PLFAs biomarkers (G+, G−, F and P) in the root rhizosphere as compared to soybean and corn MC. This increase was 12.3% and 12.3% (Table 7). Our results corroborate the findings of Li et al. 41 and Zhou et al. 4 who demonstrated that IC can enhance both bacterial and fungal populations in the plant root rhizosphere. However, these studies were not conducted on crops cultivated in cool climate production systems. To our knowledge, this is the first work demonstrating the potential of intercropping to modulate the active microbial population in podzolic soils present in cool climate boreal ecosystems. We observed that the bacterial population was the highest of the total PLFAs biomarkers present in corn and soybeans cultivated as MC or IC under cool climatic conditions in boreal ecosystem (Table 7). However, the G− population was 4% higher than the G+ bacterial population (Table 7). Higher proportion of G+ bacterial population compared to G− suggested a deficiency of organic carbon in the soil 63,64 . In contrast, the dominance of G− bacteria over G+ bacteria in the soil, and a high fungal population is characteristic of the presence of higher amount of complex organic matter in the soil 63,65 . The findings from this work demonstrate that corn soybean IC promoted the growth and diversity of the active microbial community, and as such can enhance the soil health status under cool climatic conditions in boreal ecosystems.
The significant positive correlation observed between the agronomic and soil chemical properties demonstrated that these parameters are associated with the superior FP (Figs 6 and 7; Table 8) observed when corn was IC with soybean and cultivated under cool climatic conditions in boreal ecosystems. The FP revealed a significant positive correlation with G+, G− bacteria, and the total bacterial population (Fig. 6) consistent with observations in earlier findings 12 in temperate ecosystems. The increase in above ground biomass is known to have a strong positive connection with the below ground soil health indicators 12,25,29,[66][67][68] . This was evident in the strong associations between shifts in the active microbial populations (Figs 4-6), and increase forage production (Figs 5 and 6), RS-P available (Fig. 5, Table 8) observed following intercropping forage soybean with silage corn under field conditions in cool climate boreal ecosystem. Collectively, reduction in RS-pH, increased RS-APase activity, RS-P available, as well as increased in the active fungi, protozoan and bacterial population appears to be the most important determinants of forage production (Figs 4-7), when corn and soybeans are cultivated as intercrops under cool climatic conditions in boreal environment. Similar relationships have been reported between agronomic performance and soil health indicators in various IC production systems under different climatic conditions 12,25,29,[66][67][68] . This is the first study demonstrating that similar relationships exist when silage corn is intercropped with soybeans (vine or upright varieties) and cultivated on podzols under cool climatic conditions in boreal ecosystem. Soil microbes are known to mineralize organic matter and other sources of plant nutrients located in the soil, thus making them available to the plant for uptake, growth and productivity 12 . This is consistent with the significant correlations observed between RS-A vailable , and the active soil microbial population (Table 8). Similarly, IC can  Table 8. Pearson correlation coefficients between phospholipid fatty acids (PLFAs) and rhizosphere soil acid phosphates activity (RS-APase) and rhizosphere soil available phosphorus (RS-P availlable ). ns: correlation is nonsignificant; *Correlation is significant (p ≤ 0.05); (**p ≤ 0.01); (***p ≤ 0.001). n = 66. www.nature.com/scientificreports www.nature.com/scientificreports/ stimulate the enrichment of P solubilizing soil microbes or microbial species with enhanced soil phosphatase activities, thereby increasing the RS-P available during IC 12,29,68 . Thus, it appears, under cool climatic conditions in boreal environments, that the increased forage production in the IC production system was highly dependent on RS-P available and RS-APase activity, presumably through stimulation or modification of the active microbial community structure (Fig. 5, Tables 6-8). The enhanced microbial population observed in IC (Table 7) might be more efficient in mineralizing and mobilizing P in the root rhizosphere under cool climatic conditions in boreal ecosystem. This could be the mechanism through which the improved agronomic performance observed in the tested IC system may be related to the enhanced soil health status. Further experimentation at the molecular genetics and cellular levels are needed to confirm this mechanism and are the subject of future work in our research program.
conclusion Consistent with the objective of this study. We observed that forage soybeans intercropped with silage corn resulted in significantly enhanced agronomic performance and forage production in cool climate boreal ecosystem. In general, corn intercropped with US genotypes displayed superior agronomic performance compared to when intercropped with VS.
Collectively, reduction in RS-pH, and increased RS-APase activity, RS-P available , and the active fungi, protozoan and bacterial populations appear to be the most important determinants of the soil health status and improved forage production, when silage corn and forage soybeans are cultivated as intercrops in podzols under cool climatic conditions in boreal environment.
This study is the first to demonstrate that intercropping silage corn with forage soybeans is a suitable approach to increase forage production and enhanced the soil health status under cool climatic conditions in boreal ecosystems. This work will provide significant improvement in our knowledge to better understand agriculture production in boreal ecosystems or northern climates, particularly in the context of climate change and expanding global populations in these geographic regions where food security is anticipated to be a challenge in the future 1,2 . The sowing of vine type soybean (S3) was carried out with a hand drill due to the small seed size with same line spacing (almost 1 m) as done with SAMCO system for other corn and soybean varieties. The silage corn genotypes were selected based on low corn heating units requirements 69 . There was a total of eleven treatments for the two corn and three soybean genotypes cultivated as either mono or intercrops ( Table 1).

Methods
The following seeding rates were used for MC (corn: 77,100 seeds ha −1 ; soybean: 129,200 seeds ha −1 ) and IC (60% corn + 40% soybean; total 129,200 seeds ha −1 ) during both study years. Crop nutrient requirements were fulfilled through inorganic fertilizers using the regional recommended rates for MC or IC based on the soil nutrient status prior to planting (Supplementary Table S1).
Soybean seeds were inoculated with Bradyrhizobium japonicum @ 10 g kg −1 seeds 70 before seeding. Herbicide application was carried out using roundup WeatherMax ® (Monsanto Canada Inc) during both growing seasons for weed control. The crop was harvested on October 25 th and 13 th during 2016 and 2017, respectively. The experiment was laid out in a randomized complete block design with three replications per treatment. Each experimental treatment plot was 5 m × 6 m in dimension. Weather data for both seasons were collected from a weather station located adjacent to the experimental plots and are reported in Fig. 1.

Crop agronomic performance: Chlorophyll contents, final plant height, and forage production.
The chlorophyll contents were measured using a portable chlorophyll meter (SPAD-502 Konica-Minolta, Japan) taken from the top three leaves of the corn and soybean plants at 65 and 77 days after sowing during both growing seasons. At the corn physiological maturity (R6), four plants in a transect were selected from the replicates of each experimental treatment, and the plant height measured. The same plants were then uprooted and separated into roots and shoot to measure the biomass production. Plant fresh weight was recorded, and a subsample was taken from each treatment to measure the dry matter percentage by drying in a forced air oven (Shel Lab ® ) at 65 °C for 72 h. Thereafter, the total forage production was calculated considering the dry matter percentage and total fresh biomass per treatment.
Soil health status evaluation. To quantify the effects of IC and MC on the soil health status, soil samples were collected from the root rhizosphere to measure the active microbial community composition by analyzing the microbial membrane phospholipid fatty acids (PLFA), along with soil pH, soil available P and soil acid phosphatase activities (i.e. the rhizosphere active microbial community composition, soil available P, soil pH and acid phosphatase activity were used as indicators of the soil heath status in this study). At harvest, plants were uprooted gently, and soil samples near the roots were collected, additionally roots were shaken gently to collect all soil attached to the root surface to obtain the rhizosphere soil for further analysis. The rhizosphere soil was sieved through 2 mm meshes to remove plant roots, small stones, gravel etc. Aliquots of the fresh soil (4 g) was used for PLFA analysis, and the rest of the soil stored at −20 °C for testing RS-P available , RS-APase activity and soil pH.
(METTLER TOLEDO, Canada) 71 . Briefly, 10 g of air-dried, sieved soil (2 mm) was weighed in 50 mL polypropylene centrifuge tubes, and 20 mL of 0.01 M CaCl 2 was added. The soil solution was then mixed for 30 min on an orbital shaker (Innova ™ 2300 Platform Shaker, New Brunswick Scientific, USA) at 120 rpm, then allowed to stand for 1 h before measuring the RS-pH.
Rhizosphere soil acid phosphatase activity. RS-APase activity was measured using the modified methods of Tabatabai and Bremner 72 . Briefly, 1 g of 2 mm sieved soil was weighed and extracted in 1 mL of 0.09 M (pH 4.8) citrate buffer. Polypropylene centrifuge tubes containing soil and citrate buffer were then centrifuged (Heraeus ™ Megafuge ™ 16 Centrifuge Series) at 5000 rpm for 10 min. An aliquot (50 µL) of the supernatant was collected and RS-APase activity assessed after incubating for 30 min in the oven at 37 °C with 1 mM of 4-nitrophenyl phosphate (pNP) and 50 µL citrate buffer. The reaction was terminated immediately after incubation with 20 µL of 0.5 N sodium hydroxide (NaOH). The absorbance of the mixture was recorded at 405 nm using a spectrophotometer (BioTek ™ Cytation ™ 3 imaging reader, BioTek, VT USA.) and the RS-APase activity presented in μ mole pNP g −1 soil 30 min −1 .
Rhizosphere soil available phosphorus (RS-p available ). RS-P available was analyzed using the Mehlich-3 extraction method 73 . Briefly, 2 g air dried soil was weighed in 50 mL Erlenmeyer flasks, and 20 mL of Mehlich-3 extractant solution (1:10 soil: extractant) was added. Flasks containing the mixture were shaken for 5 min on an orbital shaker at 120 rpm (Innova ™ 2300 Platform Shaker, New Brunswick Scientific, USA), and the filtrate recovered following filtration using Whatman 42 filter papers (Sigma Aldrich, ON. Canada). Aliquots of the filtrate was then analyzed using an AA3 Continuous Flow Analytical System (AA3HR, SEAL Analytical USA) to measure the phosphate content, which was then converted to total phosphate in the soil sample based on soil weight as follows: (Sigma-Aldrich, ON. Canada) into neutral lipids, glycolipids, and phospholipids using (2.5 mL) chloroform, (4 mL) acetone and (2.5 mL) methanol, respectively. The phospholipid fractions were re-dissolved into 500 µL of methyl tert-butyl ether (MTBE), and aliquots (100 µL) of the phospholipid fractions derivatized using 50 µL trimethyl sulfonium hydroxide (TMSH) in 2 mL GC vials 75 . The mixture in the vials were vortexed and incubated for 30 min at room temperature. After incubation, 10 µL of the internal standard methyl nonadecanoate (C19:0 @ 160 µg/mL) were added to the samples in the vials, and the samples analyzed via gas chromatography-mass spectrometry (GC-MS) and gas chromatography-flame ionization detection (GC-FID).

Gc-MS/fiD analysis of soil microbial pLfAs. GC-MS/FID analysis was conducted on a Thermo
Scientific Trace-1300 gas chromatography (GC) coupled to a Thermo Scientific TSQ 8000 Triple Quadrupole mass spectrometer (MS) and a flame ionization detector (FID). GC-MS was used for peak identification, while GC-FID was used for quantification. Methylated fatty acids were separated with a DB23 high resolution column (30 m × 0.25 mm × 0.2 μm; Agilent Technology, Mississauga, Canada) using helium as the carrier gas at a flow rate of 1 mL min −1 . One (1 μL) of each sample was injected in split less mode using a Tri-plus auto-sampler. The oven temperature was programed as follows: the initial oven temperature of 50 °C was held for 1 min, then programmed to increase at 20 °C min −1 to 175 °C, held for 1 min at 175 °C, then increased at 4 °C min −1 to 230 °C, where it was held for 5 min. The methylated PLFAs were identified through retention times comparison and mass spectra obtained from commercial standards (NIST database) (Thermo Scientific, ON. Canada; Supelco 37 Component FAME Mix, and Bacterial Acid Methyl Ester (BAME) Mix obtained from, Sigma Aldrich, ON, Canada). Quantification of individual PLFAs was done using standard curves prepared from the standard mixtures, and values presented as nmol g −1 soil. A total of 37 PLFAs were identified ( Table 2) and 27 of them used as biomarkers to assess different microbial groups living in the soil (active microbial community composition) at the time of sampling.
calculations and statistical analysis. The LER is the relative land area needed for MC to produce the same yield attained by IC 76 . LER was measured to evaluate the effect of IC verses MC 12,76,77 as given in following Equations: where, L C and L S are the partial LER for intercropped corn and soybean, respectively. When the LER value is greater than 1, it indicates that an advantage is gained from IC compared to monocrop cropping system in terms of the use of environmental resources for plant growth and yield (total production). When the LER is equal to 1, it means IC has no advantage over MC in the use of environmental resources; and when the LER is less than 1, it means MC use resources more efficiently than IC for plant growth and yield 12,49,77 . All statistical analyses were conducted using XLSTATS (Addinsoft Inc, Paris, France) and Statistix-10 software programs (Analytical Software, FL, USA), while graphs were created with Sigma Plot 13.0 (Systat Software Inc., San Jose, CA). All measurements of chemical parameters (PLFAs, dry matter, plant height etc.) were made in quadruplet. Analysis of variance (ANOVA) was used to determine the effects of treatments on chemical parameters. Where treatment effects were significant, the means were compared using Fisher's LSD test at α = 0.05. The effects of IC or MC treatments on plant agronomic performance (biomass production or forage yield, plant height, chlorophyll content etc.), and soil health status indicators (soil pH, RS-APase, activities et) were evaluated in this study. Redundancy analysis (RDA) and Pearson's correlation coefficients (r) were used to test the linear relationships between FP, (RS-APase) activities, RS-AP, active microbial composition and soil pH.