Intercropping Forage Mulberry Benefits Nodulation and Growth of Soybeans

1. Introduction

Soybean (Glycine max L.) is a highly versatile crop used in food products, oils, proteins, and animal feed[1], with a low yield in China[2] due to a restricted acreage of suitable arable land[3] and inherent intolerance to continuous cropping[4]. Therefore, there is a need to optimize cultivation practices to enhance soybean production and reduce its import dependency. Intercropping is a widely used agricultural practice that enhances land productivity[5,6]. Intercropping gramineous plants with legumes can reduce nitrogen depletion in soil, improve nodulation for N fixation, and increase N transfer[7], thereby decreasing the need for chemical fertilizers[8]. Mechanistic studies on the promotion of N fixation in legumes by intercropping with gramineous plants mainly focus on the absorption of nitrogen in gramineous plants to create nitrogen-depleted conditions and prompt soybean roots to release nodulation signal substances like flavonoids[7,9,10]. Intercropping also aids rhizosphere soil acidification that benefits the activation of soil phosphorus[11] and strengthens the expression of genes related to nodulation and N fixation[12]. Intercropping soybeans with corn is especially widespread in China and enables farmers to obtain an additional soybean harvest while not impacting the yield of corn[5,13]. However, during the symbiotic stage, soybeans endure considerable shade stress imposed by corn, which reduces their photosynthetic capacity and biomass[14,15] and alters the red-to-farred light ratio in their canopy[16-18], leading to their collapse [19] and reduced yield [20-22].

To mitigate the shading impact on intercropped soybeans, forage mulberry, which can be harvested 2 to 3 times annually, is usually intercropped with soybeans to obtain more time and space for sunlight absorption. Mulberry also serves as an ideal special livestock feed when it is less than 1 meter high and has unlignified tender branches and leaves. Intercropping with soybeans also enhances soil nutrient availability and the diversity of the microbial community, thereby increasing the biomass of mulberry[12]. However, the effects of intercropping forage mulberry on the growth of soybeans have not been comprehensively explored. Therefore, this study aims to investigate the impacts of intercropping soybeans with forage mulberry on soybean growth indices, yield components, photosynthesis, shading status, and root nodulation. Our findings would facilitate the formulation of a high-yield intercropping model specific to soybeans. These results will also enrich the theory of high-yield soybean cultivation, thereby offering a valuable reference point for the expansion of the soybean cultivation area.

2. Materials and Methods

Site description. A field experiment was conducted in 2021, 2022 and 2023 at the Experimental Base of Northeast Forestry University in Li Ming Town, Harbin City, Heilongjiang Province (45°44'N, 126°36'E), situated at an altitude of 95.0 m. The region experiences an annual average temperature of 5.5°C, a cumulative temperature above 10°C totalling 2700°C, a frost-free period of 115-130 days, and an annual average precipitation of 560 mm, predominantly falling between July and September, which accounts for over 59% of the total. The experimental site was characterised by a meadow soil composed of 1.42 g/kg of total nitrogen, 23.12 g/kg of total potassium, 0.62 g/kg of total phosphorus, 31.1 g/kg available matter content, and a pH of approximately 6.51.

Experimental design. The field experiment was carried out in a randomized block design in an irrigated 600 m² unobstructed arable land, with each intercropping plot measuring 6 m x 12 m. To explore the effects of intercropping with forage mulberry on soybean nodulation, the researchers deliberately excluded nitrogen fertilizer to establish a nitrogen - deficient condition. The treatments consisted of plots with sole-cropped soybean without nitrogen fertilizer (SN0) and intercrops of forage mulberry and soybeans without nitrogen fertilizer (IN0). In the intercropping system, the rows of plants were arranged in a 3:4 ratio, specifically three rows of forage mulberry followed by four rows of soybean. Both the soybean and forage mulberry were planted with a consistent row - to - row spacing of 0.65 m. The soybean was planted in rows spaced 8 cm apart to reach a final population of 190,000 plants/hm², was fertilized with 65 kg/hm² of P as triple superphosphate and 75 kg/hm² of K as KCl, whereas the forage mulberry was planted in rows 5 cm apart for a population of 300,000 plants/hm², received 35 kg/hm² of P and 60 kg/hm² of K from the same sources.

Growth performance and yield. Plant height was measured at soybean growth stages of beginning bloom (R1) as described previously[23] for both the SN0 and IN0 treatments. The diameter of the stem base and the length and width of the third compound leaf were measured using vernier callipers (DL91150, Deli, China), while the leaf area was calculated using the formula: Leaf area = length × width × 0.75, where 0.75 means leaf area coefficient of soybeans. The fresh and dry weights of the entire plant were recorded using an analytical balance (BSA223S, Sartorius, Germany). The plants were dried in an oven at 105°C for 1 h, then dried at 75°C until a constant weight was reached to determine the dry weight. The chlorophyll content was determined using the ethanol extraction colourimetric method, 0.5 grams of leaves were placed in a 50 mL centrifuge tube with 20 mL of 95% ethanol, sealed, and extracted in the dark at room temperature for 24 hours until the tissue turned white. The chlorophyll content was then measured using a UV-visible spectrophotometer at wavelengths of 665 nm and 649 nm. To determine the yield, the seeds were harvested according to the treatment plots and threshed, followed by measurement of the seed weight from each plot using an electronic balance in triplicates.

Photosynthetic gas exchange parameters. The photosynthetic parameters were measured using a Li-6400 portable photosynthetic system (LI-COR, Lincoln, USA) between 9:00 and 11:00 hours on the fully expanded functional leaves at the R1 stage before the first mowing of forage mulberry. All these measurements were repeated five times.

Measurement of nitrogenase activity. The nitrogenase activity of root nodules was indirectly determined by the acetylene reduction assay. Briefly, the roots from each plant with root nodules were collected and placed into sealed vials. Two millilitres of air were removed from each vial and replaced with an equal volume of acetylene. After incubating at 28°C for 2 h, the ethylene content was measured using gas chromatography (GC-2010, Shimadzu, Japan). The nitrogenase activity of nodules was expressed as the amount of ethylene produced in micromoles per hour per gram of sample. Each measurement was repeated five times.

Measurement of nitrogen metabolism. The third compound leaf of the soybean was collected at flowering stage (R1) in 2023. The concentrations of NH₄⁺ -N in the leaves were measured using phenol disulfonic acid colorimetry and potassium chloride extraction-indophenol blue colorimetry [24]. The activities of glutamine synthetase (GS) and glutamate synthase (GOGAT) in soybean leaves were refered to the method of Feng [25].

Transcriptome analysis. Soybean root samples were collected at the R1 stage in triplicates and stored in a cryostat. The root samples were washed with 75% ethanol to remove dust and soil and then frozen in liquid nitrogen for 1 h and stored at -80°C. The total RNA extracted from the root samples was used to construct the library and for transcriptome sequencing at the Hangzhou Kaitai Biotechnology.

Soil physical and chemical properties. The soil samples from soybean roots were collected at the R1 stages. The pH value of the soil was measured using a pH meter (delta 320, Mettler Toledo, Switzerland) with electrodes in a soil and water (1:2.5) suspension. The available nitrogen was determined using the alkali hydrolysis diffusion method [26], while the available phosphorus and potassium were determined as described by Lu [24].

Statistical analysis. All independent experiments were repeated at least five times. Statistical analyses were performed using SPSS 24.0 with a two-way analysis of variance to analyze the interaction of planting pattern and year on soil physico-chemical properties. And least significant difference (LSD) tests were used to compare differences between cropping patterns at the same year, as well as between years within the same cropping pattern differences among treatments. Canoco 4.5 was used to reveal the correlation among planting patterns, growth indices and soil properties though redundancy analysis (RDA). To establish a theoretical framework for intercropping that favours soybean cultivation and management, we employed structural equation modelling (SEM) to uncover the pathways and interaction mechanisms of various factors affecting soybean growth and nodulation traits and analyzed their direct and indirect effects. The interaction relationship was further analyzed by the partial least squares path model (PLS-PM). The path and determination coefficients (R²) in the path model were calculated by the plspm package of R 4.4.0.

3. Results

3.1. Growth Indices at the Bloom Stage of Soybean

Based on three years of data, there were no significant differences in the plant height and stem diameter of soybeans between sole-cropping and intercropping (Table 1). While, soybean intercropped with forage mulberry consistently achieved significantly higher values of leaf area, fresh weight per plant and dry weight per plant their corresponding values in sole-cropping.

Specifically, in 2021, 2022, and 2023, the leaf area of intercropped soybeans was 5.17% (P<0.05), 8.11% (P<0.05), and 5.79% (P<0.05) larger than that in sole - cropping, respectively. Similarly, the soybean fresh weight per plant of intercropped soybeans was 5.81% (P<0.05), 5.10% (P<0.05) and 4.63% (P<0.05) higher, and the dry weight was 4.19% (P<0.05), 2.69% (P<0.05), and 3.81% (P<0.05) higher than that in sole-cropping.

The results of the two-way ANOVA test showed that pattern has significant effect on soybean fresh weight per plant. While, year, pattern and year × pattern have no significant effects on plant height, stem diameters, leaf area and dry weight. Which clearly indicating intercropping with forage mulberry refrained from triggering spindly growth in soybeans and, to a certain degree, actively promoted their growth and development.

3.2. The Soybean Yield Components of Soybean

From 2021 to 2023, when compared with soybeans under sole cropping, the yield components of soybeans intercropped with forage mulberry showed significant increases (Table 2). Specifically, in 2021, 2022, and 2023, the number of pods per plant in intercropped soybeans was 5.48% (P < 0.05), 4.39% (P < 0.05), and 5.84% (P < 0.05) higher than that in sole - cropping, respectively. The number of seeds per plant in intercropped soybeans was 12.63% (P < 0.05), 11.39% (P < 0.05), and 8.47% (P < 0.05) higher than in sole - cropping. In the same vein, the grain yield per plant of intercropped soybeans was 7.41% (P < 0.05), 6.31% (P < 0.05), and 7.84% (P < 0.05) higher. The 100 - seed weight of intercropped soybeans was 12.17% (P < 0.05), 11.85% (P < 0.05), and 9.39% (P < 0.05) higher, and the overall grain yield was 8.39% (P < 0.05), 7.69% (P < 0.05), and 10.88% (P < 0.05) higher than that in sole - cropping.

The results of the two-way ANOVA test showed that the cropping pattern had a significant effect on the grain yield per plant, 100-seed weight, and fresh weight of soybean grains. In contrast, neither the year nor the interaction between the year and the cropping pattern had significant effects on these characteristics.

3.3. Root Nodule Traits at the Bloom Stage of Soybean

From 2021 to 2023, compared with soybeans in sole cropping, the number of root nodules, along with the fresh and dry weights of root nodules intercropped with forage mulberry treatment, were significantly increased (Table 3).

In particular, in 2021, 2022, and 2023, the number of root nodules per plant intercropped soybeans were 21.43% (P<0.05), 24.21% (P<0.05), and 25.79% (P<0.05) more than that in sole-cropping, respectively. Similarly, the fresh weight of nodules per plant of intercropped soybeans was 51.58% (P<0.05), 54.89% (P<0.05) and 56.78% (P<0.05) higher, and the dry weight was 40.21% (P<0.05), 44.00% (P<0.05), and 51.30% (P<0.05) higher than that in sole-cropping.

The results of the two - way ANOVA test showed that the cropping pattern had a significant effect on the number of nodules, the fresh weight of nodules, and the dry weight of nodules. Conversely, both the year and the interaction of the year and the cropping pattern had no significant impacts on these features.

3.4. Physical and Chemical Properties in Rhizosphere Soil

From 2021 to 2023, when compared with soybeans under sole cropping, the content of available nitrogen (AN) and the pH value in the treatment of soybeans intercropped with forage mulberry decreased significantly (Table 4). In contrast, the content of available phosphorus (AP), the content of available potassium (AK), and the soil water content (SWC) were higher in the intercropping system than in the sole - cropping system.

For example, in 2021, 2022, and 2023, the AN content in the soil of intercropped soybeans were 9.95% (P<0.05), 8.95% (P<0.05), and 9.31% (P<0.05) lower, respectively, and the pH value was 1.62% (P<0.05), 3.04% (P<0.05), and 2.44% (P<0.05) lower than that in sole-cropping. On the other hand, the AP content in the soil of intercropped soybeans was 16.72% (P<0.05), 16.48% (P<0.05) and 16.30% (P<0.05) higher, and the AK content was 19.77% (P<0.05), 19.53% (P<0.05), and 20.57% (P<0.05), and SWC was 13.51% (P<0.05), 11.28% (P<0.05), and 12.06% (P<0.05) higher than that in sole-cropping.

The results of the two-factor variance analysis indicated that the planting pattern had a significant influence on the pH value, AN, AP, AK and SWC in rhizosphere soil of the soybean. However, year and the interaction between the year and the planting pattern had no significant impact on the above - mentioned properties of the rhizosphere soil.

3.5. Photosythetic Parameters at Blooming Stage

Intercropping soybeans with forage mulberry significantly affected the main gas exchange parameters of soybean leaves, as presented in Figure 1. When compared to sole - cropping, the net photosynthetic rate (P_n) of soybeans in the intercropping system with forage mulberry increased by 24.86% (P < 0.01), and the intercellular CO₂ concentration (C_i) increased by 5.48% (P < 0.05). However, no significant differences were observed between the intercropping and sole - cropping systems in terms of chlorophyll content (P > 0.05), stomatal conductance (G_s) (P > 0.05), and transpiration rate (T_r) (P > 0.05).

3.6. Nitrogen Metabolism Changes in Soybean

This is example 1 of an The NH₄⁺ - N generated by nitrogen fixation in root nodules is assimilated through the GS/GOGAT cycle within the leaves of soybeans (Figure 2a). Intercropping soybeans with forage mulberry enhanced the nitrogenase activity of soybean roots (Figure 2b), along with the content of NH₄⁺ -N and the activities of glutamine synthetase (GS) in soybean leaves at the R1 stage (Figure 2c,d).

Compared to sole - cropping, intercropping with forage mulberry led to a significant increase in nitrogenase activity in nodules by 135.28% (P < 0.01), NH₄⁺ -N content in soybean leaves by 31.13% (P < 0.01), and GS activity in soybean leaves by 19.30% (P < 0.01). Conversely, the glutamate synthase activity (GOGAT) in soybean leaves within the intercropping system with forage mulberry was 24.64% (P < 0.001) lower than that in the sole - cropping system (Figure 2e).

3.7. Transcriptome Analysis

The transcriptome analysis identified 1,462 differentially expressed genes (DEGs) in soybean roots of intercropping with forage mulberry (IN0) compared to sole-cropping (SN0). Among these, 569 DGEs were up-regulated, and 893 DGEs were down-regulated, suggesting that intercropping with forage mulberry modified the expression of genes in soybean root (Figure 3a).

The DEGs were categorized into biological processes, cellular components, and molecular functions based on the gene ontology (GO) terms (Figure 3b). The DEGs enriched in the biological processes mainly included the phosphorus metabolic process, organic acid metabolic process, oxoacid metabolic process, response to organonitrogen compound and response to nitrogen compound. The enriched molecular functions mainly included oxidoreductase activity, organic acid transmembrane transporter activity and carboxylic acid transmembrane transporter activity (Figure 3b). The KEGG enrichment analysis showed that DEGs in SN0 and IN0 were highly enriched in plant-pathogen interaction and isoflavone biosynthesis pathways (Figure 3c).

3.8. Isoflavone Biosynthesis

Flavonoids serve as crucial signalling molecules in the symbiosis between soybean and rhizobium and are synthesised through the isoflavone biosynthesis pathway, which comprises the phenylpropane pathway from phenylalanine to coumarin-CoA, the flavonoid biosynthesis from p-coumarin-CoA to naringenin, and the isoflavone biosynthesis from naringenin to isoflavones (Figure 4). Compared with sole-cropping, intercropping with forage mulberry increased the expression of DEGs related to coenzyme A ligase (4-coumarate: coenzyme A ligase, 4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavone synthase (FNS), and 2-hydroxyisoflavanone dehydratase (HIDH). This indicates that intercropping with forage mulberyy activated the "isoflavone biosynthesis" pathway and enhanced the metabolism and secretion of flavonoids in the soybean roots compared to sole-cropping, thereby facilitating the biosynthesis of rhizobial nodulation factor (NF) to promote N fixation through root nodulation.

3.9. Redundancy Analysis (RDA)

Redundancy analysis (RDA) of the planting pattern, growth - related parameters, and soil properties indicated that intercropping with forage mulberry had the highest positive value along axis 1, while sole - cropping had the highest negative value along axis 1 (Figure 5). The genes related to the isoflavone pathway, nitrogenase activity, and the main soil properties, namely available phosphorus (AP), available potassium (AK), and soil water content (SWC), were positively correlated with intercropping and negatively correlated with sole - cropping. The growth indices, including grain yield per plant, 100 - seed weight, plant height, stem diameter, fresh weight per plant, dry weight per plant, leaf area, and the soil properties of available nitrogen (AN) and pH value, were positively correlated with sole - cropping and negatively correlated with intercropping.

3.10. Structural Equation Modeling (SEM)

Structural equation modeling (SEM) was used to analyze both the direct associations and indirect and causal linkages among variables for a better understanding (Figure 6). Intercropping with forage mulberry had highly significant negative effects on pH and AN with a path coefficient of -0.940 (P < 0.001) and significant positive effects on the AP, AK and SWC with path coefficients of 0.940 (P < 0.001). The contents of AP and AK had the most significant positive effects on the soybean isoflavone synthesis pathway, with path coefficients of 0.977 (P < 0.001). Moreover, isoflavone synthesis pathways were positively correlated with nodulation traits (β = 0.997, P < 0.001), while nodulation-mediated effects positively regulated nitrogen metabolism in soybean leaf (β = 0.991, P < 0.001), and the nitrogen metabolism positive correlated with the gas exchange parameters in soybean leaf (β = 0.986, P < 0.001). The results also showed a positive and insignificant correlation of the gas exchange parameters with yield (β = 0.595, P > 0.05), and a positive correlation of the intercropping with yield (β=0.960, P > 0.05).

GOF (Goodness of Fit)> 0.600 indicates a model with a good fit. The isoflavone pathway includes genes of 4CL5, 4CL7, CHI3, FNS and HIDH3, while the nodulation traits include the number, fresh weight and dry weight of nodules and nitrogenase activity. The photosynthetic parameters used include net photosynthetic rate (Pn), stomatal conductance (Gs), and transpiration rate (Tr). Other parameters used are soil pH, soil available nitrogen (AN), soil available phosphorus (AP) and soil available potassium (AK), soil water content (SWC). The nitrogen metabolism include the Nitrogenase activity, NH₄⁺-N content, glutamine synthetase activity (GS), glutamate synthase activity (GOGAT). The green arrows indicate positive effects, while the yellow arrows indicate negative effects; ∗indicates P < 0.05, and ∗∗∗indicates P < 0.001.

Figure 7. Underlying mechanisms of intercropping with forage mulberry-induced impacts on soybean growth and nodulation.

Figure 7. Underlying mechanisms of intercropping with forage mulberry-induced impacts on soybean growth and nodulation.

4. Discussion

4.1. Effects of Intercropping of Mulberry on the Growth of Soybean

Legumes are preferred for intercropping primarily due to their ability to symbiotically fix N, which they utilize and can also be used by adjacent crops [27,28]. While, intercropping with gramineous plants can reduce nitrogen (N) depletion in legumes, improve nodulation for N fixation, and increase N transfer, thereby decreasing the need for chemical fertilizers[8]. So intercropped with legumes greatly conforms to the strategy of reducing N input in agriculture [29]. However, the symbiotic N fixation is an energy-expensive process that depends greatly on shoot photosynthetic processes [30]. The light available to plants is a primary influencing factor in photosynthesis since it provides energy. Moreover, light also serves as an essential signal for biological nitrogen fixation [31]. On the other hand, shading brought by intercropping directly leads to a decline in soybean’s photosynthetic parameters [14], thinning of the stems, elongation of the internodes, increase in plant height, lodging and reduced yields [17,20]. In this study, intercropping with forage mulberry did not hinder the plant height and stem diameter of soybeans. In contrast, it enhanced the chlorophyll content and net photosynthetic rate, thereby increasing the soybean yield. This indicated that the regularly mowed forage mulberry did not exert shading effects on soybeans during the intercropping symbiosis period. To some extent, intercropping with forage mulberry also had beneficial effects on soybean growth and yield.

4.2. Effects of Intercropping on Nodulation Related Indicators

Intercropping with gramineous plants creates a nitrogen - depleted environment. This reduces the "nitrogen repression" on the root systems of leguminous plants by competing for soil nitrogen and also promotes nodulation[32]. Our previous studies found that the contents of available nitrogen content (AN) in the rhizospheric soil of the legumes intercropped with forage mulberry declined, while the number of nodules increased[12,33,34]. Similarly, in this current study, the AN content in the rhizospheric soil of soybeans intercropped with forage mulberry also declined over the three - year period. Therefore, it is likely that the forage mulberry enriched protein was mowed 2-3 times just like gramineous plants to take away a large amount of N elements from the soil, thereby creating an N-deficient environment. When the roots of leguminous plants sense the deficiency in soil N, they secrete flavonoid substances as signals into the soil, which plays the role of a key signaling molecule in symbiosis[35,36].

Intercropping can promote the expression level of CHI genes in pea roots, thereby enhancing the secretion of flavonoids in leguminous crops, which is beneficial in the nodulation of legumes[37]. In this study, intercropping with forage mulberry stimulated the expression of genes encoding the key enzymes, including 4CL, CHS, CHI, FNS, and HIDH, which are involved in the isoflavone biosynthesis in soybean roots. Correspondingly, the nodule traits of soybeans intercropped with forage mulberry are significantly superior to those of soybeans under monoculture. However, the N-fixation in leguminous plants is limited by the contents of the soil P and K[20], the legumes are stimulated to adopt potential P uptake strategies via root-secreting organic acid to meet their high P demand for N fixation[37]. In the present study, intercropping with forage mulberry decreased the rhizospheric soil pH of intercropped soybeans but increased the content of soil AP and AK. Additionally, the gene ontology analysis also showed the significant enrichment of organic acid and ion transmembrane transporter proteins in the intercropped soybeans. Therefore, these findings suggest that intercropping with forage mulberry could potentially stimulate soybean roots to secrete organic acids, consequently enhancing soil phosphorus activity. Meanwhile, ample soil moisture acts as an indispensable medium for the activated phosphorus to penetrate root. Moreover, soil water content (SWC) significantly impacts the nodulation traits of leguminous plants [38]. Specifically, inadequate SWC can impede nodule formation and reduce nodule activity. Shen et al[39] reported that intercropping with maize plants can enhance the SWC in soybean rhizosphere soils. Similarly, intercropping with forage mulberry has been shown to increase the SWC in the rhizosphere soil of soybean roots. The adequate SWC may facilitate the utilization of nitrogen fixed by root nodules, transport it to the leaves, and promote nitrogen metabolism in the leaves.

4.3. Effects of Intercropping on NH₄⁺-N Assimilation in Soybean Leaf

Previous research has demonstrated that the most vigorous stage of nitrogen fixation, transport, and assimilation in soybean leaves occurs at the flowering stage[40]. Therefore, aiming to build upon previous findings, in this experiment, all indicators associated with nitrogen metabolism were assayed during the flowering stage. And, we observed that the nitrogenase activity in soybean root nodules and the NH₄⁺ -N content in leaves in the intercropping system with forage mulberry were higher than those in the sole - cropping system. To investigate the effects of intercropping with forage mulberry on soybean nodulation, nitrogen fertilizer was deliberately omitted to create a nitrogen - limited environment. It is highly likely that the main source of ammonium ions (NH₄⁺) in soybean leaves is the nitrogen fixed via soybean nodulation. Meanwhile, the conversion of NH₄⁺ into glutamine and glutamate depends on the interaction between glutamine synthetase (GS) and glutamate synthase (GOGAT) in the GS - GOGAT cycle[41]. In this experiment, the activity of GS in soybean leaves was found to be higher in the intercropping system with forage mulberry than in the sole - cropping system. Conversely, the activity of GOGAT exhibited an opposite trend. There was no correlation found among the NH₄⁺ content, GS activity, and GOGAT activity in soybean leaves. A possible reason for this is that the function of GS/GOGAT may not be restricted to nitrogen assimilation. The products of GS/GOGAT, glutamine and glutamate, not only participate in amino acid metabolism but also potentially act as signal molecules in various pathways[42]. There are still many gaps in a better understanding of the interactions of GS/GOGAT and the concentration of NH4+ in soybean intercropped with forage mulberry.

4.4. Relationships Between Soybean Growth, Photosynthesis, Soil Physicochemical Indices, Nodulation and Yield

In recent years, there has been an exponential increase in studies examining the effects of agronomic practices on soil microbial communities[43]. However, it remains a challenge to elucidate the relationships between plant phenotypic traits, soil physicochemical properties, plant yield, and plant gene expression, as well as their impact pathways under intercropping conditions[44]. Combining RDA and SEM can help explore the correlations between variables and then delve deeply into their specific cause - effect relationships and the underlying mechanisms [45]. In this study, RDA indicated that on the same coordinate axis, the genes involved in soybean root isoflavone biosynthesis, nodule traits, along with soil AP, AK, and SWC, were positively associated with intercropping patterns but negatively correlated with sole - cropping. Meanwhile, the main soybean growth indices, along with soil pH and AN, were positively related to sole - cropping and negatively related to intercropping. Previous studies have demonstrated that intercropping soybeans with maize exerts a detrimental influence on the photosynthetic traits, soil nitrogen content, and yield of soybeans[44]. On the contrary, our SEM results indicate that the soil P and K content as the main contributors to soybean yield, had the most significant positive effects on isoflavone synthesis and nodulation traits, while the nodulation traits positively related to nitrogen metabilism and photosynthesis in leaves, which in turn affected the soybean yield. Based on the above discussion, we have elucidated the fundamental mechanisms by which forage mulberry promotes the growth of soybeans and established a theoretical framework for enhancing soybean nodulation in the intercropping system (Figure. 7).

5. Conclusions

The presented data elucidate the fundamental mechanisms of forage mulberry in facilitating the growth of soybeans, and formulate theoretical framework for enhancing nodulation in intercropping system. In this experiment, intercropped with forage mulberry did not induce shade stress on soybean. Additionally, intercropping significantly reduced the content of AN and the value of pH in the rhizosphere soil, while enhanced the levels of AP, AK and SWC. Moreover, the intercropping activated the expression of genes encoding key enzymes involved in the “isoflavone biosynthesis“ pathway, including 4CL, CHS, CHI, FNS, and HIDH. Furthermore, the intercropping of forage mulberry significantly increased the number of nodules, and nitrogenase activity in soybean roots, which collectively contributed to higher soybean yields.